ch8example9prog1.m

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

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