% sefade.m
% 多徑衰落信道仿真程序	 	

function[iout,qout,ramp,rcos,rsin]=sefade(idata,qdata,itau,dlvl,n0,itn,n1,nsamp,tstp,fd,flat)

%****************** variables *************************
% idata  input Ich data     
% qdata  input Qch data     
% iout   output Ich data
% qout   output Qch data
% ramp   : 衰落影響后的幅值
% rcos   : 衰落影響的cosine值
% rsin   : 衰落影響的sine值
% itau   : 每徑的延遲時間
% dlvl   : 每徑的相對衰落能量
% itn    : 每徑的衰落計數器
% n1     ：衰落和延遲總波數
% nsamp  : 仿真的點數      
% tstp   : 最小判定時間                    
% fd     : 最大多普勒頻率              
% no     : 每徑產生衰落的波數   
% counter  : 衰落計數器                         
% flat     : 衰落模式
% (1->flat (only amplitude is fluctuated 平衰落),0->nomal(phase and amplitude are fluctutated 相位和幅度均衰落 
%******************************************************

iout = zeros(1,nsamp);
qout = zeros(1,nsamp);

total_attn = sum(10 .^( -1.0 .* dlvl ./ 10.0));

for k = 1 : n1 

	atts = 10.^( -0.05 .* dlvl(k));

	if dlvl(k) >= 40.0 
	       atts = 0.0;
	end
	

	[itmp,qtmp] = delay ( idata , qdata , nsamp , itau(k));
	[itmp3,qtmp3,ramp,rcos,rsin] = fade (itmp,qtmp,nsamp,tstp,fd,n0(k),itn(k),flat);
	
  iout = iout + atts .* itmp3 ./ sqrt(total_attn);
  qout = qout + atts .* qtmp3 ./ sqrt(total_attn);

end
% ************************end of file***********************************