第5讲蒙特卡洛方法的应用(原稿).ppt

熊宇乐 y=zeros(1,1000); a=binornd(1,0.4,1,1000);b=binornd(1,0.4,1,1000); c=0;d=0; for i=1:1000 c=c+a(1,i).*b(1,i); y(i)=c/i; end y=y; num=1:1000; plot(num,y) liti2(1,0.4,100) liti2(1,0.4,10000) function proguji=liti5(mm) %mm 是随机实验次数 frq=0; randnum1=unifrnd(0,60,mm,1); randnum2=unifrnd(0,60,mm,1); randnum=randnum1-randnum2; proguji=0; for ii=1:mm if abs(randnum(ii,1))=20 frq=frq+1; end end proguji=frq/mm 圆的面积 单位圆的面积等于 ? 计算第一象限内的单位圆的面积，方法为：把它分成n个窄的曲边梯形，计算S大，S小，其中n可以为1000, 10000, ... 无穷级数法 arctg x =x-x3/3+x5/5-x7/7+x9/9-... ? /4 = arctg 1 =1-1/3+1/5-1/7+1/9-… 收敛太慢！ |x| 应当比