呂東澤，喻旭偉，王 俐，蔡易君

（中國電子科技集團(tuán)公司29所，成都610036）

0 引 言

數(shù)字射頻存儲器（DRFM）可以產(chǎn)生與雷達(dá)信號相干的干擾信號，對采用多普勒、脈沖壓縮等技術(shù)的相干雷達(dá)具有很好的干擾效果，因此DRFM在電子對抗領(lǐng)域得到了廣泛的關(guān)注和應(yīng)用［1］。DRFM根據(jù)其量化方式分為幅度量化和相位量化2種。其中幅度量化DRFM具有輸出信噪比高、無模糊帶寬較大等優(yōu)點(diǎn)，但是其在對干擾信號進(jìn)行相位調(diào)制時較為復(fù)雜；相位量化DRFM具有動態(tài)范圍大、實(shí)現(xiàn)結(jié)構(gòu)簡單等優(yōu)點(diǎn)，但是當(dāng)存在多信號時，會出現(xiàn)小信號抑制現(xiàn)象［2］。

為了使DRFM可根據(jù)實(shí)際情況靈活選取量化方式，本文提出了幅度和相位雙模量化的DRFM系統(tǒng)。但是傳統(tǒng)DRFM采用的模擬正交變換方法的鏡像抑制比最高只能達(dá)到30dB［3－4］。為了簡化結(jié)構(gòu)、提升性能，本文在雙模DRFM中采用數(shù)字正交變換方法形成兩正交支路。該DRFM系統(tǒng)具有結(jié)構(gòu)簡單、應(yīng)用靈活、鏡像抑制比高等優(yōu)點(diǎn)。

1 正交雙模DRFM結(jié)構(gòu)

圖1 數(shù)字正交雙模DRFM系統(tǒng)結(jié)構(gòu)

基于數(shù)字正交變換方法的雙模DRFM系統(tǒng)的結(jié)構(gòu)如圖1所示。該DRFM只需要采用1個模擬／數(shù)字（A／D）轉(zhuǎn)換器、1個現(xiàn)場可編程門陣列（FPGA）和1個數(shù)字／模擬（D／A）轉(zhuǎn)換器，利用FPGA實(shí)現(xiàn)數(shù)字正交變換、量化、雷達(dá)信號存儲、干擾信號調(diào)制、混頻等功能。該DRFM系統(tǒng)相比傳統(tǒng)DRFM系統(tǒng)具有結(jié)構(gòu)簡單、功能豐富、易于實(shí)現(xiàn)等優(yōu)點(diǎn)。

2 數(shù)字正交變換法

DRFM接收到的雷達(dá)帶通信號表示為：

式 中： I（t）＝ a（t）cos［φ（t）］；Q（t）＝a（t）sin［φ（t）］。

利用采樣頻率fs＝4f0／（2m－1）且fs≥2B對s（t）進(jìn)行采樣，其中m表示正整數(shù)，B表示信號帶寬，可得采樣序列為：

式 中： I（t）＝ a（t）cos［φ（n）］；Q（t）＝a（t）sin［φ（n）］。

常用的數(shù)字正交變換方法有Hilbert變換法、低通濾波法和多相濾波法，本文分別對其進(jìn)行分析。

2.1 基于Hilbert濾波法的DRFM

Hilbert變換的理想頻譜為［5］：

由上式可知，Hilbert變換為90°移相器，通過該變換后信號只有相位變換，幅度沒有變化。DRFM接收到的雷達(dá)實(shí)信號經(jīng)過Hilbert變換后變?yōu)榻馕鲂盘?，變換后頻譜為：

一般利用數(shù)字濾波器來實(shí)現(xiàn)Hilbert濾波器，理想Hilbert濾波器沖擊響應(yīng)為：

經(jīng)過Hilbert濾波器后，信號會產(chǎn)生（M－1）／2個采樣點(diǎn)的延遲，M為濾波器階數(shù)。由于采用偶數(shù)階Hilbert濾波器會使得延遲為非整數(shù)采樣點(diǎn)，從而增加延遲匹配的難度，因此一般Hilbert濾波器都采用奇數(shù)階數(shù)。

綜上，基于Hilbert濾波法的數(shù)字正交變換示意圖如圖2所示。將A／D轉(zhuǎn)換后的雷達(dá)信號通過Hilbert濾波器使其與原信號正交，再將另一支路進(jìn)行延時處理使2個支路時間匹配，然后對其量化。

圖2 Hilbert濾波法示意圖

2.2 基于低通濾波法的DRFM

基于低通濾波法的數(shù)字正交變換如圖3所示，將經(jīng)過A／D轉(zhuǎn)換器的雷達(dá)信號分為2路，分別利用數(shù)控振蕩器（NCO）產(chǎn)生的正交信號相乘，然后經(jīng)過低通濾波器濾除高頻分量，再對2個支路分別量化。

圖3 低通濾波法示意圖

NCO分別產(chǎn)生cos（ω0n）和sin（ω0n），與s（n）相乘（ω0＝2πf0／fs）可得［6］：

對式（6）、（7）再進(jìn)行低通濾波后，可得：

低通濾波法中由于I、Q 2個支路采用相同的低通濾波器參數(shù)，因此2個支路具有較好的一致性。

2.3 基于多相濾波法的DRFM

由式（2）中采樣后的序列可得：

對s（n）進(jìn)行奇、偶抽取和符號變換可得：

式中：I′（n）和Q′（n）分別為I（n）和Q（n）的2倍抽取，由抽取原理可知，只要I′（n）和Q′（n）數(shù)字譜小于π／2即fs≥4B，則其2倍抽取可以無失真地表示原序列［7］。它們的頻譜相差1個相位因子ejω／2，相當(dāng)于時域相差1／2個采樣時間，可以采用2個延時濾波器對其偏差進(jìn)行校正，延時校正濾波器可以采用內(nèi)插因子為4的多項(xiàng)濾波內(nèi)插器的分支濾波器［8－9］：HI1（ejω）＝1，HQ1（ejω）＝ejω／2或者 HI2（ejw）＝ejω／4，HQ2（ejω）＝ejω3／4。

因此基于多相濾波法DRFM的正交變換結(jié)構(gòu)如圖4所示，首先對經(jīng)過A／D轉(zhuǎn)換器的雷達(dá)信號進(jìn)行奇、偶抽取和符號變換獲得I、Q 2個支路，然后利用多項(xiàng)濾波器的分支濾波器實(shí)現(xiàn)不同的分?jǐn)?shù)相移從而使I、Q 2路信號時序匹配。

圖4 多相濾波法示意圖

3 仿真分析

利用Matlab對基于數(shù)字正交變換的雙模DRFM進(jìn)行仿真分析：接收到的雷達(dá)信號頻率為f0＝120MHz，幅度量化位數(shù)為8bit，相位量化位數(shù)為3bit，Hilbert濾波器、低通濾波器的階數(shù)均為33，多相濾波器的分支濾波器為8階。

處理前信號頻譜及經(jīng)過3種雙模DRFM的信號頻譜如圖5所示，圖（a）表示處理前原始信號頻譜；（b）、（c）、（d）分別表示利用3種方法經(jīng)過濾波器后，再進(jìn)行8bit幅度量化、3bit相位量化后信號頻譜。從圖5可見，處理前的信號頻譜存在干擾峰值，雙模DRFM可以有效對信號進(jìn)行幅度和相位量化，3種正交變換方法均可有效抑制干擾。

采用3種方法的DRFM隨不同頻移其鏡像抑制比（IR）的變化如圖6所示。由圖6可見，低通濾波法隨著頻移的變化其IR一直保持較好的性能；多相濾波法在頻移較小的范圍內(nèi)可以保持較好的IR性能，但隨著頻移變大，IR性能隨之大幅度下降；Hilbert法最大IR低于其他2種方法，但是相比多相濾波法可以在更廣的頻移范圍內(nèi)保持IR性能，超過范圍后隨著頻移的增大IR性能下降速度低于多相濾波法。

4 結(jié)束語

本文首先提出了雙模體制的DRFM系統(tǒng)，該DRFM結(jié)構(gòu)僅需1個A／D轉(zhuǎn)換器、1個FPGA和1個D／A轉(zhuǎn)換器，可以有效對雷達(dá)信號進(jìn)行幅度、相位量化，其結(jié)構(gòu)簡單且功能豐富。雙模DRFM中采用數(shù)字正交變換方法相比模擬方法可以大幅度提升DRFM的鏡像抑制比，其中Hilbert法可達(dá)到的最優(yōu)IR在3種方法中最低，超出一定頻移范圍后多相濾波法隨頻移變化性能下降最快，低通濾波法鏡像抑制性能隨頻移變化較小。

圖5 雷達(dá)信號及3種方法雙模DRFM處理后頻譜

圖6 鏡像抑制比隨頻移變化

［1］Schleher D C.Electronic Warefare in The Information Age［M］.Boston：Artech House，2000.

［2］李菊.基于幅相量化的DRFM 技術(shù)［J］.上海航天，2008（5）：57－60.

［3］邱兆坤.基于Hilbert變換的正交檢波器設(shè)計(jì)新方法［J］.系統(tǒng)工程與電子技術(shù)，2010，10（32）：2045－2048.

［4］曾濤，龍騰，毛二可.一種新的中頻采樣濾波器的設(shè)計(jì)與實(shí)現(xiàn)［J］.電子學(xué)報，2000，28（6）：51－54.

［5］謝小娟，楊凌云.數(shù)字正交分解算法的研究［J］.山西電子技術(shù)，2007（6）：83－87.

［6］徐新.數(shù)字接收機(jī)中I／Q解調(diào)算法的研究［D］.成都：電子科技大學(xué)，2005.

［7］藍(lán)金巧，孫曉聞，吳順君.基于多相濾波器的數(shù)字正交檢波濾波器的選擇［J］.計(jì)算機(jī)工程與應(yīng)用，2004（18）：215－224.

［8］楊小牛，樓才義，徐建良.軟件無線電原理與應(yīng)用［M］.北京：電子工業(yè)出版社，2001.

［9］Pellon L E.A double Nyquist digital product dector for quadrature sampling［J］.IEEE Transactions on Signal Processing，1992，40（7）：1670－1681.