金融MATLAB實(shí)驗(yàn)報(bào)告三答案.詳解

1、安徽財(cái)經(jīng)大學(xué)金融 證券實(shí)驗(yàn)室實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)課程名稱 金融 MATLAB開課系部金融學(xué)院班級(jí)學(xué)號(hào)姓名指導(dǎo) 教 師年月日1實(shí)驗(yàn)名稱MATLAB 金融數(shù)量分析學(xué)院學(xué)號(hào)姓名2實(shí)驗(yàn)準(zhǔn)備學(xué)會(huì)使用 MATLAB金融工具箱進(jìn)行金融數(shù)量分析，如：期權(quán)定價(jià)分析、投資組合績(jī)效分析、收益和風(fēng)險(xiǎn)的計(jì)算、有效前沿的計(jì)算、固定收益實(shí)驗(yàn)?zāi)康淖C券的久期和凸度計(jì)算、利率的期限結(jié)構(gòu)、技術(shù)指標(biāo)計(jì)算等。在下述 6 個(gè)主題中任選 3 個(gè)主題，使用 MATLAB金融工具箱進(jìn)行數(shù)量分析，數(shù)據(jù)來源自行在網(wǎng)上搜尋，要求是 2012 年之后的數(shù)據(jù)。（可參照各章的例題）1. 期權(quán)定價(jià)分析 （第 10 章）2. 收益、風(fēng)險(xiǎn)和有效前沿的計(jì)算 （第 12

2、章）3. 投資組合績(jī)效分析 （第 13 章）4. 固定收益證券的久期和凸度計(jì)算 （第 17 章）5. 利率的期限結(jié)構(gòu) （第 18 章）6. 技術(shù)指標(biāo)分析 （第 22 章）本實(shí)驗(yàn)報(bào)告不指定具體的題目，請(qǐng)大家自行設(shè)定，同學(xué)相互之間不要出實(shí)驗(yàn)設(shè)計(jì)方案現(xiàn)雷同。一、期權(quán)定價(jià)分析1.black-scholes方程求解例 1：假設(shè)歐式股票期權(quán)，六個(gè)月后到期，執(zhí)行價(jià)格 90 元，現(xiàn)價(jià)為 102 元，無股利支付，股價(jià)年化波動(dòng)率為 55%，無風(fēng)險(xiǎn)利率為 8%，計(jì)算期權(quán)價(jià)格。解： clearPrice=102;>> Strike=90;>> Rate=0.08;>> Time=6

3、/12;>> Volatility=0.55;CallDelta, PutDelta = blsprice(Price,Strike, Rate,Time, Volatility)計(jì)算結(jié)果：CallDelta =23.5648PutDelta=8.03582.期權(quán)價(jià)格與波動(dòng)率關(guān)系分析Price=102;>> Strike=90;>> Rate=0.08;>> Time=6/12;Volatility=0.08:0.01:0.5;>> N=length(Volatility)Call=zeros(1,N);Put=zeros(1,N);

4、for i=1:NCall(i), Put(i) = blsprice(Price,Strike, Rate,Time, Volatility(i);N=43endplot(Call,'b-');hold onplot(Put,'b');xlabel('Volatility')ylabel('price')legend('Call','Put')33. 計(jì)算期權(quán) Delta 。例 2. 假設(shè)歐式股票期權(quán)，六個(gè)月后到期，執(zhí)行價(jià)格 90 元，現(xiàn)價(jià)為 102 元，無股利支付，股價(jià)年化波動(dòng)率為 55%，無風(fēng)

5、險(xiǎn)利率為 8%，計(jì)算期權(quán) Delta。解： clearPrice=102;>> Strike=90;>> Rate=0.08;>> Time=6/12;>> Volatility=0.55;CallDelta,PutDelta= blsdelta(Price,Strike,Rate,Time, Volatility)計(jì)算結(jié)果： CallDelta=0.7321PutDelta=-0.26794. 利用不同的 Price 與 Time 計(jì)算 Detla 三維關(guān)系。>> Price=60:1:102;>> Strike=90;

6、Rate=0.08;4>> Time=(1:1:12)/12;>> Volatility=0.55;>> Price,Time=meshgrid(Price,Time);Calldelta,Putdelta= blsdelta(Price,Strike,Rate,Time, Volatility);>> mesh(Price, Time, Putdelta); xlabel('Stock Price '); ylabel('Time (year)'); zlabel('Delta');>>

7、;5. B-S 公式隱含波動(dòng)率計(jì)算例 3：假設(shè)歐式股票期權(quán)，一年后，執(zhí)行價(jià)格 99 元，現(xiàn)價(jià)為 105 元，無股利支付，股價(jià)年化波動(dòng)率為 40%，無風(fēng)險(xiǎn)利率為 10%，則期權(quán)價(jià)格為：解： clear>> Price=105;>> Strike=99;>> Rate=0.1;>> Time=1;>> CallValue=15;>>CallVolatility=blsimpv(Price,Strike,Rate,Time,CallValue, ,5'Call')計(jì)算結(jié)果：CallVolatility=NaN&g

8、t;> PutValue=7;>> PutVolatility= blsimpv(Price,Strike,Rate,Time,PutValue,'Put')PutVolatility=0.34556. 期權(quán)二叉樹模型的計(jì)算例：假設(shè)歐式股票期權(quán)，三個(gè)月后到期，執(zhí)行價(jià)格85 元，現(xiàn)價(jià)為 95 元，無股利支付，股價(jià)年化波動(dòng)率為 60%，無風(fēng)險(xiǎn)利率為 10%。解： clear>> Price=95;>> Strike=85;>> Rate=0.1;>> Time=4/12;>> flag=1;>> Increment=1/12;>> Volatility=0.6;>> AssetPrice,OptionValue= binprice(Price,Strike,Rate,Time,Increment,Volatility,flag)計(jì)算結(jié)果：AssetPrice=95.0000112.9654134.3283159.7312189.9379079.891795.0000112.9654134.32830067.186179.891795.0000600056.50