利用DFT做連續(xù)信號的頻譜分析 一、利用DFT計算連續(xù)信號的頻譜

1、2.2 利用利用DFT做連續(xù)信號的頻譜分析做連續(xù)信號的頻譜分析 一、利用DFT計算連續(xù)信號的頻譜分析過程采樣截短DFT)(ajX)(atx( )x n)(jeX( )( )( )NNxnx n Rn)(jNeX)(kXNNkjNeX/2)(v 從圖中可以看出，這是一次次的近似過程，首先，用離散采樣信號的 DTFT 來近似連續(xù)信號 的傅立葉變換 ，其次，將x(n) 截短，這一過程等效于用一矩形序列RN(n)與x(n)相乘，其DTFT為v 最后，再對截短的信號作DFT)(jweX)(txa)( jXa)()()(jwNjwjwNeReXeX二、頻譜分析的近似過程出現(xiàn)問題及解決辦法1、混疊（1）混疊

2、現(xiàn)象對連續(xù)信號x(t)進行數(shù)字處理前，要進行采樣采樣序列的頻譜是連續(xù)信號頻譜的周期延拓，周期為fs，如采樣率過低，不滿足采樣定理，fs2fh，則導(dǎo)致頻譜混迭，使一個周期內(nèi)的譜對原信號譜產(chǎn)生失真，無法恢復(fù)原信號，進一步的數(shù)字處理失去依據(jù)。nanTttxnTx)()()(v（2）解決辦法 在采樣前利用一模擬低通濾波器將原始信號的上限頻率 限制在采樣頻率的一半，即加一抗混疊濾波器。maxf2、泄漏(1)產(chǎn)生原因處理實際信號序列x（n）時，一般總要將它截斷為一有限長序列，長為N點，相當(dāng)于乘以一個矩形窗w(n)=RN(n)。矩形窗函數(shù)，其頻譜有主瓣，也有許多副瓣，窗口越大，主瓣越窄，當(dāng)窗口趨于無窮大時，

3、就是一個沖擊函數(shù)。我們知道，時域的乘積對應(yīng)頻域的卷積，所以，加窗后的頻譜實際是原信號頻譜與矩形窗函數(shù)頻譜的卷積，卷積的結(jié)果使頻譜延伸到了主瓣以外，且一直延伸到無窮。當(dāng)窗口無窮大時，與沖擊函數(shù)的卷積才是其本身，這時無畸變，否則就有畸變。例如，信號為，是一單線譜，但當(dāng)加窗后，線譜與抽樣函數(shù)進行卷積，原來在0處的一根譜線變成了以0為中心的，形狀為抽樣函數(shù)的譜線序列，原來在一個周期（s）內(nèi)只有一個頻率上有非零值，而現(xiàn)在一個周期內(nèi)幾乎所有頻率上都有非零值，即的頻率成份從0處“泄漏”到其它頻率處去了?？紤]各采樣頻率周期間頻譜“泄漏”后的互相串漏，卷積后還有頻譜混迭現(xiàn)象產(chǎn)生。（2）解決辦法：增加截短的長度N

4、，以減小泄漏。nTje0Tjex3、柵欄效應(yīng)（1）柵欄效應(yīng)的產(chǎn)生N點DFT是在頻率區(qū)間0，2上對信號頻譜進行N點等間隔采樣，得到的是若干個離散的頻譜點X（k），且它們限制在基頻的整數(shù)倍上，這就好像在柵欄的一邊通過縫隙看另一邊的景象一樣，只能在離散點處看到真實的景象，其余部分頻譜成分被遮擋，所以稱之為柵欄效應(yīng)。（2）減小柵欄效應(yīng)方法：尾部補零，使譜線變密，增加頻域采樣點數(shù)，原來漏掉的某些頻譜分量就可能被檢測出來。4、 DFT的分辨率 v 填補零值可以改變對DTFT的采樣密度，人們常常有一種誤解，認(rèn)為補零可以提高DFT的頻率分辨率。事實上我們通常規(guī)定DFT的頻率分辨率為 ，這里的N是指信號x(n)

5、的有效長度，而不是補零的長度。不同長度的x(n)其DTFT的結(jié)果是不同的；而相同長度的x(n)盡管補零的長度不同其DTFT的結(jié)果應(yīng)是相同的，他們的DFT只是反映了對相同的DTFT采用了不同的采樣密度。 Nfs/Nfs/參數(shù)選擇的一般原則： (1)若已知信號的最高頻率 ，為防止混疊，選定采樣頻率 ；(2)根據(jù)頻率分辯率 ，確定所需DFT的長度(3) 和N確定以后，即可確定相應(yīng)模擬信號的時間長度這里T是采樣周期。maxfmax2ffsfffNs/NTfNs/sf5、周期信號的譜分析 對于連續(xù)的單一頻率周期信號 ， 為信號的頻率。 可以得到單一譜線的DFT結(jié)果，但這是和作DFT時數(shù)據(jù)的截取長度選得是

6、否恰當(dāng)有關(guān)，截取長度N選得合理， 可完全等于 的采樣。 )2sin()(tftxaaaf)(kXN)(ajXv例例2-1 對連續(xù)的單一頻率周期信號 v按采樣頻率 采樣，截取長度N分別選N =20和N =16，觀察其DFT結(jié)果的幅度譜。 v解解 此時離散序列 v，即k=8。用MATLAB計算并作圖，函數(shù)fft用于計算離散傅里葉變換DFT，程序如下： )2sin()(tftxaavk=8; vn1=0:1:19; vxa1=sin(2*pi*n1/k); vsubplot(2,2,1) vplot(n1,xa1) vxlabel(t/T);ylabel(x(n); vxk1=fft(xa1);xk

7、1=abs(xk1); vsubplot(2,2,2) vstem(n1,xk1) vxlabel(k);ylabel(X(k); vn2=0:1:15; vxa2=sin(2*pi*n2/k); vsubplot(2,2,3) vplot(n2,xa2) vxlabel(t/T);ylabel(x(n); vxk2=fft(xa2);xk2=abs(xk2); v subplot(2,2,4) vstem(n2,xk2) vxlabel(k);ylabel(X(k);051015-1-0.500.51t/Tx(n)0510150246810kX(k)051015-1-0.500.51t/Tx(n)0510150246810kX(k)(a)(b)(c)(d)不同截取長度的正弦信號及其DFT結(jié)果v計算結(jié)果示于圖 (a)和(b)分別是N=20