gusucode.com > matlab通信工程仿真源码（张德丰等编著）程序书籍 > matlab_code/matlab通信工程仿真源码（张德丰等编著）/第9章/li9_21.m

    clear;
x=[1490 1440 1680 1610 1500 1750 1550 1420 1800 1580];  %样本
m0=1600;  %给定的期望值
n=length(x);  %样本数
xbar=mean(x);  %样本平均
s=std(x,1);  %样本标准差(有偏)
al=0.02;  %显著性水平
% 期望的假设检验
t1=tinv(1-al/2,n-1)   %自由度n-1的t分布al/2分位点
t=(xbar-m0)./(s./sqrt(n-1))   %计算统计量
h_mean=(t>abs(t1))  %判断:若拒绝,则h_mean等于1
%方差的假设检验
sig2=14400;   %给定的方差值
ch_1=chi2inv(al/2,n-1)  %1-al/2分位点
ch_2=chi2inv(1-al/2,n-1)  %al/2分位点
ch=n*s^2/(sig2)   %计算统计量
h_var=(ch<ch_1)|(ch>ch_2)  %判断:若拒绝,则h_var等于1