对海雷达小目标检测方法.pdf

《对海雷达小目标检测方法.pdf》由会员分享，可在线阅读，更多相关《对海雷达小目标检测方法.pdf（12页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201911404578.2 (22)申请日 2019.12.31 (71)申请人 扬州船用电子仪器研究所 （中国船 舶重工集团公司第七二三研究所） 地址 225001 江苏省扬州市吴州东路186号 (72)发明人 宋佳柯涛张恒王晓华 翟玉健 (74)专利代理机构 南京理工大学专利中心 32203 代理人 岑丹 (51)Int.Cl. G01S 7/41(2006.01) (54)发明名称 对海雷达小目标检测方法 (57)摘要 本发明公开了一种对海雷达小目标检测方 法， 具体包括。

2、： 对原始回波数据进行距离维FFT处 理和MTD处理； 对MTD处理后的数据进行舰速补 偿， 提取舰速补偿后的零多普勒通道数据作为固 定杂波杂波背景； 进行固定杂波对消以及地杂 波、 海杂波抑制， 并对抑制后的数据进行MTD处理 和恒虚警处理， 完成对高速目标的检测； 根据舰 速补偿后的数据多普勒频率高低分别采用检测 前跟踪处理和扫描间相关处理完成杂波下低速 小目标检测， 所述低速目标为速度低设定阈值的 目标。 本发明提高了运行效率， 从而显著提高了 对海面小目标检测性能。 权利要求书3页 说明书5页 附图3页 CN 111610501 A 2020.09.01 CN 111610501 A 。

3、1.一种对海雷达小目标检测方法， 其特征在于， 具体包括： 对原始回波数据进行距离维FFT处理和MTD处理； 对MTD处理后的数据进行舰速补偿， 提取舰速补偿后的零多普勒通道数据作为固定杂 波杂波背景； 进行固定杂波对消以及地杂波、 海杂波抑制， 并对抑制后的数据进行MTD处理和恒虚警 处理， 完成对高速目标的检测， 所述高速目标为速度高于设定阈值的目标； 根据舰速补偿后的数据多普勒频率高低分别采用检测前跟踪处理和扫描间相关处理 完成杂波下低速小目标检测， 所述低速目标为速度低设定阈值的目标。 2.根据权利要求1所述的对海雷达小目标检测方法， 其特征在于， 对地杂波、 海杂波抑 制后的数据进行。

4、MTD处理和恒虚警处理， 完成对高速目标的检测的具体方法为： 采用二阶数字对消器对固定杂波进行固定杂波抑制； 判断固定杂波抑制后的数据中是否存在海杂波， 如果存在， 则采用自适应二次对消器 进行海杂波抑制； 对抑制后的数据进行MTD处理、 恒虚警处理， 完成对高速目标进行检测。 3.根据权利要求2所述的对海雷达小目标检测方法， 其特征在于， 所述二阶数字对消器 的差分方程为： 其中，v为平台运动速度， 为天线扫描角度， fr为多普勒频率， 为雷达 信号波长， y(n)为二阶数字对消器输出， x(n)为固定杂波输入。 4.根据权利要求2所述的对海雷达小目标检测方法， 其特征在于， 判断固定杂波抑。

5、制后 的数据中是否存在海杂波的具体方法为： (1)对固定杂波抑制后的数据进行求模； (2)计算求模结果的统计平均， 按恒虚警检测信噪比要求设置初始门限TH； 其中Pf为虚警率， U0表示检测门限； (3)将求模结果与TH进行比较， 统计小于TH的模值， 计算小于TH模值的均值， 然后乘以 一个系数作为杂波定位的低门限TH1； 统计大于TH的模值， 计算大于TH模值的均值， 将大于 TH模值的均值作为杂波定位的高门限TH2； (4)设置杂波门限个数TH3， 设N为固定杂波抑制后的信号模值连续落在区间TH1,TH2 的点数， 若NTH3则判决该区域存在海杂波， 若NTH3则判决该区域不存在海杂波。。

6、 5.根据权利要求2所述的对海雷达小目标检测方法， 其特征在于， 采用自适应二次对消 器进行海杂波抑制的具体步骤为： (1)对判决为海杂波的杂波区数据进行分段处理， 每段的长度为N； (2)对海杂波区中不同的杂波段分别进行杂波的多普勒频率估计， 采用自适应二次对 消滤波器组进行海杂波抑制， 所述自适应二次对消滤波器组的系统函数为： Hk(z)1-2exp(j2 kf/fr)z-1+exp(j4 kf/fr)z-2 k0,1,2,.N-1 式中， Hk(z)为系统函数输出， 多普勒频率在0fr之间均匀分布， 按照杂波数据分段的 权利要求书 1/3 页 2 CN 111610501 A 2 每段长。

7、度N取间隔ffr/N。 6.根据权利要求15任一所述的对海雷达小目标检测方法， 其特征在于， 对抑制后的 数据进行恒虚警处理的具体步骤为： 确定参考单元， 取参考单元的核函数的平均作为检测单元的杂波概率密度函数； 根据杂波概率密函数得到分布函数， 并根据分布函数确定门限值； 若信号检测结果高于检测门限， 则检测到目标。 7.根据权利要求6所述的对海雷达小目标检测方法， 其特征在于， 采用二维的方式来选 取参考单元数， 将被检测单元D四周单元作为参考单元。 8.根据权利要求6所述的对海雷达小目标检测方法， 其特征在于， 杂波概率密度函数 为： 其中， K为核函数， an定义为： 先选定一个大于1。

8、小于n的自然数kkn， 对任何xX， 以an (x)an(x； X1,Xn)记最小的正数a， 使x-a,x+a中至少包含样本X1,Xn中的k个， X为参 考单元， x为中心单元； 对概率密度函数积分得到的分布函数为： 门限值为： 式中， Pfa为虚警概率，为的逆函数。 9.根据权利要求8所述的对海雷达小目标检测方法， 其特征在于， 对多普勒频率高于设 定阈值的数据采用检测前跟踪处理检测小目标并获得目标航迹的具体步骤为： (1)根据目标运动模型， 在需要建立起Kalman滤波运动方程的波束数内， 设置第一门 限； (2)若回波信号能量超过第一门限， 则判别回波信号能否与已有航迹关联， 若关联成 。

9、功， 则转(3)， 否则将没有关联上的回波信号， 与之前超过第一门限但未关联上的回波信号 建立Kalman滤波运动方程， 将其作为航迹建立前的航迹数据， 等待下一个波位扫描回波信 号进行处理； (3)对关联上的回波信号与航迹， 利用Kalman滤波预测下一个波束的目标回波信号， 将 下一个波束驻留时间内， 过第一门限的回波信号， 与之前过第一门限的回波信号进行相参 积累/非相参积累； (4)将不同波位驻留时间内的积累的信号能量与对应波束驻留数的检测门限进行比 较， 若超过门限， 则判决目标存在并给出对应目标航迹， 否则将其作为航迹建立前的航迹数 据， 等待下一个波位扫描回波信号进行处理。 10。

10、.根据权利要求8所述的对海雷达小目标检测方法， 其特征在于， 对多普勒频率低于 设定阈值的数据采用扫描间相关处理完成杂波下低速小目标检测， 具体步骤为： 权利要求书 2/3 页 3 CN 111610501 A 3 (1)存储经过杂波抑制的多帧目标信息； (2)利用前几帧保存的数据与当前的目标信息在对应的距离单元组成一个窗； (3)统计窗内对应距离单元目标出现的次数， 若目标出现的次数超过检测门限， 则判决 该位置是目标。 权利要求书 3/3 页 4 CN 111610501 A 4 对海雷达小目标检测方法 技术领域 0001 本发明属于雷达信号处理技术， 具体为一种对海雷达小目标检测方法。 。

11、背景技术 0002 导航雷达是舰艇的重要装备之一， 对安全航行发挥重大作用。 在常规情况下， 当其 他监视雷达如搜索雷达、 警戒雷达处于静默状态时， 导航雷达不仅负责导航任务， 还担负监 视近距离范围内的小目标， 尤其是具有攻击性的目标， 这些目标包括小冰山、 小船、 蛙人、 望 远镜等。 海杂波背景条件下的慢速小目标一直是阻碍导航雷达性能提高的难题， 其主要原 因有两个方面， 一是导航雷达自身条件。 大部分导航雷达受安装场地、 成本等方面的限制， 天线口径小， 分辨率低， 显示和信号处理功能简单， 在强海杂波环境下很难发现小目标。 二 是目标的性质。 首先雷达分辨单元内， 小目标回波强度与海。

12、杂波回波强度相比并不占优势， 对于低分辨雷达， 往往小目标回波信号强度比海杂波还弱， 因此直接进行小目标检测非常 困难； 其次由于慢速目标相对于雷达的运动速度较低， 与海杂波相对于雷达的速度差异也 很小， 使得在多普勒域上也很难找到与海杂波的差异。 发明内容 0003 本发明的目的在于提出了一种对海雷达小目标检测方法。 0004 实现本发明目的的技术解决方案为： 0005 本发明与现有技术相比， 其显著优点为： 本发明根据目标运动谱特性， 将目标分为 高速运动目标和低速运动目标， 分别选择合理的处理算法进行小目标检测； 本方法首先对 原始回波数据进行距离维FFT处理、 MTD处理和舰速补偿后得。

13、到固定杂波背景， 针对高速运 动目标， 进行固定杂波对消后， 进行抑制地杂波、 海杂波抑制处理， 并对抑制后的数据进行 MTD处理和恒虚警处理， 完成对高速目标的检测； 针对低速目标则采用检测前跟踪处理和扫 描间相关处理完成杂波下低速小目标检测， 提高了海面小目标检测运行效率， 显著提高了 对海面小目标检测性能。 0006 下面结合附图对本发明做进一步详细的描述。 附图说明 0007 图1为本发明的流程图。 0008 图2为核密度估计实现恒虚警的流程图。 0009 图3为核密度估计算法二维参考单元选取方式示意图。 0010 图4为TBD处理的流程图。 0011 图5为扫描间相关处理流程图。 具。

14、体实施方式 0012 如图1所示， 一种对海雷达小目标检测方法， 具体包括： 说明书 1/5 页 5 CN 111610501 A 5 0013 步骤1： 对原始回波数据进行距离维FFT处理和MTD处理； 0014 步骤2： 对MTD处理后的数据进行舰速补偿； 取舰速补偿后的零多普勒通道数据作 为固定杂波背景显示。 0015 步骤3： 进行固定杂波对消以及地杂波、 海杂波抑制， 并对抑制后的数据进行MTD处 理和恒虚警处理， 完成对高速目标的检测， 所述高速目标为速度高于设定阈值的目标， 在本 实施例中， 设定阈值为3m/s。 0016 采用二阶数字对消器对MTD处理的数据进行固定杂波抑制， 。

15、具体为： 将相邻三个两 两相邻的回波信号进行相减， 所得结果再进行相减， 所述二阶数字对消器单位脉冲响应为： 0017 h3(t)h2(t)-h2(t-T) 0018 (t)-3 (t-T)+ (t-2T)- (t-3T) 0019 功率增益为 0020 0021 其差分方程为 0022 0023其中，v平台运动速度， 来自平台运动信息。 天线扫描角度， 来 自天线方位信息。 0024 判断固定杂波抑制后的数据中是否存在海杂波， 如果存在， 则采用自适应二次对 消器进行海杂波抑制， 其中， 判断固定杂波抑制后的数据中是否存在海杂波的具体方法为： 0025 (1)对固定杂波抑制后的数据进行求模；。

16、 0026 (2)计算求模结果的统计平均， 按恒虚警检测信噪比要求设置初始门限TH； 0027 (3)将求模结果与TH进行比较， 统计小于TH的模值， 计算小于TH模值的均值， 然后 乘以一个系数作为杂波定位的低门限TH1； 统计大于TH的模值， 计算大于TH模值的均值， 将 大于TH模值的均值作为杂波定位的高门限TH2； 0028 (4)设置杂波门限个数TH3， 设N为信号模值连续落在区间TH1,TH2的点数， 若N TH3则判决该区域存在海杂波， 若NTH3则判决该区域不存在海杂波。 0029 采用自适应二次对消器进行海杂波抑制， 具体步骤为： 0030 对判决为海杂波的杂波区数据进行分段。

17、处理， 每段的长度为N； 0031 对运动杂波区中不同的杂波段分别进行杂波的多普勒相位估计， 海杂波频率中心 ci的估计方法采用多滤波器组探索法， 具体为： 0032 假设杂波多普勒频率在0fr之间均匀分布， 取间隔： 0033 ffr/N 0034 设计N个自适应二次对消滤波器组： 0035 Hk(z)1-2exp(j2 kf/fr)z-1+exp(j4 kf/fr)z-2 k0,1,2,.N-1 0036 将杂波数据分别通过这些滤波器组， 其输出杂波剩余功率谱为： 0037 Sk(f)|Hk(exp(j2 f/fr)|2W(f) k0,1,2,.N-1 0038 设杂波的中心频率为fd， 。

18、第k个滤波器的杂波功率输出最小， 则可认为杂波的中心 说明书 2/5 页 6 CN 111610501 A 6 频率满足： 0039 (k-1/2)ffd(k+1/2)f 0040 此时自适应权系数为WIcos(2 k/N)， WQsin(2 k/N) 0041 谱中心的最大估计误差为： 0042 e(f)fr/2fr/(2N) 0043 凹口对准ci的二次数字对消器差分方程为 0044 0045 对地杂波和海杂波抑制后的数据进行MTD处理、 恒虚警处理对高速目标进行检测。 0046 如果目标运动速度大到在一个波束驻留时间内跨距离分辨单元， 则无法使用MTD 处理。 海面上航行的目标， 目前出。

19、现这样的径向速度一般是不可能的。 因此对于运动速度较 大的目标， 可以考虑采用一个波束驻留时间内的目标回波进行MTD处理来进行检测。 0047 本发明提出一种不需要先验假设的非参量恒虚警检测方法， 采用核密度估计方法 来准确估计杂波的分布函数， 确定检测门限。 核密度估计方法， 通过核密度函数来代替所需 估计的背景杂波概率密度函数， 再通过概率密度函数与检测门限的关系确定检测门限， 从 而实现良好的恒虚警检测。 如图2是用改进的核密度估计实现恒虚警的流程图， 这里采用二 维的方式来选取参考单元数， 如图3所示， 被检测单元D四周为参考单元， 共8个。 通过改进的 核密度估计方法实时的估计杂波的。

20、分布函数， 再利用分布函数与检测门限的关系直接将参 考单元的杂波数据代入， 从而实时的求得检测门限。 而不需要知道杂波的背景函数。 0048 该核密度估计的基本原理为： 依次以每一点x为中心取核函数的平均来求待估计 的概率密度函数， 可以用下式表示 0049 0050 其中， K为核函数， an定义为： 先选定一个大于1小于n的自然数kkn， 对任何xX， 以an(x)an(x； X1,Xn)记最小的正数a， 使x-a,x+a中至少包含样本X1,Xn中的k个。 常见的核函数有以下几种： 0051均匀核： 0052高斯核： 0053Epanechnikov核： 0054四次核： 0055 本发明。

21、中采用高斯核， 对概率密度函数积分得到分布函数的表达式为： 0056 0057 由雷达检测知识， 根据以上得到的分布函数的表达式， 知道检测门限的表达式为： 说明书 3/5 页 7 CN 111610501 A 7 0058 0059式中， Pfa为虚警概率，为的逆函数。 0060 则可得到虚警概率与检测门限的关系为： 0061 0062其中， erfc函数定义为 0063 通过上述研究， 对于未知情况下杂波的分布函数， 可以通过给出的核密度估计方 法准确的估计， 从而通过给定的虚警概率确定检测门限。 如接收信号经处理后， 检测结果高 于检测门限， 则认为检测到目标。 0064 对舰速补偿后的。

22、数据， 按其多普勒频率高低分别采用检测前跟踪处理和扫描间相 关处理完成杂波下低速小目标检测， 具体为： 0065 如图4所示， 对多普勒频率高于设定阈值的数据采用检测前跟踪处理检测小目标 并获得目标航迹， 即TBD处理： 0066 探测目标为海面低速小目标时， 在一个波束驻留时间内， 尽管通过去线性调频和 相参积累， 仍然会有目标回波信号在检测端与杂波和噪声的竞争中无法通过检测门限被检 测出来的情况。 面对这样的情况， 可以采用检测前跟踪的算法来进行处理， 在使用统计判决 确定有目标回波信号的同时， 给出目标的航迹。 0067 本步骤针对海面慢速运动目标的检测跟踪问题， 提出采用动态规划和Ka。

23、lman滤波 结合的方法进行TBD处理。 即根据目标运动模型， 在需要建立起Kalman滤波运动方程的波束 数内， 首先设置第一门限， 这一门限的虚警率较高， 以确保目标回波信号能通过这一门限。 当然因为虚警率较高， 有相当的杂波也会通过第一检测门限。 按照运动方程， 建立多个可能 的目标运动的滤波方程。 然后按照动态规划的方法， 将下一个波束驻留时间内通过第一检 测门限的信号与前面信号进行相参积累或者在目标运动的情况下， 相位关系因为运动被破 坏， 故进行非相参积累， 将最大值作为动态规划的阶段目标， 回波信号是否在Kalman滤波的 关联窗内作为这一阶段优化的约束条件。 按照这样的处理流程。

24、， 将不同波束驻留时间内积 累的信号能量与对应波束驻留数较低虚警率的检测门限进行比较来确定是否有目标， 并通 过约束条件来将杂波信号形成的虚假航迹过滤。 0068 对多普勒频率低于设定阈值的数据采用扫描间相关处理完成杂波下低速小目标 检测， 具体步骤为： 0069 当目标运动速度较慢时， 相对雷达处于静止或接近静止时， 对于小目标， 可以采用 扫描间积累来提高信杂比。 假设进行扫描间积累的扫描次数为N， 则要求能进行有效的积 累， 目标在这N个扫描期间， 不移动出一个距离单元， 并且因为目标运动造成的相位变化应 该在 /2之内， 则按照条件分别有NTwvr1.5。 如果要进行相参积累， 还需要。

25、有NTwvr/ / 2。 其中N为脉冲扫描次数， Tw为天线扫描周期， vr为目标相对雷达的径向速度， 为雷达工作 波长。 能进行扫描间相参积累对目标速度的限制较大， 几乎只能对相对没有运动的目标进 行。 因此扫描间积累主要采用非相参积累。 后面的描述中， 将每个扫描期间称为一帧。 说明书 4/5 页 8 CN 111610501 A 8 0070 由于海杂波的强度和形态会随着时间的变化其相关性不高， 而目标在帧间的相关 性较强， 则可以根据这样的特征处理强的杂波， 具体实现可以在帧间进行滑窗处理： 先存储 经过前面杂波抑制的多帧目标信息； 然后利用前几帧保存的数据与当前的目标信息在对应 的单元组成一个窗； 统计窗内对应单元目标出现的次数， 当目标出现的次数超过门限时则 认为该位置是目标而不是杂波， 如果在窗内的目标出现的次数没有达到门限则认为该位置 为杂波将其剔除掉， 处理过程如图5。 说明书 5/5 页 9 CN 111610501 A 9 图1 图2 图3 说明书附图 1/3 页 10 CN 111610501 A 10 图4 说明书附图 2/3 页 11 CN 111610501 A 11 图5 说明书附图 3/3 页 12 CN 111610501 A 12 。