【正文】 %調制，產生已調信號sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt)。******** %設定信躁比(小信噪比)sn2=30。 %設定信躁比(無信噪比)db=am^2/(2*(10^(sn/10)))。 %生成高斯白躁聲nsfm=n+sfm。********enddiff_nsfmn = abs(hilbert(diff_nsfm))。diff_nsfmn1=diff_nsfmnzero。******** %抽樣間隔fs=1/ts。 %所需的頻率分辨率，用在求傅里葉變換%時，它表示FFT的最小頻率間隔%*****對調制信號m(t)求傅里葉變換*****m=am*cos(2*pi*fm*t)。if nargin==2 n1=0。endn2=length(m)。M=fft(m,n)。df1=fs/n。 %縮放，便于在頻鋪圖上整體觀察f=[0:df1:df1*(length(m)1)]fs/2。if nargin==2 n1=0。endn2=length(sfm)。U=fft(sfm,n)。df1=fs/n。 %縮放%******************************************%*****************************************%*******%***************顯示程序******************disp(39。)pause%**************figure(1)******************figure(1)subplot(3,1,1)。 %繪制調制信號的時域圖xlabel(39。)。調制信號的時域圖39。subplot(3,1,2)。 %繪制載波的時域圖xlabel(39。)。載波的時域圖39。subplot(3,1,3)。 %繪制已調信號的時域圖xlabel(39。)。已調信號的時域圖39。%******************************************disp(39。)pause%************figure(2)*********************figure(2)subplot(2,1,1)plot(f,abs(fftshift(M))) %fftshift:將FFT中的DC分量移到頻譜中心xlabel(39。)title(39。)subplot(2,1,2)plot(f,abs(fftshift(U)))xlabel(39。)title(39。)%******************************************disp(39。)pause%**************figure(3)******************figure(3)subplot(3,1,1)。 %繪制調制信號的時域圖xlabel(39。)。調制信號的時域圖39。subplot(3,1,2)。 %繪制已調信號的時域圖xlabel(39。)。無噪聲條件下已調信號的時域圖39。nsfm=sfm。enddiff_nsfmn = abs(hilbert(diff_nsfm))。diff_nsfmn1=diff_nsfmnzero。 %繪制無噪聲條件下解調信號的時域圖plot((1:length(diff_nsfmn1))./1000,diff_nsfmn1./400,39。)。時間t39。 title(39。)。按任意鍵可以看到原調制信號、小信噪比高斯白噪聲條件下已調信號和解調信號已調信號的曲線39。plot(t,mt)。時間t39。title(39。)。 %計算對應的小信噪比高斯白躁聲的方差n1=sqrt(db1)*randn(size(t))。 %生成含高斯白躁聲的已調信號（信號通%過信道傳輸）for i=1:length(t)1 %接受信號通過微分器處理 diff_nsfm1(i)=(nsfm1(i+1)nsfm1(i))./dt。 %hilbert變換，求絕對值得到瞬時幅度（包絡檢波）zero=(max(diff_nsfmn)min(diff_nsfmn))/2。subplot(3,1,2)。 %繪制含小信噪比高斯白噪聲已調信號的時域圖xlabel(39。)。含小信噪比高斯白噪聲已調信號的時域圖39。subplot(3,1,3)。r39。xlabel(39。)。含小信噪比高斯白噪聲解調信號的時域圖39。%*****************************************disp(39。)pause%**************figure(5)******************figure(5)subplot(3,1,1)。 %繪制調制信號的時域圖xlabel(39。)。調制信號的時域圖39。db1=am^2/(2*(10^(sn2/10)))。 %生成高斯白躁聲nsfm1=n1+sfm。enddiff_nsfmn1 = abs(hilbert(diff_nsfm1))。diff_nsfmn1=diff_nsfmn1zero。plot(1:length(diff_nsfm1),diff_nsfm1)。時間t39。title(39。)。 %繪制含大信噪比高斯白噪聲解調信號%的時域圖plot((1:length(diff_nsfmn1))./1000,diff_nsfmn1./400,39。)。時間t39。 title(39。