数理统计MATLAB求解.pptVIP

概率统计的MATLAB求解 随机变量及其分布 随机变量函数的分布 随机变量的数字特征 参数估计 1 随机变量及其分布 注：以后碰到命令末尾为： rnd----产生随机数X； cdf----产生分布函数F(x) pdf----产生密度函数p(x)或分布列Px=P{X=x} inv----计算x=F-1(p)→ p=F (x) 常见分布的随机数产生 专用函数计算概率密度函数表 专用函数的累积概率值函数表 常用临界值函数表 常见分布的均值和方差 3.1 随机变量及其分布 程序：》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(分布函数图) 3.1 随机变量及其分布 3.1 随机变量及其分布 程序:(1)p=normcdf(2.5,2,0.5)- normcdf(1,2,0.5) p = 0.8186 (2) x=0:0.1:4;px=normpdf(x,2,0.5); fx= normcdf(x,2,0.5); plot(x,px,+b);hold on; plot(x,fx,*r);legend(正态分布函数,正态分布密度); (3) specs=[1.5,1.9]; pp=normspec(specs,2,0.5) 3.1 随机变量及其分布 3.2 随机变量函数的分布 根据概率统计教材中的定理：如果已知随机变量X的密度fX(x),随机变量函数Y=g(X)单调，则Y的密度函数为： fY(x)= fX(h(y))|h(y)|,其中x=h(y)是y=g(x)的反函数。 如果y=g(x)不单调，则将定义域分成若干单调区间进行讨论。也可利用： 据此意思，计算随机变量函数的分布相当于编程 程序：clear; x=solve(y=exp(x)) x=log(y) dy=diff(x,y) dy=1/y fy= 1*abs(dy) fy=1/|y| 注：取值区域需要自己确定，用积分求法作为练习！ 3.3 随机变量的数字特征 3.3 随机变量的数字特征 程序：clear; x=[-1,0,2,3]; p=[1/8,1/4,3/8,1/4]; EX=sum(x.*p) 1.3750 3.3 随机变量的数字特征 程序：clear;syms a b x;fx=a+b*x^2; EX=int(x*fx,x,0,1) EX=1/4*b+1/2*a F=int(fx,x,0,1) F=a+1/3*b f1=EX-3/5;f2=f-1; [a,b]=solve(f1,f2) a=3/5,b=6/5 3.3 随机变量的数字特征 程序：clear;syms x; fx1=0.5*exp(x); fx2=0.5*exp(-x); EY=int(-x*fx1,x,-inf,0) + int(x*fx2,x,0, inf) EY= 1 3.3 随机变量的数字特征 3.3 随机变量的数字特征 程序：clear;syms x;fx=2/pi*(cos(x))^2; EX=int(x*fx,x,-pi/2,pi/2) E2X=int(x^2*fx,x,-pi/2,pi/2) DX=E2X-EX^2 3.3 随机变量的数字特征 3.3 随机变量的数字特征 3.3 随机变量的数字特征 3.4 参数估计 3.4 参数估计 3.4 参数估计 程序：clear; w=normrnd(15,2.5,50,1); 或w=15+2.5*randn(50,1); alpha=0.03; [mh,sh,mc,sc]=normfit(w,alpha) 运行一次：mh=15.1076