概率统计实验课第一讲.pptVIP

第9章?概率论与数理统计问题的计算机求?解;?§9.1 概率分布与伪随机数生?成;问题：?已知xi的值，要求F(xi),?p(xi); 或者已知某概率分?布fi=F(xi),需要求出x?i的值。 利用MATLAB语言?的统计工具箱中提供的函数可以更?容易、更精确地求出各种值。本节?将介绍几种常用的概率分布形式，?并介绍 MATLAB的求解函数?。;§9.1.2 ?常用分布的概率密度函数与分布函?数;调用格式为? y=poisspdf(x?, λ), F=poisscd?f(x, λ), x=pois?sinv(F, λ) 其中，x?为选定的一组横坐标向量，可以由?x=0:k 语句生成，y为x各?点处的概率密度函数的值。 例9?－1 试分别绘制出λ＝1，2?，5，10时，Poisson分?布的概率密度函数与分布函数曲线?。;解： x=[0:15]?; y1=[];y2=[];?lam1=[1,2,5,10]?; for i=1:le?ngth(lam1) y1=[?y1,poisspdf(x,l?am1(i))]; y2=[y?2,poisscdf(x,la?m1(i))]; end ? plot(x,y1),fi?gure; plot(x,y2?);? ;2、正态分?布;例9－2 试分别绘制出（μ?,σ2）为(-1,1),(0,?0.1),(0,1),(0,1?0),(1,1)时正态分布的概?率密度函数与分布函数曲线。 求?解： x=[-5:0.?02:5];y1=[];y2?=[]; mu1=[-1?,0,0,0,1]; sig1?=[1,0.1,1,10,1]?; sig1=sqrt(sig?1); for i=1?:length(mu1) y?1=[y1,normpdf(x?,mu1(i),sig1(i)?)]; y2=[y2,norm?cdf(x,mu1(i),si?g1(i))]; end ? plot(x,y1),fi?gure; plot(x,y2?);;3、Г分?布;MATLAB语言工具箱提供?了Г分布的概率密度函数、分布函?数及逆概率分布的语言。调用格式?为 y=gampdf(x, α?, λ), F=gamcdf(?x, α, λ), x= ga?minv(x, α, λ) 其?中，x为选定的一组横坐标向量，?y为x各点处的概率密度函数的值?。;例9－3 试分别绘?制出（α, λ)为(1,1)、?(1,0.5)、(2,1)、(?1,2)、(3,1)时Г分布的?概率密度函数与分布函数曲线。 ?求解： x=[-5:0.02:?5];y1=[];y2=[]?; a1=[1,1,2,1,3?]; lam1=[1,0.5,?1,2,1]; for i=1?:length(lam1) y?1=[y1,gampdf(x,? a1(i), lam1(i)?)]; y2=[y2, gam?cdf(x, a1(i),la?m1(i))]; end ? plot(x,y1),fi?gure; plot(x,y2?);? ;4、χ2分?布;例9－4 试分别绘制?出 k为1，2，3，4，5时χ?2分布的概率密度函数与分布函数?曲线。 求解： x=[-eps?:-0.02:-0.05,0:?0.02:2]; x=sort?(x); k1=[1,?2,3,4,5];y1=[];?y2=[]; for i?=1:length(k1) y?1=[y1,chi2pdf(x?, k1(i))]; y2=[?y2,chi2cdf(x, k?1(i))]; end ? plot(x,y1), fi?gure; plot(x,y2?);?