--正文--基于MATLAB環(huán)境下的GUI的FSK仿真演示系統(tǒng)

1、第一章 緒 論1.1課題研究背景及意義隨著社會(huì)的不斷發(fā)展，通信對(duì)我們來說越來越顯的重要。對(duì)于通信技術(shù)來講，信息傳輸?shù)挠行院涂煽啃允桥卸ㄐ阅艿臉?biāo)準(zhǔn)。而對(duì)于保證遠(yuǎn)距離傳輸信息的正確性這一方面，數(shù)字通信系統(tǒng)具有先天的優(yōu)勢(shì)。在數(shù)字載波通信中，我們采用了三種解調(diào)方式:幅移鍵控(ASK)、頻移鍵控(FSK)、相移鍵控(PSK)。其中頻移鍵控(FSK)是信息傳輸中使用得較早的一種調(diào)制方式，頻移鍵控的最簡(jiǎn)單形式是二進(jìn)制頻率鍵控2FSK。它的主要優(yōu)點(diǎn)是：實(shí)現(xiàn)方法簡(jiǎn)單，抗噪聲與抗衰減性能較好。因此在中低速數(shù)據(jù)傳輸中得到了廣泛的應(yīng)用。本課題研究FSK，并基于MATLAB環(huán)境下的GUI對(duì)FSK系統(tǒng)進(jìn)行仿真演示系統(tǒng)的

2、制作，從而通過運(yùn)用模擬的視覺化的手段來實(shí)現(xiàn)達(dá)到解調(diào)調(diào)制的目的。隨著電子通信科技領(lǐng)域的高速發(fā)展，通信系統(tǒng)復(fù)雜性隨之增加。而傳統(tǒng)的手工分析與電路板試驗(yàn)等分析設(shè)計(jì)方法己經(jīng)不能適應(yīng)發(fā)展的需要，通信系統(tǒng)計(jì)算機(jī)模擬仿真技術(shù)日益顯示出其巨大的優(yōu)越性。計(jì)算機(jī)仿真是根據(jù)被研究的真實(shí)系統(tǒng)的模型，利用計(jì)算機(jī)進(jìn)行實(shí)驗(yàn)研究的一種方法。MATLAB是矩陣實(shí)驗(yàn)室(Matrix Laboratory)的簡(jiǎn)稱，它在數(shù)學(xué)類科技應(yīng)用軟件中在數(shù)值計(jì)算方面首屈一指。MATLAB可以進(jìn)行矩陣運(yùn)算、繪制函數(shù)和數(shù)據(jù)、實(shí)現(xiàn)算法、創(chuàng)建用戶界面等，主要應(yīng)用于工程計(jì)算、控制設(shè)計(jì)、信號(hào)處理與通訊、圖像處理、信號(hào)檢測(cè)等領(lǐng)域。同時(shí)，Matlab仿真軟件也

3、是分析通信系統(tǒng)常用的工具之一。Matlab的圖形界面功能GUI (Graphical User Interface)能為仿真系統(tǒng)生成一個(gè)人機(jī)交互界面，便于仿真系統(tǒng)的操作。因此，GUI在通信系統(tǒng)的各種仿真中得到了廣泛的應(yīng)用，本文也選用該工具對(duì)數(shù)字調(diào)制系統(tǒng)進(jìn)行仿真。1.2課題設(shè)計(jì)要求設(shè)計(jì)一個(gè)“基于MATLAB環(huán)境下GUI的FSK仿真演示系統(tǒng)”，具體要求如下：1掌握FSK的原理和仿真方法；2通過腳本編程或者SIMULINK對(duì)FSK調(diào)制方式進(jìn)行仿真；3在MATLAB環(huán)境下制作一個(gè)GUI，可以由用戶輸入載波頻率、兩個(gè)調(diào)制頻率f1和f2，動(dòng)態(tài)產(chǎn)生FSK調(diào)制波形；4觀察調(diào)制信號(hào)和已調(diào)信號(hào)波形，分析調(diào)制性能。

4、第二章 系統(tǒng)組成及工作原理2.1 系統(tǒng)單元功能模塊設(shè)計(jì)此系統(tǒng)必須具備以下單元功能模塊： 原序列信號(hào)產(chǎn)生模塊； 信號(hào)加躁模塊； 調(diào)制模塊； 1、帶通濾波器模塊； 解調(diào)模塊： 2、乘法器模塊； 3、低通濾波器模塊； 4、抽樣判決模塊；2.2 數(shù)字調(diào)制解調(diào)原理數(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)制數(shù)字調(diào)制方式有二進(jìn)制振幅鍵控(2-ASK)、移頻鍵控(2-FSK)和移相鍵控(2-PSK和2-DPSK)。在此基礎(chǔ)上還可以派生出許多其他的形式。數(shù)字調(diào)制是用載波信號(hào)的某些離散狀態(tài)來表征所傳輸?shù)男畔?，在接收端也?duì)載波信號(hào)的

5、離散調(diào)制參量進(jìn)行檢測(cè)。數(shù)字調(diào)制采用正弦波調(diào)制，即信號(hào)被調(diào)制為高頻正弦波。 數(shù)字解調(diào)是一種為了傳播方便，把信息編碼傳輸?shù)姆椒āＴ趥鬏斶^程中，語音通話等會(huì)轉(zhuǎn)變成一連串的數(shù)字信號(hào)，變?yōu)橛?jì)算機(jī)二進(jìn)制代碼0和1，在信息接收處又解碼恢復(fù)到原來的語音狀態(tài)。解調(diào)是調(diào)制的逆過程，其作用是從接受的己調(diào)信號(hào)中恢復(fù)原基帶信號(hào)(即調(diào)制信號(hào))。調(diào)制的方法可分為兩類：相干解調(diào)和非相干解調(diào)。由于FSK的解調(diào)過程有相干解調(diào)法和非相干解調(diào)法，我們?cè)谶@里只介紹相干解調(diào)原理。相干解調(diào)是指利用乘法器，輸入一路與載頻相干(同頻同相)的參考信號(hào)與載頻相乘。2.3 頻移鍵控FSK原理二進(jìn)制數(shù)字頻率調(diào)制方法有模擬調(diào)制和鍵控法兩種。其中二進(jìn)制頻

6、率鍵控2FSK (frequency-shift keying)是采用最為廣泛的一種，它原理更為簡(jiǎn)單、實(shí)現(xiàn)起來也更為容易、抗噪聲和抗衰減性能好、穩(wěn)定可靠，是中低速數(shù)據(jù)傳輸?shù)淖罴堰x擇。故一般采用鍵控法來實(shí)現(xiàn)二進(jìn)制的頻率鍵控2FSK。2.3.1 2FSK調(diào)制原理一、基本原理2FSK是用原二進(jìn)制”0”/”1”序列信號(hào)調(diào)制載波的一種數(shù)字調(diào)制方式。我們將每個(gè)比特的信息轉(zhuǎn)換為一個(gè)頻率，不同的電平對(duì)應(yīng)不同頻率的載波。一般來說，“0”由較低的頻率f0表示，“1”由較高的頻率f1表示。即f0<f1。其表達(dá)式為： （式2-1）對(duì)應(yīng)典型波形為：圖 2.1 調(diào)制典型波形頻移鍵控采用的調(diào)制方式需要兩個(gè)不同頻率段的

7、載波信號(hào)，碼元為“1”和碼元為“0”時(shí)的載波頻率不同，其數(shù)學(xué)表達(dá)式為： （式2-2） （式2-3） （式2-3）二、2FSK信號(hào)的產(chǎn)生方法采用模擬調(diào)頻電路來實(shí)現(xiàn)：信號(hào)在相鄰碼元之間的相位是連續(xù)變化的。采用鍵控法來實(shí)現(xiàn)：相鄰碼元之間的相位不一定連續(xù)。鍵控法的原理框圖如圖2.2所示?；鶐盘?hào)振蕩器0（f0）振蕩器1（f1）0 鍵控開關(guān) 1e2FSK（t）圖 2.2 鍵控法原理框圖基帶信號(hào)是一個(gè)二進(jìn)制0/1信號(hào)。鍵控開關(guān)受基帶信號(hào)控制，其中基帶信號(hào)為“1”時(shí)，開關(guān)接至振蕩器1。反之，基帶信號(hào)為“0”時(shí)，開關(guān)接至振蕩器0。當(dāng)振蕩器產(chǎn)生信號(hào)為頻率分別為f0、f1的余弦信號(hào)時(shí)，調(diào)制信號(hào)即為 。其中，都是二

8、進(jìn)制序列信號(hào)。2.3.2 2FSK解調(diào)原理解調(diào)原理框圖如圖2.4所示，對(duì)應(yīng)的典型波形圖如圖2.3所示。圖2.3 解調(diào)典型波形圖帶通w0帶通w1低通低通e2fskcos w0cos w1定時(shí)脈沖抽樣判決S（t）圖2.4 解調(diào)原理框圖第三章 系統(tǒng)設(shè)計(jì)思路3.1 調(diào)制模塊設(shè)計(jì)開 始輸入信息碼w、載頻f0、f1碼長le=length(w)；設(shè)步長為0.001，每個(gè)碼元有50個(gè)步長；C1=sin(2f1t)；C0=sin(2f0t)；i=1, j=1基帶采樣信號(hào)是否為0？n=(i-1)*50+jbit(n)=1，fsk(n)= C1 (n)bit(n)=0，fsk(n)= C0 (n)j=j+1j=50

9、?i =i+1i=le?fftfsk=FTfsk, nos=awgn(fsk,20),ffsk=FTnos將bit、fsk、fftfsk、nos、ffsk畫出結(jié) 束否是否是否是圖3.1 調(diào)制模塊設(shè)計(jì)流程圖注釋：w為二進(jìn)制信息碼；bit信號(hào)是承載碼信息的二進(jìn)制單極性不歸零脈沖序列；fsk是已調(diào)制的二進(jìn)制頻移鍵控信號(hào)；nos是加噪后的fsk信號(hào)；fftfsk、ffsk分別是fsk、nos信號(hào)的傅立葉變換，即頻域函數(shù)。3.2 解調(diào)模塊設(shè)計(jì)圖3.2 解調(diào)模塊設(shè)計(jì)流程圖開 始經(jīng)過帶通濾波器0，得到以f0為中心頻率的y0(t)，及其頻譜函數(shù)k0(f)輸入nos、ffsk經(jīng)過帶通濾波器1，得到以f1為中心頻

10、率的y1(t)，及其頻譜函數(shù)k1(f)將y0與C0相乘，得到m0(t)將y1與C1相乘，得到m1(t)經(jīng)過低通濾波器0，得到含直流分量的Ym0(t)，及其頻譜函數(shù)ym0(f)經(jīng)過低通濾波器1，得到含直流分量的Ym1(t)，及其頻譜函數(shù)ym1(f)i=1Ym0(i)>Ym1(i)？g(i)=0g(i)=1將解調(diào)輸出信號(hào)g畫出結(jié) 束i=length(Ym0)？i=i+1是否否是3.3 GUI界面設(shè)計(jì)如圖3.3所示，界面一共有三個(gè)編輯框edit1edit3，兩個(gè)下拉菜單popupmenu3popupmenu4。一個(gè)pushbutton20。三個(gè)axes1axes2。他們的功能描述如表3-1所示

11、：圖3.3 GUI界面設(shè)計(jì)表3-1 GUI界面功能描述控件名功能描述edit1輸入二進(jìn)制信息碼wedit2輸入載頻f0edit3輸入載頻f1popupmenu3選擇axes2坐標(biāo)軸要顯示的圖形popupmenu4選擇axes3坐標(biāo)軸要顯示的圖形pushbutton20清除axes2、axes3坐標(biāo)軸上的圖形axes1顯示背景圖片axes2顯示popupmenu3中被選擇的函數(shù)的圖形axes3顯示popupmenu4中被選擇的函數(shù)的圖形第四章 系統(tǒng)調(diào)試結(jié)果分析4.1 調(diào)試步驟及現(xiàn)象分析4.1.1 調(diào)試步驟第一步：運(yùn)行m文件或者fig文件后，系統(tǒng)會(huì)打開“基于MATLAB環(huán)境下的GUI的FSK調(diào)制解

12、調(diào)仿真演示系統(tǒng)”。第二步：輸入任意長度的二進(jìn)制信息碼w，載頻f0及f1。（缺省值為f0=400，f1=200，w=1 0 0 0 0 ）第三步：選中任一下拉菜單中的任意選項(xiàng)，即可調(diào)用FSK調(diào)制解調(diào)函數(shù)bfsk，并且根據(jù)選中項(xiàng)的項(xiàng)號(hào)，查找到相應(yīng)已生成的函數(shù)信號(hào)，最后在相應(yīng)的坐標(biāo)軸中顯示出來。第四步：選擇另一下拉菜單中的任意項(xiàng)，并在其相應(yīng)坐標(biāo)軸中顯示所選函數(shù)信號(hào)。第五步：根據(jù)需要，使用清除圖像功能。4.1.2現(xiàn)象分析整個(gè)調(diào)制解調(diào)過程，生成了很多函數(shù)信號(hào)。仿真演示系統(tǒng)可以通過popupmenu選擇任意函數(shù)信號(hào)進(jìn)行對(duì)比觀察。下面按照整個(gè)調(diào)制解調(diào)過程，順序地選擇一些具有重要對(duì)比意義的現(xiàn)象分析系統(tǒng)性能及調(diào)

13、試結(jié)果。一、打開系統(tǒng)界面圖4.1 系統(tǒng)界面圖運(yùn)行m文件或fig文件后，顯示如圖4.1所示系統(tǒng)界面。二、輸入?yún)?shù)、選擇菜單項(xiàng)打開界面后，用戶需要輸入任意長度的信息碼w、載頻f0和f1三個(gè)參數(shù)。接著，如圖4.2所示，點(diǎn)擊下拉菜單后，共有原信號(hào)m(t)、已調(diào)信號(hào)、抽樣判決輸出信號(hào)m1(t)等16項(xiàng)函數(shù)信號(hào)可供選擇。用戶可以選擇不同的函數(shù)信號(hào)分別顯示在圖形一和圖形二坐標(biāo)軸上，這樣可以方便對(duì)比觀察。（其中，編輯框的缺省值為f0=400，f1=200，w=1 0 0 0 0 ）圖4.2 用戶輸入?yún)?shù)三、原信號(hào)&已調(diào)信號(hào)輸入信息碼1 0 0 1 0，載頻f0=50，f1=100。觀察原信號(hào)（圖一）和

14、已調(diào)信號(hào)（圖二）。如圖4.3所示。顯然，對(duì)比觀察圖中兩個(gè)信號(hào)，原信號(hào)為1時(shí)的已調(diào)信號(hào)頻率f1是原信號(hào)為0時(shí)的已調(diào)信號(hào)頻率f0的兩倍，可以確定演示結(jié)果正確。四、已調(diào)信號(hào)&已調(diào)加噪信號(hào)輸入信息碼1 0 0 1 0，載頻f0=50，f1=100。觀察已調(diào)信號(hào)（圖一）和已調(diào)加噪信號(hào)（圖二）。如圖4.4所示。當(dāng)已調(diào)信號(hào)經(jīng)過信道傳輸后，不免會(huì)有噪聲干擾。這里我們調(diào)用了MATLAB中的awgn(f,N)函數(shù)的功能產(chǎn)生噪聲（其中f指的是無噪聲信號(hào)，N指的是信噪比），令N=20。將圖4.4中的兩個(gè)圖形對(duì)比觀察，發(fā)現(xiàn)圖形二所示信號(hào)顯然已被噪聲干擾，發(fā)生了一定的變形，可以確定演示結(jié)果正確。圖4.3 原信號(hào)&

15、amp;已調(diào)信號(hào)圖4.4 已調(diào)信號(hào)&已調(diào)加噪信號(hào)五、已調(diào)加噪信號(hào)&已調(diào)加噪信號(hào)頻譜 輸入信息碼1 0 0 1 0，載頻f0=400，f1=200。觀察已調(diào)加噪信號(hào)（圖一）和已調(diào)加噪信號(hào)頻譜（圖二）。如圖4.5所示。圖4.5 已調(diào)加噪信號(hào)&已調(diào)加噪信號(hào)頻譜六、經(jīng)乘法器與相應(yīng)載波相乘后的信號(hào)m0、m1輸入信息碼1 0 0 1 0，載頻f0=400，f1=200。觀察k0與載波C0=sin(2f0t)相乘得到的信號(hào)m0（圖一）和k1與載波C1=sin(2f1t)相乘得到的信號(hào)m1（圖二）。如圖4.6所示。相乘后，若某點(diǎn)對(duì)應(yīng)的原信號(hào)應(yīng)為0，在此點(diǎn)處信號(hào)平均值>。同理，若某

16、點(diǎn)對(duì)應(yīng)的原信號(hào)應(yīng)為1，則在此點(diǎn)處信號(hào)>。從而下面可以直接用抽樣判決的方法解調(diào)出原信號(hào)。對(duì)比觀察圖中兩個(gè)信號(hào)，可以確定演示結(jié)果正確。七、經(jīng)兩個(gè)低通濾波器后的信號(hào)Ym0、Ym1輸入信息碼1 0 0 1 0，載頻f0=400，f1=200。觀察m0經(jīng)帶通濾波器0濾波后得到的頻域函數(shù)Ym0（圖一）和m1經(jīng)帶通濾波器1濾波后得到的頻域函數(shù)Ym1（圖二）。如圖4.7所示。經(jīng)過低通濾波器以后，我們將經(jīng)乘法器得到的信號(hào)m0、m1的高頻成分濾除，得到較為平緩的信號(hào)Ym0、Ym1。這樣使得抽樣判決器的抽樣值較為準(zhǔn)確，有利于提高抽樣判決器的有效性。觀察圖中兩個(gè)信號(hào)，不難看出在有些時(shí)間段m0>m1，有些時(shí)

17、間段m1>m0，且這些時(shí)間段較為嚴(yán)格地區(qū)分開來了。這恰好符合抽樣判決器的判決要求。所以可以確定演示結(jié)果正確。圖4.6 經(jīng)乘法器與相應(yīng)載波相乘后的信號(hào)m0、m1圖4.7 經(jīng)兩個(gè)低通濾波器后的信號(hào)Ym0、Ym1八、原信號(hào)m(t)&解調(diào)輸出信號(hào)m1(t)輸入信息碼1 0 0 1 0，載頻f0=400，f1=200。觀察原信號(hào)m(t)（圖一）和解調(diào)輸出信號(hào)m1(t)（圖二）。如圖4.8所示。信號(hào)在系統(tǒng)中受調(diào)制、傳輸、解調(diào)等過程是會(huì)存在一定的延遲的。而比較圖中兩個(gè)信號(hào)發(fā)現(xiàn)，m1(t)相對(duì)于m(t)有一定的延遲。所以可以確定演示結(jié)果正確。圖4.8 原信號(hào)m(t)&解調(diào)輸出信號(hào)m1(t

18、)4.2 調(diào)試故障分析系統(tǒng)調(diào)試時(shí)，多次在調(diào)制步驟出現(xiàn)故障，載波的波形出現(xiàn)頻率不正確的現(xiàn)象。經(jīng)過多次試驗(yàn)后，發(fā)現(xiàn)振蕩器的頻率f0、f1需要在200HZ以下，才能夠明顯觀察出一個(gè)碼元內(nèi)的載波個(gè)數(shù)。所以，系統(tǒng)在使用仿真時(shí)，輸入的載波頻率最好是在200HZ以下。第五章 總結(jié)自從因特網(wǎng)把我們領(lǐng)進(jìn)信息時(shí)代開始，人類的歷史翻開了璀璨的一頁。隨著信息的飛速發(fā)展，通信原理也隨之崛起。從而，使得培養(yǎng)新世紀(jì)的技術(shù)人才顯得分外重要。在學(xué)習(xí)通信原理理論基礎(chǔ)后，我們了解到FSK是信息傳輸中使用得較早的一種調(diào)制方式，也是數(shù)字通信中用得較廣的一種方式。它的主要優(yōu)點(diǎn)是：實(shí)現(xiàn)起來較容易，抗噪聲與衰減的性能較好，適合中低速數(shù)據(jù)的傳

19、輸?shù)膽?yīng)用。通常數(shù)據(jù)率在低于1200bps時(shí)使用FSK方式。在衰落信道中傳輸數(shù)據(jù)時(shí)，它也被廣泛采用。相干解調(diào)對(duì)接收設(shè)備的復(fù)雜程度比非相干解調(diào)較高。在理論學(xué)習(xí)的基礎(chǔ)上通過利用MATLAB仿真，真正的看到了通信中傳輸信息的一系列的問題。比如說要使信號(hào)不失真地傳輸?shù)浇邮斩耍鸵紤]很多的因數(shù)。在發(fā)送端要注意噪聲的加入，盡量地減少噪聲浸入信道中，以免在接受端使信號(hào)失真度過大而不能恢復(fù)成原來的信號(hào)。而在接收端，采用哪種解調(diào)方式能夠更好地恢復(fù)出原來的信號(hào)，對(duì)于不同的解調(diào)方式有相干解調(diào)和非相干解調(diào)。相干解調(diào)一般是在接收端使接收的信號(hào)通過一個(gè)相乘器，同時(shí)乘上一個(gè)與原調(diào)制信號(hào)同頻同相的載波，再通過低通濾波器濾除不

20、需要的信號(hào)。然后再經(jīng)過采樣、量化和編碼，最終可以得到原調(diào)制信號(hào)。通過這次的課程設(shè)計(jì)，進(jìn)一步了解了二進(jìn)制頻移鍵控即2FSK的基本原理及其相干解調(diào)方法。當(dāng)然在學(xué)習(xí)的過程中，也遇到了許多困難。比如參數(shù)設(shè)置的不理想時(shí)，總是會(huì)出現(xiàn)波形失真等問題。但是通過查閱資料和請(qǐng)教老師，我們很好地完成了此次設(shè)計(jì)。同時(shí)這次的設(shè)計(jì)也讓我對(duì)MATLAB環(huán)境下的GUI有了一定的了解，使我明白了用計(jì)算機(jī)仿真電子通信系統(tǒng)，具有廣泛的適應(yīng)性和極高的靈活性。心 得本次設(shè)計(jì)主要涉及到了通信原理和MATLAB的相關(guān)知識(shí)與運(yùn)用，主要有基帶信號(hào)的調(diào)制原理及方法、相干解調(diào)的原理及實(shí)現(xiàn)方法、帶通和低通濾波器的特性、抽樣判決的實(shí)現(xiàn)方法等等，加深了

21、對(duì)上述相關(guān)知識(shí)的了解，使自己更深刻理解了調(diào)制與解調(diào)的原理和實(shí)現(xiàn)方法，已經(jīng)基本掌握了MATLAB的m文件及其環(huán)境下的GUI的基本應(yīng)用。通過這次課程設(shè)計(jì)，我們學(xué)到了很多書本上沒有的知識(shí)。鍛煉了我們獨(dú)立思考問題、分析問題、解決問題的能力。而且本次設(shè)計(jì)是和同組人共同完成的，加強(qiáng)了與他人溝通的能力以及團(tuán)隊(duì)合作精神，為今后走向社會(huì)提供了很好的準(zhǔn)備。參考文獻(xiàn)1 樊昌信. 通信原理.北京：國防工業(yè)出版社，2008.2 郝文化. MATLAB圖形圖像處理應(yīng)用教程.北京：中國水利水電出版社，2002.3 徐金明.MATLAB實(shí)用教程.北京：清華大學(xué)出版社，2003.4 達(dá)新宇. 通信原理實(shí)驗(yàn)與課程設(shè)計(jì).北京：北京

22、郵電大學(xué)出版社，2003.附錄一 m文件由于有很多都是系統(tǒng)生成語句，不是函數(shù)關(guān)鍵部分，故這里刪除了許多不必要的語句，著重強(qiáng)調(diào)的還是popupmenu3、popupmenu4的回調(diào)函數(shù)及調(diào)制解調(diào)函數(shù)bfsk。function popupmenu3_Callback(hObject, eventdata, handles) axes(handles.axes2);ab=1;val=get(handles.popupmenu3,'value');switch val case 1 ab=1; case 2 ab=2; case 3 ab=3; case 4 ab=4; case 5

23、ab=5; case 6 ab=6; case 7 ab=7; case 8 ab=8; case 9 ab=9; case 10 ab=10; case 11 ab=11; case 12 ab=12; case 13 ab=13; case 14 ab=14; case 15 ab=15; case 16 ab=16; end bfsk1(ab); function popupmenu4_Callback(hObject, eventdata, handles)axes(handles.axes3);ab=1;val=get(handles.popupmenu4,'value

24、9;);switch val case 1 ab=1; case 2 ab=2; case 3 ab=3; case 4 ab=4; case 5 ab=5; case 6 ab=6; case 7 ab=7; case 8 ab=8; case 9 ab=9; case 10 ab=10; case 11 ab=11; case 12 ab=12; case 13 ab=13; case 14 ab=14; case 15 ab=15; case 16 ab=16; end bfsk1(ab); function bfsk1(ba)%-產(chǎn)生原信號(hào)（二進(jìn)制單極性不歸零脈沖序列）-w=eval(

25、get(findobj('tag','edit1'),'string');f0=eval(get(findobj('tag','edit2'),'string');f1=eval(get(findobj('tag','edit3'),'string');le=length(w);t=0:0.001:(0.05*le-0.001);c1=sin(2*pi*f1*t);c0=sin(2*pi*f0*t);bit=1:le;fsk=1:le;for i=1:

26、le for j=1:50 n=(i-1)*50+j; if w(i)=1 bit(n)=1;fsk(n)=c1(n); %原信號(hào)bir和已調(diào)信號(hào)fsk else bit(n)=0;fsk(n)=c0(n); end endendfftfsk=abs(fft(fsk); %已調(diào)信號(hào)的頻譜fftfsknos=awgn(fsk,20); %加噪已調(diào)信號(hào)nosffsk=abs(fft(nos); %加噪已調(diào)信號(hào)的頻譜ffsk%-帶通濾波器-fpl=f0-15;fpu=f0+15;fsl=f0-40;fsu=f0+40;Fs=10000;wp0=2*fpl/Fs,2*fpu/Fs;ws0=2*fsl/

27、Fs,2*fsu/Fs;rp=1;rs=50;N0,wp0=ellipord(wp0,ws0,rp,rs);B0,A0=ellip(N0,rp,rs,wp0);y0=filter(B0,A0,nos); %經(jīng)帶通濾波器0后的c0(t)-y0 k0=abs(fft(y0); %經(jīng)帶通濾波器0后的c0(f)-k0 fpl1=f1-8;fpu1=f1+8;fsl1=f1-20;fsu1=f1+20;wp1=2*fpl1/Fs,2*fpu1/Fs;ws1=2*fsl1/Fs,2*fsu1/Fs;N1,wp1=ellipord(wp1,ws1,rp,rs);B1,A1,=ellip(N1,rp,rs,wp1);y1=filter(B1,A1,nos); %經(jīng)帶通濾波器1后的c1(t)-y1k1=abs(fft(y1); %經(jīng)帶通濾波器1后的c1(f)-k1 %-相乘器-length(y0);length(c0);m0=y0.*c0; %與載波相乘后c0(t)-m0m1=y1.*c1; %與載波相乘后c1(t)-m1%-低通濾波器0-Fsm=1000; %采樣頻率wpm=f0*0.6; wsm=f0*1.2;rpm=1; rsm=20;Nm,wcm=buttord(wpm,wsm,rpm,rsm,'s');Bm,Am=butter(Nm,wcm,&