matlab实现瑞利信道需要的步骤,瑞利信道Matlab仿真程序

%%File_C7:Jakes.m

%本程序将一随机信号通过瑞利信道产生输出

%%

clear;

clc;

Ts=0.02;

fmax=2;%最大多普勒频移

Nt=400;%采样序列的长度

sig=j*ones(1,Nt);%信号

t=[0:Nt];

%设定信道仿真参数

N0=25;

D=1;

[u]=jakes_single_rayleigh(N0,D,fmax,Nt,Ts);%生成瑞利信道 RecSignal=u.*sig;

plot(20*log10(RecSignal));

%JakesRayleigh.m

%本函数用Jakes方法产生单径的符合瑞利分布的复随机过程

%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [u]=jakes_single_rayleigh(N0,D,fmax,M,Ts,Tc) % 输入参数：

% N0 频率不重叠的正弦波个数

% D 方差，可由输入功率得到

% fmax 最大多普勒频移

% M 码片数

%输出参数

%u 输出复信号

%u1 输出信号的实部

%u2 输出信号的虚部

%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%

N=4*N0+2;%Jakes仿真叠加正弦波的总个数

%计算Jakes仿真中的离散多普勒频率fi,n

f=zeros(1,N0+1);

for n=1:N0

f(n)=fmax*cos(2*pi*n/N);

end

f(N0+1)=fmax;

%计算多普勒增益ci,n

%同向分量增益c1,n

c1=zeros(1,N0+1);