matlab金融分析實(shí)驗(yàn)

matlab金融分析實(shí)驗(yàn)

12.1知識(shí)要點(diǎn)和背景：最優(yōu)投資組合及其計(jì)算

12.1.1包含無風(fēng)險(xiǎn)證券的投資組合
12.1.2無風(fēng)險(xiǎn)證券投資組合的計(jì)算
mzy1.m
【?? ra=0.12;rb=0.08;rf=0.06;rp=0.1;
Aeq=[ra,rb,rf;1 1 1];beq=[rp 1]';lb=[0 0 -100]';
x0=[1 1 1]'/3;
[x,fval]=fmincon('mzy1f',x0,[],[],Aeq,beq,lb,[],[]),
pause;x0=x
[x,fval]=fmincon('mzy1f',x0,[],[],Aeq,beq,lb,[],[]),
?x0=[6/19,20/19,-7/19]??? 】
mzy1f.m
【??? function f=mzy1f(x)
?v(1,1)=10;v(1,2)=0;v(1,3)=0;v(2,1)=0;v(2,2)=1;v(2,3)=0;v(3,:)=0;
? f=x'*v*x;????? 】
12.2?機(jī)會(huì)的價(jià)值
12.2.1 簡(jiǎn)單二項(xiàng)式模型機(jī)會(huì)價(jià)值
?以上計(jì)算方法是從期末開始向期初倒推，稱為后向式動(dòng)態(tài)規(guī)劃。Maltab程序?yàn)椋?br>mzy2.m?????????? % 一期二項(xiàng)式方法計(jì)算Call options
【??? clear;
S0=input('輸入當(dāng)前股票價(jià)格:');
E=input('輸入?yún)f(xié)定執(zhí)行價(jià)格:');
u=input('輸入上升比例u:');
d=input('輸入下降比例d:');
r=input('輸入無風(fēng)險(xiǎn)利率r:');
if any(S0 <= 0 | E<=0 | r<0 | u<=0 | d<0 | u<1+r | 1+r??? error(sprintf('輸入 S0, E, u, d 和 r > 0. 輸入 u>(1+r)>d.'))
??? break
end
fu=max(u*S0-E,0)；fd=max(d*S0-E,0)；
p=(1+r-d)/(u-d)；f=(p*fu+(1-p)*fd)/(1+r);
disp('所求的合約價(jià)值為:');f? 】

12.2.2? 兩期二項(xiàng)式模型機(jī)會(huì)價(jià)值
mzy3.m??????? (一期二項(xiàng)式方法計(jì)算Call options)
【?? clear;
SO=input('輸入當(dāng)前股票價(jià)格:');
E=input('輸入?yún)f(xié)定執(zhí)行價(jià)格:');
u=input('輸入上升比例u:');
d=input('輸入下降比例d:');
r=input('輸入無風(fēng)險(xiǎn)利率r:');
if any(SO <= 0 | E<=0 | r<0 | u<=0 | d<0 | u<1+r | 1+r?? error(sprintf('輸入 SO, E, u, d 和 r > 0. 輸入 u>(1+r)>d.'))
?? break
end
fuu=max(u*u*SO-E,0)；fud=max(u*d*SO-E,0)；
fdd=max(d*d*SO-E,0)；p=(1+r-d)/(u-d)；
fu=(p*fuu+(1-p)*fud)/(1+r)；fd=(p*fud+(1-p)*fdd)/(1+r)；
f=(p*fu+(1-p)*fd)/(1+r);
disp('所求的合約價(jià)值為:'); f?????? 】
12.2.3 觀察與思考

200多個(gè)MATLAB經(jīng)典教程和MATLAB論文請(qǐng)查看：matlab教程