同濟大學信號與系統(tǒng)第九章第4講_dft

1、信號與系統(tǒng) 同濟大學汽車學院 魏學哲 9-7 9-7 離散傅里葉變換離散傅里葉變換的應用的應用運用運用DFT方法，往往伴隨方法，往往伴隨FFT算法的實施，所謂的應算法的實施，所謂的應用幾乎成為用幾乎成為FFT應用的同意語。應用的同意語。（一）快速卷積（一）快速卷積若長度為若長度為N1的序列的序列x (n)與長度為與長度為N2的序列的序列h(n)作線作線卷積，得到：卷積，得到：mmnhmxny)()()(y (n)為長度為為長度為N1+ N2 1的有限長序列，要做的有限長序列，要做N1*N2次乘法運算，當次乘法運算，當N1=N2 =N，要做，要做N2次乘法。次乘法。信號與系統(tǒng) 同濟大學汽車學院

2、魏學哲 直接卷積與快速卷積直接卷積與快速卷積如果把求線卷積改為求圓卷積，兩序列分別補零加長如果把求線卷積改為求圓卷積，兩序列分別補零加長為為為為N1+ N2 1，則有可能減少運算次數(shù)。，則有可能減少運算次數(shù)。x(n)y(n)卷積卷積FFTFFT序列相乘序列相乘IFFTx(n)h(n)X(k)H(k)X(k)H(k)y(n)信號與系統(tǒng) 同濟大學汽車學院 魏學哲 在快速卷積過程中，共需要在快速卷積過程中，共需要2次次FFT，1次次IFFT，相當，相當于于3次次FFT運算量，在一般數(shù)字濾波中，由運算量，在一般數(shù)字濾波中，由h(n)求求H(k)是事先計算完成放在存儲器中，故只需是事先計算完成放在存儲器

3、中，故只需2次次FFT的運算的運算量，若假設量，若假設N1=N2 =N ，則補零后：，則補零后：NNN2121故需要故需要NN22log2此外此外X(k)與與H(k)兩序列相乘，還需要兩序列相乘，還需要2N次復乘，次復乘，全部復數(shù)乘法次數(shù)為：全部復數(shù)乘法次數(shù)為：次復數(shù)乘法運算次復數(shù)乘法運算)log1 (22log22222NNNNN隨著隨著N的增大，該運算量比的增大，該運算量比N2顯著減少。顯著減少。信號與系統(tǒng) 同濟大學汽車學院 魏學哲 以上分析是針對兩序列長度相近或相等的情況，如以上分析是針對兩序列長度相近或相等的情況，如果一個序列很短，而另一序列很長，則需要補很多果一個序列很短，而另一序列

4、很長，則需要補很多零，圓卷方案甚至增加運算量，可采用分段卷積的零，圓卷方案甚至增加運算量，可采用分段卷積的方法，其基本原理是將方法，其基本原理是將x(n)分成若干小段，每小段分成若干小段，每小段長度與長度與h(n)接近，將接近，將x(n)的各小段與的各小段與h(n)卷積，最卷積，最后取和，仍可發(fā)揮快速卷積的優(yōu)越性。后取和，仍可發(fā)揮快速卷積的優(yōu)越性。信號與系統(tǒng) 同濟大學汽車學院 魏學哲 重疊相加法重疊相加法 假定h(n), x(n)均為因果序列。h(n)的長度為N，如圖a x(n)長度是 現(xiàn)將 等分為若干小段，每段長M,如圖bNN11,NNh(n)圖a1N圖bMMMN1X(n)信號與系統(tǒng) 同濟大

5、學汽車學院 魏學哲 以 表示x（n）序列的第i小段（i為正整數(shù)， ），為完成 各 與h（n）之圓卷積，應將 長度補足至N+M-1，在圖c中已用虛線示意補零。MNP1,-Pi01 nxi nxi nxiN+M-1N+M-1N+M-1N+M-1 nxi圖c信號與系統(tǒng) 同濟大學汽車學院 魏學哲 輸入序列可表示為 其中 相應地，輸出序列也可分解為 10Piinxnx 為其他值nMiniMnxnxi01) 1( nhnxnhnxnyPii10信號與系統(tǒng) 同濟大學汽車學院 魏學哲 利用卷積分配律可得 由于 的長度為N+M-1，而 的有效長度只有M，故相鄰兩段的 必有N-1長度的重疊。見圖dnhnxnyny

6、nhnxnyiiPiiPii這里1010N+M-1重疊N+M-1重疊 niy圖d重疊N+M-1 nyi nyi nxi信號與系統(tǒng) 同濟大學汽車學院 魏學哲 按照上述原理，此方法的運算過程可分為兩部分： 首先求每個 與h(n)的圓卷積，樣點數(shù)為N+M-1，共需P次，求各 ，其原理仍按圖b； 然后將 取和（實際上是重疊部分的相加），即得 nxi nyi nyi ny圖e y(n)N+N1-1信號與系統(tǒng) 同濟大學汽車學院 魏學哲 有時，N1可能很長，以致趨于無限大，例如語音信號、地震波動信號、宇宙通信中產(chǎn)生的某些信號等，如果不采用分段卷積的方法將遲遲不能給出結(jié)果，而且，無法找到那樣大的存儲設備來滿足

7、N1的需要。因此，即使在分段措施改善速度不顯著的情況下，仍有可能采用這種方法。 借助FFT不僅可完成快速卷積運算，也可利用它進行解卷積運算，具體計算公式可仿照z變換解卷積方法求得。信號與系統(tǒng) 同濟大學汽車學院 魏學哲 （二）快速相關(guān)（二）快速相關(guān)相關(guān)和自相關(guān)也可借助相關(guān)和自相關(guān)也可借助FFT完成。完成。x(n)y(n)相關(guān)相關(guān)FFTFFT序列相乘序列相乘IFFTx(n)h(n)X(k)H(k)X(k)H*(k)y(n)h(n)共軛共軛H*(k)信號與系統(tǒng) 同濟大學汽車學院 魏學哲 相關(guān)運算常見于雷達和聲納系統(tǒng)中，應用該運算確相關(guān)運算常見于雷達和聲納系統(tǒng)中，應用該運算確定信號的時間延遲。定信號的

8、時間延遲。當當x(n)與與h(n)是同一信號，是同一信號， y(n) 是自相關(guān)函數(shù)，是自相關(guān)函數(shù)，而而Y(k)是是x(n)的功率譜。的功率譜。信號與系統(tǒng) 同濟大學汽車學院 魏學哲 （三）利用（三）利用DFTDFT對連續(xù)時間信號的逼近對連續(xù)時間信號的逼近hsff2hsffT211其中其中,T為抽樣間隔。為抽樣間隔?；蛘呋蛘咭灰?用用DFT計算連續(xù)時間信號的傅氏變換可能造成的誤計算連續(xù)時間信號的傅氏變換可能造成的誤差差1.混疊現(xiàn)象混疊現(xiàn)象為避免混疊為避免混疊,由抽樣定理可知，須滿足由抽樣定理可知，須滿足其中其中fs為抽樣頻率為抽樣頻率; fh為信號的最高頻率分量為信號的最高頻率分量對于時間有限信號

9、，其傅立葉變換不可能是有限帶對于時間有限信號，其傅立葉變換不可能是有限帶寬，抽樣后必然帶來混疊寬，抽樣后必然帶來混疊(aliasing)，減小抽樣間，減小抽樣間隔可減弱混疊，但總不可避免。隔可減弱混疊，但總不可避免。信號與系統(tǒng) 同濟大學汽車學院 魏學哲 若信號頻譜有限，則時間函數(shù)必然是無限的。在實若信號頻譜有限，則時間函數(shù)必然是無限的。在實際應用中，為利用際應用中，為利用FFT對信號進行分析，必須把時對信號進行分析，必須把時間截取一定范圍，也就是說，在時域?qū)π盘栠M行截間截取一定范圍，也就是說，在時域?qū)π盘栠M行截斷操作，或稱作加時間窗，亦即用時間窗函數(shù)乘以斷操作，或稱作加時間窗，亦即用時間窗函數(shù)

10、乘以信號，由卷積定理可知，時域相乘，頻域為卷積，信號，由卷積定理可知，時域相乘，頻域為卷積，時間加窗使頻譜產(chǎn)生失真，它從原有的頻率受限圖時間加窗使頻譜產(chǎn)生失真，它從原有的頻率受限圖中擴散出來，這就造成拖尾現(xiàn)象，稱之為頻譜泄漏中擴散出來，這就造成拖尾現(xiàn)象，稱之為頻譜泄漏（leakage）。）。2.2.頻譜泄漏頻譜泄漏信號與系統(tǒng) 同濟大學汽車學院 魏學哲 0n0)(1nx)(1jeXn)(2nx)(2jeXn )()(21nxnxnyjeY信號與系統(tǒng) 同濟大學汽車學院 魏學哲 pTF13.3.柵欄效應柵欄效應用用DFT計算頻譜時，只是知道為頻率計算頻譜時，只是知道為頻率的整數(shù)倍處的頻譜。在兩個譜線

11、之間的情況就不的整數(shù)倍處的頻譜。在兩個譜線之間的情況就不知道知道,這相當通過一個柵欄觀察景象一樣，故稱作這相當通過一個柵欄觀察景象一樣，故稱作柵欄效應。補零點加大周期，可使柵欄效應。補零點加大周期，可使F變小來提高變小來提高分辨力，以減少柵欄效應。分辨力，以減少柵欄效應。信號與系統(tǒng) 同濟大學汽車學院 魏學哲 例例 有一頻譜分析用的有一頻譜分析用的FFT處理器，其抽樣點數(shù)必須是處理器，其抽樣點數(shù)必須是2的整數(shù)冪。的整數(shù)冪。 假定沒有采用任何特殊的數(shù)據(jù)處理措施，已知條件為假定沒有采用任何特殊的數(shù)據(jù)處理措施，已知條件為（1）頻）頻 率分辨率為率分辨率為 ，（2） 信號的最高頻率信號的最高頻率 ，試確

12、定以下參量：（試確定以下參量：（1）最小記錄長度）最小記錄長度 ;(2) 抽樣點間的最大時間抽樣點間的最大時間間隔間隔T; (3) 在一個記錄中的最小點數(shù)在一個記錄中的最小點數(shù)N。ZH10ZkH4PT解：解： （a） 最小記錄長度最小記錄長度sFTP1 . 01011sTP1 . 0，（b）最大的抽樣時間間隔）最大的抽樣時間間隔TsffThs3310125. 01042/ 1/ 1/ 1msT1 . 0(c) 最小記錄點數(shù)最小記錄點數(shù)N1024280010/1042/2103NFfNh取信號與系統(tǒng) 同濟大學汽車學院 魏學哲 dejXtxdtetxjXtjtj21)(2.連續(xù)時間周期信號傅氏級數(shù)

13、變換對連續(xù)時間周期信號傅氏級數(shù)變換對 ktjkTTtjkejkXtxdtetxTjkXpp0002/2/01二二.DFT.DFT與連續(xù)信號傅氏變換相對數(shù)值的確定與連續(xù)信號傅氏變換相對數(shù)值的確定1.連續(xù)時間非周期信號傅氏變換對連續(xù)時間非周期信號傅氏變換對信號與系統(tǒng) 同濟大學汽車學院 魏學哲 10 ,)(1)(10 ,)()(1010NnWkXNnxNkWnxkXNknkNNnnkN3.DFT變換時變換時:信號與系統(tǒng) 同濟大學汽車學院 魏學哲 1sfT4.4.用用DFTDFT計算非周期信號的傅氏變換計算非周期信號的傅氏變換幅度電平未受到影響。幅度電平未受到影響。用用DFT計算所得的頻譜分量乘以計算

14、所得的頻譜分量乘以T， 就等于頻譜的正就等于頻譜的正常幅度電平；用常幅度電平；用IDFT計算非周期信號的傅氏反變換，計算非周期信號的傅氏反變換，再乘以再乘以fs就得到所需信號的正常幅度電平。所以就得到所需信號的正常幅度電平。所以,從從時間到頻率時間到頻率, 再從頻率到時間，整個過程總共乘了再從頻率到時間，整個過程總共乘了信號與系統(tǒng) 同濟大學汽車學院 魏學哲 設設nTdtnTt,10)()()()(NnnTjnnTjtjenTxTenTxdtetxjX用用DFT計算所得的頻譜分量乘以計算所得的頻譜分量乘以T的理由：的理由：信號與系統(tǒng) 同濟大學汽車學院 魏學哲 )()(1,.1 , 0,)()(10/210/20nxDFTTenxTNkenTxTkXNnNknjNnNfTjkns1,.1 , 0,/200NkkNfs又信號與系統(tǒng) 同濟大學汽車學院 魏學哲 0000,2,2kdNNFfFs)()(1)(2)(21)(1102021000000kXDFTfekjXNfekjXejkXnTxsNkknNjsknNTfjNkknTjks 用用IDFT計算非周期