chanmodel_sui3.m

802.16a的三径莱斯信道模型

function fade = chanmodel_sui3()
%   data bit rate = 2M bps, delay = 0.5 us, so every path delays one bit.
%   公式解析见802.16a信道模型E:\毕设相关\陆浩\xieyi\信道模型
N = 10000;
M = 256;
Dop_res = 0.1;
res_accu = 20;
OR = 4;
P = [0 -5 -10];
K = [1 0 0];
tau = [0.0 0.5 1.0];
Dop = [0.4 0.4 0.4];

Fnorm = -1.5113;  %gain normalization factor in dB

P = 10.^(P/10);   % calculate linear power
s2 = P./(K+1);    % calculate variance
m2 = P.*(K./(K+1)); % calculate constant power
m = sqrt(m2);      % calculate constant part
%rmsdel = sqrt(sum(P.*(tau.^2))/sum(P) - (sum(P.*tau)/sum(P))^2);
%fprintf('rms delay spread %6.3f us\n', rmsdel);

% creating the Ricean channel coefficients with the specified powers.
L = length(P);     % number of taps
paths_r = sqrt(1/2)*(randn(L, N) + j*randn(L, N)).*((sqrt(s2))'*ones(1,N));
paths_c = m' * ones(1,N);

% the white spectrum is shaped according to the Doppler PSD function,that
% is to mean the effect of the Doppler.
paths = paths_r + paths_c;
paths = paths * 10^(Fnorm/20);    % multiply all coefficients with F

paths_1 = paths(1,:);
paths_2 = paths(2,:);
paths_3 = paths(3,:);

fade = [paths_1;paths_2;paths_3];