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