CN110133609B - 一种机载外辐射源状态未知下移动接收源路径优化方法 - Google Patents
一种机载外辐射源状态未知下移动接收源路径优化方法 Download PDFInfo
- Publication number
- CN110133609B CN110133609B CN201910271672.9A CN201910271672A CN110133609B CN 110133609 B CN110133609 B CN 110133609B CN 201910271672 A CN201910271672 A CN 201910271672A CN 110133609 B CN110133609 B CN 110133609B
- Authority
- CN
- China
- Prior art keywords
- external radiation
- radiation source
- target
- state
- source
- 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
- 230000005855 radiation Effects 0.000 title claims abstract description 117
- 238000000034 method Methods 0.000 title claims abstract description 14
- 238000005259 measurement Methods 0.000 claims abstract description 19
- 238000005457 optimization Methods 0.000 claims abstract description 16
- 238000001914 filtration Methods 0.000 claims abstract description 6
- 230000008859 change Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 4
- 230000004888 barrier function Effects 0.000 claims description 2
- 238000011478 gradient descent method Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 claims description 2
- 230000007704 transition Effects 0.000 claims description 2
- 238000013178 mathematical model Methods 0.000 abstract 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000007499 fusion processing Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details 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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种机载外辐射源状态未知下移动接收源路径优化方法。该方法针对机载外辐射源双基雷达系统中外辐射源状态未知且运动情况,提出了一种基于序贯信息修正的扩展卡尔曼滤波跟踪算法联合估计目标的状态和外辐射源的状态。基于目标位置估计误差PCRLB指标建立接收源路径优化的数学模型,通过实时优化接收源位置,获得更精确的目标量测信息,进一步提高机载外辐射源雷达定位跟踪性能。
Description
技术领域
本发明属于传感器协同定位技术领域,具体涉及一种机载外辐射源状态未知下移动接收源路径优化方法。
背景技术
外辐射源雷达利用第三方非合作信号源作为目标的机会照射源,通过接收信号源的直达波和经过目标辐射的回波处理,实现对“无线电静默”目标定位跟踪,具有低成本、隐蔽性好、抗干扰能力强等特点。传统的外辐射源雷达系统中接收源和辐射源大多安装在地基平台,其位置固定不变的。机载外辐射源雷达将外辐射源从固定的地基平台扩展到可运动的机载平台,大幅扩展了空基预警系统的有效覆盖范围。
机载外辐射源雷达系统作为一种双/多基地结构的传感器组网系统,通过数据融合处理,实现对目标的定位跟踪。由于外辐射源大多是非合作的,外辐射源自身状态的准确性和可靠性直接影响目标的定位精度,需要对目标和外辐射源的状态进行联合估计。此外,目标定位精度受“外辐射源-目标-接收源”三者之间的几何关系影响,可以通过目标定位误差的后验克劳美罗下界(PCRLB)或目标定位精度几何解释(GDOP)指标描述。接收源通过运动规划,优化三者之间的几何关系可以获得目标位置更加精确的量测信息,从而减少定位误差。根据外辐射源和目标的联合位置估计值来实时优化接收源的位置是进一步提高机载外辐射源雷达定位跟踪性能的关键。
发明内容
本发明针对机载外辐射源双基雷达系统中外辐射源状态未知情况,设计了一种外辐射源和目标状态联合估计算法,推导了联合估计下的跟踪性能下界,并基于此设计了接收源路径优化算法。通过联合估计和路径优化提高外辐射源雷达系统目标跟踪性能。
本发明以机载外辐射源双基雷达系统为对象,通过目标/外辐射源状态联合估计和接收源路径优化提高机载外辐射源雷达系统的定位跟踪性能。本发明方法的具体步骤是:
步骤1:建立目标和外辐射源的运动模型,预测k+1时刻目标的状态和外辐射源的状态。
步骤2:计算目标和外辐射源状态联合估计误差的PCRLB指标。
步骤3:以目标和外辐射源位置误差的PCRLB为优化目标,考虑接收源运动约束,建立接收源路径优化模型。
步骤5:计算接收源在最优位置下的目标方位角、外辐射源方位角、双基距和双基距离率信息的量测信息。
步骤6:采用序贯信息修正的扩展卡尔曼滤波跟踪算法联合估计目标的状态和外辐射源的状态。
本发明的有益效果:
1.考虑机载外辐射源雷达系统中外辐射源状态未知情况,通过外辐射源和目标状态联合估计提高目标定位精度。
2.在双基外辐射源雷达系统角度/双基距/双基距率定位体制下,针对运动目标,给出了机载接收源的最优航迹规划。通过实时优化收源位置,获得更精确的目标量测信息,进一步提高机载外辐射源雷达定位跟踪性能。
具体实施方式:
步骤1:建立目标和外辐射源的运动模型,预测k+1时刻目标的状态和外辐射源的状态。
X(k+1)=F·X(k)+v(k) (1)
Xt(k+1)=Ft·Xt(k)+vt(k) (2)
其中,F和Ft分别为目标和外辐射源的状态转移矩阵;目标运动噪声为v(k),发射源运动噪声为vt(k),假设v(k)、vt(k)为零均值协方差分别为Q和Qt的高斯白噪声。
步骤1.2:基于k时刻目标的估计状态和外辐射源的估计状态,预测k+1时刻目标的状态和外辐射源的状态
步骤2:计算目标和外辐射源状态联合估计误差的后验克劳美罗下界PCRLB指标。
对目标和外辐射源的状态进行扩维,令Xa(k)=[X(k),Xt(k)]T,建立运动方程
Xa(k+1)=Fa·Xa(k)+va(k) (3)
其中,Fa=diag(F,Ft),va(k)=[v(k),vt(k)]T为协方差为Qa=diag(Q,Qt)的高斯白噪声。
计算目标和外辐射源状态联合估计误差的Fisher信息矩阵FIM
J(k+1|k)=(Qa+FaJ(k|k)-1Fa T)-1+HaR-1(Ha)T (4)
则PCRLB指标为
PCRLB(k+1|k)=J(k+1|k)-1 (5)
得到目标和外辐射源位置误差的PCRLB如下
步骤3:以目标和外辐射源位置误差的PCRLB为优化目标,考虑接收源运动约束,建立接收源路径优化模型如下
步骤4.1:获得接收源k时刻位置和速度,令i=1,设定初始障碍因子r1>0和允许误差ξ1>0。
步骤4.2:构造障碍函数B。
步骤5:计算接收源在最优位置下的目标方位角、外辐射源方位角、双基距和双基距离率信息的量测信息Z(k+1)。
其中,x(k+1)和y(k+1)表示目标在k+1时刻的真实位置,xt(k+1)和yt(k+1)表示外辐射源在k+1时刻的真实位置;rTO(k+1)为目标到外辐射的距离,rOR(k+1)为目标到接收源的距离,rTR(k+1)为外辐射源到接收源的距离;为目标到外辐射的距离变化率,为目标到接收源的距离变化率,为外辐射源到接收源的距离变化率;nθ(k+1)、nd(k+1)和分别表示目标方位角量测误差、外辐射源方位角量测误差、双基距量测误差和双基距离率量测误差。
步骤6.1.1:对目标状态估计误差协方差进行一步预测
P(k+1|k)=FP(k|k)FT+Q (12)
其中,P(k|k)为k时刻目标状态估计误差协方差,P(k+1|k)为k时刻目标状态估计误差协方差的一步预测。
对外辐射源的状态估计误差协方差进行一步预测
Pt(k+1|k)=FPt(k|k)FT+Qt (13)
其中,Pt(k|k)为k时刻外辐射源状态估计误差协方差,Pt(k+1|k)为k时刻外辐射源状态估计误差协方差的一步预测。
步骤6.1.2:引入外辐射源的状态信息,更新目标的新息协方差S
步骤6.1.3:更新目标的滤波增益K
其中,Zp(k+1)为基于目标和外辐射源一步预测值得到的量测的预测值。
步骤6.1.5:更新目标的协方差P(k+1|k+1)
步骤6.2.2:引入目标的状态信息,更新外辐射源的新息协方差St
步骤6.2.3:更新外辐射源的滤波增益Kt
步骤6.2.5:更新外辐射源的协方差
Claims (1)
1.一种机载外辐射源状态未知下移动接收源路径优化方法,其特征在于:该方法包括以下步骤:
步骤1:建立目标和外辐射源的运动模型,预测k+1时刻目标的状态和外辐射源的状态;
X(k+1)=F·X(k)+v(k) (1)
Xt(k+1)=Ft·Xt(k)+vt(k) (2)
其中,F和Ft分别为目标和外辐射源的状态转移矩阵;目标运动噪声为v(k),外辐射源运动噪声为vt(k),假设v(k)、vt(k)为零均值协方差分别为Q和Qt的高斯白噪声;
步骤1.2:基于k时刻目标的估计状态和外辐射源的估计状态,预测k+1时刻目标的状态和外辐射源的状态
步骤2:计算目标和外辐射源状态联合估计误差的后验克劳美罗下界PCRLB指标;
对目标和外辐射源的状态进行扩维,令Xa(k)=[X(k),Xt(k)]T,建立运动方程
Xa(k+1)=Fa·Xa(k)+va(k) (3)
其中,Fa=diag(F,Ft),va(k)=[v(k),vt(k)]T为协方差为Qa=diag(Q,Qt)的高斯白噪声;
计算目标和外辐射源状态联合估计误差的Fisher信息矩阵FIM
J(k+1|k)=(Qa+FaJ(k|k)-1Fa T)-1+HaR-1(Ha)T (4)
则PCRLB指标为
PCRLB(k+1|k)=J(k+1|k)-1 (5)
得到目标和外辐射源位置误差的PCRLB如下
步骤3:以目标和外辐射源位置误差的PCRLB为优化目标,考虑接收源运动约束,建立接收源路径优化模型如下
步骤4.1:获得接收源k时刻位置和速度,令i=1,设定初始障碍因子r1>0和允许误差ξ1>0;
步骤4.2:构造障碍函数B
步骤5:计算接收源在最优位置下的目标方位角、外辐射源方位角、双基距和双基距离率信息的量测信息Z(k+1);
其中,x(k+1)和y(k+1)表示目标在k+1时刻的真实位置,xt(k+1)和yt(k+1)表示外辐射源在k+1时刻的真实位置;rTO(k+1)为目标到外辐射源的距离,rOR(k+1)为目标到接收源的距离,rTR(k+1)为外辐射源到接收源的距离;为目标到外辐射源的距离变化率,为目标到接收源的距离变化率,为外辐射源到接收源的距离变化率;nθ(k+1)、nd(k+1)和分别表示目标方位角量测误差、外辐射源方位角量测误差、双基距量测误差和双基距离率量测误差;
步骤6.1.1:对目标状态估计误差协方差进行一步预测
P(k+1|k)=FP(k|k)FT+Q (12)
其中,P(k|k)为k时刻目标状态估计误差协方差,P(k+1|k)为k时刻目标状态估计误差协方差的一步预测;
对外辐射源的状态估计误差协方差进行一步预测
Pt(k+1|k)=FPt(k|k)FT+Qt (13)
其中,Pt(k|k)为k时刻外辐射源状态估计误差协方差,Pt(k+1|k)为k时刻外辐射源状态估计误差协方差的一步预测;
步骤6.1.2:引入外辐射源的状态信息,更新目标的新息协方差S
步骤6.1.3:更新目标的滤波增益K
其中,Zp(k+1)为基于目标和外辐射源一步预测值得到的量测的预测值;
步骤6.1.5:更新目标的协方差P(k+1|k+1)
步骤6.2.2:引入目标的状态信息,更新外辐射源的新息协方差St
步骤6.2.3:更新外辐射源的滤波增益Kt
步骤6.2.5:更新外辐射源的协方差
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910271672.9A CN110133609B (zh) | 2019-04-04 | 2019-04-04 | 一种机载外辐射源状态未知下移动接收源路径优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910271672.9A CN110133609B (zh) | 2019-04-04 | 2019-04-04 | 一种机载外辐射源状态未知下移动接收源路径优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110133609A CN110133609A (zh) | 2019-08-16 |
CN110133609B true CN110133609B (zh) | 2020-10-16 |
Family
ID=67569428
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910271672.9A Active CN110133609B (zh) | 2019-04-04 | 2019-04-04 | 一种机载外辐射源状态未知下移动接收源路径优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110133609B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112147587B (zh) * | 2020-09-28 | 2022-02-25 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种雷达波束方位中心海上标定方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5911331B2 (ja) * | 2012-02-24 | 2016-04-27 | 三菱電機株式会社 | モノパルス測角装置およびモノパルス測角方法 |
CN103412286B (zh) * | 2013-09-01 | 2015-04-08 | 西安电子科技大学 | 基于mimo雷达的发射极化优化doa估计方法 |
CN106446422B (zh) * | 2016-09-28 | 2019-11-19 | 西京学院 | 一种基于对数似然估计的无源定位跟踪新方法 |
CN106932759B (zh) * | 2017-01-17 | 2020-08-14 | 电子科技大学 | 一种用于主被动雷达的协同定位方法 |
CN107144827B (zh) * | 2017-04-12 | 2019-12-31 | 北京理工大学 | 一种基于联合估计crlb的分布式雷达最优构型构建方法 |
CN108490465B (zh) * | 2018-03-16 | 2020-06-16 | 中国电子科技集团公司第三十六研究所 | 基于时频差与测向的地面同频多运动辐射源跟踪方法及系统 |
-
2019
- 2019-04-04 CN CN201910271672.9A patent/CN110133609B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110133609A (zh) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109916410B (zh) | 一种基于改进平方根无迹卡尔曼滤波的室内定位方法 | |
CN109100714B (zh) | 一种基于极坐标系的低慢小目标跟踪方法 | |
CN109633592B (zh) | 运动观测站误差下外辐射源雷达时差与频差协同定位方法 | |
CN109633589A (zh) | 目标跟踪中基于多模型优化多假设的多目标数据关联方法 | |
CN113064155B (zh) | 一种空中雷达多目标跟踪下航迹关联的优化方法 | |
CN109782289B (zh) | 一种基于基线几何结构约束的水下航行器定位方法 | |
CN104077498B (zh) | 一种结合目标角度的外辐射源雷达多目标跟踪方法 | |
CN104849702B (zh) | 利用ads‑b数据的gm‑ephd滤波雷达系统误差联合估计方法 | |
CN111954153B (zh) | 一种基于uwb定位的智能召车导航定位方法 | |
CN108717174B (zh) | 基于信息论的预测快速协方差交互融合无源协同定位方法 | |
CN110749891B (zh) | 一种可估计未知有效声速的自适应水下单信标定位方法 | |
CN101627319A (zh) | 用于估计目标运动的方法和仪器 | |
CN104101876B (zh) | 外辐射源雷达中一种基于随机有限集的多目标跟踪方法 | |
CN112034445B (zh) | 基于毫米波雷达的车辆运动轨迹跟踪方法和系统 | |
CN113899369A (zh) | 一种基于自适应降噪算法的超宽带/pdr室内定位方法 | |
CN114019498A (zh) | 一种基于分数阶卡尔曼滤波的毫米波雷达单目标跟踪方法 | |
CN110133609B (zh) | 一种机载外辐射源状态未知下移动接收源路径优化方法 | |
CN110596691B (zh) | 一种考虑地球曲率的双基站三维无源定位方法 | |
CN113280821B (zh) | 基于斜率约束和回溯搜索的水下多目标跟踪方法 | |
CN117570980A (zh) | 一种基于uwb与gps融合定位算法的方法及系统 | |
CN117146830A (zh) | 一种自适应多信标航位推算和长基线的紧组合导航方法 | |
CN115421153B (zh) | 一种基于扩展卡尔曼滤波的激光雷达和uwb组合定位方法及系统 | |
CN110133608B (zh) | 合作外辐射源双基雷达协同路径规划多步优化方法 | |
CN114548159B (zh) | 信号干扰下的超宽带精确定位方法 | |
CN116679296A (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 |