sefade.m

来自「信道编码采用卷积码」· M 代码 · 共 49 行

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