第二次上機(jī)答案匯總.docx

一題第一問(wèn)首先在x,y的取值范圍內(nèi)取若干個(gè)點(diǎn)，用linspace命令，如對(duì)于x則是x1=linspace(0,1,200)括號(hào)內(nèi)前兩個(gè)數(shù)是x的取值范圍，最后一個(gè)數(shù)是分成的份數(shù)用meshgride命令將（x1,y1）從矢量轉(zhuǎn)化為x,y矩陣，具體的函數(shù)是x,y=meshgride(x1,y1)設(shè)一個(gè)函數(shù)z簡(jiǎn)寫(xiě)將要積分的函數(shù)，即z=cos（x）+sin(y)以便以后畫(huà)圖時(shí)方便，沒(méi)有這一步也是可以的最后一步就是利用mesh函數(shù)畫(huà)出三維圖形，具體命令為mesh(x1,y1,z)前面兩個(gè)是我們自己設(shè)的變量，順序可以換，最后一個(gè)是被積函數(shù)，如果前面沒(méi)有將它簡(jiǎn)寫(xiě)為z的話(huà)，可以直接寫(xiě)原函數(shù)總結(jié)：這一問(wèn)具體的操作為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完全一樣，只是在下方加上了一些線(xiàn)而構(gòu)成了矩形的方框，meshz截圖如下waterfall 和mesh功能完全一樣，只是會(huì)產(chǎn)生動(dòng)態(tài)瀑布的效果 waterfall截圖 surf用法也和mesh類(lèi)似，只是它做出的圖是黑色的，所以取點(diǎn)時(shí)不要太密，一般要少用100個(gè)，surf所做三維圖形截圖Ps:以上函數(shù)最后一行均可簡(jiǎn)寫(xiě)為函數(shù)名(z)不影響結(jié)果，如mesh(z)等價(jià)于mesh(x1,y1,z)plot用來(lái)畫(huà)三維曲線(xiàn)圖而不是曲面圖，它的功能就是將所有的x,y,z坐標(biāo)表示的點(diǎn)連起來(lái)，plot3 做的三維曲面截圖 要想畫(huà)曲面圖應(yīng)該用以上幾個(gè)一題第二問(wèn)dblquad 用于進(jìn)行數(shù)值積分，他的用法是：設(shè)一個(gè)函數(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)其中老師講的和不知道是不是版本不兼容問(wèn)題，總是做不出來(lái)一題第三問(wèn) 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對(duì)三題中所描述的曲線(xiàn)進(jìn)行繪圖如下x2=linspace(0,2*pi,1000);y2=x2.2.*sin(x2).*cos(x2);plot(x2,y2)set(gca,ytick,-15,0,8)grid on圖形為第四題在已知股票價(jià)格和行權(quán)價(jià)格，無(wú)風(fēng)險(xiǎn)利率，標(biāo)準(zhǔn)差，到期天數(shù)后，到期日現(xiàn)金流可以簡(jiǎn)化為Max0,50exp(2+4)-52n=10000