“數(shù)字信號處理”中DFT譜分析的教學探討

摘要：針對“數(shù)字信號處理”課程中學生普遍難以掌握的用離散傅里葉變換（DFT）對連續(xù)信號譜分析的問題，提出幾個關(guān)鍵的教學要點，即明確原理、分析誤差、示例驗證、計算參數(shù)、綜合問題、獲得頻譜。教學實踐表明，按此過程講授，可使學生從根本上增強對這一問題的理解，為解決DFT在工程實際中的應(yīng)用打好基礎(chǔ)。

關(guān)鍵詞：離散傅里葉變換（DFT）；譜分析；數(shù)字信號處理；教學

作者簡介：喬建華（1975-），女，山西離石人，太原科技大學電子信息工程學院，副教授；于少娟（1971-），女，江蘇漣水人，太原科技大學電子信息工程學院，教授。（山西#8194;太原#8194;030024）

基金項目：本文系山西省教育廳教學改革研究項目（項目編號：2009sjg202）、山西省高等學校教學改革項目（項目編號：J2011130）、山西省高等學校教學改革項目（項目編號：J2013064）、山西省特色專業(yè)建設(shè)項目（項目編號：201201）的研究成果。

中圖分類號：G642.0#8195;#8195;#8195;#8195;#8195;文獻標識碼：A#8195;#8195;#8195;#8195;#8195;文章編號：1007-0079（2014）14-0114-02

離散傅里葉變換（DFT）開辟了頻域離散化的道路，[1]DFT的快速算法FFT的出現(xiàn)，又使得DFT的應(yīng)用從理論走向了實踐，帶來了數(shù)字信號處理技術(shù)突飛猛進的發(fā)展，因此，DFT及其快速算法FFT具有重要的理論及現(xiàn)實意義。DFT的應(yīng)用之一——用DFT（FFT）對連續(xù)信號譜分析就是本科“數(shù)字信號處理”課程中的一個重點內(nèi)容，同時也是一個難點，學生普遍對這部分內(nèi)容理解不深入，對誤差產(chǎn)生原因不明確，對減小誤差的方法不清晰。究其原因，還是對整個譜分析過程研究不透徹的緣故。筆者經(jīng)過多年的教學改革，總結(jié)出講解這部分內(nèi)容的幾個要點，可使學生更深入明確地掌握該內(nèi)容，并加深對DFT的理解，為DFT在工程實踐中的正確應(yīng)用打好基礎(chǔ)。

一、梳理過程，明確原理

對連續(xù)信號進行譜分析，關(guān)鍵就是得到信號的傅里葉變換。[2]而要用計算機計算，就需要先進行時域采樣，得到采樣信號，再經(jīng)過傅里葉變換轉(zhuǎn)換到頻域，然后頻域采樣，而頻域采樣就可以用DFT來完成。這個過程如圖1所示，模擬信號xa（t）的譜分析，是按實線箭頭方向的過程實現(xiàn)的，虛框中的兩個步驟實際上是由虛框上方的DFT（FFT）來完成的。

按照上述過程推導可得，用DFT對連續(xù)信號譜分析的計算公式為：[2]

， （1）

即可用對連續(xù)信號采樣得到的時域離散信號x（n）的N點DFT再乘以采樣間隔T近似得到。對這個公式的講解，主要注重其原理和概念，掌握DFT在其中的應(yīng)用，并且要說明在以上分析的步驟中存在著近似。

二、提出問題，定性分析，指出誤差產(chǎn)生必然性

可以發(fā)現(xiàn)，在上述分析的過程中存在幾個問題，即圖1中向下的虛線箭頭指向的問題。一是對模擬信號采樣時，需要滿足Nyquist采樣定理，否則會產(chǎn)生頻譜混疊。然而模擬信號大多不是銳截止的，因為根據(jù)傅里葉變換理論，若信號持續(xù)時間有限長，則其頻譜無限寬；若信號的頻譜有限寬，則其持續(xù)時間無限長。[1]所以對于高頻的幅度很小的信號，只能近似看成帶限信號，因此在時域采樣時，必然會出現(xiàn)頻譜混疊現(xiàn)象。二是采樣的信號要進行DFT分析，需要是有限長的，因此要對采樣的序列進行截斷，也會產(chǎn)生誤差。三是頻譜采樣時，相當于隔著柵欄看信號，有些信號會被遮擋，這樣就有可能漏掉一些重要的成分，產(chǎn)生誤差。

要說明的是，以上這幾方面的誤差都是不可避免的，只是從工程角度看，在誤差允許的范圍內(nèi)。那么誤差對譜分析的影響有多大，如何減小這些誤差，就是接下來要重點分析的問題。

三、深入探討，示例驗證，尋求減小誤差的措施

針對以上過程產(chǎn)生的誤差，下面來逐一分析，以探求抑制誤差的措施。下表1即詳細列出了誤差的產(chǎn)生及減小誤差的方法。

表1#8194;用DFT對連續(xù)信號譜分析的誤差問題一覽表

產(chǎn)生原因誤差類型減小誤差的方法

時域采樣頻譜混疊預(yù)濾，濾除高于折疊頻率fs/2的高頻雜散分量；

要滿足Nyquist采樣定理，一般取fs≥(3～5)fc

截取有限長序列截斷效應(yīng)增大截取長度，提高頻率分辨率，減小高頻泄漏；

選擇合適的窗函數(shù)，但窄的主瓣寬度和低的旁瓣幅度只能折中選取

頻域采樣柵欄效應(yīng)原序列尾部補零；

增加截取長度

具體來說，首先時域采樣時，由于不是帶限信號，肯定會造成混疊。為了滿足Nyquist采樣定理，可以先進行預(yù)濾，濾除高于折疊頻率fs/2的高頻雜散信號，然后再按照采樣頻率fs≥2fc的原則，[3]fc為模擬信號的最高頻率，而實際應(yīng)用中fs一般?。?～5）fc。

其次要截取有限長序列，需要通過加窗來截斷，由此造成了截斷效應(yīng)，它有高頻泄露和譜間干擾兩種表現(xiàn)方式。[4]對于這個問題，先要講清窗函數(shù)的性質(zhì)，窗函數(shù)的頻譜主要關(guān)注兩點：一是主瓣。例如對于矩形窗，就是指ω=-2π/N和2π/N之間的區(qū)域，因此矩形窗的主瓣寬度就是4π/N，[5]它使得譜線展寬，即產(chǎn)生了高頻泄漏，影響的是頻率分辨率。由于每一波瓣下的面積是恒定不變的，隨著N的增大，主瓣會變窄，幅度會增加。二是旁瓣。它產(chǎn)生的波紋造成了譜間干擾，大的旁瓣可能會湮沒小的主瓣，因此旁瓣的幅度也要設(shè)法減小。但是對確定的窗函數(shù)，旁瓣的衰減是確定的。對不同的窗函數(shù)，當N一定時，主瓣寬度窄的，旁瓣幅度就高；旁瓣幅度低的，主瓣寬度就寬，所以要根據(jù)需要折中選取。掌握了窗函數(shù)的特性，也就能明白如何減小截斷效應(yīng)帶來的誤差，也即從增大截取長度N和選擇合適的窗函數(shù)兩個方面來解決這個問題。

對于頻域采樣時產(chǎn)生的柵欄效應(yīng)，一方面要減少遮擋，可以通過給原序列尾部補零，增加DFT變換的點數(shù)；另一方面是直接增大截取長度，將頻譜本身細化，頻率分辨率自然隨之提高。這個問題可以舉一些具體的例子，通過MATLAB演示來說明，[6]更能直觀地體現(xiàn)這兩種方法本質(zhì)的不同，加深學生的理解。

例如分析信號x（n）=cos（40nπ）+cos（41nπ）+sin（80nπ），則要分辨出20Hz和20.5Hz的兩條譜線，頻率分辨率至少就要達到0.5Hz。由于信號最高頻率為40Hz，采樣頻率取100Hz，則最少采樣點數(shù)為100/0.5=200，這里取256點。若點數(shù)太少，頻譜本身就是不清晰的，DFT的變換點數(shù)再多也無用。下圖2顯示了不同采樣點數(shù)及序列尾部補零時對應(yīng)的頻譜圖，從中可以清晰地看出不同情況下信號頻譜的分布。

四、回歸問題，計算參數(shù)，得到譜分析結(jié)果

總結(jié)以上過程，根據(jù)譜分析原理及誤差減小的方法，可以得到用DFT對連續(xù)信號譜分析的具體步驟：預(yù)濾，確定模擬信號最高頻率fc；時域采樣，采樣頻率fs≥2fc，一般取fs≥（3～5）fc；確定截取長度，用記錄時間Tp表示，由于頻率分辨率F=fs/N，最少記錄時間Tpmin=NT=N/fs=1/F；選擇合適的窗函數(shù)；確定最少采樣點數(shù)Nmin=fs/F；按照式（1），得到模擬信號的離散頻譜函數(shù)Xa（k）。

另外需要說明的是，根據(jù)頻率采樣定理，只有頻率采樣點數(shù)N大于等于原序列長度M，才能無失真地恢復原時域信號，所以通常將序列的截取點數(shù)也作為N。而且，在實際的應(yīng)用中，上述參數(shù)的取值一般都要大些，以保證獲得準確的頻域分析結(jié)果。一種有效的做法是，將取值按倍數(shù)擴大后，看頻譜是否出現(xiàn)新的信息，如此反復對比，直至獲得精準結(jié)果。

五、結(jié)語

“數(shù)字信號處理”是一門理論和實踐并重的課程，因此在教學中，既要講清原理，又要多以例證。對于本節(jié)內(nèi)容，這兩方面尤其要重視。本文即針對工程實踐中經(jīng)常需要做的用DFT對連續(xù)信號譜分析的問題，亦是“數(shù)字信號處理”課程中的難點之一，做了細致的分析，指出講授這部分內(nèi)容的幾個關(guān)鍵要素：先要宏觀把握，明確整個譜分析的過程；再通過定性分析指出誤差產(chǎn)生的必然性；然后可以通過不同示例尋求減小誤差的具體方法；最后一定要回到問題上，計算出所需參數(shù)，最終得到準確的頻譜，從而可使學生脫離只會依據(jù)例子計算參數(shù)的簡單層次，理清整個譜分析過程的脈絡(luò)，明確各環(huán)節(jié)之間的關(guān)系，獲得精準的頻譜分析結(jié)果。實踐表明，取得了良好的教學效果。

參考文獻：

[1]A.V.Oppenheim， R.W.Schafer.數(shù)字信號處理[M].北京:科學出版社，1980.

[2]高西全，丁玉美.數(shù)字信號處理[M].第三版.西安:西安電子科技出版社，2008.

[3]賈中云，李秀梅，董文.“數(shù)字信號處理”中采樣定理的教學探索[J].中國電力教育，2012，(29):52-53.

[4]謝海霞，孫志雄.連續(xù)非周期信號頻譜分析及Matlab實現(xiàn)[J].現(xiàn)代電子技術(shù)，2013，(11):53-56.

[5]Sophocles J.Orfanidis.信號處理導論[M].北京:清華大學出版社，

1999.

[6]黃永平，田秀麗，田曉燕.“數(shù)字信號處理”課程教學改革與實踐[J].電氣電子教學學報，2013，(2):39-40，43.

（責任編輯：王意琴）