raleighfading.m

来自「Rayleigh Fading的MATLAB程式码 相信对需要通道模拟的会有帮」· M 代码 · 共 29 行

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear; clc;

 N0 = 8;
 Nr = 4 * N0 + 2;
 fm = 83.3;
 fn = fm * cos(2*pi*(1:N0)/Nr);
 phi = pi * (1:N0) / (N0 + 1);
% Ts = 0.01 / fm;

rndt = 1000 * rand(1,1);
t=(0:0.001:500);
xc = sqrt(2) * cos(2*pi*fm*t);
xs = zeros(1,length(xc));
for k=1:N0
    xc = 2 * cos(phi(k)) * cos(2*pi*fn(k)*t) + xc;
    xs = 2 * sin(phi(k)) * cos(2*pi*fn(k)*t) + xs;
end
ag = abs(xc+ sqrt(-1) * xs)/sqrt(2*N0+1);
hist(ag, 100);
ylabel('Number of Hits');
xlabel('Magnitude');
figure(2);

plot(0.001*(1:100),ag(1:100))
xlabel('Time');
ylabel('Magnitude');