第二次上机答案汇总.docx

一题第一问首先在x,y的取值范围内取若干个点，用linspace命令，如对于x则是x1=linspace(0,1,200)括号内前两个数是x的取值范围，最后一个数是分成的份数用meshgride命令将（x1,y1）从矢量转化为[x,y]矩阵，具体的函数是[x,y]=meshgride(x1,y1)设一个函数z简写将要积分的函数，即z=cos（x）+sin(y)以便以后画图时方便，没有这一步也是可以的最后一步就是利用mesh函数画出三维图形，具体命令为mesh(x1,y1,z)前面两个是我们自己设的变量，顺序可以换，最后一个是被积函数，如果前面没有将它简写为z的话，可以直接写原函数总结：这一问具体的操作为x1=linspace(0,1,200);y1=linspace(1,2,200); [x,y]=meshgrid(x1,y1); z=cos(x)+sin(y); mesh(x1,y1,z)或者x1=linspace(0,1,200);y1=linspace(1,2,200);[x,y]=meshgrid(x1,y1); mesh(x1,y1,cos(x)+sin(y))作出图形如下meshz的作用和mesh完全一样，只是在下方加上了一些线而构成了矩形的方框，meshz截图如下waterfall 和mesh功能完全一样，只是会产生动态瀑布的效果 waterfall截图surf用法也和mesh类似，只是它做出的图是黑色的，所以取点时不要太密，一般要少用100个，surf所做三维图形截图Ps:以上函数最后一行均可简写为函数名(z)不影响结果，如mesh(z)等价于mesh(x1,y1,z)plot用来画三维曲线图而不是曲面图，它的功能就是将所有的x,y,z坐标表示的点连起来，plot3 做的三维曲面截图要想画曲面图应该用以上几个一题第二问dblquad用于进行数值积分，他的用法是：设一个函数名q ，令 q=dblquad（@（x,y）被积函数，x下限，x上限，y下限，y上限）即可算出q=dblquad(@(x,y)(cos(x)+sin(y)),0,1,1,2)q =1.7979②fun=@(x,y)cos(x)+sin(y)fun = @(x,y)cos(x)+sin(y)q=dblquad(fun,0,1,1,2)q =1.7979③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)④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)其中老师讲的④和⑤不知道是不是版本不兼容问题，总是做不出来一题第三问 n=1000;x=unifrnd(0,1,n,1);y=unifrnd(1,2,n,1);z=unifrnd(0,2,n,n);frq=0;fori=1:nfor j=1:nzt=cos(x(i))+sin(y(j));if z(i,j)ztfrq=frq+1;endendend p=frq/n^2p = 0.8969 jf=2*pjf = 1.7938第二题n=10000;frq=0;fori=1:nx1=randperm(6);x2=randperm(6);x3=randperm(6);x=x1(1)+x2(1)+x3(1);if x10frq=frq+1;endend p=frq/np =0.5036改变n的数字n=100000;frq=0;fori=1:nx1=randperm(6);x2=randperm(6);x3=randperm(6);x=x1(1)+x2(1)+x3(1);if x10frq=frq+1;endend p=frq/np = 0.5014第三题答案x=unifrnd(0,2*pi,1000,1);y=x.^2.*sin(x).*cos(x);max(y)ans = 7.9543 min(y)ans = -15.3593对三题中所描述的曲线进行绘图如下x2=linspace(0,2*pi,1000);y2=x2.^2.*sin(x2).*cos(x2);plot(x2,y2)set(gca,ytick,[-15,0,8])grid on图形为第四题在已知股票价格和行权价格，无风险利率，标准差，到期天数后，到期日现金流可以简化为Max{0,50exp(2+