关于多径信道仿真的问题求大神help

BLACK天堂

这是一个基于改进型Jakes模型的单径的平坦型瑞利衰落信道的程序，我想以此为基础改成多径型的，就是像MATLAB自带的多径信道函数
rayleighchan（fd，fs，tau，pdb）一样，可以引入tau时延变量，但不知道怎么弄，求大神帮忙啊



function [h]=rayleigh(fd,t,tau,pdb)

%该程序利用改进的jakes模型来产生单径的平坦型瑞利衰落信道 %输入变量说明：

%  fd：信道的最大多普勒频移 单位Hz     

%  t :信号的抽样时间序列，抽样间隔单位s   

%  h为输出的瑞利信道函数，是一个时间函数复序列        

N=40;  %假设的入射波数目     

wm=2*pi*fd;   

    M=N/4;  %每象限的入射波数目即振荡器数目

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

Ts=zeros(1,length(t)); %信道函数的虚部

    P_nor=sqrt(1/M);   %归一化功率系

    theta=2*pi*rand(1,1)-pi;  %区别个条路径的均匀分布随机相位     

for n=1:M

        %第i条入射波的入射角  

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

        fi_tc=2*pi*rand(1,1)-pi;    %对每个子载波而言在(-pi,pi)之间均匀分布的随机相位

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

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

        Ts=Ts+2*cos(wm*t*sin(alfa(n))+fi_ts);    %计算冲激响应函数

    end;   

   h= P_nor*(Tc+j*Ts);   %乘归一化功率系数得到传输函数