脈沖信號(hào)是現(xiàn)代雷達(dá)主要采用的信號(hào)形式，脈沖信號(hào)頻率測(cè)量是雷達(dá)偵察中不可或缺的環(huán)節(jié)，對(duì)雷達(dá)對(duì)抗起著重要的作用。數(shù)字化處理是雷達(dá)對(duì)抗系統(tǒng)發(fā)展的趨勢(shì)之一，常用的數(shù)字測(cè)頻方法包括過零點(diǎn)檢測(cè)法、相位差分法、快速傅里葉變換( FFT)法和現(xiàn)代譜估計(jì)法。其中FFT法工程可實(shí)現(xiàn)性強(qiáng)，實(shí)時(shí)性好，且適用于寬帶偵收，因此在工程中得到廣泛應(yīng)用。

　　本文以時(shí)寬較短( 0. 2～1μs)的正弦波脈沖信號(hào)為研究對(duì)象，分析了傳統(tǒng)FFT測(cè)頻法的不足之處，從工程應(yīng)用角度分析了提高測(cè)頻精度的改進(jìn)方法，并提出了基于FPGA的全數(shù)字實(shí)現(xiàn)流程。

　　1 FFT測(cè)頻

　　信號(hào)x( t)經(jīng)過數(shù)字化采樣后為x( n)，n = 0，1，2，…，N-1，為對(duì)其進(jìn)行頻譜分析，進(jìn)行離散傅里葉變換( DFT)，將信號(hào)從時(shí)域轉(zhuǎn)換到頻域，如式( 2)所示：

　　可見，F(xiàn)FT測(cè)頻的頻率分辨率只與信號(hào)時(shí)寬有關(guān)，根據(jù)譜線的最大值來換算信號(hào)的頻率，如果信號(hào)的頻率正好落在一根譜線上，得到的頻率測(cè)量結(jié)果是準(zhǔn)確的，而在多數(shù)情況下，信號(hào)頻率落在兩根譜線之間，由最大值譜線位置反映的頻率不再準(zhǔn)確，最大測(cè)頻誤差為Δf /2.

　　脈沖是雷達(dá)最常采用的信號(hào)形式，根據(jù)需要，雷達(dá)有時(shí)會(huì)采用脈內(nèi)帶調(diào)制的信號(hào)類型，例如相位編碼、線性調(diào)頻等，對(duì)于此類復(fù)雜信號(hào)可采用各種信號(hào)處理方法將其轉(zhuǎn)化為普通正弦波信號(hào)，因此正弦波脈沖的測(cè)頻方法具有通用性。根據(jù)上文分析結(jié)果，對(duì)于時(shí)寬較長的脈沖，采用FFT測(cè)頻法易于實(shí)現(xiàn)較高測(cè)頻精度，滿足設(shè)備指標(biāo)要求。但是對(duì)于短脈沖，例如一個(gè)0. 2μs寬的脈沖，根據(jù)式( 3)，理論能達(dá)到的測(cè)頻精度只有2. 5 MHz，難以滿足偵察要求。

　　2補(bǔ)零技術(shù)

　　補(bǔ)零是指在進(jìn)行FFT運(yùn)算之前在時(shí)域數(shù)據(jù)的尾部添加一些零，并使總的時(shí)域數(shù)據(jù)點(diǎn)數(shù)保持為2的冪次方。由于補(bǔ)零不增加任何新的信息，所以并不改變頻譜形狀和頻率分辨率，補(bǔ)零只是在原始點(diǎn)數(shù)的FFT結(jié)果中內(nèi)插了一些頻率分量。對(duì)于點(diǎn)數(shù)較少的FFT結(jié)果，在大多數(shù)情況下，從中找到峰值比較困難，也很難觀察到頻譜的細(xì)微結(jié)構(gòu)。而補(bǔ)零之后，功率譜的峰值位置可以較清晰的顯露出來，有助于提高對(duì)主瓣峰值頻率分量進(jìn)行精確定位的能力，由此提高測(cè)頻精度。

　　補(bǔ)零技術(shù)的缺點(diǎn)是額外增加了處理量，補(bǔ)零越多，處理時(shí)間也就越長。此外，對(duì)于存在噪聲的情況，補(bǔ)零也不能改善信噪比，存在頻譜峰值點(diǎn)定位錯(cuò)誤的可能，造成測(cè)頻誤差增大。

　　3插值FFT測(cè)頻方法分析

　　3. 1插值FFT頻率估計(jì)原理

　　插值FFT估計(jì)頻率方法利用真正的頻譜峰值兩側(cè)的2根FFT譜線，求其幅度比值，建立一個(gè)以修正頻率為變量的方程，解方程得到修正頻率值，對(duì)FFT最大譜線位置進(jìn)行校正，以實(shí)現(xiàn)對(duì)信號(hào)頻率更高精度的估計(jì)，如圖1所示。相比上節(jié)補(bǔ)零的方法，不必增加FFT的長度以及由此帶來的運(yùn)算處理量，只需從FFT結(jié)果中找出兩個(gè)點(diǎn)就足夠。

　　圖1矩形窗頻譜函數(shù)

　　在圖1中插值頻率校正即求出矩形窗譜主瓣中心與相鄰譜線的橫坐標(biāo)差，對(duì)于譜線位置x、x + 1，其矩形窗譜函數(shù)為sinc函數(shù)，表示為f( x)，頻譜值為yx、yx+1，矩形窗譜函數(shù)和頻譜值已知，可構(gòu)成一方程如下：

　　式中，α= yx /yx+1.實(shí)際應(yīng)用中，已知FFT譜峰最大值位置k1，相鄰次大值位置k2，頻率分辨率Δf，利用修正頻率值校正頻率可得：

　　當(dāng)k2 = k1 + 1時(shí)，取加號(hào); k2 = k1-1時(shí)，取減號(hào)。