第二次上机作业Word文档格式.docx
- 文档编号:17844754
- 上传时间:2022-12-11
- 格式:DOCX
- 页数:12
- 大小:181.58KB
第二次上机作业Word文档格式.docx
《第二次上机作业Word文档格式.docx》由会员分享,可在线阅读,更多相关《第二次上机作业Word文档格式.docx(12页珍藏版)》请在冰豆网上搜索。
得图如下
使用surf作图
surf(z);
使用waterfall作图
waterfall(z);
使用plot3作图
plot3(x1,y1,z);
使用meshz作图
meshz(z);
1.2进行数值积分
dblquad(@(x,y)(cos(x)+sin(y)),0,1,1,2)
ans=
1.7979
我用的是matlab7.0,所以这个可以实现
fun=@(x,y)(cos(x)+sin(y));
dblquad(fun,0,1,1,2)
1.7979
③解
fun=cos(x)+sin(y);
dblquad(@(x,y)fun,0,1,1,2)
1.8444
PS:
不知道答案为什么会不一样
④解
fun=inline('
cos(x)+sin(y)'
'
x,y'
)
fun=
Inlinefunction:
fun(x,y)=cos(x)+sin(y)
jf=dblquad(fun,0,1,1,2)
这个是抄老师的,但是我实现不了,不知道是不是抄错了
⑤解
jf=dblquad(inline('
),0,1,1,2)
这个也是,实现不了
⑥解用int
③
f=cos(x)+sin(y)
f=
cos(x)+sin(y)
g=int(f,x,0,1)
g=
sin
(1)+sin(y)
h=int(g,y,1,2)
h=
cos
(1)-cos
(2)+sin
(1)
这个是上机的时候实现的,但是不知道为什么重新做的时候就实现不了了。
。
1.3使用MC求积分
n=1000;
x=unifrnd(0,1,n,1);
y=unifrnd(1,2,n,1);
z=unifrnd(0,2,n,n);
m=0;
fori=1:
n
forj=1:
z1=cos(x(i))+sin(y(j));
ifz(i,j)<
z1
m=m+1;
end
p=m/n^2
p=
0.9015
f=2*p
1.8030
两种方法比较得出:
使用MC方法要比常规方法精确度小
2.使用MC方法得出投掷三枚骰子点数之和大于10的概率
n=10000;
m=0;
fori=1:
x=randperm(6);
y=randperm(6);
z=randperm(6);
f=x
(1)+y
(1)+z
(1);
iff>
10
p=m/n
0.4986
3.使用MC方法求x^2*sin(x)cos(x)在0<
x<
2π内的最大值和最小值
x=unifrnd(0,2*pi,n,1);
y=x.^2.*sin(x).*cos(x);
max(y)
7.9544
min(y)
-15.3595
给该题函数作图
plot(x,y)
set(gca,'
ytick'
[-15,0,8])------在-15,0,8处分别给y画线
gridon-------------------------------------------设网格-
得图如下
②x取(0,2π)内均匀值
x=linspace(0,2*pi,n);
7.9544
作图如下
[-15,0,8])
gridon
得图为
4.求期权价格
在给出各个参数具体之后,到期日期权现金流可简化为
Max{0,50*exp(2+4ε)-52}
当50*exp(2+4ε)-52>
0时期权才有价值
ε服从标准正态分布
则计算程序为:
x=normrnd(0,1,n,1);
y=50*exp(2+4*x
(1))-52;
ify>
f=m/n
1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第二次 上机 作业