2PSK仿真課程設(shè)計-2PSK調(diào)制解調(diào)器的建模與仿真

1題目2PSK調(diào)制解調(diào)器的建模與仿真專業(yè)班級指導(dǎo)教師I課程設(shè)計任務(wù)書●2PSK調(diào)制解調(diào)器的建模與仿真設(shè)計目的設(shè)計要求課程設(shè)計的目的1.通過利用matlab對數(shù)字通信相關(guān)內(nèi)容進行仿真設(shè)計，熟悉數(shù)字通信的相關(guān)理論；2.更好的了解數(shù)字通信的相關(guān)知識，提高分析問題、查閱資料、自主設(shè)計等各方面能力。課程設(shè)計的要求1.產(chǎn)生一組20位以上隨機二進制數(shù)字基帶信號(雙極性);2.產(chǎn)生載波信號，載波周期與碼元寬度相同且均為1;3.實現(xiàn)基于MATLAB仿真2PSK調(diào)制與解調(diào)；4.通過巴特沃斯低通濾波器，根據(jù)位同步信號進行抽樣判決；5.各階段信號頻譜分析、誤碼率并給出仿真結(jié)果；設(shè)計進度安排或工作計劃2014.12.30～2014.12.31:熟悉課題，查詢相關(guān)資料，完成方案選擇；2015.1.2~2015.1.6:設(shè)計模塊劃分、實現(xiàn)及各模塊調(diào)試、驗證；2015.1.7~2015.1.8:設(shè)計整體實現(xiàn)、調(diào)試及驗證，并開始撰寫報告：2015.1.9~2015.1.10:設(shè)計完成，課程設(shè)計報告撰寫并定稿，上交。其它認(rèn)真閱讀數(shù)字通信課程設(shè)計報告撰寫規(guī)范；課題小組經(jīng)協(xié)商好要指定組長并明確分工，形成良好團隊工作氛圍；基于課題基本要求，各小組與指導(dǎo)老師討論，再將課題細(xì)化、增加要求；課題小組每個成員均需各自撰寫一份課程設(shè)計報告。2PSK調(diào)制解調(diào)器的建模與仿真數(shù)字傳輸系統(tǒng)中，數(shù)字信號對高頻載波進行調(diào)制，變?yōu)轭l帶信號，通過信道傳輸，在接收端解調(diào)后恢復(fù)成數(shù)字信號，由于大多數(shù)實際信號都是帶通型的，所以必須先用數(shù)字頻帶信號對載波進行調(diào)制，形成數(shù)字調(diào)制信號再進行傳輸，因此，調(diào)制解調(diào)技術(shù)是實現(xiàn)現(xiàn)代通信的重要手段。數(shù)字調(diào)制地實現(xiàn)，促進了通信的飛速發(fā)展。研究數(shù)字通信調(diào)制理論，提供有效調(diào)制方法有著重要意義。實現(xiàn)調(diào)試解調(diào)的方法有很多種，本文應(yīng)用了鍵控法產(chǎn)生調(diào)制與解調(diào)信號。數(shù)字相位調(diào)制又稱相移鍵控記作PSK(PhaseShiftKeying),二進制相移鍵控記作2PSK,它們是利用載波振蕩相位的變化來傳送數(shù)字信號的，在二進制數(shù)字解調(diào)中，當(dāng)正弦載波的相位隨二進制數(shù)字基帶信號離散變化，則就產(chǎn)生二進制移相鍵控2PSK信號。重點介紹了2PSK的調(diào)制與解調(diào)的工作原理，以及通過MATLAB進行設(shè)計和仿真，并且分析了各階段信號的頻譜及誤碼率。關(guān)鍵詞數(shù)字調(diào)制與解調(diào)；2PSK;MATLAB仿真課程設(shè)計任務(wù)書 I Ⅱ1設(shè)計概述 11.1設(shè)計背景 11.2設(shè)計目的 11.3設(shè)計內(nèi)容及要求 12設(shè)計原理 22.12PSK數(shù)字調(diào)制過程分析 22.22PSK數(shù)字解調(diào)過程分析 32.32PSK相干解調(diào)系統(tǒng)性能分析 43設(shè)計實現(xiàn) 53.1系統(tǒng)設(shè)計框圖 53.2M文件下仿真設(shè)計思路及過程 53.2.12PSK調(diào)制部分 53.2.22PSK解調(diào)部分 6 63.3.1仿真模型建立 63.3.2參數(shù)的設(shè)置 73.3.3仿真結(jié)果 73.3.4仿真結(jié)果分析 4M文件仿真結(jié)果 94.12PSK調(diào)制、解調(diào)系統(tǒng)信號波形圖 94.22PSK調(diào)制解調(diào)系統(tǒng)頻譜分析 4.3各階段功率譜密度分析及誤碼率 參考文獻 附錄 11.1設(shè)計背景當(dāng)今社會已經(jīng)步入信息時代，在各種信息技術(shù)中，信息的傳輸及通信起著支撐作用。而對于信息的傳輸，數(shù)字通信已經(jīng)成為重要的手段。因此，數(shù)字信號的調(diào)制就顯得非常重要。數(shù)字信號的傳輸方式分為基帶傳輸和帶通傳輸，在實際應(yīng)用中，大多數(shù)信道具有帶通特性而不能直接傳輸基帶信號。為了使數(shù)字信號在帶通信道中傳輸，必須使用數(shù)字基帶信號對載波進行調(diào)制，以使信號與信道的特性相匹配。這種用數(shù)字基帶信號控制載波，把數(shù)字基帶信號變換為數(shù)字帶通信號的過程稱為數(shù)字調(diào)制。通常使用鍵控法來實現(xiàn)數(shù)字調(diào)制，比如對載波的振幅、頻率和相位進行鍵控。在這三種調(diào)制的基礎(chǔ)上為了得到更高的效果也出現(xiàn)了很多其他的調(diào)制方式，如：DPSK,MASK,MFSK,MPSK,APK。它們其中有的一些是將基本的調(diào)制方式用在多進制上或者引入了一些新的方式來解決基本調(diào)制的一些問題如相位模糊和無法提取位定時信號，另外一些是由多種基本的調(diào)制方式組合來達(dá)到更數(shù)字調(diào)制與解調(diào)系統(tǒng)性能的好壞取決于傳輸信號的誤碼率，而誤碼率不僅僅與信道、接收方法有關(guān)還和發(fā)送端所采用的調(diào)制方式有很大的關(guān)系。我們研究的ASK,FSK,PSK等就主要是發(fā)送方的基本調(diào)制方式。本設(shè)計主要對2PSK信號的原理及其相干解調(diào)系統(tǒng)性能進行了分析和仿真，這樣能讓我們對數(shù)字調(diào)制方式有一個更清楚的認(rèn)識。1.2設(shè)計目的2PSK數(shù)字調(diào)相技術(shù)由于其抗干擾能力強，實現(xiàn)簡單，而被廣泛應(yīng)用于各種通信中。本題目要求學(xué)生進行2PSK信號的產(chǎn)生及譜分析，達(dá)到以下目的：1.通過利用MATLAB對數(shù)字通信相關(guān)內(nèi)容進行仿真設(shè)計，熟悉數(shù)字通信的相關(guān)理論；2.更好的了解數(shù)字通信的相關(guān)知識，提高分析問題、查閱資料、自主設(shè)計等各方面能力。3.了解并掌握數(shù)字頻帶傳輸系統(tǒng)—2PSK的調(diào)制與解調(diào)系統(tǒng)，通過MATLAB來實現(xiàn)整個通信系統(tǒng)的設(shè)計，進一步加深對2PSK系統(tǒng)的理解，并熟練掌握MATLAB語言。1.3設(shè)計內(nèi)容及要求1.產(chǎn)生一組20位以上隨機二進制數(shù)字基帶信號(雙極性);2.產(chǎn)生載波信號，載波周期與碼元寬度相同且均為1;3.實現(xiàn)基于MATLAB仿真2PSK調(diào)制與4.通過巴特沃斯低通濾波器，根據(jù)位同步信號進行抽樣判決；5.各階段信號頻譜分析、誤碼率并給出仿真結(jié)果。22設(shè)計原理2.12PSK數(shù)字調(diào)制過程分析數(shù)字信號的傳輸方式分為基帶傳輸和帶通傳輸，在實際應(yīng)用中，大多數(shù)信道具有帶通特性而不能直接傳輸基帶信號。為了使數(shù)字信號在帶通信道中傳輸，必須使用數(shù)字基帶信號對載波進行調(diào)制，以使信號與信道的特性相匹配。這種用數(shù)字基帶信號控制載波，把數(shù)字基帶信號變換為數(shù)字帶通信號的過程稱為數(shù)字調(diào)制。數(shù)字調(diào)制技術(shù)的兩種方法：①利用模擬調(diào)制的方法去實現(xiàn)數(shù)字式調(diào)制，即把數(shù)字調(diào)制看成是模擬調(diào)制的一個特例，把數(shù)字基帶信號當(dāng)作模擬信號的特殊情況處理；②利用數(shù)字信號的離散取值特點通過開關(guān)鍵控載波，從而實現(xiàn)數(shù)字調(diào)制。這種方法通常稱為鍵控法，比如對載波的相位進行鍵控，便可獲得相移鍵控(PSK)基本的調(diào)制方式。e?psk(t)開關(guān)電路開關(guān)電路0π圖2-1模擬調(diào)制的方法圖2-2相移方法數(shù)字調(diào)相：如果兩個頻率相同的載波同時開始振蕩，這兩個頻率同時達(dá)到正最大值，同時達(dá)到零值，同時達(dá)到負(fù)最大值，它們應(yīng)處于"同相"狀態(tài)；如果其中一個開始得遲了一點，就可能不相同了。如果一個達(dá)到正最大值時，另一個達(dá)到負(fù)最大值，則稱為"反相"。一般把信號振蕩一次(一周)作為360度。如果一個波比另一個波相差半個周期，我們說兩個波的相位差180度，也就是反相。當(dāng)傳輸數(shù)字信號時，"1"碼控制發(fā)0度相位，"0"碼控制發(fā)180度相位。載波的初始相位就有了移動，也就帶上了信息。相移鍵控：利用載波的相位變化來傳遞數(shù)字信息，而振幅和頻率保持不變。在2PSK因此，上式可以改寫為3由于表示信號的兩種碼元的波形相同極性相反，故2PSK信號一般可以表述為一個雙極性全占空比矩形脈沖序列與一個正弦載波的相乘，即信號在相乘器中與本地載波相乘，然后用低通濾波器濾除高頻分量，再進行抽樣判決。判決器是按極性來判決的。即正抽樣值判為1,負(fù)抽樣值判為0。e?psx(t)bdCa圖2-32PSK相干解調(diào)系統(tǒng)框圖圖2-42PSK信號相干解調(diào)各點時間波形帶通濾波器的意義是讓有用信號(已調(diào)信號)通過，濾除一部分噪聲，所以有用信號4在a處得到信號為假設(shè)相干載波的基準(zhǔn)相位與2PSK信號的調(diào)制載波的基準(zhǔn)相位一致(通常默認(rèn)為0相載波為參考的，因此，解調(diào)時必須有與此同頻同相的同步載波。如果同步載波的相位發(fā)生2.32PSK相干解調(diào)系統(tǒng)性能分析s?(t)y;(t)y(t)x(t)輸出P。在大信噪比(r→1)條件下，上式可近似為：5低低濾通波3.1系統(tǒng)設(shè)計框圖20位二進制單極性碼單極性碼轉(zhuǎn)換為雙極性碼調(diào)制2PSK誤碼率仿真加性圖3-1系統(tǒng)設(shè)計框圖3.2M文件下仿真設(shè)計思路及過程本設(shè)計主要是通過MATLAB來實現(xiàn)整個2PSK系統(tǒng)，在用MATLAB語言編程的時候主要分為以下幾個部分。(1)通過碼型變換產(chǎn)生雙極性數(shù)字基帶信號。首先調(diào)用MATLAB的rand函數(shù)產(chǎn)生一組隨機序列，然后將離散的隨機序列變換成單極性的二進制數(shù)字基帶信號，將單極性信號乘以2然后減去1即可得到我們所需要的二進制雙極性數(shù)字基帶信號。其中單(雙)極性碼元寬度為1,每個碼元采樣250點即采樣間隔為1/250。(2)載波信號的產(chǎn)生。采用載頻為f.=1Hz的余弦載波信號cos(2*π*f*t),初始相位設(shè)為0,使得載波周期與碼元寬度相同且均為1。(3)通過隨機二進制雙極性數(shù)字基帶信號與載波信號相乘得到2PSK信號(0變1不變):psk=Bipolar_NRZ.*carrier_signal;(4)已調(diào)信號通過信道傳輸時會受到噪聲影響，我們通過MATLAB的awgn函數(shù)來模擬6高斯白噪聲與已調(diào)信號相加的過程，在這里我們只考慮加性噪聲，不考慮乘型噪聲。其中添加的高斯白噪聲是由yl=awgn(psk,10)實現(xiàn)，其中awgn為加性高斯白噪聲信道模型，信道的信噪比為20。(1)相干解調(diào)的實現(xiàn)過程是將接收機收到信號通過乘法器與載波信號相乘即y2=y1.*carrier_signal(2)用MATLAB來設(shè)計理想的巴特沃斯模擬低通濾波器，由于載波信號為1HZ,而要濾去的高頻信號為2HZ,故設(shè)通帶截止頻率1HZ,阻帶截止頻率2HZ。(3)進行抽樣判決，對每一個碼元寬度中間值作判決抽樣，增加抽樣成功率。若判決值大于0,則輸出值置為1,否則置為0。(3)求解誤碼率，利用公式(2-6)可以計算誤碼率的理論值，信噪比snrdB是從-10到10變化，將輸入的隨機二進制信號與抽樣判決值進行比較，誤碼率即為誤碼數(shù)與碼元總3.3.1仿真模型建立2PSK調(diào)制與解調(diào)及誤碼分析的總體仿真模型方框圖如下所示，上半部分為調(diào)制部分，GeneratorConverterXAWGNAWGNConyerterConverter圖3-2Simulink仿真模型7(1)正弦載波(SineWave2)參數(shù)設(shè)置：設(shè)置的振幅是1,頻率為1HZ,相位為0。(2)與載波反相正弦波(SineWavel)參數(shù)設(shè)置：設(shè)置的振幅是-1,頻率也是1hz,相位為0。(3)伯努利二進制隨機序列產(chǎn)生器(BemoulllBinaryGenerator):0和1的出現(xiàn)概率為0.5,抽樣時間為1秒，采用是基于采樣的，其幅度設(shè)置為2,周期為3,占1比為2/3。(4)碼型變化器(UnipolartoBipolarConverter)參數(shù)設(shè)置：設(shè)置依據(jù)：采用1變0不變調(diào)制，故極性設(shè)為“Positive”。(5)多路選擇器(Switch)參數(shù)設(shè)置：當(dāng)二進制序列大于0時，輸出第一路信號；當(dāng)二進制序列小于0時，輸出第二路信號。(6)帶通濾波器(DigitalFilterDesign)參數(shù)：帶通范圍為1~2HZ,載波頻率為1HZ,而基帶號帶寬為1HZ,考濾到濾波器的邊沿緩降，故設(shè)置為1~2HZ。(7)低通濾波器(DigitalFilterDesign)參數(shù)設(shè)置：截止頻率為1HZ,二進制序列的帶寬為1HZ,故取1HZ。(8)取樣判決器(Sign)參數(shù)設(shè)置：門限值取為0.5,取樣時間為1,當(dāng)大于0.5時輸出0,當(dāng)小于0.5時輸出1,能達(dá)到在1變0不變的取樣規(guī)則下正確解碼的目的。圖3-32PSK調(diào)制的波形圖中第一個圖為正相載波的波形，第二個圖為隨機產(chǎn)生的二進制序列，第三個圖為通8過碼型變換器后的波形，最后一個圖為調(diào)制后的2PSK信號。解調(diào)波形：WWVWWWWwwwwMwWwMwwMwwwwwww圖3-4加入高斯白噪聲后的解調(diào)波形圖中第一個圖為收到的2PSK波形(加入了高斯噪聲),第二個圖為與同頻同向載波相乘后的波形，第三個圖為抽樣判決后的雙極性二進制波形，最后一個圖是經(jīng)過極性變換后的波形。3.3.4仿真結(jié)果分析在Simulink下能夠清晰的看到2PSK調(diào)制和解調(diào)的波形。2PSK數(shù)字調(diào)制系統(tǒng)采用的是雙極性矩形波，所以在基帶信號產(chǎn)生后要進行極性變換。2PSK的調(diào)制是通過鍵控法實現(xiàn)的，還可以采用模擬調(diào)制的方法。2PSK的解調(diào)只有相干解調(diào)一種辦法。在加入高斯白噪聲后，依然可以得到解調(diào)的波形，但是存在一定的誤碼率。且信噪比越大，誤碼率越低。94M文件仿真結(jié)果4.12PSK調(diào)制、解調(diào)系統(tǒng)信號波形圖圖4-1二進制數(shù)字基帶信號與載波信號圖4-22PSK信號與接收機接收到的信號波形由公式(2-4)可知，通過低通濾波器以后，輸出波形近似于極性基帶信號波形，且幅度減半。圖4-4判決輸出波形個2PSK的調(diào)制、解調(diào)系統(tǒng)是正確的。4.22PSK調(diào)制解調(diào)系統(tǒng)頻譜分析數(shù)字基帶信號為二進制非歸零碼，所以它的傅里葉變換應(yīng)該是門函數(shù)，信號帶寬為1/TS,其頻譜圖-10。載波信號為正弦信號，載頻為1Hz,其頻譜如下。圖4-5數(shù)字基帶信號與載波信號頻譜2PSK信號的實質(zhì)為基帶信號與載波相乘的結(jié)果，所以2PSK的信號頻譜是基帶信號頻譜的線性搬移，中心頻率為fc.,帶寬為基帶信號帶寬的兩倍，如圖-11,2PSK信號通過信道后，加上了高斯白噪聲，其頻譜如下。圖4-62PSK信號與接收機接收的信號頻譜振幅相干解調(diào)實質(zhì)是將已調(diào)信號與提取的載波信號相乘，所以解調(diào)輸出的信號的頻率成分包括基帶信號頻率和載頻信號和噪聲如圖-12,通過低通濾波器后，高頻成分被慮掉，剩振幅下基帶信號和部分噪聲。頻率/Hz圖4-7解調(diào)輸出與低通輸出的信號頻譜圖4-8判決輸出信號頻譜與基帶信號頻譜通過比較判決輸出信號的頻譜和基帶信號頻譜，得出結(jié)論，2PSK的調(diào)制解調(diào)系統(tǒng)設(shè)計是正確的。誤碼率解調(diào)輸出信號功率譜密度判決信號功率譜密度4.3各階段功率譜密度分析及誤碼率誤碼率解調(diào)輸出信號功率譜密度判決信號功率譜密度載波信號功率譜密度載波信號功率譜密度接收機接收信號功率譜密度低通信號功率譜密度圖4-92PSK功率譜密度圖4-102PSK系統(tǒng)誤碼率由圖可以看出系統(tǒng)的誤碼率為0。通過在本次設(shè)計中的實踐明白了自己知識上的誤區(qū)，例如，在低通濾波的過程中，主要目的是濾去高頻分量，濾去載波成分，所以對于低通濾波器的截止頻率的設(shè)置較為關(guān)鍵。而在2PSK的調(diào)制與解調(diào)中所用信號為雙極性的信號，因此要將在本次設(shè)計中產(chǎn)生的單極性信號經(jīng)過碼的變化形成雙極性碼來傳輸。本次設(shè)計只是按理論上的知識結(jié)構(gòu)進行簡單地系統(tǒng)構(gòu)建，目的是明確數(shù)字基帶傳輸?shù)脑砑斑^程，而對于具體問題，例如，實際中信道噪聲一般為高斯白噪聲，本次設(shè)計為簡便并未采用而是用rand函數(shù)產(chǎn)生了隨機噪聲信號。此外，還有在2PSK實際傳輸系統(tǒng)中，在恢復(fù)載波的過程中會出現(xiàn)“倒π現(xiàn)象”即相位模糊現(xiàn)象，但是在本次設(shè)計中直接在解調(diào)時給其同頻同相的載波，所以不會出現(xiàn)此種現(xiàn)象，這是自己第一次利用Matlab編程功能實現(xiàn)通信原理中基礎(chǔ)知識系統(tǒng)地構(gòu)建和利用Matlab中Simulink模塊搭建系統(tǒng)來實現(xiàn)，通過自己親自去動手和調(diào)試我明白了實踐的重要性，尤其是對程序的調(diào)試，更需要大量的時間反復(fù)上機運行，發(fā)現(xiàn)錯誤并改正，這樣也就加強了自己對程序分析的能力，更深刻地明白了通信原理中的知識內(nèi)容，更進一步懂得了理論聯(lián)系實際的含義，同時提高了自己的思考能力，使得自己對課本里的內(nèi)容理解、記憶地更加透徹，這無論是在我以后的工作中或是學(xué)習(xí)中都是非常有用的。參考文獻[1]郝建軍，桑林.數(shù)字通信[M].北京：人民郵電出版社，2002(2):206~207[1]樊昌信，曹麗娜.通信原理[M].北京：國防工業(yè)出版社，2013(6):188~190.[3]李明明，李白萍.電子信息類專業(yè)MATLAB實驗教程[M].北京：北京大學(xué)出版社，2011(1):102~208.[4]劉學(xué)勇.詳解MATLAB/Simulink通信系統(tǒng)建模與仿真[M].北京：電子工業(yè)出版社，2011(1):附錄%數(shù)字通信課程設(shè)計2PSK調(diào)制解調(diào)器的建模與仿真i=20;j=5000;%i為隨機數(shù)個數(shù)，j為采樣點數(shù)t=linspace(0,20,5000);%在0到20之間等分成5000個點，每個碼元250點%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%功能實現(xiàn)模塊一：產(chǎn)生隨機二進制數(shù)字基帶信號%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%rem=round(rand(1,i);%調(diào)用函數(shù)產(chǎn)生0-1隨機序列，round按四舍五入取整%將離散的的隨機序列轉(zhuǎn)變成連續(xù)的單極性二進制碼元ifrem(n)<1;%若隨機數(shù)值小于1則單極性信號賦值為0form=j/i*(n-1)+1j/i*n%m=250*(n-1):250*nelse%若隨機數(shù)值等于1則單極性信號賦值為1form=j/i*(n-1)+1j/i*nBipolar_NRZ(m)=1;plot(t,Bipolar_NRZ);title(‘雙極性二進制數(shù)字基帶信號);axis([0,20,-2,2]);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%功能實現(xiàn)模塊二：產(chǎn)生載波信號%%%%%%%%%%%%%載波周期與碼元寬度相同且均為1%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%fc=1;%載波頻率為1Hzcarrier_signal=cos(2*pi*fc*t);%載波周期與碼元寬度相同且均為1subplot(212);plot(t,carrier_signal);axis([020-22]);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%功能實現(xiàn)模塊三：2psk調(diào)制部分：基帶信號與載波相乘%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%psk=Bipolar_NRZ.*carrier_signal;%基帶信號與載波相乘figure(2);subplot(211);plot(t,psk);axis([0,20,-2,2]);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%功能實現(xiàn)模塊四：調(diào)制信號通過信道傳輸產(chǎn)生加性噪聲%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%yl=awgn(psk,20);%信號psk中加入白噪聲，信噪比為10subplot(212);plot(t,y1);axis([0,20,-1,2]);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%功能實現(xiàn)模塊六：相干解調(diào)實現(xiàn)過程%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%y2=y1.*carrier_signal;%接收機收到信號輸入乘法器與載波信號相乘figure(3)subplot(211);plot(t,y2);title(解調(diào)輸出波形);grid;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%功能實現(xiàn)模塊六：通過巴特沃斯低通濾波器%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%fp=1;f=2;rp=1;rs=30;fn=300;%設(shè)計巴特沃斯低通濾波器的參數(shù)，通帶截止頻率1HZ,阻帶截止頻率2HZ,采樣率400ws=fs/(fn/2);wp=fp/(fn/2);%計算歸一化角頻率[n,wn]=buttord(wp,ws,tp,rs);%計算階數(shù)和截止頻率[b,a]=butter(n,wn);y3=filtfilt(b,a,y2);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%功能實現(xiàn)模塊七：根據(jù)位同步信號進行抽樣判決%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%y4=linspace(0,19);forn=1:20;ify3((2*n-1)*125)>0;%對每一個碼元寬度中間值作判決抽樣，增加抽樣成功率y4(n)=1;y4(n)=0;figure(4)subplot(211);stem(y4);axis([020-22]);title('抽樣后得到的離散采樣值);gridon;xlabel(時間y5=t;forn=1:20ify4(n)==1;%若當(dāng)前碼元的抽樣值為1,則該碼元的所有值(250點)為1form=j/i*(n-1)+1j/i*ny5(m)=1;form=j/i*(n-1)+1j/i*ny5(m)=-1;subplot(212)plot(t,y5);title('抽樣判決輸出基帶信號);xlabel(采樣點數(shù));%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%功能實現(xiàn)模塊八：各階段信號頻譜分析%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%fs=250;N=1024;%采樣頻率和數(shù)據(jù)點數(shù)n=0:N-1;t=n/fS;%時間序列y=ft(Bipolar_NRZ,N);%對信號進行FFT變換mag=abs(y);%求取Fourier變換的振幅plot(f(1:N/2),mag(1:N/2));xlabel(頻率/Hz);ylabel(振幅);title(基帶信號頻譜(N=1024));gridon;fs=250;N=1024;n=0:N-1;t=n/fs;y=fm(carrier_signal,N);%對信號進行FFT變換mag=abs(y);%求取Fourier變換的振幅plot(f(1:N/2),mag(1:N/2));xlabel(頻率/Hz);ylabel(振幅);title(載波信號頻譜(N=1024));gridon;fs=256;N=1024;n=0:N-1;t=n/fs;y=fh(psk,N);%對信號進行FFT變換mag=abs(y);%求取Fourier變換的振幅plot(f(1:N/2),mag(1:N/2));xlabel(頻率/Hz);ylabel(振幅);title(2PSK信號頻譜(N=1024));gridon;fs=250;N=1024;n=0:N-1;t=n/fs;y=ft(y1,N);%對信號進行FFT變換mag=abs(y);%求取Fourier變換的振幅subplot(212)plot(f(1:N/2),mag(1:N/2));xlabel(頻率/Hz);ylabel(振幅);title(接收機接收到的信號頻譜(N=1024));gridon;fs=250;N=1024;n=0:N-1;t=n/fs;y=m(y2,N);%對信號進行FFT變換mag=abs(y);%求取Fourier變換的振幅plot(f(1:N/2),mag(1:N/2));xlabel(頻率/Hz);ylabel(振幅);title(解調(diào)輸出信號頻譜(N=1024)');gridon;fs=250;N=1024;n=0:N-1;t=n/fs;y=m(y3,N);%對信號進行FFT變換mag=abs(y);%求取Fourier變換的振幅plot(f(1:N/2),mag(1:N/2));xlabel(頻率/Hz);ylabel(振幅);title('低通輸出信號頻譜(N=1024));gridon;fs=250;N=1024;n=0:N-1;t=n/fs;y=ft(y5,N);%對信號進行FFT變換mag-abs(y);%求取Fourier變換的振幅plot(f(1:N/2),mag(1:N/2));xlabel(頻率/Hz);ylabel(振幅);title(抽樣判決輸出信號頻譜(N=1024));gridon;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%功能實現(xiàn)模塊九：各階段信號功率譜密度分析%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Fs=1000:n=0:1/Fs:1;nfft=1024;window=hamming(100);%海明窗noverlap=20;%數(shù)據(jù)無重疊range=