西安交大概率论上机实验报告Word格式文档下载.docx
- 文档编号:14215255
- 上传时间:2022-10-20
- 格式:DOCX
- 页数:10
- 大小:112.96KB
西安交大概率论上机实验报告Word格式文档下载.docx
《西安交大概率论上机实验报告Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《西安交大概率论上机实验报告Word格式文档下载.docx(10页珍藏版)》请在冰豆网上搜索。
2014161114281813273125241619232617143021
1816181920221922182626132113111923182428
1311251517182216131213110915182115121713
1412161008231811162813212212081521181616
1928191214192828281321281911151824181628
1915132214162420281818281413282924281418
1818082116243216281915181810121626181933
0811182723112222132814221826181632272524
1717283316202832192318281524282916171918]
6.利用Matlab软件模拟高尔顿板钉试验。
7.自己选择一个与以上问题不同类型的概率有关的建模题目,并解决.
实验任务及结果
任务一、掌握常见分布的概率密度及分布函数的命令
分布函数
概率密度命令
分布函数命令
二项分布
binopdf(x,n,p)
binocdf(x,n,p)
泊松分布
poisspdf(x,lamda)
poisscdf(x,lamda)
均匀分布
unifpdf(x,a,b)
unifcdf(x,a,b)
正态分布
normpdf(x,mu,sigma)
normcdf(x,mu,sigma)
几何分布
geopdf(x,p)
geocdf(x,p)
超几何分布
hygepdf(x,m,k,n)
hygecdf(x,m,k,n)
t分布
tpdf(x,v)
tcdf(x,v)
F分布
fpdf(x,v1,v2)
fcdf(x,v1,v2)
任务二、利用二项分布命令计算抛硬币实验
题目分析:
掷硬币是一种简单的随机试验,服从二项分布b(n,0.5),利用MATLAB中的概率密度命令与分布函数命令,取参数n为实验次数150,参数x为计算数值45,可以直接得到结果。
程序代码:
p1=binopdf(45,150,0.5)
p2=binocdf(45,150,0.5)
任务三、使用随机数命令产生服从二项分布的随机数并验证泊松定理
1、利用binornd(n,p,m,s)可以直接产生m行s列服从b(n,p)的随机数。
2、泊松定理的内容是:
在n重贝努力试验中,事件A在每次试验中发生的概率为p,出现A的总次数K服从二项分布b(n,p),当n很大p很小,λ=np大小适中时,二项分布可用参数为λ=np的泊松分布来近似。
为了验证泊松定理,可以设置参数n、p,通过二项分布命令binopdf与泊松分布命令poisspdf分别计算出分布律,并作图对比。
%使用binornd命令产生服从二项分布b(n,p)的随机数
n=10000;
p=0.01;
binornd(n,p,1,10)%产生服从b(n,p)的随机数
x=50:
150;
y1=binopdf(x,n,p);
%利用二项分布计算分布律,用空心圈绘出
y2=poisspdf(x,n*p);
%利用泊松分布计算分布律,用星号绘出
plot(x,y1,'
o'
x,y2,'
*'
);
xlabel('
x'
ylabel('
P(X=x)'
运行结果与分析:
(1)服从二项分布的随机数
程序使用n=10000,p=0.01的二项分布,产生10个随机数结果如图,可以看出产生的10个随机数都在np=100附近。
(2)泊松定理的验证
图中空心圈为利用二项分布计算分布律结果,星号为利用泊松分布计算分布律结果,从图可以看到,两种分布计算结果几乎完全重合,即在这种条件下二项分布完全可以用泊松分布逼近,验证了泊松定理。
任务四、画出二维随机变量的概率密度函数图像
利用MATLAB命令ezsurf可以非常简单地画出二维函数图像。
ezsurf('
1/(2*pi)*exp(-(x^2+y^2)/2)'
)
任务五、调用MATLAB函数计算一组数据的均值、方差并绘直方图
MATLAB自带的mean命令可以直接计算一个向量内部数据的均值,var命令可以直接计算方差,hist命令可以直接绘出频数直方图,稍作调整就可以绘出频率直方图。
A=[16251920253324232024251715212226152322...
2014161114281813273125241619232617143021...
1816181920221922182626132113111923182428...
1311251517182216131213110915182115121713...
1412161008231811162813212212081521181616...
1928191214192828281321281911151824181628...
1915132214162420281818281413282924281418...
1818082116243216281915181810121626181933...
0811182723112222132814221826181632272524...
1717283316202832192318281524282916171918];
E=mean(A)%计算数据A的样本均值
D=var(A)%计算数据A的样本方差
[a,b]=hist(A);
bar(b,a/sum(a));
%将区间分为10个小区间绘制频率直方图
样本数据'
频率'
title('
频率直方图'
(1)均值与方差计算结果:
计算得这组数据均值为19.5176,方差为34.4025。
(2)频率分布直方图
任务六、利用MATLAB模拟高尔顿板钉实验
高尔顿钉板实验如图,每一黑点表示钉在板上的一颗钉子,它们彼此的距离均相等,上一层的每一颗的水平位置恰好位于下一层的两颗正中间。
从入口处放进一个直径略小于两颗钉子之间的距离的小圆玻璃球,当小圆球向下降落过程中,碰到钉子后皆以1/2的概率向左或向右滚下,于是又碰到下一层钉子。
如此继续下去,直到滚到底板的一个格子内为止。
把许许多多同样大小的小球不断从入口处放下,只要球的数目相当大,它们在底板将堆成近似于正态的密度函数图形。
为了模拟高尔顿钉板实验,可以通过一个向量对钉子进行模拟,数值即为钉子的横坐标,每落下一次即根据随机数判断一次运动方向,如此便可以模拟任意球个数、任意钉板层数的高尔顿钉板实验。
num=100000;
ball=zeros(1,num);
%产生num个球,横坐标为0
floor=21;
%floor层钉板
fori=1:
floor
rd=binornd(1,0.5,1,num);
%产生随机数向量
%利用随机数向量对小球运动方向判断
%随机数为0,小球向左运动;
随机数为1,小球向右运动
forj=1:
num
ifrd(j)==0
ball(j)=ball(j)-0.5;
else
ball(j)=ball(j)+0.5;
end
end
hist(ball,50)
图为100000个小球经过21层钉板后的分布图,可以很明显的看出小球分布近似于正态分布的密度函数。
任务七、利用大量随机数实验模拟会面问题
两人相约中午12时到13时在某地会面,双方约定,先到者必须等候对方15分钟,过了15分钟如果对方仍未到达则离去。
试模拟两人会面的概率
问题分析:
两人到达的时间可以认为是服从区间[0,60]上均匀分布的两个随机变量,两人会面的条件是两人到达时间之差的绝对值在15分钟以内,可以产生大量实验,在随机数的模拟值之差在一个范围内时,认为两人会面成功,计算两人会面的次数,与总实验次数之比即为两人会面概率的模拟值。
n=1000000;
k=0;
x=rand(1,n);
y=rand(1,n);
n
ifabs(x(i)-y(i))<
1/4
k=k+1;
p=k/n
两人会面的概率模拟值为0.4375,与理论计算值非常接近。
心得体会
每次做完一个需要使用MATLAB的软件,我都会更加觉得MATLAB是一个强大的软件,竟然将这些常见的分布都做成了简单方便的命令可以直接调用,因为之前在计算均值、方差时,都是自己写程序来计算,这样似乎浪费了不少时间,达到的效果还不一定好。
因此,学无止境,要想用好MATLAB,还是要多加练习。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 西安 交大 概率论 上机 实验 报告