11.txt

11.txt为产生信道冲激响应的源码

%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Channel impulse response %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%

if (Rayleigh)
    N=50;
    fm=100;
    B=20e3;
    fd=(rand(1,N)-0.5)*2*fm;
    theta=randn(1,N)*2*pi;
    c=randn(1,N);
    c=c/sum(c.^2);
    t=0:fm/B:10000*fm/B;
    Tc=zeros(size(t));
    
    Ts=zeros(size(t));
    for k=1:N
       Tc=c(k)*cos(2*pi*fd(k)*t+theta(k))+Tc;
       Ts=c(k)*sin(2*pi*fd(k)*t+theta(k))+Ts;
    end
    r=ones(Mt*Mr,1)*(Tc.^2+Ts.^2).^0.5;
    index=floor(rand(Mt*Mr,DS)*5000+1);
end