一种基于局部相关性的改进子空间海杂波抑制方法
【专利摘要】本发明公开一种基于局部相关性的改进子空间海杂波抑制方法。海杂波的抑制能够提高目标信杂比，有助于提高高频地波雷达(high?frequency?surface?wave?radar,HFSWR)对舰船目标的检测性能。原始子空间法抑制海杂波时对所有距离单元构造的协方差矩阵是相同的，造成对杂波空间估计不准从而抑制效果不佳，距离多普勒谱上反映为不同距离海杂波抑制效果类似。本发明通过计算待检测距离单元与相邻距离单元之间的相关系数，自适应的确定组成协方差矩阵的参考单元，对构成的协方差矩阵进行奇异值分解(singular?value?decomposition,SVD)分解得到海杂波所在的子空间，并最终将其抑制。该方法具有同时抑制一阶和高阶海杂波的能力，提高了目标信杂比，对所有距离单元的抑制效果都较理想。
【专利说明】一种基于局部相关性的改进子空间海杂波抑制方法
【技术领域】：
[0001] 本发明涉及一种高频地波雷达海杂波抑制方法，特别涉及一种基于局部相关性的 改进子空间海杂波抑制方法。
【背景技术】：
[0002] 高频地波雷达工作在3?30MHz频段，是超视距雷达的一种，其利用垂直极化的电 磁波沿海表面绕射传播损耗小的特点，可实现对海洋动力学参数、海上慢速目标及低空飞 机的监视与定位，近年来广泛应用于军事与民用领域。
[0003] 高频地波雷达回波情况复杂，除了包括射频干扰、瞬态干扰和电离层干扰等各类 干扰外，还包括很强的海杂波。海杂波能量主要集中在两个由雷达工作频率决定的一阶海 杂波，即Bragg峰处，实际情况时，一阶Bragg峰会有不同程度的展宽。对于舰船目标检测 来说，由于舰船目标的多普勒频率和强大的海杂波相近，海杂波才是限制检测性能的主要 因素。因此相对于飞机目标，舰船目标检测也更为困难。所以，在进行舰船目标检测时，一 个重要的思想就是要先抑制海杂波，提高信杂比，以期对后面的检测有所帮助。
[0004] 现有较常用的海杂波抑制方法主要有以下两个思路：一类是时域对消类方法， 如循环对消方法和基于幅度相位估计（amplitude and phase estimation, APES)的对 消方法。这类方法通过估计强杂波对应谐波分量的频率、幅度和相位3个参数，有选择 性的对消掉强杂波分量，使目标信号从杂波中显露出来，这类方法的性能取决于上述3 个参数的估计精度。（如文献 l:R〇〇t B.HF radar ship detection through clutter cancellation[C]//Radar Conference, 1998. RADARC0N98. Proceedings of thel998IEEE. IEEE，1998:281-286.)(如文献2 :郭欣，倪晋麟，刘国岁.短相干积累条件下天波超视距雷 达的舰船检测[J].电子与信息学报，2004, 26(4) :613-618)。另一类是基于矩阵分解的海 杂波方法，如基于子空间的方法、基于奇异值分解（singular-value decomposition, SVD) 的方法和通过多普勒频率估计杂波子空间（estimate subspace via Doppler, ESVID)的方 法，这类方法根据杂波在子空间的聚集特性来实现海杂波抑制。其中，基于子空间方法通过 参考单元数据来估计杂波子空间，（如文献3 :邢孟道，保铮，强勇.天波超视距雷达瞬态干 扰抑制[J].电子学报，2002, 30(6) :823-826.) ;SVD方法通过待检测单元来估计杂波子空 间，（如文献4 :Poon M W Y, Khan R H，Le_Ngoc S. A singular value decomposition (SVD) based method for suppressing ocean clutter in high frequency radar[J]. Signal Processing, IEEE Transactions on, 1993, 41 (3) :1421-1425.) ;ESVID 方法根据杂波子 空间与信号矢量张成子空间的等价性，通过多普勒频率估计来估计杂波子空间，（如文 献 5 :Zhiguo Z,Chen J, Zheng B.A method to estimate subspace via Doppler for ocean clutter suppression in skywave radars[C]//Radar(Radar), 2011IEEE CIE International Conference on.IEEE,2011, 1:145-148.)〇
[0005] 除了上述较常用的两大类方法外，一些新颖的方法也被提出。如中国专利公开 CN102928823B基于二次雷达应答信号相关特征的海杂波抑制方法中提出将受到海杂波影 响的探测目标计算其真实应答信号和受海杂波效应的反射信号之间的相对时间差，根据不 同相对时间差的码位值演变方式，生成与真实目标相对应的目标航管代码并保存；将真实 目标航管代码与历史探测点迹数据进行代码相关，相关成功后对探测目标航管代码进行纠 码并实现海杂波的抑制。如中国专利公开CN103645467A海杂波抑制以及海杂波背景中目 标检测的方法和系统中提出对海面搜索雷达实测的一维距离像的历程图进行二维傅里叶 变换后得到估计的海杂波一维距离像的历程图；将实测的一维距离像的历程图的图像数据 减去估计的海杂波一维距离像的历程图的图像数据，得到抑制海杂波的一维距离像的历程 图的图像数据。如中国专利公开CN102779336A基于混沌的光电图像海杂波抑制方法中提 出利用光电图像海杂波所具有的混沌特性，采用RBF神经网络方法预测与抑制海杂波。如 中国专利公开CN102169175B基于FRFT域非相参积累的海杂波抑制和动目标检测方法中通 过对FRFT域延迟数据的迭代平滑，分别采用FRFT域均值积累方法和FRFT域单极点反馈积 累方法，完成非相参积累，将FRFT域信号与积累后的数据对消，抑制海杂波。
[0006] 时域对消类海杂波抑制方法的效果很难保证，当脉冲相干积累时间较短时则频率 分辨力低，此时对Bragg峰的频率和幅度进行估计时精度下降，无法得到较好的初相估计， 造成循环对消次数增加，而随着对消的不断进行，剩余杂波将会在多普勒中扩散开来，影响 目标的检测。
[0007] 子空间的海杂波抑制方法中，对于一个进行完距离向脉冲压缩和方位向波束形成 后的三维数据块，取整个数据块进行协方差估计，即认为整个数据块的海杂波是高度相关 的，用的是同一个标准进行海杂波抑制，这显然是不合理的。海杂波确实在一定的范围内高 度相关，但这个范围并没有大到整个雷达探测的范围。对于整体用一个标准来进行海杂波 抑制，有可能出现近处海杂波抑制不干净或远处海杂波抑制过多的情况。
[0008] 那么，对于不同距离不同方位的待抑制回波数据，应该设计不同的参考单元个数 来构造与其对应的协方差矩阵，从而进行海杂波的抑制。参考单元的个数，以及协方差矩阵 的构造直接影响到海杂波抑制效果。参考单元个数取的过大，海杂波空间相关性减弱，协方 差估计不准，在很大一片区域都用同一个标准进行抑制不符合海杂波在空间分布的特点； 参考单元个数取的过小，则样本数太小也无法进行准确估计。
【发明内容】
：
[0009] 本发明为了提高高频地波雷达舰船目标检测性能，提出了基于局部相关性的改进 子空间海杂波抑制方法，对于每个特定的方位和距离单元的待抑制数据，通过计算周围单 元与待抑制单元的相关系数，自适应的调整块的大小以期获得更精确的协方差矩阵估计， 继而进行后续的海杂波抑制。
[0010] 本发明采用如下技术方案：一种基于局部相关性的改进子空间海杂波抑制方法， 其包括如下步骤：
[0011] (a)、设输入数据为距离向脉冲压缩，方位向数字波束形成后的距离-慢时间-方 位三维数据块；
[0012] (b)、每次抑制的对象为一固定方位固定距离门的一维慢时间数据，
[0013] 设当前待检测的某确定方位和距离的距离单元的K个回波慢时间数据为矢 量s。= [Sc;(o)，s。⑴，...，Sc;(K-i)]T，第i个距离单元的K个回波慢时间数据为Si =
[Si(o)，Si(i)，...，Si(K-i)]T，其中 T为转置；
[0014] (c)、遍历出与当前待检测距离单元相邻的且方位与待检测距离单元相同的距离 单元慢时间数据Sy并计算其与待检测距离单元S。的相关系数
[0015]
【权利要求】
1. 一种基于局部相关性的改进子空间海杂波抑制方法，其特征在于：包括如下步骤 (a) 、设输入数据为距离向脉冲压缩，方位向数字波束形成后的距离-慢时间-方位三 维数据块； (b) 、每次抑制的对象为一固定方位固定距离门的一维慢时间数据， 设当前待检测的某确定方位和距离的距离单元的K个回波慢时间数据为矢量 S。= [Sc;(0)，Sc;(l)，...，Sc;(K-l)]T，第i个距离单元的K个回波慢时间数据为Si = [Si(o)，Si(i)，...，Si(K-i)] T，其中 T为转置； (C)、遍历出与当前待检测距离单元相邻的且方位与待检测距离单元相同的距离单元 慢时间数据Si，并计算其与待检测距离单元S。的相关系数
其中η代表共轭转置，11112代表二范数，其中〇< | Pil < 1，1表示相关性最大，〇表 示相关性最小，设计门限a找到P彡a的P个对应S。的参考单元数据Si(i = 1，2，...P)， 和对应的相关系数P i (i = 1，2,. . . P); (d) 、按距离从近到远将所选P个参考单元数据排列，得到X = [S,, S9,. . . S丄按照对应 的顺序，构造 w = diag( P P 2,. . · P p),则S。的估计协方差矩阵戈 (e) 、将R进行SVD奇异值分解
R = V Σ VH Σ = diag ( λ丨，λ 2, · · · , λ κ)，得到K个奇异值（λ丨彡λ 2彡…彡λ κ)，V = [Vl，ν2,. . . νκ]，Vi为对应λ i的特征向量，且V为酉矩阵，所述SVD分解为矩阵的一种正交 分解，矩阵的主要分量与大奇异值相对应，定义归一化奇异值为： Η 贝U，〇<瓦引，选择一个小正数作为阈值〇，则石比0大的前r个特征值对应的特征向 量所张成的空间为信号子空间即海杂波所在空间V。= [Vl，v2, ... vj，其余特征向量所张成 的空间为噪声子空间； ⑴、将待抑制的S。投影到信号子空间V。= [Vl，v2, . . . vj上，得到滤波后的输出为
其中I为KXK维单位矩阵； (g)、重复（b)?（f)的步骤，检测当前距离-慢时间-方位三维数据块中的待检测方 位上的所有距离单元数据，接着选择其它方位，直至得出最后一个方位上的最后一个距离 单元的数据，则当前距离-慢时间-方位三维数据块海杂波抑制完成。
【文档编号】G01S7/36GK104155632SQ201410346178
【公开日】2014年11月19日 申请日期:2014年7月18日 优先权日:2014年7月18日
【发明者】赵慧俊, 周涛, 张弓, 徐坤, 胡文, 王鑫海, 万明雪, 张璐, 季秀霞 申请人:南京航空航天大学