安徽财经大学金融证券实验室实验报告实验课程名称 《金融MATLAB》开课系部 金融学院 班 级 学 号 姓 名 指导教师 年月日实验名称MATLAB金融数量分析学院学号姓名实验准备实 验 目 的学会使用MATLAB金融工具箱进行金融数量分析，如：期权定价分析、 投资组合绩效分析、收益和风险的计算、有效前沿的计算、固定收益 证券的久期和凸度计算、利率的期限结构、技术指标计算等实 验 设 计 方 案在下述6个主题中任选3个主题，使用MATLAB金融工具箱进行数量分析， 数据来源自行在网上搜寻，要求是2012年之后的数据可参照各章的例题）1. 期权定价分析（第10章）2. 收益、风险和有效前沿的计算（第12章）3•投资组合绩效分析（第13章）4•固定收益证券的久期和凸度计算（第17章）5•利率的期限结构（第18章）6•技术指标分析（第22章）本实验报告不指定具体的题目，请大家自行设定，冋学相互之间不要出现 雷同一、期权定价分析l.black-scholes 方程求解例1：假设欧式股票期权，六个月后到期，执行价格90元，现价为102元，无股利支付，股 价年化波动率为55%，无风险利率为8%，计算期权价格。
解：clearPrice=102;>> Strike=90;>> Rate=0.0 8;>> Time=6/12;>> VOlatility=O.55;[CallDelta, PutDelta] = blsprice(Price, Strike, Rate, Time, Volatility)计算结果：CallDelta =23.5648PutDelta =8.03582•期权价格与波动率关系分析Price=102;>> Strike=90;>> Rate=0.0 8;>> Time=6/12;Vblatility=0.08:0.01:0.5;>> N=length(Volatility)Call=zeros(1,N); Put=zeros(1,N);for i=1:N[Call(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')^Figure 1 □旦区3. 计算期权Delta。
例2.假设欧式股票期权，六个月后到期，执行价格90元，现价为102元，无股利支付，股 价年化波动率为55%，无风险利率为8%,计算期权Delta解：clearPrice=102;>> Strike=90;>> Rate=0.0 8;>> Time=6/12;>> Volatility=0.55;[CallDel ta, Put Del ta] = blsdel ta(Price, St rike, Rate, Time, Vola tility)计算结果：CallDelta =0.7321PutDelta =-0.26794. 利用不同的Price与Time计算Det la三维关系>> Price=60:1:102;>> Strike=90;Rat e=0.08;>> Time=(l:l:12)/12;>> Volatility=0.55;>> [Price,Time]二meshgrid(Price,Time);[Calldel ta. Put del ta] = blsdel ta(Price, St rike, Rate, Time, Vola tility); >> mesh(Price, Time, Putdelt a);xlabel('S tock Price '); ylabel('Time (year)');zlabel('De lt a');>>股价年化[],[],5. B-S公式隐含波动率计算例3：假设欧式股票期权，一年后，执行价格99元，现价为105元，无股利支付, 波动率为40%，无风险利率为10%,则期权价格为：解：clear>> Price=105;>> Strike=99;>> Ra te=0.1;>> Time=1;>> CallValue=15;>> CallVolatility 二 blsimpv(Price, Strike, Rate, Time, CallValue,[],{'Call'})计算结果：CallVolatility 二NaN>> PutValue=7;>> PutVolatility 二 blsimpv(Price, Strike, Rate, Time, PutValue, [], [], [], {'Put'}) PutVolatility 二0.34556. 期权二叉树模型的计算例：假设欧式股票期权，三个月后到期，执行价格85元，现价为95元，无股利支付，股价 年化波动率为60%,无风险利率为10%。
解：clear>> Price=95;>> Strike=85;>> Ra te=0.1;>> Time=4/12;>> flag=1;>> Increment二1/12;>> Vola tilit y=0.6;>> [AssetPrice, OptionValue] = binprice(Price, Strike, Rate, Time, Increment, Vola tility, flag)计算结果：AssetPrice =95.0000112.9654134.3283159.7312189.9379079.891795.0000112.9654134.32830067.186179.891795.000000056.501267.1861。