大學(xué)數(shù)學(xué)實驗之蒙特卡洛方法

《數(shù)學(xué)實驗》報告班級：序號：姓名：問題描述=1\*ROMANI、用蒙特卡羅方法計算以下函數(shù)在區(qū)間上的積分，并改變隨機(jī)點數(shù)目觀察對結(jié)果的影響。y=1/(1+x),0=<x=<1;y=(exp(3*x))*sin(2*x),0=<x=<2;y=(1+x^2)^0.5,0=<x=<2;y=(1/(2*pi)^0.5)*exp(-x(i)^2/2),0=<x=<2;y=exp(x(i)/2)*(sin(x(i)))^2,0=<x=<2*pi;f(x,y)=exp(-x^2-y^2)0=<x=<pi,0=<y=<sin(x);=2\*ROMANII、用蒙特卡羅法求解全局最優(yōu)化及約束問題并通過圖形做出評論，求下列函數(shù)的最大值。f(x)=(1-x.^2).*sin(3*x),-2*pi=<x=<2*pi;maxf(x)=x1*x2*x3,s.t.:-x1+2x2+2x3>=0,x1+2x2+2x3<=72,10<=x2<=20,x1-x2=10;f(x,y)=(X.^2+2*(Y.^2)+X.*Y).*exp(-X.^2-Y.^2),abs(x)<1.5,abs(y)<1.5;問題分析與實驗過程=1\*ROMANI、(1)使用均值估計法程序：functionp=shell1(a,b,n)z=0;x=unifrnd(a,b,1,n);fori=1:nu=(x(i)+1)^(-1);z=z+u;endp=(b-a)*z/n;運行結(jié)果：p=shell1(0,1,1000)p=0.6975>>p=shell1(0,1,10000)p=0.6922>>p=shell1(0,1,100)p=0.7001>>p=shell1(0,1,500)p=0.6890結(jié)果分析：改變了四次隨機(jī)點數(shù)，結(jié)果都趨近于0.69，說明積分值約等于0.69，但是點數(shù)越多，值越接近。=1\*ROMANI、(2)使用均值估計法程序：functionp=shell2(a,b,n)z=0;x=unifrnd(a,b,1,n);fori=1:nu=(exp(3*x(i)))*sin(2*x(i));z=z+u;endp=(b-a)*z/n;運行結(jié)果：>>p=shell2(0,2,1000)p=-24.4911>>p=shell2(0,2,100)p=-43.8720>>p=shell2(0,2,10000)p=-30.8699>>p=shell2(0,2,500)p=-23.2955>>p=shell2(0,2,100000)p=functionf81(n)x=unifrnd(-2*pi,2*pi,1,n);y=(1-x.^2).*sin(3*x);max(y)x=-2*pi:0.001:2*pi;y=(1-x.^2).*sin(3*x);plot(x,y)xlabel('x');ylabel('y');運行結(jié)果：>>f81(1000)ans=32.3293>>f81(10000)ans=32.4002>>f81(100000)ans=32.4006做出函數(shù)的圖像，并且標(biāo)出最高點的值結(jié)果分析：可以看到，蒙特卡羅法求出的最大值接近于32.4，而從圖中可以看出最大值是32.33，求出的結(jié)果比較符合。=2\*ROMANII、(2)使用均值估計法分析：由于x1=x2+10,所以可以消元，使其變?yōu)閮蓚€自變量x2和x3。x2,x3在它們被允許的范圍內(nèi)生成多個隨機(jī)的數(shù)值，利用max函數(shù)可以近似地求出結(jié)果。然后做出圖像，進(jìn)行結(jié)果的比較。程序：functionf82(n)x2=unifrnd(10,20,1,n);x1=10+x2;x3=unifrnd(-10,20,1,n);fori=1:nif-x1(i)+2*x2(i)+2*x3(i)>=0ifx1(i)+2*x2(i)+2*x3(i)<=72y(i)=(x1(i))*(x2(i))*(x3(i));endendendmax(y)x2=10:0.1:20;x3=-5:21/100:16;[X,Y]=meshgrid(x2,x3);err1=X+2*Y<10;err2=3*X+2*Y>62;X(err1)=nan;Y(err2)=nan;Z=X.*Y.*(X+10);surf(X,Y,Z)運行結(jié)果：>>f82(1000)ans=3.3889e+03>>f82(10000)ans=3.4357e+03>>f82(100)ans=3.3726e+03>>f82(100000)ans=3.4441e+03結(jié)果分析：可以看到，蒙特卡羅法求出的最大值接近于3400，而從圖中可以看出最大值是3437，求出的結(jié)果比較符合。=2\*ROMANII、(3)使用蒙特卡羅法分析：x,y在它們被允許的范圍內(nèi)生成多個隨機(jī)的數(shù)值，利用max函數(shù)可以近似地求出結(jié)果。然后做出圖像，進(jìn)行結(jié)果的比較。程序：functionf83(n)x=unifrnd(-1.5,1.5,1,n);y=unifrnd(-1.5,1.5,1,n);z=(x.^2+2*(y.^2)+x.*y).*exp(-x.^2-y.^2);max(z)x=-1.5:0.1:1.5;y=-1.5:0.1:1.5;[X,Y]=meshgrid(x,y);Z=(X.^2+2*(Y.^2)+X.*Y).*exp(-X.^2-Y.^2);surf(X,Y,Z)運行結(jié)果：>>f83(1000)ans=0.8105>>f83(10000)ans=0.8117作出函數(shù)圖，并且標(biāo)出最大值結(jié)果分析：可以看到，蒙特卡羅法求出的最大值接近于0.81，而從圖中可以看出最大值是0.8025，求出的結(jié)果比較符合。3.實驗總結(jié)和實驗感悟這次蒙特