概率论实验

实验一随机变量的概率分布。

一实验目的。

1. 掌握计算随机变量分布律或概率密度值的matlab命令；

2. 掌握计算分布函数的matlab命令；

3. 学习常见分布的随机变量的模拟与应用。

二实验背景知识介绍。

1. 随机变量及其概率分布。

随机变量是定义在样本空间上的实函数，按其取值情况常见有两类：离散型与连续型。

设是随机变量，给定任意实数，记。

则称函数为随机变量的概率分布函数，简称分布函数。分布函数能完整地描述随机变量的统计规律性。

若已知随机变量的分布函数为，则对于任意的实数有。

若为连续型随机变量，是的分布函数，则存在非负函数，对任意实数，有。

称为的概率密度函数或密度函数。

在概率与统计中，常用的分布有：二项分布、几何分布、泊松分布、正态分布、指数分布、均匀分布、分布、分布、分布等。

2. 统计工具箱与常见命令介绍。

为了便于研究概率与统计的计算问题，matlab提供了专门的统计工具箱(stastoolbox),其概率计算的主要功能有：计算相应分布的概率、分布函数、逆分布函数和产生相应分布的随机数。工具箱的统计计算主要功能有：

统计量的数字特征、统计图形的绘制、参数估计、假设检验、方差分析等。

表1：常见分布名称。

在统计工具箱中，matlab为每一种分布提供了5类命令函数，其命令字符分别为：pdf表示概率密度；cdf表示概率分布函数（累积概率）；inv表示逆概率分布函数；stat表示均值与方差；rnd表示生成相应分布的随机数。这样，当需要一种分布的某一类命令函数时，只要将表6.

1中的分布名字符后缀命令函数字符并输入命令参数即可。如，binopdf(x,n,p)表示计算服从参数为二项分布的随机变量在的概率；normcdf(x,) 表示计算服从参数为的正态分布的随机变量在的分布函数；expstas()表示计算服从参数为的指数分布的随机变量的期望与方差，等等。

三实验内容。

1. 常见分布的概率计算。

例1 掷硬币100次，其中正面出现的概率为0.5，这100次中正面出现的次数记为：

1）试计算的概率和的概率；

2）绘制分布函数图形和概率分布律图形。

分析这是100重的贝努力试验，～，计算二项分布的分布函数的命令为binocdf(x,n,p),计算概率分布律的命令为binopdf(x,n,p)。

解设计程序并保存为m文件“

p2=binocdf(45,100,0.5)

x=1:100;

p=binocdf(x,100,0.5);

px=binopdf(x,100,0.5);

subplot(1,2,1)

plot(x,p,'r*')

subplot(1,2,2)

plot(x,px)

程序运行结果：p1 = 0.0485，p2 = 0.1841。输出图形为。

分布函数分布律。

例2 设～，1）当时，求，，；

2）当时，若求；

3）分别绘制时的概率密度函数图形。

分析本题是关于正态分布的有关概率计算问题，只要调用正态分布(norm)的有关命令就能实现其计算。这些命令分别是分布函数命令normcdf()；概率密度命令normpdf()；逆分布函数命令norminv()。

解设计程序并保存为m文件“

p1=normcdf(2.9,2,0.5)-normcdf(1.8,2,0.5)

p2=1-normcdf(-3,2,0.5)

p3=1+normcdf(0.5,2,0.5)-normcdf(3.5,2,0.5)

px=0.95;

x0=norminv(px,2,0.5)

x=-2:0.05:5;

y1=normpdf(x,2,0.2);

y2=normpdf(x,2,0.5);

y3=normpdf(x,2,0.9);

plot(x,y1,x,y2,x,y3)

gtext('标准差为0.2')

gtext('标准差为0.5')

gtext('标准差为0.9')

title('正态概率密度曲线')

程序运行结果为p1 = 0.6195，p2 = 1，3 = 0.0027，x0 = 2.8224。

2 常见分布随机数的产生与计算机模拟。

在matlab的统计工具箱中，提供了产生满足常用分布的随机数命令，这些命令的名字构成为：分布名字符+rnd.如unifrnd(a,b),表示产生均匀分布的随机数。

的随机数还可以用命令rand产生，如rand(m,n), 表示产生[0,1]区间上的维的随机数，还可以通过变换将[0,1]区间的随机数转化为任意区间的随机数。

例3 投掷硬币的计算机模拟。投掷硬币1000次，试模拟掷硬币的结果。

分析设硬币是质地均匀的，掷硬币一次的结果为正面与反面。用命令unifrnd(0,1)或rand(1,n)产生[0,1]区间均匀分布的随机数，将区间二等分。若，则对应于硬币的正面，否则，对应于硬币的反面，掷硬币的结果可以是统计正面与反面出现的频数与频率。

解设计程序并保存为m文件“glst1201li03”

n=1000;

m=0;k=0;

for i=1:n

x(i)=unifrnd(0,1);

if x(i)<=0.5

m=m+1;

else k=k+1;

endend

p1=m/n

p2=k/n

运行结果为 p1 = 0.5370，p2 = 0.4630。

例4个人中至少有两人生日相同的概率是多少？通过计算机模拟此结果。

假设每个人在一年（按365天计算）内每一天出生的可能性都相同，现在随机的选取个人，则他们中至少有两人生日相同的概率为。

经计算，得到表1：

表1个人中至少有两人生日相同的概率。

由表1可以看出，在40人左右的人群里，两人或两人以上生日相同这一事件十有**会发生。

下面通过matlab软件模拟这个结果。

实验过程先在区间上产生个随机数，用365乘以这个随机数，这相当于是在区间上产生个随机数，再给这个随机数取整，这相当于是在区间上随机抽取的个人的生日，通过对这个数的比较，来决定是否有人生日相同。下面的matlab程序，模拟了这个结果。

n=10000; %实验的次数。

for k=10:5:50 %分别抽取的人数。

s=0;w=0;

for i=1:n

x=sort(fix(365*rand(1,k)))产生个不同人的生日。

for j=1:(k-1)

y=x(j+1)-x(j); 比较两人的生日。

if y==0

w=1; %作出两人生日相同的判断。

break;

elsew=0;

endend

s=w+s; %统计此实验中生日相同的次数。

endp=s/n; %计算生日相同的频率。

fprintf('k=%.0f,p=%.5f',k,p)

end输出结果。

k=10,p=0.11694

k=15,p=0.25287

k=20,p=0.41044

k=25,p=0.56922

k=30,p=0.70843

k=35,p=0.81595

k=40,p=0.89100

k=45,p=0.94112

k=50,p=0.97008

与表1相比，接近程度相当好。我们还画出了随着人数的变化，频率的变化曲线图。

图3 个人的生日相同的频率图。

3 如何产生具有的分布律的随机数？

离散型随机数的模拟：设随机变量的分布律为，

令，将作为区间的分点。

若随机变量～，有。令。则有。

据此，可产生具有的分布律的随机数的具体执行过程为：每产生一个区间上均匀分布的随机数，如果，则令的取值为。

例如离散型随机变量的分布律为。

设是上均匀分布的随机数，令。

则是具有的分布律的随机数。

例5.一列火车从甲站经过乙站开往丙站，某人每天赶往乙站乘这趟火车。已知火车从甲站到乙站运行时间为均值30分钟，标准差为2分钟的正态随机变量。

火车大约下午1点离开甲站，火车离开时刻的频率分布和这个人到达乙站的频率分布如表2所示，问他能赶上火车的概率有多大？

表2 火车出发时刻和人到达乙站的时刻及频率分布。

这是一个概率问题，但如果用概率论的方法来求解，将是十分困难的。因为它涉及此人到达时刻、火车离开甲站的时刻、火车运行时间这几个随机变量。计算时应将此人赶火车事件分解为以下几个事件：：

此人1：28分到达火车站并赶上火车；：此人1：

30分到达火车站并赶上火车；：此人1：32分到达火车站并赶上火车；：

此人1：34分到达火车站并赶上火车。。为了计算事件的概率，需要将火车出发的时间进行分解：

火车1：00从甲站出发且在1：28分以后到达乙站；火车1：

05从甲站出发且在1：28分以后到达乙站；火车1：10从甲站出发且在1：

28分以后到达乙站。其余的也可类似进行分解。又因为火车运行时间是随机的，没有有效的方法来进行概率的计算。

下面我们用随机模拟的方法来解决此问题。

用随机模拟的方法解决此问题就是通过多次随机试验来得到此人能赶上火车的概率有多大。为此，首先引入以下随机变量：：火车从甲站开出的时刻；火车从甲站运行到乙站所需要的时间；此人到达乙站的时刻。

此人能赶上火车的条件是。取下午1：00为时间0点，则取值为0，5，10，其概率分别为0.

7，0.2，0.1；服从；取值为28，30，32，34，其概率分别为0.

3，0.4，0.2，0.1.

概率论实验

概率论实验

概率论实验

概率论实验

其他用户还读了