【正文】 。[N,wn]=buttord(wp/pi,ws/pi,Rp,As)。[b,a]=butter(N,wn)。y_out_end =filter(b,a,y_out_end)。%濾波y_out_end =2* y_out_end。%恢復(fù)幅度%原信號的頻譜K=fft(SignalInput)。%DSB調(diào)制后信號的頻譜L=fft(y_in)。%y_out的頻譜（含包絡(luò)）M=fft(y_out)。%最終解調(diào)的頻譜N=fft(y_out_end)。%輸出 figure(1)。subplot(4,2,1)。plot(SignalInput(delay(6)+1:LengthOfSignal))。axis([0,3000,2,2])。title(39。原始輸入信號39。)。subplot(4,2,2)。plot(abs(fftshift(K)))。 axis([4900,5100,0,6000])。title(39。原始輸入信號的頻譜 39。)。subplot(4,2,3)。plot(y_in(delay(6)+1:LengthOfSignal))。axis([0,3000,2,2])。 %去除時延造成的空白信號title(39。 進(jìn)入瑞利信道前，DSB調(diào)制后的信號39。)。subplot(4,2,4)。plot(abs(fftshift(L)))。 axis([3500,6500,0,3000])。title(39。進(jìn)入瑞利信道前，DSB調(diào)制后的信號的頻譜 39。)。subplot(4,2,5)。plot(y_out(delay(6)+1:LengthOfSignal))。axis([0,3000,])。 %去除時延造成的空白信號title(39。經(jīng)瑞利信道后，DSB解調(diào)前的信號39。)。subplot(4,2,6)。plot(abs(fftshift(M)))。axis([3500,6500,0,100])。title(39。經(jīng)瑞利信道后，DSB解調(diào)前的信號的頻譜39。)。subplot(4,2,7)。plot(y_out_end(delay(6)+1:LengthOfSignal))。axis([0,3000,])。%去除時延造成的空白信號title(39。最終解調(diào)后的信號39。)。subplot(4,2,8)。plot(abs(fftshift(N)))。axis([4900,5100,0,200])。title(39。最終解調(diào)后的信號的頻譜 39。)。figure(2)。subplot(3,1,1)。hist(r,256)。%繪制直方圖title(39。瑞利信道的幅度分布39。)subplot(3,1,2)。hist(angle(r0))。title(39。瑞利信道的相位分布39。)。subplot(3,1,3)。plot(Sf1)。title(39。多普勒濾波器的頻響特性39。)。% 參考【1】f=1:2*fm1。 %通頻帶長度y=./((1((ffm)/fm).^2).^(1/2))/pi/fm。 %多普勒功率譜(基帶)圖3Sf=zeros(1,LengthOfSignal)。Sf1=y。Sf(fcfm+1:fc+fm1)=y。%(把基帶映射到載波頻率)x1=randn(1,LengthOfSignal)。x2=randn(1,LengthOfSignal)。nc=ifft(fft(x1+1i*x2).*sqrt(Sf))。 %同相分量%首先產(chǎn)生獨立的復(fù)高斯噪聲的樣本，并經(jīng)過FFT后形成頻域的樣本，%然后與S（f）開方后的值相乘，以獲得滿足多普勒頻譜特性要求的信號，%經(jīng)IFFT后變換成時域波形，再經(jīng)過平方，將兩路的信號相加并進(jìn)行開方運算后，形成瑞利衰落的信號r(t)x3=randn(1,LengthOfSignal)。x4=randn(1,LengthOfSignal)。ns=ifft(fft(x3+1i*x4).*sqrt(Sf))。 %正交分量r0=(real(nc)+1i*real(ns))。 %瑞利信號r=abs(r0)。 %瑞利信號幅值（nc、ns分別為窄帶高斯過程的同相和正交支路的基帶信號） 模擬仿真及結(jié)果分析圖5多普勒濾波器的頻響圖6 多普勒濾波器的統(tǒng)計特性圖7 信道的時域/頻域輸入/輸出波形次實驗主要是通過MATLAB仿真瑞利衰落信道的傳輸過程，通過雙邊帶調(diào)幅的調(diào)制與解調(diào)實現(xiàn)信號的傳輸。正如右圖所示：圖中第一、第二個波形是在進(jìn)入瑞利衰落前，第三、第四個波形是在進(jìn)入瑞利衰落后，有明顯的噪聲的存在。由第一個圖輸入，第四個圖輸出，信號的傳遞在存在干擾的情況下基本實現(xiàn)。第一個波形到第二個波形是實現(xiàn)了抑制載波的雙邊帶調(diào)幅；第三個波形到第四個波形是運用同步解調(diào)或相干解調(diào)實現(xiàn)對載波信號的解調(diào)功能。 圖8 信道的時域輸入/輸出波形圖9 信道的頻域變化分析圖9中的第二、第四幅圖可以發(fā)現(xiàn)輸出信號的頻譜圖上一段頻率內(nèi)出現(xiàn)了多余的小頻率，但總體頻率沒有很大變化。說明信道對輸入波形存在影響，但本題中，輸出的波形還是可以基本反映輸入情況的。如圖6中顯示出了瑞利信道的幅度分布和相位分布情況，而在簡介中提到包絡(luò) r 服從瑞利分布，θ在0～2π內(nèi)服從均勻分布。瑞利分布的概率分布密度如圖1，經(jīng)過分析對照可以發(fā)現(xiàn)，輸出的多普勒濾波器的統(tǒng)計特性完全符合要求。 小結(jié)與體會這次實習(xí)是自由定題，當(dāng)聽到這個時讓我產(chǎn)生一種不知所措感，因為，我從沒有真正講學(xué)習(xí)的東西運用到稍微大一點的實際中。不過，越是這樣才能讓我們鍛煉的更多，在用程序編出想要的波形時，我們不僅能體會成功的快樂，而且對通信的傳輸原理有了更加清晰的認(rèn)識，這次設(shè)計要求我們學(xué)會綜合運用我們所學(xué)的通信知識及用Matlab實現(xiàn)，可以上網(wǎng)查資料，也可以圖書館里找。設(shè)計一開始的時候我充滿了恐懼感，面對好多的書，好多的代碼，好混亂的腦中知識……總之，一萬個不想做，但是想到如果不實踐是不會有進(jìn)步的，只得硬著頭皮來，然而隨著一點一點進(jìn)行，發(fā)現(xiàn)一切都不像想象的那么艱巨，慢慢的我完成了我的課程設(shè)計——瑞利衰落信道的仿真。其實做什么事都是那樣，不開始，永遠(yuǎn)害怕打破那層膜，只用勇敢的邁出第一步，一切問題都會在面臨的時候解決。這次設(shè)計最深刻的感受就是：只要開始，一切都可以變得簡單。 參考文獻(xiàn)[1]胡宴如，耿蘇燕主編，高頻電子線路，北京：高等教育出版社，2009。[2]沈衛(wèi)康，宋宇飛，宋紅梅，數(shù)字信號處理，北京：清華大學(xué)出版社，2011。[3]樊昌信，曹麗娜編著，通信原理，北京：國防工業(yè)出版社，2010。[4]潘子宇，Matlab通信仿真設(shè)計指導(dǎo)書，南京工程學(xué)院，2011。10