matlab连续型随机变量的分布分析.docx

matlab连续型随机变量的分布分析.docx

《matlab连续型随机变量的分布分析.docx》由会员分享，可在线阅读，更多相关《matlab连续型随机变量的分布分析.docx（11页珍藏版）》请在冰豆网上搜索。

matlab连续型随机变量的分布分析

连续型随机变量的分布及其数字特征

一、基本概念

设随机变量X的分布函数为F（x），若存在非负函数f（x），使对任意实数x，有

≤

则称X为连续型随机变量，并称f（x）为X的概率密度，它满足以下性质：

1f（x）≥0，-∞＜x＜+∞;②

;③P{a

;④P{x=a}=0．

二、常见的三种连续型随机变量的概率分布

常用的三种连续型随机变量的概率分布是均匀分布、指数分布和正态分布．

（1）均匀分布

若连续型随机变量X的概率密度为

MATLAB提供的有关均匀分布的函数如下：

unifpdf（X,A,B）均匀分布的密度函数

unifcdf（X,A,B）均匀分布的累积分布函数

unifinv（P,A,B）均匀分布的逆累积分布函数

unirnd（A,B,m,n）均匀分布的随机数发生器

unifstat（A,B）均匀分布的数学期望与方差

其中X为随机变量，P为概率值，A，B为均匀分布参数，m和n为生成随机数矩阵的行数和列数．

（2）指数分布

如果随机变量X的概率密度为

其中

为常数，则称X服从参数为

的指数分布，记作X～e（

）．

MATLAB提供的有关指数分布的函数如下：

exppdf（X,L）指数分布的密度函数

expcdf（X,L）指数分布的累积分布函数

expinv（P,L）指数分布的逆累积分布函数

exprnd（X,L,m,n）产生服从指数分布的随机数

expstat（L）求指数分布的数学期望与方差

其中X为随机变量，L为参数

，P为显著概率，m和n为随机数矩阵的行数和列数．

绘制指数分布密度函数和累积分布函数图形的程序如下

x=-0.1:

0.001:

0.4;

y=exppdf（x,0.05）;z=expcdf（x,0.05）;

subplot（1,2,1）;plot（x,y,'k'）;

axis（[-0.1,0.4,-0.1,21]）;

subplot（1,2,2）;plot（x,z,'k'）;

-0.1

0

0.1

0.2

0.3

0

5

10

15

20

-0.1

0

0.1

0.2

0.3

0

0.5

1

axis（[-0.1,0.4,-0.1,1.1]）;

指数分布的密度函数及累积分布函数图

（3）标准正态分布

如果随机变量X的概率密度为：

其中

和

均为常数，且

>0，则称X服从参数为

和

2的正态分布，记作X～N（

2）．当

=0，

=1时，称X服从标准正态分布，记作X～N（0,1）．

MATLAB提供的有关正态分布的函数如下：

normpdf（X,M,C）正态分布的密度函数

normcdf（X,M,C）正态分布的累积分布函数

norminv（P,M,C）正态分布的逆累积分布函数

normrnd（M,C,m,n）产生服从正态分布的随机数

normstat（M,C）求正态分布的数学期望和方差

其中X为随机变量，M为正态分布参数

，C为参数

，P为显著概率，m和n为随机矩阵的行数和列数．

绘制标准正态分布的密度函数及累积分布函数图和一般正态分布的密度函数及累积分布函数图的程序如下：

-4

-2

0

2

4

0

0.2

0.4

-4

-2

0

2

4

0

0.5

1

0

5

10

15

0

0.2

0.4

0.6

0

5

10

15

0

0.5

1

x=-4:

0.01:

4;

y=normpdf（x,0,1）;z=normcdf（x,0,1）;

subplot（2,2,1）;plot（x,y,'k'）;

axis（[-4,4,-0.1,0.5]）;

subplot（2,2,2）;plot（x,z,'k'）;

axis（[-4,4,-0.1,1.1]）;

x=-4:

0.01:

16;

y1=normpdf（x,6,1）;z1=normcdf（x,6,1）;

y2=normpdf（x,6,4）;z2=normcdf（x,6,4）;

y3=normpdf（x,6,0.6）;z3=normcdf（x,6,0.6）;

subplot（2,2,3）;plot（x,y1,'k',x,y2,'k',x,y3,'k'）;

axis（[-4,16,-0.1,0.8]）;

subplot（2,2,4）;plot（x,z1,'k',x,z2,'k',x,z3,'k'）;

axis（[-4,16,-0.1,1.1]）;

三、求解方法

（1）通用函数介绍.

Pdf计算已选函数的概率密度函数,调用格式为:

Y=Pdf（name,X,A）

Y=Pdf（name,X,A,B）

Y=Pdf（name,X,A,B,C）

Name为上表中取stat后的字符，如beta、bino、chiz、exp等。

（2）利用专用函数.

Betapdf（X1，A1，B）

Binopaf（X,N,P）

四、例题

绘制卡方分布密度函数在n分别等于1，5，15时的值

>>x=0:

0.2:

30

y1=chi2pdf（x,1）

plot（x,y1,'+'）

holdon

y2=chi2pdf（x,5）

plot（x,y2,'+'）

y2=chi2pdf（x,15）

plot（x,y2,'o'）

axis（[0,30,0,0.2]）

Columns1through13

00.20000.40000.60000.80001.00001.20001.40001.60001.80002.00002.20002.4000

Columns14through26

2.60002.80003.00003.20003.40003.60003.80004.00004.20004.40004.60004.80005.0000

Columns27through39

5.20005.40005.60005.80006.00006.20006.40006.60006.80007.00007.20007.40007.6000

Columns40through52

7.80008.00008.20008.40008.60008.80009.00009.20009.40009.60009.800010.000010.2000

Columns53through65

10.400010.600010.800011.000011.200011.400011.600011.800012.000012.200012.400012.600012.8000

Columns66through78

13.000013.200013.400013.600013.800014.000014.200014.400014.600014.800015.000015.200015.4000

Columns79through91

15.600015.800016.000016.200016.400016.600016.800017.000017.200017.400017.600017.800018.0000

Columns92through104

18.200018.400018.600018.800019.000019.200019.400019.600019.800020.000020.200020.400020.6000

Columns105through117

20.800021.000021.200021.400021.600021.800022.000022.200022.400022.600022.800023.000023.2000

Columns118through130

23.400023.600023.800024.000024.200024.400024.600024.800025.000025.200025.400025.600025.8000

Columns131through143

26.000026.200026.400026.600026.800027.000027.200027.400027.600027.800028.000028.200028.4000

Columns144through151

28.600028.800029.000029.200029.400029.600029.800030.0000

y1=

Columns1through13

Inf0.80720.51640.38150.29900.24200.19990.16740.14170.12090.10380.08950.0776

Columns14through26

0.06740.05880.05140.04500.03950.03480.03060.02700.02380.02110.01860.01650.0146

Columns27through39

0.01300.01150.01030.00910.00810.00720.00640.00570.00510.00460.00410.00360.0032

Columns40through52

0.00290.00260.00230.00210.00180.00170.00150.00130.00120.00110.00090.00090.0008

Columns53through65

0.00070.00060.00050.00050.00040.00040.00040.00030.00030.00030.00020.00020.0002

Columns66through78

0.00020.00010.00010.00010.00010.00010.00010.00010.00010.00010.00010.00010.0000

Columns79through91

0.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.0000

Columns92through104

0.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.0000

Columns105through117

0.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.0000

Columns118through130

0.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.0000

Columns131through143

0.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.0000

Columns144through151

0.00000.00000.00000.00000.00000.00000.00000.0000

y2=

Columns1through13

00.01080.02750.04580.06380.08070.09590.10940.12090.13060.13840.14440.1489

Columns14through26

0.15190.15360.15420.15370.15230.15010.14730.14400.14020.13600.13150.12690.1220

Columns27through39

0.11710.11210.10720.10220.09730.09250.08780.08320.07870.07440.07020.06620.0623

Columns40through52

0.05860.05510.05170.04850.04550.04260.03990.03730.03490.03260.03040.02830.0264

Columns53through65

0.02460.02290.02130.01980.01840.01710.01590.01480.01370.01270.01180.01090.0101

Columns66through78

0.00940.00870.00800.00740.00690.00640.00590.00540.00500.00460.00430.00390.0036

Columns79through91

0.00340.00310.00290.00260.00240.00220.00210.00190.00170.00160.00150.00140.0013

Columns92through104

0.00120.00110.00100.00090.00080.00080.00070.00060.00060.00050.00050.00050.0004

Columns105through117

0.00040.00040.00030.00030.00030.00020.00020.00020.00020.00020.00020.00010.0001

Columns118through130

0.00010.00010.00010.00010.00010.00010.00010.00010.00010.00010.00010.00000.0000

Columns131through143

0.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.00000.0000

Columns144through151

0.00000.00000.00000.00000.00000.00000.00000.0000

y2=

Columns1through13

00.00000.00000.00000.00000.00000.00000.00000.00000.00010.00010.00020.0003

Columns14through26

0.00040.00060.00080.00110.00150.00200.00260.00330.00410.00500.00600.00720.0085

Columns27through39

0.00990.01140.01310.01490.01680.01880.02090.02310.02540.02770.03020.03260.0351

Columns40through52

0.03760.04010.04260.04510.04750.04990.05230.05460.05680.05890.06100.06290.0647

Columns53through65

0.06650.06810.06950.07090.07210.07320.07420.07500.07570.07630.07670.07700.0772

Columns66through78

0.07720.07720.07700.07670.07630.07580.07530.07460.07380.07300.07200.07100.0700

Columns79through91

0.06890.06770.06650.06520.06390.06250.06120.05980.05840.05690.05550.05400.0526

Columns92through104

0.05110.04970.04820.04680.04530.04390.04250.04110.03970.03840.03700.03570.0344

Columns105through117

0.03320.03200.03070.02960.02840.02730.02620.02520.02410.02310.02220.02120.0203

Columns118through130

0.01940.01860.01780.01700.01620.01550.01480.01410.01340.01280.01220.01160.0110

Columns131through143

0.01050.01000.00950.00900.00860.00810.00770.00730.00700.00660.00630.00590.0056

Columns144through151

0.00530.00500.00480.00450.00430.00400.00380.0036