单径的平坦型瑞利衰落信道rayleigh()的matlab函数实现_matlab rayleigh函数

        在通信系统中有时我们不仅需要用到高斯信道，有时也会用到瑞利信道，但是matlab并没有内置相应的函数块，只有关于瑞利分布函数Rayleigh()，所以接下来，我们来看一看rayleigh()具体在matlab中如何实现。

        在我们需要使用的时候，进行该函数的调用即可。

function [h]=rayleigh(fd,t)
%输入变量说明：
%  fd：信道的最大多普勒频移 单位Hz     
%  t :信号的抽样时间序列，抽样间隔单位s  
%  h为输出的瑞利信道函数，是一个时间函数复序列 
 
    %假设的入射波数目
    N=40; 
    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 );

注：思路来源于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