2.带你入门matlab数理统计常见分布的概率密度函数（matlab程序）

1.简述

      计算概率分布律及密度函数值

matlab直接提供了通用的计算概率密度函数值的函数，它们是pdf 和namepdf函数，使用方式如下： Y=pdf(‘name’，K，A，B)或者：namepdf (K，A，B) 上述函数表示返回在X=K处、参数为A、B、C的概率值或密度值，对于不同的分布，参数个数是不同；name为分布函数名，使用时需要按照对应分布进行改动。函数名总结如下表：

name的取值    函数说明 ‘beta’ 或 ‘Beta’    Beta分布 ‘bino’ 或 ‘Binomial’    二项分布 ‘chi2’ 或 ‘Chisquare’    卡方分布 ‘exp’ 或 ‘Exponential’    指数分布 ‘f’ 或 ‘F’    F分布 ‘gam’ 或 ‘Gamma’    GAMMA分布 ‘geo’ 或 ‘Geometric’    几何分布 ‘hyge’ 或 ‘Hypergeometric’    超几何分布 ‘logn’ 或 ‘Lognormal’    对数正态分布 ‘nbin’ 或 ‘Negative Binomial’    负二项式分布 ‘ncf’ 或 ‘Noncentral F’    非中心F分布 ‘nct’ 或 ‘Noncentral t’    非中心t分布 ‘ncx2’ 或 ‘Noncentral Chi-square’    非中心卡方分布 ‘norm’ 或 ‘Normal’    正态分布 ‘poiss’ 或 ‘Poisson’    泊松分布 ‘rayl’ 或 ‘Rayleigh’    瑞利分布 ‘t’ 或 ‘T’    T分布 ‘unif’ 或 ‘Uniform’    连续均匀分布 ‘unid’ 或 ‘Discrete Uniform’    离散均匀分布 ‘weib’ 或 ‘Weibull’    Weibull分布  

2.代码及运行结果

%%  二项分布的密度函数 clear all; x=1:20; y=binopdf(x,200,0.06);  figure; plot(x,y,'r*'); title('二项分布（n=200，p=0.06）');

%%  泊松分布的密度函数 clear all; x=1:20; y=poisspdf(x,20);       %泊松分布 figure; plot(x,y,'r+'); title('泊松分布');

%%   几何分布 clear all; x=1:10; y=geopdf(x,0.4);        %几何分布 figure; plot(x,y,'rx'); title('几何分布');

%%  均匀分布（离散） clear all; n=10; x=1:n; y=unidpdf(x,n);             %均匀分布（离散） figure; plot(x,y,'ro'); title('均匀分布（离散）');

%%  均匀分布（连续） clear all; x=-2:0.1:15; y=unifpdf(x,0,6);           %均匀分布（连续）   0到6之间 figure; plot(x,y,'r:'); title('均匀分布（连续）');

%%  指数分布 clear all; x=0:0.1:10; y=exppdf(x,2);      %指数分布 figure; plot(x,y,'r:'); title('指数分布');

%%  正态分布 clear all; x=-5:0.1:5; y1=normpdf(x,0,1);          %标准正态分布 y2=normpdf(x,3,3);          %非标准正态分布 figure; plot(x,y1,x,y2,':'); legend('标准正态分布','非标准正态分布'); x1=-5:0.1:5; y3=normpdf(x1,3,1);         %SIGMA=1 y4=normpdf(x1,3,2);         %SIGMA=2 y5=normpdf(x1,3,3);         %SIGMA=3

figure; plot(x1,y3,'r-',x1,y4,'b:',x1,y5,'k--'); legend('SIGMA=1','SIGMA=2','SIGMA=3'); y6=normpdf(x1,0,2);                   %MU=0 y7=normpdf(x1,2,2);                   %MU=2 y8=normpdf(x1,4,2);                   %MU=4

figure; plot(x1,y6,'r-',x1,y7,'b:',x1,y8,'k--'); legend('MU=0','MU=2','MU=4');

%%  三大抽样分布的概率密度函数 %%  卡方分布 clear all; x=0:0.1:15; y1=chi2pdf(x,2);    %卡方分布n=2 y2=chi2pdf(x,3);    %卡方分布n=3 figure; hold on; plot(x,y1); plot(x,y2,':'); legend('n=2','n=3'); title('卡方分布');

%%  t分布 clear all; x=-5:0.1:5; y1=tpdf(x,2);            %t分布（n=2） y2=tpdf(x,10);           %t分布（n=10）   figure; plot(x,y1,'r:',x,y2,'b-'); legend('n=2','n=10'); title('t分布');

%% F分布 clear all; x=0.1:0.1:5; y=fpdf(x,2,5);             %F分布 figure; plot(x,y,'r:'); title('F分布（m=2,n=5）');