瑞利信道Matlab仿真程序.doc

《瑞利信道Matlab仿真程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、%%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方法产生单径的符合瑞利分布的复随

2、机过程%%%%%%%%%%%%%%%%%%%%%%%%%%%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,nf=zeros(1,N0+1);forn=1:N0f(n

3、)=fmax*cos(2*pi*n/N);endf(N0+1)=fmax;%计算多普勒增益ci,n%同向分量增益c1,nc1=zeros(1,N0+1);forn=1:N0c1(n)=D*(2/sqrt(N))*2*cos(pi*n/N0);endc1(N0+1)=D*(2/sqrt(N))*sqrt(2)*cos(pi/4);%正交分量增益c2,nc2=zeros(1,N0+1);forn=1:N0c2(n)=D*(2/sqrt(N))*2*sin(pi*n/N0);endc2(N0+1)=D*(2/sqrt(N))*sqrt(2)

4、*sin(pi/4);%插入随机相移ph_i,解决Jakes方法的广义平稳问题n=(1:N0+1);U=rand(size(n));[x,k]=sort(U);ph_i=2*pi*n(k)/(N0+1);%计算复包络u1=zeros(1,M);%Rc(t)u2=zeros(1,M);%Rs(t)u=zeros(1,M);%R(t)k=0;%计算Rc(t)k=0;fort=0:Ts:(M-1)*Ts;w2=cos(2*pi*f*t+ph_i);ut2=c2*w2.';k=k+1;u2(k)=ut2;end%计算u(t)k=0;fort=

5、0:Ts:(M-1)*Tsk=k+1;u(k)=u1(k)-j*u2(k);end%程序结束