ch8example10prog1.m

通信程序源代码和模型通信程序源代码和模型通信程序源代码和模型

% ch8example10prog1.m
clear;clc;
N=100;         % 样本数
mu=2; sigma=1; % 参数：均值和标准差
r=normrnd(mu,sigma,N,1); % 随机数产生
[f,x] = ecdf(r);         % 计算随机数样本r的经验分布
muhat=mean(r);           % 为K-S检验估计分布参数
sigmahat=std(r);
P = normcdf(x,muhat,sigmahat);% 正态分布的理论概率分布函数
Kn=sqrt(N)*max(abs(f-P))      % 计算统计量
u_alpha=1.3581                % 显著性水平为0.05时候的分位点
H=(Kn>u_alpha)                % 判决
stairs(x,f);hold on;plot(x,P);% 经验分布和理论分布曲线对比
[Hks,Pks,KSSTAT,CVks] = kstest(r,[x,P],0.05) % K-S法
[Hli,Pli,LSTAT,CVli] = lillietest(r,0.05)    % Lilliefors法
[Hjb,Pjb,JBSTAT,CVjb] = jbtest(r,0.05)       % Jarque-Bera法