數(shù)字調(diào)制的系統(tǒng)設(shè)計(jì)

1、一、設(shè)計(jì)的目的與要求數(shù)字調(diào)制是通信系統(tǒng)中最為重要的環(huán)節(jié)之一，數(shù)字調(diào)制技術(shù)的改進(jìn)也是通信系統(tǒng)性能提高的重要途徑。通過分析了解數(shù)字調(diào)制系統(tǒng)的幾種基本調(diào)制解調(diào)方法，Matlab用于仿真，分析和修改，應(yīng)用圖形界面功能GUI能為仿真系統(tǒng)生成一個人機(jī)交互界面，便于仿真系統(tǒng)的操作，因此采用Matlab對數(shù)字系統(tǒng)進(jìn)行仿真。通過對系統(tǒng)的仿真，我們可以更加直觀的了解數(shù)字調(diào)制系統(tǒng)的性能及影響性能的因素，從而便于改進(jìn)系統(tǒng)，獲得更佳的傳輸性能。二、系統(tǒng)設(shè)計(jì)思路與原理1.基礎(chǔ)知識1.1數(shù)字調(diào)制的意義 數(shù)字調(diào)制是指用數(shù)字基帶信號對載波的某些參量進(jìn)行控制，使載波的這些參量隨基帶信號的變化而變化。根據(jù)控制的載波參量的

2、不同，數(shù)字調(diào)制有調(diào)幅、調(diào)相和調(diào)頻三種基本形式，并可以派生出多種其他形式。由于傳輸失真、傳輸損耗以及保證帶內(nèi)特性的原因，基帶信號不適合在各種信道上進(jìn)行長距離傳輸。為了進(jìn)行長途傳輸，必須對數(shù)字信號進(jìn)行載波調(diào)制，將信號頻譜搬移到高頻處才能在信道中傳輸。因此，大部分現(xiàn)代通信系統(tǒng)都使用數(shù)字調(diào)制技術(shù)。因此，對數(shù)字通信系統(tǒng)的分析與研究越來越重要，數(shù)字調(diào)制作為數(shù)字通信系統(tǒng)的重要部分之一，對它的研究也是有必要的。 1.2 Matlab在通信系統(tǒng)仿真中的應(yīng)用 Matlab是一種交互式的、以矩陣為基礎(chǔ)的軟件開發(fā)環(huán)境,它用于科學(xué)和工程的計(jì)算與可視化。Matlab的編程功能簡單,并且很容易擴(kuò)展和創(chuàng)造

3、新的命令與函數(shù)。應(yīng)用Matlab可方便地解決復(fù)雜數(shù)值計(jì)算問題。Matlab具有強(qiáng)大的Simulink動態(tài)仿真環(huán)境,可以實(shí)現(xiàn)可視化建模和多工作環(huán)境間文件互用和數(shù)據(jù)交換。用戶可以在Matlab和Simulink兩種環(huán)境下對自己的模型進(jìn)行仿真、分析和修改。用于實(shí)現(xiàn)通信仿真的通信工具包(Communication toolbox,也叫Commlib,通信工具箱)是Matlab語言中的一個科學(xué)性工具包,提供通信領(lǐng)域中計(jì)算、研究模擬發(fā)展、系統(tǒng)設(shè)計(jì)和分析的功能,可以在Matlab環(huán)境下獨(dú)立使用,也可以配合Simulink使用。另外Matlab的圖形界面功能GUI（Graphical Us

4、er Interface）能為仿真系統(tǒng)生成一個人機(jī)交互界面，便于仿真系統(tǒng)的操作。因此，Matlab在通信系統(tǒng)仿真中得到了廣泛應(yīng)用。2.系統(tǒng)設(shè)計(jì)原理數(shù)字調(diào)制系統(tǒng)的相關(guān)原理 數(shù)字調(diào)制可以分為二進(jìn)制調(diào)制和多進(jìn)制調(diào)制，多進(jìn)制調(diào)制是二進(jìn)制調(diào)制的推廣，所以本文主要討論二進(jìn)制的調(diào)制與解調(diào)，最后簡單討論一下多進(jìn)制調(diào)制中的差分相位鍵控調(diào)制（M-DPSK）。最常見的二進(jìn)制數(shù)字調(diào)制方式有二進(jìn)制振幅鍵控（2-ASK）、移頻鍵控（2-FSK）和移相鍵控（2-PSK和2-DPSK）。下面是這幾種調(diào)制方式的相關(guān)原理。 2.1 二進(jìn)制幅度鍵控（2-ASK）幅度鍵控可以通過乘法器和開關(guān)電

5、路來實(shí)現(xiàn)。載波在數(shù)字信號1或0的控制下通或斷，在信號為1的狀態(tài)載波接通，此時傳輸信道上有載波出現(xiàn)；在信號為0的狀態(tài)下，載波被關(guān)斷，此時傳輸信道上無載波傳送。那么在接收端我們就可以根據(jù)載波的有無還原出數(shù)字信號的1和0。2-ASK信號功率譜密度的特點(diǎn)如下：  （1）由連續(xù)譜和離散譜兩部分構(gòu)成；連續(xù)譜由傳號的波形g(t)經(jīng)線性調(diào)制后決定，離散譜由載波分量決定； （2）已調(diào)信號的帶寬是基帶脈沖波形帶寬的二倍。2.2 二進(jìn)制頻移鍵控（2-FSK） 頻移鍵控是利用兩個不同頻率f1和f2的振蕩源來代表信號1和0，用數(shù)字信號的1和0去控制兩個獨(dú)立的振蕩源交替輸

6、出。對二進(jìn)制的頻移鍵控調(diào)制方式，其有效帶寬為B=2xF+2Fb,xF是二進(jìn)制基帶信號的帶寬也是FSK信號的最大頻偏，由于數(shù)字信號的帶寬即Fb值大，所以二進(jìn)制頻移鍵控的信號帶寬B較大，頻帶利用率小。2-FSK功率譜密度的特點(diǎn)如下： (1) 2FSK信號的功率譜由連續(xù)譜和離散譜兩部分構(gòu)成,離散譜出現(xiàn)在f1和f2位置；(2) 功率譜密度中的連續(xù)譜部分一般出現(xiàn)雙峰。若兩個載頻之差|f1 -f2|fs，則出現(xiàn)單峰。2.3二進(jìn)制相移鍵控（2-PSK） 在相移鍵控中，載波相位受數(shù)字基帶信號的控制，如在二進(jìn)制基帶信號中為0時，載波相位為0或，為1時載波相位為或

7、0。載波相位和基帶信號有一一對應(yīng)的關(guān)系，從而達(dá)到調(diào)制的目的。2-PSK信號的功率密度有如下特點(diǎn)：  (1) 由連續(xù)譜與離散譜兩部分組成；(2) 帶寬是絕對脈沖序列的二倍；  (3) 與2ASK功率譜的區(qū)別是當(dāng)P1/2時，2PSK無離散譜，而2ASK存在離散譜。2.4 多進(jìn)制數(shù)字調(diào)制 上面所討論的都是在二進(jìn)制數(shù)字基帶信號的情況，在實(shí)際應(yīng)用中，我們常常用一種稱為多進(jìn)制（如4進(jìn)制，8進(jìn)制，16進(jìn)制等）的基帶信號。多進(jìn)制數(shù)字調(diào)制載波參數(shù)有M種不同的取值，多進(jìn)制數(shù)字調(diào)制比二進(jìn)制數(shù)字調(diào)制有兩個突出的優(yōu)點(diǎn)：一是有于多進(jìn)

8、制數(shù)字信號含有更多的信息使頻帶利用率更高；二是在相同的信息速率下持續(xù)時間長，可以提高碼元的能量，從而減小由于信道特性引起的碼間干擾?，F(xiàn)實(shí)中用得最多的一種調(diào)制方式是多進(jìn)制相移鍵控（MPSK）。 多進(jìn)制相移鍵控又稱為多相制，因?yàn)榛鶐盘栍蠱種不同的狀態(tài)，所以它的載波相位有M種不同的取值，這些取值一般為等間隔。在多相制移鍵控有絕對移相和相對移相兩種，實(shí)際中大多采用四相絕對移相鍵控（4PSK，有稱QPSK），四相制的相位有0、/2、3/2四種，分別對應(yīng)四種狀態(tài)11、01、00、10。三、軟件實(shí)現(xiàn)過程與結(jié)果分析1.2-ASK通常，二進(jìn)制振幅鍵控信號（2-ASK）的產(chǎn)生方法（調(diào)制方法）有兩種，如

9、圖3.1所示：圖3.1 2-ASK信號產(chǎn)生的兩種方法2-ASK解調(diào)的方法也有兩種相應(yīng)的接收系統(tǒng)組成方框如圖3.2所示：圖3.2 2-ASK信號接收系統(tǒng)組成框圖根據(jù)3.1（a）所示方框圖產(chǎn)生2-ASK信號，并用圖3.2（b）所示的相干解調(diào)法來解調(diào)，設(shè)計(jì)2-ASK仿真模型如圖3.3所示：圖3.3 2-ASK模型在該模型中，調(diào)制和解調(diào)使用了同一個載波，目的是為了保證相干解調(diào)的同頻同相，雖然這在實(shí)際運(yùn)用中是不可能實(shí)現(xiàn)的，但是作為仿真，這樣能獲得更理想的結(jié)果。 主要模塊參數(shù)設(shè)置如下： 1Bernoulli Binary Generator的參數(shù)設(shè)置為： 

10、Probability of a zero ：0.5 Initial seed ：67 Sample time：1 2. 載波頻率設(shè)為：50（可調(diào)） 3. Sample and Decide模塊是一個子系統(tǒng)，其內(nèi)部結(jié)構(gòu)由抽樣和判決兩部分組成，其中，抽樣由同步?jīng)_激信號（Sychronizing signal）完成，其參數(shù)perio（sec）設(shè)置和信號源的參數(shù)Sample time保持一致。Sample and 

11、;Decide 模塊內(nèi)部結(jié)構(gòu)如圖3.4所示：圖3.4 Sample and Decide子系統(tǒng)內(nèi)部結(jié)構(gòu)4. Error Rate Calculation的參數(shù)設(shè)置：   Receive delay：2 Output data： Work space Variable name：ErrorVec仿真結(jié)果時域分析 設(shè)信息源發(fā)出的是由二進(jìn)制符號0、1組成的序列，且假定0符號出現(xiàn)的概率為P，1符號出現(xiàn)的概率為1-P，他們彼此獨(dú)立。則2ASK信號的時間表示

12、式為: （3.1）（3.2）將圖3.3中各示波器的值輸出到Work space中做統(tǒng)一處理（處理程序見附錄2），各環(huán)節(jié)波形如圖3.5所示；圖3.5 2-ASK各環(huán)節(jié)波形示意圖從圖3.5中可以看出，經(jīng)過調(diào)制后的信號波形在符號1持續(xù)時間內(nèi)是載波的波形，在符號0持續(xù)時間內(nèi)無波形，這與式（3.2）是完全吻合的。最后經(jīng)過解調(diào)和抽樣判決出來的信號與源信號波形大體一致，只是有兩個碼元的延遲，這說明如果將Error Rate Calculation的Receive delay參數(shù)設(shè)置為2，則此模型最后的誤碼率為0。這個值與理論值有些出入，原因是我們在仿真時為了便于觀察信

13、號的波形，將信號源發(fā)送的碼元數(shù)設(shè)定為20個（碼元速率為1，仿真時間20秒），這大大低于現(xiàn)實(shí)中的傳碼率，所以在只傳送20個碼元的情況下，誤碼率為0是可能的。 仿真結(jié)果頻域分析 由于二進(jìn)制的隨機(jī)脈沖序列是一個隨機(jī)過程,所以調(diào)制后的二進(jìn)制數(shù)字信號也是一個隨機(jī)過程，因此在頻率域中只能用功率譜密度表示。2ASK的功率譜密度為：（3.3）又因?yàn)镚(f)=TsSa(fTs), 當(dāng)概率P=0.5時，2ASK的功率譜密度可進(jìn)一步整理為:（3.4）由上式可知，2-ASK信號的中心頻譜被搬移到了載波頻率fc上。對圖3.5中各環(huán)節(jié)數(shù)據(jù)做1024點(diǎn)快速傅立葉可得頻域波形，如圖3.6所示：圖3.6 2

14、-ASK各環(huán)節(jié)頻譜圖從圖3.6中可以看到，源信號中心頻率經(jīng)調(diào)制后搬移到了載波頻率上，這與公式（3.4）是相符的。最后經(jīng)過抽樣判決后的頻譜與源信號頻譜也大體一致，說明該2-ASK仿真模型是成功的、符合理論的。2.2-FSK如果信號源同2-ASK一樣的假設(shè)，那么，2-FSK信號便是0符號對應(yīng)于載波1，而1符號則對應(yīng)于2（與1不同的另一載波）的已調(diào)波形，而且1與2之間的改變是瞬間完成的。2-FSK信號的產(chǎn)生如圖3.7所示：圖3.7 2-FSK信號產(chǎn)生方法2-FSK信號最常用的解調(diào)方法是采用的相干檢測法，如圖3.8所示：圖3.8 2-FSK相干解調(diào)的方法利用Simulink通信工具箱中的FSK調(diào)制解調(diào)

15、模塊及信號源與信道即可。設(shè)計(jì)的2-FSK仿真模型如圖3.9圖3.9 2-FSK仿真模型重要模塊參數(shù)設(shè)置如下：     1.信號源參數(shù)設(shè)置同2-ASK； 2. M-FSK Modulator Passband及M-FSK Demodulator Passband：M-ary number：2 Symbol period (s)：1（與Bernoulli Binary Generator/ Sample tim

16、e一致）Frequency separation (Hz)：1（可調(diào)）Carrier frequency (Hz)：30（可調(diào)）Carrier initial phase (rad)：0   Input sample time(s)：1/100 Output sample time(s)：1/1003. AWGN Channel： Initial seed：120（與Bernoulli Binar

17、y Generator/ Initial seed不同）；  Mode： SNR（dB）SNR (dB)：10（可調(diào)）4. Error Rate Calculation的參數(shù)設(shè)置：  Receive delay：3 Output data： Work space Variable name：ErrorVec2 仿真結(jié)果時域分析：根據(jù)上述2-FSK信號產(chǎn)生原理，已調(diào)信號的時間表達(dá)式可表示為：由式（3

18、.5）可看出2-FSK信號是由兩個2-ASK信號相加而成的，將圖3.9中各示波器的值輸出到Work space中做統(tǒng)一處理（處理程序見附錄3），其中源信號、調(diào)制后信號及解調(diào)后信號波形如圖3.10所示：(a)源信號波形（b）調(diào)制后信號波形（c）解調(diào)后信號波形圖3.10 2-FSK源信號、調(diào)制后信號及解調(diào)后信號波形由圖3.12可知，調(diào)制后信號波形由兩種頻率不同的波形組成，且兩種頻率分別對應(yīng)解調(diào)后信號的符號0和符號1，即2-FSK信號波形可以看作是由兩個2-ASK信號相加而成的，這與式（3.5）完全相符。仿真結(jié)果頻域分析 改變Frequency separation&#

19、160;(Hz)和Carrier frequency (Hz)兩個參數(shù)的值單獨(dú)觀察調(diào)制后的頻譜，獲得圖3.11中的兩個頻譜圖（a）載波差值：1 載波：30（b）載波差值：5 載波：20圖3.11 2-FSK調(diào)制后頻譜對比圖3.13（a）和 (b)可知，當(dāng)兩個載波差值很小時，已調(diào)信號的頻譜呈現(xiàn)單峰如（a）圖；當(dāng)兩個載波差值較大時，已調(diào)信號的頻譜呈現(xiàn)雙峰如（b）圖,這與2.2節(jié)中闡述的2-FSK頻譜的特點(diǎn)完全相符。 仿真結(jié)果的分析說明該2-FSK仿真模型是可行的。3.2-PSK和2-FSK一樣我們也用Simulink通信工具箱提供的現(xiàn)成DPSK調(diào)制解調(diào)模塊來

20、構(gòu)建仿真模型，并由M文件編制程序?qū)Ψ抡娼Y(jié)果進(jìn)行統(tǒng)一處理。2-PSK仿真模型如圖3.12所示：圖3.12 2-PSK仿真模型主要模塊參數(shù)設(shè)置如下： 將信號源的Sample time設(shè)為1/2，仿真觀察40個碼元，調(diào)制解調(diào)模塊中的Symbol period (s)也相應(yīng)設(shè)成1/2，其余參數(shù)可參照2-FSK，兩者參數(shù)類似。仿真結(jié)果分析：2-PSK信號的時間表達(dá)式為：若在某一碼元持續(xù)時間Ts內(nèi)觀察時，可以簡寫為： （3.7）或以相反的形式。將圖3.12中各示波器數(shù)據(jù)做統(tǒng)一處理，得到各環(huán)節(jié)時域頻域?qū)Ρ葓D如圖3.13所示，仿真結(jié)果處理程序見附錄4圖3.13 2-PS

21、K各環(huán)節(jié)時域頻域波形 從上圖可以看出，調(diào)制后的信號波形由兩種相位不同的波形組成，而且兩種波形是反相的，即相位相差180度，這與式（3.7）一致。解調(diào)后的時域波形和源信號相比，不僅有一個碼元的延遲，而且第一個碼元由1變成了0，出現(xiàn)了誤碼，由誤碼計(jì)算儀的計(jì)算數(shù)據(jù)可知，該系統(tǒng)在傳送40個碼元的情況下誤碼率為0.025，這是一個理論上和現(xiàn)實(shí)中都可以接受的值。 將式（3.6）與式（3.2）比較可見，2-ASK和2-PSK時間表達(dá)式形式完全相同，所不同的只是an的取值，因此，兩者的頻域波形也相似。將圖3.13中的2-PSK信號調(diào)制后頻譜與圖3.6中的2-ASK信號調(diào)制后頻譜比較也可得出

22、相同結(jié)論。4. GUI用戶界面設(shè)計(jì)為了使整個數(shù)字調(diào)制仿真系統(tǒng)便于操作，我們利用Matlab 中的GUI圖形用戶界面將各種調(diào)制模型集合到一個操作界面上，使該仿真系統(tǒng)成為一個整體。 GUI全稱Graphical User Interface即，圖形用戶界面，就是通過窗口、選單、按鈕、文字說明等對象構(gòu)成一個美觀的界面，提供用戶利用鼠標(biāo)或鍵盤方便地實(shí)現(xiàn)操作。GUI是Matlab自帶的用于生成人機(jī)交互界面的工具，Matlab設(shè)計(jì)圖形用戶界面的方法有兩種：使用可視化的界面環(huán)境和通過編寫程序。用戶可根據(jù)自己的需要生成符合用戶要求的圖形界面，并通過編寫各控件的回調(diào)函數(shù)（

23、callback function）使控件和該控件所要實(shí)現(xiàn)的功能聯(lián)系起來，從而實(shí)現(xiàn)圖形界面的操作功能。圖3.14 數(shù)字調(diào)制仿真系統(tǒng)操作界面 只要在Matlab命令窗口輸入仿真系統(tǒng)的名稱：ddmod即可進(jìn)入如上圖所示的仿真操作界面。從上圖可以看出，該仿真系統(tǒng)操作界面主要由時間域波形顯示框、頻率波形顯示框、調(diào)制方式單選項(xiàng)、設(shè)置參數(shù)按鈕、打開模型按鈕、仿真按鈕、關(guān)閉按鈕及包含所要觀察的各環(huán)節(jié)的下拉式菜單等控件組成。通過該操作界面，我們可以使各個仿真模型及相關(guān)函數(shù)在后臺運(yùn)行，同時將我們需要的結(jié)果在界面上顯示出來，使系統(tǒng)的操作方便簡潔。另一方面，可以通過“參數(shù)設(shè)置”按鈕對仿真模型的主

24、要參數(shù)進(jìn)行設(shè)置，以觀察在不同參數(shù)下仿真的結(jié)果，“參數(shù)設(shè)置”對話框如圖3.15所示：圖3.15 M-DPSK參數(shù)設(shè)置對話框如果需要對更多的參數(shù)進(jìn)行設(shè)置或者需要修改仿真模型的某些模塊，可以通過單擊“打開模型”按鈕來打開所需要的模型。當(dāng)各種參數(shù)設(shè)置好后，可以在下拉式菜單中選擇一個要觀察的仿真環(huán)節(jié)，然后單擊仿真按鈕就可以在時間域波形顯示框、頻率波形顯示框中觀察到時間域波形和頻率波形，當(dāng)然，這一切的前提是在調(diào)制方式單選項(xiàng)上選擇了其中一項(xiàng)，否則將會出現(xiàn)如圖3.16所示的警告：圖3.16 警告對話框另外，為了方便查看與仿真相關(guān)的程序，F(xiàn)ile中菜單中設(shè)計(jì)了OpenMfile一項(xiàng)，如圖3.17所示，圖3.17

25、 File菜單 仿真結(jié)束后可單擊“關(guān)閉”按鈕或File/Close關(guān)閉仿真系統(tǒng)。四、設(shè)計(jì)總結(jié)與心得體會課程設(shè)計(jì)是培養(yǎng)學(xué)生綜合運(yùn)用所學(xué)的理論知識，發(fā)現(xiàn)、提出、分析和解決實(shí)際問題，鍛煉實(shí)踐能力的重要環(huán)節(jié)，是對學(xué)生實(shí)際工作能力的具體訓(xùn)練和考察過程。通過動手去做，理論與實(shí)踐相結(jié)合使我更加明白什么是數(shù)字調(diào)制系統(tǒng)。通過兩周的的努力，我們的課程設(shè)計(jì)終于完成了，這意味著這一階段艱苦的熱烈的付出即將告一段落。在整個課程設(shè)計(jì)中，我在學(xué)習(xí)上和思想上都有很大收獲，這除了自身的努力外，與指導(dǎo)的老師、同學(xué)和朋友的關(guān)心、支持和鼓勵是分不開。理論應(yīng)用到實(shí)際是需要一個過程，遇到自己解決不了的問題，通過查閱資料，請教同學(xué)和請教老

26、師等方法來解決。而且可以發(fā)現(xiàn)自己的不足之處，對以前所學(xué)知識進(jìn)一步掌握，也可以通過課程設(shè)計(jì)過程提高實(shí)際工作能力和獨(dú)立思考的能力。由于所學(xué)知識水平和能力的有限，在課程設(shè)計(jì)過程中，我遇到了很多問題，有些通過查閱相關(guān)資料得到了解決，也有一部分是請教同學(xué)，而更重要的是老師的答疑給我們一個系統(tǒng)的框架和結(jié)構(gòu)，使我們在總體方案的設(shè)計(jì)中有了正確的方法和方向。本次課程設(shè)計(jì)按是以小組為單位進(jìn)行的，小組中的每個人既有不同的任務(wù)也可以通過小組討論來提出更好的解決辦法和不同的思路。在這個過程中，獨(dú)立思考很重要，團(tuán)隊(duì)合作也是必不可少的?？傊舜握n程設(shè)計(jì)使我受益匪淺！參考文獻(xiàn) 1 樊昌信. 通

27、信原理教程.北京：電子工業(yè)出版社，2005 2 張輝，曹麗娜現(xiàn)代通信原理與技術(shù) M西安電子科學(xué)出版社，2001 3 張威MATLAB入門與基礎(chǔ) M西安電子科技大學(xué)出版社，2004 4 施曉紅，周佳精通GUI圖形界面編程 M北京大學(xué)出版社，2003 5 JDGibsonPrinciples of Digital and Analog Communications1990 6 F Xiong， Dig

28、ital Modulation TechniquesArtech House， 2000  7 鄧華MATLAB通信仿真及應(yīng)用實(shí)例詳解 M北京：人郵電出版社，20038 李建新現(xiàn)代通信系統(tǒng)分析與仿真一MATLAB通信工具箱 M.西安電子科技大學(xué)出版社，2000附錄附錄1  2-ASK仿真結(jié)果處理程序function ask1() sim('dm'); t1=ScopeData5(:,1); y1=ScopeData5(:,2);&

29、#160;t2=ScopeData2(:,1); y2=ScopeData2(:,2);t3=ScopeData(:,1);y3=ScopeData(:,2);t5=ScopeData1(:,1);y5=ScopeData1(:,2); t6=ScopeData7(:,1);y6=ScopeData7(:,2); t7=ScopeData3(:,1); y7=ScopeData3(:,2);   %引入示波器的數(shù)據(jù)subplot(3,2,5); plot(t1,y1),grid; xlabel('

30、濾波后的波形')subplot(3,2,4);plot(t2,y2),grid; xlabel('相干解調(diào)后的波形')subplot(3,2,1); stairs(t3,y3),axis(0 20 -0.5 2.5),grid;title('ASK調(diào)制與解調(diào)各環(huán)節(jié)波形') xlabel('源信號')subplot(3,2,2); plot(t5,y5),grid; xlabel('調(diào)制后的波形')subplot(3,2,3);plot(t6,y

31、6),grid; xlabel('經(jīng)過高斯信道后的波形') subplot(3,2,6); stairs(t7,y7),axis(0 20 -0.5 2.5),grid;xlabel('抽樣判決后的波形') figure(2) n1=1024; f1=(-n1/2:1:n1/2-1);s=fft(y3,n1); subplot(3,2,1); plot(f1,abs(fftshift(s),%axis(-30 30 0 7)

32、xlabel('源信號頻譜') n2=1024; f2=(-n2/2:1:n2/2-1); s=fft(y5,n2); subplot(3,2,2); plot(f2,abs(fftshift(s),%axis(-400 400 0 300) xlabel('調(diào)制后的頻譜') s=fft(y2,n2);subplot(3,2,3); plot(f2,abs(fftshift(s),%axis(-400 400 0 300)x

33、label('相干解調(diào)后的頻譜')s=fft(y1,n2); subplot(3,2,4); plot(f2,abs(fftshift(s),%axis(-400 400 0 300) xlabel('濾波后的頻譜') s=fft(y7,n2);subplot(3,2,5); plot(f2,abs(fftshift(s),%axis(-400 400 0 300)xlabel('抽樣判決后的頻譜')附錄2  2-FSK仿真結(jié)果

34、處理程序 function bfsk2() sim('bfsk'); t1=ScopeData(:,1);y1=ScopeData(:,2);t2=ScopeData1(:,1);y2=ScopeData1(:,2); t3=ScopeData2(:,1); y3=ScopeData2(:,2);t4=ScopeData3(:,1); y4=ScopeData3(:,2);subplot(4,2,1); stairs(t1,y1),axis(0 20 -1/2 3/2)

35、,grid;xlabel('源信號')subplot(4,2,3); plot(t2,y2),grid; xlabel('調(diào)制后的波形') subplot(4,2,5); plot(t3,y3),grid; xlabel('經(jīng)過高斯信道后的波形')subplot(4,2,7); stairs(t4,y4),axis(0 20 -1/2 3/2),grid; xlabel('相干解調(diào)后的波形') n1=1024; 

36、f1=(-n1/2:1:n1/2-1);s=fft(y1,n1); subplot(4,2,2); plot(f1,abs(fftshift(s),%axis(-30 30 0 7)xlabel('源信號頻譜')n2=1024; f2=(-n2/2:1:n2/2-1);s=fft(y2,n2); subplot(4,2,4); plot(f2,abs(fftshift(s),%axis(-400 400 0 300) xlabel('調(diào)制后的頻譜

37、9;) f2=(-n2/2:1:n2/2-1);s=fft(y3,n2); subplot(4,2,6); plot(f2,abs(fftshift(s),%axis(-400 400 0 300)xlabel('經(jīng)過高斯信道后的頻譜')f2=(-n1/2:1:n1/2-1); s=fft(y4,n1); subplot(4,2,8); plot(f2,abs(fftshift(s),%axis(-400 400 0 300) xlabel('調(diào)制后的頻譜')r=ErrorVec2(1); w=sprintf('符號誤差率= %dn',r); uiwait(msgbox(w,'符號誤差率','modal');附錄3 2-PSK仿真結(jié)果處理程序 function psk() sim('msk');