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

1、貴州大學(xué)明德學(xué)院高頻電子線路課程設(shè)計(jì)報(bào)告題目模擬角度調(diào)制系統(tǒng)學(xué)院：明德學(xué)院專業(yè)：電子信息工程班級：學(xué)號：姓名：周科遠(yuǎn)指導(dǎo)老師：2012年1月1日高頻電子線路課程設(shè)計(jì)任務(wù)書一、課程設(shè)計(jì)的目的高頻電子線路課程設(shè)計(jì)是專業(yè)實(shí)踐環(huán)節(jié)之一， 是學(xué)習(xí)完 高頻電子線路 課程后進(jìn) 行的一次全面的綜合練習(xí)。其目的讓學(xué)生掌握高頻電子線路的基本原理極其構(gòu)造和運(yùn) 用，特別是理論聯(lián)系實(shí)踐，提高學(xué)生的綜合應(yīng)用能力。二、課程設(shè)計(jì)任務(wù)課程設(shè)計(jì)一、高頻放大器課程設(shè)計(jì)二、高頻振蕩器課程設(shè)計(jì)三、模擬線性調(diào)制系統(tǒng)課程設(shè)計(jì)四、模擬角度調(diào)制系統(tǒng)課程設(shè)計(jì)五、數(shù)字信號的載波傳輸課程設(shè)計(jì)六、通信系統(tǒng)中的鎖相環(huán)調(diào)制系統(tǒng)共 6 個(gè)課題選擇，學(xué)生任選

2、一個(gè)課題為自己的課程設(shè)計(jì)題目，獨(dú)立完成；具體內(nèi)容按方向分別進(jìn)行，不能有雷同；任務(wù)包括原理介紹、系統(tǒng)仿真、波形分 析等；要求按學(xué)校統(tǒng)一的課程設(shè)計(jì)規(guī)范撰寫一份設(shè)計(jì)說明書。三、課程設(shè)計(jì)時(shí)間課程設(shè)計(jì)總時(shí)間 1 周（ 5個(gè)工作日）四、課程設(shè)計(jì)說明書撰寫規(guī)范1、在完成任務(wù)書中所要求的課程設(shè)計(jì)作品和成果外，要撰寫課程設(shè)計(jì)說明書1 份。課程設(shè)計(jì)說明書須每人一份，獨(dú)立完成。2、設(shè)計(jì)說明書應(yīng)包括封面、任務(wù)書、目錄、摘要、正文、參考文獻(xiàn)（資料） 等內(nèi)容，以及附圖或附件等材料。3、題目字體用小三，黑體，正文字體用五號字，宋體，小標(biāo)題用四號及小四， 宋體，并用 A4 紙打印。目錄摘要IABSTRACT II一課程設(shè)計(jì)的

3、目的與要求 . 11.1課程設(shè)計(jì)的目的 11.2課程設(shè)計(jì)的要求 1二.FM調(diào)制解調(diào)系統(tǒng)設(shè)計(jì) 22.1 FM調(diào)制模型的建立 32.2 調(diào)制過程分析 32.3 FM解調(diào)模型的建立 42.4 解調(diào)過程分析 52.5 高斯白噪聲信道特性 62.6調(diào)頻系統(tǒng)的抗噪聲性能分析 9三仿真實(shí)現(xiàn) 103.1 MATLA源代碼 113.2 仿真結(jié)果 15四心得體會(huì) 18五參考文獻(xiàn) 192 / 25 第 6 頁 共 25 頁摘要FM在通信系統(tǒng)中的使用非常廣泛。FM廣泛應(yīng)用于高保真音樂廣 播、電視伴音信號的傳輸、衛(wèi)星通信和蜂窩電話系統(tǒng)等。本設(shè)計(jì)主要是利用MATLABS成環(huán)境下的M文件，編寫程序來實(shí) 現(xiàn)FM調(diào)制與解調(diào)過程

4、，并分別繪制出基帶信號，載波信號，已調(diào)信 號的時(shí)域波形； 再進(jìn)一步分別繪制出對已調(diào)信號疊加噪聲后信號， 相 干解調(diào)后信號和解調(diào)基帶信號的時(shí)域波形；最后繪出FM基帶信號通過上述信道和調(diào)制和解調(diào)系統(tǒng)后的誤碼率與信噪比的關(guān)系， 并通過與 理論結(jié)果波形對比來分析該仿真調(diào)制與解調(diào)系統(tǒng)的正確性及噪聲對 信號解調(diào)的影響。在課程設(shè)計(jì)中，系統(tǒng)開發(fā)平臺(tái)為 Windows Vista ， 使用工具軟件為MATLAB 7.0在該平臺(tái)運(yùn)行程序完成了對 FM調(diào)制和 解調(diào)以及對疊加噪聲后解調(diào)結(jié)果的觀察。 通過該課程設(shè)計(jì)， 達(dá)到了實(shí) 現(xiàn)FM信號通過噪聲信道，調(diào)制和解調(diào)系統(tǒng)的仿真目的。關(guān)鍵詞 FM;調(diào)制；解調(diào)；MATLAB 7

5、.0 ;噪聲ABSTRACTKEYWORDSFM: in the communication system is widely used. FM is widely used in high fidelity music radio, television signal transmission, satellite communication and cellular telephone system.The main design is the use of MATLAB integrated environment M documents, write a program to achi

6、eve FM modulation and demodulation process, and draw the baseband signal, a carrier signal, modulated signal time-domain waveform; further are plotted on a modulated signal superimposed noise signal, the coherent demodulated signal and the demodulated baseband signal in the time domain waveform; fin

7、ally draw FM baseband signal through the channel and modulation and demodulation system bit error rate and signal to noise ratio relationship, and through the theoretical results and the simulation waveform comparison of modulation and demodulation system is correct and the influence of noise on sig

8、nal demodulation. In the course of design, system development platform for Windows Vista, the use of tools for MATLAB 7. In the platform to run the program completed the FM modulation and demodulation of superimposed noise after the demodulation results of observation. Through the curriculum design,

9、 to achieve the realization of the FM signal through a noisy channel, modulation and demodulation system simulation purpose.課程設(shè)計(jì)的目的與要求1.1 課程設(shè)計(jì)的目的通過FM調(diào)制解調(diào)系統(tǒng)設(shè)計(jì)與防真的課程設(shè)計(jì)，掌握通信原理中模擬信 號的調(diào)制和解調(diào)、數(shù)字基帶信號的傳輸、 數(shù)字信號的調(diào)制和解調(diào)， 模擬信號的抽 樣、量化和編碼與信號的最佳接收等原理。應(yīng)用原理設(shè)計(jì) FM 調(diào)制解調(diào)系統(tǒng)，并 對其進(jìn)行防真。1.2 課程設(shè)計(jì)的要求要求能夠熟練應(yīng)用MATLA語言編寫基本的通信系統(tǒng)的應(yīng)用程序

10、，進(jìn)行模擬 調(diào)制系統(tǒng)， 數(shù)字基帶信號的傳輸系統(tǒng)的建模、 設(shè)計(jì)與仿真。 所有的仿真用 MATLAB 程序?qū)崿F(xiàn)（即只能用代碼的形式，不能用SIMULINK實(shí)現(xiàn)），系統(tǒng)經(jīng)過的信道都假 設(shè)為高斯白噪聲信道。模擬調(diào)制要求用程序畫出調(diào)制信號，載波，已調(diào)信號、解 調(diào)信號的波形，數(shù)字調(diào)制要求畫出誤碼率隨信噪比的變化曲線。FM調(diào)制解調(diào)系統(tǒng)設(shè)計(jì)通信的目的是傳輸信息。通信系統(tǒng)的作用就是將信息從信息源發(fā)送到一個(gè)或 多個(gè)目的地。對于任何一個(gè)通信系統(tǒng)，均可視為由發(fā)送端、信道和接收端三大部 分組成（如圖1所示）0圖1通信系統(tǒng)一般模型信息源（簡稱信源）的作用是把各種信息轉(zhuǎn)換成原始信號。根據(jù)消息的種類 不同信源分為模擬信源和數(shù)

11、字信源。發(fā)送設(shè)備的作用產(chǎn)生適合傳輸?shù)男盘枺词?發(fā)送信號的特性和信道特性相匹配，具有抗噪聲的能力，并且具有足夠的功率滿 足原距離傳輸?shù)男枨?。信息源和發(fā)送設(shè)備統(tǒng)稱為發(fā)送端。發(fā)送端將信息直接轉(zhuǎn)換得到的較低頻率的原始電信號稱為基帶信號。通?；鶐盘柌灰酥苯釉谛诺乐袀鬏敗Ｒ虼?，在通信系統(tǒng)的發(fā)送端需將基帶信號的頻譜 搬移（調(diào)制）到適合信道傳輸?shù)念l率范圍內(nèi)進(jìn)行傳輸。這就是調(diào)制的過程。信號通過信道傳輸后，具有將信號放大和反變換功能的接收端將已調(diào)制的信 號搬移（解調(diào)）到原來的頻率范圍，這就是解調(diào)的過程。信號在信道中傳輸?shù)倪^程總會(huì)受到噪聲的干擾， 通信系統(tǒng)中沒有傳輸信號時(shí) 也有噪聲，噪聲永遠(yuǎn)存在于通信系統(tǒng)中。

12、由于這樣的噪聲是疊加在信號上的， 所 以有時(shí)將其稱為加性噪聲。噪聲對于信號的傳輸是有害的，它能使模擬信號失真。 在本仿真的過程中我們假設(shè)信道為高斯白噪聲信道。調(diào)制在通信系統(tǒng)中具有十分重要的作用。一方面，通過調(diào)制可以把基帶信號 的頻譜搬移到所希望的位置上去，從而將調(diào)制信號轉(zhuǎn)換成適合于信道傳輸或便于 信道多路復(fù)用的已調(diào)信號。另一方面，通過調(diào)制可以提高信號通過信道傳輸時(shí)的 抗干擾能力，同時(shí)，它還和傳輸效率有關(guān)。具體地講，不同的調(diào)制方式產(chǎn)生的已 調(diào)信號的帶寬不同，因此調(diào)制影響傳輸帶寬的利用率?？梢姡{(diào)制方式往往決定一個(gè)通信系統(tǒng)的性能。在本仿真的過程中我們選擇用調(diào)頻調(diào)制方法進(jìn)行調(diào)制。調(diào)制過程是一個(gè)頻譜搬

13、移的過程，它是將低頻信號的頻譜搬移到載頻位置。 而解調(diào)是將位于載頻的信號頻譜再搬回來，并且不失真地恢復(fù)出原始基帶信號。 在本仿真的過程中我們選擇用非相干解調(diào)方法進(jìn)行解調(diào)。10 / 25第20頁共25頁2.1 FM調(diào)制模型的建立圖2 FM調(diào)制模型其中，m(t)為基帶調(diào)制信號，設(shè)調(diào)制信號為m(t) Acos(2 fnt)設(shè)正弦載波為c(t) cos(2 fct)信號傳輸信道為高斯白噪聲信道，其功率為2 02.2 調(diào)制過程分析在調(diào)制時(shí)，調(diào)制信號的頻率去控制載波的頻率的變化， 載波的瞬時(shí)頻偏隨調(diào)制信號m(t)成正比例變化，即d dtK fm (t)式中，Kf為調(diào)頻靈敏度rad.(s?V)這時(shí)相位偏移為

14、(t)Kfm( )d則可得到調(diào)頻信號為$M(t) Acos ctKf m( )ddt=0.001;t=0:dt:1.5;am=15; fm=15;mt=am*cos(2*pi*fm*t);fc=50;%生成調(diào)制信號%設(shè)定載波頻率-可更改調(diào)制信號產(chǎn)生的M文件:%設(shè)定時(shí)間步長%產(chǎn)生時(shí)間向量 %設(shè)定調(diào)制信號幅度J可更改%設(shè)定調(diào)制信號頻率-可更改ct=cos(2*pi*fc*t);kf=10;in t_mt(1)=0;for i=1:le ngth(t)-1in t_mt(i+1)=i nt_mt(i)+mt(i)*dt; endsfm=am*cos(2*pi*fc*t+2*pi*kf*i nt_mt

15、);%生成載波%設(shè)定調(diào)頻指數(shù)%對mt進(jìn)行積分%調(diào)制，產(chǎn)生已調(diào)信號載波的時(shí)域圖時(shí)間t已調(diào)信號的時(shí)域圖圖3 FM調(diào)制2.3 FM解調(diào)模型的建立調(diào)制信號的解調(diào)分為相干解調(diào)和非相干解調(diào)兩種。相干解調(diào)僅僅適用于窄帶 調(diào)頻信號，且需同步信號，故應(yīng)用范圍受限；而非相干解調(diào)不需同步信號，且對 于NBFM信號和WBFM信號均適用，因此是FM系統(tǒng)的主要解調(diào)方式。在本仿 真的過程中我們選擇用非相干解調(diào)方法進(jìn)行解調(diào)。圖4 FM解調(diào)模型非相干解調(diào)器由限幅器、鑒頻器和低通濾波器等組成，其方框圖如圖5所示。 限幅器輸入為已調(diào)頻信號和噪聲，限幅器是為了消除接收信號在幅度上可能出現(xiàn) 的畸變；帶通濾波器的作用是用來限制帶外噪聲，

16、使調(diào)頻信號順利通過。鑒頻器 中的微分器把調(diào)頻信號變成調(diào)幅調(diào)頻波， 然后由包絡(luò)檢波器檢出包絡(luò)，最后通過 低通濾波器取出調(diào)制信號。2.4解調(diào)過程分析設(shè)輸入調(diào)頻信號為tS(t) SFM(t) Acos(ct Kf m( )d )微分器的作用是把調(diào)頻信號變成調(diào)幅調(diào)頻波。微分器輸出為dSFM (t)dtc Kfm(t) sin( cttKfm( )d )包絡(luò)檢波的作用是從輸出信號的幅度變化中檢出調(diào)制信號。包絡(luò)檢波器輸出為So(t)Kdc Km(t)Kd c KdKm(t)Kd稱為鑒頻靈敏度(VHz)，是已調(diào)信號單位頻偏對應(yīng)的調(diào)制信號的幅度，經(jīng)低通濾波器后加隔直流電容，隔除無用的直流，得mo(t) KdK

17、fm(t)微分器通過程序?qū)崿F(xiàn)，代碼如下：for i=1:le ngth(t)-1%接受信號通過微分器處理diff_ nsfm(i)=( nsfm(i+1)-nsfm(i)./dt;enddiff_nsfmn = abs(hilbert(diff_nsfm);%hilbert 變換，求絕對值得到瞬時(shí)幅度(包絡(luò)檢波) 通過M文件繪制出兩種不同信噪比解調(diào)的輸出波形如下：調(diào)制信號的時(shí)域圖50.511.5時(shí)間t無噪聲條件下已調(diào)信號的時(shí)域圖5 0 -/% JF-5 無噪聲條件下解調(diào)信號的時(shí)域圖00.511.5時(shí)間t圖5 FM解調(diào)2.5高斯白噪聲信道特性設(shè)正弦波通過加性高斯白噪聲信道后的信號為r(t) Ac

18、os( ct) n(t)其中，白噪聲n(t)的取值的概率分布服從高斯分布。MATLAB本身自帶了標(biāo)準(zhǔn)高斯分布的內(nèi)部函數(shù) randn。randn函數(shù)產(chǎn)生的隨 機(jī)序列服從均值為m 0，方差21的高斯分布。正弦波通過加性高斯白噪聲信道后的信號為r(t) Acos( ct) n(t)故其有用信號功率為A2S A2噪聲功率為N 2信噪比S N滿足公式則可得到公式B 10logio(SN)AB2?1010我們可以通過這個(gè)公式方便的設(shè)置高斯白噪聲的方差。在本仿真過程中，我們選擇了 10db和30db兩種不同信噪比以示區(qū)別，其時(shí) 域圖如圖7和圖&5-5050-5調(diào)制信號的時(shí)域圖00.511.5時(shí)間t無噪聲條件

19、下已調(diào)信號的時(shí)域圖051.50.51時(shí)間t無噪聲條件下解調(diào)信號的時(shí)域圖000.511.5時(shí)間t-5圖6無噪聲條件下已調(diào)信號的時(shí)域圖50-550000 -5000200-2050-550000 -500050-5調(diào)制信號的時(shí)域圖00.511.5時(shí)間t含小信噪比高斯白噪聲已調(diào)信號的時(shí)域圖050010001500時(shí)間t含小信噪比高斯白噪聲解調(diào)信號的時(shí)域圖圖7含小信噪比高斯白噪聲已調(diào)信號的時(shí)域圖調(diào)制信號的時(shí)域圖00.511.5時(shí)間t含大信噪比高斯白噪聲已調(diào)信號的時(shí)域圖050010001500時(shí)間t含大信噪比高斯白噪聲解調(diào)信號的時(shí)域圖00.511.5時(shí)間t圖8含大信噪比高斯白噪聲已調(diào)信號的時(shí)域圖2.6調(diào)

20、頻系統(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)的抗噪聲性能，其分析模型如圖 9所示。廠1解調(diào)器圖9調(diào)頻系統(tǒng)抗噪聲性能分析模型圖中帶通濾波器的作用是抑制信號帶寬以外的噪聲。n(t)是均值為零，單邊功率譜密度為no的高斯白噪聲，經(jīng)過帶通濾波器后變?yōu)檎瓗Ц咚乖肼昻i(t)。限幅器是為了消除接收信號在幅度上可能出現(xiàn)的畸變設(shè)調(diào)頻信號為Sfm (t)Acos(tct Kfm( )d )故其輸入功率為SiA

21、22輸入噪聲功率為NinoBFM因此輸入信噪比為SA2Ni2Bfm在大信噪比條件下，信號和噪聲的相互作用可以忽略， 這時(shí)可以把信號和噪聲分開來算，這里，我們可以得到解調(diào)器的輸出信噪比So3A2K：m2(t)No 8 2nof；上式中，A為載波的振幅，Kf為調(diào)頻器靈敏度，fm為調(diào)制信號m(t)的最高頻率，no為噪聲單邊功率譜密度。我們?nèi)缛艨紤]m(t)為單一頻率余弦波時(shí)的情況，可得到解調(diào)器的制度增益為G FM2mfA2 2n ofm考慮在寬帶調(diào)頻時(shí)，信號帶寬為BFM2(mf 1)fm 2( ffm)則可以得到2Gfm 3mf (mf 1)可以看出，大信噪比時(shí)寬帶調(diào)頻系統(tǒng)的信噪比增益是很高的，它與調(diào)

22、頻指數(shù)的立方成正比?？梢?，加大調(diào)頻指數(shù) mf，可使調(diào)頻系統(tǒng)的抗噪聲性能迅速改善。圖10程序流程圖3.1 MATLAB源代碼%FM調(diào)制解調(diào)系統(tǒng).m%頻率調(diào)制與解調(diào)的 Matlab演示源程序%可以任意改原調(diào)制信號函數(shù)m(t)%信息工程陳麗丹07323202. * * * * * * * * * * * * * * * .%*初始化*echo off close all clear all clc %*. * * * * * * * * * * * * * * * .*FM調(diào)制*dt=0.001;t=0:dt:1.5;%設(shè)定時(shí)間步長%產(chǎn)生時(shí)間向量% 設(shè)定調(diào)制信號幅度 % 設(shè)定調(diào)制信號頻率% 生成調(diào)制

23、信號% 設(shè)定載波頻率% 生成載波% 設(shè)定調(diào)頻指數(shù)%求信號m（t）的積分% 調(diào)制，產(chǎn)生已調(diào)信號% 調(diào)制信號%*sn1=10;sn2=30;添加高斯白噪聲 *sn=0;db=amA2/(2*(10A(s n/10);n=sqrt(db)*randn(size(t);nsfm=n+sfm;%設(shè)定信躁比 （小信噪比 ）% 設(shè)定信躁比 （大信噪比 ）% 設(shè)定信躁比 （無信噪比 ） %計(jì)算對應(yīng)的高斯白躁聲的方差 %生成高斯白躁聲 %生成含高斯白躁聲的已調(diào)信號（信號通 % 過信道傳輸）am=5;fm=5; mt=am*cos(2*pi*fm*t);fc=50;ct=cos(2*pi*fc*t);kf=10;

24、int_mt(1)=0;for i=1:length(t)-1int_mt(i+1)=int_mt(i)+mt(i)*dt; end sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt); %* % . * * * * * * * * * * * * * * * .17 / 25第 24 頁 共 25 頁%*% * * * * * * * * * * * * * * *FM解調(diào)*%接受信號通過微分器處理for i=1:length(t)-1diff_nsfm(i)=(nsfm(i+1)-nsfm(i)./dt;enddiff_nsfmn = abs(hilbert(diff

25、_nsfm);%hilbert 變換，求絕對值得到瞬時(shí)幅度（包絡(luò)檢波）zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn-zero;%*% * * * * * * * * * * * * * * * fs=1/ts;%* ts=0.001; fs=1/ts; df=0.25;時(shí)域到頻域轉(zhuǎn)換 *%*對調(diào)制信號m(t)求傅里葉變換* m=am*cos(2*pi*fm*t);% 抽樣間隔%抽樣頻率% 所需的頻率分辨率， 用在求傅里葉變換%時(shí)，它表示 FFT 的最小頻率間隔%原調(diào)信號if nargin=2 n1=0;elsen1

26、=fs/df;endn2=length(m);n=2A(max (n extpow2( n1),n extpow2( n2);M=fft(m,n);m=m,zeros(1,n-n2);df1=fs/n;%以上程序是對調(diào)制后的信號 u 求傅里變換M=M/fs;%縮放，便于在頻鋪圖上整體觀察f=0:df1:df1*(length(m)-1)-fs/2;%時(shí)間向量對應(yīng)的頻率向量%*對已調(diào)信號 u 求傅里變換 *fs=1/ts;if nargin=2 n1=0;elsen1=fs/df;end n2=length(sfm);n=2A(max(nextpow2(n1),nextpow2(n2);U=ff

27、t(sfm,n);u=sfm,zeros(1,n-n2);% 以上是對已調(diào)信號 u 求傅里變換% 縮放df1=fs/n;U=U/fs;%* %*% . * * * * * * * * * * * * * * * .%*顯示程序*disp（ 按任意鍵可以看到原調(diào)制信號、載波信號和已調(diào)信號的曲線）pause*figure(1)*figure(1)subplot(3,1,1);plot(t,mt);xlabel( 時(shí)間 t);title( 調(diào)制信號的時(shí)域圖 );subplot(3,1,2);plot(t,ct);xlabel( 時(shí)間 t);title( 載波的時(shí)域圖 );subplot(3,1,3)

28、;plot(t,sfm);xlabel( 時(shí)間 t);title( 已調(diào)信號的時(shí)域圖 );%*% 繪制調(diào)制信號的時(shí)域圖% 繪制載波的時(shí)域圖% 繪制已調(diào)信號的時(shí)域圖disp（ 按任意鍵可以看到原調(diào)制信號和已調(diào)信號在頻域內(nèi)的圖形） pause%* *figure(2) *figure(2) subplot(2,1,1) plot(f,abs(fftshift(M) %fftshift: 將 FFT 中的 DC 分量移到頻譜中心 xlabel( 頻率 f)title( 原調(diào)制信號的頻譜圖 )subplot(2,1,2) plot(f,abs(fftshift(U) xlabel( 頻率 f) tit

29、le( 已調(diào)信號的頻譜圖 )%*disp( 按任意鍵可以看到原調(diào)制信號、無噪聲條件下已調(diào)信號和解調(diào)信號的曲線) pause %* *figure(3) *figure(3)%繪制調(diào)制信號的時(shí)域圖%繪制已調(diào)信號的時(shí)域圖%接受信號通過微分器處理subplot(3,1,1);plot(t,mt);xlabel( 時(shí)間 t);title( 調(diào)制信號的時(shí)域圖 );subplot(3,1,2);plot(t,sfm);xlabel( 時(shí)間 t);title( 無噪聲條件下已調(diào)信號的時(shí)域圖 ); nsfm=sfm;for i=1:length(t)-1diff_nsfm(i)=(nsfm(i+1)-nsfm

30、(i)./dt; enddiff_nsfmn = abs(hilbert(diff_nsfm);%hilbert 變換，求絕對值得到瞬時(shí)幅度（包絡(luò)檢波）zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;%繪制無噪聲條件下解調(diào)信號的時(shí)域圖diff_nsfmn1=diff_nsfmn-zero;subplot(3,1,3);plot(1:length(diff_nsfmn1)./1000,diff_nsfmn1./400,r); xlabel( 時(shí)間 t);title( 無噪聲條件下解調(diào)信號的時(shí)域圖 );%* disp（按任意鍵可以看到原調(diào)制信號、小信噪比高斯白噪聲條

31、件下已調(diào)信號和解調(diào)信號已調(diào) 信號的曲線 ）pause*figure(4)*figure(4)subplot(3,1,1);plot(t,mt);xlabel( 時(shí)間 t);title( 調(diào)制信號的時(shí)域圖 );db1=amA2/(2*(10A(s n1/10);%繪制調(diào)制信號的時(shí)域圖n1=sqrt(db1)*randn(size(t);nsfm1=n1+sfm;%計(jì)算對應(yīng)的小信噪比高斯白躁聲的方差%生成高斯白躁聲%生成含高斯白躁聲的已調(diào)信號（信號通%過信道傳輸)for i=1:length(t)-1%接受信號通過微分器處理diff_nsfm1(i)=(nsfm1(i+1)-nsfm1(i)./d

32、t;enddiff_nsfmn1 = abs(hilbert(diff_nsfm1); %hilbert 變換，求絕對值得到瞬時(shí)幅度(包絡(luò)檢波) zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn1-zero;subplot(3,1,2);plot(1:length(diff_nsfm),diff_nsfm);%繪制含小信噪比高斯白噪聲已調(diào)信號的時(shí)域圖xlabel( 時(shí)間 t);title( 含小信噪比高斯白噪聲已調(diào)信號的時(shí)域圖 );subplot(3,1,3);%繪制含小信噪比高斯白噪聲解調(diào)信號的時(shí)域圖plot(1:l

33、ength(diff_nsfmn1)./1000,diff_nsfmn1./400,r);xlabel( 時(shí)間 t);title( 含小信噪比高斯白噪聲解調(diào)信號的時(shí)域圖 );%*disp(按任意鍵可以看到原調(diào)制信號、大信噪比高斯白噪聲條件下已調(diào)信號和解調(diào)信號已調(diào) 信號的曲線 )pause%* *figure(5) * figure(5) subplot(3,1,1);plot(t,mt); xlabel( 時(shí)間 t);title( 調(diào)制信號的時(shí)域圖 );db1=amA2/(2*(10A(s n2/10); n1=sqrt(db1)*randn(size(t); nsfm1=n1+sfm;過信道傳輸) for i=1:length(t)-1diff_nsfm1(i)=(nsfm1(i+1)-nsfm1(i)./dt; end