ch8example10prog1.m

来自「清华大学出版社 邵玉斌编写的《通信系统建模与仿真实例分析》一书的所有MATLAB」· M 代码 · 共 16 行

% 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法