- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
华中科技大学《概率论与数理统计》课程概率统计实验报告(综合性实验.PDF
华中科技大学《概率论与数理统计》课程
概率统计实验报告(综合性实验)
高尔顿订板实验及其应用
班级 电气0707 电气0707 电气07 中英班
学号 U200712024 U200712174 U200712200
姓名 张锐 王鑫 杨光
Email zhangruitiancai@163.com
自我评分 7
1、 问题概述和分析
Galton 钉板运用
高尔顿 GALTON, Francis 1822.2 .16—1911.1.17。英国的遗传学家。E.达尔文的孙子, C.达
尔文的表弟。 1882 年创建人体测量研究室,1904 年建立优生学研究所。1907 年该所移交给伦敦大学,
与此同时高尔顿( Galton )担任了第一代教授。最早在遗传学研究中采用数学方法。是生物统计学、优生
学和人类遗传学的创始人。而 Galton 钉板是他数学成就中最有名,最直接,最被大家广为认识
的实验。而我们小组想通过这个实验加深对MATLAB 的应用,以及运用这个程序能做一些
比如小抽奖的概率论运用。
主要结合本实验分析正态分布的图形。以及GALTON 钉板实验。
最主要对Matlab 的较为深入的学习,以掌握Matlab 语言为目的而开展的中型实验。
2、实验设计总体思路*
2.1、引论*
这是一个给我们团队学习MATLAB 实验程序的机会。这个程序在大学理工科中广泛应用,
我们希望更早得接触它,了解他,从而为自己以后的数学,专业课程学习与分析中抢得先机。
同时,也增加我们团队的默契和谐了寝室氛围。同时我们都很尊敬伟大的遗传学家,数学家。
高尔顿,我们希望通过他给我们留下的伟大财富,高尔顿钉板来加深我们对 MATLAB,对
概率论的学习。
2.2、实验主题部分*
2.2.1、实验设计思路
1、理论分析
做一个动画模拟高尔顿钉板实验。进而得到其的正态分布区间图,最后同化到一个小抽奖实
验,得出抽奖期望等结果
2、实现方法
首先用 MATLAB 做一个关于高尔顿钉板的实验动画从而得到近似的正态分布区间图。
(MATLAB 中可以做出动态模型从而得到高尔顿钉板动画)
然后将最外部赋值最高,越靠近中间赋值越低,得出近似抽奖模型。(小抽奖活动很多都是
依据次来得到的,而所得分布图近似于正态分布,从而使两边的奖金最高,让中间的奖金最
低,从而得到抽奖模型)
最后通过MATLAB 计算出得奖期望从而得到得奖情况。(根据MATLAB 中得到的期望设计
最好的抽奖方案)
2.2.2、实验结果及分析
近似于正态分布的高尔顿钉板实验。
这可以近似于抽奖实验中得奖概率分布图。
得出的得奖期望:
2.2.3、程序及其说明*
得到动画的函数:
clear
clf
m=100;n=5;y0=2; %设置参数。
ballnum=zeros(1,n+1);
p=0.5;q=1-p;
for i=n+1:-1:1 %创建钉子的坐标x,y 。
x(i,1)=0.5*(n-i+1);y(i,1)=(n-i+1)+y0;
for j=2:i
x(i,j)=x(i,1)+(j-1)*1;y(i,j)=y(i,1);
end
end
mm=moviein(m); % 动画开始,模拟小球下落路径。
for i=1:m
s=rand(1,n); %产生n 个随机数。
xi=x(1,1);yi=y(1,1);k=1;l=1; % 小球遇到第一个钉子。
for j=1:n
plot(x(1:n,:),y(1:n,:),o,x(n+1,:),y(n+1,:),.-), % 画钉子的位置 。
axis([-2 n+2 0 y0+n+1]),hold on
k=k+1; % 小球下落一格 。
if s(j)p
l=l+0; %小球左移 。
文档评论(0)