船载高频地波雷达海杂波时域抑制方法与流程

本发明涉及杂波抑制领域，具体涉及船载高频地波雷达海杂波时域抑制方法。

背景技术：

对于高频地波雷达来讲，其主要作用就是检测和跟踪超视距的海面目标和低空目标，从而实现早期预警的作用。岸基高频地波雷达主要作用在于海岸监测和低空防卫，随着岸基高频地波雷达发展成熟，船载高频地波雷达继承了岸基雷达一些特点，同样也具备自身的优势，例如具备更好的机动性和生存能力。然而随着平台的运动，它也面临一些新的挑战。最大的问题是：平台运动使得在岸基雷达中较窄的两根一阶海杂波谱线展宽明显，当平台慢速运动时，两根一阶海杂波频率谱线展宽成两个通带；当平台快速运动时，两根一阶海杂波频率谱线融合成一个通带。船载高频地波雷达的主要探测目标为海面上的舰船，这种目标的速度相对较慢，对应的多普勒频率很可能位于展宽的一阶海杂波频谱之内，而这种船目标的检测又相对困难。

传统的船载高频地波雷达信号处理方法是二维傅里叶变换加数字波束形成(FFT-DBF)，之后进行目标检测，由于没有经过杂波抑制的过程，所以FT-DBF输出信杂噪比较低，目标检测性能较差；而其他方法主要是先经过海杂波抑制再进行目标检测，输出的信杂噪比有所提高，目标检测性能得到改善，但是却增加了算法复杂度，耗费了更多的计算时间。目前，海杂波抑制方面的研究比较有限。其中，正交投影算法OW是一种基于一阶海杂波空时分布的比较有效的海杂波抑制方法。然而，当目标和一阶海杂波位于同一波束主瓣之内时，目标增益将会损失；而船载平台的雷达天线孔径有限，方位分辨率较差，使得目标和一阶海杂波位于同一波束主瓣的情况普遍存在。斜投影算法OP的原理与OW类似，均是在空域进行海杂波抑制，即通过在特定频率海杂波的空间方位形成凹口实现的；OP的提出克服了正交投影存在目标增益损失的这一缺陷，在特定的多普勒频率位置，斜投影算法在空域一阶海杂波位置形成的凹口更窄，解决了目标增益损失的问题。但是，当一阶海杂波和目标的方位不在同一波束主瓣之内时，相比OW，OP的海杂波抑制性能较差。传统的空时自适应处理方法是二维联合处理方法，从算法原理上分析，海杂波抑制性能较OW和OP有很大改善，然而，其计算量相当大。为了解决传统空时自适应处理算法计算量大的问题，降维的空时自适应处理算法JDL被提出，但是由于船载高频地波雷达方位分辨率较差和杂波协方差矩阵估计误差较大，所以特定方位的海杂波抑制性能依然是很有限的，而且目标方位估计误差大。

技术实现要素：

本发明的目的是为了克服传统方法FFT-DBF输出信杂噪比低，目标检测性能差的问题，而提出船载高频地波雷达海杂波时域抑制方法。

船载高频地波雷达海杂波时域抑制方法具体过程为：

步骤一：利用不含船目标的实测船载高频地波雷达三维数据，针对待检测距离单元，对阵元数据进行数字波束形成(DBF)处理，得到不含有船目标的待检测距离单元和方位的时域数据x_D；

所述高频为3-30MHz；

步骤二：利用经过步骤一得到的不含有船目标的待检测距离单元和方位的时域数据对建立的RBF神经网络进行训练；得到训练好的RBF神经网络；

RBF神经网络即径向基函数神经网络；

步骤三：根据含有船目标的船载高频地波雷达三维数据，利用已经训练好的RBF神经网络对待检测距离单元和方位的时域数据进行预测、对消，得到对消结果z'[k₀]；

步骤四：对对消结果z'[k₀]进行处理，得到待检测距离单元和方位的多普勒结果z'_fft。

本发明的有益效果为：

本发明提出了一种在时域进行海杂波抑制的方法PCRBF，PCRBF使用的数据是三维的，第一维度代表天线阵元，第二维度代表时域数据，第三维度代表距离单元。PCRBF首先利用含有船目标的船载高频地波雷达三维数据，对待检测距离单元数据进行数字波束形成处理，得到待检测距离单元和方位的时域数据；然后利用已经训练好的RBF神经网络对待检测距离单元和方位的时域数据进行预测和对消；而后将对消结果进行傅里叶变换得到待检测距离单元和方位的频域数据。

本发明是一种建立在RBF预测网络之上的设计方法，相对于其它方法，该发明有以下几点优势：第一，本发明时域预测对消方法PCRBF与现有技术二维傅里叶变换加数字波束形成FFT-DBF相比，本发明方法具有很好海杂波抑制效果，输出信杂噪比提高了大约5dB，改善了目标检测性能；第二，本发明与正交投影方法OW、斜投影方法OP、和JDL相比，计算复杂度降低，需要较少的运算时间。

附图说明

图1为船载均匀线性阵列示意图；

图2为信号处理框图；

图3为RBF网络框架图；

图4为FFT-DBF和PCRBF两种算法的频域性能比较示意图，SCNR_out为输出信杂噪比。

具体实施方式

具体实施方式一：本实施方式的船载高频地波雷达海杂波时域抑制方法具体过程为：

步骤一：利用不含船目标的实测船载高频地波雷达三维数据，针对待检测距离单元，对阵元数据进行数字波束形成(DBF)处理，得到不含有船目标的待检测距离单元和方位的时域数据x_D；

所述高频为3-30MHz；

步骤二：利用经过步骤一得到的不含有船目标的待检测距离单元和方位的时域数据对建立的RBF神经网络进行训练；得到训练好的RBF神经网络；

RBF神经网络即径向基函数神经网络；

步骤三：根据含有船目标的船载高频地波雷达三维数据，利用训练好的RBF神经网络对待检测距离单元和方位的时域数据进行预测、对消，得到对消结果z′[k₀]；

步骤四：对对消结果z′[k₀]进行处理，得到待检测距离单元和方位的多普勒结果z'_fft。

具体实施方式二：本实施方式与具体实施方式一不同的是：所述步骤一中利用不含船目标的实测船载高频地波雷达三维数据，针对待检测距离单元，对阵元数据进行数字波束形成处理，得到不含有船目标的待检测距离单元和方位的时域数据x_D；具体过程为：

船载高频地波雷达一阶海杂波的空时分布

船载均匀线性阵列天线示意图如图1，V_p代表雷达平台运动速度，d代表相邻的两个天线阵元的距离，φ₀是雷达回波的入射方向和雷达平台运动方向夹角。

理想条件下，一阶海杂波空时分布如下：

f_d＝f_dpcosφ₀±f_B

其中，f_d是运动平台和一阶海杂波共同产生的多普勒频率，f_dp＝2V_p/λ是运动平台产生的多普勒频率，λ是雷达波长，是一阶海杂波布拉格频率，f_c代表载波频率(MHz)。

时域海杂波抑制原理

对于PCRBF算法，针对待检测距离单元的二维数据，雷达信号处理框图如图2。

对不含有船目标的实测船载高频地波雷达地波三维数据中第三维待检测距离单元和第一维天线阵元n∈{1,2，...，N}，定义第二维时域数据x_n＝[x_n1…x_nM]^T；

其中N＝7是天线阵元个数，下标n是天线阵元标号，上标T代表转置，M是船载高频地波雷达第二维时域数据长度(即相参积累脉冲个数)；

因此，不含有船目标的实测船载高频地波雷达待检测距离单元的二维数据表示为：

X＝[x₁,x₂,…,x_N]^T (1)

对(1)进行如下的数字波束形成处理，得到不含有船目标的船载高频地波雷达待检测距离单元和方位的时域数据；

其中，j²＝-1，j是虚数单位，φ₁是待检测方位，x_D为不含有船目标的船载高频地波雷达待检测距离单元和方位的时域数据。

其它步骤及参数与具体实施方式一相同。

实施方式三：本实施方式与具体实施方式一或二不同的是：所述步骤二中利用经过步骤一得到的不含有船目标的待检测距离单元和方位的时域数据对建立的RBF神经网络进行训练；得到训练好的RBF神经网络；具体过程为：

设定RBF神经网络输入向量为[x_D[k],x_D[k+τ],…,x_D[k+(L-1)τ]]；

所述，x_D[k]是不含有船目标的船载高频地波雷达待检测距离单元和方位的时域数据x_D的第k个元素，其中k＝1,2,…,M-(L-1)τ-1；

不含有船目标的RBF神经网络输出向量为

所述，y[k₀]是x_D[k₀]的预测值，即不含有船目标的RBF神经网络输出数据；x_D[k₀]为真值，即不含有船目标的船载高频地波雷达待检测距离单元和方位的时域数据x_D的第k₀个元素；其中k₀＝k+(L-1)τ+1，i＝1,2,…,p，w_i为相应的权值，p为中间层神经元个数，为正整数；为径向基函数(高斯函数)；

由于k₀不断变化，y是一个预测向量，目标向量为[x_D[1+(L-1)τ+1],x_D[2+(L-1)τ+1],…,x_D[M]]，当预测向量和目标向量的均方误差小于等于设定值时(例如0.01倍的x_D中元素的幅度平均值)，此刻得到的权值w_i作为最优权值，即可得到训练好的RBF神经网络。

RBF神经网络被用作预测模型，RBF网络框架如图3。

说明：x_D[k]是向量x_D的第k个元素，其中k＝1,2,…,M-(L-1)τ-1。两个连续的输入元素间隔为τ，每个输入向量的元素个数为L，径向基函数为高斯函数其中i＝1,2,…,p，相应的权值是w_i，中间层神经元个数是p，网络输出是x_D[k₀]的预测值，其中k₀＝k+(L-1)τ+1。

实际应用中，向量x_D是复向量，然而RBF神经网络的输入向量是实向量，因此预测需要将实部和虚部分别进行，之后将预测结果的实部和虚部进行复数重构。尽管神经网络是非线性的，但是只要实部和虚部预测误差足够小，将两部分的预测结果进行组合重构是合理的。

利用不含船目标的实测船载高频地波雷达三维数据(相比其他噪声，实测数据中海杂波必须占优)，经过步骤一得到待检测距离单元和方位的时域数据，应用此时域数据对步骤二中构建的RBF神经网络进行训练，具体过程为：

训练函数应用MATLAB软件工具箱中的newrb函数，用法为net＝newrb(p₀,t₀,goal,spread,m₀)，p₀代表图3中的输入向量，t₀代表输出结果，训练网络时，t₀的取值为x_D[k₀]，网络训练需要达到的均方误差水平goal取值为0.01倍的向量x_D各元素的平均幅度，中间层高斯函数展宽参数spread取值为100，中间层神经元最大数目m₀取值为7。

其它步骤及参数与具体实施方式一或二相同。

实施方式四：本实施方式与具体实施方式一至三之一不同的是：所述步骤三中根据含有船目标的船载高频地波雷达三维数据，利用已经训练好的RBF神经网络对待检测距离单元和方位的时域数据进行预测、对消，得到对消结果z′[k₀]；具体过程为：

对含有船目标的船载高频地波雷达三维数据中第三维待检测距离单元和第一维天线阵元n∈{1,2，…，N}，定义第二维时域数据x′_n＝[x′_n1…x′_nM]^T，

因此，含有船目标的船载高频地波雷达待检测距离单元的二维数据表示为：

X′＝[x₁′,x′₂,…,x′_N]^T (3)

对(1)进行如下的数字波束形成处理，得到含有船目标的船载高频地波雷达待检测距离单元和方位的时域数据；

其中，x′_D为含有船目标的船载高频地波雷达待检测距离单元和方位的时域数据，

利用已经训练好的RBF神经网络对得到的待检测距离单元和方位的时域数据进行预测和对消处理，具体过程为：

预测过程参照图3所示，含有船目标的RBF神经网络输出数据为预测结果，对消过程如下：

z′[k₀]＝x′_D[k₀]-y′[k₀] (5)

式中z′[k₀]为对消结果，y′[k₀]是含有船目标的RBF神经网络输出数据；x′_D[k₀]为含有船目标的船载高频地波雷达待检测距离单元和方位的时域数据x′_D的第k₀个元素；

其它步骤及参数与具体实施方式一至三之一相同。

具体实施方式五：本实施方式与具体实施方式一至四之一不同的是：所述步骤四中对对消结果z′[k₀]进行处理，得到待检测距离单元和方位的多普勒结果z'_fft；具体过程为：

为了保持快速傅里叶变换FFT之后的频率分辨率，向量z'中的前(L-1)τ+1个元素被赋予零值；快速傅里叶变换过程如下：

z'_fft＝FFT(z′) (6)

其中，z'_fft是待检测距离单元和方位的多普勒结果，FFT是快速傅里叶变换符号。

其它步骤及参数与具体实施方式一至四之一相同。

具体实施方式六：本实施方式与具体实施方式一至五之一不同的是：所述是船载高频地波雷达第二维时域数据长度(即相参积累脉冲个数)M取值为1024。

其它步骤及参数与具体实施方式一至五之一相同。

采用以下实施例验证本发明的有益效果：

实施例一：

本实施例船载高频地波雷达海杂波时域抑制方法具体是按照以下步骤制备的：

为了实现海杂波抑制性能和计算复杂度之间的平衡，L的值设定为5，p的值取决于设定的均方误差水平，为了实现输入向量中相邻元素之间的独立性和相关性的折中，τ的值设定为4。本专利中用于训练RBF网络的数据是不含船目标的实测船载高频地波雷达三维数据，用于目标检测的数据是通过向实测的船载高频地波雷达三维数据加入一个仿真船目标得到的，在船目标的仿真中不考虑船目标回波强度随距离的衰减，船目标的仿真方法如下：

S(r)＝am[r]×a_s×S_st (7)

S_st＝s_s(φ_0t)×s_t(f_0t)^T (8)

a_s＝a_cn×10^SCNR/20 (11)

am＝FT(ph⊙h) (13)

其中，S是仿真船目标最终需要的三维数据块，am[r]是向量am中的第r个元素，a_s是所加船目标的平均幅度，S_st是所加船目标的空时导向矢量，f_0t和φ_0t分别是所加船目标的频率和方位。另外，T_r是脉冲重复间隔，a_cn是海杂波加噪声的平均幅度，SCNR代表输入的信杂噪比，am是在距离维进行归一化的复幅度，FT代表傅里叶变换，h是汉宁窗，ph是所加船目标在距离维上的相位信息，⊙代表哈达玛乘积，r₀是所加船目标的距离单元，r_max是距离单元个数。

本专利中所用数据的参数如下：d是14m/s，T_r是0.26s，λ是56.786m，M是1024，N是7，r_max是256。实测的船载高频地波雷达海杂波数据来自1998年进行的船载雷达平台试验，雷达平台运动速度是5m/s。本专利使用的三维数据是将频率和方位分别为-0.2044Hz和130°的仿真船目标加到实测的船载高频地波雷达三维数据的第17个距离单元而得到。

为了比较两种算法的性能，我们利用输出信杂噪比的概念，计算方法如式(15)：

其中，K是频域数据长度，A₀代表所关心位置的幅度，A_n代表所关心位置的周围幅度。

实验中，加入仿真船目标的信杂噪比SCNR为-15dB，FFT-DBF和PCRBF算法的多普勒结果如图4，纵坐标是归一化幅度取dB。

利用式(15)计算算法的输出信杂噪比SCNR_out，由图4可以看出：FFT-DBF算法目标位置为22.50dB，PCRBF算法目标位置的输出信杂噪比为27.05dB。如果门限设定为25dB，PCRBF算法能检测到船目标，而FFT-DBF算法不能；如果门限设定为13dB，两种算法均能检测到船目标，但是此时FFT-DBF算法存在虚警，所以PCRBF算法能够起到一定的海杂波抑制效果，性能较好。

本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。