前言因?yàn)槎鄰胶团矂优_運(yùn)動等影響要素，使得挪動信道對傳輸信號在時間、頻次和角度上造成了色散，如時間色散、頻次色散、角度色散等等，所以多徑信道的特征對通訊質(zhì)量有著至關(guān)重要的影響，而多徑信道的包絡(luò)統(tǒng)計(jì)特征成為我們研究的焦點(diǎn)。依據(jù)不一樣無線環(huán)境，接收信號包絡(luò)一般聽從幾種典型散布，如瑞利散布、萊斯散布和Nakagami-m散布。在本文中，特意針對聽從瑞利散布的多徑信道進(jìn)行模擬仿真，進(jìn)一步加深對多徑信道特征的認(rèn)識。仿真原理1、瑞利散布簡介環(huán)境條件：往常在離基站較遠(yuǎn)、反射物許多的地域，發(fā)射機(jī)和接收機(jī)之間沒有直射波路徑，存在大批反射波；抵達(dá)接收天線的方向角隨機(jī)且在（0~2π）平均散布；各反射波的幅度和相位都統(tǒng)計(jì)獨(dú)立。幅度、相位的散布特征：包絡(luò)r聽從瑞利散布，θ在0～2π內(nèi)聽從平均散布。瑞利散布的概率散布密度如圖1所示：圖1瑞利散布的概率散布密度2、多徑衰敗信道基本模型依據(jù)ITU-RM.1125標(biāo)準(zhǔn)，失散多徑衰敗信道模型為N(t)y(t)rk(t)x(tk)k1(1)此中,rk(t)復(fù)路徑衰敗，聽從瑞利散布;k是多徑時延。多徑衰敗信道模型框圖如圖2所示：圖2多徑衰敗信道模型框圖3、產(chǎn)生聽從瑞利散布的路徑衰敗r(t)利用窄帶高斯過程的特征，其振幅聽從瑞利散布，即r(t)n(t)2n(t)2cs（2）上式中，nc(t)、ns(t)分別為窄帶高斯過程的同相和正交支路的基帶信號。第一產(chǎn)生獨(dú)立的復(fù)高斯噪聲的樣本，并經(jīng)過FFT后形成頻域的樣本，而后與S（f）開方后的值相乘，以獲取知足多普勒頻譜特征要求的信號，經(jīng)IFFT后變換成時域波形，再經(jīng)過平方，將兩路的信號相加并進(jìn)行開方運(yùn)算后，形成瑞利衰敗的信號r(t)。以下列圖3所示:圖3瑞利衰敗的產(chǎn)生表示圖此中，1.5S(f)（3）ffc)2fm1(fm4、產(chǎn)生多徑延時k多徑/延時參數(shù)如表1所示：表1多徑延時參數(shù)TapRelativedelay(ns)Averagepower(dB)1002310-1.03710-9.041090-10.051730-15.062510-20.0仿真框架依據(jù)多徑衰敗信道模型（見圖2），利用瑞利散布的路徑衰敗r(t)（見圖3）和多徑延時參數(shù)k（見表1），我們能夠獲取多徑信道的仿真框圖，如圖4所示；圖4多徑信道的仿真框圖仿真結(jié)果1、多普勒濾波器的頻響圖5多普勒濾波器的頻響2、多普勒濾波器的統(tǒng)計(jì)特征圖6多普勒濾波器的統(tǒng)計(jì)特征3、信道的時域輸入/輸出波形圖7信道的時域輸入/輸出波形小組分工程序編寫：吳溢升報告撰寫：譚世恒仿真代碼%main.mclc;LengthOfSignal=10240;

%信號長度

(最好大于兩倍

fc)fm=512;%最大多普勒頻移fc=5120;%載波頻次t=1:LengthOfSignal;%SignalInput=sin(t/100);SignalInput=sin(t/100)+cos(t/65);%信號輸入power=[0-1-9-10-15-20];

%dBy_in=[zeros(1,delay(6))SignalInput];

%為時移補(bǔ)零y_out=zeros(1,LengthOfSignal);

%用于信號輸出fori=1:6Rayl;y_out=y_out+r.*y_in(delay(6)+1-delay(i):delay(6)+LengthOfSignal-delay(i))*10^(power(i)/20);end;figure(1);subplot(2,1,1);plot(SignalInput(delay(6)+1:LengthOfSignal));%去除時延造成的空白信號title('SignalInput');subplot(2,1,2);plot(y_out(delay(6)+1:LengthOfSignal));%去除時延造成的空白信號title('SignalOutput');figure(2);subplot(2,1,1);hist(r,256);title('AmplitudeDistributionOfRayleighSignal')subplot(2,1,2);hist(angle(r0));title('AngleDistributionOfRayleighSignal');figure(3);plot(Sf1);title('TheFrequencyResponseofDopplerFilter');%Rayl.mf=1:2*fm-1;

%通頻帶長度y=0.5./((1-((f-fm)/fm).^2).^(1/2))/pi;

%多普勒功率譜

(基帶)Sf=zeros(1,LengthOfSignal);Sf1=y;%多普勒濾波器的頻響Sf(fc-fm+1:fc+fm-1)=y;

%(把基帶映照到載波頻次

)x1=randn(1,LengthOfSignal);x2=randn(1,LengthOfSignal);nc=ifft(fft(x1+i*x2).*sqrt(Sf));

%同相重量x3=randn(1,LengthOfSignal);x4=