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