推荐概率统计与计算机仿真.ppt

* * * 概率统计与计算机仿真 1 概率统计问题的MATLAB实现 一. 随机数的产生 二项分布的随机数据的产生 命令 参数为N，P的二项随机数据 函数 binornd 正态分布的随机数据的产生 R = normrnd(MU,SIGMA) 或者 randn 通用函数求各分布的随机数据 命令 求指定分布的随机数 函数 random rn = random(‘Normal’,0,1,2,4) %产生一个2×4的数组，服从 %均值为0,方差为1的正态分布 rp = random(‘Poisson’,1:6,1,6) %产生一个1×6的数字，在1 ％到6上服从泊松分布， 产生均匀分布随机数 命令 在[0 1]上均匀分布的随机数 函数 rand 交互式随机数工具：randtool 二.统计数据的图形处理 （1）饼形图和条形图 例 全国不同教育水平人群吸烟率: 序号 1 2 3 4 教育水平 大学 中学 小学 文盲 男性 44.72 53.75 66.47 67.83 女性 3.61 1.93 6.10 11.56 x=[44.72 53.75 66.47 67.83; 3.61 1.93 6.10 11.56]; x=x; colormap jet bar(x,1.2,group) legend(Male,Female,-1) （2）直方图 将数据分成若干个小组；统计出数据落在每 组内的个数，称为频数，形成一个频数表；最后 以数据的取值为横坐标轴，频数为纵坐标轴，在 每一组上竖起一个矩形，这样的图形就称为直方图. Matlab作直方图的命令是hist [n, x]=hist(a, k) 将数组a所在的区间分为k份， x返回各个小区间的中点，n返回a在各小区间的频数. hist(a, k) 绘出数组a的直方图，k的含义同上， k的默认值为10. 直方图常用于显示数据的统计分布规律， （1）通用函数计算概率密度函数值 命令 通用函数计算概率密度函数值 函数 pdf （2）专用函数计算概率密度函数值 如：binopdf ，poisspdf 三 概率密度函数值和积累概率 （3）通用函数计算累积概率值 命令 通用函数cdf用来计算随机变量X （4）专用函数计算累积概率值 如：normcdf binocdf （5）逆累积分布函数值的计算 MATLAB中的逆累积分布函数是已知 求x。 通用函数计算逆累积分布 命令 icdf x=icdf(norm,0.975,0,1) 专用函数*inv计算逆累积分布函数 例 公共汽车门的高度是按成年男子与车门顶碰头 的机会不超过1%设计的。设男子身高X服从正态分 布N（175，36），求车门的最低高度。 h=icdf(norm,0.99,175,6) h=norminv(0.99, 175, 6) 2 计算机仿真及其matlab实现 一、 什么是计算机仿真 计算机仿真是利用计算机对一个实际系统的结构和行为进行动态演示，以评价或预测该系统的行为效果。它是解决较复杂的实际问题的一条有效途径。 计算机仿真在航空、机电、冶金等工业部门及社会经济、交通运输、生态系统等方面有着广泛的应用，已成为分析、研究和设计各种系统的重要手段。 二、 为什么要进行计算机仿真 1)在一个实际系统还没有建立起来之前，要对系统的行为或结果进行分析研究时，计算机仿真是一种行之有效的方法。 2)在有些真实系统上做实验会影响系统的正常运行，如在生产中任意改变工艺参数可能会导致废品，在经济活动中随意将一个决策付诸行动可能会引起经济混乱。 3)当人是系统的一部分时，他的行为往往会影响实验的效果，这时最好对系统进行仿真研究。 4)在实际系统上做实验时，很难保证每次操作的条件相同，因而对实验结果的好坏很难作出正确的判断。 5)在实际系统上做实验可能会花费很长的时间，费用太高或者有危险，使得实验不易进行。 6)有些系统一旦建立起来之后就无法复原。例如要投资建立一个大型企业，要分析它建成之后的经济效益和社会效益，不能用建立起来试试看的办法，因为建成后就无法回到原来的状态了。 1)难以用数学公式描述的系统，或者没有求解的有效数学方法； 2)虽然可以用解析方法解决问题，但数学的分析与计算过于复杂，这时计算机仿真可能提供简单可行的求解方法； 3)希望能在较短的时间内观察到系统发展的全过程，以估计参数对系统行为的影响； 4)难以在实际环境中进行实验和观察时，计算机仿真