fm系統(tǒng)的設(shè)計(jì)(全)(編輯修改稿)
2025-07-26 19:19 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 指數(shù)，可使調(diào)頻系統(tǒng)的抗噪聲性能迅速改善。參考文獻(xiàn)【1】樊昌信，：國(guó)防工業(yè)出版社【2】【3】附錄1 Matlab源代碼%%頻率調(diào)制與解調(diào)的Matlab演示源程序%可以任意改原調(diào)制信號(hào)函數(shù)m(t)%*****************初始化******************echo off close allclear allclc%*****************************************%***************%****************FM調(diào)制*******************dt=。 %設(shè)定時(shí)間步長(zhǎng)t=0:dt:。 %產(chǎn)生時(shí)間向量am=5。 %設(shè)定調(diào)制信號(hào)幅度fm=5。 %設(shè)定調(diào)制信號(hào)頻率mt=am*cos(2*pi*fm*t)。 %生成調(diào)制信號(hà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)1 int_mt(i+1)=int_mt(i)+mt(i)*dt。 %求信號(hào)m(t)的積分end %調(diào)制，產(chǎn)生已調(diào)信號(hào)sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt)。 %調(diào)制信號(hào)%*****************************************%***************%*************添加高斯白噪聲**************sn1=10。 %設(shè)定信躁比(小信噪比)sn2=30。 %設(shè)定信躁比(大信噪比)sn=0。 %設(shè)定信躁比(無信噪比)db=am^2/(2*(10^(sn/10)))。 %計(jì)算對(duì)應(yīng)的高斯白躁聲的方差n=sqrt(db)*randn(size(t))。 %生成高斯白躁聲nsfm=n+sfm。 %生成含高斯白躁聲的已調(diào)信號(hào)（信號(hào)通%過信道傳輸）%*****************************************%***************%****************FM解調(diào)*******************for i=1:length(t)1 %接受信號(hào)通過微分器處理 diff_nsfm(i)=(nsfm(i+1)nsfm(i))./dt。enddiff_nsfmn = abs(hilbert(diff_nsfm))。 %hilbert變換，求絕對(duì)值得到瞬時(shí)幅度（包絡(luò)檢波）zero=(max(diff_nsfmn)min(diff_nsfmn))/2。diff_nsfmn1=diff_nsfmnzero。%*****************************************%***************%**************時(shí)域到頻域轉(zhuǎn)換**************ts=。 %抽樣間隔fs=1/ts。 %抽樣頻率df=。 %所需的頻率分辨率，用在求傅里葉變換%時(shí)，它表示FFT的最小頻率間隔%*****對(duì)調(diào)制信號(hào)m(t)求傅里葉變換*****m=am*cos(2*pi*fm*t)。 %原調(diào)信號(hào)fs=1/ts。if nargin==2 n1=0。else n1=fs/df。endn2=leng
研究報(bào)告相關(guān)推薦
鄂ICP備17016276號(hào)-1