rayleigh_chan.m

非常好用

function [data_out,Gain]=rayleigh_chan(Ts,data_in,A,LOS)


    gain_linear=10^(A/10);
    Gain=gain_linear*(randn(1,1)+j*randn(1,1))/2;
    data_temp=data_in*Gain;   
    
    x=data_temp;
    L=20;
    Tao=1e-6;
    N=length(data_in);
    for kk=1:N    
        Mid=round(Tao/Ts);
        Low=Mid-L/2;
        High=Mid+L/2;
        tempii=0;
        for ii=Low:High
            temp=pi*(ii-Tao/Ts+1e-30);
            if ((kk-ii)>0)&&((kk-ii)<N+1)
            tempii=tempii+ x(kk-ii)*sin(temp)/(temp);
            end
        end        
        yy(kk)=tempii;
    end     
    
    if LOS==1
        data_out=data_in+yy;
    else
        data_out=yy;
    end