高尔顿钉板
- 格式:doc
- 大小:225.00 KB
- 文档页数:9
《概率论与数理统计》的课堂思政设计—以大数定律为例摘要:本文以大数定律教学内容为基础,以课程思政为导向,挖掘概率论与数理统计教学内容中蕴含的思政元素,把思政元素融入教学内容、教学方法和教学活动。
通过引入历史两个著名试验,追溯大数定律的发展和演变,剖析大数定律的内涵和意义三个方面阐述课程思政教学的有效运行,结合具体教学内容为课程思政的有效运行提供切实可行的方法。
关键词:大数定律、思政元素一、引言概率论与数理统计是研究与随机现象相关的数量规律的学科,也是高等学校理工科专业的通识必修课之一。
概率论与数理统计应用广泛,几乎遍及学术研究和日常生活的方方面面。
通过学习,学员可以提升运用概率论的思想观察、处理随机事件的能力。
大数定律在概率论与数理统计课程的教学内容中具有承上启下的重要意义,既是前面概率论内容的一个补充,又为数理统计提供了理论基础。
大数定律是概率论的基本理论,在理论研究和应用中起着重要的作用。
大数定律是概率论中讨论随机变量序列的算术平均值向常数收敛的定理。
从数学上严格地解释了频率稳定于概率,平均值稳定于数学期望。
本文以大数定律的教学内容为基础,深挖案例中蕴含的“思政元素”及所承载的思想政治教育功能,将思政元素有机的融入课堂教学,把“知识传授”与“价值引领”有机统一起来,做到立德树人,培养具有社会主义核心价值观的有用人才。
二、课堂思政设计(一)通过引入试验,透过实验现象看本质高尔顿钉板(Galton board),是弗朗西斯高尔顿以验证中心极限定理的试验。
从漏斗形上口掉落的小球会遇上一系列排列成三角形的“钉子”。
每当小球从正上方下落到一个“钉子”上时,它总是会有50%的概率跑到左边,有50%的概率跑到右边。
在经过数次这样随机的“左右选择”之后,小球掉落到下方的格子中。
如图1所示。
图1 高尔顿钉板试验引入高尔顿钉板试验,可以从直观上看到无数的随机因素共同作用的结果即每一个因素或多或少都起到一点作用,但都没有起到很大的甚至决定性的作用也就是说每种因素的微小差异对总的影响作用不是很大,最终综合在一起就形成了正态分布。
二、高尔顿钉板的理论解释及仿真
高尔顿钉板中(见教材第一章图1—2)每一黑点表示钉在板上的一颗钉子,它们彼此的距离均相等,上一层的每一颗的水平位置恰好位于下一层的两颗中间.从入口处放进一个直径略小于两颗钉子之间的距离的小圆玻璃球,当小圆球向下降落过程中,碰到钉子后均以2
1的概率向左或右滚下,于是又碰到下一层钉子.如此继续下去,直到滚到底板的一个格子内为止.把许许多多同样大小的小球不断从入口处放下,只要球的数目相当大,它们在底板将堆成近似于正态),0(n N 的密度函数图形,其中n 为钉子层数。
解释如下:
令k X 表示某个小球在第k 次碰钉子后向左或向右落下这一随机现象联系的
随机变量,k X 的概率分布为(...,2,1=k )k X 1-1
k
p 2121令∑==n i k n X Y 1,其中
)...21(n k X k ,,=相互独立,则n Y 表示这个小球第n 次
碰钉后的位置。
由中心极限定理,有⎰∞--∞→=⎭
⎬⎫⎩⎨⎧<x t n n dt e x n Y P 2221lim π故n Y ~),0(n N 。
下面进行计算机仿真。
%“高尔顿钉板问题”模拟程序
function y=f2
n=input('请输入球数n=');
m=input('请输入层数m=');
x(n)=0;yy(m+1)=0;
for i=1:n
for j=1:m
if rand>0.5
x(i)=x(i)+1;
else
x(i)=x(i)-1;
end
end。
高尔顿钉板试验模拟(程序)...这是我2005年12的课程设计中程序的核心部分,写完后自己非常得意,等着老师表扬。
等啊等,等待现在也没等到:em16:现将它献给大家...(若有版权那遵守BSD吧)注1:程序以前是用Matlab写的现用Java重写注2:原程序中galton返回值为int[] grid 、没有“输出结果”部分public void galton(int sumOfGrid, int sumOfBall){int[] grid = new int[sumOfGrid];int number = 0; //一个小球从顶端落下过程中向右偏移的总次数int rand ; //随机数,取值范围为{0,1},为0、为1的概率相等for( int counter_ball = 1; counter_ball <= sumOfBall; counter_ball++ ){//<核心>// (sumOfGrid - 1)为钉板的层数for( int times = 1; times <= ( sumOfGrid - 1 ); times++ ){rand = (int)( Math.random()*2 );number += rand;}grid[number]++;number = 0;//</核心>}//输出结果System.out.println( "小球的总数为"+sumOfBall+"\t格子的个数为"+sumOfGrid );for( int index = 0; index < grid.length; index++ )System.out.println( (index+1)+"号格子中的小球数为:\t"+grid[index] );}}//end of metod galton补充:(谢谢2楼提醒:-D )高尔顿钉板试验:自板上端放入一小球, 任其自由落下.在下落过程中, 当小球碰到钉子时, 从左边落下与从右边落下的机会相等.碰到下一排钉子也是如此.自板上端放入n(n自行输入)个小球, 观察小球落下后呈现曲线并统计小球落入各个格子的频率.高尔顿钉板试验可见《概率论》(复旦大学李贤平)当小球数量少时分布无明显特征,当小球数量多时(>100)分布近似正态分布。
2024年秋季普通高中11月份高三年级阶段性联考数学本试卷共4页,19题.全卷满分150分.考试用时120分钟.★祝考试顺利★注意事项:1.答题前,先将自己的姓名、准考证号、考场号、座位号填写在试卷和答题卡上,并将准考证号条形码粘贴在答题卡上的指定位置.2选择题的作答:每小题选出答案后,用2B 铅笔把答题卡上对应题目的答案标号涂黑.写在试卷、草稿纸和答题卡上的非答题区域均无效.3.非选择题的作答:用黑色签字笔直接答在答题卡上对应的答题区域内.写在试卷、草稿纸和答题卡上的非答题区域均无效.4.考试结束后,请将答题卡上交.一、选择题:本题共8小题,每小题5分,共40分.在每小题给出的四个选项中,只有一项是符合题目要求的.1.在复平面内,复数对应的点位于( )A.第一象限 B.第二象限C.第三象限 D.第四象限2.已知,则的值为( )A.B. C.D.3.已知,且,则与的夹角为( )A.B. C. D.4.已知曲线在点处的切线在轴上的截距为,则的值为( )A.1B.0C.D.5.暑假期间某校5名学生计划去黄冈旅游,体验黄冈的风俗与文化.现有黄梅东山问梅村、罗田天堂寨、黄州的东坡赤壁三个景区可供选择若每名学生只去一个景区,且恰有2人前往黄梅东山问梅村,则不同的游览方案种数为( )A.40B.90C.80D.16011i+π1cos 33α⎛⎫-=- ⎪⎝⎭πsin 6α⎛⎫+ ⎪⎝⎭1313-(),2a b == ()2a a b ⊥+ a bπ32π33π45π6ln ay x x=+()1,a y 3-a 1-2-6.已知函数的最小正周期为,将的图象向右平移个单位后得到函数的图象,若为偶函数,则正实数的最小值为( )A.B. C. D.7英国生物统计学家高尔顿设计了高尔顿钉板来研究随机现象.如图是一个高尔顿钉板的设计图,每一黑点表示钉在板上的一颗钉子,它们彼此的距离均相等,上一层的每一颗钉子恰好位于下一层两颗打子的正中间,小球每次下落,将随机的向两边等概率的下落.数学课堂上,老师向学生们介绍了高尔顿钉板放学后,爱动脑的小明设计了一个不一样的“高尔顿钉板”,它使小球在从钉板上一层的两颗钉子之间落下后砸到下一层的钉子上时,向左下落的概率为向右下落的概率的2倍.当有大量的小球依次滚下时,最终都落入钉板下面的5个不同位置.若一个小球从正上方落下,经过5层钉板最终落到4号位置的概率是()A.B. C. D.8.是定义在上的函数,为的导函数,若方程在上至少有3个不同的解,则称为上的“波浪函数”.已知定义在上的函数为“波浪函数”,则实数的取值范围是()A. B.C. D.二、多选题:本题共3小题,每小题6分,共18分.在每小题给出的四个选项中,有多项是符合题目要求.全部选对的得6分,部分选对得部分分,有错选的得0分.9.下列结论中正确的有( )A.已知,若,则;B.某学生8次考试的数学成绩分别为:101、108、109、120、132、135、141、141,则这8次数学成绩的第75百分位数为135;C.已知的平均值为8,则的平均值为7;D.已知为两个随机事件,若,则.()()cos 0f x x x ωωω=->π()f x ϕ()g x ()g x ϕπ12π6π32π3881168124813281()f x [],a b ()f x '()f x ()()f x f x ='[],a b ()f x [],a b []4,3-()3228f x x x mx =+++m 5675m -<- (56)45m -<- (56)45m -< (74)m -<-…()24,X N σ~()50.1P X =…()340.4P X =……128,,,,11,13x x x 128,,,x x x A B 、()()()0.4,0.3,0.2P A P B P AB ===∣()0.15P B A =∣10.已知正实数满足,下列结论中正确的是()A.的最大值是B.的最小值是C.的最小值是3D.的最小值为11.高斯被誉为“数学王子”,是世界上伟大数学家.用他名字定义的函数(表示不超过的最大整数)称为高斯函数.已知正项数列的前项和为,且,令,则下列结论正确的有( )A.B.C.D.三、填空题:本题共3小题,每小题5分,共15分.12.已知函数,则__________.13.已知的角的对边分别为,且,若,则__________.14.已知函数在区间上存在零点,则的取值范围为__________.四、解答题:本题共5小题,共77分.解答应写出文字说明、证明过程或演算步骤.15.(本题满分13分)已知,函数.(1)求的单调递减区间;(2)在中,若,求和长.16.(本题满分15分)已知是公差不为0的等差数列,,且成等比数列,数列满足:,且.,a b 23a b ab +=ab 982a b +832a b +1b a-3-()[]f x x =[]x x {}n a n n S 112n n n S a a ⎛⎫=+ ⎪⎝⎭21n n n b S S +=+()*n a n n =∈N)*n S n =∈N []12636b b b +++= 1210011118S S S ⎡⎤+++=⎢⎥⎣⎦ ()()2222ln f x x f x x -'=+()2f '=ABC A B C 、、a b c 、、sin a C =π6A =22b c bc+=()()()()13e 0xf x a x b a =-++≠[]1,3-3b a+()π,cos ,cos ,sin 2m x x n x x ⎫⎛⎫=-= ⎪⎪⎝⎭⎭()32f x m n =-⋅()f x ABC ()0,ABC f A BC S ===AC AB {}n a 421a =125,,a a a {}n b 143n n b b +=-1121b a =-(1)求和的通项公式;(2)若为数列的前项和,求.17.(本题满分15分)东风学校有甲乙两个食堂,学校后勤服务中心为了调查学生对两个食堂的满意度,随机调査300名学生.设表示事件“学生喜欢去甲食堂”,表示事件“调査的学生是男生”.若.调查的是男生调查的是女生合计喜欢去甲食堂喜欢去乙食堂合计(1)完成上列列联表,并根据小概率值的独立性检验,判断学生喜欢去哪个食堂与性别是否有关?(2)为了答谢参与调查的学生,学校后勤服务中心从参与调查的300名学生中按性別分层抽样的方法选15名幸运学生参与抽奖活动,并为他们准备了15张奖券,其中一等奖奖券有3张,二等奖奖券有5张,三等奖奖券有7张,每人抽取一张.设15名幸运学生中男生抽中一等奖的人数为,写出的分布列,并计算.附0.10.050.010.0050.0012.7063.8416.6357.87910.82818.(本题满分17分)已知函数.(1)讨论的单调性;(2)当时,恒成立,求实数的取值范围;(3.19.(本题满分17分)马尔科夫链是一种随机过程,它具有马尔科夫性质,也称为“无记忆性”,即一个系统在某时刻的状态仅{}n a{}n b n T1n n a b ⎧⎫⎨⎬-⎩⎭n n T M N ()()()457|,|,7815P M N P N M P N ===22⨯0.001α=X X ()E X ()()()()22():ad bc na b c d a c b d χ-⋅=++++αax ()1ln f x x a x x=--()f x 1x …()0f x …a ()ln 1n ++>+与前一时刻的状态有关.为了让学生体验马尔科夫性质,数学老师在课堂上指导学生做了一个游戏.他给小明和小美各一个不透明的箱子,每个箱子中都有个红球和1个白球,这些球除了颜色不同之外,其他的物质特征完全一样规定“两人同时从各自的箱子中取出一个球放入对方的箱子中”为一次操作,假设经过次操作之后小明箱子里的白球个数为随机变量,且.(1)求的值;(2)求;(3)证明:为定值.x n n X ()1518P X ==x ()1n P X =()n E X2024年秋季普通高中11月份阶段性联考高三数学试卷参考答案一、选择题:本题共8小题,每小题5分,共40分.1.D2.B3.B4.C5.C6.B7.A8.D8.【解析】,显然不满足上式,所以,令,则,在,且,画出的图像,可知:.二、选择题(多选)【有错选得0分,全对得6分,部分对得部分分.两解题,每答对一个得3分,三解题,每答对一个得2分】9.ACD 10.BCD11.BCD10.解析:(1)(当时取等号);(2)(当时取等号);()()()32481f x f x x x x m x '=⇒--+=-1x =32481,1x x x x m x--+≠=-()32481x x x g x x --+=-()()()22221(1)x x g x x '-+=--()g x ∴[)(4,1,1,2,2,3⎤⎤⎡-↑↑↓⎦⎣⎦()()()564,24,375g g g -=-=-=-[)7,4m ∈--8329ab a b ab =+≥⇒≥⇒≥24,33a b ==8233a b ab +=≥24,33a b ==(3)(当时取等号);(4)(当时取等号).11.解析:(1)当时,,又A 错,B 对;(2),.故C 对;(3),当时,,,;故D对;三、填空题:12.13.14.14.【解析】,令,在,在,()()212122233,3225923a b a b ab a b a b a b b a b a b a ⎛⎫+=⇒+=∴+=++=++≥⇒+≥ ⎪⎝⎭1a b ==132233b b b b a b b --=-=+-≥-b =11,2n n nS a a ⎛⎫=+∴ ⎪⎝⎭2n ≥2211112,1n n n n n n n S S S S S S S ---=-+⇒-=-11111,02n S a a a ⎛⎫=+> ⎪⎝⎭211;n n n a S n S a ⇒=∴=⇒==∴()1263211176,722n n n b b b b S S +===-∴+++=+-∈+ []12636b b b ∴+++= 12n S =>=]1210011122118;S S S ⎡⎤∴+++>+++=->⎣⎦2n ≥12n S =<=-]121001111212119S S S ⎡⎤∴+++<++++=+-=⎣⎦1210011118S S S ⎡⎤∴+++=⎢⎥⎣⎦ 3-21,2e e ⎡⎤-⎢⎥⎣⎦()()()03e 1;x f x b a x =⇔+=-310,e x b x a a +-≠∴= ()()12,e ex x x x g x g x --=='()g x ∴()1,2-↓()2,3↑作出的图像,可知:.四、解答题:本题共5小题,共77分.解答应写出文字说明、证明过程或演算步骤.15.(本题满分13分)解:(1)由减区间为(2),或.16.(本题满分15分)解:(1)设的公差为,又(2),两式相减,得:17.(本题满分15分)()g x 2132e e b a+-≤≤()23π3cos cos sin sin 222f x x x x x x x ⎛⎫=---=- ⎪⎝⎭()311π1cos21cos2sin 21,2226x x x x x ⎫⎛⎫=--=--=--+⎪ ⎪⎪⎝⎭⎭πππππ2π22πππ,26263k x k k x k -+≤-≤+⇒-+≤≤+()f x ∴()*πππ,π63k k k ⎡⎤-++∈⎢⎥⎣⎦N ()ππ0sin 21,,63f A A A ⎛⎫=⇒-== ⎪⎝⎭6,ABC S AB AC =⇒⋅= 227,BC AB AC AB AC =⇒+-⋅=2,3AB AC ∴==3,2,AB AC ==⋅{}n a ()()()221520,,21321(212)6d d a a a d d d d ≠=∴-+=-⇒= ()14133,16 3.n a a d a a n d n ∴=-==+-=-()1143141,n n n n b b b b ++=-⇒-=-111215,14,b a b =-=-=()*1441n n n n b b n ∴-=⇒=+∈N 6314n nn a n b -=-2323411633915631391563;;4444444444nn n n n n k n n n T T +=---==++++∴=++++∑2341336666635165;4444444334n n n n n n n T T +-+=+++++-⇒=-⋅解:(1)被调查的学生中男生有140人,女生有160人.男生中喜欢去乙食堂的有80人,喜欢去甲食堂的有60人..被调查的学生中喜欢去甲食堂的有160人.调查的是男生调查的是女生合计喜欢去甲食堂60100160喜欢去乙食堂8060140合计140160300零假设:假设学生喜欢去哪个食堂与性别无关.,根据小概率值的独立性检验,我们推断不成立,即认为学生喜欢去哪个食堂与性别有关.此推断犯错误的概率不大0.001.(2)根据男女生人数之比可知,被抽取的15人中男生7人,女生8人.,,X 的分布列为:X 0123p,18.(本题满分17分)解(1)定义域为;..当时,恒成立,;()77,300140,1515P N =⨯=∴44(),14080,77P M N =⨯=∴∣533()(),60160,888P N M P N M =⇒=÷=∴∣∣0H 220.001(606010080)30011.5810.828160140160140χχ⨯-⨯⨯=≈>=⨯⨯⨯0.001α=0H 0,1,2,3X =()()()()615243712312312312777715151515C C C C C C C 8282450,1,2,3C 65C 65C 65C 65P X P X P X P X ============86528652465113()82824570123656565655E X =⨯+⨯+⨯+⨯=()0,∞+()()22211,Δ4,f x x ax a x=-+=-⋅'0122a -≤≤2Δ0,10x ax ≤-+≥()()0,f x f x ≥↑'.当时,有两根,但两根均为负数,当时,.当时,有两正根,当时,;当时,;当时;综上所述:.当时,增区间为;.当时,增区间为和;减区间为.(2),令,则在,若,则,与题意相符;若,则,所以必存在,使得当时,,从而使得当时,,与题意相矛盾;综上:.(3)证明:由(2)知,当时,(仅当时取等号),,令;,得证.19.(本题满分17分)解:(1)(2)022a<-2Δ0,10x ax >-+=()0,x ∞∈+()()0,;f x f x '≥↑32a >2Δ0,10x ax>-+=1x =2x =()10,x x ∈()()0,f x f x >↑'()12,x x x ∈()()0,f x f x <↓'()2,x x ∞∈+()(),0,f x f x >'↑012a ≤()f x ()0,∞+022a >()f x ⎛ ⎝∞⎫+⎪⎪⎭()11f x x a x x ⎛⎫=+- ⎪⎝⎭'()1g x x a x =+-()()()22110,g x x g x x =-≥∴'[)()1,,12g a ∞+↑=-2a ≤()()()()()()10,0,,10g x g f x f x f x f ≥≥≥↑≥='2a >()120g a =-<01x >()01,x x ∈()()()0,0,g x f x f x <'<↓()01,x x ∈()()10f x f <=2a ≤1x ≥()12ln 0f x x x x=--≥1x =12ln x x x∴-≥x =11ln ln n n n n ++>=⇒>()2341ln ln ln ln ln 1123n n n +>+++=+ ()111513;11118x x P x x x x x x ==⋅+⋅=⇒=++++()()()()()()()11111010111212n n n n n n n n n n P x P x P x x P x P x x P x P x x ++++===⋅==+=⋅==+=⋅==∣∣∣,又,.(3),令,则而,..得证.()()()()()()11331111510120122244442282n n n n n n P x P x P x P x P x P x ⎛⎫==⋅+=⋅⨯+⨯+=⋅==+=+= ⎪⎝⎭()()()0121n n n P x P x P x =+=+==()()()()()()11151141411111,11,2882787n n n n n n P x P x P x P x P x P x ++⎡⎤⎡⎤∴==-=+===+⇒=-==-⎣⎦⎢⎥⎣⎦()()()114543431314311,11;78756756878778n n nn n P x P x P x -⎛⎫⎛⎫⎛⎫=-=-=∴=-=⨯=⨯⇒==+⨯ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭ ()()()()()1112020121n n n n n n n P x P x P x x P x P x x +++===⋅==+=⋅==∣∣()()1222n n n P x P x x ++=⋅==∣()()()1311913122162214828n n n n P x P x P x +⎛⎫==+===++ ⎪⎝⎭()()()()111131391339228248214214148141414n n n n n n n P x P x P x P x ++++⎡⎤⎛⎫⎡⎤⎡⎤⇒=-==-+⇒=-=⨯=-+ ⎪⎢⎥⎢⎥⎢⎥⎣⎦⎝⎭⎣⎦⎣⎦()38214n n n a P x ⎡⎤==-⎢⎥⎣⎦1193344,141414n n n n a a a a ++⎛⎫=+⇒+=+ ⎪⎝⎭()113333338280141414161414a P x ⎡⎤⎡⎤+==-+=-+=⎢⎥⎢⎥⎣⎦⎣⎦()()3333310820214141414148n n n n n a P x P x ⎡⎤∴+=⇒=-+=⇒==-⨯⎢⎥⎣⎦()()()()43133100112212177814148n n n n n n E X P x P x P x ⎡⎤⎡⎤=⨯=+⨯=+⨯==⨯+⨯+⨯-⨯=⎢⎥⎢⎥⎣⎦⎣⎦。
高尔顿钉板数学解释
高尔顿钉板是一种几何工具,也叫高尔顿板或者布朗尼安运动。
它由一块平面上的钉板组成,其上方有许多等间距的突起物(通常是钉子),形成一系列固定的排列。
高尔顿钉板的数学解释可以通过落球模型来理解。
假设在钉板上方有一个小孔,并且从该小孔中往下自由下落一个小球。
每当小球碰到钉板时,它会按照一定的规则反射或者改变运动方向,最终最终落到一个特定的位置。
通过观察这些球的运动轨迹,我们可以观察到一些有趣的现象。
数学上,高尔顿钉板的解释可以通过物理学中的弹性碰撞和反射来进行分析。
当小球经过碰撞后,其运动方向和速度会发生改变,但是总能量是守恒的。
通过数学建模和计算,我们可以推导出小球的轨迹和最终的落点。
高尔顿钉板模型在数学研究中被广泛应用于混沌理论、动力系统和分形几何等领域。
它展示了一种看似无序但却是有规律的运动,引发了人们对随机性、确定性和无序性的思考。
总结而言,高尔顿钉板的数学解释可以通过物理学中的弹性碰撞和反射来理解。
通过观察球的运动轨迹,我们可以观察到一些有趣的现象,这在数学研究中具有重要的意义。
高尔顿(Galton)钉板实验
一、问题描述
Galton钉板试验是英国生物统计学家Galton设计的。
在一板上钉有n排钉子,如图示,其中n=5。
右图中15个圆
点表示15颗钉子,在钉子的下方有n+1个各子,分别编号为
0,1,2,…,n。
从Galton钉板的上方扔进一个小球任其自
由下落,在下落的过程中当小球碰到钉子时,从左边落下与从
右边落下的机会相等。
碰到下一排钉子时又是如此。
最后落入
底板中的某一个格子,图中用一条折线显示小球下落的一条轨迹。
二、高尔顿钉板试验中的相关问题
1、小球落入各个格子中的概率与频数
做一个小球的高尔顿钉板试验,其落入第i个格子的概率正好满足二项分布。
设高尔顿钉板有n行钉,第n行铁钉共有n个,有(n+1)个空。
把这(n+1)个空由左到右依次编号为i=0,1,2,…,n共(n+1)个空。
观察i=0这个空,小球从这个空落下的条件是:小球从第一次与铁钉碰撞后必须连续向
左落下,即连续n次选择向左落下,所以落入第i=0个空的概率为P(i=0)=C0
n (
2
1
)n(
2
1
)0。
观察i=1这个空,小球从这个空落下的条件是:小球从第一次与铁钉碰撞后连续n次碰撞落下过程中,有且只有一次选择向右落下,其余都只能是向左落下,所以落入第i=1个空
的概率为P(i=1)=C1
n (
2
1
)n-1(
2
1
)1。
小球从第一次与铁钉碰撞后连续n次碰撞落下过程中,有i次选择向右落下,其余都选择向左落下,所以落入第i个空的概率为
P(i)= C i
n (
1
2
)n-i(
2
1
)i(i=0,1,2,…,n)。
故,当一个一个从顶部放入k个小球,低槽中各格的理论频数为:
h(i)=k×P(i),(i=0,1,2,…,n).
2、程序运行
2.1 基本功能
①输入小球数k、概率p;
②计算高尔顿钉板n=4时,放入k个小球后,落入底槽各格中的实验小球数;
③计算高尔顿钉板n=4时,放入k个小球后,落入底槽各格中的理论小球数;
④动画演示每个小球下落路径及底槽各格小球数频率增长情况;
④画出落入底槽各格中的实验小球数频率的柱状图;
⑤画出落入底槽各格中的实验小球数、落入底槽各格中的理论实验小球数的频率曲线图;
⑥关闭。
2.2 图形界面
2.3 程序使用方法
➢打开ex01.fig;
➢填写“输入球数”、“输入p值”,点击求值运行;
➢实验值一栏显示为落入底槽各格中的实验小球数,理论值一栏显示为落入底槽各格中的理论小球数;
➢界面动画演示小球下落路径;
➢红线代表落入底槽各格中的实验小球频率分布,绿线代表落入底槽各格中的理论小球频率分布,便于比较;
➢多次点击“求值”按钮,可得到多组实验数据;
➢单击“关闭”按钮,可退出程序窗口。
2.4 运行结果
(1)动画演示:小球数k=500 p=0.3
(2)小球数k=500 p=0.5
由运行结果我们可以看出,所做实验中落入底槽各格中的实验小球频率与落入底槽各格中的理论小球频率基本符合,但存在较大误差。
(3)小球数k=5000 p=0.5
由运行结果我们可以看出,所做实验中落入底槽各格中的实验小球频率与落入底槽各格中的理论小球频率基本符合,误差相对(1)有所减小。
(4)小球数k=50000 p=0.5
由(2)、(3)、(4)运行结果我们可以看出,随着小球数的增加,实验值与理论值
越来越符合,小球落入底槽各格中的概率也愈加符合二项分布、各格小球数符合正态分布。
(5)小球数k=500 p=0.1
(6)小球数k=500 p=0.3
2.5卡方检验
由2.4(2)运行结果知:
p=[0.0625,0.25,0.375,0.25,0.0625];
v=[31.25,125,187.5,125,31.25];
function kf
p=[0.0625,0.25,0.375,0.25,0.0625];
V=[31.25,125,187.5,125,31.25];
for i=1:5
V=0;
V(i)=(v(1,i)-k*p(1,i))^2/(k*p(1,i));
V=V+V(i);
end
运行结果与9.488相比,可得到V<9.488,因此验证所实验分布服从B(4,p)。
3、matlab模拟Poisson分布与几何分布
3.1 Poisson分布
3.1.1高尔顿钉板模拟Poisson分布
➢数学原理
对于参数n,p二项分布,当n很大p很小,近似于λ=n*p为参数的泊松分布。
➢取n=50000,p=0.3
3.1.2时间间隔序列模拟Poisson分布
➢数学原理
在实际事例中,当一个随机事件,例如某电话交换台收到的呼叫、来到某
共汽车站的乘客、某放射性物质发射出的粒子、显微镜下某区域中的白血球等等,以固定的平均瞬时速率λ(或称密度)随机且独立地出现时,那么这个事件在单位时间(面积或体积)内出现的次数或个数就近似地服从泊松分布。
由具有指数分布的时间间隔序列模拟泊松过程,T n是由random (‘exponential’,λ)生成,T n间相互独立。
➢程序内容
function Poisson
λ=2;T max=50;
i=1;
T(1)=random('exponential',λ); %产生服从参数为λ的指数分布随机数T(1) while(T(i)<T max)
T(i+1)=T(i)+random('exponential',λ);
i=i+1;
end
T(i)=T max; x=0:1:i;w(1)=0;
for p=1:i
w(p+1)=T(p);
end
stairs(w,x);
➢运行结果
05101520253035404550
05
10
15
20
25
3.2模拟几何分布
➢ 数学原理
几何分布可定义为n 次伯努利实验中,实验k 次才得到第一次成功的概率。
在高尔顿钉板试验中,我们可以认为是第m (m=1,2,3,…,n )个小球第一次落入第i (i=0,1,2,3,4)个底槽中的概率。
四、参考资料与任务分配
1、参考资料
【1】《数学实验》 焦光虹主编 科学技术出版社
【2】《MATLAB 语言与数学建模》 曾建军 李世行等编著 安徽大学出版社
【3】《概率论与数理统计》 哈尔滨工业大学数学系组编 科学出版社
【4】《高尔顿(Galton )》 数理统计
【5】《高尔顿钉板实验的算法实现与分析》 聂燕著 中国民航飞行学院学报 May.2008 Vol.19 No.3
五、总结
我们所编高尔顿钉板matlab程序基本满足实验要求,可以通过改变小球数k值、概率p值,从而得到不同情况下底槽各格中小球的分布情况,并与理论值相比较,对实验结果的正确与否进行验证。
同时,程序也给出了底槽各格中实验小球数与理论小球数的频率分布曲线,从而可以对其有更直观地认识。
但我们的程序也有明显的缺点,如程序只给出了试验所要求的高尔顿钉板的层数n(程序中n=4),不能进一步模拟得到更多的数据,程序给出小球下落轨迹的动画演示但当小球数很多时程序运行时间很长,程序编学较为繁琐等。
通过该次实验,我们对《数学实验》这一课程有了进一步的了解与掌握,对于matlab 程序的编学、matlab语言的应用有了更加深刻地体会,同时我们在实际问题在应用数学知识的能力得到加强。