《第二次上机答案汇总.docx》由会员分享，可在线阅读，更多相关《第二次上机答案汇总.docx（8页珍藏版）》请在金锄头文库上搜索。
1、一题第一问首先在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的话，可以直接写原函数总结：这一问具体的操作
2、为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类似
3、，只是它做出的图是黑色的，所以取点时不要太密，一般要少用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.7979fun=(x
4、,y)cos(x)+sin(y)fun = (x,y)cos(x)+sin(y)q=dblquad(fun,0,1,1,2)q =1.7979f=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(inl
5、ine(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;for i=1:nfor j=1:nzt=cos(x(i)+sin(y(j);if z(i,j) p=frq/n2p = 0.8969 jf=2*pjf = 1.7938第二题n=10000;frq=0;for i=1:nx1=randperm(6);x2=randperm(6);x3=randperm(6);x=x1(1)+x2
6、(1)+x3(1);if x10frq=frq+1;endend p=frq/np =0.5036改变n的数字n=100000;frq=0;for i=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图形为第四题在已知股票价格和行权价格，无风险利率，标准差，到期天数后，到期日现金流可以简化为Max0,50exp(2+4)-52n=10000000;x=randn(n,1);frq=0;for i=1:nx1=x(i);f=max(0,50*exp(2+4*x1)-52);frq=frq+f;endp=frq/np = 9.8911e+005
