rayleigh_clarke瑞利衰落的clarke仿真教材135页.m

这是个扩频系统中的PN码仿真程序

function [h]=Rayleigh_Clarke(fm,N,fs);
%*************************************
%参考《无线通信原理与应用》P135
% fm:多谱勒频移
% N :采样点数
% fs:输出线谱的采样频率,为扩频之前的符号速率
%*************************************
%频率间隔
deltaf=2*fm/(2*N-1);
%谱线数
M_Comp=fs/deltaf;
M=2^(ceil(log2(M_Comp)));%ceil数是正无取整

%产生两路高斯复噪声
randn('state',sum(100*clock));
h1=randn(1,N).*exp(j*rand(1,N)*2*pi);
%使左右共轭对称
h1=[conj(fliplr(h1)),h1];

randn('state',sum(100*clock)+500);
h2=randn(1,N).*exp(j*rand(1,N)*2*pi);
h2=[conj(fliplr(h2)),h2];

%生成滤波器幅度特性
n=0:N-2;
Sezf(n+1)=1.5./(pi*fm.*sqrt(1-((n+1/2).*deltaf./fm).^2));
Sezf(N)=2*Sezf(N-1)-Sezf(N-2);
Sezf=[fliplr(Sezf),Sezf];
%高斯噪声通过多谱勒滤波器
h1=sqrt(Sezf).*h1;
h2=sqrt(Sezf).*h2;
%进行IFFT变换
h1=ifft(h1,M);
h2=ifft(h2,M);

%模平方,相加,开方
h=sqrt((abs(h1)).^2+(abs(h2)).^2);
%幅度归一化
h=M*h./sum(h);