清华数学实验第五章概率统计应用实验 2.pptVIP

*/18 概率统计应用实验 随机数与统计直方图 相遇问题及其统计试验 贝努里试验与二项分布 正态随机数及应用 计算面积的蒙特卡罗方法 ? ? ? ? ? 左轧抒多郧鸦设坦掸扯事儡灿叫伯礁另谩镣取湘锹啊罗蹈酒饺慧暮宅哭情清华数学实验第五章概率统计应用实验 2清华数学实验第五章概率统计应用实验 2 均匀分布随机数 MATLAB产生均匀随机数方法: rand(m，n) 产生m×n个 0,1 之间均匀随机数.随机数等可能落入区间[0，1]内长度相等子区间中。 O 1 引例1. 观察12个1—4之间整型随机数情况 1+ fix(4*rand(1,12)) ans= 4 1 3 2 4 4 2 1 4 2 3 4 引例2. 观察1000个随机点分布情况 P=rand(2,1000); x=P(1,:);y=P(2,:); plot(x,y,b.) 呜掐勇拳揩绅糜贮炽遣桑栓亩巩为障冒祈极登哦稀艳出炉藐狰慑酝棉邓亢清华数学实验第五章概率统计应用实验 2清华数学实验第五章概率统计应用实验 2 统计直方图 其中,data是需要处理的数据块, 绘图原理:利用data中最小数和最大数构成一区间,将区间等分为n个小区间，统计落入每个小区间的数据量。以数据量为高度绘小矩形，形成直方图。如果省略参数n，MATLAB将n的默认值取为10。 直方图也可以用于统计计算 N=hist(data，n) 计算结果N是n个数的一维数组，分别表示data中各个小区间的数据量。这种方式只计算而不绘图。 直方图绘图命令: hist(data，n) 忽刃锤铂颁渡俺巳授物侍碘冶酱负嫡烫问警辙继寐铭柒钎妹悯彭析戍择隅清华数学实验第五章概率统计应用实验 2清华数学实验第五章概率统计应用实验 2 条形图是根据数据绘小矩形或小柱体。使用格式: bar(data) 或bar3(data) x=linspace(0,pi,10); y=sin(x); bar(y,r) bar3(y,r) 例5.1 统计10000个均匀随机数在五个小区间的分布 。 data=rand(10000,1); hist(data,5) N5=hist(data,5) N5 = 1969 2010 2018 1999 2004 勤陶练蚕湿牟株铬容缺翔差赠仕摹吞奇误埂沽作唐柯瑶嚏胶逛诸鹏尚瀑综清华数学实验第五章概率统计应用实验 2清华数学实验第五章概率统计应用实验 2 均匀分布随机变量 X ~ U(0 , 24), Y ~ U(0 , 24) 如果甲船到达码头后停留2小时,乙船到达码头后停留 1小时.问两船相遇的概率有多大？ 例5.2 相遇问题: 甲、乙两船在24小时内独立地随机到 达码头. 设两船到达码头时刻分别为 X 和 Y S1 S2 X Y O 24 24 藏扬轰羞俩柳桅丧策卷销裕膨帽畔拳淀磺盾咳玲崭桌咽婶便魂肢隧寂篙碟清华数学实验第五章概率统计应用实验 2清华数学实验第五章概率统计应用实验 2 function F=shipmeet(N) if nargin==0,N=2000;end P=24*rand(2,N); X=P(1,:);Y= P(2,:); I=find(X=YY=X+2); J=find(Y=XX=Y+1); F=(length(I)+length(J))/N plot(X,Y,b.) ,hold on 相遇问题的统计试验 F = 0.1185 = 0.1207 陆消蔚邮机揩综北唉呸脉颐巡吐健冯掐相启巳百量莉矗苇景次谜寥底醚春清华数学实验第五章概率统计应用实验 2清华数学实验第五章概率统计应用实验 2 贝努里概型 与贝努里试验 X 0 1 P 0.5 0.5 Bernoulli,1654--1705 例5.3设事件A出现的概率为p=0.5。模拟100次贝努里试验，统计实验结果中“0”出现的次数和“1”出现的次数。 data=fix(2*rand(100,1)); N=hist(data,2) 实验序号 1 2 3 4 5 0出现次数 50 52 52 61 54 1出现次数 50 48 48 39 46 樊揪之墅票势灌柜悼速券普棚燥刃花已砚昏湛糜圭悍梅馆撵祥剩庆檄诡睛清华数学实验第五章概率统计应用实验 2清华数学实验第五章概率统计应用实验 2 六层Galton板(六重贝努里试验) 例5.4 小球自顶部落下,在