fade.m

OFDM系统仿真

%fade.m
%瑞利衰落仿真程序
function[iout,qout,ramp,rcos,rsin]=fade(idata,qdata,nsamp,tstp,fd,no,counter,flat)
%******************variables*************************
%idata:输入I信道的数据
%qdata:输入Q信道的数据
%iout:输出I信道的数据
%qout:输出Q信道的数据
%ramp:衰落影响后的幅度
%rcos:衰落影响的cosine值
%rsin:衰落影响的sine值
%nsamp:仿真的点数
%tstp:最小判定时间
%fd:最大doppler频率
%no:产生衰落的波数
%counter:衰落计数器
%flat:衰落模式
%(1->平衰落,0->相位和幅度均衰落)
%******************************************************

if fd~=0.0
ac0=sqrt(1.0./(2.0.*(no+1)));
as0=sqrt(1.0./(2.0.*no));
ic0=counter;

pai=3.14159265;
wm=2.0.*pai.*fd;
n=4.*no+2;
ts=tstp;
wmts=wm.*ts;
paino=pai./no;
xc=zeros(1,nsamp);
xs=zeros(1,nsamp);
ic=[1:nsamp]+ic0;

for nn=1:no
  cwn=cos(cos(2.0.*pai.*nn./n).*ic.*wmts);
  xc=xc+cos(paino.*nn).*cwn;
  xs=xs+sin(paino.*nn).*cwn;
end

cwmt=sqrt(2.0).*cos(ic.*wmts);
xc=(2.0.*xc+cwmt).*ac0;
xs=2.0.*xs.*as0;

ramp=sqrt(xc.^2+xs.^2);
rcos=xc./ramp;
rsin=xs./ramp;
if flat==1
iout=sqrt(xc.^2+xs.^2).*idata(1:nsamp);
qout=sqrt(xc.^2+xs.^2).*qdata(1:nsamp);
else
iout=xc.*idata(1:nsamp)-xs.*qdata(1:nsamp);
qout=xs.*idata(1:nsamp)+xc.*qdata(1:nsamp);
end

else
iout=idata;
qout=qdata;
end
%************************end of file***********************************