概率论与数理统计实验.pptx

实验目的;;;调用格式：; 当只知道一个随机变量取值在（a，b）内，但不知道（也没理由假设）它在何处取值的概率大，在何处取值的概率小，就只好用U（a，b）来模拟它。;例1、产生U(2,8)上的一个随机数，10个随机数， 2行5列的随机数。;2、正态分布随机数; 当研究对象视为大量相互独立的随机变量之和，且其中每一种变量对总和的影响都很小时，可以认为该对象服从正态分布。;3、指数分布随机数;结果 ( (2) y2=[0.1465 0.0499 0.0722 0.0115 0.0272 0.0784 0.3990 0.0197 0.0810 0.0485 0.0233] (3) ； 0.0207 0.1974 0.1616 0.1301 0.4182 0.0809];排队服务系统中顾客到达率为常数时的到达间隔、故障率为常数时零件的寿命都服从指数分布。;4、二项分布随机数;20、其他分布随机数的产生方法;（一）直接抽样法（反函数法）;设分布律为P(X=xi)= pi , i=1, 2, ...，其分布函数为F(x) ;（二）变换抽样法;例5、 设X分布函数为;生成n=20的1行10列随机数;例6 生成单位圆上均匀分布的1行10000列随机数，并画经验分布函数曲线。;第20页/共48页;例6 频率的稳定性;掷一枚硬币，记录掷硬币试验中频率P*的波动情况。;在Matlab中编辑.m文件输入以下命令：;在Matlab命令行中输入以下命令：;在Matlab命令行中输入以下命令：;在Matlab命令行中输入以下命令：;1、表示位置的统计量—平均值和中位数;2、表示变异（离散）程度的统计量 —方差、标准差、极差 ;例7：产生5组正态分布随机数，每组100个，计算均值、中位数、标准差、极差、方差。;3.表示分布形状的量—偏度和峰度 ;命令：y=skewness(x) 若x为矢量，返回样本偏度，若x为矩阵，结果为行向量，每个元素对应x中列元素的样本偏度。;分布峰度 ;命令：y=kurtosis(x) 若x为矢量，返回x的样本峰度，若x为矩阵，结果为行向量，每个元素对应x中列元素的样本峰度。;例8：产生5组正态分布随机数，每组100个，计算样本偏度和峰度。;三、分布函数的近似求法（直方图）;;1、整理资料： 把样本值x1，x2，…，xn进行分组，先将它们依大小次序排列，得 ;2、求出各组的频数和频率：;3、作图： ;(1)、描绘数组data的频数直方图的命令为： hist(data, k); 例10 一道工序用自动化车床连续加工某种零件，由于刀具损坏等会出现故障.故障是完全随机的，并假定生产任一零件时出现故障机会均相同.工作人员是通过检查零件来确定工序是否出现故障的.现积累有100次故障纪录，故障出现时该刀具完成的零件数如下：作频数直方图及经验分布函数图;459 362 624 542 509 584 433 748 815 505 612 452 434 982 640 742 565 706 593 680 926 653 164 487 734 608 428 1153 593 844 527 552 513 781 474 388 824 538 862 659 775 859 755 49 697 515 628 954 771 609 402 960 885 610 292 837 473 677 358 638 699 634 555 570 84 416 606 1062 484 120 447 654 564 339 280 246 687 539 790 581 621 724 531 512 577 496 468 499 544 645 764 558 378 765 666 763 217 715 310 851;解;hist(x,10) histfit(x,10);第45页/共48页;作经验分布函数图;例11、产生参数为10的指数分布随机数500个， 并画出直方图和经验分布函数图。;感谢您的观看！