CN107561527A - 一种舰载雷达海目标航向航速高精度补偿计算方法 - Google Patents

一种舰载雷达海目标航向航速高精度补偿计算方法 Download PDF

Info

Publication number
CN107561527A
CN107561527A CN201710621893.5A CN201710621893A CN107561527A CN 107561527 A CN107561527 A CN 107561527A CN 201710621893 A CN201710621893 A CN 201710621893A CN 107561527 A CN107561527 A CN 107561527A
Authority
CN
China
Prior art keywords
ship
compensation
speed
point mark
plane
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
CN201710621893.5A
Other languages
English (en)
Other versions
CN107561527B (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.)
724th Research Institute of CSIC
Original Assignee
724th Research Institute of CSIC
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 724th Research Institute of CSIC filed Critical 724th Research Institute of CSIC
Priority to CN201710621893.5A priority Critical patent/CN107561527B/zh
Publication of CN107561527A publication Critical patent/CN107561527A/zh
Application granted granted Critical
Publication of CN107561527B publication Critical patent/CN107561527B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)
  • Navigation (AREA)

Abstract

一种舰载雷达海目标航向航速高精度补偿计算方法,本发明涉及舰载雷达海目标航向航速计算,主要是针对在舰载雷达中本舰运动情况下的海目标的航向航速计算。舰载雷达观测到海目标形成了历史关联点迹,利用本舰航向航速信息将历史关联点迹由相对位置补偿为绝对位置,再采用最小二乘法,可以计算出海目标的绝对航向航速。对于某一个历史周期的关联点迹,补偿原则是使用从当前时间到该历史周期的时间段内所有本舰航向航速信息,因此,历史周期越久的航迹点,对应补偿的本舰航向航速信息也多。经分析,在本舰机动情况下由补偿带来的每秒位置误差为0.0874m,远远小于传统算法的误差0.8721m,大大提高了舰载雷达海目标航向航速计算的精确性。

Description

一种舰载雷达海目标航向航速高精度补偿计算方法
技术领域
本发明属于雷达数据处理技术。
背景技术
随着雷达技术的发展,目标跟踪算法的不断研究更新,为了实现主动探测系统对目标高精度的实时跟踪,其任务不仅是准确测量目标的位置,而且还包括精确计算目标的航向航速。目前舰载雷达由于其本舰运动的特点,实际是匀速运动的海目标,由雷达相对看来,其目标的运动形式较为复杂,往往是非线性变化的。这是因为雷达本身可能在做非线性变化,目标的运动状态是由雷达运动状态和目标实际运动状态合成导致的。传统的方法只关注滤波方法的选择,而未对可利用的资源作全局分析取舍,这样不能很好的抓住问题的本质,只能产生事倍功半的结果,最典型的计算方法是通过先计算相对航向航速再与本舰航向航速合成得到,其中对相对航向航速的估计方法,采用kalman滤波或alpha-beta滤波等算法,或在此基础上做的适应性改进,其研究的滤波算法虽然能在一定程度上达到效果,但并未能满足实际使用的需求,传统方法已然不适用。
本发明提出了通过将本舰航向航速信息补偿目标相对位置,再通过最小二乘法计算目标绝对航向航速,这样通过本舰航向航速信息补偿计算出来的结果误差较小,再经过最小二乘法滤波可以进一步减少由于目标位置突变带来的误差,经理论分析和实践检验,方法可行,目前在雷达目标跟踪航向航速参数计算应用方面有着良好的应用前景。
发明内容
实现本发明目的的技术解决方案为:舰载雷达观测到海目标形成了历史关联点迹,根据历史关联点迹和最小二乘法原理,可以得到海目标的相对航向航速;而在舰上,由于雷达平台是运动的,求目标的绝对航向航速与本舰航向航速有关。本发明将本舰航向航速先补偿目标历史关联点迹得到目标绝对位置再计算目标的绝对航向航速,而不是通过目标历史关联点迹先计算目标相对航向航速再叠加本舰航向航速的方法。历史周期越久的关联点迹,对应补偿的本舰航向航速信息也多,从而积累的误差也大。即历史周期为N的关联点迹,补偿原则是从当前时间到该历史周期时间段内的所有本舰航向航速信息。
本发明的显著优点为:目标绝对航向航速精度更高,能适应本舰的机动运动,计算成本低,其工程实现简单。本发明研究设计的这种新的数据处理技术能更好的满足作战使用要求,具有很高的推广应用价值。
附图说明
图1是本舰在1s内实际走过的路线示例图;
图2是补偿间隔为1s下不同历史周期的关联点迹对应的本舰运动误差Es;
图3是补偿间隔为1s下不同历史周期的关联点迹对应的补偿误差Ea;
图4是补偿间隔为10s下不同历史周期的关联点迹对应的本舰运动误差Es;
图5是补偿间隔为10s下不同历史周期的关联点迹对应的补偿误差Ea;
图6是方法二实现流程图。
具体实施方式
下面结合实施例和附图对本发明技术方案具体说明。
舰载雷达观测海目标,形成了历史关联点迹,设由历史关联点迹分解到x,y上的分量和时间t构成的数据为(xi,yi,ti,i=1,2,…,n),i=0为当前周期关联点迹,i=1为上个周期关联点迹,依次类推,i越大历史越久。由于该雷达观测的是海目标,其数据率ft不高,本发明实施例设为10s。
在舰上,由于雷达平台是运动的,设本舰历史航向航速分解到x,y上的分量和时间t构成的数据为(Vxj,Vyj,tj,j=0,1,2,…,m),j=0为当前时刻本舰航向航速信息,j=1为上个时刻本舰航向航速信息,依次类推,j越大历史越久。该数据率fs相对于ft要高,本发明实施例设为1s。
由于历史关联点迹是相当于本舰的相对点迹,根据它计算出来的航向航速是相对量,如果能把它直接转变为绝对量,就可以直接得到绝对航向航速结果。本发明实施例使用本舰航向航速信息和目标历史关联点迹采用如下公式,得到目标绝对位置分别为:
该转换过程称之为本舰航向航速信息对历史关联点迹作补偿。由此可见,历史周期越久的关联点迹,对应补偿的本舰航向航速信息也多。目标航向航速的计算方法采用了最小二乘法,需要若干周期(本发明实施例假设为N个周期)的历史关联点迹,因此转换的绝对量与该最近的N个历史关联点迹有关,而与更早的历史关联点迹无关。根据最小二乘法原理求得经验公式分别为x=axt+bx和y=ayt+by,将补偿后的目标绝对位置代入经验公式用最小二乘法计算可得到目标绝对航向航速。有:
sx=∑[xxi-(axti+bx)]2=min (3)
sy=∑[yyi-(ayti+by)]2=min (4)
可求得ax,ay,则绝对航向c=arc tan(ay/ax),绝对航速
由于本发明实施例1s使用一次本舰航向航速信息,因此,补偿得到的目标绝对位置存在误差,其误差与本舰航向航速信息有关。对于本舰机动情况下,该误差表现的更加突出。传统的方法(下面简称方法一)是先求目标的相对航向航速,再叠加本舰的航向航速,其补偿数据率更低,因此误差较大。本发明实施例方法(下面简称方法二)是将本舰航向航速先补偿得到的目标绝对位置再计算目标的绝对航向航速,其补偿数据率较高,因此误差较小。下面理论计算并分析比较两种方法的误差大小。
在上述假设的历史关联点迹和本舰航向航速信息的基础上,设本舰转向时航速为dv=10m/s,航向变化率dc为1度/s,如图1中所示S0是本舰1s内实际走过的路线,S’是补偿算法近似路线,R为转向半径,dx和dy分别为切向和法向的每秒位置误差,每秒位置误差为dd,则有:
经计算,有:R=573.2484m;dx=-0.0046m;dy=0.0873m;dd=0.0874m。
设航迹历史周期为N个,每个扫描周期为t=10s,补偿间隔为1s下不同历史周期的关联点迹对应本舰运动误差Es和对应的补偿误差Ea分别如图2和3所示。
而方法一先求目标的相对航向航速,再叠加本舰航向航速,相当于一个扫描周期做一次补偿,根据上面假定其补偿间隔为10s,则误差计算如下:
dx=0.04565m;dy=0.87089m;dd=0.8721m
方法一中补偿间隔为10s,不同历史周期的关联点迹对应本舰运动误差Es和对应的补偿误差Ea分别如图4和5所示。
根据两种方法的误差数据,本发明选择方法二,实现方法如下,具体流程图如图6。
(1)保存历史关联点迹和历史本舰航向航速信息;
(2)将历史关联点迹时间和当前时间范围内的历史本舰航向航速挑选出来,利用公式 (1)和(2)将每个历史关联点迹由相对位置补偿为绝对位置;
(3)将补偿后的历史关联点迹输入最小二乘法系统;
(4)根据最小二乘法求得拟合参数,绝对航向航速信息。

Claims (3)

1.一种舰载雷达海目标航向航速高精度补偿计算方法,其特征在于:保存本舰航向航速信息和雷达探测海目标若干周期的历史关联点迹,通过本舰航向航速信息把历史关联点迹由相对量补偿为绝对量;采用最小二乘法计算得到目标绝对航向航速。
2.根据权利要求1所述的一种舰载雷达海目标航向航速高精度补偿计算方法,其特征在于所述航向航速的计算方法包括以下步骤:
(1)保存历史关联点迹和历史本舰航向航速信息;
(2)将历史关联点迹时间和当前时间范围内的历史本舰航向航速挑选出来,利用公式(1)和公式(2)将每个历史关联点迹由相对位置补偿为绝对位置,其中公式(1)和公式(2)如下;
(3)将补偿后的历史关联点迹输入最小二乘法系统;
(4)根据最小二乘法求得拟合参数,绝对航向航速信息。
3.根据权利要求1或权利要求2所述的一种舰载雷达海目标航向航速高精度补偿计算方法,其特征在于:补偿原则是历史周期越久的关联点迹,对应的需要补偿的本舰航向航速信息越多。
CN201710621893.5A 2017-07-27 2017-07-27 一种舰载雷达海目标航向航速高精度补偿计算方法 Active CN107561527B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710621893.5A CN107561527B (zh) 2017-07-27 2017-07-27 一种舰载雷达海目标航向航速高精度补偿计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710621893.5A CN107561527B (zh) 2017-07-27 2017-07-27 一种舰载雷达海目标航向航速高精度补偿计算方法

Publications (2)

Publication Number Publication Date
CN107561527A true CN107561527A (zh) 2018-01-09
CN107561527B CN107561527B (zh) 2020-11-10

Family

ID=60974994

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710621893.5A Active CN107561527B (zh) 2017-07-27 2017-07-27 一种舰载雷达海目标航向航速高精度补偿计算方法

Country Status (1)

Country Link
CN (1) CN107561527B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2939904A1 (fr) * 2008-12-16 2010-06-18 Thales Sa Methode de traitement de poursuite multi-capteur a temps de latence reduit
CN102176010A (zh) * 2011-01-21 2011-09-07 西安电子科技大学 基于多发单收的无源雷达定位跟踪系统及定位跟踪方法
CN104280723A (zh) * 2014-10-22 2015-01-14 四川大学 雷达杂波区域新建系统航迹处理方法
CN105353368A (zh) * 2015-11-09 2016-02-24 中国船舶重工集团公司第七二四研究所 一种基于策略判决的自适应变结构雷达对海目标跟踪方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2939904A1 (fr) * 2008-12-16 2010-06-18 Thales Sa Methode de traitement de poursuite multi-capteur a temps de latence reduit
CN102176010A (zh) * 2011-01-21 2011-09-07 西安电子科技大学 基于多发单收的无源雷达定位跟踪系统及定位跟踪方法
CN104280723A (zh) * 2014-10-22 2015-01-14 四川大学 雷达杂波区域新建系统航迹处理方法
CN105353368A (zh) * 2015-11-09 2016-02-24 中国船舶重工集团公司第七二四研究所 一种基于策略判决的自适应变结构雷达对海目标跟踪方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
PHANI CHAVALI ET AL.: "Concurrent Particle Filtering and Data Association Using Game Theory for Tracking Multiple Maneuvering Targets", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 *
于振华等: "一种雷达航迹融合技术的设计原理及实现", 《舰船电子对抗》 *

Also Published As

Publication number Publication date
CN107561527B (zh) 2020-11-10

Similar Documents

Publication Publication Date Title
CN108957466B (zh) 移动机器人的雷达数据补偿方法、装置、设备和存储介质
CN106871924B (zh) 一种基于卫星定位信号的车辆行驶里程计算方法
US11300663B2 (en) Method for predicting a motion of an object
CN108759838A (zh) 基于秩卡尔曼滤波器的移动机器人多传感器信息融合方法
CN103268407B (zh) 一种基于拉格朗日插值及卡尔曼滤波的轨道数据插值方法
US9864067B2 (en) Method for determining a current position of a motor vehicle in a geodetic coordinate system and motor vehicle
CN101511022A (zh) 一种机载视频压缩与目标跟踪联合实现方法
CN104035083A (zh) 一种基于量测转换的雷达目标跟踪方法
CN110727923B (zh) 一种短期潮位预报方法
CN109143192A (zh) 多目标有效起始确认管理方法、系统及计算机可读存储介质
CN110673148A (zh) 一种主动声纳目标实时航迹解算方法
Chen et al. A novel fusion methodology to bridge GPS outages for land vehicle positioning
CN110231620A (zh) 一种噪声相关系统跟踪滤波方法
CN109581302A (zh) 一种车载雷达数据跟踪方法及系统
CN104318586B (zh) 基于自适应形态学滤波的运动模糊方向估计方法及装置
CN106153045B (zh) 一种抑制gnss信息异常的滤波增益动态调整方法
CN114779204A (zh) 一种基于雷达目标幅度最小二乘跟踪处理方法及系统
CN107765240A (zh) 一种运动状态的判断方法、装置和电子设备
CN107561527A (zh) 一种舰载雷达海目标航向航速高精度补偿计算方法
CN104155669B (zh) 基于北斗二代卫星导航系统的接收机深组合实现方法
CN104020468A (zh) 一种基于二次雷达的近距离目标更新周期的计算方法
CN110007298B (zh) 一种目标超前预测跟踪方法
CN109188422A (zh) 一种基于lu分解的卡尔曼滤波目标跟踪方法
CN114488237A (zh) 定位方法、装置、电子设备及智能驾驶方法
CN100399048C (zh) 变换线性调频尺度成像算法因子的可复用计算方法及装置

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