第二次上機答案匯總.docx

一題第一問首先在x,y的取值范圍內(nèi)取若干個點，用linspace命令，如對于x則是x1=linspace(0,1,200)括號內(nèi)前兩個數(shù)是x的取值范圍，最后一個數(shù)是分成的份數(shù)用meshgride命令將（x1,y1）從矢量轉(zhuǎn)化為x,y矩陣，具體的函數(shù)是x,y=meshgride(x1,y1)設(shè)一個函數(shù)z簡寫將要積分的函數(shù)，即z=cos（x）+sin(y)以便以后畫圖時方便，沒有這一步也是可以的最后一步就是利用mesh函數(shù)畫出三維圖形，具體命令為mesh(x1,y1,z)前面兩個是我們自己設(shè)的變量，順序可以換，最后一個是被積函數(shù)，如果前面沒有將它簡寫為z的話，可以直接寫原函數(shù)總結(jié)：這一問具體的操作為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完全一樣，只是在下方加上了一些線而構(gòu)成了矩形的方框，meshz截圖如下waterfall 和mesh功能完全一樣，只是會產(chǎn)生動態(tài)瀑布的效果 waterfall截圖 surf用法也和mesh類似，只是它做出的圖是黑色的，所以取點時不要太密，一般要少用100個，surf所做三維圖形截圖Ps:以上函數(shù)最后一行均可簡寫為函數(shù)名(z)不影響結(jié)果，如mesh(z)等價于mesh(x1,y1,z)plot用來畫三維曲線圖而不是曲面圖，它的功能就是將所有的x,y,z坐標表示的點連起來，plot3 做的三維曲面截圖 要想畫曲面圖應(yīng)該用以上幾個一題第二問dblquad 用于進行數(shù)值積分，他的用法是：設(shè)一個函數(shù)名q ，令 q=dblquad（x,y）被積函數(shù)，x下限，x上限，y下限，y上限）即可算出 q=dblquad(x,y)(cos(x)+sin(y),0,1,1,2)q =1.7979fun=(x,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(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;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(1)+x3(1);if x10frq=frq+1;endend p=frq/np =0.5036改變n的數(shù)字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圖形為第四題在已知股票價格和行權(quán)價格，無風(fēng)險利率，標準差，到期天數(shù)后，到期日現(xiàn)金流可以簡化為Max0,50exp(2+4)-52n=10000