陳 旭,王金鋒,朱 平

(中國船舶集團(tuán)有限公司第七二三研究所,江蘇 揚(yáng)州 225101)

0 引 言

在雷達(dá)研制、調(diào)試和試驗(yàn)的各個(gè)階段,都需要對(duì)雷達(dá)性能進(jìn)行測試,尤其是雷達(dá)的主要戰(zhàn)術(shù)指標(biāo)如探測威力、探測精度等。而傳統(tǒng)的外場試驗(yàn)成本高,周期長,試驗(yàn)效率低,因此需要采用現(xiàn)代化的便捷測試方法,能夠在試驗(yàn)室或近場對(duì)雷達(dá)研制全周期進(jìn)行性能考核評(píng)估,提高雷達(dá)的研制效率。其中以電磁藍(lán)軍模擬為核心需求的雷達(dá)目標(biāo)回波模擬器,成為當(dāng)前目標(biāo)模擬仿真領(lǐng)域的一個(gè)研究熱點(diǎn)。

本文以高分辨率雷達(dá)探測需求為背景,根據(jù)雷達(dá)發(fā)射信號(hào)建立滿足要求的高分辨率擴(kuò)展目標(biāo)模型,給出一種基于目標(biāo)布點(diǎn)方式的擴(kuò)展目標(biāo)模擬方法,并以數(shù)字儲(chǔ)頻單元(DRFM)為核心模塊研制了一套雷達(dá)目標(biāo)回波模擬系統(tǒng)[1],最后采用中頻數(shù)據(jù)采集、逆合成孔徑雷達(dá)(ISAR)成像處理的方法對(duì)模擬的目標(biāo)數(shù)據(jù)進(jìn)行分析測試,測試結(jié)果滿足高分辨率雷達(dá)目標(biāo)模擬的要求。

1 擴(kuò)展目標(biāo)模型

常規(guī)的目標(biāo)模擬方法以點(diǎn)目標(biāo)為主,僅考慮目標(biāo)的雷達(dá)截面(RCS)特性,這一特性由雷達(dá)電磁波波長、極化方式和目標(biāo)的姿態(tài)角等參數(shù)決定。而對(duì)于高分辨率雷達(dá)來說,這種點(diǎn)目標(biāo)模擬方式不能全面、合理地模擬艦船、飛機(jī)等復(fù)雜目標(biāo)。這類目標(biāo)的散射特性可表現(xiàn)為多個(gè)強(qiáng)散射點(diǎn)。本文采用擴(kuò)展目標(biāo)模型進(jìn)行描述并模擬,分別計(jì)算各個(gè)散射點(diǎn)中心目標(biāo)特性與雷達(dá)發(fā)射信號(hào)的時(shí)域卷積,得到單個(gè)散射點(diǎn)的回波信號(hào),再將各散射點(diǎn)的回波信號(hào)進(jìn)行時(shí)域疊加,得到擴(kuò)展目標(biāo)的回波信號(hào)[2-3]。

設(shè)雷達(dá)發(fā)射信號(hào)為線性調(diào)頻信號(hào),其基帶信號(hào)用s(t)表示為:

(1)

式中:TP為發(fā)射脈沖寬度;Kr為調(diào)頻斜率;t為采樣時(shí)間。

則其中任一個(gè)散射點(diǎn)在ta時(shí)刻的回波信號(hào)即為該點(diǎn)在雷達(dá)發(fā)射信號(hào)基礎(chǔ)上增加時(shí)延、幅度調(diào)制和多普勒調(diào)制,用下式表示為:

(2)

式中:σ為該散射點(diǎn)的幅度值;R(ta)為該散射點(diǎn)與雷達(dá)間距離;λ為電磁波波長;c為電磁波傳播速度,為一常數(shù)。

擴(kuò)展目標(biāo)回波可以看成多個(gè)散射點(diǎn)的集合,即表示為各散射點(diǎn)的矢量和,所以擴(kuò)展目標(biāo)回波模型可用下式表示:

(3)

式中:M為散射點(diǎn)數(shù)目;σi為第i個(gè)點(diǎn)的幅度值;Ri(ta)為第i個(gè)散射點(diǎn)與雷達(dá)間距離;其它參數(shù)定義同式(2)。

2 擴(kuò)展目標(biāo)模擬方法

本文設(shè)計(jì)并模擬實(shí)現(xiàn)的擴(kuò)展目標(biāo)包括一維距離像目標(biāo)、二維面目標(biāo)和三維體目標(biāo)。根據(jù)上一節(jié)所建立的擴(kuò)展目標(biāo)模型,以三維體目標(biāo)為例,詳細(xì)描述其具體的模擬實(shí)現(xiàn)方法。

三維體目標(biāo)即模擬的擴(kuò)展目標(biāo)各散射點(diǎn)坐標(biāo)值包括x、y、z三維信息,設(shè)置雷達(dá)坐標(biāo)為(xr,yr,zr),目標(biāo)的基準(zhǔn)點(diǎn)在雷達(dá)坐標(biāo)系下的位置為(x0,y0,z0),目標(biāo)的某個(gè)散射點(diǎn)P在以基準(zhǔn)點(diǎn)為原點(diǎn)的目標(biāo)坐標(biāo)系下的位置為(xi,yi,zi)。設(shè)目標(biāo)在(x0-xr,y0-yr,z0-zr)方向上平動(dòng)的速度為V,與Z軸方向的夾角為θ,與X軸方向的夾角為φ;目標(biāo)在ta時(shí)刻的姿態(tài)信息即橫滾、俯仰和偏航表示為(γ,?,ψ),如圖1所示。

圖1 三維體目標(biāo)運(yùn)動(dòng)示意圖

(1) 計(jì)算t時(shí)刻由目標(biāo)姿態(tài)變化引起的位置變化:

(4)

(2) 計(jì)算t時(shí)刻由目標(biāo)平動(dòng)所引起的位置變化:

(5)

式中:(V,θ,φ)根據(jù)目標(biāo)運(yùn)動(dòng)態(tài)勢得到,根據(jù)系統(tǒng)需求定時(shí)更新。

(3) 計(jì)算散射點(diǎn)P與雷達(dá)的瞬時(shí)斜距為:

(6)

(4) 計(jì)算t時(shí)刻散射點(diǎn)P的回波信號(hào),即式(2)。

(5) 計(jì)算t時(shí)刻整個(gè)擴(kuò)展目標(biāo)的回波信號(hào),即式(3)。

3 基于擴(kuò)展目標(biāo)的雷達(dá)目標(biāo)回波模擬系統(tǒng)設(shè)計(jì)

根據(jù)上節(jié)建立的擴(kuò)展目標(biāo)模型以及目標(biāo)模擬方法,本文基于數(shù)字儲(chǔ)頻法(數(shù)字射頻存儲(chǔ)器(DRFM)產(chǎn)生法)構(gòu)建技術(shù)方案:采用射頻接收機(jī)接收雷達(dá)信號(hào);利用下變頻電路將雷達(dá)射頻信號(hào)變?yōu)楣β氏鄬?duì)恒定的中頻信號(hào),然后利用高速模數(shù)轉(zhuǎn)換器(ADC)對(duì)雷達(dá)中頻信號(hào)進(jìn)行高速采樣、存儲(chǔ),并在數(shù)字邏輯電路的控制下,將存儲(chǔ)的數(shù)字中頻基帶信號(hào)經(jīng)延遲、多普勒調(diào)制、疊加等,利用高速數(shù)模轉(zhuǎn)換器(DAC)產(chǎn)生滿足要求的擴(kuò)展目標(biāo)中頻信號(hào);然后經(jīng)過上變頻、功率放大后輸出模擬的射頻信號(hào)[4]。實(shí)現(xiàn)原理如圖2所示。

圖2 基于擴(kuò)展目標(biāo)的雷達(dá)目標(biāo)回波模擬系統(tǒng)原理框圖

數(shù)字儲(chǔ)頻單元為本系統(tǒng)的核心部分,實(shí)現(xiàn)對(duì)雷達(dá)信號(hào)的采集存儲(chǔ),以及對(duì)目標(biāo)回波的重構(gòu)以產(chǎn)生所需要的擴(kuò)展目標(biāo)。該單元基于高速信號(hào)數(shù)字采樣和實(shí)時(shí)存儲(chǔ)技術(shù),采用現(xiàn)場可編程門陣列(FPGA)和數(shù)字信號(hào)處理器(DSP)作為核心處理架構(gòu),并由高速AD采樣器、DA轉(zhuǎn)換器、大容量存儲(chǔ)器及控制接口等組成。硬件組成見圖3。

圖3 數(shù)字儲(chǔ)頻單元硬件總體框圖

如圖3所示,數(shù)字儲(chǔ)頻單元以1片Virtex-6系列 FPGA芯片 XC6VSX315T和1片DSP芯片TMS320C6455作為核心處理架構(gòu)。其中DSP芯片實(shí)現(xiàn)對(duì)外部數(shù)據(jù)包括系統(tǒng)配置參數(shù)、雷達(dá)參數(shù)等的解析,生成擴(kuò)展目標(biāo)調(diào)制參數(shù),包括本振參數(shù)、速度、距離、功率等,并發(fā)送給FPGA芯片;FPGA芯片對(duì)采樣后的雷達(dá)中頻數(shù)據(jù)進(jìn)行抽取,然后根據(jù)DSP發(fā)送過來的參數(shù)完成回波信號(hào)的延時(shí)、擴(kuò)展點(diǎn)卷積、多普勒調(diào)制、功率調(diào)制、多個(gè)散射點(diǎn)目標(biāo)疊加、插值、上混頻、幅度補(bǔ)償后輸出。此外,數(shù)字儲(chǔ)頻單元其它主要電路還包括1片高速ADC,分辨率10 bit,最高采樣率可達(dá)3 GSPS,用于對(duì)輸入的雷達(dá)信號(hào)進(jìn)行數(shù)字采樣;1片DAC轉(zhuǎn)換器,最高頻率2.5 GSPS,用于將產(chǎn)生的擴(kuò)展目標(biāo)信號(hào)轉(zhuǎn)換成模擬信號(hào)輸出;1片Virtex-5系列 FPGA芯片 XC5VLX30T用于外部接口控制、板內(nèi)數(shù)據(jù)傳輸?shù)?。?shù)字儲(chǔ)頻單元信號(hào)處理流程如圖4所示。

圖4 數(shù)字儲(chǔ)頻單元信號(hào)處理流程

4 應(yīng)用舉例及驗(yàn)證

本節(jié)舉例說明對(duì)研制的擴(kuò)展目標(biāo)模擬系統(tǒng)的測試方法和測試結(jié)果。測試框圖見圖5。

圖5 系統(tǒng)測試框圖

如圖5所示,由雷達(dá)信號(hào)模擬器產(chǎn)生寬帶雷達(dá)線性調(diào)頻信號(hào);控制計(jì)算機(jī)產(chǎn)生需要模擬的擴(kuò)展目標(biāo)的各種參數(shù),包括所有散射點(diǎn)位置、散射點(diǎn)功率、目標(biāo)運(yùn)動(dòng)速度、目標(biāo)轉(zhuǎn)動(dòng)速度等;目標(biāo)數(shù)據(jù)采集系統(tǒng)則實(shí)現(xiàn)對(duì)擴(kuò)展目標(biāo)中頻數(shù)據(jù)的采集,以進(jìn)行后續(xù)數(shù)據(jù)分析。

本文采用基于擴(kuò)展目標(biāo)布點(diǎn)的方法模擬目標(biāo),然后對(duì)錄取的數(shù)據(jù)進(jìn)行逆合成孔徑雷達(dá)(ISAR)成像處理及檢驗(yàn)。因此首先由控制計(jì)算機(jī)產(chǎn)生擴(kuò)展目標(biāo)的參數(shù),包括質(zhì)心位置、各散射點(diǎn)位置、擴(kuò)展目標(biāo)的運(yùn)動(dòng)速度、目標(biāo)轉(zhuǎn)速等。圖6為產(chǎn)生的一個(gè)三維擴(kuò)展目標(biāo)0°時(shí)的散射點(diǎn)布置模型。

圖6 擴(kuò)展目標(biāo)布點(diǎn)設(shè)計(jì)示意圖

擴(kuò)展目標(biāo)模擬系統(tǒng)工作參數(shù)設(shè)置:

(1) 擴(kuò)展目標(biāo)參數(shù)設(shè)置:散射點(diǎn)數(shù)量21個(gè),散射點(diǎn)位置相對(duì)基準(zhǔn)點(diǎn)坐標(biāo)為(-9,-18,1.5)、(-6,-18,1.5)、(-3,-18,1.2)、(0,-18,1.2)、(3,-18,1.2)、(6,-18,1.0)、(9,-18,1.0)、(0,-15,1.0)、(0,-12,0.6)、(0,-9,0.6)、(0,-6,0.6)、(0,-3,0.3)、(0,0,0)、(-6,1.5,-0.3)、(6,1.5,-0.3)、(-3,3,-0.6)、(0,3,-1)、(3,3,-1)、(0,6,-1)、(0,9,-1.2)、(0,12,-1.2),目標(biāo)運(yùn)動(dòng)速度V=20 m/s,目標(biāo)轉(zhuǎn)動(dòng)速度A=30°/s。

(2) 雷達(dá)信號(hào)參數(shù):采用線性調(diào)頻信號(hào),信號(hào)帶寬Br=400 MHz,脈寬Tp=100 μs,重頻fPRF=5 kHz;

(3) 目標(biāo)模擬系統(tǒng)AD采樣率fs=1 200 MHz。

對(duì)擴(kuò)展目標(biāo)模擬系統(tǒng)產(chǎn)生的目標(biāo)數(shù)據(jù),錄取后進(jìn)行ISAR成像處理[5]。處理流程如圖7所示。

圖7 錄取數(shù)據(jù)ISAR成像處理流程

(1) 去斜處理

對(duì)錄取的擴(kuò)展目標(biāo)回波數(shù)據(jù)首先進(jìn)行去斜處理,針對(duì)上節(jié)所建立的目標(biāo)回波模型,設(shè)置去斜參考信號(hào)為:

(7)

式中:Rref為參考距離,即為雷達(dá)到目標(biāo)質(zhì)心的距離,其它參數(shù)定義同式(3)。

將回波信號(hào)與參考信號(hào)的共軛相乘,并進(jìn)行傅里葉變換,得到下式:

(8)

(2) 一維距離成像

對(duì)去斜后的回波信號(hào)進(jìn)行快速傅里葉變換可得到一維距離像,如圖8所示。

(3) 包絡(luò)對(duì)齊

包絡(luò)對(duì)齊的目的是消除由于目標(biāo)運(yùn)動(dòng)相對(duì)于雷達(dá)的平動(dòng)造成的相鄰回波在距離向上的錯(cuò)位,通過對(duì)齊調(diào)整,把積累時(shí)間內(nèi)的回波信號(hào)在距離向上調(diào)整為一致,作為后面初相校正和方位壓縮的基礎(chǔ)。本文采用相鄰相關(guān)法實(shí)現(xiàn)包絡(luò)對(duì)齊。

(4) 相位補(bǔ)償

在進(jìn)行包絡(luò)對(duì)齊后,各次回波中還存在由于目標(biāo)平動(dòng)所造成的相位變化,因此需要通過相位補(bǔ)償來進(jìn)一步校正。常用的相位補(bǔ)償法包括特顯點(diǎn)相位補(bǔ)償法、參數(shù)估計(jì)相位補(bǔ)償法和基于圖像準(zhǔn)則的相位補(bǔ)償法等。本文采用計(jì)算量相對(duì)小、易于工程實(shí)現(xiàn)的特顯點(diǎn)法。

(5) 成像處理

在完成包絡(luò)對(duì)齊和相位補(bǔ)償后,可進(jìn)行成像處理。本次錄取數(shù)據(jù)方位向取256個(gè)脈沖,成像處理結(jié)果如圖9所示。

對(duì)比圖6設(shè)計(jì)模型、目標(biāo)轉(zhuǎn)動(dòng)角度和圖9成像結(jié)果,本文所設(shè)計(jì)的雷達(dá)目標(biāo)回波模擬系統(tǒng)滿足要求,能夠模擬典型復(fù)雜擴(kuò)展目標(biāo)。

5 結(jié)束語

本文以三維體目標(biāo)為例,詳細(xì)介紹了基于目標(biāo)布點(diǎn)方式的擴(kuò)展目標(biāo)模擬方法,即將擴(kuò)展目標(biāo)分解為多個(gè)強(qiáng)散射點(diǎn),分別計(jì)算各散射點(diǎn)的回波信號(hào),再將各散射點(diǎn)的回波信號(hào)進(jìn)行時(shí)域疊加,得到擴(kuò)展目標(biāo)的回波信號(hào)。通過推導(dǎo)得到了具體的計(jì)算公式,然后基于數(shù)字儲(chǔ)頻原理研制了雷達(dá)目標(biāo)回波模擬系統(tǒng),并對(duì)模擬的目標(biāo)采用ISAR成像處理進(jìn)行測試驗(yàn)證。測試結(jié)果表明:研制的雷達(dá)目標(biāo)回波模擬系統(tǒng)能夠有效模擬復(fù)雜擴(kuò)展目標(biāo),為高分辨率雷達(dá)調(diào)試、試驗(yàn)提供了一種行之有效的手段。