clarke.m

这里有2组图

% clarke.m
% Clarke的一个参考模型
hold off;
clear;                                  % 清理工作现场
N=1000;                                 % 多径数目
tstep=0.00005;                          % 时间步长
t=0.0001:tstep:0.25;                    % 时间取值区间
fc=900*10^6;                            % 载波频率
wc=2*pi*fc;                             % 载波角频率
v=120/3600;	                            % 接收机速度(km/h)   
c=3*10^5;                               % 光速
wm=wc*(v/c);                            % 最大多普勒角频移
fm=wm/(2*pi);                           % 最大多普勒频移
R=unifrnd(0,2*pi,1,N);                  % 1×N阶的[0,2*pi]内的随机数

for n=1:N
    wn(n)=wm*cos(2*pi*n/N);             % 多径波的多普勒角频移
end
for m=1:length(t)                       % m个时间取样点
    Tsum(m)=0;
    for n=1:N                           % N个多径
        Tc(n)=cos(wn(n)*t(m)+R(n));     % 同相分量Tc
        Ts(n)=sin(wn(n)*t(m)+R(n));     % 正交分量Ts
        T(n)=Tc(n)-j*Ts(n);             % 复包络 
        Tsum(m)=Tsum(m)+T(n);           % 某个时间点上的多径相加
    end
end

rt=abs(Tsum)*sqrt(2/N);                 % 计算E场包络值

hist(rt,300);                           % 统计分布状况
title('用Clarke法仿真Rayleigh衰落- hist统计 ');
pause;

p=hist(rt,300);                         % hist的包络
plot(p);                                % 绘制E场分布图
title('用Clarke法仿真Rayleigh衰落-用hist的包络统计出的概率密度');
ylabel('概率密度p');
pause;

rtdB=10*log10(rt);
plot(rtdB,'r');
title('用Clarke法得到的Rayleigh衰落仿真，接收机速度速率120km/hr');
ylabel('接收E场包络(dB)');
xlabel('time(0～250ms分成5000组,步长0.05ms)');
grid on;