概率论实验报告 -赵恒伟.docx

概率论与数理统计实验报告 概率论与数理统计实验 电信学院 信息24班 2120502099 赵恒伟：1871036262 2120502079 尚 悦2120502089 鲁 龙2014-6-15 XJTU  【实验选题】 4，7，10，15，19题 【实验目的】 熟练掌握MATLAB软件的关于概率论与数理统计的基本操作 应用概率论与数理统计知识解决实际问题 问题一： 4.俗语说“三个臭皮匠顶一个诸葛亮”。假设诸葛亮成功解决某一棘手问题的概率为 0.9；0.95；0.99。试就三个“臭皮匠”解决问题的能力不同，验证这句俗语。树形图分析。 【分析与解答】 设事件：一个臭皮匠的成功处理一个问题 树形图分析如下： 事件 处理问题 未能处理问题 由树形图可知，三个“臭皮匠”处理事件的概率： 因为三??臭皮匠处理问题能力不同，这里随机数法来决定他们的处理问题能力，随机数采用正态分布，因为一般人不可能每次都能处理问题，而且不按照处理能力均匀分布，故可设定平均处理能力=0.5，而且个人处理能力不超过1，所以根据原则，即约为0.5，。此时当能力时，位于之外，几乎可以忽略。 计算三个臭皮匠处理问题的能力： Matlab代码： p=0.9; t=0; for i=1:1:1000 a1=random(norm,0.5,0.33,1,1); %臭皮匠1能成功解决问题的概率。 a2=random(norm,0.5,0.33,1,1); %臭皮匠2能成功解决问题的概率。 a3=random(norm,0.5,0.33,1,1); %臭皮匠3能成功解决问题的概率。 a=1-((1-a1)*(1-a2)*(1-a3)); %三个臭皮匠能解决问题的概率。 if ap %三个臭皮匠顶一个诸葛亮一次成功。 t=t+1; %记录一次。 end end x=t/1000 %1000次试验后三个臭皮匠顶一个诸葛亮的概率。 运行结果：x = 0.5570 同上，当p=0.95,运行结果x =0.4130；当p=0.99，运行结果x = 0.2400。 【结果与讨论】 由该结果可知，诸葛亮的处理问题能力越强，“三个臭皮匠胜过诸葛亮”这一结论成立的可能性越低；当“臭皮匠”平均处理问题能力越高，结论成立可能越大。 问题二： 7. 在非血缘群体中，骨髓配型的概率为十万分之一。今有一白血病患者希望骨髓移植。若该患者能成功找到配型的概率为0.5；0.8；0.99；问骨髓库至少该有多少份志愿捐献者的资料？ 【分析与解答】 该问题有两种解题思路，一种是利用古典概型，较为简单，另一种是结合中心极限定理。 解法一： 设Ai表示第i个志愿者能成功配型，i=1,2,…n，则P()=1/100000，利用独立性有，P(患者能成功找到配型)=P(A1A2…An)=1-(99999/100000)^n根据上述方程Matlab求解n。 程序设计： syms p n [n,p]=solve(1-p=(99999/100000)^n,p=0.5); n1=vpa(n,10) [n,p]=solve(1-p=(99999/100000)^n,p=0.8); n2=vpa(n,10) [n,p]=solve(1-p=(99999/100000)^n,p=0.99); n3=vpa(n,10) 实验结果： n1 =69314.37148 n2 =160942.9865 n3 =460514.7160 所以为了使配型概率为0.5， 0.8， 0.99.至少需要志愿者分别： n1 =69314.37148， n2 =160942.9865， n3 =460514.7160 解法二： 一个志愿捐献者配型这一事件， 用表示成功配型的人数。 用q表示配型成功的概率，q=0.5, 0.8, 0.99; ,为运算方便取等式。 利用代码： syms p q q=0.5; p=1/100000; x=norminv(1-q,0,1) 可求得：q=0.5, 0.8, 0.99时， 对应x=0， -0.8416， -2.3263，在将x带入， 下面程序： [n,p,q]=solve(0=(1-