CN110907907A - 一种海杂波多普勒谱特性分析对比方法 - Google Patents

一种海杂波多普勒谱特性分析对比方法 Download PDF

Info

Publication number
CN110907907A
CN110907907A CN201910996586.4A CN201910996586A CN110907907A CN 110907907 A CN110907907 A CN 110907907A CN 201910996586 A CN201910996586 A CN 201910996586A CN 110907907 A CN110907907 A CN 110907907A
Authority
CN
China
Prior art keywords
frequency shift
sea
sea clutter
wave
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910996586.4A
Other languages
English (en)
Other versions
CN110907907B (zh
Inventor
许心瑜
张金鹏
张玉石
李清亮
李慧明
张浙东
万晋通
余运超
黎鑫
尹雅磊
朱秀芹
尹志盈
赵鹏
夏晓云
李善斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Research Institute of Radio Wave Propagation CRIRP
Original Assignee
China Research Institute of Radio Wave Propagation CRIRP
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Research Institute of Radio Wave Propagation CRIRP filed Critical China Research Institute of Radio Wave Propagation CRIRP
Priority to CN201910996586.4A priority Critical patent/CN110907907B/zh
Publication of CN110907907A publication Critical patent/CN110907907A/zh
Application granted granted Critical
Publication of CN110907907B publication Critical patent/CN110907907B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/414Discriminating targets with respect to background clutter

Abstract

本发明公开了一种海杂波多普勒谱特性分析对比方法,包括如下步骤:步骤1,海杂波平均多普勒谱估计及频移谱宽计算;步骤2,船载平台的海杂波频移和谱宽修正;步骤3,海洋环境参数匹配及测试参数影响分析对比。本发明针对同一雷达在多平台测试的海杂波多普勒谱特性分析对比中存在的平台差异影响和测试参数条件复杂多样的现象,提供了一种海杂波多普勒谱特性分析对比方法,修正船载平台运动带来的影响,并按雷达参数和海洋环境参数影响的大小差异设计一定的分析次序和步骤,便于快速得到同一雷达在不同平台测试的多种海洋环境条件下的多普勒谱特性差异。

Description

一种海杂波多普勒谱特性分析对比方法
技术领域
本发明属于杂波多普勒特性研究领域,特别涉及该领域中的一种多平台测试的海杂波多普勒谱特性分析对比方法。
背景技术
海上雷达监视系统在工作中难免会受到来自海面的雷达回波的随机信号即海杂波的影响,时间相关性的分析有助于设计检测算法区分目标和海杂波。海杂波时间相关性通常称脉间相关性,常用杂波功率谱来描述。由于海面是运动的,海杂波谱将产生多普勒频移,所以也称为多普勒谱。在雷达动目标检测等处理中,为了有效抑制杂波,需要准确估计杂波多普勒谱中心位置和谱宽,但是受雷达参数和海洋参数影响,海杂波谱通常会呈现不同变化趋势,不能单一地用一种条件或几种条件概括。需要研究不同海洋环境下的海面运动特性对多普勒谱的影响。
为了获取更多海洋参数条件下的海杂波数据,一种是时间上的积累,进行长时间的观测,这种方式比较适合岸基固定平台,另一种是空间上的变化,在不同海域进行测试,但是有些海域受多种条件限制,很难找到合适的平台搭建地点进行长期观测,而搭载某些周期性航行大型船只,可进行多次短期测试。但是这样带来的问题就是不同平台模式下,即使雷达相同,也会由于平台运动状态的影响,引入附加的多普勒谱特性差异。在分析这种移动平台测试的海杂波谱特性时,需要尽量去除平台运动引入的一些非海洋环境参数的影响,避免在分析变化趋势时带来误差。
此外,由于不同测试环境下试验条件的差异,可能获取的海洋环境参数类型有一定差异,如何在不对等的环境参数中挑选对海杂波多普勒谱特性影响更突出的因素进行分析,需要根据理论分析和实测数据处理经验设计一定的分析对比原则,从而在多样的雷达参数和海洋环境参数组合中更迅速的找出多普勒谱特性变化规律。
发明内容
本发明所要解决的技术问题就是提供一种多平台测试的海杂波多普勒谱特性分析对比方法。
本发明采用如下技术方案:
一种海杂波多普勒谱特性分析对比方法,其改进之处在于,包括如下步骤:
步骤1,海杂波平均多普勒谱估计及频移谱宽计算:
对船载平台和固定平台下的雷达采集获取的海杂波数据,首先采用功率谱密度函数估计算法得到各组数据在每个距离门上的平均多普勒谱,然后根据每个距离门的平均多普勒谱的最大谱峰位置动态设置一定的区间范围,利用积分法估计多普勒谱的谱中心和谱展宽,得到频移和谱宽的初始值;
功率谱密度函数估计算法采用Welch或Burg谱估计方法,得到平均多普勒谱曲线Pf(fi),其中f为多普勒频率,i=1,2...N,N为数据点个数,在最高谱峰对应的频率值fm两边设定积分区间宽度[fm-fk,fm+fk],其中fk根据实测数据的平均多普勒曲线覆盖范围设定,然后由下式计算频移和谱宽的初始值fd和bw
Figure BDA0002239917330000021
Figure BDA0002239917330000022
步骤2,船载平台的海杂波频移和谱宽修正:
读取船载平台的姿态记录仪同步采集的经纬度值,按时间分段,根据每个时间段起始和结束的经纬度位置估算实时的船速大小和航行方向,与海杂波数据测试时间进行匹配,然后根据船体航行方向与雷达天线波束指向的几何关系,计算船体运动引起的频移和展宽值,对步骤1中的频移和谱宽的初始值进行修正;
船体运动引起的频移f和展宽值b由以下计算公式得到:
Figure BDA0002239917330000023
式中:vs为船速大小;α为船体航行方向与雷达天线波束指向的夹角;φaz为雷达天线方位波束宽度;λ为雷达波长;
通过下式对频移和谱宽的初始值fd和bw进行修正,得到船载平台运动修正后的频移fmd和谱宽bmw
Figure BDA0002239917330000024
步骤3,海洋环境参数匹配及测试参数影响分析对比:
以测试时间作为联系进行匹配,将海杂波数据与测试区域的海洋环境参数测试数据进行匹配,然后按照浪向一致性判读、随擦地角变化趋势分析及平稳区域确定、海情和浪向影响对比、距离分辨率影响对比的顺序依次分析海杂波频移和谱宽的变化规律,具体步骤如下:
步骤31,浪向一致性判读及海情等级划分:
计算雷达天线波束指向与浪向间的夹角,判断相对浪向是逆浪、顺浪还是侧浪,并增加浪向一致性参数φcst,结合步骤2计算得到的频移修正值进行判读,对于符合逆浪条件下频移为正或顺浪条件下频移为负的数据将其对应的φcst赋值为1,否则赋值为0,侧浪条件对应的数据统一将φcst赋值为1,在以下分析过程中,优先选择φcst为1的海杂波数据;根据道格拉斯海情等级表,结合对应浪高值估算海情,标注每组海杂波数据对应海情等级,海情等级与相对浪向两两组合形成L种参数组合以及每种参数组合下对应的数据集;
步骤32,随擦地角变化趋势分析及平稳区域确定:
根据雷达距离采样单元长度将海杂波数据中的距离门序号换算为实际距离,结合雷达架设高度计算每一个距离门对应的擦地角,分别绘制L种参数组合下频移和谱宽随擦地角变化的散点分布图,然后利用曲线拟合分别得到L种参数组合下的频移和谱宽随擦地角的变化曲线,对每一种参数组合下的多组数据的频移值和谱宽值求平均,得到该参数组合下频移和谱宽均值随擦地角的变化趋势,截取曲线平稳区域对应的擦地角区间;
步骤33,海情和浪向影响对比:
首先求取需要分析的所有参数组合下的擦地角区间的交集,提取各组海杂波数据对应该擦地角区间交集内的频移和谱宽计算值,按海情等级和相对浪向形成的L种参数组合分类,将每种参数组合下的频移和谱宽计算值分别进行合并,绘制在逆浪、顺浪、侧浪情况下频移和谱宽随海情等级变化的散点分布图,对同一海情同一相对浪向的数据求取频移和谱宽的均值,记录均值随海情和浪向的变化趋势;
步骤34,距离分辨率影响对比:
对于L种参数组合的每一种组合,若具有不同距离分辨率的海杂波数据,在擦地角区间交集中,分别将每种距离分辨率下对应的频移值和谱宽值在距离向求自相关函数,并计算相关长度比较大小,然后统计每种距离分辨率下的谱宽值的概率密度分布,并分别采用高斯分布和伽马分布进行拟合,比较更接近哪一种分布,得到不同距离分辨率下海杂波在距离向的多普勒谱特性变化剧烈程度对比结果。
本发明的有益效果是:
本发明针对同一雷达在多平台测试的海杂波多普勒谱特性分析对比中存在的平台差异影响和测试参数条件复杂多样的现象,提供了一种海杂波多普勒谱特性分析对比方法,修正船载平台运动带来的影响,并按雷达参数和海洋环境参数影响的大小差异设计一定的分析次序和步骤,便于快速得到同一雷达在不同平台测试的多种海洋环境条件下的多普勒谱特性差异。
与现有技术相比,本发明所公开海杂波多普勒谱特性分析对比方法的有益效果,一是通过修正不同海杂波测量平台的运动影响,扩大了数据来源,便于分析对比多种测试参数下的海杂波多普勒谱特性,二是在海杂波数据对应的海洋环境参数复杂多样的情况下,根据理论分析和实测数据分析经验设计了参数影响的分析对比次序和步骤,便于快速有序地找出多种海洋环境条件下的海杂波多普勒谱特性变化规律,为进一步深入分析指明了方向。
附图说明
图1是本发明实施例1所公开方法的流程示意图;
图2(a)是S波段5级海情的逆浪、顺浪和侧浪条件下的频移随擦地角的变化曲线及曲线拟合结果示意图;
图2(b)是S波段5级海情的逆浪、顺浪和侧浪条件下的谱宽随擦地角的变化曲线及曲线拟合结果示意图;
图3(a)是逆浪条件下2到5级海情的海杂波频移的散点分布及均值计算结果示意图;
图3(b)是逆浪条件下2到5级海情的海杂波谱宽的散点分布及均值计算结果示意图;
图4(a)是4级海情逆浪条件下四种分辨率对应的频移的距离相关函数曲线示意图;
图4(b)是4级海情逆浪条件下四种分辨率对应的谱宽的距离相关函数曲线示意图;
图5(a)是4级海情逆浪条件下60m分辨率的海杂波数据谱宽的概率密度函数及采用高斯分布和伽马分布的拟合曲线示意图;
图5(b)是4级海情逆浪条件下30m分辨率的海杂波数据谱宽的概率密度函数及采用高斯分布和伽马分布的拟合曲线示意图;
图5(c)是4级海情逆浪条件下15m分辨率的海杂波数据谱宽的概率密度函数及采用高斯分布和伽马分布的拟合曲线示意图;
图5(d)是4级海情逆浪条件下7.5m分辨率的海杂波数据谱宽的概率密度函数及采用高斯分布和伽马分布的拟合曲线示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图和实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例1,如图1所示,本实施例公开了一种海杂波多普勒谱特性分析对比方法,包括如下步骤:
步骤1,海杂波平均多普勒谱估计及频移谱宽计算:
对船载平台和固定平台下的雷达采集获取的海杂波数据,首先采用功率谱密度函数估计算法得到各组数据在每个距离门上的平均多普勒谱,然后根据每个距离门的平均多普勒谱的最大谱峰位置动态设置一定的区间范围,利用积分法估计多普勒谱的谱中心和谱展宽,得到频移和谱宽的初始值;
功率谱密度函数估计算法采用Welch或Burg谱估计方法,得到平均多普勒谱曲线Pf(fi),其中f为多普勒频率,i=1,2...N,N为数据点个数,在最高谱峰对应的频率值fm两边设定积分区间宽度[fm-fk,fm+fk],其中fk根据实测数据的平均多普勒曲线覆盖范围设定,然后由下式计算频移和谱宽的初始值fd和bw
Figure BDA0002239917330000051
Figure BDA0002239917330000052
步骤2,船载平台的海杂波频移和谱宽修正:
读取船载平台的姿态记录仪同步采集的经纬度值,按时间分段,根据每个时间段起始和结束的经纬度位置估算实时的船速大小和航行方向,与海杂波数据测试时间进行匹配,然后根据船体航行方向与雷达天线波束指向的几何关系,计算船体运动引起的频移和展宽值,对步骤1中的频移和谱宽的初始值进行修正;
船体运动引起的频移f和展宽值b由以下计算公式得到:
Figure BDA0002239917330000053
式中:vs为船速大小;α为船体航行方向与雷达天线波束指向的夹角;φaz为雷达天线方位波束宽度;λ为雷达波长;
通过下式对频移和谱宽的初始值fd和bw进行修正,得到船载平台运动修正后的频移fmd和谱宽bmw
Figure BDA0002239917330000061
步骤3,海洋环境参数匹配及测试参数影响分析对比:
以测试时间作为联系进行匹配,将海杂波数据与测试区域的海洋环境参数测试数据进行匹配,然后按照浪向一致性判读、随擦地角变化趋势分析及平稳区域确定、海情和浪向影响对比、距离分辨率影响对比的顺序依次分析海杂波频移和谱宽的变化规律,具体步骤如下:
步骤31,浪向一致性判读及海情等级划分:
计算雷达天线波束指向与浪向间的夹角,判断相对浪向是逆浪、顺浪还是侧浪,并增加浪向一致性参数φcst,结合步骤2计算得到的频移修正值进行判读,对于符合逆浪条件下频移为正或顺浪条件下频移为负的数据将其对应的φcst赋值为1,否则赋值为0,侧浪条件对应的数据统一将φcst赋值为1,在以下分析过程中,优先选择φcst为1的海杂波数据;根据道格拉斯海情等级表,结合对应浪高值估算海情,标注每组海杂波数据对应海情等级,海情等级与相对浪向两两组合形成L种参数组合以及每种参数组合下对应的数据集;
步骤32,随擦地角变化趋势分析及平稳区域确定:
根据雷达距离采样单元长度将海杂波数据中的距离门序号换算为实际距离,结合雷达架设高度计算每一个距离门对应的擦地角,分别绘制L种参数组合下频移和谱宽随擦地角变化的散点分布图,然后利用曲线拟合分别得到L种参数组合下的频移和谱宽随擦地角的变化曲线,对每一种参数组合下的多组数据的频移值和谱宽值求平均,得到该参数组合下频移和谱宽均值随擦地角的变化趋势,截取曲线平稳区域对应的擦地角区间;
步骤33,海情和浪向影响对比:
首先求取需要分析的所有参数组合下的擦地角区间的交集,提取各组海杂波数据对应该擦地角区间交集内的频移和谱宽计算值,按海情等级和相对浪向形成的L种参数组合分类,将每种参数组合下的频移和谱宽计算值分别进行合并,绘制在逆浪、顺浪、侧浪情况下频移和谱宽随海情等级变化的散点分布图,对同一海情同一相对浪向的数据求取频移和谱宽的均值,记录均值随海情和浪向的变化趋势;
步骤34,距离分辨率影响对比:
对于L种参数组合的每一种组合,若具有不同距离分辨率的海杂波数据,在擦地角区间交集中,分别将每种距离分辨率下对应的频移值和谱宽值在距离向求自相关函数,并计算相关长度比较大小,然后统计每种距离分辨率下的谱宽值的概率密度分布,并分别采用高斯分布和伽马分布进行拟合,比较更接近哪一种分布,得到不同距离分辨率下海杂波在距离向的多普勒谱特性变化剧烈程度对比结果。
以S波段某雷达海杂波数据为例,分别搭载了船载平台和岸基平台进行测试,按照本实施例步骤对海杂波数据进行以下处理:
步骤1,海杂波平均多普勒谱估计及频移谱宽计算:
基于获取的海杂波数据,首先采用Welch估计算法得到各组数据在每个距离门上的平均多普勒谱,然后搜索每个距离门的平均多普勒谱的最大值对应的频率值fmax,结合S波段实测数据多普勒谱曲线分布特点,设置fk=200Hz,积分区间则为[fm-200,fm+200],然后利用积分法估计得到频移和谱宽的初始值fd和bw
步骤2,船载平台的海杂波频移和谱宽修正:
读取船载平台的姿态记录仪在海杂波测试时同步采集的经纬度值,本实施例按30秒间隔分段,提取每个时间段开始和结束处的船体所在的经纬度值,估算船速大小vs和航行方向
Figure BDA0002239917330000073
将其与雷达海杂波数据按时间信息进行匹配,读取该组数据雷达参数中的雷达天线波束指向即方位角,统一以正北方向为0°,计算雷达天线波束指向与船体航行方向的夹角,然后采用下式计算船体运动引起的在海面自身运动基础上附加的频移f和展宽值b
Figure BDA0002239917330000071
式中:α为船体航行方向
Figure BDA0002239917330000074
与雷达天线波束指向的夹角;φaz为雷达天线方位波束宽度;λ为雷达波长。进而对步骤1中的fd和bw进行修正,得到修正后的频移fmd和谱宽bmw
Figure BDA0002239917330000072
步骤3,海洋环境参数匹配及测试参数影响分析对比:
以测试时间作为联系进行匹配,本实施例中首先将海杂波测试时间和海洋环境参数记录时间统一采用“年+月+日+时+分”的形式记录,对于每组海杂波数据,寻找与其时间相隔最近且小于一定的预设时间间隔阈值的一组海洋环境参数,依次循环完成每组海杂波数据与海洋环境参数的匹配,然后按照以下步骤分析测试参数对海杂波谱特性的影响:
步骤31,浪向一致性判读及海情等级划分:
计算每组数据雷达天线波束指向与浪向间的夹角,判断相对浪向是逆浪、顺浪还是侧浪,并增加浪向一致性参数φcst,结合步骤2计算得到的频移修正值fmd进行判读,一般逆浪条件下频移为正,顺浪条件下频移为负,对于符合以上情况的数据将其对应的φcst赋值为1,否则赋值为0,侧浪条件对应的数据统一将φcst赋值为1,在以下分析过程中,优先选择φcst为1的海杂波数据;根据道格拉斯海情等级表,结合对应浪高值估算海情,标注每组海杂波数据对应海情等级,海情等级与相对浪向两两组合形成L种参数组合,每种参数组合下对应一个数据集。本实施例中海情等级和相对浪向情况如下:
海情等级:2、3、4、5级;相对浪向:逆浪、顺浪、侧浪,且在每一种海情下均有三种浪向的海杂波数据;因此共计12种参数组合。
步骤32,随擦地角变化趋势分析及平稳区域确定:
根据雷达距离采样单元长度将海杂波数据中的距离门序号换算为实际距离,本实施例中需首先找出该S波段雷达泄漏信号所在位置,将其作为第一个距离门,然后将变换后的距离门序号乘以距离采样长度60m,得到每个距离门和雷达发射天线的真实距离R,结合雷达架设高度H,利用arcsin(H/R)得到每个距离门对应的擦地角。分别绘制12种参数组合下频移和谱宽随擦地角变化的散点分布图,然后进行曲线拟合,本实施例中频移散点图采用多项式模型进行拟合,谱宽散点图采用指数模型进行拟合,如图2(a),2(b)所示为S波段5级海情的逆浪、顺浪和侧浪条件三种参数组合下的海杂波数据的频移和谱宽随擦地角的散点分布图及曲线拟合结果。对每一种参数组合下的多组数据的频移值和谱宽值求平均,得到该参数组合下频移均值和谱宽均值随擦地角的变化趋势,截取曲线平稳区域对应的擦地角区间,对于不同参数组合,情况会有所差异,有时在整个区间内比较平缓,可作为一个区间,有时根据变化趋势急缓程度分为几个区间。
步骤33,海情和浪向影响对比:
首先求取需要分析的所有参数组合下的擦地角区间的交集[θ12],对于本实施例中的12种参数组合,这个擦地角交集对于频移值为[0.5°,1°],对于谱宽值为[0.7°,1°]。分别提取各组海杂波数据对应其擦地角交集内的频移和谱宽计算值,按12种参数组合对数据分类,将每类下所有数据的频移和谱宽计算值分别进行合并,绘制在逆浪、顺浪、侧浪情况下频移和谱宽随海情等级变化的散点分布图,对同一海情同一浪向的数据求取频移和谱宽的均值,记录均值随海情和浪向的变化趋势。如图3(a),3(b)所示为逆浪条件下2到5级海情的海杂波频移和谱宽的散点分布图,图中“□”为频移和谱宽的均值。可以清晰看到,随着海情的升高,频移和谱宽逐步增高。
步骤34,距离分辨率影响对比:
对于本实施例的12种参数组合,均具有四种不同距离分辨率(60m,30m,15m,7.5m)的测试数据,在某一参数组合下,提取各组数据擦地角区间交集[θ12]内的频移值和谱宽值,在距离向求取自相关函数,并计算相关长度比较大小,如图4(a),4(b)所示为4级海情逆浪条件下四种分辨率对应的频移和谱宽的距离相关函数,频移的距离相关长度分别为75m、45m、22.5m和15m,谱宽的距离自相关长度分别为60m、37.5m、22.5m和15m,相关长度越小,距离向的频移和谱宽变化程度越剧烈。
然后统计每种距离分辨率下的谱宽值的概率密度分布,并分别采用高斯分布和伽马分布进行拟合,比较更接近哪一种分布,如图5(a),5(b),5(c),5(d)所示为4级海情逆浪条件下四种分辨率的海杂波频移和谱宽与高斯分布和伽马分布的拟合结果,由图可见,在4级海情逆浪条件下,整体来看四种分辨率都更接近伽马分布,而且分辨率越高越接近伽马分布,说明随着分辨率的提高不同距离门之间的变化剧烈程度会越大。

Claims (3)

1.一种海杂波多普勒谱特性分析对比方法,其特征在于,包括如下步骤:
步骤1,海杂波平均多普勒谱估计及频移谱宽计算:
对船载平台和固定平台下的雷达采集获取的海杂波数据,首先采用功率谱密度函数估计算法得到各组数据在每个距离门上的平均多普勒谱,然后根据每个距离门的平均多普勒谱的最大谱峰位置动态设置一定的区间范围,利用积分法估计多普勒谱的谱中心和谱展宽,得到频移和谱宽的初始值;
步骤2,船载平台的海杂波频移和谱宽修正:
读取船载平台的姿态记录仪同步采集的经纬度值,按时间分段,根据每个时间段起始和结束的经纬度位置估算实时的船速大小和航行方向,与海杂波数据测试时间进行匹配,然后根据船体航行方向与雷达天线波束指向的几何关系,计算船体运动引起的频移和展宽值,对步骤1中的频移和谱宽的初始值进行修正;
步骤3,海洋环境参数匹配及测试参数影响分析对比:
以测试时间作为联系进行匹配,将海杂波数据与测试区域的海洋环境参数测试数据进行匹配,然后按照浪向一致性判读、随擦地角变化趋势分析及平稳区域确定、海情和浪向影响对比、距离分辨率影响对比的顺序依次分析海杂波频移和谱宽的变化规律,具体步骤如下:
步骤31,浪向一致性判读及海情等级划分:
计算雷达天线波束指向与浪向间的夹角,判断相对浪向是逆浪、顺浪还是侧浪,并增加浪向一致性参数φcst,结合步骤2计算得到的频移修正值进行判读,对于符合逆浪条件下频移为正或顺浪条件下频移为负的数据将其对应的φcst赋值为1,否则赋值为0,侧浪条件对应的数据统一将φcst赋值为1,在以下分析过程中,优先选择φcst为1的海杂波数据;根据道格拉斯海情等级表,结合对应浪高值估算海情,标注每组海杂波数据对应海情等级,海情等级与相对浪向两两组合形成L种参数组合以及每种参数组合下对应的数据集;
步骤32,随擦地角变化趋势分析及平稳区域确定:
根据雷达距离采样单元长度将海杂波数据中的距离门序号换算为实际距离,结合雷达架设高度计算每一个距离门对应的擦地角,分别绘制L种参数组合下频移和谱宽随擦地角变化的散点分布图,然后利用曲线拟合分别得到L种参数组合下的频移和谱宽随擦地角的变化曲线,对每一种参数组合下的多组数据的频移值和谱宽值求平均,得到该参数组合下频移和谱宽均值随擦地角的变化趋势,截取曲线平稳区域对应的擦地角区间;
步骤33,海情和浪向影响对比:
首先求取需要分析的所有参数组合下的擦地角区间的交集,提取各组海杂波数据对应该擦地角区间交集内的频移和谱宽计算值,按海情等级和相对浪向形成的L种参数组合分类,将每种参数组合下的频移和谱宽计算值分别进行合并,绘制在逆浪、顺浪、侧浪情况下频移和谱宽随海情等级变化的散点分布图,对同一海情同一相对浪向的数据求取频移和谱宽的均值,记录均值随海情和浪向的变化趋势;
步骤34,距离分辨率影响对比:
对于L种参数组合的每一种组合,若具有不同距离分辨率的海杂波数据,在擦地角区间交集中,分别将每种距离分辨率下对应的频移值和谱宽值在距离向求自相关函数,并计算相关长度比较大小,然后统计每种距离分辨率下的谱宽值的概率密度分布,并分别采用高斯分布和伽马分布进行拟合,比较更接近哪一种分布,得到不同距离分辨率下海杂波在距离向的多普勒谱特性变化剧烈程度对比结果。
2.根据权利要求1所述的海杂波多普勒谱特性分析对比方法,其特征在于:步骤1中的功率谱密度函数估计算法采用Welch或Burg谱估计方法,得到平均多普勒谱曲线Pf(fi),其中f为多普勒频率,i=1,2...N,N为数据点个数,在最高谱峰对应的频率值fm两边设定积分区间宽度[fm-fk,fm+fk],其中fk根据实测数据的平均多普勒曲线覆盖范围设定,然后由下式计算频移和谱宽的初始值fd和bw
Figure FDA0002239917320000021
Figure FDA0002239917320000022
3.根据权利要求1所述的海杂波多普勒谱特性分析对比方法,其特征在于:步骤2中的船体运动引起的频移f和展宽值b由以下计算公式得到:
Figure FDA0002239917320000023
式中:vs为船速大小;α为船体航行方向与雷达天线波束指向的夹角;φaz为雷达天线方位波束宽度;λ为雷达波长;
通过下式对频移和谱宽的初始值fd和bw进行修正,得到船载平台运动修正后的频移fmd和谱宽bmw
Figure FDA0002239917320000031
CN201910996586.4A 2019-10-19 2019-10-19 一种海杂波多普勒谱特性分析对比方法 Active CN110907907B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910996586.4A CN110907907B (zh) 2019-10-19 2019-10-19 一种海杂波多普勒谱特性分析对比方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910996586.4A CN110907907B (zh) 2019-10-19 2019-10-19 一种海杂波多普勒谱特性分析对比方法

Publications (2)

Publication Number Publication Date
CN110907907A true CN110907907A (zh) 2020-03-24
CN110907907B CN110907907B (zh) 2022-06-14

Family

ID=69815475

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910996586.4A Active CN110907907B (zh) 2019-10-19 2019-10-19 一种海杂波多普勒谱特性分析对比方法

Country Status (1)

Country Link
CN (1) CN110907907B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112068087A (zh) * 2020-10-17 2020-12-11 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种岸基多通道雷达仿机载海杂波测量方法
CN112068086A (zh) * 2020-10-17 2020-12-11 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于外定标试验数据的岸基多通道雷达幅相校正方法
CN112147587A (zh) * 2020-09-28 2020-12-29 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种雷达波束方位中心海上标定方法
CN112213709A (zh) * 2020-10-17 2021-01-12 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于相参脉冲串的仿机载数据分组方法
CN112213708A (zh) * 2020-10-14 2021-01-12 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种岸基海杂波观测方法
CN112255600A (zh) * 2020-10-17 2021-01-22 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于收发共置多通道雷达仿机载运动方法
CN112255601A (zh) * 2020-10-17 2021-01-22 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种岸基多通道雷达仿机载数据诊断方法
CN112415478A (zh) * 2020-12-16 2021-02-26 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种凝视状态下窄波束雷达杂波的数据处理方法
CN112526464A (zh) * 2020-11-14 2021-03-19 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于实测数据估计岸基多通道雷达方位波束宽度的方法
CN112782690A (zh) * 2021-01-14 2021-05-11 中国科学院国家空间科学中心 星载雷达高度计的近岸非海洋波形检测分类方法及系统
CN114488107A (zh) * 2022-04-13 2022-05-13 南方海洋科学与工程广东省实验室(广州) 一种海杂波时空分布及影响分级产品制作的方法及装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102540162A (zh) * 2011-12-12 2012-07-04 中国船舶重工集团公司第七二四研究所 基于海杂波的低空电磁波传播特性估计方法
CN105259537A (zh) * 2015-11-10 2016-01-20 武汉大学 基于频移迭代的多普勒谱中心频率估计方法
CN105388465A (zh) * 2015-12-17 2016-03-09 西安电子科技大学 基于海浪谱模型的海杂波模拟方法
CN107678003A (zh) * 2017-09-15 2018-02-09 国家海洋局第海洋研究所 一种地波雷达海杂波背景下的目标检测方法和装置
CN111337895A (zh) * 2020-01-14 2020-06-26 北京理工大学 多通道海杂波空时相关性分析方法
WO2020136258A1 (fr) * 2018-12-27 2020-07-02 Thales Dispositif de génération d'un ensemble de données simulées du fouillis de mer, procédé et programme d'ordinateur associés

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102540162A (zh) * 2011-12-12 2012-07-04 中国船舶重工集团公司第七二四研究所 基于海杂波的低空电磁波传播特性估计方法
CN105259537A (zh) * 2015-11-10 2016-01-20 武汉大学 基于频移迭代的多普勒谱中心频率估计方法
CN105388465A (zh) * 2015-12-17 2016-03-09 西安电子科技大学 基于海浪谱模型的海杂波模拟方法
CN107678003A (zh) * 2017-09-15 2018-02-09 国家海洋局第海洋研究所 一种地波雷达海杂波背景下的目标检测方法和装置
WO2020136258A1 (fr) * 2018-12-27 2020-07-02 Thales Dispositif de génération d'un ensemble de données simulées du fouillis de mer, procédé et programme d'ordinateur associés
CN111337895A (zh) * 2020-01-14 2020-06-26 北京理工大学 多通道海杂波空时相关性分析方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张玉石 等: "海杂波测量定标的姿态修正数据处理方法", 《电子与信息学报》 *
张金鹏 等: "基于不同散射机制特征的海杂波时变多普勒谱模型", 《物理学报》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112147587A (zh) * 2020-09-28 2020-12-29 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种雷达波束方位中心海上标定方法
CN112147587B (zh) * 2020-09-28 2022-02-25 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种雷达波束方位中心海上标定方法
CN112213708A (zh) * 2020-10-14 2021-01-12 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种岸基海杂波观测方法
CN112213708B (zh) * 2020-10-14 2022-04-22 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种岸基海杂波观测方法
CN112068086B (zh) * 2020-10-17 2022-03-01 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于外定标试验数据的岸基多通道雷达幅相校正方法
CN112213709B (zh) * 2020-10-17 2022-03-04 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于相参脉冲串的仿机载数据分组方法
CN112255601A (zh) * 2020-10-17 2021-01-22 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种岸基多通道雷达仿机载数据诊断方法
CN112068086A (zh) * 2020-10-17 2020-12-11 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于外定标试验数据的岸基多通道雷达幅相校正方法
CN112255600B (zh) * 2020-10-17 2022-03-22 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于收发共置多通道雷达仿机载运动方法
CN112255600A (zh) * 2020-10-17 2021-01-22 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于收发共置多通道雷达仿机载运动方法
CN112255601B (zh) * 2020-10-17 2022-02-01 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种岸基多通道雷达仿机载数据诊断方法
CN112213709A (zh) * 2020-10-17 2021-01-12 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于相参脉冲串的仿机载数据分组方法
CN112068087A (zh) * 2020-10-17 2020-12-11 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种岸基多通道雷达仿机载海杂波测量方法
CN112068087B (zh) * 2020-10-17 2022-03-01 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种岸基多通道雷达仿机载海杂波测量方法
CN112526464A (zh) * 2020-11-14 2021-03-19 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于实测数据估计岸基多通道雷达方位波束宽度的方法
CN112415478A (zh) * 2020-12-16 2021-02-26 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种凝视状态下窄波束雷达杂波的数据处理方法
CN112782690A (zh) * 2021-01-14 2021-05-11 中国科学院国家空间科学中心 星载雷达高度计的近岸非海洋波形检测分类方法及系统
CN112782690B (zh) * 2021-01-14 2023-10-13 中国科学院国家空间科学中心 星载雷达高度计的近岸非海洋波形检测分类方法及系统
CN114488107A (zh) * 2022-04-13 2022-05-13 南方海洋科学与工程广东省实验室(广州) 一种海杂波时空分布及影响分级产品制作的方法及装置

Also Published As

Publication number Publication date
CN110907907B (zh) 2022-06-14

Similar Documents

Publication Publication Date Title
CN110907907B (zh) 一种海杂波多普勒谱特性分析对比方法
Lund et al. Marine radar ocean wave retrieval’s dependency on range and azimuth
CN110208806B (zh) 一种航海雷达图像降雨识别方法
CN110837078B (zh) 基于相关特性的阵列地波雷达海杂波背景下目标检测方法
JP4825183B2 (ja) 表層流速推定方法、装置並びにプログラム
CN103308910B (zh) 利用高频地波雷达检测海上非航行状态船只目标的方法
Xie et al. Measuring ocean surface wind field using shipborne high-frequency surface wave radar
CN110398722A (zh) 基于随机矩阵有限谱的扩展目标回波检测方法
Lipa et al. Developments in compact HF-radar ocean wave measurement
CN111751797B (zh) 一种基于方位角的高频地波雷达一阶和二阶回波谱分界线确定方法
EP2278355A1 (en) Radar system
CN113514833B (zh) 一种基于海浪图像的海面任意点波向的反演方法
CN112924959B (zh) 基于带宽合成的分布式高频地波雷达目标定位和测速方法
CN111796288B (zh) 一种基于杂波频谱补偿技术的三坐标雷达动目标处理方法
CN113820680B (zh) 一种基于协方差的多帧海陆雷达回波分割方法
CN112836638B (zh) 高海况sar船只检测方法及应用
Rino et al. X-band low-grazing-angle ocean backscatter obtained during LOGAN 1993
CN113687321B (zh) 雷达目标探测距离评估方法和装置
Crisp et al. Comparison of ship detectors for polarimetric SAR imagery
CN113093174B (zh) 基于phd滤波雷达起伏微弱多目标的检测前跟踪方法
JP2005140687A (ja) 目標類別装置
KR20180130282A (ko) 해류 정보 측정 방법과 프로그램 및 이를 이용한 해양 레이더 시스템
Srivastava et al. Analysis of Sea Waves and Ship Wake Detection
Wang et al. Elimination of the impact of vessels on ocean wave height inversion with X-band wave monitoring radar
Wang et al. Seafloor terrain detection from acoustic images utilizing the fast two-dimensional CMLD-CFAR

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant