matlab应用之Galton抽奖模型的奖品设计概述.pptVIP

一、实验目的 　 用二项分布描述Galton钉板模型，模拟Galton钉板试验，观察和体会概率分布列的意义； 　　用Galton钉板实验刻画抽奖模型，解决实际生活中的问题。 　　 所有抽奖模型都是要赚钱的，没有人想花费精力却一无所获。以Galton钉板实验为例，如在各个格子中设计适当价值的奖品用以获取利润，当然是抽奖主办者所关心的事情。 　　Galton钉板模型可以看成一个有趣的游戏。小球自上方落下，经过个钉子，每经过一个钉子时只有两种结果：向左或向右，这是一个具有两个结果的随机试验。将向右视为成功，成功的概率为，向左视为失败，失败的概率为，小球碰到一个钉子下落一格，相当于进行了一次试验。小球自顶端落下，碰到个钉子，最终落在某个格子的过程，相当于一个重伯努利试验。重伯努利试验的成功次数正好是小球向右移动的次数，它是一个随机变量。根据概率论的结果有。 　　 　　 所有抽奖模型都是要赚钱的，没有人想花费精力却一无所获。以Galton钉板实验为例，如在各个格子中设计适当价值的奖品用以获取利润，当然是抽奖主办者所关心的事情。 　　Galton钉板模型可以看成一个有趣的游戏。小球自上方落下，经过个钉子，每经过一个钉子时只有两种结果：向左或向右，这是一个具有两个结果的随机试验。将向右视为成功，成功的概率为，向左视为失败，失败的概率为，小球碰到一个钉子下落一格，相当于进行了一次试验。小球自顶端落下，碰到个钉子，最终落在某个格子的过程，相当于一个重伯努利试验。重伯努利试验的成功次数正好是小球向右移动的次数，它是一个随机变量。根据概率论的结果有。 　　 　 现考虑5层钉板，这时有6个格子，首先对所给的，必须了解大量投球后小球的堆积形状，考虑一种设计方案，给定两个参数和，在0~10元中取值，在0~1之间取值。各个格子的奖品价值为，即编号为0的格子设置价值为的奖品，编号为1的格子设置价值为奖品，依此类推。 　　对于不同的和，根据服从二项分布的随机变量的分布列，计算数学期望，它是和的二元函数。如果规定抽奖者在每次抽奖中付出一元的代价，则就是主办者在一次抽奖中所获得的平均利润。注意主办者获利的前提是有大量的人参加抽奖，所以一方面希望利润较大，另一方面也需对抽奖者具有一定的吸引力。试对二者进行权衡，为该抽奖模型配置适当的奖品，以获得较高的利润。 三. 实验步骤 　　　一、动画模拟Galton钉板试验 1) 确定钉子的位置。将钉子的横、纵坐标存储在一个矩阵中； 2) 模拟了小球从顶端随机地落入某一格子的过程。设向右的概率为，向左的概率为；将[0,1]分成两段，区间[0,p]和(p,1]。利用rand[]产生一个介于0和1之间的随机数，如果随机数，让小球落向左边，否则落向右边；将这一过程重复n次，并用直线连接小球落下时所经过的点。 3) 模拟小球堆积的形状。输入扔球次数，计算落在第个格子的小球数在总球数中所占的比例，这样当模拟结束时，就得到了频率，用频率反映小球堆积的形状。 4) 利用movie完成动画。 二、用二项分布描述Galton钉板模型 　　 根据概率论的结果，它服从二项分布，即,其取值与模拟模型的对应关系为： 　　　 　　　 三、求出主办者的获利，即随机变量的数学期望 　　求出主办者的获利，即随机变量的数学期望。 　　给定两个参数和，在0~10元中取值，在0~1之间取值。各个格子的奖品价值为 　　利用概率论知识，二项随机变量的分布列为： 　　　 　　则 五．实验程序 程序一：动画模拟Galton钉板试验 clear; clc; clf; m=input(请输入小球的个数:); n=6; y0=2; ballnum=zeros(1,n+1); p= input(请输入p的值:); q=1-p; for i=n+1:-1:1 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); 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 y