肖 雷，周秀珍

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

0 引 言

隨著現(xiàn)代軍事工業(yè)的迅速發(fā)展,雷達(dá)高度計(jì)在導(dǎo)彈的制導(dǎo)過程中起著重要作用，現(xiàn)在的雷達(dá)高度計(jì)主要分為脈沖雷達(dá)高度計(jì)和調(diào)頻連續(xù)波(FMCW)雷達(dá)高度計(jì)兩大類[1]。根據(jù)雷達(dá)發(fā)射波束角的寬度，脈沖雷達(dá)高度計(jì)可以劃分為寬波束與窄波束2種。在平臺(tái)穩(wěn)定、指向不變的情況下，窄波束雷達(dá)相比較寬波束雷達(dá)來說，測(cè)量精度更高，因此在實(shí)際應(yīng)用中，窄波束雷達(dá)高度計(jì)相較寬波束雷達(dá)高度計(jì)而言，其應(yīng)用更加廣泛。但是由于彈載雷達(dá)平臺(tái)不穩(wěn)定，在這種條件下，窄波束雷達(dá)高度計(jì)不能保證波束能夠準(zhǔn)確指向雷達(dá)天底點(diǎn)，因而在彈載雷達(dá)領(lǐng)域，須采用寬波束雷達(dá)高度計(jì)。

與點(diǎn)目標(biāo)信號(hào)的模擬不同，寬波束雷達(dá)高度計(jì)的目標(biāo)模型為地面。在此背景下，介紹了3種雷達(dá)測(cè)高算法，并設(shè)計(jì)仿真流程完成了對(duì)3種雷達(dá)測(cè)高算法的仿真實(shí)現(xiàn)，對(duì)仿真生成的回波進(jìn)行測(cè)高，將測(cè)高結(jié)果與高程圖的真實(shí)高度進(jìn)行對(duì)比，最終對(duì)3種算法的測(cè)高精度進(jìn)行了比較。

1 脈沖雷達(dá)回波仿真模型

回波功率仿真和相干視頻回波仿真是現(xiàn)代雷達(dá)回波仿真中比較常見的2種方法，回波功率仿真只關(guān)注功率大小,而相干視頻回波仿真主要關(guān)注的是信號(hào)的相位信息。

1.1 相干視頻回波模型

如果雷達(dá)的發(fā)射信號(hào)表示為：

ψT(t)=VT(t)ej2πfct

(1)

式中：VT(t)表示為矩形包絡(luò)。

而ψT(t)又可以表示為：

(2)

式中：Tt為發(fā)射信號(hào)脈寬。

如果從雷達(dá)發(fā)射信號(hào)到接收回波的這段時(shí)間里天線指向、地面散射特性、位置沒有變化，且不忽略頻率的多普勒變化，則由雷達(dá)方程可推出點(diǎn)目標(biāo)的回波方程為：

(3)

(4)

γ和目標(biāo)散射截面積σ的關(guān)系可表示為：

(5)

將式(1)、式(4)、式(5)代入式(3)中，可得到雷達(dá)的單個(gè)散射單元回波信號(hào)為：

(6)

由于大部分自然地形都是粗糙的表面，散射元回波相位并不十分精確，需要在式(6)基礎(chǔ)上加上一個(gè)相位的隨機(jī)分量φ，φ的概率密度為：

(7)

最后得到的散射單元相干視頻回波表達(dá)式如下：

(8)

1.2 功率疊加回波模型

由雷達(dá)方程可推出單個(gè)散射目標(biāo)的回波功率如下：

(9)

經(jīng)研究[2]，在粗糙地面的假設(shè)條件下，雷達(dá)地面回波總功率的期望值等于照射區(qū)域內(nèi)各個(gè)散射單元回波功率的期望值之和：

(10)

地面回波的功率模型可由以上這個(gè)原理計(jì)算得出。

2 雷達(dá)測(cè)高算法

2.1 回波前沿雷達(dá)測(cè)高半功率點(diǎn)算法

回波前沿的雷達(dá)測(cè)高半功率點(diǎn)算法是把雷達(dá)天底點(diǎn)位置反射的回波和回波前沿的半功率點(diǎn)位置相對(duì)應(yīng)。由此可以通過半功率點(diǎn)相對(duì)于回波上升沿起點(diǎn)的時(shí)間延時(shí)τH來解算出天底點(diǎn)的高度：

(11)

式中：c為光速。

2.2 偏差重心雷達(dá)測(cè)高跟蹤算法

偏差重心(OCOG)雷達(dá)測(cè)高跟蹤算法的基本思路是首先須定義一個(gè)矩形跟蹤窗口，然后用一個(gè)與窗內(nèi)回波重心相同、面積相同的矩形和窗內(nèi)的回波進(jìn)行擬合，通過調(diào)整跟蹤窗使波形中心和矩形中心重合，最后矩形寬度的一半加上跟蹤窗中心就是跟蹤點(diǎn)。假設(shè)矩形的幅度為A，矩形的寬度為W，E為誤差量，由此得出它們的計(jì)算公式：

(12)

(13)

(14)

(15)

通過誤差量E來調(diào)整跟蹤窗的位置，直至誤差等于零，這時(shí)跟蹤窗的中心位置與W/2之和就是跟蹤天底點(diǎn)的時(shí)延。

2.3 平均高度雷達(dá)測(cè)高算法

平均高度雷達(dá)測(cè)高算法的基本思路是：某時(shí)刻的回波功率可以通過一定空間距離之內(nèi)所有散射元在這一時(shí)刻回波功率的疊加來表示，通過計(jì)算雷達(dá)的平均距離和這些散射元就可以得到平均高度。

可計(jì)算得出平均高度：

(16)

式中：PR(i)代表回波功率的采樣。

3 3種雷達(dá)測(cè)高算法實(shí)現(xiàn)

對(duì)上面介紹的3種雷達(dá)測(cè)高算法,下面分別給出了其實(shí)現(xiàn)的流程圖,如圖1、圖2、圖3所示。

圖1 半功率點(diǎn)雷達(dá)測(cè)高算法流程圖

圖2 OCOG雷達(dá)測(cè)高算法流程圖

4 仿真結(jié)果及分析

選取出50個(gè)雷達(dá)的航跡點(diǎn)，分別通過上面介紹的3種雷達(dá)測(cè)高算法可以計(jì)算出高度，之后再與真實(shí)的地形數(shù)據(jù)進(jìn)行比較，計(jì)算出3種雷達(dá)測(cè)高算法的誤差，仿真結(jié)果如圖4所示。

圖3 平均高度雷達(dá)測(cè)高算法流程圖

圖4 3種雷達(dá)測(cè)高算法誤差高度

由圖4可以觀察出寬波束雷達(dá)高度計(jì)在復(fù)雜地形的回波，其中半功率點(diǎn)雷達(dá)測(cè)高算法會(huì)引起比較大的誤差，而平均高度雷達(dá)測(cè)高算法和偏差重心雷達(dá)測(cè)高跟蹤算法能更好地估計(jì)出雷達(dá)高度，其中平均高度雷達(dá)測(cè)高算法的估計(jì)精確度最好。