CN108761442B - 一种无源雷达目标定位方法及装置 - Google Patents
一种无源雷达目标定位方法及装置 Download PDFInfo
- Publication number
- CN108761442B CN108761442B CN201810557096.XA CN201810557096A CN108761442B CN 108761442 B CN108761442 B CN 108761442B CN 201810557096 A CN201810557096 A CN 201810557096A CN 108761442 B CN108761442 B CN 108761442B
- Authority
- CN
- China
- Prior art keywords
- radiation source
- external radiation
- error
- target
- time difference
- 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
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000005855 radiation Effects 0.000 claims abstract description 74
- 239000011159 matrix material Substances 0.000 claims description 57
- 239000013598 vector Substances 0.000 claims description 14
- 238000005259 measurement Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/42—Simultaneous measurement of distance and other co-ordinates
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种无源雷达目标定位方法及装置,本发明是在考虑外辐射源位置误差的条件下进行设计的,定位结果为外辐射源存在误差时的最优估计结果;而且可以在不需要目标位置和速度初值先验的情况下实现联合时差和频差的目标无源定位;同时本发明由于是关于目标位置和速度参数代数闭式解的形式,计算复杂度低,不存在收敛性问题。
Description
技术领域
本发明属于无源定位技术领域,具体涉及一种无源雷达目标定位方法及装置。
背景技术
无源雷达本身不辐射源电磁波,而是通过接收和处理目标对环境中现有非合作照射源的反射或散射信号,来探测和定位目标。与传统有源雷达相比,无源雷达系统具有结构简单、成本低、隐蔽性高、不占用频谱资源、抗电子干扰等优点。多年来,一直是国际雷达领域的研究热点。
目前,可用的外辐射源较多,根据其用途,主要可以分为两类:一类为民用辐射源信号,一类为军用辐射源信号。民用辐射源的位置一般可以较为准确的获取,但是受限于信号波形、功率等因素的限制,基于民用辐射源的无源雷达系统的探测距离和定位精度往往较低。而军用辐射源由于信号本身即为目标探测而设计,因此基于军用辐射源的无源雷达系统具有更高的探测距离和定位精度,近年来受到学术界关注。
但由于采用的军用辐射源(例如目前关注较多的预警机信号)是非合作,甚至敌方的,因此其位置往往无法准确获得,仅能通过ESM系统进行估计,得到的外辐射源位置是不准确的,因而常常造成无源雷达系统的定位精度低。
发明内容
本发明的目的在于提供一种无源雷达目标定位方法及装置,用以解决现有技术中无源雷达系统的定位精度不高的问题。
为解决上述技术问题,本发明的技术方案为:
本发明提供了一种无源雷达目标定位方法,包括如下方法方案:
方法方案一,包括如下步骤:
构造无源雷达系统的时差观测方程和频差观测方程;
借助辅助向量以及构造的时差观测方程和频差观测方程,构造联合时差和频差的关于目标位置和速度的第一伪线性方程组;以第一加权矩阵W1为加权矩阵,采用加权最小二乘法对所述第一伪线性方程组进行求解,得到第一最小二乘解θ1;其中,所述第一加权矩阵W1为单位阵;其中,u为目标位置,Rt为目标到外辐射源之间的距离;
根据第一最小二乘解θ1,构造第二加权矩阵W2,以第二加权矩阵W2为加权矩阵,采用加权最小二乘法对所述第一伪线性方程组进行求解,得到第二最小二乘解θ2;
根据辅助向量中的辅助参数Rt、与目标位置参数u、之间的关系,构造第二伪线性方程组;以第三加权矩阵W3为加权矩阵,采用加权最小二乘法对所述第二伪线性方程组求解,得到最终最小二乘解θ3,所述最终最小二乘解为预测的目标位置和速度。
方法方案二,在方法方案一的基础上,所述时差观测方程和频差观测方程分别为:
rn=Rt+Rnr+Δrn
将时差观测方程和频差观测方程表示成矩阵形式其中,r=[r1 r2... rN]T,测量误差为Δr=[Δr1 Δr2 ... ΔrN]T,协方差矩阵为Qα=E[ΔαΔαT];将时差和频差的观测方程中外辐射源位置和速度的真实值分别替换为和并移项、平方、整理,得到联合时差和频差的关于目标位置和速度的第一伪线性方程组;其中,为含有误差的外辐射源位置,st为外辐射源位置,Δst为误差外辐射源位置误差,为含有误差的外辐射源速度,为外辐射源速度,为外辐射源速度误差。
方法方案七,在方法方案六的基础上,所述第三加权矩阵为:
本发明还提供了一种无源雷达目标定位装置,包括如下装置方案:
装置方案一,包括处理器,所述处理器用于执行指令实现如下方法:
构造无源雷达系统的时差观测方程和频差观测方程;
借助辅助向量以及构造的时差观测方程和频差观测方程,构造联合时差和频差的关于目标位置和速度的第一伪线性方程组;以第一加权矩阵W1为加权矩阵,采用加权最小二乘法对所述第一伪线性方程组进行求解,得到第一最小二乘解θ1;其中,所述第一加权矩阵W1为单位阵;其中,u为目标位置,Rt为目标到外辐射源之间的距离;
根据第一最小二乘解θ1,构造第二加权矩阵W2,以第二加权矩阵W2为加权矩阵,采用加权最小二乘法对所述第一伪线性方程组进行求解,得到第二最小二乘解θ2;
根据辅助向量中的辅助参数Rt、与目标位置参数u、之间的关系,构造第二伪线性方程组;以第三加权矩阵W3为加权矩阵,采用加权最小二乘法对所述第二伪线性方程组求解,得到最终最小二乘解θ3,所述最终最小二乘解为预测的目标位置和速度。
装置方案二,在装置方案一的基础上,所述时差观测方程和频差观测方程分别为:
rn=Rt+Rnr+Δrn
将时差观测方程和频差观测方程表示成矩阵形式其中,r=[r1 r2... rN]T,测量误差为Δr=[Δr1 Δr2 ... ΔrN]T,协方差矩阵为Qα=E[ΔαΔαT];将时差和频差的观测方程中外辐射源位置和速度的真实值分别替换为和并移项、平方、整理,得到联合时差和频差的关于目标位置和速度的第一伪线性方程组;其中,为含有误差的外辐射源位置,st为外辐射源位置,Δst为误差外辐射源位置误差,为含有误差的外辐射源速度,为外辐射源速度,为外辐射源速度误差。
装置方案七,在装置方案六的基础上,所述第三加权矩阵为:
本发明的有益效果:
本发明的无源雷达目标定位方法及装置,在考虑外辐射源位置误差的条件下,有针对性的对无源雷达目标进行定位,从而解决了因外辐射源位置误差的存在导致无源雷达系统的定位结果不准确的问题。该定位结果为外辐射源存在误差时的最优估计结果;而且,可以在不需要目标位置和速度初值先验的情况下实现联合时差和频差的目标无源定位。同时,本发明由于是关于目标位置和速度参数代数闭式解的形式,计算复杂度低,不存在收敛性问题。
附图说明
图1是本发明的方法流程图;
图2是实验仿真中外辐射源和接收站几何位置示意图;
图3是目标位置估计误差随测量误差变化的仿真对比图;
图4是目标速度估计误差随测量误差变化的仿真对比图;
图5是目标位置估计误差随外辐射源位置误差变化的仿真对比图;
图6是目标速度估计误差随外辐射源位置误差变化的仿真对比图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚,下面结合附图及实施例,对本发明作进一步的详细说明。
本发明提供了一种无源雷达目标定位装置,该装置包括处理器,所述处理器用于执行指令实现本发明的无源雷达目标定位方法。下面对该方法做详细、具体的介绍。
首先,构造无源雷达系统的时差观测方程和频差观测方程。
目标到外辐射源之间的距离为Rt=||u-st||,目标到接收站n之间的距离为Rnr=||u-snr||。那么,考虑实际测量中不可避免地存在误差,得到无源雷达系统的双基地(对应于时差)及其变化率(对应于频差)的观测方程为:
然后,将时差观测方程和频差观测方程,即式(1)表示成矩阵形式:
G1θ′=h1+Δh1 (3)
因此,θ1的加权最小二乘解为:
在第一次求解时,以单位阵W1=IN×N(即第一加权矩阵)为加权矩阵W,对方程组(3)采用加权最小二乘法进行求解,求解得到第一加权最小二乘解:
然后,利用第一最小二乘解θ1中目标位置和速度的估计值,构造一个新的加权矩阵,即第二加权矩阵:
以该第二加权W2为加权矩阵,重新对方程组(3)采用加权最小二乘法进行求解,求解得到第二加权最小二乘解:
该步骤中,重复迭代1~2次即可,尽可能的提高目标位置和速度的估计精度。
G2θ″=h2+Δh2 (8)
构造第三加权矩阵W3,以第三加权矩阵为加权矩阵,对方程组(8)进行求解,得到最终最小二乘解:
其中,第三加权矩阵W3为:
最终最小二乘解即为预测的目标位置和速度。
现举一个具体的实例来对本发明的方法的效果做进一步的说明。
图2无源雷达系统和目标的几何位置示意图,对本发明进行模拟实验仿真。图3和图4分别展示了在外辐射源位置误差一定的情况下,本发明目标位置和速度估计误差随测量误差变化的仿真对比,可以看出本发明考虑外辐射源位置误差后系统对目标位置和速度的估计精度显著提高;图5和图6分别展示了在测量误差一定的情况下,本发明目标位置和速度的估计误差随外辐射源位置误差变化的仿真对比,结果表明本发明可在外辐射源位置误差增大时仍能保持较高的定位精度,估计性能明显优于忽略外辐射源位置误差的定位算法。
Claims (6)
1.一种无源雷达目标定位方法,其特征在于,包括如下步骤:
构造无源雷达系统的时差观测方程和频差观测方程;
借助辅助向量以及构造的时差观测方程和频差观测方程,构造联合时差和频差的关于目标位置和速度的第一伪线性方程组;以第一加权矩阵W1为加权矩阵,采用加权最小二乘法对所述第一伪线性方程组进行求解,得到第一最小二乘解θ1;其中,所述第一加权矩阵为单位阵;其中,u为目标位置,Rt为目标到外辐射源之间的距离;
根据第一最小二乘解θ1,构造第二加权矩阵W2,以第二加权矩阵W2为加权矩阵,采用加权最小二乘法对所述第一伪线性方程组进行求解,得到第二最小二乘解θ2;
根据辅助向量中的辅助参数Rt、与目标位置参数u、之间的关系,构造第二伪线性方程组;以第三加权矩阵W3为加权矩阵,采用加权最小二乘法对所述第二伪线性方程组求解,得到最终最小二乘解θ3,所述最终最小二乘解为预测的目标位置和速度;
所述时差观测方程和频差观测方程分别为:
rn=Rt+Rnr+Δrn
其中,Rt=||u-st||,Rnr=||u-snr||,u=[x y z]T,st为外辐射源位置,st=[xt yt zt]T,snr为接收站的位置,snr=[xnr ynr znr]T,为外辐射源速度, 为接收站速度,Rnr为目标到接收站之间的距离;
所述第一伪线性方程组为:G1θ′=h1+Δh1;
所述第二伪线性方程组为:G2θ″=h2+Δh2;
4.一种无源雷达目标定位装置,其特征在于,包括处理器,所述处理器用于执行指令实现如下方法:
构造无源雷达系统的时差观测方程和频差观测方程;
借助辅助向量以及构造的时差观测方程和频差观测方程,构造联合时差和频差的关于目标位置和速度的第一伪线性方程组;以第一加权矩阵W1为加权矩阵,采用加权最小二乘法对所述第一伪线性方程组进行求解,得到第一最小二乘解θ1;其中,所述第一加权矩阵W1为单位阵;其中,u为目标位置,Rt为目标到外辐射源之间的距离;
根据第一最小二乘解θ1,构造第二加权矩阵W2,以第二加权矩阵W2为加权矩阵,采用加权最小二乘法对所述第一伪线性方程组进行求解,得到第二最小二乘解θ2;
根据辅助向量中的辅助参数Rt、与目标位置参数u、之间的关系,构造第二伪线性方程组;以第三加权矩阵W3为加权矩阵,采用加权最小二乘法对所述第二伪线性方程组求解,得到最终最小二乘解θ3,所述最终最小二乘解为预测的目标位置和速度;
所述时差观测方程和频差观测方程分别为:
rn=Rt+Rnr+Δrn
其中,Rt=||u-st||,Rnr=||u-snr||,u=[x y z]T,st为外辐射源位置,st=[xt yt zt]T,snr为接收站的位置,snr=[xnr ynr znr]T,为外辐射源速度, 为接收站速度,Rnr为目标到接收站之间的距离;
所述第一伪线性方程组为:G1θ′=h1+Δh1;
所述第二伪线性方程组为:G2θ″=h2+Δh2;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810557096.XA CN108761442B (zh) | 2018-06-01 | 2018-06-01 | 一种无源雷达目标定位方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810557096.XA CN108761442B (zh) | 2018-06-01 | 2018-06-01 | 一种无源雷达目标定位方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108761442A CN108761442A (zh) | 2018-11-06 |
CN108761442B true CN108761442B (zh) | 2021-06-11 |
Family
ID=64001911
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810557096.XA Active CN108761442B (zh) | 2018-06-01 | 2018-06-01 | 一种无源雷达目标定位方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108761442B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112763980B (zh) * | 2020-12-28 | 2022-08-05 | 哈尔滨工程大学 | 一种基于方位角及其变化率的目标运动分析方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8866672B2 (en) * | 2012-04-05 | 2014-10-21 | L-3 Communications Integrated Systems Lp | Cooperative systems and methods for TDOA-based emitter location |
CN106353720A (zh) * | 2016-09-04 | 2017-01-25 | 中国人民解放军海军航空工程学院 | 基于tdoa/groa的多站连续定位模型 |
CN107526073A (zh) * | 2017-08-22 | 2017-12-29 | 哈尔滨工程大学 | 一种运动多站无源时差频差联合定位方法 |
-
2018
- 2018-06-01 CN CN201810557096.XA patent/CN108761442B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8866672B2 (en) * | 2012-04-05 | 2014-10-21 | L-3 Communications Integrated Systems Lp | Cooperative systems and methods for TDOA-based emitter location |
CN106353720A (zh) * | 2016-09-04 | 2017-01-25 | 中国人民解放军海军航空工程学院 | 基于tdoa/groa的多站连续定位模型 |
CN107526073A (zh) * | 2017-08-22 | 2017-12-29 | 哈尔滨工程大学 | 一种运动多站无源时差频差联合定位方法 |
Non-Patent Citations (3)
Title |
---|
Closed-Form Localization for Distributed MIMO Radar Systems Using Time Delay Measurements;Chee-Hyun Park等;《IEEE Transactions on Wireless Communications》;20160229;第15卷(第2期);第1481-1484页 * |
使用单校正源的多站时频差定位算法;吴魏等;《信息工程大学学报》;20160229;第17卷(第1期);第7-9页 * |
基于TDOA和TOA的定位技术研究;赵海霞;《中国优秀硕士学位论文全文数据库 信息科技辑》;20160315(第03期);第44-49页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108761442A (zh) | 2018-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104076348B (zh) | 一种雷达超视距基线无源协同定位方法 | |
CN106324597B (zh) | 基于pfa的大转角isar雷达的平动补偿和成像方法 | |
CN107831480B (zh) | 弹载雷达和差通道稳健自适应杂波抑制方法 | |
CN104749570B (zh) | 一种移不变机载双基合成孔径雷达目标定位方法 | |
CN104730520B (zh) | 基于子孔径合成的圆周sar后向投影自聚焦方法 | |
CN102411136A (zh) | 一种扩展基线解模糊的相位干涉仪测向方法 | |
CN104101876B (zh) | 外辐射源雷达中一种基于随机有限集的多目标跟踪方法 | |
CN106405533A (zh) | 基于约束加权最小二乘的雷达目标联合同步与定位方法 | |
CN103852759A (zh) | 扫描雷达超分辨成像方法 | |
CN106918810B (zh) | 一种存在阵元幅相误差时的微波关联成像方法 | |
CN109031190B (zh) | 一种高重频脉冲信号无源时差定位方法 | |
CN105005036A (zh) | 一种用于近程mimo成像的传播损耗补偿方法 | |
CN104077498A (zh) | 一种结合目标角度的外辐射源雷达多目标跟踪方法 | |
CN105445718A (zh) | 一种基于阵列重构的分布式多载舰超视距雷达的doa估计方法 | |
CN104678370B (zh) | 用于估计和补偿极化校准二面角反射器双站散射影响的方法 | |
CN105738887A (zh) | 基于多普勒通道划分的机载雷达杂波功率谱的优化方法 | |
CN106932759A (zh) | 一种用于主被动雷达的协同定位方法 | |
CN104931923A (zh) | Grid Iterative ESPRIT,一种可扩展的用于均匀圆阵二维到达角的快速估计算法 | |
CN110471029B (zh) | 一种基于扩展卡尔曼滤波的单站无源定位方法及装置 | |
CN108761442B (zh) | 一种无源雷达目标定位方法及装置 | |
CN109521418B (zh) | 基于干涉场的地基雷达测角方法 | |
CN104931965B (zh) | 一种固定发射站双基正侧视sar成像方法 | |
CN106291547A (zh) | 基于天线方向图辅助的多普勒模糊分量自适应抑制方法 | |
CN109633596A (zh) | 脉冲多普勒雷达体制下的多目标运动参数联合估计方法 | |
CN1601298A (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 |