大学数学实验之蒙特卡洛方法.pdfVIP

大学数学实验之蒙特卡洛方法 《数学实验》报告 班级： 序号： 姓名： 1．问题描述 I 、用蒙特卡罗方法计算以下函数在 区间上的积分，并改变随机点数目观 察对结果的影响。 （1） y=1/(1+x), 0=x=1; （2 ） y= (exp(3*x))*sin(2*x), 0=x=2 ; （3 ） y=(1+x^2)^0.5, 0=x=2; （4 ） y=(1/(2*pi)^0.5)*exp(-x( i)^2/2), 0=x=2; （5 ） y=exp(x(i)/2)*(sin(x (i)))^2, 0=x=2*pi; （6 ） f(x,y)=exp(-x^2-y^2) 0=x=pi, 0=y=sin(x); II 、用蒙特卡罗法求解全局最优化及 约束问题并通过图形做出评论， 求下 列函数的最大值。 （1） f(x)=(1-x.^2).*sin(3*x), -2*pi=x=2*pi; （2 ） max f(x)=x1*x2*x3,s.t.:-x1+2x2+ 2x3=0,x1+2x2+2x3=72,10=x 2=20,x1-x2=10; （3 ） f(x,y)=(X.^2+2*(Y.^2)+X. *Y).*exp(-X.^2-Y.^2), abs(x)1.5,abs(y)1.5; 2 ．问题分析与实验过程 I 、(1) 使用均值估计法 程序： function p=shell1(a,b,n) z=0; x=unifrnd(a,b,1,n); for i=1:n u=(x(i)+1)^(-1); z=z+u; end p=(b-a)*z/n; 运行结果： 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 结果分析：改变了四次随机点 数，结果都趋近于 0.69 ，说明积分值 约等于 0.69 ，但是点数越多，值越接 近。 I 、(2) 使用均值估计法 程序： function p=shell2(a,b,n) z=0; x=unifrnd(a,b,1,n); for i=1:n u=(exp(3*x(i)))*sin(2*x(i )); z=z+u; end p=(b-a)*z/n; 运行结果： 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