概率统计的MATLAB求解详解.ppt

第10讲 MATLAB求解概率统计问题 一 随机变量及其分布 命令1：Fx=hygecdf(x,M,N,K) 功能：计算超几何分布的累积概率，总共M件产品，其中次品N 件，抽取K件检查，计算发现次品不多于x件的概率Fx=P{次品数X≤x}=F(x) 命令2：x=hygeinv(p,M, N,K) 功能：在已知参数M、N 、 K和p的情况下计算随机量x，使得p=P{0≤次品数X≤x} 命令3：X=hygernd(M,N,K,m,n) 功能：在已知参数M,N ,K的情况下产生m*n维符合超几何分布的随机数矩阵X 命令4：Px=hygepdf(x,M, N, K) 功能：总共M件产品，其中次品N 件，抽取K件检查，计算发现恰好x件次品的概率Px=P{X=x} 命令5：stairs(x,Px) 功能：绘制以 x为横坐标，Px为纵坐标的阶梯平面图；当Px是分布列(或密度)时，绘制概率密度分布图；当Px是累积分布时，绘制概率分布函数图 注：以后碰到命令末尾为： rnd----产生随机数X； cdf----产生分布函数F(x) pdf----产生密度函数p(x)或分布列Px=P{X=x} inv----计算x=F-1(p)→ p=F (x) 命令1：Fx=binocdf(x,n,p) 功能：计算二项分布的累积概率Fx=P{X≤x}=F(x) 命令2：x=binoinv(y, n,p) 功能：计算随机量x，使得y=P{X≤x} 命令3：X=binornd(n,p,M,N) 功能：产生M*N维符合二项分布的随机数矩阵X 命令4：Px=binopdf(x,n, p) 功能：计算试验中事件恰好发生x次的概率Px=P{X=x} 命令1：Fx=poisscdf(x,lambda) 功能：计算累积概率Fx=P{X≤x}=F(x) 命令2：x=poissinv(p, lambda) 功能：计算随机量x，使得p=P{X≤x} 命令3：X=poissrnd(lambda,M,N) 功能：产生M*N维随机数矩阵X 命令4：Px=poisspdf(x,lambda) 功能：计算概率Px=P{X=x} 命令1：Fx=normcdf(x, mu,sigma) 功能：计算累积概率Fx=P{X≤x}=F(x) 命令2：x=norminv(p, mu,sigma) 功能：计算随机量x，使得p=P{X≤x} 命令3：X=normrnd(mu,sigma,M,N) 功能：产生M*N维随机数矩阵X 命令4：Px=normpdf(x, mu,sigma) 功能：计算分布密度p(x)在x的值 补充：randn()---标准正态分布随机数 命令1：Fx=expcdf(x, lambda) 功能：计算累积概率Fx=P{X≤x}=F(x) 命令2：x=expinv(p, lambda) 功能：计算随机量x，使得p=P{X≤x} 命令3：X=exprnd(lambda,M,N) 功能：产生M*N维随机数矩阵X 命令4：Px=exppdf(x, lambda) 功能：计算分布密度p(x)在x的值 命令1：Fx=unifcdf(x, a,b) 功能：计算累积概率Fx=P{X≤x}=F(x) 命令2：x=unifinv(p, a,b) 功能：计算随机量x，使得p=P{X≤x} 命令3：X=unifrnd(a,b,M,N) 功能：产生M*N维随机数矩阵X 命令4：Px=unifpdf(x, a,b) 功能：计算分布密度p(x)在x的值 补充：rand()---(0,1)均匀分布随机数 命令:gamcdf(x, a, lambda), gaminv(p, a, lambda) gampdf(x, a,lambda), gamrnd(a, lambda,m,n) 命令:chi2cdf(x, k), chi2inv(p, k),chi2pdf(x, k) chi2rnd(k,m,n) 命令:tcdf(x, k), tinv(p, k),tpdf(x, k) trnd(k,m,n) 命令:fcdf(x, p,q), finv(F,p,q),fpdf(x, p,q) frnd(p,q,m,n) 程序：》clear; px=binopdf(45,100,0.5) % 计算x=45的概率 px = 0.0485 fx=binocdf(45,100,0.5) % 计算x≤45的概率 fx =0.1841 》x=1:100;p1=binocdf(x,100,0.5);plot(x,p1,+); title(分布函数图) 程序:(1)p=normcdf(2.5,2,0.5)- normcdf(1,2,0.5) p = 0.8186 (2) x=0