function [h]=Jakes_gen(v,fc,t,seed)
%函數(shù)[h]=myjakesmodel(v,t,fc,k)，
%該程序利用改進的jakes模型來產(chǎn)生單徑的平坦型瑞利衰落信道
%Yahong R.Zheng and Chengshan Xiao "Improved Models for 
%the Generation of Multiple Uncorrelated Rayleigh Fading Waveforms" 
%IEEE Commu letters, Vol.6, NO.6, JUNE 2002
%輸入變量說明：
%  fc：載波頻率 單位Hz
%  v:移動臺速度  單位m/s     
%  t :信道持續(xù)時間  單位s
%  seed: 產(chǎn)生衰落的隨機數(shù)種子
%  h為輸出的瑞利信道函數(shù)，是一個時間函數(shù)復(fù)序列 
%作者：婁文科   日期：05.3.13
		%產(chǎn)生第k條徑的隨機數(shù)
		rand('state',seed);
    %假設(shè)的入射波數(shù)目
    N=40; 
    %電磁波傳播速度即光速
    c=3*10^8; 
    %最大多普勒頻移
    wm=2*pi*fc*v/c;
    %每象限的入射波數(shù)目即振蕩器數(shù)目
    N0=N/4;
    %信道函數(shù)的實部
    Tc=zeros(1,length(t));
    %信道函數(shù)的虛部
    Ts=zeros(1,length(t));
    %歸一化功率系數(shù)
    P_nor=sqrt(1/N0);
    %區(qū)別個條路徑的均勻分布隨機相位
    sita=2*pi*rand(1,1)-pi;
    for I=1:N0
          %第i條入射波的入射角 
            alfa(I)=(2*pi*I-pi+sita)/N;
            %對每個子載波而言在(-pi,pi)之間均勻分布的隨機相位
            fi_tc=2*pi*rand(1,1)-pi;
            fi_ts=2*pi*rand(1,1)-pi;
            %計算沖激響應(yīng)函數(shù)
            Tc=Tc+cos(cos(alfa(I))*wm*t+fi_tc);
            Ts=Ts+cos(sin(alfa(I))*wm*t+fi_ts);
    end;
    %乘歸一化功率系數(shù)得到傳輸函數(shù)
   h=P_nor*(Tc+j*Ts );