基于最小二乘逼近的一阶海杂波检测方法
【专利摘要】基于最小二乘逼近的一阶海杂波检测方法，本发明涉及一阶海杂波检测方法。本发明是为了解决现有技术查找到的布拉格峰位置不准确及虚警率高的问题。本发明方法首先由雷达的工作频率f0，求得理论布拉格频率fB，并在每个海元上划定噪声区，求得每个海元的阈值；根据雷达探测海域海态数据，确定洋流径向流速的最大值Vm，并求得布拉格频偏的最大值Δfm；根据Δfm和fB确定搜索范围，进而确定布拉格峰的疑似点，并对疑似点的幅度值进行阈值确认，判断当前海元是否有效；其次对于有效海元，取适当个数疑似点的临近点构成点集，利用最小二乘逼近法求得该点集的逼近函数；最后取逼近函数的峰值点作为布拉格峰；本发明应用于高频地波雷达探测领域。
【专利说明】
基于最小二乘逼近的一阶海杂波检测方法
技术领域
[0001] 本发明涉及一阶海杂波检测方法，涉及高频地波雷达探测领域。
【背景技术】
[0002] 高频地波雷达，可以突破地球曲率的限制，探测到视线以下的目标。在雷达回波中 往往掺杂着大量的干扰和噪声。高频雷达一阶海杂波是雷达回波中最主要的干扰之一，它 的能量往往很高，严重影响目标的检测，应该尽量抑制海杂波;但从另一方面来讲，海杂波 包含着大量的海态信息，是海态反演的主要检测对象。海杂波的产生机理可以由布拉格谐 振原理来解释，即当高频电磁波照射到粗糙的海面上时，电磁波会与海面相互作用发生强 烈的散射作用，从而产生一阶海杂波。理论上，一阶海杂波在雷达回波RD谱上会在频率轴上 形成一对关于零频对称的能量尖峰，即为Bragg峰。Bragg峰对于其存在频率上的海面目标 来说是非常大的干扰，将直接导致该频率范围内目标的淹没，在有洋流，电离层杂波和强风 等情况下，Bragg峰还会分裂展宽，扩大目标检测的干扰范围。而通过测量Bragg峰位置的偏 移，左右Bragg峰能量的比值，一阶Bragg峰与二阶Bragg峰能量比值又能得到海面浪场、风 场、流场等十分丰富的海态信息。所以，无论是海态信息的提取还是对海面目标的跟踪监 测，海杂波的检测都是研究的重点。
[0003] 经典的一阶海杂波检测方法有局部峰值检测法和基于特征识别法等，经典检测方 法的共同特点是将检测范围内的最大值点视为Bragg (布拉格)峰，并未对Bragg峰附近的谱 做任何处理。但相参积累时间过高会带来频谱分裂，这会使Bragg峰分裂为两个峰;高频噪 声干扰会使Bragg峰附近的谱有毛刺，这两种情况下，都会使得查找到的Bragg峰位置不准 确，虚警率非常高。
【发明内容】
[0004] 本发明是为了解决现有技术查找到的Bragg峰位置不准确及虚警率高的问题，而 提出的基于最小二乘逼近的一阶海杂波检测方法。
[0005] 基于最小二乘逼近的一阶海杂波检测方法按以下步骤实现：
[0006] 步骤一：由雷达的工作频率fQ，求得理论布拉格频率fB;
[0007] 步骤二:根据步骤一得到的fB在每个海元上划定噪声区，并求得每个海元的阈值；
[0008] 步骤三:根据雷达探测海域海态数据，确定洋流径向流速的最大值Vm，并求得布拉 格频偏的最大值A fm;
[0009] 步骤四：根据步骤三中求得的布拉格频偏的最大值A fm与步骤一中得到的理论布 拉格频率fB，确定布拉格峰的搜索范围；
[0010] 步骤五:在步骤四确定的搜索范围内取最大值点作为布拉格峰疑似点；
[0011] 步骤六:将步骤五确定的疑似点与步骤二得到的阈值作比较，若疑似点的幅度值 大于阈值，则认为海元数据有效，为有效海元，执行步骤七;否则，检测下一海元；
[0012] 步骤七:在步骤六中确定的有效海元内，取以布拉格峰疑似点为中心的点做点集， 求取逼近函数y;
[0013]步骤八:取步骤七中得到的逼近函数y的峰值点记为布拉格峰值点。
[0014]发明效果：
[0015]本发明相对传统检测算法，求取了 Bragg峰附近谱的逼近函数，相当于对Bragg峰 附近的谱进行了平滑处理，减弱了频谱分裂和高频噪声带来的干扰，可有效提高一阶海杂 波检测的准确率。
[0016]本发明利用最小二乘逼近法求得高频雷达回波多普勒谱的逼近函数，对Bragg峰 附近的谱进行了平滑处理，降低了由相参积累时间过高带来的频谱分裂和高频噪声等干扰 的影响，解决了一阶海杂波检测中Bragg峰位置识别不准确的问题。
【附图说明】
[0017]图1为距离多普勒谱图；
[0018]图2为多普勒谱图；
[0019]图3为利用最小二乘逼近法处理后的RD谱结果图；
[0020]图4为第4个波束在10km处结果图；
[0021]图5为第4个波束在50km处结果图；
[0022]图6为第4个波束在100km处结果图；
[0023]图7为第4个波束在150km处结果图；
[0024]图8为本发明流程图。
【具体实施方式】
[0025]【具体实施方式】一:如图8所示，基于最小二乘逼近的一阶海杂波检测方法包括以下 步骤：
[0026]原始雷达回波数据经两次FFT运算，一次DBF运算得到回波数据的方位-距离-多普 勒谱(ARD谱）。一阶海杂波的识别即是要在ARD谱上识别出一阶Bragg峰的位置。
[0027]步骤一：由雷达的工作频率fo,求得理论布拉格频率ft;
[0028] 步骤二:根据步骤一得到的fB在每个海元上划定噪声区，并求得每个海元的阈值；
[0029] 步骤三:根据雷达探测海域海态数据，确定洋流径向流速的最大值Vm，并求得布拉 格频偏的最大值A fm;
[0030] 步骤四：根据步骤三中求得的布拉格频偏的最大值A fm与步骤一中得到的理论布 拉格频率fB，确定布拉格峰的搜索范围；
[0031 ]步骤五:在步骤四确定的搜索范围内取最大值点作为布拉格峰疑似点；
[0032]步骤六:将步骤五确定的疑似点与步骤二得到的阈值作比较，若疑似点的幅度值 大于阈值，则认为海元数据有效，为有效海元，执行步骤七;否则，检测下一海元；
[0033]步骤七:在步骤六中确定的有效海元内，取以布拉格峰疑似点为中心的点做点集， 求取逼近函数y;
[0034] 步骤八:取步骤七中得到的逼近函数y的峰值点记为布拉格峰值点。
[0035] 算法是基于最小二乘逼近法的，在此对最小二乘逼近法进行介绍。
[0036] 在实际工程应用中，通常得到的并不是连续的函数f(x)，而是由实验或观测得到 的离散的点集(Xi，fi)，i = 0,1，...，m。当需要求点集中未给出的变量的函数值时，就必须做 出数值逼近。设/WeC[a,6],%(x)屬(x),...,仍,(x)为c[a，b]上n + 1个线性无关函数，用 {供。(x)肩(x),.名(x)}表示由灼(x)肩(x),…仍,(x)张成的线性子空间，则对任意的 p p 中；. n
[0037] 二 JXw, (A.) (1) j=Q
[0038] 最小二乘逼近问题就是要求一个我功使 tn
[0039] ||/-列|; = Z A. [乂' -沪⑷]2 =瞻||/-炉| (2 ) s：-i) 1
[0040] 其中 Pi(i = 0，l，…，m)为权值。
[0041] 求贫X)的问题等价于求多元函数 - :~\2 m n
[0042] F(a0,al,---,an) = Yipi ^{xt)-j] (3) /-0 戶0
[0043]的极小值，其法方程为： n 一
[0044] 2](多?;.=(./'，色以=() U.-" C4)
[0045] 其中， 咏，死）=Z A % (A M (A )，
[0046] H (5) - m , 仁o
[0047] 求解法方程即可得到函数故X),并且函数贫X)存在且唯一。
【具体实施方式】 [0048] 二:本实施方式与一不同的是:所述步骤一中求得理 论布拉格频率fB的公式为：
[0049] i#=±0.102^/T (6X
[0050]【具体实施方式】三:本实施方式与【具体实施方式】一或二不同的是:所述步骤二中求 得每个海元的阈值具体为：
[0051 ]每个海元的阈值为^与5冊的和，所述NP为基底噪声，SNR(根据实际情况给定)为设 定的最小信噪比。
[0052]本次算法将多普勒频率在[_2fB，2fB]范围外的区域划分为噪声区（噪声区的范围 可根据实际情况调整）。
[0053]【具体实施方式】四：本实施方式与【具体实施方式】一至三之一不同的是：所述步骤三 中求得布拉格频偏的最大值△ fm的公式为：
[0054]根据经验确定洋流径向流速的最大值Vm(本算法取lm/s，可根据实际情况调整）， 依此求出频偏的最大值，公式如下： (7)
[0056]其中所述A为雷达发射电磁波的波长。
[0057]【具体实施方式】五:本实施方式与【具体实施方式】一至四之一不同的是：所述步骤四 中确定布拉格峰的搜索范围的具体为：[_fB_ A fm~-fB+ A fm]~[fB- A fm~fB+ A fm]。
[0058]【具体实施方式】六:本实施方式与【具体实施方式】一至四之一不同的是:所述步骤七 中取以布拉格峰疑似点为中心的点做点集，求取逼近函数y的具体过程为：
[0059] 取以布拉格峰疑似点为中心的9~15个点做点集（点集个数根据相参积累时间调 整)用来求取逼近函数，并去除奇异点，其中所述奇异点为幅度值小于左右两侧点的幅度值 的点；
[0060] -阶Bragg峰的形状类似抛物线，因此取①= Span{l，z，'"，zn}，Pi = l(i = l，2,…， m);其中①为线性子空间，1，z，"_，zn为n+1个选定的线性无关函数，Pi为权函数，m代表点的 个数；
[0061]确定逼近函数的形式为公式(8)
[0062] y = ax2+bx+c (8)
[0063] 其中x为每个点对应的多普勒频率，y为每个点对应的幅度值；
[0064]采用公式(9)求解系数a、b、c: ni
[0065] \ ，=l (9) m (v,%)=Y.y^,)
[0066] 其中％ =_r'(y = 0,1,2)，济=Z(& = 0,1,2)和與为选定的线性无关函数，Xi为i 点对应的多普勒频率，点对应的幅度值；由公式(9)求解出逼近函数的系数a、b、c后， 即得到了逼近函数。
[0067] 实施例一：
[0068]实验中采用的数据来源于哈尔滨工业大学电子所某雷达实验站于2015年8月采集 的数据，表1给出雷达的相关参数。
[0069] 表1雷达系统参数
[0072]根据相关要求，部分参数不予给出，以*代替。
[0073]原始雷达回波数据经两次FFT运算、一次DBF运算后得到回波数据的方位-距离-多 普勒谱(ARD谱），由ARD谱可得到每个波束的距离多普勒谱(RD谱)和每个海元的多普勒谱， 如图1和图2所示。图1中两条相对零频对称的亮带是一阶海杂波；图2中相对零频对称的两 个尖峰是一阶海杂波。
[0074]图3给出利用最小二乘逼近法处理之后的RD图，并用紫色圆圈标出了查找到的 Bragg峰。从图中可看出，紫色圆圈基本覆盖了一阶海杂波的两条亮带，说明查找结果与原 RD谱中的一阶海杂波吻合程度很高，也即利用该方法准确的找到了左右Bragg峰的位置。 [0075]图4至图7给出了利用最小二乘逼近法处理之后的多普勒回波谱，图中用竖线标注 理论Bragg峰的位置，用*点标注算法识别出的Bragg峰的位置。从图中可看出，在不同海元 内，本算法都准确的找到了Bragg峰的位置，验证了算法的适用性和正确性。且图中右上方 给出了 Bragg峰的左侧频偏和右侧频偏及频偏差值。理论上，Bragg峰的左侧频偏和右侧频 偏应该相等。因此，识别的Bragg峰的左右频偏差越小，说明算法效果越好。表2给出采用最 小二乘逼近法时不同距离范围内的左右频偏差的均值和均方根误差表。
[0076]表2左右频偏差均值和均方根误差表
[0078]从表2可看出，由最小二乘逼近法得到的左右频偏差的均值和均方根误差都较小， 其中，左右频偏差的均值控制在0.0057Hz以内，左右频偏差的均方根误差控制在0.0112Hz 以内，相对于传统算法局部峰值法的识别结果，左右频偏差的均值提高了9.5%，左右频偏 差的均方根误差提高了 1.75%。充分说明最小二乘逼近法的效果较好，大大提高了一阶海 杂波识别结果的准确率，能够准确识别一阶海杂波。
【主权项】
1. 基于最小二乘逼近的一阶海杂波检测方法，其特征在于，所述基于最小二乘逼近的 一阶海杂波检测方法包括以下步骤： 步骤一：由雷达的工作频率f O，求得理论布拉格频率fB ; 步骤二:根据步骤一得到的f B在每个海元上划定噪声区，并求得每个海元的阈值； 步骤三:根据雷达探测海域海态数据，确定洋流径向流速的最大值￥"，并求得布拉格频 偏的最大值Afm; 步骤四：根据步骤三中求得的布拉格频偏的最大值A ^与步骤一中得到的理论布拉格 频率fB，确定布拉格峰的搜索范围； 步骤五:在步骤四确定的搜索范围内取最大值点作为布拉格峰疑似点； 步骤六:将步骤五确定的疑似点与步骤二得到的阈值作比较，若疑似点的幅度值大于 阈值，则认为海元数据有效，为有效海元，执行步骤七;否则，检测下一海元； 步骤七:在步骤六中确定的有效海元内，取以布拉格峰疑似点为中心的点做点集，求取 逼近函数y; 步骤八:取步骤七中得到的逼近函数y的峰值点记为布拉格峰值点。2. 根据权利要求1所述的基于最小二乘逼近的一阶海杂波检测方法，其特征在于所述 步骤一中求得理论布拉格频率f B的公式为：(6) ο3. 根据权利要求2所述的基于最小二乘逼近的一阶海杂波检测方法，其特征在于所述 步骤二中求得每个海元的阈值具体为： 每个海元的阈值为^与5顺的和，所述Np为基底噪声，SNR为设定的最小信噪比。4. 根据权利要求3所述的基于最小二乘逼近的一阶海杂波检测方法，其特征在于所述 步骤三中求得布拉格频偏的最大值A匕的公式为：(7): 其中所述λ为雷达发射电磁波的波长。5. 根据权利要求4所述的基于最小二乘逼近的一阶海杂波检测方法，其特征在于所述 步骤四中确定布拉格峰的搜索范围的具体为：[_fB_ Δ fm~-fB+ Δ fm]~[fB- Δ fm~fB+ Δ fm]。6. 根据权利要求5所述的基于最小二乘逼近的一阶海杂波检测方法，其特征在于所述 步骤七中取以布拉格峰疑似点为中心的点做点集，求取逼近函数y的具体过程为： 取以布拉格峰疑似点为中心的9~15个点做点集，并去除奇异点，其中所述奇异点为幅 度值小于左右两侧点的幅度值的点； 取Φ = SpanU，z，…，zn}，Pi=l(i = l，2,…，m);其中 Φ为线性子空间，I，z，…，zn为n+1 个选定的线性无关函数，P1为权函数，m代表点的个数； 确定逼近函数的形式为公式(8) y = ax2+bx+c (8) 其中X为每个点对应的多普勒频率； 采用公式(9)求解系数a、b、c: (9) 其中％ = .r;( / = 0,1,2)，仍.=V 0 = 0,1,2)，以为i点对应的多普勒频率，yi为i点对应的 幅度值；由公式(9)求解出逼近函数的系数a、b、c后，即得到了逼近函数。
【文档编号】G01S7/292GK105891787SQ201610206443
【公开日】2016年8月24日
【申请日】2016年4月5日
【发明人】董英凝, 赵东阳, 张贺磊, 邓维波
【申请人】哈尔滨工业大学