第4章離散傅里葉變換培訓(xùn)教材

離散傅里葉變換(DiscreteFourierTransform，簡(jiǎn)稱DFT)是數(shù)字信號(hào)處理中非常重要的一種數(shù)學(xué)變換，本章主要討論DFT的定義、物理意義及基本性質(zhì)、頻率采樣和DFT應(yīng)用舉例。主要內(nèi)容包括：?周期序列的離散傅里葉級(jí)數(shù)定義及其性質(zhì)?有限長(zhǎng)度序列的離散傅里葉變換推導(dǎo)?離散傅里葉變換的基本性質(zhì)?頻率域采樣理論?離散傅里葉變換的基本應(yīng)用

第3章離散傅里葉變換2025/1/2423.1離散傅里葉變換的定義3.2離散傅里葉變換的性質(zhì)3.3頻率域采樣3.4DFT的應(yīng)用舉例2025/1/2433.1離散傅里葉變換定義

3.1.1傅里葉變換的幾種可能形式 傅里葉變換就是建立以時(shí)間為自變量的“信號(hào)”與以頻率為自變量的“頻譜函數(shù)”之間的某種變換關(guān)系。所以，當(dāng)自變量“時(shí)間”或“頻率”取連續(xù)或離散值時(shí)，就形成各種不同形式的傅里葉變換對(duì)。1、連續(xù)時(shí)間、連續(xù)頻率——傅里葉變換就是連續(xù)時(shí)間非周期信號(hào)x(t)的傅里葉變換關(guān)系,所得到的是連續(xù)的非周期的頻譜函數(shù)。2025/1/245圖3-1連續(xù)的非周期信號(hào)及其非周期、連續(xù)的頻譜密度2025/1/2462、連續(xù)時(shí)間、離散頻率——傅里葉級(jí)數(shù)

設(shè)x(t)代表一個(gè)周期為T0的周期性連續(xù)時(shí)間函數(shù)，可展成傅里葉級(jí)數(shù)，其傅里葉級(jí)數(shù)的系數(shù)為X(jkΩ0)，是離散頻率的非周期函數(shù)，和X(jkΩ0)組成變換對(duì)，表示為其中Ω0=2πF=2π/T0為離散頻譜相鄰兩譜線之間的角頻率間隔，k為諧波序號(hào)。

3-33-42025/1/247圖3-2連續(xù)的周期信號(hào)及其非周期離散譜線

這一變換對(duì)的示意圖如圖3-2所示?？梢钥闯?，時(shí)域連續(xù)函數(shù)造成頻域是非周期的頻譜函數(shù)，而頻域的離散頻譜就與時(shí)域的周期時(shí)間函數(shù)相對(duì)應(yīng)。2025/1/2483、離散時(shí)間、連續(xù)頻率——序列的傅里葉變換

這正是第二章中討論過(guò)的序列（離散時(shí)間信號(hào)）的傅里葉變換，即3-53-6上面討論的三種傅里葉變換對(duì)都不適合在計(jì)算機(jī)上運(yùn)算，因?yàn)樗鼈冎辽僭谝粋€(gè)域(時(shí)域或頻域)中函數(shù)是連續(xù)的。因而從數(shù)字計(jì)算角度出發(fā)，我們感興趣的是時(shí)域及頻域都是離散的情況，這就是我們這里要研究的離散傅里葉變換。

2025/1/249離散傅立葉變換—有限長(zhǎng)序列的離散頻域表示

離散傅里葉變換的定義設(shè)x(n)是一個(gè)長(zhǎng)度為N的有限長(zhǎng)序列，則定義x(n)的N點(diǎn)離散傅里葉變換為3-313-32(3-31)和(3-32)分別是有限長(zhǎng)序列的離散傅里葉正變換和反變換。x(n)與X(k)是一個(gè)有限長(zhǎng)序列的離散傅里葉變換對(duì)。已知其中的一個(gè)序列，就能唯一確定另一序列。離散傅里葉變換記作DFT，IDFT為逆變換。2025/1/24103.2離散傅里葉變換的性質(zhì)3.2.1線性性質(zhì)如果和是兩個(gè)有限長(zhǎng)序列，長(zhǎng)度分別為N1和N2。若式中a，b為任意常數(shù)，取N=max[N1,N2]，則的N點(diǎn)DFT為3-35其中X1(k)和X2(k)分別為x1(n)和x2(n)的N點(diǎn)DFT。2025/1/24113.2.2循環(huán)移位性質(zhì)

1、序列的循環(huán)移位設(shè)為x(n)有限長(zhǎng)序列，長(zhǎng)度為N，則的循環(huán)移位定義為3-36式(3-36)表明，將x(n)以N為周期進(jìn)行周期延拓得到再將左移m位得到，最后取的主值序列，則得到有限長(zhǎng)序列的循環(huán)移位序列y(n)。x(n)及其循環(huán)移位過(guò)程如圖3-7所示。顯然y(n)仍是長(zhǎng)度為N有限長(zhǎng)序列。觀察圖3-7可見(jiàn)，循環(huán)移位的實(shí)質(zhì)是將x(n)左移m位，而左側(cè)移出主值區(qū)的序列值又依次從右側(cè)進(jìn)入主值區(qū)。2025/1/2412圖3-7循環(huán)移位過(guò)程示意圖

2025/1/24132、時(shí)域循環(huán)移位定理設(shè)x(n)是長(zhǎng)度為N的有限序列，y(n)為x(n)的循環(huán)移位，即則3-37其中證明：2025/1/2414令n+m=n’，則有由于上式中求和項(xiàng)以N為周期，所以對(duì)其在任一周期上的求和結(jié)果相同。將上式的求和區(qū)間改在主值區(qū)則得2025/1/24153、頻域循環(huán)移位定理如果X(k)=DFT[x(n)]，0≤k≤N-1Y(k)=X((k+l))NRN(k)則3-38式(3-38)的證明方法與時(shí)域循環(huán)定理類似，直接對(duì)Y(k)=X((k+l))NRN(k)進(jìn)行IDFT即得證。2025/1/24163.2.3循環(huán)卷積定理有限長(zhǎng)序列x1(n)和x2(n)，長(zhǎng)度分別為N1和N2，N=max[N1，N2]。x1(n)和x2(n)的N點(diǎn)DFT分別為

X1(k)=DFT[x1(n)]，X2(k)=DFT[x2(n)]如果X(k)=X1(k)·X2(k)則3-39或一般稱(3-39)式所表示的運(yùn)算為和的循環(huán)卷積。2025/1/2417證明：直接對(duì)(3-39)式兩邊進(jìn)行DFT令n-m=n’，則有2025/1/2418因?yàn)樯鲜街幸訬為周期，所以對(duì)其在任一周期上的求和結(jié)果不變。因此式(3-39)的循環(huán)卷積過(guò)程如圖3-8所示。2025/1/2419圖3-8循環(huán)卷積過(guò)程示意圖2025/1/2420由于循環(huán)卷積過(guò)程中，要求對(duì)x2(m)的循環(huán)反轉(zhuǎn)，循環(huán)移位，特別是兩個(gè)N長(zhǎng)的序列的循環(huán)卷積長(zhǎng)度仍為N。顯然與一般的線性卷積不同，故稱之為循環(huán)卷積，記為由于所以即循環(huán)卷積亦滿足交換率。2025/1/2421作為習(xí)題請(qǐng)讀者證明頻域卷積定理如果x(n)=x1(n)·x2(n)則3-40或2025/1/24223.2.4復(fù)共軛序列的DFT設(shè)x*(n)是x(n)的復(fù)共軛序列，長(zhǎng)度為NX(k)=DFT[x(n)]則

DFT[x*(n)]=X*(N-k)，0≤k≤N-1且

X(N)=X(0)證明：根據(jù)DFT的唯一性，只要證明(3-41)式右邊等于左邊即可。3-412025/1/2423證又由X(k)的隱含周期性有X(N)=X(0)用同樣的方法可以證明

DFT[x*(N-n)]=X*(k)3-422025/1/24243.2.5DFT的共軛對(duì)稱性

(1)如果x(n)=xr(n)+jxi(n)由DFT的線性性質(zhì)即可得3-493-50其中X(k)的共軛對(duì)稱分量X(k)的共軛反對(duì)稱分量(2)如果x(n)=xep(n)+xop(n)3-513-52其中2025/1/2425綜上所述：如果序列x(n)的DFT為X(k)，則x(n)的實(shí)部和虛部（包括j）的DFT分別為X(k)的共軛對(duì)稱分量和共軛反對(duì)稱分量；而x(n)的共軛對(duì)稱分量和共軛反對(duì)稱分量的DFT分別為X(k)的實(shí)部和虛部乘以j。設(shè)x(n)是長(zhǎng)度為N的實(shí)序列，且X(k)=DFT[x(n)]，則(1)X(k)共軛對(duì)稱，即

X(k)=X*(N-k)，0≤k≤N-1(2)如果x(n)=x(N-n)，則X(k)實(shí)偶對(duì)稱，即

X(k)=X(N-k)(3)如果x(n)=-x(N-n)，則X(k)純虛奇對(duì)稱，即

X(k)=-X(N-k)3-533-543-552025/1/24263.3頻率域采樣設(shè)任意序列x(n)存在Z變換且X(z)的收斂域包含單位圓(即x(n)存在傅里葉變換)。在單位圓上對(duì)X(z)等間隔采樣N點(diǎn)得到3-56顯然，式(3-56)表示在區(qū)間[0，2π]上對(duì)x(n)的傅里葉變換X(ejω)的N點(diǎn)等間隔采樣。如果將X(k)看成長(zhǎng)度為N的有限長(zhǎng)序列xN(n)的DFT，即

xN(n)=IDFT[X(k)]，0≤k≤N-12025/1/2427下面推導(dǎo)序列xN(n)與原序列x(n)之間的關(guān)系，并導(dǎo)出頻域采樣定理。

由DFT和DFS的關(guān)系可知，X(k)是xN(n)以N為周期的周期延拓序列的離散傅立葉系數(shù)的主值序列，即將式(3-56)代入上式得2025/1/2428式中

所以3-573-58式(3-58)說(shuō)明X(z)在單位圓上的N點(diǎn)等間隔采樣X(jué)(k)的IDFT為原序列x(n)以N為周期的周期延拓序列的主值序列。2025/1/2429頻域采樣定理

如果序列x(n)的長(zhǎng)度為M，則只有當(dāng)頻域采樣點(diǎn)數(shù)N≥M時(shí)，才有

xN(n)=IDFT[X(k)]=x(n)

即可由頻域采樣X(jué)(k)恢復(fù)原序列，否則產(chǎn)生時(shí)域混疊現(xiàn)象。這就是所謂的頻域采樣定理。下面推導(dǎo)用頻域采樣X(jué)(k)表示X(z)的內(nèi)插公式和內(nèi)插函數(shù)。設(shè)序列x(n)的長(zhǎng)度為M，在頻域0到之間等間隔采樣N點(diǎn)，N≥M，則有2025/1/2430式中

將上式代入X(z)的表達(dá)式中得上式中，因此3-592025/1/2431令

則3-603-61式(3-61)稱為用X(k)表示X(z)的內(nèi)插公式，φk(z)稱為內(nèi)插函數(shù)。當(dāng)z=ejω時(shí)，(3-60)式和(3-61)式就成為x(n)的傅里葉變換變換X(ejω)的內(nèi)插函數(shù)和內(nèi)插公式，即2025/1/2432進(jìn)一步化簡(jiǎn)可得

在數(shù)字濾波器的結(jié)構(gòu)與設(shè)計(jì)中，我們將會(huì)看到，頻域采樣理論及有關(guān)公式可提供一種有用的濾波器結(jié)構(gòu)和濾波器設(shè)計(jì)途徑。3-623-632025/1/24333.4DFT的應(yīng)用舉例DFT的快速算法FFT的出現(xiàn)，使DFT在數(shù)字通信、語(yǔ)音信號(hào)處理、圖象處理、功率譜估計(jì)、仿真、系統(tǒng)分析、雷達(dá)理論、光學(xué)、醫(yī)學(xué)、地震以及數(shù)值分析等各個(gè)領(lǐng)域都得到廣泛應(yīng)用。然而，各種應(yīng)用一般都以卷積和相關(guān)運(yùn)算的具體處理為依據(jù)，或者以DFT作為連續(xù)傅里葉變換的近似為基礎(chǔ)。本節(jié)主要介紹用DFT計(jì)算卷積和相關(guān)系數(shù)的基礎(chǔ)原理以及用DFT對(duì)連續(xù)信號(hào)和序列進(jìn)行譜分析等最基本的應(yīng)用。只要掌握了這兩種基本應(yīng)用的原理，就為用DFT解決數(shù)字濾波和系統(tǒng)分析等問(wèn)題打下了基礎(chǔ)。2025/1/24343.4.1用DFT計(jì)算線性卷積如果且

則由時(shí)域循環(huán)卷積定理有0≤k≤L-1由此可見(jiàn)，循環(huán)卷積既可在時(shí)域直接計(jì)算，也可以按照?qǐng)D3-10所示的計(jì)算框圖，在頻域計(jì)算。由于DFT有快速算法FFT，當(dāng)N很大時(shí)，在頻域計(jì)算的速度快得多，因而常用DFT(FFT)計(jì)算循環(huán)卷積。2025/1/2435圖3-10用DFT計(jì)算循環(huán)卷積

在實(shí)際應(yīng)用中，為了分析時(shí)域離散線性非時(shí)變系統(tǒng)或者對(duì)序列進(jìn)行濾波處理等，需要計(jì)算兩個(gè)序列的線性卷積。和計(jì)算循環(huán)卷積一樣，為了提高運(yùn)算速度，也希望用DFT(FFT)計(jì)算線性卷積。而DFT只能用來(lái)計(jì)算循環(huán)卷積，為此導(dǎo)出線性卷積和循環(huán)卷積之間的關(guān)系以及循環(huán)卷積和線性卷積相等的條件。2025/1/2436假設(shè)h(n)和x(n)都是有限長(zhǎng)序列，長(zhǎng)度分別是N和M。它們的線性卷積和循環(huán)卷積分別表示如下：其中，L≥max[N，M]，，所以3-643-652025/1/2437對(duì)照式(3-64)可以看出，上式中

即3-66式(3-66)說(shuō)明，yc(n)等于yl(n)以L為周期的周期延拓序列的主值序列。我們知道，yl(n)長(zhǎng)度為N+M-1，因此只有當(dāng)循環(huán)卷積長(zhǎng)度L≥N+M-1時(shí)，yl(n)以L為周期進(jìn)行周期延拓才混疊現(xiàn)象。此時(shí)取主值序列顯然滿足yc(n)=yl(n)。由此證明循環(huán)卷積等于線性卷積的條件是L≥N+M-1。圖3-11中畫出了h(n)、x(n)、h(n)*x(n)和L分別取6、8、10時(shí)h(n)×x(n)的波形。由于h(n)長(zhǎng)度N=4，x(n)長(zhǎng)度M=4，N+M-1=8，所以只有L≥8時(shí)，h(n)×x(n)的波形才與h(n)*x(n)相同。2025/1/2438圖3-11線性卷積與循環(huán)卷積

2025/1/2439如果取L=N+M-1，則可用DFT(FFT)計(jì)算線性卷積，計(jì)算框圖如圖3-12所示。其中DFT和IDFT通常用快速算法(FFT)來(lái)實(shí)現(xiàn)，故常稱其為快速卷積。

圖3-12用DFT計(jì)算線性卷積框圖

2025/1/2440實(shí)際上，經(jīng)常遇到兩個(gè)序列的長(zhǎng)度相差很大的情況，例如M>>N時(shí)，直接套用上述方法是不行的。解決這個(gè)問(wèn)題的方法是將長(zhǎng)序列分段處理計(jì)算。這種分段處理法有重疊相加法和重疊保留法兩種。下面只介紹重疊相加法。設(shè)序列h(n)長(zhǎng)度為，x(n)為無(wú)限長(zhǎng)序列。將均勻分段，每段長(zhǎng)度取，則式中xk(n)=x(n)·RM(n-kM)2025/1/2441于是，h(n)與x(n)的線性卷積可表示為式中

3-67式(3-67)說(shuō)明，計(jì)算h(n)與x(n)的線性卷積時(shí)，可先進(jìn)行分段線性卷積yk(n)=h(n)*xk(n)，然后再把分段卷積結(jié)果疊加起來(lái)即可。如圖3-13所示。每一分段卷積yk(n)的長(zhǎng)度為N+M-1，因此yk(n)與yk+1(n)有N-1個(gè)點(diǎn)重疊，必須把重疊部分的yk(n)與yk+1(n)相加，才能得到完整的卷積序列y(n)。2025/1/24西安建筑科技大學(xué)信于與控制學(xué)42圖3-13重疊相加法卷積示意圖

2025/1/24433.4.2用DFT對(duì)信號(hào)進(jìn)行譜分析所謂信號(hào)的譜分析就是計(jì)算信號(hào)的傅里葉變換。連續(xù)信號(hào)與系統(tǒng)的傅里葉分析顯然不便于直接用計(jì)算機(jī)進(jìn)行計(jì)算，使其應(yīng)用受到限制，而DFT是一種時(shí)域和頻域均離散化的變換，適合數(shù)值運(yùn)算，成分分析離散信號(hào)和系統(tǒng)的有力工具。1.用DFT對(duì)連續(xù)信號(hào)進(jìn)行譜分析工程實(shí)際中，經(jīng)常遇到的連續(xù)信號(hào)xa(t)，其頻譜函數(shù)Xa(jΩ)也是連續(xù)函數(shù)。2025/1/2444設(shè)連續(xù)信號(hào)xa(t)持續(xù)時(shí)間和Tp，最高頻率為fc，如圖3-14所示。xa(t)的傅里葉變換為以間隔T≤1/2fc(即fs=1/T≥2fc)采樣xa(t)得=xa(nT)。設(shè)共采樣N點(diǎn)，并對(duì)Xa(jf)作零階近似(t=nT，dt=T)得顯然，Xa(jf)仍是f的連續(xù)周期函數(shù)，和如圖3-14(b)所示。對(duì)在區(qū)間等間隔采樣N點(diǎn)，采樣間隔為F，如圖3-14(c)所示。2025/1/2445參數(shù)fc、Tp、N和F滿足如下關(guān)系式：

由于NT=Tp，所以3-683-69將f=kF和(3-68)式代入中可得Xa(jf)的采樣令Xa(k)=，x(n)=xa(nT)，則3-702025/1/2446用同樣方法，由

可以推出3-71(3-70)式說(shuō)明，連續(xù)信號(hào)的頻譜特性可以通過(guò)對(duì)連續(xù)信號(hào)采樣并進(jìn)行DFT再乘以T的近似方法得到。時(shí)域采樣信號(hào)可由(3-71)式得出。對(duì)持續(xù)時(shí)間有限的帶限信號(hào)，在滿足時(shí)域采樣定理時(shí)，上述分析方法不丟失信息。2025/1/2447但直接由分析結(jié)果Xa(k)看不到Xa(jf)的全部頻譜特性，而只能看到N個(gè)離散采樣點(diǎn)的譜特性，這就是所謂的柵欄效應(yīng)。如果的持續(xù)時(shí)間無(wú)限長(zhǎng)，上述分析中要進(jìn)行截?cái)嗵幚?，所以?huì)產(chǎn)生頻率混疊和泄露現(xiàn)象，從而使譜分析產(chǎn)生誤差。下面將討論上述問(wèn)題產(chǎn)生的原因及改進(jìn)措施。2025/1/2448圖3-14用DFT計(jì)算連續(xù)信號(hào)頻譜原理2025/1/2449理想低能濾波器的單位沖擊響應(yīng)ha(t)及其頻響函數(shù)Ha(if)如圖3-15(a)、(b)所示。圖中現(xiàn)在用DFT來(lái)分析ha(t)的頻率響應(yīng)特性。由于ha(t)的持續(xù)時(shí)間為無(wú)窮長(zhǎng)，故要截取一段Tp，假設(shè)Tp=8s，采樣間隔T=0.25s(即采樣頻率fs=4Hz)，采樣點(diǎn)數(shù)N=Tp/T=32。此時(shí)頻域采樣間隔F=1/NT=0.125Hz。則

H(k)=T·DFT[h(n)],0≤k≤31其中h(n)=ha(nT)R32(n)3-712025/1/2450圖3-15用DFT計(jì)算理想低通濾波器頻響曲線2025/1/2451H(k)的曲線如圖3-15(c)所示。由圖可見(jiàn)，低頻部分近似理想低通頻響特性，而高頻誤差較大，且整個(gè)頻響都有波動(dòng)。這些差別就是由于對(duì)ha(t)截?cái)嗨a(chǎn)生的。為減少這種截?cái)嗾`差，可適當(dāng)加長(zhǎng)Tp，增加采樣點(diǎn)數(shù)N或用窗函數(shù)處理后在進(jìn)行DFT。有關(guān)窗函數(shù)的內(nèi)容將在FIR數(shù)字濾波器設(shè)計(jì)中詳細(xì)敘述。2025/1/2452連續(xù)信號(hào)進(jìn)行譜分析對(duì)信號(hào)進(jìn)行譜分析主要關(guān)心兩個(gè)問(wèn)題，就是譜分析范圍和頻率分辨率。譜分析范圍受采樣速率fs的限制。為了不產(chǎn)生頻率混疊失真，通常要求信號(hào)的最高頻率fs<fs/2。頻率分辨率用采樣間隔F來(lái)描述，表示譜分析中能夠分辨的兩個(gè)頻譜分量的最小間隔。顯然F越小譜分析就越接近Ha(jf)，所以F較小時(shí)，我們稱頻率分辨率較高。下面討論用DFT對(duì)連續(xù)信號(hào)譜分析的參數(shù)選擇原則。2025/1/2453信號(hào)譜分析的參數(shù)選擇原則

在已知信號(hào)最高頻率fc（即譜分析范圍）時(shí)，為了避免在DFT運(yùn)算中發(fā)生頻率混疊現(xiàn)象，要求采樣速率fs滿足下式：

fs>2fc

按照式(3-69)，譜分辨率F=fs/N，如果保持采樣點(diǎn)數(shù)N不變，要提高譜的分辨率，必須降低采樣速率，采樣速率的降低會(huì)引起譜分析范圍減少；如維持fs不變，為提高譜的分辨率可以增加采樣點(diǎn)數(shù)N，因?yàn)镹T=Tp，T=1/fs-1，只有增加對(duì)信號(hào)的觀察時(shí)間，才能增加N，Tp和N可按照下式選擇：3-723-733-742025/1/2454例3-1對(duì)實(shí)信號(hào)進(jìn)行譜分析，要求譜分辨率F≤10Hz，信號(hào)最高頻率fc=2.5kHz，試確定最小記錄時(shí)間TPmin，最大的采樣間隔Tmax，最少的采樣點(diǎn)數(shù)Nmin。如果fc不變，要求譜分辨率增加一倍，最少的采樣點(diǎn)數(shù)和最小的記錄時(shí)間是多少?解：因此TPmin=0.1s，因?yàn)橐骹s≥2fc，所以2025/1/24552.用DFT對(duì)序列進(jìn)行譜分析已知道單位圓上的Z變換就是序列傅里葉變換，即為使頻率分辨率提高一倍，F(xiàn)=5Hz，要求X(ejω)是ω的連續(xù)周期函數(shù)。如果對(duì)序列x(n)進(jìn)行N點(diǎn)DFT，得到X(k)，X(k)是在區(qū)間[0，2π]上對(duì)X(ejω)的N點(diǎn)等間隔采樣。因此序列的傅里葉變換可利用DFT(FFT)來(lái)計(jì)算。2025/1/2456對(duì)周期為N的周期序列，其頻譜函數(shù)為其中

由DFT的隱含周期性知道，截取的主值序列x(n)=RN(n)，并進(jìn)行N點(diǎn)DFT得到所以可用X(k)表示的頻譜結(jié)構(gòu)。2025/1/2457如果截取長(zhǎng)度M等于的整數(shù)個(gè)周期，即M=mN，m為正整數(shù)，則令n=n’+rN，r=0，1，...，m-1，n’=0，1，...，N-1，則2025/1/2458由此可見(jiàn)，XM(k)也能表示的頻譜結(jié)構(gòu)，只是在k=rm時(shí)，XM(rm)=，表示的r次諧波譜線，其幅度擴(kuò)大m倍，而其它k值時(shí)，XM(k)=0。當(dāng)然，X(r)與XM(rm)對(duì)應(yīng)點(diǎn)頻率是相等的。所以，只要截取的整數(shù)個(gè)周期進(jìn)行DFT，就可得到它的頻譜結(jié)構(gòu)，達(dá)到譜分析的目的。因?yàn)樗?025/1/2459如果的周期預(yù)先不知道，可先截取M點(diǎn)進(jìn)行DFT，即

再將截取長(zhǎng)度擴(kuò)大一倍，截取比較XM(k)和X2M(k)，如果二者的主譜差別滿足分析誤差要求，則已以XM(k)或X2M(k)近似表示的頻譜，否則，繼續(xù)將截取長(zhǎng)度加倍，直至前后兩次分析所得主譜頻率差別滿足誤差要求。設(shè)最后截取長(zhǎng)度為rM，則XrM(k0)表示點(diǎn)的譜線強(qiáng)度。2025/1/2460實(shí)際應(yīng)用中，并非整個(gè)單位圓上的頻譜都很有意義。

例如，對(duì)于窄帶信號(hào)，往往只希望對(duì)信號(hào)所在的一段頻段進(jìn)行譜分析，這時(shí)便希望采樣能密集地在著段頻帶內(nèi)進(jìn)行，而帶外部分可完全不管。有時(shí)希望采樣點(diǎn)不局限于單位圓上。例如，語(yǔ)音信號(hào)處理中，常常需要知道系統(tǒng)極點(diǎn)所對(duì)應(yīng)的頻率，如果極點(diǎn)位置偏離單位圓較遠(yuǎn)，則其單位圓上的頻譜就很平滑，如圖3-16(a)所示，這時(shí)很難從中識(shí)別出極點(diǎn)對(duì)應(yīng)的頻率。如果使采樣點(diǎn)軌跡接近這些極點(diǎn)的弧線或圓周進(jìn)行，則采樣結(jié)果將會(huì)在極點(diǎn)對(duì)應(yīng)的頻率上出現(xiàn)明顯的尖峰，如圖3-16(b)所示，這樣就能準(zhǔn)確地測(cè)定出極點(diǎn)頻率。2025/1/2461圖3-16單位圓與非單位圓采樣2025/1/2462對(duì)均勻分布在以原點(diǎn)為圓心的任何圓上的N點(diǎn)頻譜采樣，可用DFT（FFT）計(jì)算，而沿螺旋弧線采樣，則要用線性調(diào)頻Z變換（Chirp-Z變換，簡(jiǎn)稱CZT）計(jì)算。例如，要求計(jì)算序列在半徑為r的圓上的頻譜，那么N個(gè)等間隔采樣點(diǎn)為k=0,1,2,…,N-1，zk點(diǎn)的頻譜分量為2025/1/2463令，則上式表明，要計(jì)算x(n)在半徑為r的圓上的等間隔頻譜分量，可以先對(duì)x(n)乘以r-n，再計(jì)算N點(diǎn)DFT(FFT)即可得到。若要求x(n)分布在該圓的有限角度內(nèi)N點(diǎn)等間隔頻譜分量，可以通過(guò)尾部補(bǔ)零的方法，仍按式(3-75)用DFT分析整個(gè)圓上的等間隔頻譜，最后只取所需角度內(nèi)的頻譜分量即可。顯然這種方法的計(jì)算量大，下面要介紹的Chirp-Z變換可使這種譜分析的運(yùn)算量大大減少。3-752025/1