编辑: 霜天盈月祭 | 2017-09-16 |
y=normpdf(x);
z=normpdf(x,0,2);
>> plot(x,y,x,z),gtext('N(0,1)'),gtext('N(0,2^2)') (请同学通过help gtext了解gtext命令) 练习2: 请同学用仿照练习1画出如下图所示的卡方分布和F分布概率密度函数图形 2)计算正态分布的累积分布函数 语法:Y=normcdf(X,mu sigma) 说明:根据相应的均值mu 和方差sigma计算X中每个值的正态分布的累积分布函数值. 练习:求服从标准正态分布的样本值落在区间[-2,2]上的概率 >> P=normcdf(2)-normcdf(-2) (均值mu =0和方差sigma=1时可缺省) 3)计算正态分布的逆累积分布函数 语法:X=norminv(P,mu sigma) 说明:根据相应的,mu 和sigma计算正态分布中累积分布概率值为P的正态分布对应点.P中的值必须位于[0,1]区间上. 练习1:>> x=norminv(0.5,0,1) (标准正态分布累积分布概率为0.5对应的点) 练习2:找一个区间,使它包含95%的标准正态分布值, >> x=norminv([0.025 0.975],0,1) 4)二项分布均值和方差 语法:[m,v]=binostat (N,P) 说明:返回二项分布的均值m和方差v 练习:>> [m,v]=binostat(500,0.01) (试验500次,每次事件发生概率0.01,返回二项分布的均值m和方差v) 5)生成服从正态分布的随机数 语法:R=normrnd(mu,sigma,m,n) 说明:生成m*n形式的正态分布的随机矩阵. 练习:生成均值u=70,方差v=25,30行1列的正态随机数组. >> R=normrnd(70,25,30,1) ****卡方分布和F分布图命令 >> x=0:0.01:20;
y=chi2pdf(x,5);
z=chi2pdf(x,10);
>> plot(x,y,x,z),gtext('chi2(5)'),gtext('chi2(10)') >> x=0:0.01:3;
y=fpdf(x,10,50);
z=fpdf(x,10,5);
>> plot(x,y,x,z),gtext('F(10,50)'),gtext('F(10,5)')