基于MATLAB的FSK調制解調1(1)

1、 目錄TOC o 1-5 h z HYPERLINK l bookmark2 o Current Document 1引言3 HYPERLINK l bookmark4 o Current Document 1.1課程設計目的3 HYPERLINK l bookmark6 o Current Document 1.2課程設計要求3 HYPERLINK l bookmark8 o Current Document 1.3課程設計步驟31.4設計平臺簡介4 HYPERLINK l bookmark12 o Current Document 2設計原理5 HYPERLINK l bookmark14

2、 o Current Document FSK基本原理5 HYPERLINK l bookmark10 o Current Document FSK調制原理5 HYPERLINK l bookmark16 o Current Document FSK解調原理6 HYPERLINK l bookmark18 o Current Document 3仿真實現(xiàn)過程7 HYPERLINK l bookmark20 o Current Document FSK信號的產生7 HYPERLINK l bookmark22 o Current Document FSK的解調8 HYPERLINK l book

3、mark24 o Current Document 3.3疊加噪聲的FSK解調9 HYPERLINK l bookmark28 o Current Document 4遇到問題及解決辦法12 HYPERLINK l bookmark30 o Current Document 4.1編程中問題解答12 HYPERLINK l bookmark32 o Current Document 4.2程序解調過程中的問題12 HYPERLINK l bookmark34 o Current Document 結束語13 HYPERLINK l bookmark36 o Current Document 參

4、考文獻14基于MATLAB的FSK調制解調摘要：本課程設計利用MATLAB集成環(huán)境下的M文件，編寫程序來實現(xiàn)FSK的調制解調，并繪制出解調前后的時域和頻域波形及疊加噪聲時解調前后的時頻波形，并觀察解調前后頻譜有何變化以加深對FSK信號解調原理的理解。對信號疊加噪聲，并進行解調，繪制出解調前后信號的時頻波形，改變噪聲功率進行解調，根據(jù)運行結果和波形來分析該解調過程的正確性及信道對信號傳輸?shù)挠绊?。完成整個FSK的調制解調過程。程序開發(fā)平臺為MATLAB7.1,使用其自帶的M文件實現(xiàn)。運行平臺為Windows2000。關鍵詞：程序設計；FSK；調制解調；MATLAB7.1；M文件1引言本課程設計是利

5、用MATLAB集成環(huán)境下的M文件，編寫程序來實現(xiàn)FSK的調制解調，并繪制出解調前后的時域和頻域波形及疊加噪聲時解調前后的時頻波形，根據(jù)運行結果和波形來分析該解調過程的正確性及信道對信號傳輸?shù)挠绊憽?.1課程設計目的此次課程設計的目的是熟悉MATLAB中M文件的使用方法，編寫M文件實現(xiàn)FSK的調制和解調，繪制出FSK信號解調前后在時域和頻域中的波形，觀察調解前后頻譜的變化，再對信號進行噪聲疊加后解調同樣繪制解調前后的信號時頻波形，最后改變噪聲功率進行調解，分析噪聲對信號傳輸造成的影響，加深對FSK信號解調原理的理解。1.2課程設計要求熟悉MATLAB中M文件的使用方法，并在掌握FSK調制解調原理

6、的基礎上,編寫出FSK調制解調程序。在M文件環(huán)境下運行程序繪制出FSK信號解調前后在時域和頻域中的波形，觀察波形在解調前后的變化，對其作出解釋，同時對信號加入噪聲后解調，得到解調后的時頻波形，分析噪聲對信號傳輸造成的影響。解釋所得到的結果。1.3課程設計步驟本課程設計采用M文件編寫的方法實現(xiàn)二進制的FSK的調制與解調，然后在信號中疊加高斯白噪聲。一，調用dmode函數(shù)實現(xiàn)FSK的解調，并繪制出FSK信號調制前后在時域和頻域中的波形，兩者比較。二，調用ddemod函數(shù)解調，繪制出FSK信號解調前后在時域和頻域中的波形，兩者比較。三，調用awgn函數(shù)在新海中疊加不同信噪比的噪聲，繪制在各種噪聲下的

7、時域頻域圖。最后分析結果。1.4設計平臺簡介Matlab是美國MathWorks公司開發(fā)的用于概念設計，算法開發(fā)，建模仿真,實時實現(xiàn)的理想的集成環(huán)境。是目前最好的科學計算類軟件。作為和MathematicalMaple并列的三大數(shù)學軟件。其強項就是其強大的矩陣計算以及仿真能力。Matlab的由來就是Matrix+Laboratory=Matlab，這個軟件在國內也被稱作矩陣實驗室。Matlab提供了自己的編譯器：全面兼容C+以及Fortran兩大語言。Matlab7.1于2005.9最新發(fā)布-完整版,提供了MATLAB，SIMULINK的升級以及其他最新的75個模塊的升級。7.1版本提高了產品

8、質量,同時也提供了新的用于數(shù)據(jù)分析、大規(guī)模建模、固定點開發(fā)、編碼等。所以Matlab是工程師，科研工作者手上最好的語言，最好的工具和環(huán)境。2設計原理2.1FSK基本原理頻移鍵控是利用載波的頻率變化來傳遞數(shù)字信息。在2FSK中載波的頻率隨二進制基帶信號在fl和f2兩個頻率點間變化。故其表達式為e2fsk（t）二Acos（31t+申n）發(fā)送“1時；e2fsk（t）二Acos（32t+0n）發(fā)送“0時。一個2FSK可以看成兩個不同載頻的2ASK信號的疊加。因此2FSK信號時域表達式又可以寫成e2fsk(t)=工ang(t-nTs)cos(3lt+申n)+工ang(t-nTs)cos(32+0n)n(

9、2-1)式中：g（t）為單個矩形脈沖，脈寬為Ts；an=1時，概率為P;an=0時，概率為1-P;an是an的反碼，若an=1，則an=0；若若an=0，則an=1,于是an=1時，概率為l-p；an=0時，概率為P;n和On分別是第n個信號碼元（1或0）的初始相位。在移頻鍵控中n和On不攜帶信息，通常和令n和On為零。因此2FSK信號表達式可簡化為:e2fsk(t)二s1(t)cos31t+s2(t)cos32t(2-2)s1(t)二工ang(t-nTs)s2(t)二ang(t-nTs)(2-3)(2-4)2.2FSK調制原理在二進制頻移鍵控中，幅度恒定不變的載波信號的頻率隨著輸入碼流的變化

10、而切換（稱為高音和低音，代表二進制的1和0）。產生FSK信號最簡單的方法是根據(jù)輸入的數(shù)據(jù)比特是0還是1,在兩個獨立的振蕩器中切換。采用這種方法產生的波形在切換的時刻相位是不連續(xù)的，因此這種FSK信號稱為不連續(xù)FSK信號。）其實現(xiàn)如圖一所示：2.3FSK解調原理對于FSK信號的解調方式很多：相干解調、濾波非相干解調、正交相乘非相干解調。而FSK的非相干解調一般采用濾波非相干解調，解調原理是將2FSK信號分為上下兩路2ASK信號分別進行解調，然后判決，這里的抽樣判決是直接比較兩路信號抽樣值的大小，可以不專門設置門限。判決則應與調制規(guī)則相呼應,調制時若規(guī)定“1”符號為對應載波頻率fl,則接受時上支路

11、的樣值大小，應判為“1”反之則判為“0”。圖2-22FSK信號非相干解調原理圖3仿真實現(xiàn)過程3.1FSK信號的產生調用domde函數(shù)實現(xiàn)FSK的調制：y二dmod(x,Fc,Fd,Fs,fsk,M,df);(3T)其中y為已調信號；x為調制信號；Fc=10為載波頻率；Fs=40系統(tǒng)采樣頻率；Fd=l碼元速率。此程序輸入的是二進制的FSK信號且為數(shù)字信號，所謂為連續(xù)的二進制圖形如圖三所示：O765432100.0.0.0.0.0.0.-gnl_-dE10.90.8二進制隨機序列146im圖3-1輸入的二進制隨即序列對該信號進行頻譜分析在M文件下得到了調制后的頻域和時域圖設定好各個的參數(shù)調制后的信

12、號505n-O.apnl_-duu圖3-2調制后的信號圖3-3調制后的頻譜圖通過調制，從所得的圖形可知輸入的二進制隨機信號沒受到噪聲干擾的時候是很規(guī)則的信號調制后的信號3.2FSK的解調對調制好的FSK信號進行非相干解調。輸入的FSK中頻信號分別經過中心頻率為fH、fL的帶通濾波器，然后分別經過包絡檢波，包絡檢波的輸出在t=kTb。時抽樣（其中k為整數(shù)），并且將這些值進行比較。根據(jù)包絡檢波器輸出的大小,比較器判決數(shù)據(jù)比特是1還是0。:553E551?5?iiE55=|E02468101214Time非相干解調后的信號aAAaAmAaeeX1-.30ei原輸入二進制隨機序列02468101214

13、Elp_l_-duuTime 圖3-4解調前后的頻譜圖比較從所得到的結果來看，經過解調后，在沒有噪聲的干擾下，解調出來的信號在頻譜上和原信號十分吻合，基本沒有誤差。3.3疊加噪聲的FSK解調在信號中如果存在了不需要的電信號，我們稱之為噪聲。而在通信系統(tǒng)中，我們一般需要處理的是高斯白噪聲。一般在物理上把它翻譯成白噪聲(whitenoise)。白噪聲是指功率譜密度在整個頻域內均勻分布的噪聲。所有頻率具有相同能量的隨機噪聲稱為白噪聲。從我們耳朵的頻率響應聽起來它是非常明亮的“咝”聲。白噪聲或白雜訊，是一種功率頻譜密度為常數(shù)的隨機信號或隨機過程。換句話說，此信號在各個頻段上的功率是一樣的，由于白光是由

14、各種頻率(顏色)的單色光混合而成，因而此信號的這種具有平坦功率譜的性質被稱作是“白色的”，此信號也因此被稱作白噪聲。相對的，其他不具有這一性質的噪聲信號被稱為有色噪聲。理想的白噪聲具有無限帶寬，因而其能量是無限大，這在現(xiàn)實世界是不可能存在的。實際上，我們常常將有限帶寬的平整訊號視為白噪音，因為這讓我們在數(shù)學分析上更加方便。然而，白噪聲在數(shù)學處理上比較方便，因此它是系統(tǒng)分析的有力工具。一般，只要一個噪聲過程所具有的頻譜寬度遠遠大于它所作用系統(tǒng)的帶寬，并且在該帶寬中其頻譜密度基本上可以作為常數(shù)來考慮，就可以把它作為白噪聲來處理。我們要分析的是在加了高斯白噪聲后的FSK信號。調用awgn函數(shù)來實現(xiàn)給

15、FSK信號的噪聲疊加：yl二awgn(y,20);(3-2)yl為加了噪聲后的信號，y為需要疊加噪聲的信號。調用前面已經產生的調制信號，我們設定噪聲比為20，則可以得到加入噪聲后的已調信號和原二進制隨機序列的比較。如圖6所示調制后的信號1.51111二進制隨機序列46810Time加入高斯白噪聲后的已調信號5cdsnl一-dEoapnl-dE505o.-O.Bpnl=duuLnu51214468101214Time圖3-5原始信號調制后信號和加入高斯白噪聲后的已調信號比較從圖中可知，加入噪聲之后信號的波動明顯增大，幅度超出了“1”“-1”的范圍，信號的失真情況很明顯，但對信號的波動沒有很大的改

16、變。由此說明，疊加噪聲之后信號傳輸效果會受到影響，但正確性不會改變。3QC20000.20.40.6o.e1-2141.S1e2圖3-6調制后的頻譜圖300200100Ja.21110.SUP11.21.41.61.823-7加入高斯白噪聲解調后的頻譜圖（snr=10）3-8加入高斯白噪聲解調后的頻譜圖（snr=3）加入噪聲后，信號的變化很明顯。改變噪聲的大小：當信噪比為10時，失真不是非常明顯，與原波形基本上沒大的變化。當信噪比為3時，噪聲大了之后對信號的影響非常大。說明FSK信號具有一定的抗噪聲能力。pnl_Qu10.80.60.40.2原輸入二進制隨機序列281012144非相干WiSf

17、e的信號_eie0.80.63-9加入噪聲后非相干解調后,加J入噪聲后,對信號進行非E相干解調,與原信號:殳有變化，任然能正確的解調出原垃始信號，這是因為雖然加并且我們取的數(shù)據(jù)組數(shù)有限，不能很好的反應其真實的抗噪能力。通過上述的操作之后，信道中噪聲對信號的傳輸有一定的影響，噪聲越大影口入了噪聲，但數(shù)字信號的抗噪,性能較模擬信號而0響越明顯，掩蓋了原始信號，所以我們在平常信號傳輸中，盡量減少信道中噪聲的干擾，保證信號的傳輸效率。4遇到問題及解決辦法4.1編程中問題解答在此次課程實際過程中，碰到了很多的問題，而解決這些問題使我學得了很多的新知識，同時也重拾以前遺忘的知識，加深了對現(xiàn)有知識的鞏固。在

18、同學們和老師的幫助下，將這些問題一一解決。現(xiàn)將問題做如下總結。首先是在編寫程序時，雖然在網上找到了相關程序，但這些程序是網上程序員編寫的，其中有許多的函數(shù)是我們的MATLAB雖然有，但是我們所謂涉及過的，因而無法讀懂這些函數(shù)對其作出解釋。而我要做的是將這些函數(shù)變?yōu)樽约含F(xiàn)有的程序在MATLAB的函數(shù)庫中現(xiàn)有的我們接觸過的能自己組織語言描述的。例如：yl二awgn(y,SNR-10*logl0(0.5)-10*loglO(N),measured,dB);此函數(shù)在以前的課程中并未看到過，在問過老師之后，才真正明白其含義：次句程序的意思就和我們的yl=awgn(y,20);一樣，在信號中添加高斯白噪聲

19、，對于這一系列的函數(shù)和程序用語問題，不懂的經過老師和同學的解答，都得到了很好的解釋。也是我懂得在平時學習中不能只局限于課本知識，對于課后的各類文章也要去多多留意。4.2程序解調過程中的問題對于程序中做出FSK解調后的頻譜圖的程序我們可以單獨作為一個M文件在MATLAB中運行，這樣就減少了程序在各個M文件中存在的問題，同樣也就較少了程序過于繁瑣的問題。對于矩陣n二O:l:length(y)T;w二0:1:1000-1*pi/500;)在沒有用1000-1之前，程序一直運行有錯誤。在請教了我們班的同學之后，按照矩陣乘法的規(guī)則第一個矩陣的列數(shù)要和第二個矩陣的行數(shù)相等，否則不能相乘。所以在定義1000個