盧文良

(海軍701工廠 研發(fā)部,北京　100015)

?

兩種針對(duì)小信號(hào)的脈壓優(yōu)化算法及實(shí)現(xiàn)

盧文良

(海軍701工廠 研發(fā)部,北京100015)

摘要結(jié)合某雷達(dá)項(xiàng)目,針對(duì)項(xiàng)目LFM和Taylor編碼兩種小時(shí)寬帶寬積(TB)信號(hào)脈沖壓縮效果不理想的問題,通過對(duì)兩種不同小TB信號(hào)的實(shí)際回波波形和特性進(jìn)行了分析,針對(duì)性地使用了倒推法和循環(huán)迭代法,獲得了優(yōu)化的脈壓系數(shù)。試驗(yàn)得出結(jié)論,較之傳統(tǒng)的匹配濾波脈沖壓縮系數(shù),采用優(yōu)化脈壓系數(shù)的LFM和Taylor編碼小TB信號(hào)脈沖壓縮主副瓣比均能提高十幾dB。

關(guān)鍵詞脈沖壓縮;倒推法;循環(huán)迭代法

Two Small Signal Pulse Compression Algorithms and Their Implementation

LU Wenliang

(R & D Department,Navy 701 Factory,Beijing 100015,China)

AbstractBased on a radar project,in which the result of two small TB signal (the LFM and Taylor) pulse compression methods are not good enough,the article firstly analyses the shape and characteristics of the two small TB signal.The optimum compressed coefficient is obtained by the reverse algorithm and the circulation algorithm respectively.Experiments show that the main lobe to sidelobe ratio can be improved by more than 10 dB compared with that by the traditional compression algorithms.

Keywordspulse compression;reverse algorithm;circulation algorithm

眾所周知,脈沖壓縮體制的雷達(dá)可在發(fā)射時(shí)采用寬脈沖保證足夠的作用距離,接收時(shí)通過脈沖壓縮獲得窄脈沖,將目標(biāo)信息從雜波背景中提取出來,因而能較好地解決作用距離和分辨能力之間的矛盾[1]。隨著高速大規(guī)模集成電路的發(fā)展,數(shù)字脈沖壓縮系統(tǒng)受到了廣泛的重視并普遍應(yīng)用于各種雷達(dá)設(shè)計(jì)中。它利用一組匹配濾波脈壓系數(shù)模擬脈沖壓縮系統(tǒng),并通過對(duì)脈壓系數(shù)進(jìn)行一系列變換來改進(jìn)脈沖壓縮的性能,提高脈沖壓縮的主副瓣比?？蓪?duì)于小時(shí)寬帶寬積信號(hào)而言,即便進(jìn)行加窗脈沖壓縮,得到的主副瓣比仍較小,小目標(biāo)不易被檢測(cè)出來。且在實(shí)際處理過程中,由于回波與發(fā)射信號(hào)不能夠完全匹配、雷達(dá)噪聲影響較大等因素,進(jìn)一步降低了小時(shí)寬帶寬積信號(hào)脈沖壓縮的主副瓣比,加大了雷達(dá)檢測(cè)小目標(biāo)的困難程度。因此常常會(huì)出現(xiàn)通過理論計(jì)算和計(jì)算機(jī)仿真設(shè)定的脈沖壓縮技術(shù)指標(biāo)與實(shí)際應(yīng)用結(jié)果相差甚遠(yuǎn),不容易滿足要求的情況出現(xiàn)。為解決這個(gè)問題,提高雷達(dá)檢測(cè)性能,人們進(jìn)行了諸多努力[2-6]。

結(jié)合某雷達(dá)項(xiàng)目,針對(duì)信號(hào)處理中LFM和Taylor編碼小時(shí)寬帶寬積信號(hào)脈沖壓縮效果不理想的問題,本文提出了使用倒推法和循環(huán)迭代法分別對(duì)兩種小時(shí)寬帶寬積信號(hào)產(chǎn)生優(yōu)化的脈沖壓縮系數(shù),以提高脈沖壓縮的主副瓣比。通過大量試驗(yàn),表明倒推法對(duì)于LFM信號(hào)效果較好,循環(huán)迭代法則對(duì)Taylor編碼信號(hào)效果顯著。本文介紹了算法的硬件實(shí)現(xiàn)和使用方法,即采用FPGA+DSP芯片組合方式的信號(hào)處理板,具有并行、高速、大容量、高可靠性的優(yōu)點(diǎn)。由FPGA芯片進(jìn)行下變頻獲得基帶信號(hào)作為算法輸入,根據(jù)相應(yīng)算法得到優(yōu)化的脈壓系數(shù),再將脈壓系數(shù)存儲(chǔ)到DSP芯片中進(jìn)行脈壓運(yùn)算。最后將兩種脈壓算法的脈壓結(jié)果表示出來,大量試驗(yàn)表明這兩種算法較傳統(tǒng)匹配濾波脈壓算法主副瓣比性能有大幅提高。目前上述兩種算法已成功得到應(yīng)用。

1背景介紹

理想情況下的脈沖壓縮系統(tǒng)是一個(gè)匹配濾波系統(tǒng),其要求其相頻特性與發(fā)射信號(hào)實(shí)現(xiàn)相位共軛匹配[7]。

假設(shè)輸入信號(hào)為ui(t),其頻域表達(dá)式為式(1)所示

(1)

若使用傳統(tǒng)數(shù)字頻域脈壓,脈壓系數(shù)應(yīng)為輸入信號(hào)頻譜的共軛轉(zhuǎn)置,即如式(2)所示

(2)

通常會(huì)對(duì)脈壓系數(shù)在時(shí)域或頻域加窗來提高脈壓性能。以LFM小時(shí)寬帶寬積信號(hào)為例,加窗后脈沖壓縮的主副瓣比較之加窗前有較大提高,仿真證明,時(shí)域加窗結(jié)果比頻域加窗結(jié)果效果更為顯著。在此選擇時(shí)域加窗,則時(shí)域加窗可表示為

h(t)=ui(t)×hamming(t)

(3)

hamming(t)為漢明窗函數(shù)。

將h(t)變換到頻域可表示為

(4)

頻域加窗系數(shù)為輸入信號(hào)頻譜的共軛轉(zhuǎn)置,即

(5)

頻域脈壓輸出為

Uo(ω)=Ui(ω)Hhamming(ω)

(6)

可得匹配濾波脈沖壓縮結(jié)果為

uo(t)=ifft[Uo(ω)]

(7)

由于匹配濾波脈沖壓縮將信號(hào)能量集中于一點(diǎn),而隨機(jī)噪聲能量維持不變,所以能提高輸出信號(hào)信噪比,將信號(hào)從噪聲中分離出來。脈沖壓縮輸出信號(hào)的信噪比與輸入信號(hào)信噪比之比提高了D倍,D等于時(shí)寬帶寬積,即

(8)

由此可看出,在相同輸入信號(hào)條件下,大時(shí)寬帶寬積信號(hào)較之小時(shí)寬帶寬積信號(hào)有更好的脈壓性能,更易被檢測(cè)出來。

必須指出在實(shí)際實(shí)現(xiàn)時(shí)通常不可能得到完全匹配,迫使系統(tǒng)工作在一定程度的“失配”狀態(tài)下,所得脈沖壓縮性能要遠(yuǎn)低于式(8)的結(jié)論。因此對(duì)于小時(shí)寬帶寬積信號(hào)來說,如何提高脈壓性能,達(dá)到實(shí)際需求就顯得尤為重要。

本項(xiàng)目中用到兩種信號(hào)形式:LFM信號(hào)和Taylor四相編碼信號(hào)。

LFM信號(hào)的數(shù)學(xué)表達(dá)式[8]為

(9)

其中,τ為信號(hào)時(shí)寬;μ為信號(hào)的調(diào)頻斜率;ω0為信號(hào)中心頻率;A為信號(hào)幅度。

Taylor編碼信號(hào)的一般表達(dá)式為

(10)

其中,vk代表第k位Taylor編碼,其含義為

vk=js(k-1)Wk

(11)

式中,s是固定值,值為+1或-1;Wk表示第k位二相碼,是+1(0°)或-1(180°)。

式(10)中,p(t)是周期為4T的半余弦波形

(12)

這里T為子脈沖寬度。

2改進(jìn)方法

2.1倒推法

倒推法通過提取匹配濾波脈沖壓縮輸出信號(hào)的相位和幅度信息,合成理想的脈沖壓縮結(jié)果,并對(duì)脈沖壓縮過程做逆變換,獲得優(yōu)化的脈沖壓縮系數(shù)。

由前面脈沖壓縮的基本原理和推導(dǎo)公式可知,脈沖壓縮的輸出可表示如式(7)所示,其最大值可以表示為

(13)

表征脈沖壓縮性能的主要指標(biāo)主副瓣比通過下式表征

(14)

為獲得理想脈沖壓縮結(jié)果的理想波形,首先選擇輸出信號(hào)主副瓣比波形中信號(hào)主瓣以外的區(qū)域,乘以一個(gè)主副瓣比降低因子N,將副瓣和噪聲區(qū)域的分貝值大幅降低,如下式所示

sdbdream(t)=sdb×N

(15)

結(jié)合式(14)和式(15)推導(dǎo)得到uabs(t)的表達(dá)式為

(16)

脈沖壓縮輸出信號(hào)uo(t)的相位可表示為

α(t)=angle(uo(t))

(17)

則理想脈沖壓縮的輸出信號(hào)udream(t)可表示為

udream(t)=uabs(t)ejα(t)/π

(18)

將其變換到頻域,理想脈沖壓縮結(jié)果為

Udream(ω)=fft(udream(t))

(19)

可得脈壓系數(shù)的頻域表達(dá)形式為

(20)

使用式(20)得到的理想頻域脈壓系數(shù)Hdream(ω)進(jìn)行脈沖壓縮可提高輸出信號(hào)的主副瓣比十幾dB。通過改變N的大小,可相應(yīng)改善脈沖壓縮主副瓣比的大小,但N達(dá)到一定值后,改善作用便不能提高,這與實(shí)際信號(hào)的相關(guān)性和信號(hào)純凈度有關(guān)。

2.2循環(huán)迭代法

循環(huán)迭代法以脈沖壓縮輸出信號(hào)的積分副瓣電平為指標(biāo),通過對(duì)脈沖壓縮過程的不斷循環(huán)迭代,最終獲得積分副瓣電平最小時(shí)的脈沖壓縮系數(shù)。

假設(shè)長(zhǎng)度為N的Taylor信號(hào)第k個(gè)采樣點(diǎn)可表示為sk,k=0…N-1。將信號(hào){sk}拓展,得到長(zhǎng)度為P(P≥N)的脈沖壓縮輸入信號(hào)

S=sT=[0…0s0s1…sN-10…0]

(21)

假設(shè)脈沖壓縮數(shù)字濾波器系數(shù)可表示為

hT=[h0h1…h(huán)N-1hN…h(huán)p-1]

(22)

信號(hào)通過脈沖壓縮數(shù)字濾波器后輸出可表示為

y=hHX

(23)

(24)

為增大脈沖壓縮輸出信號(hào)的主副瓣比,采用積分副電平(ISL)作為循環(huán)收斂的指標(biāo),如下式所示

ISL=yHWy

(25)

W=diag(w1…wp-1wp…w2p-1)=diag(1…101…1)

(26)

同時(shí)滿足在主瓣方向的輸出為一常數(shù),即

(27)

求解,得

h=B-1S(SHBS)-1N

(28)

其中,B=XWXH。由于(SHBS)-1N為一常數(shù),從而

h=B-1S

(29)

式(27)是針對(duì)積分副瓣電平而言,實(shí)際希望輸出峰值副瓣電平最小,這可調(diào)整式(26)中權(quán)矢量的構(gòu)成,對(duì)高副瓣賦以大的權(quán)值。調(diào)整公式為

(30)

通過調(diào)整加權(quán)系數(shù)矩陣W以及脈沖壓縮數(shù)字濾波器階數(shù)P,采用循環(huán)迭代算法便可得到較為滿意的濾波器設(shè)計(jì)結(jié)果,其過程如圖1所示。

圖1　循環(huán)迭代法流程圖

3硬件實(shí)現(xiàn)

上述兩種脈壓系數(shù)產(chǎn)生算法由包含F(xiàn)PGA+DSP芯片組合的信號(hào)處理板來實(shí)現(xiàn),FPGA采用Altera公司的CycloneⅡ系列的EP2C484C6芯片[9],完成下變頻功能,DSP采用AD公司的TS201芯片,進(jìn)行脈沖壓縮運(yùn)算,FPGA和DSP之間采用Link口通信。整個(gè)系統(tǒng)具有大存儲(chǔ)量、大運(yùn)算量和并行處理的能力,支持高速I/O和級(jí)聯(lián)特性,其處理板結(jié)構(gòu)圖如圖2所示。算法實(shí)現(xiàn)方法如下:

首先信號(hào)處理板接收頻綜產(chǎn)生的中頻輸入信號(hào),通過A/D變換后,進(jìn)入FPGA芯片實(shí)現(xiàn)下變頻,將中頻信號(hào)轉(zhuǎn)換為基帶信號(hào),并經(jīng)過抽取,通過Link口進(jìn)入DSP芯片,DSP芯片將采到的基帶信號(hào)導(dǎo)出,作為上述兩種算法的輸入信號(hào)。計(jì)算機(jī)使用Matlab對(duì)兩種算法進(jìn)行編程,并使用基帶信號(hào)產(chǎn)生我們所需要的優(yōu)化脈壓系數(shù)。將產(chǎn)生的脈壓系數(shù)燒入DSP芯片內(nèi)作為脈壓系數(shù)存儲(chǔ),重新運(yùn)行脈沖壓縮程序,此時(shí)信號(hào)處理板就能夠得到更好的脈沖壓縮主副瓣比,脈壓結(jié)果可導(dǎo)出至Matlab進(jìn)行觀察和比較。

圖2信號(hào)處理板結(jié)構(gòu)圖

4實(shí)驗(yàn)結(jié)果

圖3是采集到的解調(diào)后實(shí)際LFM小時(shí)寬帶寬積回波信號(hào)波形,時(shí)寬為20μs,帶寬2MHz,采樣頻率8MHz。圖4是采集到的解調(diào)后Taylor編碼小時(shí)寬帶寬積回波實(shí)、虛部信號(hào)波形,其子脈沖時(shí)寬為0.5μs,碼元個(gè)數(shù)為40,采樣頻率為4MHz?？梢钥吹?LFM和Taylor信號(hào)的波形均已發(fā)生形變,與計(jì)算機(jī)仿真波形相差甚遠(yuǎn)。

圖3　LFM小信號(hào)實(shí)際波形圖

圖4　Taylor小信號(hào)波形圖

下面是與匹配濾波加權(quán)(Hamming)法進(jìn)行比較的結(jié)果圖,由于傳統(tǒng)的不加權(quán)匹配濾波法脈壓結(jié)果只有約10dB,本文為了比較結(jié)果清晰,不將不加權(quán)匹配濾波結(jié)果表示在圖上。

圖5表示對(duì)LFM小時(shí)寬帶寬積實(shí)際信號(hào)分別用匹配濾波加權(quán)法和倒推法進(jìn)行脈壓的結(jié)果,可以看到,匹配濾波加權(quán)算法的脈壓主副瓣比為22dB,在主副瓣比增加因子N=4的情況下,倒推法的脈壓主副瓣比為36dB,較匹配濾波提高了14dB。

圖5　 LFM信號(hào)匹配濾波加權(quán)法與倒推法脈壓結(jié)果的比較

圖6表示對(duì)Taylor小時(shí)寬帶寬積實(shí)際信號(hào)分別用匹配濾波加權(quán)法和循環(huán)迭代法進(jìn)行脈壓的結(jié)果,可以看到,匹配濾波加權(quán)的脈壓主副瓣比為17dB,在脈沖壓縮系數(shù)長(zhǎng)度P=300的情況下,倒推法的脈壓主副瓣比為34dB,較匹配濾波提高了17dB。

圖6　Taylor信號(hào)匹配濾波加權(quán)法與循環(huán)迭代法脈壓結(jié)果的比較

從實(shí)驗(yàn)結(jié)果可知,這兩種算法較傳統(tǒng)匹配濾波較算法對(duì)于提高小時(shí)寬帶寬積信號(hào)的主副瓣比有較大提高,但是這兩種方法會(huì)受到不同信號(hào)質(zhì)量的影響,主副瓣比的提高有限。

5結(jié)束語

在雷達(dá)信號(hào)處理中,對(duì)小時(shí)寬帶寬積信號(hào)進(jìn)行匹配濾波脈沖壓縮的結(jié)果主副瓣比通常不夠理想,本文針對(duì)LFM和Taylor兩種小時(shí)寬帶寬積信號(hào)采用了不同的優(yōu)化算法,試驗(yàn)證明這兩種算法可以將脈沖壓縮的主副瓣比提高十幾dB,大幅提高了雷達(dá)的探測(cè)能力。除了算法原理,本文還介紹了算法實(shí)現(xiàn)的硬件平臺(tái),即一個(gè)由FPGA和DSP芯片組成的高速、大容量、高速I/O信號(hào)處理板,能滿足優(yōu)化脈壓系數(shù)的產(chǎn)生和應(yīng)用。無論是對(duì)于倒推法或循環(huán)迭代法而言,基帶信號(hào)的純凈度和參數(shù)的選擇對(duì)脈沖壓縮的性能有較大影響,需要考慮到波形完整性、噪聲大小、運(yùn)算量、主瓣寬度和旁瓣長(zhǎng)度等因素。

參考文獻(xiàn)

[1]MerrillISkolnik.Radarhandbook[M].3版.北京:電子工業(yè)出版社,2010.

[2]陶廣源,廖桂生,劉宏偉.多相碼信號(hào)數(shù)字脈壓濾波器設(shè)計(jì)[J].電波科學(xué)學(xué)報(bào),2003,18(2):143-146,152.

[3]王萱芬,李巍,王子立.一種實(shí)現(xiàn)超低旁瓣的脈沖壓縮算法[J].現(xiàn)代電子技術(shù),2009,32(23):5-7.

[4]王麗萍,蘇濤.一種迭代加權(quán)最小二乘旁瓣抑制濾波器設(shè)計(jì)[J].火控雷達(dá)技術(shù),2008,37(1):92-96.

[5]徐景明,朱燦焰.一種復(fù)合調(diào)制信號(hào)的脈沖壓縮及DSP實(shí)現(xiàn)[J].雷達(dá)科學(xué)與技術(shù),2010,8(5):448-451.

[6]張宏偉,任新濤,湯宮民,等.編碼信號(hào)的脈沖壓縮技術(shù)及實(shí)現(xiàn)[J].計(jì)算機(jī)測(cè)量與控制,2012,20(12):3348-3350.

[7]丁鷺飛,耿富錄.雷達(dá)原理[M].西安:西安電子科技大學(xué)出版社,2001.

[8]向敬成,張明友.雷達(dá)系統(tǒng)[M].北京:電子工業(yè)出版社,2001.

[9]AlteraCorporation.CycloneⅡdevicehand-book[M].MAUSA:AlteraCorporation,2005.

中圖分類號(hào)TN957.51

文獻(xiàn)標(biāo)識(shí)碼A

文章編號(hào)1007-7820(2016)03-061-04

doi:10.16180/j.cnki.issn1007-7820.2016.03.015

作者簡(jiǎn)介:盧文良(1981—),男,工程師。研究方向:雷達(dá)信號(hào)處理。

收稿日期:2015- 07- 23