【正文】 n2=length(m)。 %************************************************************ %*.*.*.*.*.*.*.*.*.*.*.*.*. %***********************顯示程序 ***************************** disp(39。載波的時域圖 39。原調(diào)制信號的頻譜圖 39。 %繪制已調(diào)信號的時域圖 xlabel(39。)。 echo off end echo on diff_nsfmn1=abs(hilbert(diff_nsfm1))。)。 %生成高斯白噪聲 nsfm1=n1+sfm。)。)。調(diào)制信號的時域圖 39。r39。 db1=am^2/(2*(10^(sn1/10)))。r39。 title(39。)。 %繪制載波的時域圖 xlabel(39。 U=fft(sfm,n)。 %抽樣頻率 df=。 for i=1:length(t)1 int_mt(i+1)=int_mt(i)+mt(i)*dt。 n(t)是均值為零，單邊功率譜密度為 n0的高斯白噪聲，經(jīng)過帶通濾波器后變?yōu)檎瓗Ц咚乖肼?ni(t)。 %調(diào)制，產(chǎn)生已調(diào)信號 圖 3 FM調(diào)制 FM解調(diào)模型的建立 調(diào)制信號的解調(diào)分為相干解調(diào)與非相干解調(diào)兩種。在本仿真的過程中我們選擇用調(diào)頻調(diào)制方式進行調(diào)制。另一方面，通過調(diào)制可以提高信號通過信道傳輸時的抗干擾能力，同時，它還和傳輸效率有關(guān)。 %設(shè)定載波頻率 ← 可更改 ct=cos(2*pi*fc*t)。 在本仿真中，選擇了 10dB和 30dB兩種不同信噪比以示區(qū)別，其時域圖如圖 6和圖 7 0 0 . 5 1 1 . 5505時間 t調(diào)制信號的時域圖0 0 . 5 1 1 . 5505時間 t無噪聲條件下已調(diào)信號的時域圖0 0 . 5 1 1 . 5 1 0010時間 t無噪聲條件下解調(diào)信號的時域圖姓名：黃婷 班級： 13 級電子信息工程 1 班 學(xué)號： 1308421103 圖 6 含小信噪比高斯白噪聲已調(diào)信號的時域圖 圖 7 含大信噪比高斯白噪聲已調(diào)信號的時域圖 2。 %設(shè)定載波頻率←可更改 ct=cos(2*pi*fc*t)。 diff_nsfmn1=diff_nsfmnzero。 else n1=fs/df。)。)。 %繪制調(diào)制信號的時域圖 xlabel(39。 姓名：黃婷 班級： 13 級電子信息工程 1 班 學(xué)號： 1308421103 diff_nsfmn1=diff_nsfmnzero。) title(39。)。時間 t39。)。) title(39。 %hilbert變換，求絕對值得到瞬時幅度（包絡(luò)檢波） zero=(max(diff_nsfmn)min(diff_nsfm))/2。)。 subplot(3,1,2)。) %********************************************* disp(39。無噪聲條件下已調(diào)信號的時域圖 39。) title(39。 plot(t,sfm)。plot(t,mt)。 m=[m,zeros(1,nn2)]。 %計算對應(yīng)的高斯白噪聲的方差 n=sqrt(db)*randn(size(t))。 初始化程序 FM 調(diào)制程序 添加高斯白噪聲程序 FM 解調(diào)程序 繪制頻譜圖程序 顯示程序 姓名：黃婷 班級： 13 級電子信息工程 1 班 學(xué)號： 1308421103 程序： %FM調(diào)制解調(diào)系統(tǒng) .m %頻率調(diào)制與解調(diào)的 %Matlab演示源程序 %可以任意改變原調(diào)制信號函數(shù) m(t) %*.*.*.*.*.*.*.*.*. %********************初始化 ************************ echo off close all clear all clc %**************************************************** %*.*.*.*.*.*.*.*.*.*.*.*.*. %********************FM調(diào)制 **************************** dt=。微分器輸出為 Sd(t) = dSi(t)dt = dSFM(t)dt = ?[wc + Kfm(t)sin(wc (t) +Kf ∫ m(τ)dτt?∞ 包絡(luò)檢波的作用是從輸出信號的幅度變化中檢出調(diào)制信號。 調(diào)制過程分析 在調(diào)制時，調(diào)制信號的頻率去控制載波的頻率的變化，載波的瞬時頻偏隨調(diào)制信號 m(t)成正比例變化，即 dφ(t)dt = Kfm(t) 式中， kf為調(diào)頻靈敏度 (rad (s ? v))? 這時相位偏移為 φ(t) = Kfm(τ)dτ 則可得到調(diào)頻信號為