第二次上机作业我自己留的

第二次上机答案1.1做三维图形
①解
>> n=100;
>>x1=linspace(0,1,n);
>> y1=linspace(1,2,n);
>> [x,y]=meshgrid(x1,y1);
>> z=cos(x)+sin(y);
>> mesh(x1,y1,z);
②解
>> n=100;
>>x1=linspace(0,1,n);
>> y1=linspace(1,2,n);
>> [x,y]=meshgrid(x1,y1);
>>mesh(cos(x)+sin(y));
得图如下
使用surf作图
>> n=100;
>>x1=linspace(0,1,n); >> y1=linspace(1,2,n); >> [x,y]=meshgrid(x1,y1); >> z=cos(x)+sin(y);
>>surf(z);
得图如下
使用waterfall作图>> n=100;
>>x1=linspace(0,1,n); >> y1=linspace(1,2,n); >> [x,y]=meshgrid(x1,y1); >> z=cos(x)+sin(y);
>>waterfall(z);
得图如下
使用plot3作图
>> n=100;
>>x1=linspace(0,1,n); >> y1=linspace(1,2,n); >> [x,y]=meshgrid(x1,y1); >> z=cos(x)+sin(y);
>> plot3(x1,y1,z);
得图如下
使用meshz作图
>> n=100;
>>x1=linspace(0,1,n); >> y1=linspace(1,2,n); >> [x,y]=meshgrid(x1,y1); >> z=cos(x)+sin(y);
>>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)
ans =
1.7979
③解
>>fun=cos(x)+sin(y);
>>dblquad(@(x,y)fun,0,1,1,2)
ans =
1.8444
PS：不知道答案为什么会不一样
④解
>>fun=inline('cos(x)+sin(y)','x,y')
fun =
Inline function:
fun(x,y) = cos(x)+sin(y)
>>jf=dblquad(fun,0,1,1,2)
这个是抄老师的，但是我实现不了，不知道是不是抄错了⑤解
>>jf=dblquad(inline('cos(x)+sin(y)','x,y'),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
for j=1:n
z1=cos(x(i))+sin(y(j));
if z(i,j)<z1
m=m+1;
end
end
end
>> p=m/n^2
p =
0.9015
>> f=2*p
f =
1.8030
两种方法比较得出：使用MC方法要比常规方法精确度小2．使用MC方法得出投掷三枚骰子点数之和大于10的概率>> n=10000;
m=0;
fori=1:n
x=randperm(6);
y=randperm(6);
z=randperm(6);
f=x(1)+y(1)+z(1);
if f>10
m=m+1;
end
end
>> p=m/n
p =
0.4986
3.使用MC方法求x^2*sin(x)cos(x)在0<x<2π内的最大值和最小值>> n=10000;
>> x=unifrnd(0,2*pi,n,1);
>> y=x.^2.*sin(x).*cos(x);
>> max(y)
ans =
7.9544
>>min(y)
ans =
-15.3595
给该题函数作图
>>plot(x,y)
>> set(gca,'ytick',[-15,0,8]) ------在-15,0,8处分别给y画线>> grid on -------------------------------------------设网格- 得图如下
②x取（0,2π）内均匀值>> n=10000;
>> x=linspace(0,2*pi,n); >> y=x.^2.*sin(x).*cos(x); >> max(y)
ans =
7.9544
>>min(y)
ans =
-15.3595
作图如下
>>plot(x,y)
>>set(gca,'ytick',[-15,0,8]) >>grid on
得图为
4.求期权价格
在给出各个参数具体之后，到期日期权现金流可简化为Max{0,50*exp（2+4ε）-52}
当50*exp（2+4ε）-52>0时期权才有价值
ε服从标准正态分布
则计算程序为：
>> n=10000;
>> x=normrnd(0,1,n,1);
>> m=0;
>>fori=1:n
y=50*exp(2+4*x(1))-52;
if y >0
m=m+1;
end
end
>> f=m/n
f =
1
不知道答案是什么，总感觉不太对。