瑞利信道Matlab仿真程序

瑞利信道Matlab仿真程序瑞利信道Matlab仿真程序瑞利信道Matlab仿真程序xxx公司瑞利信道Matlab仿真程序文件編號(hào)：文件日期：修訂次數(shù)：第1.0次更改批準(zhǔn)審核制定方案設(shè)計(jì)，管理制度%%File_C7:%本程序?qū)⒁浑S機(jī)信號(hào)通過(guò)瑞利信道產(chǎn)生輸出%%clear;clc;Ts=;fmax=2;%最大多普勒頻移Nt=400;%采樣序列的長(zhǎng)度sig=j*ones(1,Nt);%信號(hào)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));%%本函數(shù)用Jakes方法產(chǎn)生單徑的符合瑞利分布的復(fù)隨機(jī)過(guò)程%%%%%%%%%%%%%%%%%%%%%%%%%%%function[u]=jakes_single_rayleigh(N0,D,fmax,M,Ts,Tc)%輸入?yún)?shù)：%N0頻率不重疊的正弦波個(gè)數(shù)%D方差，可由輸入功率得到%fmax最大多普勒頻移%M碼片數(shù)%輸出參數(shù)%u輸出復(fù)信號(hào)%u1輸出信號(hào)的實(shí)部%u2輸出信號(hào)的虛部%%%%%%%%%%%%%%%%%%%%%%%%%%%%%N=4*N0+2;%Jakes仿真疊加正弦波的總個(gè)數(shù)%計(jì)算Jakes仿真中的離散多普勒頻率fi,nf=zeros(1,N0+1);forn=1:N0f(n)=fmax*cos(2*pi*n/N);endf(N0+1)=fmax;%計(jì)算多普勒增益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)*sin(pi/4);%插入隨機(jī)相移ph_i,解決Jakes方法的廣義平穩(wěn)問(wèn)題n=(1:N0+1);U=rand(size(n));[x,k]=sort(U);ph_i=2*pi*n(k)/(N0+1);%計(jì)算復(fù)包絡(luò)u1=zeros(1,M);%Rc(t)u2=zeros(1,M);%Rs(t)u=zeros(1,M);%R(t)k=0;%計(jì)算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%計(jì)算u(t