期权定价实验报告M黄清霞.doc
- 文档编号:2502997
- 上传时间:2023-05-03
- 格式:DOC
- 页数:9
- 大小:237.50KB
期权定价实验报告M黄清霞.doc
《期权定价实验报告M黄清霞.doc》由会员分享,可在线阅读,更多相关《期权定价实验报告M黄清霞.doc(9页珍藏版)》请在冰点文库上搜索。
广东金融学院实验报告
课程名称:
金融工程
实验编号
及实验名称
期权定价模型及数值方法综合实验
系别
应用数学系
姓名
黄清霞
学号
101613110
班别
1016131
实验地点
实验日期
2013-06-01
实验时数
指导教师
张学奇
其他成员
黄冬璇、马燕纯
成绩
一、实验目的及要求
1.实验目的
(1)通过期权定价模型与数值方法综合实验,使学生加深对BSM期权模型的理解;
(2)熟练掌握运用Matlab计算欧式期权价格实际应用方法;
(3)熟练掌握运用Matlab软件计算美式期权价格的有限差分法、蒙特卡罗模拟法。
(4)培养学生运用软件工具解决期权定价问题的应用和动手能力。
2.实验要求
实验以个人形式进行,要求每位实验人员按照实验指导书,在实验前做好实验原理复习工作,实验软件的熟悉工作。
实验报告要包括:
实验原理、实验工具、实验程序与实验结论。
实验内容要详实和规范,实验过程要完整和可验证,实验结果要准确。
二、实验环境及相关情况
实验设备:
实验中心和个人计算机
实验软件:
Matlab软件。
实验资料:
期权定价模型及数值方法综合实验指导书。
三、实验内容、步骤及结果
(一)基于Matlab的无收益资产的欧式期权定价实验
A.实验原理
1.参量与符号
(1):
标的资产的价格;
(2):
行权价格;(3):
到期期限;
(4):
标的资产价格波动率;(5):
连续复利的无风险利率;
2.无收益资产欧式期权定价公式
无收益欧式看涨期权的定价公式
无收益资产欧式看跌期权的定价公式
其中,
B.实验算例
算例:
股票的价格为100,股票的波动率标准差为0.25,无风险利率为5.47%,期权的执行价格为100,存续期为0.25年,试计算该股票欧式期权的价格。
C.实验过程
在Matlab中计算欧式期权价格的函数为blsprice。
调用方式为
[Call,Put]=blsprice(Price,Strike,Rate,Time,Volatility,Yield)
输入参数
Price%标的资产价格
Strike%执行价格
Rate%无风险利率
Time%距离到期日的时间,即期权的存续期
Volatility%标的资产的标准差
Yield%标的资产红利率
输出参数
Call%欧式看涨期权价格
Put%欧式看跌期权价格
MATLAB中的计算过程和结果如下:
(请将运算过程和结果粘贴下面)
>>[Call,Put]=blsprice(100,100,0.0547,0.25,0.25,0)
Put=
4.3001
Call=
5.6583
所以此算例中欧式看涨期权价格为5.6583,欧式看跌期权价格为4.3001。
(二)基于Matlab的期权定价二叉树方法实验
A.实验原理
1.二叉树模型结构
对于多时段二叉树模型,在时刻,证券价格有中可能,他们可用符号表示为,其中。
应用多时段二叉树模型来表示证券价格变化的树型结构如下图所示。
2.二叉树模型中参数的确定
;;
3.无收益欧式期权二叉树模型定价公式
(1)对于无收益欧式看涨期权,节点的期权价值为
;,
最后一列节点的期权价值为
(2)对于无收益欧式看跌期权,节点的期权价值为
;,
最后一列节点的期权价值为
B.实验算例
算例:
考虑一个不分红利5个月欧式看涨期权:
股票价格为50,执行价格为50,无风险利率为10%,波动率为40。
试构造二叉树模型,确定期权的价格并与求解公式所得解进行比较。
C.实验过程
在Matlab中可以直接利用二叉树定价函数确定期权价格,函数名称为binprice。
调用方式为
[AssetPrice,OptionValue]=binprice(Price,Strike,Rate,Time,Increment,Volatility,Flag,DividendRate,Dividend,ExDiv)
输入参数
Price%标的资产价格
Strike%执行价格
Rate%无风险利率
Time%距离到期日的时间,即期权的存续期
Increment%时间增量
Volatility%标的资产的标准差
Flag%确定期权种类,看涨期权为1,看跌期权为0
DividendRate%红利发放率(可选项)
Dividend%标的资产价外红利金额(可选项)
ExDiv%标的资产除息日期(可选项)
输出参数
Price%二叉树每个节点价格
Option%期权在每个节点的现金流
MATLAB中的计算过程和结果如下:
(请将运算过程和结果粘贴下面)
1.二叉树定价函数确定期权价格
>>[AssetPrice,OptionValue]=binprice(50,50,0.1,5/12,1/12,0.4,1)
AssetPrice=
50.000056.120062.989270.699179.352889.0656
044.547450.000056.120062.989270.6991
0039.689444.547450.000056.1200
00035.361139.689444.5474
000031.504935.3611
0000028.0692
OptionValue=
6.35959.873414.859721.525629.767739.0656
02.84934.90668.248113.404120.6991
000.77941.54913.07916.1200
000000
000000
000000
2.求解公式确定期权价格
q=
0.4927
p=
0.5073
a=
1.0084
d=
0.8909
u=
1.1224
BinTree=
50.000056.120062.989270.699179.352889.0656
044.547450.000056.120062.989270.6991
0039.689444.547450.000056.1200
00035.361139.689444.5474
000031.504935.3611
0000028.0692
BinPrice=
6.35959.873414.859721.525629.767739.0656
02.84934.90668.248113.404120.6991
000.77941.54913.07916.1200
000000
000000
000000
所以可见二叉树定价函数确定期权价格与求解公式确定期权价格所得解是一样的。
(三)基于Matlab的期权定价的显性有限差分法实验
A.实验原理
1.BSM期权价格微分方程
2.看跌期权定价的显性差分法
(1)BSM微分方程中偏导数的差分近似
;;
(2)差分方程
其中
(3)边界条件
T时刻看跌期权的价值为
下边界上期权价值为
上边界上期权价值为
(4)期权的价值
对于差分方程和边界条件
;
解出每个的期权值,然后再与每个格点的期权内在价值进行比较,判断是否提前执行,从而得到时刻每个格点的期权价格。
依此类推,可以计算出,当等于初始资产价格时,该格点对应的就是所求的期权价值。
B.实验算例
算例:
已知股票价格为50,美式看跌期权执行价格为50,到期日为5个月,无风险利率为10%,波动率标准差为0.4。
试用有限差分法确定期权的价格。
C.实验过程
1.用显性差分法求解美式看跌期权
s0=50;k=50;r=0.1;sigma=0.4;T=5/12;
dt=T/10;ds=5;
Smax=100;
M=round(Smax/ds);%对Smax/ds取整运算
N=round(T/dt);
ds=Smax/M;%重新确定股票价格步长
dt=T/N;%重新确定时间步长
S=0:
ds:
Smax
S=S'
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
veti=1:
N;
vetj=1:
M+1;
a=1/(1+r*dt)*(-1/2*r*vetj*dt+1/2*sigma^2*vetj.^2*dt);
b=1/(1+r*dt)*(1-sigma^2*vetj.^2*dt);
c=1/(1+r*dt)*(1/2*r*vetj*dt+1/2*sigma^2*vetj.^2*dt);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
L=zeros(M-1,M+1);
forj=1:
M-1
L(j,j)=a(j);L(j,j+1)=b(j);L(j,j+2)=c(j);
end
f=zeros(M+1,N+1);
f(:
N+1)=max(k-S,0);
f(1,:
)=k
f(M+1,:
)=0
%递推求解期权价格
fori=N:
-1:
1
F=L*f(:
i+1);
f(2:
M,i)=max([F,k-S(2:
M)],[],2)
end
MATLAB中的计算过程和结果如下:
(请将运算过程和结果粘贴下面)
程序:
functionaa_3
s0=50;%股价
k=50;%执行价
r=0.1;%无风险利率
T=5/12;%存续期
sigma=0.3;%股票波动率
Smax=100;%确定股票价格最大价格
ds=2;%确定股价离散步长
dt=5/1200;%确定时间离散步长
M=round(Smax/ds);%计算股价离散步数,对Smax/ds取整运算
ds=Smax/M;%计算股价离散实际步长
N=round(T/dt);%计算时间离散步数
dt=T/N;%计算时间离散实际步长
matval=zeros(M+1,N+1);
vets=linspace(0,Smax,M+1);%将区间[0,Smax]分成M+1段
veti=0:
N;
vetj=0:
M;
matval(:
N+1)=max(k-vets,0);%建立偏微分方程终值条件
matval(1,:
)=k*exp(-r*dt*(N-veti));%建立偏微分方程边界条件
matval(M+1,:
)=0;
a=0.5*dt*(sigma^2*vetj-r).*vetj;%确定迭代矩阵L的元素
b=1-dt*(sigma^2*vetj.^2+r);
c=0.5*dt*(sigma^2*vetj+r).*vetj;
L=zeros(M-1,M+1)
fori=2:
M
%建立递推关系
L(i-1,i-1)=a(i);L(i-1,i)=b(i);L(i-1,i+1)=c(i);
end
fori=N:
-1:
1
matval(2:
M,i)=L*matval(:
i+1);
end
%插值求期权价格
jdown=floor(s0/ds);
jup=ceil(s0/ds);
ifjdown==jup
price=matval(jdown+1,1)+(s0-jdown*ds)*(matval(jup+1,1)-matval(jup+1,1))/ds
end
price=
2.8288
结果:
所以此算例中用有限差分法确定美式看跌期权的价格为2.8。
(四)基于Matlab的期权定价的蒙特卡洛模拟实验
A.实验原理
期权定价蒙特卡罗模拟的步骤:
1.从初始时刻的标的资产价格开始,直到到期为止,为取一条在风险中性世界中跨越整个有效期的随机路径。
(1)将变化路径的有效期等分为个长度为的时间段,由公式
或
计算时刻的值。
(2)接着计算时的值,这样通过个正态分布的随机数就可以组件一条资产价格的蒙特卡洛模拟样本路径。
2.计算出这条路径下期权的回报以及估计期权价值的标准差。
3.重复第一步和第二步,得到许多样本结果,即风险中性世界中期权回报的大量可能取值。
4.计算这些样本回报的均值,得到风险中性世界中预期的期权回报值。
5.用无风险利率贴现,得到这个期权的估计价值。
B.实验算例
算例:
假设股票价格服从几何布朗运动,股票现价为50美元,欧式期权执行价格为52美元,股票价格的波动率为0.4,无风险利率为0.1,到期日是5个月后,在这期间,该标的股票不进行任何的分红,试用蒙特卡洛模拟(分别模拟1000次和10000次)计算该期权的价格。
C.实验过程
在Matlab中用蒙特卡洛模拟计算该期权价格的程序为:
functioneucall=blsmc(s0,K,r,T,sigma,Nu)
randn('seed',0);%定义随机数发生器种子是0,这样可以保证每次模拟结果相同
nuT=(r-0.5*sigma^2)*T;
sit=sigma*sqrt(T);
discpayoff=exp(-r*T)*max(0,s0*exp(nuT+sit*randn(Nu,1))-K);%期权到期时现金流
[eucall,varprice,ci]=normfit(discpayoff)
MATLAB中的计算过程和结果如下:
(请将运算过程和结果粘贴下面)
>>blsmc(50,52,0.1,5/12,0.4,1000)
eucall=
5.4445
varprice=
9.1361
ci=
4.8776
6.0115
%eucall:
欧式看涨期权的价格
%varprice:
模拟期权价格的方差
%ci:
95%概率保证的期权价格区间
>>blsmc(50,52,0.1,5/12,0.4,10000)
ci=
4.9587
5.3089
eucall=
5.1338
varprice=
8.9335
所以蒙特卡洛模拟分别模拟1000次和10000次计算出该期权的价格分别为5.4445和5.1338。
四、实验总结(请将实验心得体会和实验改进意见填写在下面)
这次实验让我学会了用Matlab软件二叉树法、有限差分法、蒙特卡罗模拟法计算期权价格,在实验的过程中我遇到了许多问题,但是经过努力都一一解决了。
通过这次期权定价模型与数值方法综合实验,使我加深对BSM期权模型的理解,同时培养了我运用软件工具解决期权定价问题的应用和动手能力。
五、教师评语
评语
评语等级
优
良
中
及格
不合格
1.实验态度认真,实验目的明确
2.实验方案、程序设计合理
3.实验过程(实验步骤详细,记录完整,数据合理)
4.实验结论正确,分析透彻
5.实验报告独立完成,无抄袭现象,并按时提交,格式规范,文字叙述流畅,逻辑性强
综合评定:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 期权 定价 实验 报告 黄清霞