瑞利信道的生成的具体问题

wm=2*pi*fd; 

%每象限的入射波数目即振荡器数目 

N0=N/4; 

%信道函数的实部 

Tc=zeros(1,length(t)); 

%信道函数的虚部 

Ts=zeros(1,length(t)); 

%归一化功率系数 

P_nor=sqrt(1/N0); 

%区别个条路径的均匀分布随机相位 

theta=2*pi*rand(1,1)-pi; 

for ii=1:N0 

%第i条入射波的入射角 

alfa(ii)=(2*pi*ii-pi+theta)/N; 

%对每个子载波而言在(-pi,pi)之间均匀分布的随机相位 

fi_tc=2*pi*rand(1,1)-pi; 

fi_ts=2*pi*rand(1,1)-pi; 

%计算冲激响应函数 

Tc=Tc+cos(cos(alfa(ii))*wm*t+fi_tc); 

Ts=Ts+cos(sin(alfa(ii))*wm*t+fi_ts); 

end; 

%乘归一化功率系数得到传输函数 

h=P_nor*(Tc+j*Ts ); 

这段您写的代码中，Tc与Ts的计算公式是什么呢？为什么cos里面还有一个cos？wm为什么等于2pi×多普勒频移呢？能否请您说说这段代码的理论基础？