CN107561527B - 一种舰载雷达海目标航向航速高精度补偿计算方法 - Google Patents
一种舰载雷达海目标航向航速高精度补偿计算方法 Download PDFInfo
- Publication number
- CN107561527B CN107561527B CN201710621893.5A CN201710621893A CN107561527B CN 107561527 B CN107561527 B CN 107561527B CN 201710621893 A CN201710621893 A CN 201710621893A CN 107561527 B CN107561527 B CN 107561527B
- Authority
- CN
- China
- Prior art keywords
- ship
- target
- course
- historical
- course speed
- 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.)
- Active
Links
Images
Landscapes
- Navigation (AREA)
- Radar Systems Or Details Thereof (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)
由于本发明实施例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 (2)
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 CN107561527A (zh) | 2018-01-09 |
CN107561527B true 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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102176010A (zh) * | 2011-01-21 | 2011-09-07 | 西安电子科技大学 | 基于多发单收的无源雷达定位跟踪系统及定位跟踪方法 |
CN104280723A (zh) * | 2014-10-22 | 2015-01-14 | 四川大学 | 雷达杂波区域新建系统航迹处理方法 |
FR2939904B1 (fr) * | 2008-12-16 | 2015-08-07 | Thales Sa | Methode de traitement de poursuite multi-capteur a temps de latence reduit |
CN105353368A (zh) * | 2015-11-09 | 2016-02-24 | 中国船舶重工集团公司第七二四研究所 | 一种基于策略判决的自适应变结构雷达对海目标跟踪方法 |
-
2017
- 2017-07-27 CN CN201710621893.5A patent/CN107561527B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2939904B1 (fr) * | 2008-12-16 | 2015-08-07 | 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)
Title |
---|
Concurrent Particle Filtering and Data Association Using Game Theory for Tracking Multiple Maneuvering Targets;Phani Chavali et al.;《IEEE TRANSACTIONS ON SIGNAL PROCESSING》;20131015;第61卷(第20期);第4943-4948页 * |
一种雷达航迹融合技术的设计原理及实现;于振华等;《舰船电子对抗》;20150831;第38卷(第4期);第36-39页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107561527A (zh) | 2018-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108362281B (zh) | 一种长基线水下潜艇匹配导航方法及系统 | |
CN106950562B (zh) | 一种基于预测值量测转换的状态融合目标跟踪方法 | |
CN109059964B (zh) | 一种基于“重力峰”的惯性导航与重力测量双校准方法 | |
CN109143224B (zh) | 一种多目标关联方法和装置 | |
CN103791902A (zh) | 适用于高机动载体的星敏感器自主导航方法 | |
CN103047982A (zh) | 基于角度信息的目标自适应跟踪方法 | |
CN110672871B (zh) | 一种基于组合惯导信息和光电转塔视频跟踪的运动目标测速测向方法 | |
CN109738902B (zh) | 一种基于同步信标模式的水下高速目标高精度自主声学导航方法 | |
CN106767760A (zh) | 基于多维度的多源船舶目标融合方法 | |
CN110673148A (zh) | 一种主动声纳目标实时航迹解算方法 | |
CN110702110A (zh) | 一种基于无迹卡尔曼滤波的舰船升沉运动测量方法 | |
CN111307160B (zh) | 一种自主水下航行器长航时航迹的修正方法及装置 | |
CN110738275A (zh) | 基于ut-phd的多传感器序贯融合跟踪方法 | |
CN109581302A (zh) | 一种车载雷达数据跟踪方法及系统 | |
JP2012154647A (ja) | 物標運動推定装置 | |
CN107561527B (zh) | 一种舰载雷达海目标航向航速高精度补偿计算方法 | |
CN113049052B (zh) | 一种大含沙量河流的流量测量方法 | |
Xu et al. | Optimal design of damping network based on DVL velocity and IMU | |
Yuan et al. | Reaearch on underwater integrated navigation system based on SINS/DVL/magnetometer/depth-sensor | |
CN111008549A (zh) | 基于样本熵和ifoa-grnn的uuv平台dvl信号失真重构方法 | |
Kim et al. | Optimal guidance for autonomous underwater vehicle navigation within undersea areas of current disturbances | |
CN113032994B (zh) | 断面多点潮流调和常数提取计算方法、装置及电子设备 | |
Bo et al. | GPS/INS/speed log integrated navigation system based on MAKF and priori velocity information | |
Park et al. | An effective target selection algorithm for ASM (Anti-ship Missile) | |
CN110057383B (zh) | 一种auv推位导航系统杆臂误差标校方法 |
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 |