課程設(shè)計(jì)報(bào)告(psk調(diào)制解調(diào))要點(diǎn)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上長(zhǎng)沙理工大學(xué)通信原理課程設(shè)計(jì)報(bào)告 李柳 學(xué) 院 計(jì)算機(jī)與通信工程 專(zhuān) 業(yè) 通信工程 班 級(jí) 通信1203 學(xué) 號(hào) 8 學(xué)生姓名 李柳 指導(dǎo)教師 吳志敏 課程成績(jī) 完成日期 2015年1 月 4日課程設(shè)計(jì)成績(jī)?cè)u(píng)定學(xué) 院 計(jì)算機(jī)與通信工程 專(zhuān) 業(yè) 通信工程 班 級(jí) 通信1203 學(xué) 號(hào) 8 學(xué)生姓名 李柳 指導(dǎo)教師 吳志敏 課程成績(jī) 完成日期2015年1 月 4日指導(dǎo)教師對(duì)學(xué)生在課程設(shè)計(jì)中的評(píng)價(jià)評(píng)分項(xiàng)目?jī)?yōu)良中及格不及格課程設(shè)計(jì)中的創(chuàng)造性成果學(xué)生掌握課程內(nèi)容的程度課程設(shè)計(jì)完成情況課程設(shè)計(jì)動(dòng)手能力文字表達(dá)學(xué)習(xí)態(tài)度規(guī)范要求課程設(shè)計(jì)論文的質(zhì)量指導(dǎo)教師對(duì)課程設(shè)計(jì)的評(píng)定意見(jiàn)綜合成績(jī)

2、指導(dǎo)教師簽字 年 月 日課程設(shè)計(jì)任務(wù)書(shū)計(jì)算機(jī)與通信工程學(xué)院 通信工程專(zhuān)業(yè) 課程名稱(chēng)通信原理課程設(shè)計(jì)時(shí)間2014/2015學(xué)年第一學(xué)期1819周學(xué)生姓名李柳指導(dǎo)老師吳志敏題 目基于MATLAB的PSK調(diào)制解調(diào)實(shí)現(xiàn)主要內(nèi)容：利用MATLAB集成環(huán)境下的M文件，編寫(xiě)程序來(lái)實(shí)現(xiàn)PSK的調(diào)制解調(diào)，并繪制出解調(diào)前后的時(shí)域和頻域波形及疊加噪聲時(shí)解調(diào)前后的時(shí)頻波形，根據(jù)運(yùn)行結(jié)果和波形來(lái)分析該解調(diào)過(guò)程的正確性及信道對(duì)信號(hào)傳輸?shù)挠绊?。要求?）熟悉MATLAB中M文件的使用方法，并在掌握PSK調(diào)制解調(diào)原理的基礎(chǔ)上，編寫(xiě)出PSK調(diào)制解調(diào)程序。2）繪制出PSK信號(hào)解調(diào)前后在時(shí)域和頻域中的波形，并觀(guān)察解調(diào)前后頻譜有何變

3、化以加深對(duì)PSK信號(hào)解調(diào)原理的理解。3）對(duì)信號(hào)疊加噪聲，并進(jìn)行解調(diào)，繪制出解調(diào)前后信號(hào)的時(shí)頻波形，改變?cè)肼暪β蔬M(jìn)行解調(diào)，分析噪聲對(duì)信號(hào)傳輸造成的影響。 4）在老師的指導(dǎo)下，要求獨(dú)立完成課程設(shè)計(jì)的全部?jī)?nèi)容，并按要求編寫(xiě)課程設(shè)計(jì)學(xué)年論文，能正確闡述和分析設(shè)計(jì)和實(shí)驗(yàn)結(jié)果。應(yīng)當(dāng)提交的文件：（1）課程設(shè)計(jì)學(xué)年論文。（2）課程設(shè)計(jì)附件?；贛ATLAB的PSK調(diào)制解調(diào)實(shí)現(xiàn) 學(xué)生姓名：李柳 指導(dǎo)老師：吳志敏摘 要 利用MATLAB集成環(huán)境下的M文件，編寫(xiě)程序來(lái)實(shí)現(xiàn)PSK的調(diào)制解調(diào)，并繪制出解調(diào)前后的時(shí)域和頻域波形及疊加噪聲時(shí)解調(diào)前后的時(shí)頻波形，根據(jù)運(yùn)行結(jié)果和波形來(lái)分析該解調(diào)過(guò)程的正確性及信道對(duì)信號(hào)傳輸?shù)挠绊?/p>

4、。 關(guān)鍵詞 數(shù)字調(diào)制解調(diào) 分析與仿真 Matlab 1 引 言現(xiàn)代通信的發(fā)展趨勢(shì)為數(shù)字化，隨著現(xiàn)代通信技術(shù)的不斷開(kāi)發(fā)，數(shù)字調(diào)制技術(shù)已日趨成熟，在各個(gè)領(lǐng)域都得到了廣泛的應(yīng)用和認(rèn)同。數(shù)字通信系統(tǒng)較模擬通信系統(tǒng)而言，具有抗干擾能力強(qiáng)、便于加密、易于實(shí)現(xiàn)集成化、便于與計(jì)算機(jī)連接等優(yōu)點(diǎn)。因而，數(shù)字通信更能適應(yīng)對(duì)通信技術(shù)的高要求。MATLAB是一種使用簡(jiǎn)便的、特別適用于科學(xué)研究和工程計(jì)算的高級(jí)語(yǔ)言，與其他計(jì)算機(jī)語(yǔ)言相比，它的特點(diǎn)是簡(jiǎn)潔和智能化，具有極高的編程和調(diào)試效率。通過(guò)使用MATLAB工具箱函數(shù)對(duì)數(shù)字調(diào)制進(jìn)行仿真，更能直觀(guān)徹底的掌握數(shù)字通信，數(shù)字調(diào)制的原理。有助于我們的學(xué)習(xí)和研究，加深對(duì)知識(shí)的理解和運(yùn)

5、用。1.1課程設(shè)計(jì)目的要求熟練使用MATLAB集成環(huán)境下的M文件，編寫(xiě)程序來(lái)實(shí)現(xiàn)PSK的調(diào)制解調(diào)。通過(guò)實(shí)現(xiàn)PSK調(diào)制解調(diào)的課程設(shè)計(jì)，更加熟練的掌握MATLAB的編程語(yǔ)言，熟練使用條件控制，循環(huán)等語(yǔ)句以及函數(shù)的調(diào)用等關(guān)于MATLAB的知識(shí)。同時(shí)加深對(duì)PSK調(diào)制與解調(diào)的理論知識(shí)的理解，增強(qiáng)matlab的軟件的使用方面和各種編碼解碼的能力，對(duì)以后其他的軟件的使用和其他方式的調(diào)制解調(diào)能更好的理解等。1.2課程設(shè)計(jì)要求能夠理解并掌握PSK調(diào)制解調(diào)的原理，能夠熟練的利用MATLAB中的M文件編寫(xiě)PSK的調(diào)制解調(diào)程序并繪制出PSK信號(hào)解調(diào)前后在時(shí)域和頻域的波形，觀(guān)察其頻譜的變化；在此基礎(chǔ)上對(duì)信號(hào)進(jìn)行疊加噪聲

6、，繪制出解調(diào)前后的時(shí)域波形，改變?cè)肼暪β蔬M(jìn)行解調(diào)，分析噪聲對(duì)信號(hào)傳輸?shù)挠绊憽?1.3課程設(shè)計(jì)步驟 在M文件中編寫(xiě)程序產(chǎn)生數(shù)字基帶信號(hào)，對(duì)信號(hào)進(jìn)行數(shù)字調(diào)制，再已調(diào)信號(hào)進(jìn)行解調(diào)，繪制調(diào)制解調(diào)前后信號(hào)在時(shí)域的波形和頻譜 ，然后在已調(diào)信號(hào)上疊加噪聲并解調(diào)，改變?cè)肼暪β剩^(guān)察解調(diào)前后時(shí)域波形，分析噪聲對(duì)信號(hào)傳輸?shù)挠绊憽＃?） 利用正弦信號(hào)產(chǎn)生方波，作為基帶信號(hào)。（2） 基帶信號(hào)與載波相乘得到PSK信號(hào)。（3） 通過(guò)帶通濾波器，相乘器，低通濾波器的到解調(diào)信號(hào)。（4） 疊加噪聲信號(hào)再次進(jìn)行解調(diào)，改變?cè)肼暪β试龠M(jìn)行解調(diào)。（5） 分析對(duì)比解調(diào)前后的時(shí)域波形，解釋噪聲對(duì)信號(hào)傳輸?shù)挠绊憽? PSK調(diào)制解調(diào)原理調(diào)制和

7、解調(diào)的基本原理是利用信號(hào)與系統(tǒng)的頻域分析和傅里葉變換的基本性質(zhì)，將信號(hào)的頻譜進(jìn)行搬移，使之滿(mǎn)足一定需要，從而完成信號(hào)的傳輸或處理。調(diào)制與解調(diào)又分模擬和數(shù)字兩種，在現(xiàn)代通信中，調(diào)制器的載波信號(hào)幾乎都是正弦信號(hào)，數(shù)字基帶信號(hào)通過(guò)調(diào)制器改變正弦載波信號(hào)的幅度、頻率或相位，產(chǎn)生幅度鍵控（ASK）、相位鍵控（PSK）、頻率鍵控（FSK）信號(hào)，或同時(shí)改變正弦載波信號(hào)的幾個(gè)參數(shù)，產(chǎn)生復(fù)合調(diào)制信號(hào)。相移鍵控(PSK)：一種用載波相位表示輸入信號(hào)信息的調(diào)制技術(shù)。移相鍵控分為絕對(duì)移相和相對(duì)移相兩種。以未調(diào)載波的相位作為基準(zhǔn)的相位調(diào)制叫作絕對(duì)移相。以二進(jìn)制調(diào)相為例，取碼元為“1”時(shí)，調(diào)制后載波與未調(diào)載波同相；取碼元

8、為“0”時(shí)，調(diào)制后載波與未調(diào)載波反相；“1”和“0”時(shí)調(diào)制后載波相位差180°。2.1 PSK調(diào)制原理PSK信號(hào)用載波相位的變化來(lái)表征被傳輸信息的狀態(tài)，通常規(guī)定0相位載波和相位載波分別表示傳“1”和傳“0”。設(shè)二進(jìn)制單極性碼為an,其對(duì)應(yīng)的雙極性二進(jìn)制碼為bn,則2PSK信號(hào)的一般時(shí)域信號(hào)可以表示為： S2psk（t）= bn g(t-nTs)cosct式中 bn=-1（當(dāng)an=0時(shí)，概率為P） bn=1（當(dāng)an=1時(shí)，概率為1-P）則時(shí)域信號(hào)可以變?yōu)? S2psk（t）= g(t-nTs)cos（ct+）， 當(dāng)an=0時(shí) S2psk（t）= g(t-nTs)cos（ct+0） 當(dāng)a

9、n=1時(shí)由此可知2PSK信號(hào)是一種雙邊帶信號(hào)，功率譜為：P2PSK（）= s=P（1-P）|G（+ s）|2+|G（- s）|2 + s2（1-P）2|G（0）|2(+ s)+ (- s)2PSK信號(hào)的帶寬為B2PSK=（c+Rs）-（c-Rs）= 2Rs 式中Rs為碼元速率。碼變換相乘S（t）載波eo（t）雙極性不歸零2PSK的調(diào)制方框圖 圖2.12.2 PSK解調(diào)原理 2PSK信號(hào)的解調(diào)通常采用相干解調(diào)法，其原理框圖如下：定時(shí)脈沖帶通濾波器低通濾波器抽樣判決器輸出 圖 2.2相干解調(diào)也叫同步檢波，它適用于所有線(xiàn)性調(diào)制信號(hào)的解調(diào)。實(shí)現(xiàn)相干解調(diào)的關(guān)鍵是接收端要恢復(fù)出一個(gè)與調(diào)制載波嚴(yán)格同步的相干

10、載波。相干解調(diào)是指利用乘法器，輸入一路與載頻相干（同頻同相）的參考信號(hào)與載頻相乘。 3 仿真實(shí)現(xiàn)過(guò)程3.1 PSK信號(hào)的產(chǎn)生利用正弦波產(chǎn)生基帶數(shù)字信號(hào)（雙極性），再與載波信號(hào)直接相乘的方法。x = cos(wfc_t); x(x>0)=1; x(x<0)=-1; 利用正弦波產(chǎn)生雙極性基帶數(shù)字信號(hào)；car=sin(2*pi*fc*t); 載波 ； psk=x.*car; 載波與基帶信號(hào)相乘進(jìn)行載波調(diào)制；如圖3.1為基帶信號(hào)與已調(diào)制信號(hào)。 圖3.13.2 PSK解調(diào)實(shí)現(xiàn) （1）構(gòu)建帶通濾波器，將已調(diào)信號(hào)通過(guò)帶通濾波器得到濾波信號(hào)。Hz=0.135*(z+1).*(z-1)./(z.2-

11、(p); 利用此系統(tǒng)響應(yīng)函數(shù)構(gòu)建帶通濾波器。faskn=filter(b,a,ask); 利用filter函數(shù)得到濾波信號(hào)；（2）將通過(guò)帶通濾波器的信號(hào)與載波信號(hào)相乘?？梢杂靡幌鲁绦蛘Z(yǔ)句進(jìn)行操作，從而實(shí)現(xiàn)乘法器的功能。cm=faskn.*car; 解調(diào) （3）將經(jīng)過(guò)乘法器的信號(hào)通過(guò)低通濾波器進(jìn)行濾波。p=0.72; gain1=0.14;%gain=(1-p)/2 Hz1=gain1*(z+1)./(z-(p); 利用此系統(tǒng)響應(yīng)函數(shù)構(gòu)建低通濾波器。so=filter(b1,a1,cm); 利用filter函數(shù)得到濾波信號(hào)； （4）設(shè)置判決電平，設(shè)置為最佳判決電平0，將信號(hào)通過(guò)判決器得到解調(diào)輸出

12、信號(hào)。下圖是未疊加噪聲的調(diào)制解調(diào)圖形： 圖3.2 對(duì)比基帶信號(hào)，已調(diào)信號(hào)，解調(diào)輸出信號(hào)，參考通信原理教程，已調(diào)信號(hào)的帶寬等于2倍的基帶信號(hào)帶寬（即fbw=2fm），由得出的各個(gè)信號(hào)的時(shí)域波形可知滿(mǎn)足要求，對(duì)比基帶信號(hào)和解調(diào)輸出信號(hào)的頻譜圖可以看出，兩者都出現(xiàn)兩個(gè)邊帶的頻譜。 圖3.3分析：由上圖可知基帶信號(hào)與解調(diào)信號(hào)基本相同，并未出現(xiàn)失真，只是存在一定延時(shí)。說(shuō)明相干解調(diào)的方法能實(shí)現(xiàn)PSK信號(hào)的解調(diào)無(wú)失真輸出。3.3疊加噪聲的PSK解調(diào)在實(shí)際的信號(hào)傳輸過(guò)程中，系統(tǒng)外部干擾以電磁波或經(jīng)電源串進(jìn)系統(tǒng)內(nèi)部而引起噪聲。如電氣設(shè)備，天體放電現(xiàn)象等引起的噪聲。同時(shí)由光和電的基本性質(zhì)所引起的噪聲，電器的機(jī)械運(yùn)

13、動(dòng)產(chǎn)生的噪聲，器材材料本身引起的噪聲，系統(tǒng)內(nèi)部設(shè)備電路也能引起噪聲。這些噪聲有些是不可避免的，所以我們?cè)谘芯啃盘?hào)的傳輸是一定要疊加上噪聲，這樣才有實(shí)際意義。下面是對(duì)疊加上大信噪比和小信噪比之后的信號(hào)的調(diào)制與解調(diào)做分析。 信號(hào)疊加上大信噪比的噪聲之后，信號(hào)在幅度上有微弱變化，在經(jīng)過(guò)帶通濾波器之后，與上圖為疊加噪聲的對(duì)比可以發(fā)現(xiàn)，幾乎已經(jīng)一樣，在通過(guò)乘法器和低通濾波器之后仍能恢復(fù)出原基帶信號(hào)。疊加小噪聲之后的圖形： 圖3.4疊加上大信噪比的噪聲之后，由于疊加的上的噪聲信號(hào)由randn(size(t)函數(shù)產(chǎn)生，是完全隨機(jī)的信號(hào)，疊到已調(diào)信號(hào)上會(huì)把信號(hào)的波形完全打亂。如圖3.5噪聲信號(hào)和加噪聲之后的已

14、調(diào)信號(hào)： 圖3.5 下圖3.6是疊加大噪聲情況下的基帶信號(hào)經(jīng)調(diào)制解調(diào)之后的信號(hào)的時(shí)域頻域波形： 圖 3.64 遇到的問(wèn)題及解決辦法 在此次課程設(shè)計(jì)中，我認(rèn)為我的難點(diǎn)在于帶通與低通濾波器的構(gòu)建上，在構(gòu)建帶通濾波器的過(guò)程中因?yàn)閰?shù)的設(shè)置不正確算出的系統(tǒng)函數(shù)不正確，導(dǎo)致調(diào)制信號(hào)通過(guò)帶通濾波器，乘法器和低通濾波器之后無(wú)法得到無(wú)失真的解調(diào)信號(hào)。解決方案：調(diào)整濾波器的參數(shù)，帶通濾波器的中心頻率：fc=2×10e-5;所以可以求的通帶截至頻率為：f1=fc-fm=1.8×10e-5和f2=fc+fm=2.2×10e5；由這些參數(shù)可以求的的低通濾波器的系統(tǒng)函數(shù)為：Hz=0.135

15、×(z+1)(z-1)./(z.2+0.8547) ;調(diào)整之后的解調(diào)系統(tǒng)能夠符合要求。另外，就是各種函數(shù)的調(diào)用的問(wèn)題，因?yàn)橐婚_(kāi)始沒(méi)有仔細(xì)的查閱資料，對(duì)MATLAB中的已有可直接調(diào)用的函數(shù)了解不多，致使做了許多繁重的工作，使整個(gè)設(shè)計(jì)顯得復(fù)雜。5 結(jié)束語(yǔ)通過(guò)這段時(shí)間的親身經(jīng)歷，我感覺(jué)自己學(xué)到了：收集、整理資料、共同協(xié)作、分析及處理問(wèn)題等許多方面的知識(shí)。在作設(shè)計(jì)的過(guò)程中，我有許多不懂得地方，在老師的指導(dǎo)下我一步步的解決問(wèn)題完成論文，在完成過(guò)程中老師指導(dǎo)我去怎么選擇資料，如何去利用網(wǎng)絡(luò)資源，在這個(gè)學(xué)習(xí)的過(guò)程中，我了解到MATLAB的實(shí)用價(jià)值，更深的理解數(shù)字調(diào)制技術(shù)的調(diào)制原理。對(duì)三種鍵控方式的

16、調(diào)制原理也理解的更加透徹。同時(shí)也深感用計(jì)算機(jī)仿真電子通信系統(tǒng)，具有廣泛的適應(yīng)性和極高的靈活性。在硬件試驗(yàn)中改變硬件的參數(shù)設(shè)置就意味著重做硬件，而在軟件中只需對(duì)相應(yīng)的參數(shù)進(jìn)行重新設(shè)置，同時(shí)利用MATLAB簡(jiǎn)單編程的特點(diǎn)，可以實(shí)現(xiàn)較為復(fù)雜的系統(tǒng)，因此MATLAB在通信系統(tǒng)仿真方面具有強(qiáng)大的功能和優(yōu)越性！在課程設(shè)計(jì)的過(guò)程中明顯覺(jué)得自己對(duì)MATLAB軟件的不熟練，很多函數(shù)不會(huì)使用，還需要查閱參考書(shū)才能正確使用。而且再做到最后的時(shí)候才發(fā)現(xiàn)matlab中的y=dmod(x,Fc,Fd,Fs,method)函數(shù)可以直接進(jìn)行PSK的調(diào)制，而z=ddemod(y,Fc,Fd,Fs,method)函數(shù)也可直接對(duì)已

17、調(diào)信號(hào)進(jìn)行解調(diào)，不再需要設(shè)置低通，帶通濾波器，也不需要設(shè)置判決器就可直接調(diào)制解調(diào)，而自己的設(shè)計(jì)的解調(diào)做了很多繁重的工作，由于時(shí)間緊迫，所以沒(méi)有進(jìn)行優(yōu)化。在以后的學(xué)習(xí)或者其他設(shè)計(jì)中，應(yīng)該提前做好充足的準(zhǔn)備，對(duì)要做的東西有足夠的了解之后再開(kāi)始進(jìn)行，其次以后在工作中要腳踏實(shí)地，在學(xué)術(shù)上要嚴(yán)謹(jǐn)，在思維上要活躍，在學(xué)業(yè)上要勤奮刻苦，不能懵懂的就以為自己已經(jīng)會(huì)了。最后我真誠(chéng)感謝這期間老師給予我的全力幫助，細(xì)心指導(dǎo)以及對(duì)我的嚴(yán)格要求，是她在我遇到問(wèn)題時(shí)，不辭辛苦幫我解決，感謝她在設(shè)計(jì)和任務(wù)安排上長(zhǎng)時(shí)間的指導(dǎo)。參考文獻(xiàn)1 黃文梅, 熊桂林, 楊勇.信號(hào)分析與處理MATLAB語(yǔ)言及應(yīng)用. 長(zhǎng)沙: 國(guó)防科技大學(xué)出

18、版社, 2000.2唐向宏, 岳恒立, 鄭雪峰. MATLAB及在電子信息類(lèi)課程中的應(yīng)用.北京: 電子工業(yè)出版社, 2006, 8.3 鄧華. MATLAB通信仿真及應(yīng)用實(shí)例詳解. 人民郵電出版社, 2003.4 樊昌信 徐炳祥等.通信原理(第5版).北京.國(guó)防工業(yè)出版社.20055 程佩青 數(shù)字信號(hào)處理教程（第四版）清華大學(xué)出版社6王興亮.數(shù)字通信原理與技術(shù)（第二版）.西安電子科技大學(xué)出版社.2000附錄一：未疊加噪聲PSK調(diào)制解調(diào)程序清單%程序名稱(chēng)：PSK.m%程序功能：實(shí)現(xiàn)PSK調(diào)制解調(diào)%程序作者：李柳%=clear all; close all; fs=8e5;%抽樣頻率 fm=20e

19、3;%基帶頻率 n=2*(6*fs/fm); final=(1/fs)*(n-1); fc=2e5; % 載波頻率 t=0:1/fs:(final); Fn=fs/2;%耐奎斯特頻率 %用正弦波產(chǎn)生方波 wfc_t=2*pi*fm*t; x = cos(wfc_t ); % 方波 x(x>0)=1 x(x<0)=-1; figure(1) subplot(221); plot(t,x); axis(0 2e-4 -2 2); title('基帶信號(hào)'); grid on car=sin(2*pi*fc*t);%載波 ask=x.*car;%載波調(diào)制subplot(2

20、22); plot(t,ask); axis(0 200e-6 -2 2); title('PSK信號(hào)'); grid onfBW=40e3; %編碼后的PSK信號(hào)頻率f=0:3e3:4e5; w=2*pi*f/fs; z=exp(w*j); BW=2*pi*fBW/fs; %通帶帶寬 a=.8547;%BW=2(1-a)/sqrt(a) p=(j2*a2); Hz=0.135*(z+1).*(z-1)./(z.2-(p); %帶通濾波器的系統(tǒng)函數(shù)subplot(223); plot(f,abs(Hz); title('帶通濾波器');grid on;Hz(Hz

21、=0)=10(8);%避免 log(0) subplot(224); plot(f,20*log10(abs(Hz); grid on;title('Receiver -3dB Filter Response'); axis(1e5 3e5 -3 1); %濾波器系數(shù) a=1 0 0.7305;%1 0 p b=0.135 0 -0.135;%0.135*1 0 -1 faskn=filter(b,a,ask); figure(2) subplot(321); plot(t,faskn); axis(0 100e-6 -2 2); title('通過(guò)帶通濾波后輸出

22、9;); grid on;cm=faskn.*car;%解調(diào) subplot(322); plot(t,cm); axis(0 100e-6 -2 2); grid on; title('通過(guò)相乘器后輸出'); %低通濾波器 %= p=0.72; gain1=0.14;%gain1=(1-p)/2 Hz1=gain1*(z+1)./(z-(p); subplot(323); Hz1(Hz1=0)=10(-8);%避免 log(0) plot(f,20*log10(abs(Hz1); grid on;title('LPF -3dB response'); axis

23、(0 5e4 -3 1);%濾波器系數(shù) a1=1 -0.72;%(z-(p) b1=0.14 0.14;%gain*1 1 so=filter(b1,a1,cm); so=so*10;%add gain so=so-mean(so);%removes DC component subplot(324); plot(t,so); axis(0 8e-4 -3.5 3.5); title('通過(guò)低通濾波器后輸出'); grid on; %Comparator %= High=1; Low=-1; vt=0;%判決電平 len1=length(so); for ii=1:len1 i

24、f so(ii) >= vt Vs(ii)=High; else Vs(ii)=Low; end end Vo=Vs; subplot(325); plot (t,Vo); title('解調(diào)后輸出信號(hào)'); axis(0 2e-4 -2 2) ; grid on; xlabel('時(shí)間 (s)'), ylabel('幅度(V)');figure(3)F1=fft(x);subplot(211); plot (t,F1); title('基帶信號(hào)頻譜');F2=fft(so);subplot(212); plot (t,F2

25、); title('解調(diào)后輸出信號(hào)頻譜');附錄二：疊加大信噪比噪聲PSK調(diào)制解調(diào)程序清單%程序名稱(chēng)：PSK.m%程序功能：實(shí)現(xiàn)PSK調(diào)制解調(diào)%程序作者：李柳%=clear all; close all; fs=8e5;%抽樣頻率 fm=2e4;%基帶頻率 n=2*(6*fs/fm); final=(1/fs)*(n-1); fc=2e5; % 載波頻率 t=0:1/fs:(final); %用正弦波產(chǎn)生方波 wfc_t=2*pi*fm*t; x = cos(wfc_t ); % 方波 x(x>0)=1; x(x<0)=-1; figure(1) subplot(3

26、21); plot(t,x); axis(0 2e-4 -2 2); title('基帶信號(hào)'); grid on car=sin(2*pi*fc*t);%載波 ask=x.*car;%載波調(diào)制subplot(323); plot(t,ask); axis(0 200e-6 -2 2); title('PSK信號(hào)'); grid onvn=0.1; noise=vn*(randn(size(t);%產(chǎn)生噪音 subplot(322); plot(t,noise); grid on; title('噪音信號(hào)'); axis(0 .2e-3 -1 1

27、); askn=(ask+noise);%調(diào)制后加噪 subplot(324); plot(t,askn); axis(0 200e-6 -2 2); title('加噪后信號(hào)'); grid on;fBW=4e4; %已調(diào)信號(hào)帶寬f=0:3e3:4e5; w=2*pi*f/fs; z=exp(w*j); BW=2*pi*fBW/fs; a=.8547;%BW=2(1-a)/sqrt(a) p=(j2*a2); Hz=0.135*(z+1).*(z-1)./(z.2-(p); subplot(325); plot(f,abs(Hz); title('帶通濾波器'

28、);grid on;Hz(Hz=0)=10(8);%avoid log(0) subplot(326); plot(f,20*log10(abs(Hz); grid on;title('Receiver -3dB Filter Response'); axis(1e5 3e5 -3 1); %濾波器系數(shù) a=1 0 0.7305;%1 0 p b=0.135 0 -0.135;%gain*1 0 -1 faskn=filter(b,a,askn); figure(2) subplot(322); plot(t,faskn); axis(0 100e-6 -2 2); title('通過(guò)帶通濾波后輸出'); grid on;cm=faskn.*car;%解調(diào) subplot(323); plot(t,cm); axis(0 100e-6 -2 2); grid on; title('通過(guò)相乘器后輸