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

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

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

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

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

5、nication system is widely used. FM is widelyused 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 achieve FM modulation and demodulation proc

6、ess, 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; finally draw FM baseband signal through th

7、e 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 signal demodulation. In the course of desi

8、gn, 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, to achieve the realization of the FM s

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

10、LAB 程序?qū)崿F(xiàn)（即只能用代碼的形式， 不能用 SIMULINK實現(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è)計通信的目的是傳輸信息。 通信系統(tǒng)的作用就是將信息從信息源發(fā)送到一個或多個目的地。對于任何一個通信系統(tǒng)，均可視為由發(fā)送端、信道和接收端三大部分組成（如圖 1 所示）。發(fā)送端信道接收端信息源發(fā)送設(shè)備信 道接受設(shè)備信息源噪聲源圖 1 通信系統(tǒng)一般模型信息源（簡稱信源）的作用是把各種信息轉(zhuǎn)換成原始信號。根據(jù)消息的種類不同信源分為模擬信源和數(shù)字信源。發(fā)送設(shè)備的作用產(chǎn)

11、生適合傳輸?shù)男盘枺词拱l(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ù)倪^程總會受到噪聲的干擾， 通信系統(tǒng)中沒有傳輸信號時也有噪聲，噪聲永遠(yuǎn)存在于通信系統(tǒng)中。 由于這樣的噪聲是疊加在信號

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

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

14、可得到調(diào)頻信號為sFM(t)AcosctKfm( )d調(diào)制信號產(chǎn)生的M 文件：dt=0.001;%設(shè)定時間步長t=0:dt:1.5;am=15;%產(chǎn)生時間向量%設(shè)定調(diào)制信號幅度可更改fm=15;mt=am*cos(2*pi*fm*t);%設(shè)定調(diào)制信號頻率可更改%生成調(diào)制信號fc=50;%設(shè)定載波頻率可更改ct=cos(2*pi*fc*t);%生成載波kf=10;%設(shè)定調(diào)頻指數(shù)int_mt(1)=0;%對 mt 進(jìn)行積分for i=1:length(t)-1int_mt(i+1)=int_mt(i)+mt(i)*dt;endsfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt)

15、;%調(diào)制，產(chǎn)生已調(diào)信號10調(diào)制信號的時域圖0-1000.511.5時 間 t載波的時域圖10-100.511.5時 間 t已調(diào)信號的時域圖100-1000.511.5時 間 t圖3 FM調(diào)制2.3FM 解調(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)。圖 4FM 解調(diào)模型非相干解調(diào)器由限幅器、 鑒頻器和低通濾波器等組成， 其方框圖如圖 5 所示。限幅器輸入為已調(diào)頻信號

16、和噪聲， 限幅器是為了消除接收信號在幅度上可能出現(xiàn)的畸變；帶通濾波器的作用是用來限制帶外噪聲，使調(diào)頻信號順利通過。 鑒頻器中的微分器把調(diào)頻信號變成調(diào)幅調(diào)頻波， 然后由包絡(luò)檢波器檢出包絡(luò)， 最后通過低通濾波器取出調(diào)制信號。2.4解調(diào)過程分析設(shè)輸入調(diào)頻信號為tSt(t)SFM(t)Acos(ct Kfm( )d )微分器的作用是把調(diào)頻信號變成調(diào)幅調(diào)頻波。微分器輸出為ddS(t)dS(t)S (t)iFMdtdtKfm(t) sin(tcctKf m( )d )包絡(luò)檢波的作用是從輸出信號的幅度變化中檢出調(diào)制信號。包絡(luò)檢波器輸出為So(t)KdcKmf(t)KdcKdKmf (t)Kd稱為鑒頻靈敏度（

17、VHz），是已調(diào)信號單位頻偏對應(yīng)的調(diào)制信號的幅度，經(jīng)低通濾波器后加隔直流電容，隔除無用的直流，得m o (t )K d K f m (t )微分器通過程序?qū)崿F(xiàn)，代碼如下：for i=1:length(t)-1diff_nsfm(i)=(nsfm(i+1)-nsfm(i)./dt;%接受信號通過微分器處理enddiff_nsfmn = abs(hilbert(diff_nsfm);%hilbert變換，求絕對值得到瞬時幅度（包絡(luò)檢波）通過 M 文件繪制出兩種不同信噪比解調(diào)的輸出波形如下：調(diào)制信號的時域圖50-500.511.5時 間 t無噪聲條件下已調(diào)信號的時域圖50-500.511.5時 間

18、t無噪聲條件下解調(diào)信號的時域圖50-500.511.5時 間 t圖5 FM解調(diào)2.5高斯白噪聲信道特性設(shè)正弦波通過加性高斯白噪聲信道后的信號為r (t )A cos( 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)故其有用信號功率為S噪聲功率為A222N信噪比 S N 滿足公式B 10log10(SN)則可得到公式2 A 2B210 10我們可以通過這個公式方便的

19、設(shè)置高斯白噪聲的方差。在本仿真過程中， 我們選擇了 10db 和 30db 兩種不同信噪比以示區(qū)別， 其時域圖如圖 7 和圖 8。調(diào)制信號的時域圖50-500.511.5時 間 t無噪聲條件下已調(diào)信號的時域圖50-500.511.5時 間 t無噪聲條件下解調(diào)信號的時域圖50-500.511.5時 間 t圖 6無噪聲條件下已調(diào)信號的時域圖調(diào)制信號的時域圖50-500.511.5時 間 t含小信噪比高斯白噪聲已調(diào)信號的時域圖50000-5000050010001500時 間 t含小信噪比高斯白噪聲解調(diào)信號的時域圖200-2000.511.5時 間 t圖 7 含小信噪比高斯白噪聲已調(diào)信號的時域圖調(diào)制

20、信號的時域圖50-500.511.5時 間 t含大信噪比高斯白噪聲已調(diào)信號的時域圖50000-5000050010001500時 間 t含大信噪比高斯白噪聲解調(diào)信號的時域圖50-500.511.5時 間 t圖 8含大信噪比高斯白噪聲已調(diào)信號的時域圖2.6調(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)的抗噪聲性能，其分析模型如圖 9 所示。圖 9調(diào)頻系統(tǒng)抗噪聲性能分析模型圖中帶通濾波器的作

21、用是抑制信號帶寬以外的噪聲。n(t ) 是均值為零，單邊功率譜密度為 n0 的高斯白噪聲，經(jīng)過帶通濾波器后變?yōu)檎瓗Ц咚乖肼昻i(t )。限幅器是為了消除接收信號在幅度上可能出現(xiàn)的畸變。設(shè)調(diào)頻信號為tSFM (t)Acos( ctKfm( )d )故其輸入功率為S iA22輸入噪聲功率為NinoBFM因此輸入信噪比為SiA2Ni2BFM在大信噪比條件下， 信號和噪聲的相互作用可以忽略， 這時可以把信號和噪聲分開來算，這里，我們可以得到解調(diào)器的輸出信噪比o3A2K2 m2(t )SfNo8 2 nof m3上式中， A 為載波的振幅， Kf 為調(diào)頻器靈敏度，fm 為調(diào)制信號 m(t) 的最高頻率，

22、 no 為噪聲單邊功率譜密度。我們?nèi)缛艨紤] m(t ) 為單一頻率余弦波時的情況， 可得到解調(diào)器的制度增益為SA22o3G FMN om2Si2f n o f mN i考慮在寬帶調(diào)頻時，信號帶寬為B2(m 1) fm2( ffm)FMf則可以得到GFM3m2f ( mf1)可以看出，大信噪比時寬帶調(diào)頻系統(tǒng)的信噪比增益是很高的， 它與調(diào)頻指數(shù)的立方成正比。可見，加大調(diào)頻指數(shù) mf ，可使調(diào)頻系統(tǒng)的抗噪聲性能迅速改善。三仿真實現(xiàn)初始化程序FM調(diào)制程序添加高斯白噪聲程序FM解調(diào)程序繪制頻譜圖程序顯示程序圖 10程序流程圖3.1MATLAB源代碼%FM 調(diào)制解調(diào)系統(tǒng).m%頻率調(diào)制與解調(diào)的Matlab

23、演示源程序%可以任意改原調(diào)制信號函數(shù)m(t)%信息工程陳麗丹07323202%* *%*初始化 *echooffcloseallclearallclc%*%* *%*FM調(diào)制 *dt=0.001;t=0:dt:1.5;%設(shè)定時間步長%產(chǎn)生時間向量am=5;%設(shè)定調(diào)制信號幅度fm=5;%設(shè)定調(diào)制信號頻率mt=am*cos(2*pi*fm*t);% 生成調(diào)制信號fc=50;%設(shè)定載波頻率ct=cos(2*pi*fc*t);%生成載波kf=10;% 設(shè)定調(diào)頻指數(shù)int_mt(1)=0;for i=1:length(t)-1int_mt(i+1)=int_mt(i)+mt(i)*dt;% 求信號 m(

24、t) 的積分end% 調(diào)制，產(chǎn)生已調(diào)信號sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt);% 調(diào)制信號%*%* *%*添加高斯白噪聲 *sn1=10;%設(shè)定信躁比 (小信噪比 )sn2=30;% 設(shè)定信躁比 (大信噪比 )sn=0;% 設(shè)定信躁比 (無信噪比 )db=am2/(2*(10(sn/10);%計算對應(yīng)的高斯白躁聲的方差n=sqrt(db)*randn(size(t);%生成高斯白躁聲nsfm=n+sfm;% 生成含高斯白躁聲的已調(diào)信號（信號通%過信道傳輸）%*%* *%*FM解調(diào) *for i=1:length(t)-1%接受信號通過微分器處理diff_ns

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

26、號fs=1/ts;if nargin=2n1=0;elsen1=fs/df;endn2=length(m);n=2(max(nextpow2(n1),nextpow2(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;%縮放，便于在頻鋪圖上整體觀察%時間向量對應(yīng)的頻率向量%*對已調(diào)信號u 求傅里變換*fs=1/ts;if nargin=2n1=0;elsen1=fs/df;endn2=length(sfm);n=2(max(nextpow2(n

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

28、3,1,3);plot(t,sfm);%繪制已調(diào)信號的時域圖xlabel( 時間 t);title( 已調(diào)信號的時域圖);%*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)title( 已調(diào)信號的頻譜圖)%*disp( 按任意鍵可以看到原

29、調(diào)制信號、無噪聲條件下已調(diào)信號和解調(diào)信號的曲線)pause%*figure(3)*figure(3)subplot(3,1,1);plot(t,mt);%繪制調(diào)制信號的時域圖xlabel( 時間t);title( 調(diào)制信號的時域圖);subplot(3,1,2);plot(t,sfm);%繪制已調(diào)信號的時域圖xlabel( 時間t);title( 無噪聲條件下已調(diào)信號的時域圖);nsfm=sfm;for i=1:length(t)-1% 接受信號通過微分器處理diff_nsfm(i)=(nsfm(i+1)-nsfm(i)./dt;enddiff_nsfmn = abs(hilbert(diff

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

31、re(4)subplot(3,1,1);plot(t,mt);%繪制調(diào)制信號的時域圖xlabel( 時間t);title( 調(diào)制信號的時域圖);db1=am2/(2*(10(sn1/10);n1=sqrt(db1)*randn(size(t);nsfm1=n1+sfm;%計算對應(yīng)的小信噪比高斯白躁聲的方差%生成高斯白躁聲% 生成含高斯白躁聲的已調(diào)信號（信號通%過信道傳輸）for i=1:length(t)-1diff_nsfm1(i)=(nsfm1(i+1)-nsfm1(i)./dt;%接受信號通過微分器處理enddiff_nsfmn1 = abs(hilbert(diff_nsfm1);%h

32、ilbert變換，求絕對值得到瞬時幅度（包絡(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)信號的時域圖xlabel( 時間t);title( 含小信噪比高斯白噪聲已調(diào)信號的時域圖);subplot(3,1,3);%繪制含小信噪比高斯白噪聲解調(diào)信號的時域圖plot(1:length(diff_nsfmn1)./1000,diff_nsfmn1./400,r);xlabel( 時間 t);title( 含小信噪比高斯白噪聲解調(diào)信號的時域圖);%*disp( 按任意鍵可以看到原調(diào)制信號、大信噪比高斯白噪聲條件下已調(diào)信號和解調(diào)信號已調(diào)信號的曲線)pause%*figure(5)*figure(5)subplot(3,1,1);plot(t,mt);%繪制調(diào)制信號的時域圖xlabel( 時間t);title( 調(diào)制信號的時域圖);db1=am2/(2*(10(sn2/10);n1=sqrt(db1)*randn(size(t);nsfm1=n1+sfm;%計算