資源描述：
《瑞利信道Matlab仿真程序.doc》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、%%File_C7:Jakes.m%本程序?qū)⒁浑S機(jī)信號通過瑞利信道產(chǎn)生輸出%%clear;clc;Ts=0.02;fmax=2;%最大多普勒頻移Nt=400;%采樣序列的長度sig=j*ones(1,Nt);%信號t=[0:Nt];%設(shè)定信道仿真參數(shù)N0=25;D=1;[u]=jakes_single_rayleigh(N0,D,fmax,Nt,Ts);%生成瑞利信道RecSignal=u.*sig;plot(20*log10(RecSignal));%JakesRayleigh.m%本函數(shù)用Jakes方法產(chǎn)生單徑的符合瑞利分布的復(fù)隨
2、機(jī)過程%%%%%%%%%%%%%%%%%%%%%%%%%%%function[u]=jakes_single_rayleigh(N0,D,fmax,M,Ts,Tc)%輸入?yún)?shù)：%N0頻率不重疊的正弦波個數(shù)%D方差，可由輸入功率得到%fmax最大多普勒頻移%M碼片數(shù)%輸出參數(shù)%u輸出復(fù)信號%u1輸出信號的實(shí)部%u2輸出信號的虛部%%%%%%%%%%%%%%%%%%%%%%%%%%%%%N=4*N0+2;%Jakes仿真疊加正弦波的總個數(shù)%計算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);%插入隨機(jī)相移ph_i,解決Jakes方法的廣義平穩(wěn)問題n=(1:N0+1);U=rand(size(n));[x,k]=sort(U);ph_i=2*pi*n(k)/(N0+1);%計算復(fù)包絡(luò)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%程序結(jié)束
