同濟大學(xué)信號及系統(tǒng)第九章第3講2-dft

一、有限長序列Z變換的抽樣若有限長序列x(n)的長度為N，其Z變換的表示式為：若有限長序列滿足絕對可和條件，則其收斂域包括單位圓在內(nèi)，在單位圓上N個等間距點計算Z變換，令：§9-5離散傅里葉變換與Z變換的關(guān)系1weixzh@在Z平面單位圓上，取輻角為的等間距的第k個點，計算其Z變換。因此，有限長序列的DFT可解釋為它的Z變換在單位圓上的均勻抽樣。1weixzh@二.以如何從頻域抽樣恢復(fù)原序列1.兩種抽樣時域抽樣:

對一個頻帶有限的信號，根據(jù)抽樣定理對其進行抽樣，所得抽樣信號的頻譜是原帶限信號頻譜的周期延拓，因此，完全可以由抽樣信號恢復(fù)原信號。頻域抽樣:對一有限序列(時間有限序列)進行DFT所得X(k)就是序列傅氏變換的采樣。所以DFT就是頻域抽樣。1weixzh@2.由頻域抽樣恢復(fù)序列由于x(n)絕對可和，故其傅氏變換存在且連續(xù)，也即其Z變換收斂域包括單位圓。這樣，對X(Z)在單位圓上N等份抽樣，就得到X(k)。一個絕對可和的非周期序列x(n)的Z變換為1weixzh@對X(k)進行反變換,并令其為xN(n),則1weixzh@是非周期序列x(n)的周期延拓，也就是說，頻域抽樣造成時域周期延拓?？梢姡傻玫降闹芷谛蛄衅渲?weixzh@當(dāng)x(n)不是有限長時，無法周期延拓；當(dāng)x(n)的長度M,只有NM時,才能不失真的恢復(fù)信號,即3.頻域抽樣不失真的條件1weixzh@1.由X(k)恢復(fù)X(Z)序列x(n),（0nN-1）的Z變換為二.由X(k)表達X(Z)-內(nèi)插公式1weixzh@1weixzh@稱作內(nèi)插函數(shù)。上式就是由X(k)恢復(fù)X(Z)的內(nèi)插公式,其中1weixzh@。。。。。。。將內(nèi)插函數(shù)寫成如下式:2.內(nèi)插函數(shù)的特性1weixzh@令分子為零，得:所以有N個零點。為一階極點，Z=0為(N-1)階極點。但是極點與一零點相消。這樣只有(N-1)個零點,抽樣點稱作本抽樣點。因此說,內(nèi)插函數(shù)僅在本抽樣點處不為零,其他(N-1)個抽樣點均為零。令分母為零，得1weixzh@代入4.內(nèi)插函數(shù)的頻率特性3.頻率響應(yīng)單位圓上的Z變換即為頻響,1weixzh@的函數(shù)又是k的函數(shù),其可表示為當(dāng)k=0時,則有可見,既是1weixzh@所以時時這說明在本抽樣點為1，在其他抽樣點為0。1weixzh@其中，相位特性幅度特性當(dāng)N=5時,的1weixzh@N=51weixzh@1weixzh@而在其他抽樣點上這就是說,內(nèi)插函數(shù)在本抽樣點由于i與k均為整數(shù),所以時1weixzh@就精確等于X(k)。即的特性可知,在每個抽樣點上其值為1,故由于與X(k)的關(guān)系1weixzh@而在抽樣點之間等于加權(quán)的內(nèi)插函數(shù)值疊加而得。1weixzh@§9-7離散傅里葉變換的應(yīng)用運用DFT方法，往往伴隨FFT算法的實施，所謂的應(yīng)用幾乎成為FFT應(yīng)用的同意語。（一）快速卷積若長度為N1的序列x(n)與長度為N2的序列h(n)作線卷積，得到：y(n)為長度為N1+N2–1的有限長序列，要做N1*N2次乘法運算，當(dāng)N1=N2=N，要做N2次乘法。1weixzh@直接卷積與快速卷積如果把求線卷積改為求圓卷積，兩序列分別補零加長為為N1+N2–1，則有可能減少運算次數(shù)。x(n)y(n)卷積FFTFFT序列相乘IFFTx(n)h(n)X(k)H(k)X(k)H(k)y(n)1weixzh@在快速卷積過程中，共需要2次FFT，1次IFFT，相當(dāng)于3次FFT運算量，在一般數(shù)字濾波中，由h(n)求H(k)是事先計算完成放在存儲器中，故只需2次FFT的運算量，若假設(shè)N1=N2=N，則補零后：故需要此外X(k)與H(k)兩序列相乘，還需要2N次復(fù)乘，全部復(fù)數(shù)乘法次數(shù)為：次復(fù)數(shù)乘法運算隨著N的增大，該運算量比N2顯著減少。1weixzh@以上分析是針對兩序列長度相近或相等的情況，如果一個序列很短，而另一序列很長，則需要補很多零，圓卷方案甚至增加運算量，可采用分段卷積的方法，其基本原理是將x(n)分成若干小段，每小段長度與h(n)接近，將x(n)的各小段與h(n)卷積，最后取和，仍可發(fā)揮快速卷積的優(yōu)越性。1weixzh@（二）快速相關(guān)相關(guān)和自相關(guān)也可借助FFT完成。x(n)y(n)相關(guān)FFTFFT序列相乘IFFTx(n)h(n)X(k)H(k)X(k)H*(k)y(n)h(n)共軛H*(k)1weixzh@相關(guān)運算常見于雷達和聲納系統(tǒng)中，應(yīng)用該運算確定信號的時間延遲。當(dāng)x(n)與h(n)是同一信號，y(n)是自相關(guān)函數(shù)，而Y(k)是x(n)的功率譜。1weixzh@（三）利用DFT對連續(xù)時間信號的逼近其中,T為抽樣間隔。或者一.用DFT計算連續(xù)時間信號的傅氏變換可能造成的誤差1.混疊現(xiàn)象為避免混疊,由抽樣定理可知，須滿足其中fs為抽樣頻率;fh為信號的最高頻率分量對于時間有限信號，其傅立葉變換不可能是有限帶寬，抽樣后必然帶來混疊(aliasing)，減小抽樣間隔可減弱混疊，但總不可避免。1weixzh@若信號頻譜有限，則時間函數(shù)必然是無限的。在實際應(yīng)用中，為利用FFT對信號進行分析，必須把時間截取一定范圍，也就是說，在時域?qū)π盘栠M行截斷操作，或稱作加時間窗，亦即用時間窗函數(shù)乘以信號，由卷積定理可知，時域相乘，頻域為卷積，時間加窗使頻譜產(chǎn)生失真，它從原有的頻率受限圖中擴散出來，這就造成拖尾現(xiàn)象，稱之為頻譜泄漏（leakage）。2.頻譜泄漏1weixzh@0n0nn1weixzh@3.柵欄效應(yīng)用DFT計算頻譜時，只是知道為頻率的整數(shù)倍處的頻譜。在兩個譜線之間的情況就不知道,這相當(dāng)通過一個柵欄觀察景象一樣，故稱作柵欄效應(yīng)。補零點加大周期，可使F變小來提高分辨力，以減少柵欄效應(yīng)。1weixzh@[例]有一頻譜分析用的FFT處理器，其抽樣點數(shù)必須是2的整數(shù)冪。假定沒有采用任何特殊的數(shù)據(jù)處理措施，已知條件為（1）頻率分辨率為，（2）信號的最高頻率，試確定以下參量：（1）最小記錄長度;(2)抽樣點間的最大時間間隔T;(3)在一個記錄中的最小點數(shù)N。解：（a）最小記錄長度（b）最大的抽樣時間間隔T(c)最小記錄點數(shù)N1weixzh@2.連續(xù)時間周期信號傅氏級數(shù)變換對二.DFT與連續(xù)信號傅氏變換相對數(shù)值的確定1.連續(xù)時間非周期信號傅氏變換對1weixzh@3.DFT變換時:1weixzh@4.用DFT計算非周期信號的傅氏變換幅度電平未受到影響。用DFT計算所得的頻譜分量乘以T，就等于頻譜的正常幅度電平；用IDFT計算非周期信號的傅氏反變換，再乘以fs就得到所需信號的正常幅度電平。所以,從時間到頻率,再從頻率到時間，整個過程總共乘了1weixzh@設(shè)用DFT計算所得的頻譜分量乘以T的理由：1weixzh@1weixzh@mai