基于WDFT的頻率偏移估計

1、基于WDFT的頻率偏移估計摘要：WDFT提供了在不增加采樣數(shù)目N的情況下，對任意選定的頻譜區(qū)域增加頻譜精度?？梢詫⒋幚硇盘柕闹攸c頻段的頻譜精度大幅提高，又能保持信號非重點頻段的頻譜精度的基本要求。該文重點研究了利用WDFT算法來估計高頻載波信號的一個比較小的頻率偏移。在估計精度和計數(shù)復(fù)雜程度上,與傳統(tǒng)的DFT和NDFT進(jìn)行了比較，突出了WDFT的優(yōu)越性，并用數(shù)例充分說明了這個比較。 關(guān)鍵詞：DFT; NDFT; WDFT; Computational complexity Estimation of frequency offset using warped discrete fourier

2、 transformAbstract: WDFT provided without increasing the number of samples in the case of N, for any selected frequency spectrum to increase the precision region. Signal processing can be the focus of a substantial increase in precision frequency spectrum, while preserving the signal frequencies of

3、the spectrum of non-focus accuracy of the basic requirements. This paper focuses on the use of WDFT algorithm to estimate the high-frequency carrier signal of a relatively small frequency offset. Count in the estimation accuracy and complexity, and the traditional DFT and NDFT compared, highlighting

4、 the advantages of WDFT and fully illustrated with several cases of this comparison.1.引言隨著數(shù)字技術(shù)與計算機技術(shù)的發(fā)展，數(shù)字信號處理技術(shù)已深入到各學(xué)科領(lǐng)域?？焖俑盗⑷~變換（FFT）技術(shù)的提出，大大減少了離散傅立葉變換（DFT）的計算量，使之得到廣泛的應(yīng)用。但是DFT只能給出均勻分布點處Z變換的值，要提高頻譜精度，必須增加采樣點數(shù)目N，這就導(dǎo)致計算量迅速增加。為克服DFT算法的上述缺陷，我們期望獲得更為一般的采樣點處的Z變換值，于是提出了使用FFT來計算單位圓上非均勻分布采樣點處Z變換的算法，從而可以用不均勻

5、的頻譜精度來進(jìn)行頻譜的計算分析，這就是彎曲離散傅立葉變換( Warped Discrete Fourier Transform , WDFT )。WDFT提供了在不增加采樣數(shù)目N的情況下，對任意選定的頻譜區(qū)域增加頻譜精度的一種良好選擇。WDFT與被分析信號的頻率特性相結(jié)合，通過選擇全通彎曲函數(shù)及其彎曲參數(shù)，可以將待處理信號的重點頻段的頻譜精度大幅提高，又能保持信號非重點頻段的頻譜精度的基本要求，因而可作為數(shù)字信號處理中非常有用的工具。我們知道，信號處理領(lǐng)域的信號都和一定的頻段有關(guān)系，例如在分析語音信號時，由于人耳的生理心理模型的作用，我們更關(guān)注語音信號低頻段的情況，而較少關(guān)注其高頻段的情況。這

6、就剛好和某一類AWF及彎曲參數(shù)相吻合，因此，可以使用WDFT技術(shù)，不增加采樣點數(shù)目N就可以得到更好的語音信號分析精度。應(yīng)用WDFT的理論可以估計射頻信號的一個很小的頻率偏移,并且在估計精度和計算次數(shù)上,與DFT和NDFT相比,有很大的優(yōu)越性。在很多地方（如接收器或蜂窩系統(tǒng)的基站），都需要估計在一個大載波信號附近的頻率偏移。典型的，如信號 (1)是有限長時間信號，經(jīng)過頻率為的抽樣信號抽樣后就轉(zhuǎn)化成了一個N點離散序列，這時可表示為： ，n=0，2，,N-1 和 (2) Xn中的頻率偏移就是用WDFT的理論來估計的。 該論文的結(jié)構(gòu)如下：第一部分簡要介紹了WDFT和NDFT的理論，第二部分討論了WDF

7、T估計頻率偏移的應(yīng)用，第三部分將WDFT在估計精度和計算復(fù)雜度上與DFT、NDFT進(jìn)行了比較，第四部分做出了總結(jié)。1.1 DFT DFT在分析離散信號頻譜上的應(yīng)用很廣泛。對于給定長度為N的序列xn，其離散傅立葉變換Xk就定義為z平面單位圓上均勻分布的N個點處z變換的值。 (3)在估計頻率偏移的問題上，xn的頻率抽樣必須達(dá)到高頻，由于DFT只能提供固定精度，要估計頻率偏移，其長度必須很大，因此考慮到了采用WDFT或NDFT理論來估計。1.2 NDFT有限長度N序列x(n)的非均勻離散傅立葉變換定義為： k=0,1,N-1 (4)其中,是Z平面上任意分布的N個不同點，式(4)可以寫成如下矩陣形式：

8、 X=Dx (5) 其中： , (6) (7)注意，NDFT矩陣D是范德蒙特矩陣（Vandermonde），完全由N個點來決定。對于頻率偏移的估計，NDFT的頻率響應(yīng)中的w可以達(dá)到，因此可以估計頻偏。通常情況下，NDFT的計算涉及到了由x(n)組成的長度為N的矢量D的乘法。其抽樣點的靈活性可以使頻譜精度達(dá)到我們所期望的任何值，然而，隨著計算的深入，需要個復(fù)雜的乘法運算。在某些特殊情況下，例如在單位圓上，NDFT可以應(yīng)用Goertzel算法達(dá)到計算量的簡化。1.3 WDFTWDFT是一般非均勻離散傅立葉變換(NDFT)的特殊形式，與被分析信號的頻率特性相結(jié)合，通過選擇全通彎曲函數(shù)及其彎曲參數(shù)，全

9、通彎曲函數(shù)(AWF)將頻率坐標(biāo)彎曲了，在平面單位圓上均勻分布的點被映射到z平面單位圓上的不均勻分布點?？梢詫⒋幚硇盘柕闹攸c頻段的頻譜精度大幅提高，又能保持信號非重點頻段的頻譜精度的基本要求，它與DFT最大的不同是避免了單純靠提高采樣點數(shù)目N來提高頻譜精度的問題。為方便起見，我們記長度為N的序列xn的N點WDFT 等于采用下列變換將Xz的修正Z變換在N個均勻分布點的頻譜采樣值： (8)將映射應(yīng)用到 (9) 可得 (10) 若定義 (11) 其中是的鏡像多項式，即=，即得到 (12) WDFT就定義為在處的值 (13)2. 應(yīng)用WDFT估計頻率偏移對于Xn的頻率響應(yīng)，須在頻率附近抽樣密集，才能有

10、效的估計(2)式中的頻率偏移，的最大值頻率抽樣可以幫助我們確定附近的頻率偏移。需要選擇的參數(shù)如下：1 WDFT的長度2 全通函數(shù)的階數(shù)M3 的系數(shù)參數(shù)選定后，接下來就討論WDFT的計算理論。2.1 WDFT的計算上面已經(jīng)定義了WDFT的值， 再定義 = (14) 其中是的第i個系數(shù)，同樣地，定義 (15)和的階次都是N-1，因為 (16) 則WDFT的計算可以從式(13)而化為 (17) 令，分別為從和的系數(shù)得到的長度為N的序列的N點DFT ， (18) 則 (19) 最后，WDFT的系數(shù)可以寫成 (20) 其中 (21)接下來我們用一個例子來說明這個過程，并比較分別由WDFT,DFT,NDF

11、T計算的結(jié)果。2.2 例子 在(1)中，令，并且序列Xn的長度N為64。若使用傳統(tǒng)的DFT，其頻率精度為，因此用64點DFT來估計的頻率偏移的最大誤差為 (22) 如果我們假設(shè)，那就意味著在附近用64點DFT來估計的頻率偏移的最大誤差為 (23) 現(xiàn)在我們用2階全通函數(shù)，64點WDFT來估計頻率偏移。令，此時頻率映射為： ，a，b為適當(dāng)?shù)膶嵪禂?shù) (24) 令，可以得到 (25) 由 (26) (27) 由(25)式比較兩邊的實部和虛部，可以得到以下的頻率映射 (28) 3三種算法的比較三種算法的比較主要是從估計頻偏的精度和計算的復(fù)雜程度兩方面來對比的。3.1 精度在WDFT的頻率范圍a，b內(nèi)，

12、我們?nèi)山M值，一組為的正向偏差頻率點，一組為的負(fù)向偏差頻率，可以獲得比直接在附近取單向偏差更高的精度。對于一個64點的WDFT，當(dāng)a=2.176,b=-1.166和a=-2.087,b=-1.633,在WDFT抽樣幅度取最大的情況下,可以較精確的計算兩組值的頻率偏差。MATLAB可以很好的驗證上述結(jié)論，此時在附近的頻率偏移量的范圍為-0.05rad/sample到+0.05rad/sample。(22)式表明，64點DFT不能估計此范圍的頻偏，而64點WDFT可以估計此范圍的頻偏，且最大誤差為： (29) 若和取(23)的值，此時 =35Hz (30) 我們可以看出利用WDFT理論需要64點就

13、可以有效的估計，但是使用DFT需要1024個點。利用64點NDFT的理論來估計上述例子，在的頻率范圍內(nèi)分均勻取樣，可以得到最大誤差 (31) (32) 顯然，對于M=2，與利用WDFT的理論相比，利用NDFT的理論來估計頻偏可獲得更高的精度，但是這是以犧牲計算量來獲得的。如果使用高階彎曲函數(shù)，并且選擇適當(dāng)?shù)膹澢鷧?shù)，可以提高WDFT的精度。3.2 計算復(fù)雜度令x為N維復(fù)數(shù)輸入向量，N點的DFT要求次復(fù)數(shù)乘和次復(fù)數(shù)加，假設(shè)一次復(fù)數(shù)相乘包括四次實數(shù)乘和二次實數(shù)加，一次復(fù)數(shù)加包括兩次實數(shù)加，即等于次實數(shù)乘及次實數(shù)加。WDFT要求次實數(shù)乘和次實數(shù)加.對于N=64，則需要5120次實數(shù)乘和9344次實數(shù)

14、加，而1024點FFT需要20480次實數(shù)乘和30720次實數(shù)加。如果應(yīng)用NDFT理論，我們靈活選取抽樣點，可以獲得更高的精度，但是同時也增加了計算的復(fù)雜度。從x直接計算NDFT，要求計算x與階復(fù)矩陣相乘，即次實數(shù)乘法和次實數(shù)加。對于N=64，則需要16384次實數(shù)乘和16256次實數(shù)加。事實上NDFT是在Z平面的單位圓上進(jìn)行抽樣的，利用Goertzel算法可以減少計算量。二階Goertzel算法需要次實數(shù)乘和次實數(shù)加。對于N=64，則要求8704次實數(shù)乘和16768次實數(shù)加。 對于N=64，不同算法的運算次數(shù)總結(jié)如下表，當(dāng)N增加時，其計算次數(shù)也會增加。N=64 實數(shù)乘 實數(shù)加DFT 768 1152 WDFT 5120 9344NDFT 16384 162564 總結(jié)在這篇論文中，我們討論了應(yīng)用WDFT估計頻率偏移。選擇適當(dāng)?shù)膮?shù)，可以獲得較高的精度和較理想的計算量。特別是對于N較大時，WDFT比NDFT更有效率，計算次數(shù)少，雖然后者在選取采樣點時更隨便。盡管這篇論文的例子都是應(yīng)用二階全通翹曲函數(shù)，然而高階全通翹曲函數(shù)更方便。隨著階數(shù)的的增高，估計精度也會