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

1.1FM調(diào)制模型的建立1.2調(diào)制過程分析制信號m(t)成正比例變化，即=Kfm(t)則可得到調(diào)頻信號為%設(shè)定時間步長mt=am*cos(2*pi*fm*t);ct=cos(2*pi*fc*t);int_mt(i+1)=int_mt(i)+mt(i)*dt;sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt);%產(chǎn)生時間向量%設(shè)定調(diào)制信號幅度←可更改%設(shè)定調(diào)制信號頻率←可更改%生成調(diào)制信號%設(shè)定載波頻率←可更改%設(shè)定調(diào)頻指數(shù)%調(diào)制，產(chǎn)生已調(diào)信號01000.511.5010-100.5001000.501.3FM解調(diào)模型的建立1.4解調(diào)過程分析dHzK稱為鑒頻靈敏度(VdHzt(包絡(luò)檢波)5000.505000.505000.501.5高斯白噪聲信道特性設(shè)正弦波通過加性高斯白噪聲信道后的信號為A22N=2N信噪比S滿足公式N則可得到公式2=B5000.505000.505000.50圖6無噪聲條件下已調(diào)信號的時域圖5000.50聲已調(diào)信號的時域圖0050010001500聲解調(diào)信號的時域圖000.50圖7含小信噪比高斯白噪聲已調(diào)信號的時域圖5000.511.5聲已調(diào)信號的時域圖0050010001500聲解調(diào)信號的時域圖5000.50圖8含大信噪比高斯白噪聲已調(diào)信號的時域圖初始化程序初始化程序FM調(diào)制程序添加高斯白噪聲程序FM解調(diào)程序繪制頻譜圖程序3.1MATLAB源代碼FM統(tǒng).m%*****************初始化******************offlearall%*****************************************%****************FM調(diào)制*******************%設(shè)定時間步長生時間向量%設(shè)定調(diào)制信號幅度mt=am*cos(2*pi*fm*t);ct=cos(2*pi*fc*t);intmti+1)=int_mt(i)+mt(i)*dt;sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt);%設(shè)定調(diào)制信號頻率%生成調(diào)制信號%設(shè)定載波頻率載波定調(diào)頻指數(shù)%調(diào)制，產(chǎn)生已調(diào)信號調(diào)制信號%*****************************************%*************添加高斯白噪聲**************%設(shè)定信躁比(小信噪比)%設(shè)定信躁比(大信噪比)%設(shè)定信躁比(無信噪比)%計算對應(yīng)的高斯白躁聲的方差%生成高斯白躁聲%生成含高斯白躁聲的已調(diào)信號(信號通%*****************************************%****************FM解調(diào)*******************forilengtht%接受信號通過微分器處理diff_nsfmn=abs(hilbert(diff_nsfm));%hilbert變換，求絕對值得到瞬時幅度(包絡(luò)檢波)%*****************************************%**************時域到頻域轉(zhuǎn)換**************%*****對調(diào)制信號m(t)求傅里葉變換*****m=am*cos(2*pi*fm*t);抽樣間隔信號M=fft(m,n);換M=M/fs;%以上程序是對調(diào)制后的信號u求傅里變%縮放，便于在頻鋪圖上整體觀察%時間向量對應(yīng)的頻率向量%************對已調(diào)信號u求傅里變換**********%******************************************%*****************************************%***************顯示程序******************disp('按任意鍵可以看到原調(diào)制信號、載波信號和已調(diào)信號的曲線')%**************figure(1)******************title('調(diào)制信號的時域圖');title('載波的時域圖');title('已調(diào)信號的時域圖');%繪制調(diào)制信號的時域圖%繪制載波的時域圖%繪制已調(diào)信號的時域圖%******************************************disp('按任意鍵可以看到原調(diào)制信號和已調(diào)信號在頻域內(nèi)的圖形')%************figure(2)*********************title('原調(diào)制信號的頻譜圖')title('已調(diào)信號的頻譜圖')%fftshift:將FFT中的DC分量移到頻譜中心%******************************************disp('按任意鍵可以看到原調(diào)制信號、無噪聲條件下已調(diào)信號和解調(diào)信號的曲線')%**************figure(3)******************title('調(diào)制信號的時域圖');title('無噪聲條件下已調(diào)信號的時域圖');%繪制調(diào)制信號的時域圖%繪制已調(diào)信號的時域圖%接受信號通過微分器處理diff_nsfmn=abs(hilbert(diff_nsfm));%hilbert變換，求絕對值得到瞬時幅度(包絡(luò)檢波)subplot噪聲條件下解調(diào)信號的時域圖title('無噪聲條件下解調(diào)信號的時域圖');%*****************************************disp('按任意鍵可以看到原調(diào)制信號、小信噪比高斯白噪聲條件下已調(diào)信號和解調(diào)信號已調(diào)信號的曲線')%**************figure(4)******************title('調(diào)制信號的時域圖');%繪制調(diào)制信號的時域圖%計算對應(yīng)的小信噪比高斯白躁聲的方差%生成高斯白躁聲%生成含高斯白躁聲的已調(diào)信號(信號通forilengtht)-1%接受信號通過微分器處理diff_nsfmn1=abs(hilbert(diff_nsfm1));%hilbert變換，求絕對值得到瞬時幅度(包絡(luò)檢波)plotlengthdiffnsfmdiffnsfm噪比高斯白噪聲已調(diào)信號的時域圖title('含小信噪比高斯白噪聲已調(diào)信號的時域圖');subplot繪制含小信噪比高斯白噪聲解調(diào)信號的時域圖title('含小信噪比高斯白噪聲解調(diào)信號的時域圖');%*****************************************disp('按任意鍵可以看到原調(diào)制信號、大信噪比高斯白噪聲條件下已調(diào)信號和解調(diào)信號已調(diào)信號的曲線')%**************figure(5)******************title('調(diào)制信號的時域圖');輸)%繪制調(diào)制信號的時域圖%計算對應(yīng)的大信噪比高斯白躁聲的方差%生成高斯白躁聲%生成含高斯白躁聲的已調(diào)信號(信號通%接受信號通過微分器處理%繪制含大信噪比高斯白噪聲已調(diào)信號title('含大信噪比高斯白噪聲已調(diào)信號的時域