概率统计随机数应用实验.pptxVIP

随机数应用实验随机数与统计直方图相遇问题与保险问题平面多边形填充图积分计算蒙特卡罗方法????

均匀分布随机数MATLAB产生均匀随机数方法:rand(m，n)产生m×n个0,1之间均匀随机数.随机数等可能落入区间[0，1]内长度相等子区间中。O1引例1.观察12个1—4之间整型随机数情况1+fix(4*rand(1,12))ans=413244214234引例2.观察1000个随机点分布情况P=rand(2,1000);x=P(1,:);y=P(2,:);plot(x,y,b.)

例2.观察1000个随机数在[0,0.5],[0.5,1]分布情况functionF=myrand(n)ifnargin==0,n=1000;endX=rand(1,n);Index=find(X0.5);f1=length(Index);F=[f1,n-f1];第一次实验:490510第二次实验:497503第三次实验:508492第四次实验:511489

统计直方图其中,data是需要处理的数据块,绘图原理:利用data中最小数和最大数构成一区间,将区间等分为n个小区间，统计落入每个小区间的数据量。以数据量为高度绘小矩形，形成直方图。如果省略参数n，MATLAB将n的默认值取为10。直方图也可以用于统计计算N=hist(data，n)计算结果N是n个数的一维数组，分别表示data中各个小区间的数据量。这种方式只计算而不绘图。直方图绘图方法:hist(data，n)

N5=19692010201819992004例5.1统计10000个均匀随机数在五个小区间的分布。data=rand(10000,1);figure(1),hist(data,5)N5=hist(data,5)figure(2),bar(N5,r)即观察10000个随机数在[0,0.2],[0.2,0.4],[0.4,0.6],[0.6,0.8],[0.8,1]分布情况

例3.观察1000个平面随机点在单位正方形内的分布情况P(x,y)的坐标均是[0,1]上均匀随机数,functionF=myrand2(n)ifnargin==0,n=1000;endP=rand(n,2);x=P(:,1);y=P(:,2);I1=find(x0.5y0.5);I2=find(x=0.5y0.5);I3=find(x=0.5y=0.5);I4=find(x0.5y=0.5);F(1,1)=length(I1);F(2,1)=length(I2);F(2,2)=length(I3);F(1,2)=length(I4);bar3(F,c)ans=244233259264

引例3.实验观察10个1—4之间随机数情况1+3*rand(12,1)一般区间[a，b]上的均匀随机数产生方法R=a+(b-a)*rand3.57153.82602.29641.57883.25882.78131.52963.66803.05073.2386第1次?3.17601.95103.37111.15782.90263.36701.53391.06431.53491.5648第2次?2.71973.91602.87901.88923.54161.77361.42281.84191.16861.8978第3次?

01产生“0”和“1”随机数:fix(2*rand)产生“1”到“100”整型随机数:1+fix(100*rand)rand(m,n)产生区间(0,1)上均匀分布的m×n个随机数.0203产生整型随机数方法随机数注记

均匀分布随机变量X~U(0,24),Y~U(0,24)如果甲船到达码头后停留2小时,乙船到达码头后停留1小时.问两船相遇的概率有多大？例5.2相遇问题:甲、乙两船在24小时内独立地随机到达码头.设两船到达码头时刻分别为X和YS1S2XYO2424

相遇问题的统计试验functionF=shipmeet(N)ifnargin==0,N=2000;endP=24*rand(2,N);X