fm調(diào)制解調(diào)系統(tǒng)設計與仿真(編輯修改稿)
2024-09-01 01:26 本頁面
【文章內(nèi)容簡介】 83。***********%****************FM調(diào)制*******************dt=。 %設定時間步長t=0:dt:。 %產(chǎn)生時間向量am=5。 %設定調(diào)制信號幅度fm=5。 %設定調(diào)制信號頻率mt=am*cos(2*pi*fm*t)。 %生成調(diào)制信號fc=50。 %設定載波頻率ct=cos(2*pi*fc*t)。 %生成載波kf=10。 %設定調(diào)頻指數(shù)int_mt(1)=0。for i=1:length(t)1 int_mt(i+1)=int_mt(i)+mt(i)*dt。 %求信號m(t)的積分end %調(diào)制，產(chǎn)生已調(diào)信號sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt)。 %調(diào)制信號%*****************************************%***************%*************添加高斯白噪聲**************sn1=10。 %設定信躁比(小信噪比)sn2=30。 %設定信躁比(大信噪比)sn=0。 %設定信躁比(無信噪比)db=am^2/(2*(10^(sn/10)))。 %計算對應的高斯白躁聲的方差n=sqrt(db)*randn(size(t))。 %生成高斯白躁聲nsfm=n+sfm。 %生成含高斯白躁聲的已調(diào)信號（信號通%過信道傳輸）%*****************************************%***************%****************FM解調(diào)*******************for i=1:length(t)1 %接受信號通過微分器處理 diff_nsfm(i)=(nsfm(i+1)nsfm(i))./dt。enddiff_nsfmn = abs(hilbert(diff_nsfm))。 %hilbert變換，求絕對值得到瞬時幅度（包絡檢波）zero=(max(diff_nsfmn)min(diff_nsfmn))/2。diff_nsfmn1=diff_nsfmnzero。%*****************************************%***************%**************時域到頻域轉(zhuǎn)換**************ts=。 %抽樣間隔fs=1/ts。 %抽樣頻率df=。 %所需的頻率分辨率，用在求傅里葉變換%時，它表示FFT的最小頻率間隔%*****對調(diào)制信號m(t)求傅里葉變換*****m=am*cos(2*pi*fm*t)。 %原調(diào)信號fs=1/ts。if nargin==2 n1=0。else n1=fs/df。endn2=length(m)。n=2^(max(nextpow2(n1),nextpow2(n2)))。M=fft(m,n)。m=[m,zeros(1,nn2)]。df1=fs/n。 %以上程序是對調(diào)制后的信號u求傅里變換M=M/fs。 %縮放，便于在頻鋪圖上整體觀察f=[0:df1:df1*(length(m)1)]fs/2。 %時間向量對應的頻率向量%************對已調(diào)信號u求傅里變換**********fs=1/ts。if nargin==2 n1=0。else n1=fs/df。endn2=length(sfm)。n=2^(max(nextpow2
環(huán)評公示相關推薦
鄂ICP備17016276號-1