实验六1_蒙特卡罗方法.pptVIP

38/wlxt/ncourse/mathlab/web/shiyan/ * */10 * */10 ? ? ? ? 实验六:蒙特卡罗方法实验 面积、体积计算问题 冰淇淋锥的体积计算 思考题与练习题 二维填充图绘制方法 rand 产生一个0到 1之间均匀随机数 rand(m,n) 产生m×n个0到 1之间均匀随机数 X=rand(10000,1); hist(X) 蒙特卡罗方法——利用随机试验做近似计算 10000个随机数较均匀地分布在各个小区间上，随机变量X落入小区间的概率仅与小区间长度有关,而与小区间位置无关 unifrnd(a,b) 产生一个a到 b之间均匀随机数 X=unifrnd(2,10) Y=2+(10-2)*rand 例1 计算两条抛物线y = x2, x = y2所围面积. 在正方形[0,1]×[0,1]区域投入2000个均匀随机点 则随机点落入抛物线所围区域的概率为所求面 积与正方形面积之比 function S=area1(N) if nargin==0,N=2000;end X=rand(N,1);Y=rand(N,1); II=find(Y=sqrt(X)Y=X.^2); m=length(II);S=m/N; x1=0:0.01:1;x2=1:-0.01:0; y1=sqrt(x1);y2=x2.^2; fill([x1,x2],[y