ch8example9prog1.m

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

% ch8example9prog1.m
clear;
r(1)=13;    % 随机种子（为奇数）
N=100;      % 样本数
for n=2:N
    r(n)=mod(7^5*r(n-1), 2^31);% 乘同余法产生随机数
end
r=r./2^31; % 随机数归一化
[f,x] = ecdf(r);   % 计算随机数样本r的经验分布
P = unifcdf(x,0,1);% 均匀分布的理论概率分布函数
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);% 经验分布和理论分布曲线对比
[H,P,KSSTAT,CV] = kstest(r,[x,P],0.05)