基于matlab的2fsk的課程設計

1、<p><b>　　課程設計</b></p><p><b>　　摘 要</b></p><p>　　本文主要是利用MATLAB7.0來實現(xiàn)2FSK數(shù)字調制系統(tǒng)解調器的設計。該設計模塊包含信源調制、發(fā)送濾波器模塊、信道、接受濾波器模塊、解調以及信宿，并未各個模塊進行相應的參數(shù)設置。在此基礎上熟悉MATLAB的功能及操作，最后通過觀

2、察仿真圖形進行波形分析及系統(tǒng)的性能評價。</p><p>　　關鍵詞：2FSK MATLAB 調制解調</p><p><b>　　目 錄</b></p><p><b>　　引 言1</b></p><p>　　1. 背景知識介紹2</p><p>　　

3、1.2MATLBE介紹4</p><p>　　2. 2FSK的調制與解調原理6</p><p>　　2.1 2FSK的產(chǎn)生6</p><p>　　2.2 2FSK濾波器的解調和抗噪聲性能9</p><p>　　3仿真系統(tǒng)模型設計12</p><p>　　3.1仿真思路12</p><p&

4、gt;　　3.2仿真結果12</p><p>　　3.3 Simulink仿真模型圖15</p><p>　　3.4結果分析18</p><p>　　3.4.1Matlab仿真結果分析18</p><p>　　3.4.2 simulink仿真結果分析19</p><p>　　3.5心得體會19</p&

5、gt;<p><b>　　4參考文獻21</b></p><p>　　附錄1 MATLAB仿真程序22</p><p><b>　　引 言</b></p><p>　　2FSK信號的產(chǎn)生方法主要有兩種：一種是調頻法，一種是開關法。這兩種方法產(chǎn)生的2FSK信號的波形基本相同，只有一點差異，即由調頻器產(chǎn)

6、生的2FSK信號在相鄰碼元之間的相位是連續(xù)的，而開關法產(chǎn)生的2FSK信號則分別由兩個獨立的頻率源產(chǎn)生兩個不同頻率的信號，故相鄰碼元之間的相位不一定是連續(xù)的。本設計采用后者--開關法。2FSK信號的接收也分為相干和非相干接收兩種，非相干接收方法不止一種，它們都不利用信號的相位信息。故本設計采用相干解調法。</p><p><b>　　1.背景知識介紹</b></p><p&

7、gt;<b>　　1.1通信介紹</b></p><p>　　通信就是克服距離上的障礙，從一地向另一地傳遞和交換消息。消息是信息源所產(chǎn)生的，是信息的物理表現(xiàn)，例如，語音、文字、數(shù)據(jù)、圖形和圖像等都是消息。消息有模擬消息（如語音、圖像等）以及數(shù)字消息（如數(shù)據(jù)、文字等）之分。所有消息必須在轉換成電信號（通常簡稱為信號）后才能在通信系統(tǒng)中傳輸。所以，信號是傳輸消息的手段，信號是消息的物質載體。&l

8、t;/p><p>　　相應的信號可分為模擬信號和數(shù)字信號，模擬信號的自變量可以是連續(xù)的或離散的，但幅度是連續(xù)的，如電話機、電視攝像機輸出的信號就是模擬信號。數(shù)字信號的自變量可以是連續(xù)的或離散的，但幅度是離散的，如電船傳機、計算機等各種數(shù)字終端設備輸出的信號就是數(shù)字信號。</p><p>　　通信的目的是傳遞消息，但對受信者有用的是消息中包含的有效內容，也即信息。消息是具體的、表面的，而信息是抽

9、象的、本質的，且消息中包含的信息的多少可以用信息量來度量。</p><p>　　通信技術，特別是數(shù)字通信技術近年來發(fā)展非常迅速，它的應用越來越廣泛。通信從本質上來講就是實現(xiàn)信息傳遞功能的一門科學技術，它要將大量有用的信息無失真，高效率地進行傳輸，同時還要在傳輸過程中將無用信息和有害信息抑制掉。當今的通信不僅要有效地傳遞信息，而且還有儲存、處理、采集及顯示等功能，通信已成為信息科學技術的一個重要組成部分。</

10、p><p>　　通信系統(tǒng)就是傳遞信息所需要的一切技術設備和傳輸媒質的總和，包括信息源、發(fā)送設備、信道、接收設備和信宿(受信者) ，它的一般模型如圖1所示。</p><p><b>　　↑</b></p><p>　　圖1 通信系統(tǒng)一般模型</p><p>　　通信系統(tǒng)可分為數(shù)字通信系統(tǒng)和模擬通信系統(tǒng)。數(shù)字通信系統(tǒng)是利用數(shù)字信

11、號來傳遞消息的通信系統(tǒng)，其模型如圖2所示，</p><p><b>　　↑</b></p><p>　　圖2 數(shù)字通信系統(tǒng)模型</p><p>　　模擬通信系統(tǒng)是利用模擬信號來傳遞消息的通信系統(tǒng)，其模型如圖3所示。</p><p>　　圖3 模擬通信系統(tǒng)模型</p><p>　　數(shù)字通信系統(tǒng)較模擬

12、通信系統(tǒng)而言，具有抗干擾能力強、便于加密、易于實現(xiàn)集成化、便于與計算機連接等優(yōu)點。因而，數(shù)字通信更能適應對通信技術的越來越高的要求。近二十年來，數(shù)字通信發(fā)展十分迅速，在整個通信領域中所占比重日益增長，在大多數(shù)通信系統(tǒng)中已代替模擬通信，成為當代通信系統(tǒng)的主流。</p><p>　　在當前飛速發(fā)展的信息時代，隨著數(shù)字通信技術計算機技術的發(fā)展，以及通信網(wǎng)絡與計算機網(wǎng)絡的相互融合，信息技術已成為21世紀社會國際化的強大動

13、力。信息作為一種資源，只有通過廣泛的傳播與交流，才能產(chǎn)生利用價值，促進社會成員之間的合作，推動社會生產(chǎn)力的發(fā)展，創(chuàng)造巨大的經(jīng)濟效益。而信息的傳播和交流，是依靠各種通信方式和技術來實現(xiàn)的。學習和掌握現(xiàn)代通信理論和技術是信息社會每一位成員，尤其是未來通信工作著的迫切要求。</p><p>　　1.2MATLBE介紹</p><p>　　在70年代中期，Cleve Moler博士和其同事在美國國

14、家科學基金的資助下開發(fā)了調用EISPACK和LINPACK的FORTRAN子程序庫.EISPACK是特征值求解的FOETRAN程序庫，LINPACK是解線性方程的程序庫.在當時,這兩個程序庫代表矩陣運算的最高水平。到70年代后期,身為美國New Mexico大學計算機系系主任的Cleve Moler，在給學生講授線性代數(shù)課程時。想教學生使用EISPACK和LINPACK程序庫，但他發(fā)現(xiàn)學生用FORTRAN編寫接口程序很費時間。于是他開始

15、自己動手,利用業(yè)余時間為學生編寫EISPACK和LINPACK的接口程序，Cleve Moler給這個接口程序取名為MATLAB，該名為矩陣(matrix)和實驗室(labotatory)兩個英文單詞的前三個字母的組合。在以后的數(shù)年里，MATLAB在多所大學里作為教學輔助軟件使用。并作為面向大眾的免費軟件廣為流傳。1983年春天，Cleve Moler到Standford大學講學，MATLAB深深地吸引了工程師John Little.J

16、ohn Little敏銳地覺察到MATLAB在工程領域的廣闊前</p><p>　　在MATLAB進入市場前，國際上的許多軟件包都是直接以FORTRANC語言等編程語言開發(fā)的。這種軟件的缺點是使用面窄，接口簡陋，程序結構不開放以及沒有標準的基庫，很難適應各學科的最新發(fā)展，因而很難推廣。MATLAB的出現(xiàn)，為各國科學家開發(fā)學科軟件提供了新的基礎。在MATLAB問世不久的80年代中期，原先控制領域里的一些軟件包紛紛被

17、淘汰或在MATLAB上重建。時至今日，經(jīng)過MathWorks公司的不斷完善，MATLAB已經(jīng)發(fā)展成為適合多學科，多種工作平臺的功能強大大大型軟件。在國外，MATLAB已經(jīng)經(jīng)受了多年考驗。在歐美等高校，MATLAB已經(jīng)成為線性代數(shù)，自動控制理論，數(shù)理統(tǒng)計，數(shù)字信號處理，時間序列分析，動態(tài)系統(tǒng)仿真等高級課程的基本教學工具；成為攻讀學位的大學生，碩士生，博士生必須掌握的基本技能。在設計研究單位和工業(yè)部門，MATLAB被廣泛用于科學研究和解決各

18、種具體問題。在國內，特別是工程界，MATLAB一定會盛行起來?？梢哉f，無論你從事工程方面的哪個學科，都能在MATLAB里找到合適的功能。</p><p>　　2. 2FSK的調制與解調原理</p><p>　　二進制頻率調制是用二進制數(shù)字信號控制正弦波的頻率隨二進制數(shù)字信號的變化而變化。由于二進制數(shù)字信息只有兩個不同的符號，所以調制后的已調信號有兩個不同的頻率f1和f2，f1對應數(shù)字信息“

19、1”，f2對應數(shù)字信息“0”。二進制數(shù)字信息及已調載波如圖4所示。</p><p>　　圖4 2FSK信號 </p><p>　　2.1 2FSK的產(chǎn)生</p><p>　　在2FSK信號中，當載波頻率發(fā)生變化時，載波的相位一般來說是不連續(xù)的，這種信號稱為不連續(xù)2FSK信號。相位不連續(xù)的2FSK通常用頻率選擇法產(chǎn)生，如圖5所示：</p><p&g

20、t;　　振蕩器1 c 選通開關 e</p><p>　　a基帶信號 反相器 b 相加器 g</p><p>　　振蕩器2 d 選通開關 f</p><p>　　圖5 2FSK調制器</p><p>　　兩個獨立的振蕩器作為兩個頻率發(fā)生器，他們受控于輸

21、入的二進制信號。二進制信號通過兩個選通開關電路，控制其中的一個載波通過。調制器各點波形如圖6所示：</p><p>　　圖6 2FSK調制器各點波形</p><p>　　由圖6可知，波形g是波形e和f的疊加。所以，二進制頻率調制信號2FSK可以看成是兩個載波頻率分別為f1和f2的2ASK信號的和。由于“1”、“0”統(tǒng)計獨立，因此，2FSK信號功率譜密度等于這兩個2ASK信號功率譜密度之和，

22、即</p><p><b>　?。?-1）</b></p><p>　　2FSK信號的功率譜如圖7所示：</p><p>　　圖7、2FSK信號的功率譜</p><p>　　由圖7看出，2FSK信號的功率譜既有連續(xù)譜又有離散譜，離散譜位于兩個載波頻率f1和f2處，連續(xù)譜分布在f1和f2附近，若取功率譜第一個零點以內的成分

23、計算帶寬，顯然2FSK信號的帶寬為</p><p><b>　　（2-2）</b></p><p>　　為了節(jié)約頻帶，同時也能區(qū)分f1和f2，通常取，因此2FSK信號的帶寬為</p><p><b>　?。?-3）</b></p><p>　　當時，圖2-4中2FSK的功率譜由雙峰變成單峰，此時帶寬

24、為</p><p><b>　?。?-4）</b></p><p>　　對于功率譜是單峰的2FSK信號，可采用動態(tài)濾波器來解調。此處介紹功率譜為雙峰的2FSK信號的解調。</p><p>　　2.2 2FSK濾波器的解調和抗噪聲性能</p><p>　　2FSK信號的解調也有相干解調和非相干解調（包絡解調）兩種。由于2F

25、SK信號可看做是兩個2ASK信號之和，所以2FSK解調器由兩個并聯(lián)的2ASK解調器組成。如圖8為2FSK的相干解調和非相干解調： </p><p>　　圖8 2FSK信號調解器</p><p>　　相干2FSK抗噪聲性能的分析方法和相干2ASK很相似。現(xiàn)將收到的2FSK信號表示為</p><p><b>　?。?-5）</b></p&g

26、t;<p>　　當發(fā)送數(shù)字信息為“1”時，2FSK信號的載波頻率為f1，信號能通過上支路的帶通濾波器。上支路帶通濾波器的輸出是信號和窄帶噪聲(t)的疊加(噪聲中的下標1表示上支路窄帶高斯噪聲)，即</p><p><b>　?。?-6）</b></p><p>　　此信號與同步載波cos2πf1t相乘，再經(jīng)低通濾波器濾除其中的高頻成分，送給取樣判決器的信

27、號為</p><p><b>　?。?-7）</b></p><p>　　上式中未計入系數(shù)1/2。與此同時，頻率為f1的2FSK信號不能通過下支路中的帶通濾波器，因為下支路中的帶通濾波器的中心頻率為f2，所以下支路帶通濾波器的輸出只有窄帶高斯噪聲，即</p><p><b>　　（2-8）</b></p>&

28、lt;p>　　此噪聲與同步載波cos2πf2t相乘，再經(jīng)低通濾波器濾波后輸出為</p><p><b>　?。?-9）</b></p><p>　　上式中未計入系數(shù)1/2。定義</p><p><b>　　（2-10）</b></p><p>　　取樣判決器對x(t)取樣，取樣值為</p

29、><p><b>　?。?-11）</b></p><p>　　其中，nI1、 nI2都是均值為0、方差為的高斯隨機變量，所以x是均值為a、方差為的高斯隨機變量，x的概率密度函數(shù)為</p><p><b>　?。?-12）</b></p><p>　　概率密度曲線如圖9所示：</p>&l

30、t;p>　　圖9、 判決值的函數(shù)示意圖</p><p>　　判決器對x進行判決，當x>0時，判發(fā)送信息為“1”，此判決是正確的； 當x<0時，判決發(fā)送信息為“0”，顯然此判決是錯誤的。由此可見，x<0的概率就是發(fā)“1”錯判成“0”的概率，即</p><p><b>　?。?-13）</b></p><p>　　當發(fā)送數(shù)字

31、信號“0”時，下支路有信號，上支路沒有信號。用與上面分析完全相同的方法，可得到發(fā)“0” 碼時錯判成“1”碼的概率P(1/0)，容易發(fā)現(xiàn)，此概率與上式表示的P(0/1)相同，所以解調器的平均誤碼率為</p><p>　　Pe=P(1)P(0/1)+P(0)P(1/0)=P(0/1)［P(1)+P(0)］=P(0/1) (2-14)</p><p><b>　　所以<

32、/b></p><p><b>　　(2-15)</b></p><p><b>　　式中</b></p><p>　　注意，式中無需“1”、“0”等概這一條件。</p><p><b>　　3仿真系統(tǒng)模型設計</b></p><p><b&

33、gt;　　3.1仿真思路</b></p><p>　　先確定采樣頻率fs和兩個載波頻率的值f1，f2，產(chǎn)生2FSK信號，然后進行調制與解調。</p><p>　?。?）寫出輸入已知信號的表達式S(t)。由于S(t)中有反碼的存在，則需要將信號先反轉后在從原信號和反轉信號中進行抽樣。寫出已調信號的表達式S(t)。</p><p>　?。?）在2FSK的解調

34、過程中，信號首先通過帶通濾波器，設置帶通濾波器的參數(shù)，后用一維數(shù)字濾波函數(shù)filter對信號S(t)的數(shù)據(jù)進行濾波處理。輸出經(jīng)過帶通濾波器后的信號波形。由于已調信號中有兩個不同的載波（ω1,ω2）,則經(jīng)過兩個不同頻率的帶通濾波器后輸出兩個不同的信號波形H1,H2。</p><p>　?。?）經(jīng)過帶通濾波器后的2FSK信號再經(jīng)過相乘器（cosω1，cosω2），兩序列相乘的MATLAB表達式y(tǒng)=x1.*x2→SW=

35、Hn.*Hn，輸出得到相乘后的兩個不同的2FSK波形h1,h2。</p><p>　?。?）經(jīng)過相乘器輸出的波形再通過低通濾波器，設置低通濾波器的參數(shù)，用一維數(shù)字濾波韓式filter對信號的數(shù)據(jù)進行新的一輪的濾波處理。輸出經(jīng)過低通濾波器后的兩個波形（sw1,sw2）。</p><p>　?。?）將信號sw1和sw2同時經(jīng)過抽樣判決器，分別輸出st1,st2。其抽樣判決器輸出的波形為最后的輸

36、出波形st。對抽樣判決器經(jīng)定義一個時間變量長度i，當st1(i)>=st2(i)5時，則st=0，否則st=st2(i).其中st=st1+st2。</p><p>　?。?）使用MATLAB編程（m文件）完成系統(tǒng)的仿真。</p><p><b>　　3.2仿真結果</b></p><p><b>　　仿真程序見附錄1。<

37、/b></p><p><b>　　仿真圖像如下</b></p><p>　　圖3-1 MATLAB仿真圖</p><p>　　圖3-2 MATLAB仿真圖</p><p>　　圖3-3 MATLAB仿真圖</p><p>　　圖3-4 MATLAB仿真圖</p><p&

38、gt;　　圖3-5 MATLAB仿真圖</p><p>　　3.3 Simulink仿真模型圖</p><p>　　2FSK信號的Simulink仿真模型圖如下：</p><p>　　圖3-6 Simulink仿真模型圖</p><p>　　其中sin wave和sin wave1是兩個頻率分別為f1和f2的載波，Pulse Generato

39、r模塊是信號源，NOT實現(xiàn)反相，再經(jīng)過相乘器和相加器生成2FSK信號，然后接帶通濾波器與低通濾波器完成調制與解調，又接入Error Rate calculation實現(xiàn)誤碼率的計算。</p><p><b>　　參數(shù)設置如下</b></p><p>　　表3-1 Simulink仿真各模型的參數(shù)設置</p><p>　　經(jīng)過以上參數(shù)的設定后就可

40、以進行系統(tǒng)的仿真，其各點的時間波形如下：</p><p>　　圖3-7 Simulink仿真各點的波形</p><p>　　由仿真系統(tǒng)中的誤碼率計算可知：此系統(tǒng)的誤碼率為0.1818</p><p>　　提高噪聲影響，將Gaussian Noise Generator中的Variance“1”改為“10”,得到波形如下圖所示：</p><p>

41、;　　圖3-8 改變參數(shù)后各點的波形</p><p>　　由仿真系統(tǒng)中的誤碼率計算可知：此系統(tǒng)的誤碼率為0.3636</p><p>　　再次提高噪聲影響，將Gaussian Noise Generator中的Variance“1”改為“100”得到波形如下圖所示</p><p>　　圖3-9 修改參數(shù)后得到的各點的波形</p><p>　　

42、由仿真系統(tǒng)中的誤碼率計算可知：此系統(tǒng)的誤碼率為0.4545</p><p><b>　　3.4結果分析</b></p><p>　　3.4.1Matlab仿真結果分析</p><p>　　2FSK信號的調制解調原理是通過帶通濾波器將2FSK信號分解為上下兩路2FSK信號后分別解調，然后進行抽樣判決輸出信號。本實驗對信號2FSK采用相干解調進行解

43、調。對于2FSK系統(tǒng)的抗噪聲性能，本實驗采用同步檢測法。設“1”符號對應載波頻率f1，“0”符號對應載波頻率f2。在原理圖中采用兩個帶通濾波器來區(qū)分中心頻率分別為f1和f2的信號。中心頻率為f1的帶通濾波器之允許中心頻率為f1的信號頻譜成分通過，濾除中心頻率為f2的信號頻譜成分。</p><p>　　接收端上下支路兩個帶通濾波器的輸出波形中H1,H2。在H1,H2波形中在分別含有噪聲n1,n2，其分別為高斯白噪聲

44、ni經(jīng)過上下兩個帶通濾波器的輸出噪聲——窄帶高斯噪聲，其均值同為0，方差同為(σn)2,只是中心頻率不同而已。</p><p>　　其抽樣判決是直接比較兩路信號抽樣值的大小，可以不專門設置門限。判決規(guī)制應與調制規(guī)制相呼應，調制時若規(guī)定“1”符號對應載波頻率f1，則接收時上支路的抽樣較大，應判為“1”，反之則判為“0”。</p><p>　　在（0，Ts）時間內發(fā)送“1”符號（對應ω1），則

45、上下支路兩個帶通濾波器輸出波形H1,H2。H1,H2分別經(jīng)過相干解調（相乘—低通）后，送入抽樣判決器進行判決。比較的兩路輸入波形分別為上支路st1=a+n1,下支路st2=n2，其中a為信號成分；n1和n2均為低通型高斯噪聲，其均值為零，方差為。當st1的抽樣值st1(i)小于st2的抽樣值st2(i),判決器輸出“0”符號，造成將“1”判為“0”的錯誤。</p><p>　　3.4.2 simulink仿真結果

46、分析</p><p>　　本次課程設計實現(xiàn)了2FSK的調制與解調過程。通過誤碼率為0的分析，可能是系統(tǒng)自身的原因造成的或者系統(tǒng)沒有多大的誤差，但根據(jù)示波器的波形看出初始的時候有一些失真原因可能是初始信號不穩(wěn)定造成的。通過改變信噪比，系統(tǒng)的誤碼率升高，噪聲干擾越大，相應的波形失真越嚴重，誤碼率越高。本實驗用一個邏輯比較器替代抽樣判決模塊，簡化了模塊。</p><p><b>　　3

47、.5心得體會</b></p><p>　　本次課程設計主要涉及到了通信原理和MATLB的相關知識與運用，主要有基帶信號的調制原理及方法、低通和帶通濾波器等等，加深了對上述相關知識的了解，使自己更深刻理解了調制與解調的原理和實現(xiàn)方法，以及基本掌握了MATLAB的基本應用。因為是以所學理論為基礎，所以在課程設計的過程中，我又重溫2FSK的調制與解調等知識，更加熟悉了MATLB里面的Simulink工具箱，

48、學會了獨立建立模型，分析調制與解調結果，和加入噪聲之后的情況，通過自己不斷的調試，更好的理解加入噪聲對信道的影響。</p><p>　　這次課程設計對我的自身能力有了進一步了解。第一點，這進一步端正了我的學習態(tài)度，學會了實事求是、嚴謹?shù)淖黠L，提高了動手能力。也要對自己嚴格要求，不能夠一知半解，要力求明明白白。浮躁的性格對于搞設計來說或者學習是致命的，一定要靜下心來，踏實的做事。第二點，我覺得動手之前，頭腦里必須清

49、楚應該怎么做，這一點是很重要的，所謂三思而后行。第三點，在這兩周的時間里，深深感覺到我對所學的知識的了解還是很局限、很膚淺的，在以后的時間里還要繼續(xù)加深對這門課程的理解，加強知識儲備。</p><p>　　在這次課程設計中我們遇到了許多的困難，由于粗心大意出了一些簡單的錯誤，浪費了許多時間去改正。還好在同學和老師的幫助下，給我指出了錯誤的原因以及改正的方法，我們組才順利的完成了本次課程設計。通過這次課程設計，我學

50、到了很多書本上沒有的知識。鍛煉了我們獨立思考問題、分析問題、解決問題的能力。而且本次設計有自己和本組成員共同完成。加強了和別人溝通的能力以及團隊精神，對我們走向社會是個很好的鍛煉。這個課程設計完成倉促，在編程過程中，我發(fā)現(xiàn)自己的程序還有很多地方可以完善，其中若有不足之前，請老師指出，我將及時改正。</p><p><b>　　4參考文獻</b></p><p>　　[

51、1]徐明遠，邵玉斌編著《MATLAB仿真在通信與電子工程中的應用》，西安電子科技大學出版社，2005</p><p>　　[2]孫屹，吳磊編著《Simulink通信仿真開發(fā)手冊》,國防工業(yè)出版社,2003</p><p>　　[3]樊昌信，曹麗娜編著《通信原理》，國防工業(yè)出版社</p><p>　　[4]張德豐等編著《MATLAB數(shù)值計算方法》，機械工業(yè)出版社<

52、;/p><p>　　附錄1 MATLAB仿真程序</p><p>　　fs=2000; %抽樣頻率</p><p><b>　　dt=1/fs;</b></p><p>　　f1=20; %載波1頻率</p><p>　　f2=50; %載波2頻率</p><p>　　a=ro

53、und(rand(1,10)); %隨機信號</p><p><b>　　g1=a;</b></p><p>　　g2=~a; %信號反轉，和g1反向</p><p>　　g11=(ones(1,2000))'*g1; %抽樣</p><p>　　g1a=g11(:)';</p><p

54、>　　g21=(ones(1,2000))'*g2;</p><p>　　g2a=g21(:)';</p><p>　　t=0:dt:10-dt;</p><p>　　t1=length(t);</p><p>　　fsk1=g1a.*cos(2*pi*f1.*t);</p><p>　　fsk2

55、=g2a.*cos(2*pi*f2.*t);</p><p>　　fsk=fsk1+fsk2; %產(chǎn)生的信號</p><p>　　no=0.01*randn(1,t1); %噪聲</p><p>　　sn=fsk+no;</p><p>　　subplot(311);</p><p>　　plot(t,no); %噪聲

56、波形</p><p>　　title('噪聲波形')</p><p>　　ylabel('幅度')</p><p>　　subplot(312);</p><p>　　plot(t,fsk);</p><p>　　title('產(chǎn)生的波形')</p><

57、;p>　　ylabel('幅度')</p><p>　　subplot(313);</p><p>　　plot(t,sn);</p><p>　　title('將要通過濾波器的波形')</p><p>　　ylabel('幅度的大小')</p><p>　　xla

58、bel('t')</p><p>　　figure(2) %FSK解調</p><p>　　b1=fir1(101,[1/800 20/800]);</p><p>　　b2=fir1(101,[40/800 60/800]); %設置帶通參數(shù)</p><p>　　H1=filter(b1,1,sn);</p>

59、<p>　　H2=filter(b2,1,sn); %經(jīng)過帶通濾波器后的信號</p><p>　　subplot(211);</p><p>　　plot(t,H1);</p><p>　　title('經(jīng)過帶通濾波器f1后的波形')</p><p>　　ylabel('幅度')</p>

60、<p>　　subplot(212)</p><p>　　plot(t,H2);</p><p>　　title('經(jīng)過帶通濾波器f2后的波形')</p><p>　　ylabel('幅度')</p><p>　　xlabel('t')</p><p>　　s

61、w1=H1.*H1;</p><p>　　sw2=H2.*H2; %經(jīng)過相乘器</p><p><b>　　figure(3)</b></p><p>　　subplot(211);</p><p>　　plot(t,sw1);</p><p>　　title('經(jīng)過相乘器h1后的波形&#

62、39;)</p><p>　　ylabel('幅度')</p><p>　　subplot(212);</p><p>　　plot(t,sw2);</p><p>　　title('經(jīng)過相乘器h2后的波形')</p><p>　　ylabel('幅度')</p&g

63、t;<p>　　xlabel('t')</p><p>　　bn=fir1(101,[2/800 10/800]); %經(jīng)過低通濾波器</p><p><b>　　figure(4)</b></p><p>　　st1=filter(bn,1,sw1);</p><p>　　st2=filt

64、er(bn,1,sw2);</p><p>　　subplot(211);</p><p>　　plot(t,st1);</p><p>　　title('經(jīng)過低通濾波器sw1后的波形')</p><p>　　ylabel('幅度')</p><p>　　subplot(212);<

65、;/p><p>　　plot(t,st2);</p><p>　　title('經(jīng)過低通濾波器sw2后的波形')</p><p>　　ylabel('幅度')</p><p>　　xlabel('t') %判決</p><p>　　for i=1:length(t)</

66、p><p>　　if(st1(i)>=st2(i))</p><p><b>　　st1(i)=0;</b></p><p>　　else st1(i)=st2(i);</p><p><b>　　end</b></p><p><b>　　end</b&g

67、t;</p><p><b>　　figure(5)</b></p><p>　　st=st1+st2;</p><p>　　subplot(211);</p><p>　　plot(t,st);</p><p>　　title('經(jīng)過抽樣判決器后的波形')</p>

68、<p>　　ylabel('幅度')</p><p>　　subplot(212);</p><p>　　plot(t,sn);</p><p>　　title('原始的波形')</p><p>　　ylabel('幅度')</p><p>　　xlabel(&