1用蒙特卡洛方法估计积分,和的值,并将估计值与真值进行比较。=1;fort=1:1:50n=1000;a=(pi)/2*rand(1,n);fori=1:1:nb(i)=(pi)/2*a(i)*sin(a(i));endc(t)=mean(b);endcmean((c-true).*(c-true))结论当N为200时计算的值分别为c======-004再用莱姆大为1的指数分布估计clearclctrue=1;fort=1:1:50n=1000;b=exprnd(1,1,n);fori=1:1:nif(b(i)>0&b(i)<pi/2)c(i)=b(i)*sin(b(i))*exp(b(i));elsec(i)=0;endendmean(c)a(t)=mean(b);endamean((a-true).*(a-true))当N为1000时计算的值分别为a=
用蒙特卡洛计算积分 来自淘豆网www.taodocs.com转载请标明出处.