suzuki_fading.m

suzuki衰落M文件

function envelope=suzuki_fading(k_user,n_multi)
%%已知用户和路径计算乘性噪声
%%其中参数multipath_ID=(k_user，n_multi),envelope为（1，7040）表示乘性噪声
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%加载全局变量condition_multipath(N_multi,4,K_user),condition_rayleigh(N_multi,2N-1,K_user)
%%condition_lognormal(N_multi,N,K_user)及多普勒频率系数矩阵;
% load condition_multipath;
global condition_rayleigh;
global condition_lognormal;
global d_n1;
%d_n1=25;
global Ts;%修改了！必须用这个采样频率
%Ts=0.0977*10^(-1);%;采样间隔,为1/8chip
%k_user=5;
%n_multi=3;
%Ts=1/1280/10;
load f_c1;
load f_c2;
load f_c3;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
p=condition_rayleigh(n_multi,:,k_user);
p1=p(1:d_n1);
p2=p((d_n1+1):(2*d_n1-1));
p3=condition_lognormal(n_multi,:,k_user);%分别表示三个随即过程的初相
ff1=2*f_c1(1,:)'*(1:7040)*Ts*pi+p1'*ones(1,7040);
ff2=2*f_c2(1,:)'*(1:7040)*Ts*pi+p2'*ones(1,7040);
ff3=2*f_c3(1,:)'*(1:7040)*Ts*pi+p3'*ones(1,7040);%表示相应的7040个相位;
ff1=cos(ff1);ff2=cos(ff2);ff3=cos(ff3);
ff1=ff1.*(f_c1(2,:)'*ones(1,7040));
ff2=ff2.*(f_c2(2,:)'*ones(1,7040));
ff3=ff3.*(f_c3(2,:)'*ones(1,7040));%家多普勒系数
u1=sum(ff1);
u2=sum(ff2);
u3=sum(ff3);%余弦叠加逼近正态
u3=exp(u3);
u_rayleigh=(u1.^2+u2.^2).^0.5;%瑞利分布
u_rayleigh1=u1+j*u2;
envelope=u_rayleigh.*u3;%信道乘性噪声
%注意：这个地方的影响十分严重，我们在这里暂时不考虑
%相位的影响吧！
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%上面采样点数少,时间短,不易进行噪声估计
%包含相位随机与非随机两种形式
%增加采样点数
%M=1280000;%采样点数
%m=M*Ts;%祯数量
%u1=zeros(1,M);
%u2=zeros(1,M);
%p1=rand(m,d_n1)*2*pi;
%p2=rand(m,d_n1-1)*2*pi;
%p3=rand(m,25)*2*pi;%作为每振随即相位
%for i1=1:d_n1
%    for j1=1:m
%    ff1((j1-1)/Ts+1:j1/Ts)=2*f_c1(1,i1)*(1:1/Ts)*Ts*pi+p1(j1,i1);
%    end;
%    ff1=f_c1(2,i1)*cos(ff1);
%    u1=u1+ff1;
%end;
%for i1=1:d_n1-1
    %ff2=2*f_c2(1,i1)*(1:m)*Ts*pi+p2(i1);
    %ff2=2*f_c2(1,i1)*(1:m)*Ts*pi+rand(1,1)*2*pi;
%    for j1=1:m
%    ff2((j1-1)/Ts+1:j1/Ts)=2*f_c2(1,i1)*(1:1/Ts)*Ts*pi+p2(j1,i1);
%    end;
%    ff2=f_c2(2,i1)*cos(ff2);
%    u2=u2+ff2;
%end;
%u_rayleigh=(u1.^2+u2.^2).^0.5;%瑞利幅值分布
%u_rayleigh1=u1+j*u2;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%u3=zeros(1,M);
%for i1=1:25
%    for j1=1:m
%      ff3((j1-1)/Ts+1:j1/Ts)=2*f_c3(1,i1)*(1:1/Ts)*Ts*pi+p3(j1,i1);
%    end;
%      ff3=f_c3(2,i1)*cos(ff3);
%      u3=u3+ff3;
%end;
%u3=exp(u3);
%envelope=u_rayleigh1.*u3;