数学实验6-2概率论与数理统计问题求解.pptVIP

第六讲 概率论与数理统计问题的求解（下）; 8.1概率分布与伪随机数生成 8.1.1 概率密度函数与分布函数概述;通用函数计算概率密度函数值 ;例: 计算正态分布N（0，1）的随机变量X在点0.6578的密度函数值。 解： pdf(norm,0.6578,0,1) ans = 0.3213 例:自由度为8的卡方分布，在点2.18处的密度函数值。 解： pdf(chi2,2.18,8) ans = 0.0363; 随机变量的累积概率值(分布函数值) ;例: 求标准正态分布随机变量X落在区间(-∞，0.4)内的概率。 解： cdf(norm,0.4,0,1) ans = 0.6554 例：求自由度为16的卡方分布随机变量落在[0，6.91]内的概率。 解： cdf(chi2,6.91,16) ans = 0.0250;随机变量的逆累积分布函数 ;例:在标准正态分布表中，若已知F=0.6554,求X 解： icdf(norm,0.6554,0,1) ans = 0.3999 例：公共汽车门的高度是按成年男子与车门顶碰头的机会不超过1%设计的。设男子身高X（单位：cm）服从正态分布N（175，6），求车门的最低高度。 解：设h为车门高度，X为身高。 求满足条件 F{Xh}=0.99，即 F{Xh}=0.01故 h=icdf(norm,0.99, 175, 6) h = 188.9581;8.1.2 常见分布的概率密度函数与分布函数 8.1.2.1 Poisson分布;其中：x为选定的一组横坐标向量， y为x各点处的概率密度函数值。;例：绘制 l =1,2,5,10 时 Poisson 分布的概率密度函数与概率分布函数曲线。 x=[0:15]; y1=[]; y2=[]; lam1=[1,2,5,10]; for i=1:length(lam1) y1=[y1,poisspdf(x,lam1(i))]; y2=[y2,poisscdf(x,lam1(i))]; end plot(x,y1), figure; plot(x,y2) ;8.1.2.2 正态分布;例： x=[-5:.02:5]; y1=[]; y2=[]; mu1=[-1,0,0,0,1]; sig1=[1,0.1,1,10,1]; sig1=sqrt(sig1); for i=1:length(mu1) y1=[y1,normpdf(x,mu1(i),sig1(i))]; y2=[y2,normcdf(x,mu1(i),sig1(i))]; end plot(x,y1), figure; plot(x,y2);8.1.2.3 分布;例： x=[-0.5:.02:5]‘; ％x=[-eps:-0.02:-0.5,0:0.02:5]; x=sort(x’);替代 y1=[]; y2=[]; a1=[1,1,2,1,3]; lam1=[1,0.5,1,2,1]; for i=1:length(a1) y1=[y1,gampdf(x,a1(i),lam1(i))]; y2=[y2,gamcdf(x,a1(i),lam1(i))]; end plot(x,y1), figure; plot(x,y2);8.1.2.4 分布（卡方分布）;例： x=[-eps:-0.02:-0.5,0:0.02:2]; x=sort(x); k1=[1,2,3,4,5]; y1=[]; y2=[]; for i=1:length(k1) y1=[y1,chi2pdf(x,k1(i))]; y2=[y2,chi2cdf(x,k1(i))];end plot(x,y1), figure; plot(x,y2);8.1.2.5 分布;例： x=[-5:0.02:5]; k1=[1,2,5,10]; y1=[]; y2=[]; for i=1:length(k1) y1=[y1,tpdf(x,k1(i))]; y2=[y2,tcdf(x,k1(i))]; end plot(x,y1), figure; plot(x,y2);8.1.2.6 Rayleigh分布;例： x=[-eps:-0.02:-0.5,0:0.02:5]; x=sort(x); b1=[.5,1,3,5]; y1=[]; y2=[]; for i=1:length