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