FM信號的仿真分析

1、 武漢理工大學matlab課程設(shè)計 Matlab課程設(shè)計任務(wù)書 學生姓名： 王曉丹 專業(yè)班級： 通信1103班 指導教師： 許建霞 工作單位： 信息工程學院 題 目: FM信號的仿真分析 初始條件：理論方面電路分析基礎(chǔ)、模擬電子、高頻電子線路、通信原理、數(shù)字信號處理等專業(yè)知識設(shè)備方面微型計算機軟件方面matlab仿真平臺場地方面鑒主15樓通信實驗室1要求完成的主要任務(wù): 1. 理論設(shè)計與分析，包括理論的推導和計算過程。2. 制作相應(yīng)的matlab編制仿真程序，并與理論結(jié)果對比，最后得出結(jié)論。3. 完成課程設(shè)計報告，包括必要的公式、文字/圖表說明和仿真結(jié)果。時間安排：本課程設(shè)計自任務(wù)書發(fā)布之日開

2、始，分散進行，答辯時間統(tǒng)一。指導教師簽名： 年 月 日系主任（或責任教師）簽名： 年 月 日1目 錄摘要1Abstract21緒 論31.1課題背景31.2課程設(shè)計目的與意義42 FM調(diào)制系統(tǒng)設(shè)計原理52.1調(diào)頻信號的定義及產(chǎn)生52.2 噪聲6 2.3調(diào)制信號仿真圖.73 FM解調(diào)信號模型113.1解調(diào)模型的建立113.2解調(diào)過程分析114 高斯白噪聲信道特性144.1 高斯白噪聲信道特性理論144.2 輸入輸出信噪比關(guān)系曲線145 調(diào)頻系統(tǒng)的抗噪聲性能分析187 心得體會208 參考文獻21附錄 MATLAB源程序22摘要 調(diào)制在通信系統(tǒng)中有十分重要的作用。通過調(diào)制，不僅可以進行頻譜搬移，把

3、調(diào)制信號的頻譜搬移到所希望的位置上，從而將調(diào)制信號轉(zhuǎn)換成適合于傳播的已調(diào)信號，而且它對系統(tǒng)的傳輸有效性和傳輸?shù)目煽啃杂兄艽蟮挠绊懀{(diào)制方式往往決定了一個通信系統(tǒng)的性能。FM信號的調(diào)制屬于頻譜的非線性搬移，它的解調(diào)也有相干和非相干解調(diào)兩種方式。本課程設(shè)計使用的仿真軟件為MATLAB 7.0，利用MATLAB 集成環(huán)境下的M文件，編寫程序來實現(xiàn)FM信號的仿真分析，并分別繪制出基帶信號、載波信號、已調(diào)信號的時域波形和頻域波形；再進一步分別繪制出對已調(diào)信號疊加噪聲后信號、同步解調(diào)前信號和解調(diào)后基帶信號的時域波形；最后繪出 FM 基帶信號調(diào)制和解調(diào)系統(tǒng)后的輸入輸出信噪比的關(guān)系，并通過與理論結(jié)果波形對比

4、來分析該仿真調(diào)制與解調(diào)過程的正確性及噪聲對FM 信號解調(diào)的影響。在課程設(shè)計中，系統(tǒng)開發(fā)平臺為Windows XP，在該平臺運行 MATLAB 程序完成了對FM調(diào)制和解調(diào)以及對疊加噪聲后解調(diào)結(jié)果的觀察和分析。 關(guān)鍵詞 ：FM 調(diào)制 解調(diào) MATLAB 噪聲 1AbstractModulation in a communication system has extremely important role. Through the modulation, not only can frequency shift, the modulation signal spectrum is moved to

5、 want, thus will position modulation signal is converted into suitable for transmission of the signal, and it has attune to the transmission effectiveness and reliability of transferring having very big effect, modulation method often determines a communication system performance. FM signals modulat

6、ion of nonlinear shift, belong to the spectrum of its demodulation also have coherent and incoherent demodulation in two ways. The course is designed to use simulation software for MATLAB 7.0, use of MATLAB integration environment of M files, write a program to realize FM signals of the simulation a

7、nalysis, and draw the baseband signal respectively, carrier signal, already adjustable signal and the time domain waveform frequency domain waveform; Further respectively to map out already adjustable signal after adding noise signal, synchronous demodulation signal demodulation before and after bas

8、eband signal temporal profile, Finally draw FM baseband signal modulation and demodulation system after the input/output SNR, and the relationship with the theoretical results wave contrast to analyze the simulation modulation and demodulation process accuracy and noise on FM signal demodulation of

9、the influence. In the curriculum design, system developing platform for Windows XP, on this platform operation MATLAB finished to FM modulation and demodulation of adding noise and after the observation and analysis results demodulation.Key words ：FM Modulation Demodulation MATLAB 7.0 noise1 緒論1.1課題

10、背景在當代社會中，信息的交換日益頻繁，隨著通信技術(shù)和計算機技術(shù)的發(fā)展及它們的密切結(jié)合，通信能克服對空間和實踐的限制，大量的、遠距離的信息傳遞和存取已成為可能。而調(diào)制能夠使通信的質(zhì)量獲得提高，因而有必要對調(diào)制技術(shù)進行研究。通信系統(tǒng)的一般模型可用圖1來概括：信息源接收設(shè)備信 道受信者發(fā)送設(shè)備（發(fā)送端）（接收端）噪聲源 圖1通信系統(tǒng)的一般模型信息源（簡稱信源）的作用是把各種信息轉(zhuǎn)換成原始信號。根據(jù)消息的種類不同信源分為模擬信源和數(shù)字信源。發(fā)送設(shè)備的作用產(chǎn)生適合傳輸?shù)男盘?，即使發(fā)送信號的特性和信道特性相匹配，具有抗噪聲的能力，并且具有足夠的功率滿足原距離傳輸?shù)男枨蟆Ｐ畔⒃春桶l(fā)送設(shè)備統(tǒng)稱為發(fā)送端。發(fā)送端

11、將信息直接轉(zhuǎn)換得到的較低頻率的原始電信號稱為基帶信號。通?；鶐盘柌灰酥苯釉谛诺乐袀鬏?。因此，在通信系統(tǒng)的發(fā)送端需將基帶信號的頻譜搬移（調(diào)制）到適合信道傳輸?shù)念l率范圍內(nèi)進行傳輸。這就是調(diào)制的過程。信號通過信道傳輸后，具有將信號放大和反變換功能的接收端將已調(diào)制的信號搬移（解調(diào)）到原來的頻率范圍，這就是解調(diào)的過程。信號在信道中傳輸?shù)倪^程總會受到噪聲的干擾，通信系統(tǒng)中沒有傳輸信號時也有噪聲，噪聲永遠存在于通信系統(tǒng)中。由于這樣的噪聲是疊加在信號上的，所以有時將其稱為加性噪聲。噪聲對于信號的傳輸是有害的，它能使模擬信號失真。在本仿真的過程中我們假設(shè)信道為高斯白噪聲信道。調(diào)制在通信系統(tǒng)中具有十分重要的作用

12、。一方面，通過調(diào)制可以把基帶信號的頻譜搬移到所希望的位置上去，從而將調(diào)制信號轉(zhuǎn)換成適合于信道傳輸或便于信道多路復用的已調(diào)信號。另一方面，通過調(diào)制可以提高信號通過信道傳輸時的抗干擾能力，同時，它還和傳輸效率有關(guān)。具體地講，不同的調(diào)制方式產(chǎn)生的已調(diào)信號的帶寬不同，因此調(diào)制影響傳輸帶寬的利用率?？梢?，調(diào)制方式往往決定一個通信系統(tǒng)的性能。在本仿真的過程中我們選擇用調(diào)頻調(diào)制方法進行調(diào)制。調(diào)制過程是一個頻譜搬移的過程，它是將低頻信號的頻譜搬移到載頻位置。而解調(diào)是將位于載頻的信號頻譜再搬回來，并且不失真地恢復出原始基帶信號。在本仿真的過程中我們選擇用非相干解調(diào)方法進行解調(diào)。1.2課程設(shè)計目的與意義通過FM系

13、統(tǒng)仿真的電路設(shè)計，掌握通信原理中模擬信號的調(diào)制與解調(diào)、數(shù)字基帶信號的傳輸、數(shù)字信號的調(diào)制與解調(diào)、模擬信號的抽樣、量化和編碼與信號的最佳接收等原理。應(yīng)用原理設(shè)計FM調(diào)制解調(diào)系統(tǒng)，并對其進行仿真。并能夠熟練應(yīng)用matlab語言編寫基本的通信系統(tǒng)程序。2 FM調(diào)制系統(tǒng)設(shè)計原理2.1調(diào)制信號的定義及產(chǎn)生 所謂頻率調(diào)制，指載波振幅不變，用調(diào)制信號m(t)去控制載波的瞬時頻率來實現(xiàn)其調(diào)制的一種方法。已調(diào)信號的瞬時頻率隨著調(diào)制信號的大小變化，更明確一點說是載波的瞬時頻率偏移隨m(t)成正比變化。調(diào)頻信號的產(chǎn)生主要有兩種方法：直接法和間接法。直接法師采用壓控振蕩器作為產(chǎn)生FM信號的調(diào)制器，使壓控振蕩器的瞬時頻

14、率隨調(diào)制信號的變化而呈線性變化。間接法也稱為倍頻法它是先用調(diào)制信號產(chǎn)生一個窄帶信號，然后將窄帶信號通過倍頻器得到帶寬信號。其兩種調(diào)制模型如圖2、圖3所示：m(t)SFM(t)FM調(diào)制器 圖2 FM直接調(diào)頻信號SFM(t)m(t)積分器PM調(diào)制器 圖3 FM間接調(diào)頻信號根據(jù)調(diào)制原理：設(shè)正弦型載波為 （2.1-1）是基帶調(diào)制信號，調(diào)制信號為： （2.1-2） 信號傳輸信道為高斯白噪聲信道，其功率為。 在調(diào)制時，調(diào)制信號的頻率控制載波的頻率變化，載波的瞬時頻偏隨調(diào)制信道信號成正比例變化，即： （）（）（2.1-3）9）（）（099是式中，為調(diào)頻靈敏度。此時相位偏移為 位偏移90 （2.1-4）S此時

15、可以得到調(diào)頻信號為 （2.1-5）式式中：為調(diào)頻指數(shù)式中 （2.1-6） 表示最大的相位偏移；其中的，為最大角偏移；，為最大頻偏。調(diào)頻波的有效帶寬為 FFM信號的頻譜理論值無窮大，但可根據(jù)調(diào)頻指數(shù)分為寬帶調(diào)頻和窄帶調(diào)頻。2.2 噪聲我們將信道中存在的不需要的電信號統(tǒng)稱為噪聲。通信系統(tǒng)中的噪聲是疊加在信號上的，沒有傳輸信號時通信系統(tǒng)中也有噪聲，噪聲永遠存在于通信系統(tǒng)中。噪聲可以看成是信道中的一種干擾，也稱為加性干擾，因為它是疊加在信號之上的。噪聲對于信號的傳輸是有害的，它能使模擬信號失真，是數(shù)字信號發(fā)生錯碼，并限制著信息的傳輸速率。按照來源分類，噪聲可以分為人為噪聲和自然噪聲兩大類。人為噪聲是有

16、人類的活動產(chǎn)生的，自然噪聲是自然界中存在的各種電磁波輻射，此外還有一種很重要的自然噪聲，即熱噪聲。熱噪聲來自一切電子型元器件中電子的熱運動。由于在一般的通信系統(tǒng)的工作頻率范圍內(nèi)熱噪聲的頻譜是均勻分布的，好像白光的頻譜在可見光的頻譜范圍內(nèi)均勻分布那樣，所以熱噪聲又常稱為白噪聲。由于熱噪聲是由大量自由電子的運動產(chǎn)生的，其統(tǒng)計特性服從高斯分布，故常將熱噪聲稱為高斯白噪聲。在本仿真的過程中我們假設(shè)信道為高斯白噪聲信道。2.3調(diào)制信號仿真圖在MATLAB程序中利用自建的積分函數(shù)實現(xiàn)對調(diào)制信號的間接調(diào)頻，其下為調(diào)制信號與已調(diào)信號分別在=0.5、=1和=3的情況下時域與頻域的圖形：圖4 調(diào)制信號與載波信號的

17、波形 圖5 =0.5時正弦波調(diào)制的各種波形 圖6 =0.5時三角波已調(diào)制的各種波形 圖7 =1時正弦波調(diào)制的各種波形 圖8 =1時三角波調(diào)制的各種波形 圖9 =3時正弦波已調(diào)制的各種波形圖10 =3時三角波已調(diào)制的各種波形3 FM解調(diào)信號模型3.1解調(diào)模型的建立調(diào)制信號的解調(diào)分為相干解調(diào)和非相干解調(diào)。其中相干解調(diào)是同步解調(diào)，僅僅適用于窄帶調(diào)頻信號，且需同步信號，因此其范圍受限；而非相干解調(diào)對NBFM信號和WBFM信號同時使用，因此是FM系統(tǒng)的主要調(diào)制方式。由于本課設(shè)要求使用同步解調(diào)，因此仿真過程針對窄帶信號進行調(diào)制與解調(diào)。NBFM相干解調(diào)原理圖如圖11所示：LPF微分BPF 圖11 NBFM信

18、號的相干解調(diào) 3.2解調(diào)過程分析根據(jù)式，設(shè)窄帶調(diào)頻信號 9090 （3.2-1）并設(shè)相干波，則相乘器的輸出為 （3.2-2）則通過低通濾波器取出其低頻分量 （3.2-3）再經(jīng)微分器，即得解調(diào)輸出 （3.2-4） 由此可見，相干解調(diào)可以恢復原調(diào)制信號。這種解調(diào)方法與線性中的相干解調(diào)一樣，要求本地載波與調(diào)制載波同步，否則將使解調(diào)信號失真。 matlab仿真結(jié)果如下所示：圖12 =0.5時調(diào)制下各種解調(diào)波形 圖13 =1時調(diào)制下各種解調(diào)波形圖14 =3時調(diào)制下各種解調(diào)波形4 高斯白噪聲信道特性4.1高斯白噪聲信道特性理論設(shè)正弦波通過加性高斯白噪聲到后的信號為 （4.1-1）其中，白噪聲其中，白噪聲的

19、取值的概率分布服從高斯分布。MATLAB本身自帶了標準高斯分布的內(nèi)部函數(shù)。函數(shù)產(chǎn)生的隨機序列服從均值為，方差的高斯分布。因為設(shè)正弦波通過加性高斯白噪聲信道后的信號如式（4.1-1），故其有用信號功率為 （4.1-2）噪聲功率為 （4.1-3）信噪比滿足公式 （4.1-4）則可得到公式 （4.1-5）通過此公式方便的設(shè)置高斯白噪聲的方差。4.2 輸入輸出信噪比關(guān)系曲線在本仿真過程中，選擇了5db、10db、15db、20db、25db五種不同信噪比以示區(qū)別，輸入輸出信噪比關(guān)系曲線圖如下所示：以下為調(diào)制信號與解調(diào)信號分別在=0.5、=1和=3的情況下，分別在有高斯噪聲和無高斯噪聲時的時域與頻域的圖

20、形：圖15 =0.5時正弦波調(diào)制下輸入輸出信噪比關(guān)系曲線圖16 =0.5時三角波調(diào)制下輸入輸出信噪比關(guān)系曲線圖17 =1時正弦波調(diào)制下輸入輸出信噪比關(guān)系曲線圖18 =1時三角波調(diào)制下輸入輸出信噪比關(guān)系曲線圖19 =3時正弦波調(diào)制下輸入輸出信噪比關(guān)系曲線圖20 =3時三角波調(diào)制下輸入輸出信噪比關(guān)系曲線5 調(diào)頻系統(tǒng)的抗噪聲性能分析從前面的分析可知，調(diào)頻信號的解調(diào)有相干解調(diào)和非相干解調(diào)兩種。相干解調(diào)僅適用于窄帶調(diào)頻信號，且需同步信號；而非相干解調(diào)適用于窄帶和寬帶調(diào)頻信號，而且不需同步信號，因而是FM系統(tǒng)的主要解調(diào)方式，所以這里僅僅討論非相干解調(diào)系統(tǒng)的抗噪聲性能，其分析模型如圖23所示。鑒頻器BPF限

21、幅器LPF解調(diào)器 圖23 調(diào)頻系統(tǒng)抗噪聲性能分析模型圖中帶通濾波器的作用是抑制信號帶寬以外的噪聲。是均值為零，單邊功率譜密度為的高斯白噪聲，經(jīng)過帶通濾波器后變?yōu)檎瓗Ц咚乖肼?。限幅器是為了消除接收信號在幅度上可能出現(xiàn)的畸變。設(shè)調(diào)頻信號為 （5-1）故其輸入功率為 （5-2） 輸入噪聲功率為 （5-3） 因此輸入信噪比為 （5-4）在大信噪比條件下，信號和噪聲的相互作用可以忽略，這時可以把信號和噪聲分開來算，這里，我們可以得到解調(diào)器的輸出信噪比 （5-5）上式中，為載波的振幅，為調(diào)頻器靈敏度，為調(diào)制信號的最高頻率，為噪聲單邊功率譜密度。我們?nèi)缛艨紤]為單一頻率余弦波時的情況，可得到解調(diào)器的制度增益

22、為 （5-6） 考慮在寬帶調(diào)頻時，信號帶寬為 （5-7） 則可以得到 （5-8） 可以看出，大信噪比時寬帶調(diào)頻系統(tǒng)的信噪比增益是很高的，它與調(diào)頻指數(shù)的立方成正比。可見，加大調(diào)頻指數(shù)，可使調(diào)頻系統(tǒng)的抗噪聲性能迅速改善。6 心得體會通過本次通信原理課程設(shè)計，我對信號的產(chǎn)生、調(diào)制與解調(diào)理論知識有了更深一步了解，此次課程設(shè)計鍛煉了我實踐能力，也懂得了很多將理論與實踐相結(jié)合的重要性，對matlab的操作更加熟悉，通過查找相關(guān)書籍與網(wǎng)絡(luò)，終于克服重重困難，懂得了如何用這個軟件將自己想要編寫的程序通過其本身自帶的函數(shù)表達出來，在這次設(shè)計過程中，發(fā)現(xiàn)我理論知識學的不夠扎實，通過這次實踐后，我對這部分知識掌握更

23、牢固了。 通過這次課程設(shè)計，我發(fā)現(xiàn)了許多以前從沒有發(fā)現(xiàn)的問題，并且學到了更多在現(xiàn)實運用中需要注意的東西，這些寶貴的經(jīng)驗都是無法從課本中獲得的，只有通過具體的實踐操作，才能印象深刻。一直以來，無論是老師還是學長都強調(diào)一定要把理論知識與實際相結(jié)合，本來以為是很簡單的事情，真正做起來卻讓我感到并沒有那么容易。雖然在制作此課程設(shè)計中遇到許多的困難，但是我學會了認真耐心分析，主要還是自己對所學的知識并沒有完全的消化，以至對FM原理都是似懂非懂，浪費了較多的時間，另外通過此次的課程設(shè)計，對MATLAB軟件的使用也有了更深一步的了解與掌握，特別是如何對基帶信號如何進行調(diào)制，時域到頻域變換，及加入噪聲方面有了

24、深入的了解，也為以后工作設(shè)計積累了經(jīng)驗，在這種實際操作中不斷熟悉。本次MATLAB課程設(shè)計，我通過查閱相關(guān)書籍、網(wǎng)上資料和一些仿真實例，實現(xiàn)了基本要求。通過這次的課程設(shè)計，使我更清楚的了解了理論與實際相結(jié)合的重要性，在以后的學習和生活中，我會更注重培養(yǎng)自身的實踐能力以及自學能力，并加強專業(yè)知識的系統(tǒng)性，努力扎實自己的專業(yè)基礎(chǔ)，進而提高自身的實際操作以及應(yīng)用能力。7 參考文獻1通信原理（第六版） 樊昌信等。 北京：國防工業(yè)出版社。2MATLAB7.0基礎(chǔ)教程 孫祥等。北京：清華大學出版社。3MATLAB在電子信息工程中的應(yīng)用 張德豐等。 北京：電子工業(yè)出版社4MATLAB7.X系統(tǒng)建模與仿真 張

25、亮等。 北京：人民郵電出版社附錄 Matlab源程序%*FM 調(diào)制*%dt=0.00001; %設(shè)定步長N=1000; %FFT長度Fs=1/dt; %設(shè)定采樣頻率n=0:N-1;t=n/Fs; %截止時間fc=30000; %設(shè)定載波頻率為要求的30KHZct=cos(2*pi*fc*t); %生成載波am=0.5; %設(shè)定調(diào)制信號幅度為0.5fm=300; %設(shè)定調(diào)制信號頻率為要求的300HZmt_sin=am*cos(2*pi*fm*t); %生成正弦波調(diào)制信號mt_tri=sawtooth(2*pi*fm*t,0.5);%生成三角波調(diào)制信號ma=0.5; %設(shè)定調(diào)頻指數(shù)（可更改為1或3

26、）int_mt_sin(1)=0; %定義int_mt_sin為mt_sin信號在時域上的積分for i=1:length(t)-1 int_mt_sin(i+1)=int_mt_sin(i)+mt_sin(i)*dt; %對mt_sin進行積分endint_mt_tri(1)=0; %定義int_mt_tri為mt_tri信號在時域上的積分for i=1:length(t)-1 int_mt_tri(i+1)=int_mt_tri(i)+mt_tri(i)*dt; %對mt_tri進行積分endsfm_sin=am*cos(2*pi*fc*t+ma*2*pi*fm*int_mt_sin);%

27、對正弦波調(diào)制，產(chǎn)生已調(diào)信號sfm_tri=am*cos(2*pi*fc*t+ma*2*pi*fm*int_mt_tri);%對三角波調(diào)制，產(chǎn)生已調(diào)信號figure(1);subplot(3,1,1);plot(t,mt_sin);%繪制正弦波調(diào)制信號的時域圖xlabel(時間t);title(正弦波調(diào)制信號的時域圖);subplot(3,1,2);plot(t,mt_tri);%繪制三角波調(diào)制信號的時域圖xlabel(時間t);title(三角波調(diào)制信號的時域圖);subplot(3,1,3);plot(t,ct);%繪制載波的時域圖xlabel(時間t);title(載波信號的時域圖);fi

28、gure(2);subplot(2,1,1);plot(t,sfm_sin);%繪制未加入噪聲時正弦波已調(diào)信號的時域圖xlabel(時間t);title(未加入噪聲時正弦波已調(diào)信號的時域圖);axis(0,0.001,-0.5,0.5);subplot(2,1,2);plot(t,sfm_tri);xlabel(時間t);title(未加入噪聲時三角波已調(diào)信號的時域圖);%繪制未加入噪聲時三角波已調(diào)信號的時域圖axis(0,0.001,-0.5,0.5); %* 添加高斯白噪聲*%yn_sin=awgn(sfm_sin,5);%在正弦波已調(diào)信號中加入高斯白噪聲yn_tri=awgn(sfm_t

29、ri,5);%在三角波已調(diào)信號中加入高斯白噪聲figure(1); %繪制加入噪聲前后正弦波已調(diào)信號的時域圖subplot(2,1,1);plot(t,sfm_sin);xlabel(時間t);title(未加入噪聲時正弦波已調(diào)信號的時域圖);axis(0,0.001,-0.5,0.5); subplot(2,1,2);plot(t,yn_sin);xlabel(時間t);title(加入高斯白噪聲時正弦波已調(diào)信號的時域圖);axis(0,0.001,-0.5,0.5); figure(2); %繪制加入噪聲前后三角波已調(diào)信號的時域圖subplot(2,1,1);plot(t,sfm_tri)

30、;xlabel(時間t);title(未加入噪聲時三角波已調(diào)信號的時域圖);axis(0,0.001,-0.5,0.5);subplot(2,1,2);plot(t,yn_tri);xlabel(時間t);title(加入高斯白噪聲時三角波已調(diào)信號的時域圖);axis(0,0.001,-0.5,0.5); %* FM解調(diào)*%yy_sin=demod(sfm_sin,fc,Fs,fm);yy_tri=demod(sfm_tri,fc,Fs,fm); yyn_sin=demod(yn_sin,fc,Fs,fm);yyn_tri=demod(yn_tri,fc,Fs,fm); figure(1)su

31、bplot(2,2,1);plot(t,yy_sin);xlabel(時間t);title(未加噪聲時正弦波解調(diào)信號時域圖);subplot(2,2,2);plot(t,yy_tri);xlabel(時間t);title(未入噪聲時三角波解調(diào)信號時域圖); subplot(2,2,3);plot(t,yyn_sin);xlabel(時間t);title(加噪聲時正弦波解調(diào)信號時域圖);subplot(2,2,4);plot(t,yyn_tri);xlabel(時間t);title(加噪聲時三角波解調(diào)信號時域圖); %*求傅里葉變換*%ff1=fft(mt_sin,N);%正弦調(diào)制信號傅里葉變換

32、mag1=abs(ff1); %取模f1_sin=(0:length(ff1)-1)*Fs/length(ff1);%頻率轉(zhuǎn)換ff2=fft(sfm_sin,N);%正弦已調(diào)信號傅里葉變換mag2=abs(ff2);f2_sin=(0:length(ff2)-1)*Fs/length(ff2);ff3=fft(mt_tri,N); %三角波調(diào)制信號傅里葉變換mag3=abs(ff3); %取模f3_tri=(0:length(ff3)-1)*Fs/length(ff3);%頻率轉(zhuǎn)換ff4=fft(sfm_tri,N); %三角波已調(diào)信號傅里葉變換mag4=abs(ff4);f4_tri=(0:

33、length(ff4)-1)*Fs/length(ff4); figure(1);subplot(2,1,1);plot(f1_sin,mag1);axis(0,500,0,600);xlabel(頻率HZ);title(正弦波調(diào)制信號頻譜圖);subplot(2,1,2);plot(f2_sin,mag2);axis(25000,35000,0,500);xlabel(頻率HZ);title(正弦波已調(diào)信號頻譜圖); figure(2)subplot(2,1,1);plot(f3_tri,mag3);axis(0,500,0,600);xlabel(頻率f);title(三角波調(diào)制信號頻譜圖);subplot(2,1,2);plot(f4_tri,mag4);axis(25000,35000,0,500);xlabel(頻率f);title(三角波已調(diào)信號頻譜圖); %*解調(diào)器抗噪聲性能分析*% yn_sin0=awgn(sfm_sin,5);yn_sin1=awgn(sfm_sin,10);yn_sin2=awgn(sfm_sin,15);yn_sin3=awgn(sfm_sin,20);yn_sin4=