- 1、本文档共13页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
概率论实验报告 -赵恒伟
概率论与数理统计实验报告
概率论与数理统计实验
电信学院 信息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-
文档评论(0)