最新课件MATLAB数值计算第9章随机数.docx

《最新课件MATLAB数值计算第9章随机数.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、MATLAB数值计算(读书日记及程序编写)第九章随机数2第九章随机数如果电脑不能访问外部设备，是不能出现随机行为的，所以只能是伪随机数。#随机数生成过程（相乘取模算法）有三个整数参数a,c,m,以及一个称为种子的初始值x0,其中操作modm的意思是除以m后取余数当a=13,c=0,m=31,x0=1时…这个序列以m-1为一个周期Matlab均匀随机数函数rand也是一个乘法同余生成器采用的参数是：C=0#产生二维随机的点n=20000s=0.02x=rand(n,1)-0.5;y=rand(n,1)-0.5

2、;h=plot(x,y,'.');可以通过随机数来点名，rand(1)产生一个[0,1]间的一个随机数，乘以学生总数即可。rand(1)*76蒙特卡罗方法取正方形，边长为1,以1为半径，做一个1/4的圆，一个人随机向正方形中投掷飞镖，总共投了n次，进行计数，在圆内的就是Sqrt(x^2+y^2)≤1，假如圆点数为m,圆外的就是sqrt(x^2+y^2)>1的，点数就为n-m，总的点数是n。投中的点数和其面积是存在比例关系的：1/4圆的面积为1/4*π*1^2=π/4,正方形的面积为1即(π/4)/1=m/n

3、,即π=4*m/n#通过随机数求圆周率πk=1000000;n=0;fori=1:kRx=rand(1);Ry=rand(1);ifsqrt(Rx^2+Ry^2)<1n=n+1;endendnPI=4*n/k#randguihelprandguiRANDGUIMonteCarlocomputationofpi.Generaterandompointsinacubeandcounttheportionthatarealsointheinscribedsphere.Theratioofthevolumeofth

4、espheretothevolumeofthecubeispi/6.RANDGUIwithnoargumentsorRANDGUI('rand')usesMATLAB'sbuilt-inrandomnumbergenerator.RANDGUI('randtx')usesourtextbookversionofthebuilt-ingenerator.RANDGUI('randmcg')andRANDGUI('randssp')usedifferentLehmercongruentialgenerators

5、,onewithgoodparametersandonewiththeparametersusedyearsagobyIBM's"RANDU"function.得到的显示就是π的值。