FM調(diào)制解調(diào)系統(tǒng)設(shè)計與仿真lin.

1、課程設(shè)計任務(wù)書學(xué)生姓名：學(xué)生姓名： 專業(yè)班級：專業(yè)班級： 指導(dǎo)教師：指導(dǎo)教師： 工作單位：工作單位： 題題 目目: FM 信號的仿真分析 初始條件：初始條件：調(diào)制信號：分別為 300Hz 正弦信號和三角波信號；載波頻率：30kHz；解調(diào)方式：同步解調(diào)。要求完成的主要任務(wù)要求完成的主要任務(wù): 要求能夠熟練應(yīng)用 MATLAB 語言編寫基本的通信系統(tǒng)的應(yīng)用程序，進(jìn)行模擬調(diào)制系統(tǒng)，數(shù)字基帶信號的傳輸系統(tǒng)的建模、設(shè)計與仿真。所有的仿真用 MATLAB 程序?qū)崿F(xiàn)，系統(tǒng)經(jīng)過的信道都假設(shè)為高斯白噪聲信道。 畫出以下三種情況下調(diào)制信號、已調(diào)信號、解調(diào)信號的波形、頻譜以及解調(diào)器輸入 輸出信噪比的關(guān)系曲線；（調(diào)制指

2、數(shù)；調(diào)制指數(shù)=1；調(diào)制指數(shù)=3）時間安排：時間安排：1、 2013 年 12 月 19 日，布置課設(shè)具體實施計劃與課程設(shè)計報告格式的要求說明。 2、 2013 年 12 月 19 日至 2013 年 12 月 20 日，方案選擇和電路設(shè)計。 3、 2013 年 12 月 21 日至 2013 年 12 月 25 日，電路調(diào)試和設(shè)計說明書撰寫。 4、 2014 年 1 月 8 日，上交課程設(shè)計成果及報告，同時進(jìn)行答辯。指導(dǎo)教師簽名：指導(dǎo)教師簽名： 年年 月月 日日系主任（或責(zé)任教師）簽名：系主任（或責(zé)任教師）簽名： 年年 月月 日日目 錄摘摘 要要 I IABSTRACTABSTRACT III

3、I一通信系統(tǒng)介一通信系統(tǒng)介 1二二FMFM 調(diào)制解調(diào)系統(tǒng)設(shè)計調(diào)制解調(diào)系統(tǒng)設(shè)計 3 32.2 FM 調(diào)制模型的建立 32.3 FM 調(diào)制仿真結(jié)果 62.4 FM 解調(diào)模型的建立 62.5 解調(diào)過程分析 7高斯白噪聲信道特性8 92.8 調(diào)頻系統(tǒng)的抗噪聲性能分析 10三仿真實現(xiàn)三仿真實現(xiàn)12123.1 MATLAB 源代碼123.2MATLAB 仿真結(jié)果及分析 12四心得體會四心得體會1414五參考文獻(xiàn)五參考文獻(xiàn)1414摘 要 調(diào)制在通信系統(tǒng)中有十分重要的作用。通過調(diào)制，不僅可以進(jìn)行頻譜搬移，把調(diào)制信號的頻譜搬移到所希望的位置上，從而將調(diào)制信號轉(zhuǎn)換成適合于傳播的已調(diào)信號，而且它對系統(tǒng)的傳輸有效性和

4、傳輸?shù)目煽啃杂兄艽蟮挠绊?，調(diào)制方式往往決定了一個通信系統(tǒng)的性能。FM 信號的調(diào)制屬于頻譜的非線性搬移，它的解調(diào)也有相干和非相干解調(diào)兩種方式。 本課程設(shè)計使用的仿真軟件為 MATLAB，利用 MATLAB 集成環(huán)境下的 M 文件，編寫程序來實現(xiàn) FM 信號的仿真分析，并分別繪制出基帶信號、載波信號、已調(diào)信號的時域波形和頻域波形；再進(jìn)一步分別繪制出對已調(diào)信號疊加噪聲后信號、同步解調(diào)前信號和解調(diào)后基帶信號的時域波形；最后繪出 FM 基帶信號調(diào)制和解調(diào)系統(tǒng)后的輸入輸出信噪比的關(guān)系，并通過與理論結(jié)果波形對比來分析該仿真調(diào)制與解調(diào)過程的正確性及噪聲對 FM 信號解調(diào)的影響。在課程設(shè)計中，系統(tǒng)開發(fā)平臺為

5、Windows7，在該平臺運(yùn)行 MATLAB 程序完成了對FM 調(diào)制和解調(diào)以及對疊加噪聲后解調(diào)結(jié)果的觀察和分析。關(guān)鍵詞: FM 調(diào)制 解調(diào) MATLAB 噪聲Abstract Modulation has a very important role in the communication system . Modulation, can be used to not only move the spectrum , the spectrum of the modulated signal is moved to a desired position so as to convert the

6、modulated signal suitable for transmission of modulated signals,but also it significantly affect the validity of the transmission system and the transmission reliability . Often determines the performance of the modulation scheme of a communication system . FM modulated signal is Non-linear movement

7、 of spectrum . it also has the coherent demodulation and non-coherent demodulation of two ways. Simulation Software This course is designed for use as MATLAB, the use of M-files MATLAB integrated environment , write a program to achieve the simulation of FM signals , and were drawn out of the baseba

8、nd signal , the carrier signal, modulated signal domain waveform and frequency domain waveforms ; were drawn further modulated signal superimposed on the noise signal, a synchronization signal and demodulates the first time domain waveform of the demodulated baseband signal ; final output SNR plot F

9、M input baseband signal after the modulation and demodulation system of relations and, through comparison with the theoretical results to analyze waveforms and noise affect the correctness of the simulation process of modulation and demodulation of the FM signal demodulation . In the course design ,

10、 system development platform for Windows7, the platform running MATLAB program completed and FM modulation and demodulation noise superimposed on the demodulation results after observation and analysis .Keywords : FM Modulation Demodulation MATLAB Noise一通信系統(tǒng)簡介通信的目的是傳輸信息。通信系統(tǒng)的作用就是將信息從信息源發(fā)送到一個或多個目的地。對

11、于任何一個通信系統(tǒng)，均可視為由發(fā)送端、信道和接收端三大部分組成（如圖1 所示） 。信息源發(fā)送設(shè)備信 道接受設(shè)備信息源噪聲源發(fā)送端發(fā)送端接收端接收端信道信道圖 1 通信系統(tǒng)一般模型信息源（簡稱信源）的作用是把各種信息轉(zhuǎn)換成原始信號。發(fā)送設(shè)備的作用產(chǎn)生適合傳輸?shù)男盘?，即使發(fā)送信號的特性和信道特性相匹配，具有抗噪聲的能力，并且具有足夠的功率滿足原距離傳輸?shù)男枨蟆Ｐ畔⒃春桶l(fā)送設(shè)備統(tǒng)稱為發(fā)送端。發(fā)送端將信息直接轉(zhuǎn)換得到的較低頻率的原始電信號稱為基帶信號。通?；鶐盘柌灰酥苯釉谛诺乐袀鬏?。因此，在通信系統(tǒng)的發(fā)送端需將基帶信號的頻譜搬移（調(diào)制）到適合信道傳輸?shù)念l率范圍內(nèi)進(jìn)行傳輸。這就是調(diào)制的過程。信號通過信

12、道傳輸后，具有將信號放大和反變換功能的接收端將已調(diào)制的信號搬移（解調(diào)）到原來的頻率范圍，這就是解調(diào)的過程。信號在信道中傳輸?shù)倪^程總會受到噪聲的干擾，通信系統(tǒng)中沒有傳輸信號時也有噪聲，噪聲永遠(yuǎn)存在于通信系統(tǒng)中。噪聲對于信號的傳輸是有害的，它能使模擬信號失真。在本仿真的過程中我們假設(shè)信道為高斯白噪聲信道。調(diào)制在通信系統(tǒng)中具有十分重要的作用。一方面，通過調(diào)制可以把基帶信號的頻譜搬移到所希望的位置上去，從而將調(diào)制信號轉(zhuǎn)換成適合于信道傳輸或便于信道多路復(fù)用的已調(diào)信號。另一方面，通過調(diào)制可以提高信號通過信道傳輸時的抗干擾能力，同時，它還和傳輸效率有關(guān)。具體地講，不同的調(diào)制方式產(chǎn)生的已調(diào)信號的帶寬不同，因此

13、調(diào)制影響傳輸帶寬的利用率?？梢姡{(diào)制方式往往決定一個通信系統(tǒng)的性能。在本仿真的過程中我們選擇用調(diào)頻調(diào)制方法進(jìn)行調(diào)制。調(diào)制過程是一個頻譜搬移的過程，它是將低頻信號的頻譜搬移到載頻位置。而解調(diào)是將位于載頻的信號頻譜再搬回來，并且不失真地恢復(fù)出原始基帶信號。在本仿真的過程中我們選擇用相干解調(diào)方法進(jìn)行解調(diào)。角度調(diào)制及解調(diào)電路不同于頻譜線性搬移電路。它是用低頻信號去調(diào)制高頻振蕩的相角，或是從已調(diào)波中解出調(diào)制信號所進(jìn)行的頻譜變換，這種變換不是線性變換，而是非線性變換。因此，我們把角度調(diào)制及調(diào)角波的解調(diào)電路稱為頻譜非線性變換電路。FM在通信系統(tǒng)中的使用非常廣泛。FM 廣泛應(yīng)用于高保真音樂廣播、電視伴音信號的

14、傳輸、衛(wèi)星通信和蜂窩 系統(tǒng)等。本次課設(shè)中使用功能強(qiáng)大的 MATLAB 仿真軟件對 FM 信號進(jìn)行仿真分析。二FM 調(diào)制解調(diào)系統(tǒng)設(shè)計2.1 MATBLAB 簡介MATLAB 的名稱源自 Matrix Laboratory，它的首創(chuàng)者是在數(shù)值線性代數(shù)領(lǐng)域頗有影6max)(Kdttmf響的 Cleve Moler 博士，他也是生產(chǎn)經(jīng)營 MATLAB 產(chǎn)品的美國 Mathworks 公司的創(chuàng)始人之一。MATLAB 是一種科學(xué)計算軟件，專門以矩陣的形式處理數(shù)據(jù)。MATLAB 將高性能的數(shù)值計算和可視化集成在一起，并提供了大量的內(nèi)置函數(shù)，從而使其被廣泛地應(yīng)用于科學(xué)計算、控制系統(tǒng)、信息處理等領(lǐng)域的分析、仿真

15、和設(shè)計工作中。MATLAB 與其它計算機(jī)高級語言如 C，C+等相比，MATLAB 語言編程要簡潔得多，編程語句更加接近數(shù)學(xué)描述，可讀性好，其強(qiáng)大的圖形功能和可視化數(shù)據(jù)處理能力也是其它高級語言望塵莫及的。對于具有任何一門高級語言基礎(chǔ)的學(xué)生來說，學(xué)習(xí) MATLAB 十分容易。MATLAB 使得人們擺脫了常規(guī)計算機(jī)編程的繁瑣，讓人們能夠?qū)⒋蟛糠志ν度氲窖芯繂栴}的數(shù)學(xué)建模上?？梢哉f，應(yīng)用 MATLAB 這一數(shù)學(xué)計算和系統(tǒng)仿真的強(qiáng)大工具，可以使科學(xué)研究的效率得以成百倍的提高。 2.2 FM 調(diào)制模型的建立從頻率調(diào)制的相位與頻率關(guān)系可以看出，調(diào)頻信號可通過直接調(diào)頻和間接調(diào)頻兩種方法得到，所謂間接調(diào)頻就是

16、先對調(diào)制信號積分再調(diào)相而得到。同樣，調(diào)相信號也可以通過直接調(diào)相和間接調(diào)相兩種方法得到，間接調(diào)相就是先對調(diào)制信號進(jìn)行微分再進(jìn)行頻率調(diào)制。根據(jù)調(diào)制后已調(diào)信號的瞬時相位偏移的大小，可將角度調(diào)制分為寬帶調(diào)制(寬帶調(diào)頻和寬帶調(diào)相)和窄帶調(diào)制(窄帶調(diào)頻和窄帶調(diào)相)。如果調(diào)頻信號或調(diào)相信號的最大瞬時相位偏移保持在很小的范圍內(nèi)，一般小于 30即滿足條件： 時，則稱為窄帶調(diào)頻或窄帶調(diào)相。當(dāng)上述條件不滿足時，就稱為寬帶調(diào)頻或?qū)拵д{(diào)相。圖 2 FM 間接調(diào)制模型其中，m(t)為基帶調(diào)制信號，設(shè)調(diào)制信號為m(t)=cos(2*pi*fm*t) 設(shè)正弦載波為c(t)=cos(2*pi*fc*t) 在調(diào)制時，調(diào)制信號的頻

17、率去控制載波的頻率的變化，載波的瞬時頻偏隨調(diào)制信號 m(t)成正比例變化，即:dttmktttf00)()(式中，為調(diào)頻靈敏度。fK這時相位偏移為dttmktttf)()(00 則可得到調(diào)頻信號為 )(cos()(u00FMdttmktUttfcFM 的頻譜的計算： 可以看出 FM 的頻譜與)(fnmJ的值有關(guān)。其信號帶寬為 ) 1(2fmB調(diào)制信號產(chǎn)生的部分 M 文件如下：fc=30000; %載波頻率ct=cos(2*pi*fc*t); %載波信號am=1; %調(diào)制信號幅度mmfm1204202cos)(2)(4cos)(22cos)(2)()sincos(nfnffffftnmJmJtm

18、JtmJmJtm0)12(531) 12sin()(25sin)(23sin)(2sin)(2)sinsin(nfnfffftnmJtmJtmJtmJtmtnmJAtVktAtcnfnfcFM)(cos)(sincos)(U_fm=300; %調(diào)制信號頻率mt_sin=am*cos(2*pi*fm*t); %正弦波mt_tri=sawtooth(2*pi*fm*t,0.5); %三角波mf=1; %調(diào)制指數(shù)，可調(diào)int_mt_sin(1)=0; %正弦波積分 for i=1:length(t)-1 int_mt_sin(i+1)=int_mt_sin(i)+mt_sin(i)*dt; end

19、int_mt_tri(1)=0;for i=1:length(t)-1 %三角波積分 int_mt_tri(i+1)=int_mt_tri(i)+mt_tri(i)*dt;endsfm_sin=am*cos(2*pi*fc*t+mf*2*pi*fm*int_mt_sin);sfm_tri=am*cos(2*pi*fc*t+mf*2*pi*fm*int_mt_tri);2.3 FM 調(diào)制仿真結(jié)果圖 3 FM 調(diào)制2.4 FM 解調(diào)模型的建立 由前述公式知道輸入調(diào)頻信號為)(cos()(u00FMdttmktUttfc相干載波為c(t)=cos(2*pi*fc*t) 乘法器的作用是把調(diào)頻信號變成有

20、多種頻率的波的混合，乘法器輸出為 )2cos1 ()(212sin21)(tdttmKttscfcp經(jīng)低通濾波器后取出器低頻分量為 dttmKtsfd)(21)(在經(jīng)過微分器，即得出解調(diào)出的基帶信號： )(21)(0tmKtmf相干解調(diào)可以恢復(fù)出原來的基帶信號，而且要求本地載波與調(diào)制載波同步，否則會使解調(diào)信號失真。2.5 解調(diào)過程分析 在這次的解調(diào)設(shè)計中，采用 MATLAB 自帶的解調(diào)函數(shù) demod 對已調(diào)信號進(jìn)行解調(diào)。該函數(shù)使用形式如下所示：Z = demod(Y,Fc,Fs,Method)其中，Y 為輸入的調(diào)頻波函數(shù)，F(xiàn)c、Fs 分別為載波頻率和抽樣頻率，而且要滿足，Method 為調(diào)制

21、方式。Fc2Fs FM 解調(diào)的 M 文件部分代碼如下：y_sin=demod(sfm_sin,fc,Fs,fm);y_tri=demod(sfm_tri,fc,Fs,fm);yyn_sin=demod(yn_sin,fc,Fs,fm);yyn_tri=demod(yn_tri,fc,Fs,fm);通過 M 文件繪制出大、小信噪比解調(diào)的輸出波形如下：圖 4 小信噪比時的 FM 解調(diào)圖 5 大信噪比時的 FM 解調(diào)2.6 高斯白噪聲信道特性 設(shè)正弦波通過加性高斯白噪聲信道后的信號為( )cos()( )cr tAtn t 其中，白噪聲的取值的概率分布服從高斯分布。( )n tMATLAB 本身自帶

22、了標(biāo)準(zhǔn)高斯分布的內(nèi)部函數(shù) awgn。awgn 函數(shù)產(chǎn)生的隨機(jī)序列服從均值為，方差的高斯分布。0m 21正弦波通過加性高斯白噪聲信道后的信號為( )cos()( )cr tAtn t故其有用信號功率為22AS 噪聲功率為2N信噪比滿足公式SN1010log ()SBN則可得到公式22102 10BA我們可以通過這個公式方便的設(shè)置高斯白噪聲的方差。 在本仿真過程中，選擇了 1,3,5 三種不同調(diào)制指數(shù)以示區(qū)別，并繪制輸入輸出信噪比關(guān)系曲線，如下圖 6 所示。 mf=1 mf=3 mf=5圖 6 正弦波輸入輸出信噪比2.8 調(diào)頻系統(tǒng)的抗噪聲性能分析從前面的分析可知，調(diào)頻信號的解調(diào)有相干解調(diào)和非相干解

23、調(diào)兩種。相干解調(diào)僅適用于窄帶調(diào)頻信號，且需同步信號；而非相干解調(diào)適用于窄帶和寬帶調(diào)頻信號，而且不需同步信號，因而是 FM 系統(tǒng)的主要解調(diào)方式，所以這里僅僅討論非相干解調(diào)系統(tǒng)的抗噪聲性能，其分析模型如圖所示。圖 7 調(diào)頻系統(tǒng)抗噪聲性能分析模型圖中帶通濾波器的作用是抑制信號帶寬以外的噪聲。是均值為零，單邊功率譜密( )n t度為的高斯白噪聲，經(jīng)過帶通濾波器后變?yōu)檎瓗Ц咚乖肼?。限幅器是為了消除接0n( )in t收信號在幅度上可能出現(xiàn)的畸變。設(shè)調(diào)頻信號為( )cos( )tFMcfStAtKmd故其輸入功率為22iAS 輸入噪聲功率為ioFMNn B因此輸入信噪比為22iiFMSANB在大信噪比條

24、件下，信號和噪聲的相互作用可以忽略，這時可以把信號和噪聲分開來算，這里，我們可以得到解調(diào)器的輸出信噪比 222233( )8ofoomA K m tSNn f上式中，為載波的振幅，為調(diào)頻器靈敏度，為調(diào)制信號的最高頻率，AfKmf( )m t為噪聲單邊功率譜密度。on我們?nèi)缛艨紤]為單一頻率余弦波時的情況，可得到解調(diào)器的制度增益為( )m t22232AooFMfio miSNGmSn fN考慮在寬帶調(diào)頻時，信號帶寬為2(1)2()FMfmmBmfff 則可以得到23(1)FMffGmm可以看出，大信噪比時寬帶調(diào)頻系統(tǒng)的信噪比增益是很高的，它與調(diào)頻指數(shù)的立方成正比?？梢?，加大調(diào)頻指數(shù)，可使調(diào)頻系統(tǒng)

25、的抗噪聲性能迅速改善。fm三仿真實現(xiàn) 本次設(shè)計是在 MATLAB 平臺上，在前述 FM 通信系統(tǒng)的分析基礎(chǔ)上，進(jìn)行程序代碼的編寫，從而得到相關(guān)仿真波形，由此進(jìn)行分析討論。FM調(diào)制程序添加高斯白噪聲程序FM解調(diào)程序繪制頻譜圖程序顯示程序初始化程序圖 8 程序流程圖3.1 MATLAB 源代碼根據(jù)設(shè)計要求所編寫的源代碼已在在附頁中給出，詳情請看附頁 I。在 MATLAB 平臺上運(yùn)行程序得到各仿真圖，由于篇幅有限，已在附頁 II 給出。這里針對部分仿真結(jié)果進(jìn)行分析。（1）系統(tǒng)生成滿足條件的信號波形，包括正弦波和三角波。（2）在未加入高斯白噪聲時，調(diào)制系統(tǒng)輸出理想的 FM 已調(diào)信號，已調(diào)信號的頻率隨著

26、原信號的幅度大小而變。解調(diào)系統(tǒng)也能理想地解調(diào)出原始輸入信號，解調(diào)后的波形光滑，與原信號一致。（3）在加入較弱的高斯白噪聲時，調(diào)制系統(tǒng)輸出較為理想的 FM 已調(diào)信號，已調(diào)信號的頻率隨著原信號的幅度大小而變但并不是理想的。解調(diào)系統(tǒng)也能正確地解調(diào)出原始輸入信號，解調(diào)后的波形稍帶毛刺，與原信號基本一致。（4）在加入較強(qiáng)的高斯白噪聲時，調(diào)制系統(tǒng)輸出不太理想的 FM 已調(diào)信號，已調(diào)信號的頻率與原信號的幅度大小關(guān)系不明顯。解調(diào)后的波形有失真現(xiàn)象。（5）輸入輸出信噪比呈正向線性關(guān)系，輸入信噪比越大，輸出信噪比也越大?？偸觯篎M 調(diào)制解調(diào)系統(tǒng)設(shè)計成功。結(jié)果表明：信道噪聲很大程度上影響著調(diào)制解調(diào)系統(tǒng)性能，實際運(yùn)用

27、中應(yīng)盡量降低信道噪聲；為了少失真地得到原始信號，在提高系統(tǒng)性能的基礎(chǔ)上，還應(yīng)增大輸入信噪比。四心得體會課程設(shè)計是培養(yǎng)學(xué)生綜合運(yùn)用所學(xué)知識,發(fā)現(xiàn),提出,分析和解決實際問題,鍛煉實踐能力的重要環(huán)節(jié),是對學(xué)生實際工作能力的具體訓(xùn)練和考察過程. 回顧起此次通信原理課程設(shè)計，至今我仍感慨頗多，的確，從選題到定稿，從理論到實踐，在整整一星期的日子里，可以說得是苦多于甜，但是可以學(xué)到很多很多的的東西，同時不僅可以鞏固了以前所學(xué)過的知識，而且學(xué)到了很多在書本上所沒有學(xué)到過的知識。通過這次課程設(shè)計使我懂得了理論與實際相結(jié)合是很重要的，只有理論知識是遠(yuǎn)遠(yuǎn)不夠的，只有把所學(xué)的理論知識與實踐相結(jié)合起來，從理論中得出結(jié)

28、論，才能真正為社會服務(wù)，從而提高自己的實際動手能力和獨立思考的能力。在設(shè)計的過程中遇到問題，可以說得是困難重重，這畢竟第一次做的，難免會遇到過各種各樣的問題，同時在設(shè)計的過程中發(fā)現(xiàn)了自己的不足之處，對以前所學(xué)過的知識理解得不夠深刻，掌握得不夠牢固通過這次課程設(shè)計之后，一定把以前所學(xué)過的知識重新溫故。這次課程設(shè)計終于順利完成了，在設(shè)計中遇到了很多問題，最后在不懈的努力下終于迎刃而解。感謝這次旅程，我又成長不少！ 五參考文獻(xiàn)1通信原理（第六版） 樊昌信等。 北京：國防工業(yè)出版社。2MATLAB7.0 在數(shù)字信號處理中的應(yīng)用 羅軍輝等。北京：機(jī)械工業(yè)出版社。3MATLAB 程序設(shè)計教程 劉衛(wèi)國等。

29、北京：中國水利水電出版社。附頁 I仿真波形圖 附頁 II源代碼%*本程序主要實現(xiàn) FM 調(diào)制與解調(diào)系統(tǒng)設(shè)計*%*以正弦波和三角波進(jìn)行設(shè)計*%*并考慮信道中高斯白噪聲的影響*%*初始化*dt=0.00001; %積分步長N=1000; %采樣點數(shù)Fs=1/dt; %抽樣頻率 Fs=100000hz,滿足抽樣定理n=0:N-1;t=n/Fs; %*信號設(shè)定*fc=30000; %載波頻率ct=cos(2*pi*fc*t); %載波信號am=1; %調(diào)制信號幅度fm=300; %調(diào)制信號頻率mt_sin=am*cos(2*pi*fm*t); %正弦波mt_tri=sawtooth(2*pi*fm*t

30、,0.5); %三角波mf=3; %調(diào)制指數(shù)%*調(diào)制模塊*int_mt_sin(1)=0; %正弦波積分 for i=1:length(t)-1 int_mt_sin(i+1)=int_mt_sin(i)+mt_sin(i)*dt; end int_mt_tri(1)=0;for i=1:length(t)-1 %三角波積分 int_mt_tri(i+1)=int_mt_tri(i)+mt_tri(i)*dt;endsfm_sin=am*cos(2*pi*fc*t+mf*2*pi*fm*int_mt_sin);sfm_tri=am*cos(2*pi*fc*t+mf*2*pi*fm*int_mt

31、_tri);%*加高斯白噪聲處理*yn_sin=awgn(sfm_sin,50);yn_tri=awgn(sfm_tri,50);yn_sin1=awgn(sfm_sin,10);yn_tri1=awgn(sfm_tri,10);yn_sin2=awgn(sfm_sin,60);yn_tri2=awgn(sfm_tri,60);yn_sin3=awgn(sfm_sin,80);yn_tri3=awgn(sfm_tri,80);yn_sin4=awgn(sfm_sin,40);yn_tri4=awgn(sfm_sin,40);%*解調(diào)模塊*y_sin=demod(sfm_sin,fc,Fs,fm

32、);y_tri=demod(sfm_tri,fc,Fs,fm);yyn_sin=demod(yn_sin,fc,Fs,fm);yyn_tri=demod(yn_tri,fc,Fs,fm);%*求正弦信號頻譜*ff1=fft(mt_sin,N);mag1=abs(ff1);f1_sin=(0:length(ff1)-1)*Fs/length(ff1);%*求正弦已調(diào)信號頻譜*ff2=fft(sfm_sin,N);mag2=abs(ff2);f2_sin=(0:length(ff2)-1)*Fs/length(ff2);%*求三角波信號頻譜*ff3=fft(mt_tri,N);mag3=abs(f

33、f3);f3_tri=(0:length(ff3)-1)*Fs/length(ff3);%*求三角波已調(diào)信號頻譜*ff4=fft(sfm_tri,N);mag4=abs(ff4);f4_tri=(0:length(ff4)-1)*Fs/length(ff4);%*加噪解調(diào)*yyn_sin1=demod(yn_sin1,30000,Fs,fm);yyn_sin2=demod(yn_sin2,30000,Fs,fm);yyn_sin3=demod(yn_sin3,30000,Fs,fm);yyn_sin4=demod(yn_sin4,30000,Fs,fm);yyn_tri1=demod(yn_t

34、ri1,30000,Fs,fm);yyn_tri2=demod(yn_tri2,30000,Fs,fm);yyn_tri3=demod(yn_tri3,30000,Fs,fm);yyn_tri4=demod(yn_tri4,30000,Fs,fm);%*輸入輸出信噪比*%*不同噪聲強(qiáng)度下的輸入輸出比*dy1=yn_sin-sfm_sin; snr1=var(sfm_sin)/var(dy1)/Fs/10;%snr1=var(sfm_sin)/var(dy1);dy2=yyn_sin-y_sin;snr2=var(y_sin)/var(dy2)/Fs*400;%snr2=var(y_sin)/v

35、ar(dy2);dy11=yn_sin1-sfm_sin;snr11=var(sfm_sin)/var(dy11)/Fs/10;%snr11=var(sfm_sin)/var(dy11);dy21=yyn_sin1-y_sin;snr21=var(y_sin)/var(dy21)/Fs*400;%snr21=var(y_sin)/var(dy21);dy12=yn_sin2-sfm_sin;%snr12=var(sfm_sin)/var(dy12);snr12=var(sfm_sin)/var(dy12)/Fs/10;dy22=yyn_sin2-y_sin;%snr22=var(y_sin)/

36、var(dy22);snr22=var(y_sin)/var(dy22)/Fs*400;dy13=yn_sin3-sfm_sin;%snr13=var(sfm_sin)/var(dy13);snr13=var(sfm_sin)/var(dy13)/Fs/10;dy23=yyn_sin3-y_sin;%snr23=var(y_sin)/var(dy23);snr23=var(y_sin)/var(dy23)/Fs*400;dy14=yn_sin4-sfm_sin;snr14=var(sfm_sin)/var(dy14)/Fs/10;%snr14=var(sfm_sin)/var(dy14);dy

37、24=yyn_sin4-y_sin;snr24=var(y_sin)/var(dy24)/Fs*400;%snr24=var(y_sin)/var(dy24);in=snr1,snr11,snr12,snr13,snr14;out=snr2,snr21,snr22,snr23,snr24;%*顯示模塊*%*調(diào)制信號及載波的顯示*figure(1)subplot(3,1,1);plot(t,mt_sin)xlabel(time);title(正弦波調(diào)制信號時域圖);subplot(3,1,2);plot(t,mt_tri)xlabel(time);title(三角波調(diào)制信號時域圖);subplo

38、t(3,1,3);plot(t,ct)xlabel(time);title(載波時域圖);%*已調(diào)波的顯示*figure(2)subplot(2,1,1);plot(t,sfm_sin)xlabel(time);title(未加入噪聲正弦波已調(diào)信號時域圖);axis(0 0.007 -2 2);figure(2)subplot(2,1,2);plot(t,sfm_tri)xlabel(time);title(未加入噪聲三角波已調(diào)信號時域圖);axis(0 0.007 -2 2);%*加入高斯白噪聲已調(diào)波的顯示*figure(3)subplot(2,1,1);plot(t,yn_sin)xlabel(time);title(加入高斯白噪聲正弦波已調(diào)信號時域圖);