《瑞利信道Matlab仿真程序》由会员分享，可在线阅读，更多相关《瑞利信道Matlab仿真程序（2页珍藏版）》请在装配图网上搜索。
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(NO,D,fmax,Nt,Ts);%生成瑞利信道RecSignal=u.*sig;plot(20*log10(RecSignal);%JakesRayleigh.m%本函数用Jakes方法产生单径的符合瑞利分布的复随机过程%functionu=jakes_single_rayleigh(N
2、0,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)=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);end
3、c1(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)*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=0:Ts:(M-1)*Tsk=k+1;u(k)=u1(k)-j*u2(k);end%程序结束
瑞利信道Matlab仿真程序
