CN116148829A - 一种基于非合作外辐射源的多基雷达定位方法 - Google Patents
一种基于非合作外辐射源的多基雷达定位方法 Download PDFInfo
- Publication number
- CN116148829A CN116148829A CN202310214713.7A CN202310214713A CN116148829A CN 116148829 A CN116148829 A CN 116148829A CN 202310214713 A CN202310214713 A CN 202310214713A CN 116148829 A CN116148829 A CN 116148829A
- Authority
- CN
- China
- Prior art keywords
- external radiation
- radiation source
- estimation
- radar
- matrix
- 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.)
- Pending
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 75
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000005259 measurement Methods 0.000 claims description 31
- 239000011159 matrix material Substances 0.000 claims description 29
- 238000004422 calculation algorithm Methods 0.000 claims description 12
- 239000013598 vector Substances 0.000 claims description 12
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims description 2
- 230000004807 localization Effects 0.000 claims 5
- 101100001674 Emericella variicolor andI gene Proteins 0.000 claims 2
- 238000001514 detection method Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012876 topography 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/003—Bistatic radar systems; Multistatic radar systems
-
- 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/006—Theoretical aspects
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
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
技术领域
本发明属于无线电定向技术领域,涉及无源雷达的数据处理方法,具体涉及一种基于非合作外辐射源的多基雷达定位方法。
背景技术
多基外辐射源雷达是指不主动发射电磁信号,通过外部已经存在的第三方非合作信号源,例如手机通讯信号、电视广播信号或敌方雷达信息等,对检测目标的反射回波进行探测、定位和跟踪的雷达装备。由于其本身不辐射电磁信号,可以实现静默探测,具有作用距离远、被动探测、隐蔽性强,生存能力强等优点,近三十年来受到国内外学者的广泛关注。多基外辐射源雷达定位系统从目标散射的外辐射源信号和直达波信号中提取到达时间延迟,并基于该参数估计目标的位置,完成探测、定位和跟踪的任务。
目前,多基外辐射源的目标定位方法主要应用在一些特殊的简化场景,在大部分研究中假设外辐射源位置固定且发射信号时间同步。然而在实际的非合作外辐射源定位环境中,量测信息存在极大的不确定性,主要的误差产生原因为以下三种:
①外辐射源为第三方信号源,且与雷达接收站非合作,外辐射源的位置未知。
②雷达接收站无法知晓外辐射源的信号发射时间,通过信号中的已知结构获得的时间延迟受未知时间偏移的影响。
③当雷达接收站部署在运动平台,例如飞机或海面上时,其位置受到风力或海浪的影响,雷达接收站的位置不能准确获知。
由于这三种不确定情况的存在,严重影响了多基外辐射源雷达定位系统在复杂空间地貌环境中对目标进行定位的性能,限制了非合作外辐射源多基雷达定位的发展。
发明内容
针对现有技术的不足,本发明提出了一种基于非合作外辐射源的多基雷达定位方法,同时考虑了外辐射位置未知、时钟偏差和雷达接收站站址误差三类不确定性因素,提出一个新的联合定位模型和估计算法,实现误差校正和目标精确定位。
一种基于非合作外辐射源的多基雷达定位方法,具体包括以下步骤:
步骤1、雷达接收站接收来自第三方外辐射源的直达波信号与目标散射后的第三方外辐射源信号,获得雷达接收站与第三方外辐射源之间的直接距离d和间接距离r,建立非线性量测方程。
步骤2、考虑量测误差、雷达接收站站址误差和时钟偏差,引入与目标位置、外辐射源位置和时钟偏差相关的辅助变量a0(1)、a0(2)、a0(3),将步骤1建立的非线性量测方程转化为伪线性方程,建立定位模型。
本发明具有以下有益效果:
1、同时考虑时钟误差、接收站站址误差以及外辐射源未知位置对定位精度的影响,在量测方程高度非线性的情况下通过新的闭式解算法,联合估计目标位置、外辐射源位置和时钟偏差,提高了雷达定位系统在不确定、多源误差下处理目标定位问题的能力,取得了较高的定位精度。
2、针对强非线性量测方程,通过引入辅助变量,建立目标位置、外辐射源位置和时钟偏差之间的代数方程,将非线性量测模型转化为伪线性估计模型,在保证估计性能的前提下降低了非线性估计的复杂度。
3、考虑辅助变量与目标位置、外辐射源位置和时钟偏差的关联性,设计两步加权最小二乘算法,进一步提高目标定位估计精度,将多基外辐射源雷达定位系统的应用拓展到更为复杂的空间地貌环境中,具有实际意义。
附图说明
图1基于非合作外辐射源的多基雷达系统示意图。
图2为实施例中仿真实验得到目标辐射源位置参数估计值。
具体实施方式
以下结合附图对本发明作进一步的解释说明。
一种基于非合作外辐射源的多基雷达定位方法,该方法具体包括以下步骤:
步骤1、如图1所示,在K维空间的多基外辐射源雷达定位系统中,包括第三方外辐射源t、M个接收站s和固定目标u。定义第i个雷达接收站的真实位置为i=1,2,…,M。第三方外辐射源的真实位置为to=[xt,yt]T,固定目标的真实位置为uo=[x,y]T,上标T表示矩阵转置。
在非合作的情况下,第三方外辐射源发射信号的时间也是未知的,因此量测方程具有未知的时间偏移量是一个常量,且对于系统中的所有雷达接收站都相同。考虑量测噪声和接收站站址误差,将M个雷达接收站的量测方程整合成矩阵形式,建立非线性量测方程:
步骤2、定义辅助变量a0(1)、a0(2)、a0(3):
同时引入量测误差、雷达接收站站址误差和时钟偏差:
将公式(4)所示的辅助变量和公式(5)所述的三种误差代入非线性方程式(1)、(2)中,将其转化为伪线性方程:
构造联合未知向量θ和辅助变量的未知向量ψ=[θT,ao(1),ao(2),ao(3)]T。整合M个量测方程,建立联合目标位置、外辐射源位置、时钟偏差和辅助变量的定位模型:
h1-G1ψ=B1n+D1ΔV (7)
B1=diag(B1r,B1d)、D1=diag(D1r,D1d),diag(*)表示对角矩阵,其对角项为*中的元素。ΔV=[Δv1 T,Δv2 T,…,ΔvM T]表示雷达接收站站址误差,接收站位置的协方差矩阵Qs=E[ΔVΔVT]:
其中,0k表示长度为k的零向量。
步骤3.1、首先将B1和D1设为单位矩阵,利用公式(8)获得初始解,再通过初始解更新B1和D1:
步骤4、针对步骤3的估计结果,考虑辅助变量ao(1)、ao(2)、ao(3)与目标位置、外辐射源位置和时钟偏差的耦合性,并考虑步骤3得到的估计值的误差构建加权最小二乘估计模型并采用加权最小二乘算法优化求解;
步骤4.2、整合线性估计方程,得到关联线性估计模型:
h2-G2θ=B2Δψ (13)
步骤4.4、使用步骤3中的计算方法,首先将W2中的矩阵B2设为单位矩阵,利用式(14)获得初始解,再进一步通过初始解对矩阵B2进行更新,采用加权最小二乘估计算法不断迭代更新,得到最终的目标位置估计值:代表向量的第1项到第K项。
为了证明本方法的有效性,本实施例提出以下仿真情景,模拟外辐射源位置未知、时钟偏差和雷达接收站站址误差三类不确定性误差下目标定位过程:
设置在二维定位场景下,使用4个雷达接收站和单个位置未知的外辐射源组成多基外辐射源雷达定位系统,其中,四个接收站的位置分别为s1=[0,1000]T、s2=[1000,0]T、s3=[-1000,0]T、s4=[0,-1000]T,用于求解未知的目标位置、外辐射源位置和时钟偏差。设定目标和外辐射源的真实位置分别为u=[2000,5000]T、t=[3000,2000]T,时钟偏差bt=500。设置蒙特卡罗模拟计算仿真次数L=5000次,u(l)代表第l次蒙特卡罗模拟得到的目标辐射源位置参数估计值。
本实施例使用均方误差准则(MSE)来评估仿真实验的结果:
将MSE(u)绘制为与不同等级量测噪声方差相关的函数,并与克拉美罗下界CRLB进行比较,如图2所示,展示了在恒定接收站站址误差σs=0.1条件下,随着噪声功率σ2的增加,本实施例所得目标辐射源位置参数估计值的均方误差(MSE)估计精度。实验结果表明:基于外辐射位置未知、时钟偏差和雷达接收站站址误差三种不确定因素下,在联合量测噪声误差σ2不超过30dB时,本方法可以达到克拉美罗下界,实现对目标的精确定位,验证了本方法在复杂空间地貌环境下的实用性和有效性。值得注意的是,本方法还能够在较高的噪声水平下依然保持良好的定位效率和定位精度。
Claims (6)
1.一种基于非合作外辐射源的多基雷达定位方法,其特征在于:具体包括以下步骤:
步骤1、雷达接收站接收来自第三方外辐射源的直达波信号与目标散射后的第三方外辐射源信号,获得雷达接收站与第三方外辐射源之间的直接距离d和间接距离r,考虑未知的时间偏移量和量测噪声,建立非线性量测方程;
步骤2、考虑量测误差、雷达接收站站址误差和时钟偏差,引入与目标位置、外辐射源位置和时钟偏差相关的辅助变量a0(1)、a0(2)、a0(3),将步骤1建立的非线性量测方程转化为伪线性方程,建立定位模型;
2.如权利要求1所述一种基于非合作外辐射源的多基雷达定位方法,其特征在于:所述直接距离d和间接距离r根据接收的第三方外辐射源信号上的时间戳确定,第i个雷达接收站的量测方程为:
其中,表示第i个雷达接收站与第三方外辐射源之间的直接距离,表示第i个雷达接收站与第三方外辐射源之间的间接距离,表示第i个雷达接收站的真是位置;i=1,2,…,M,M为多基雷达定位系统中雷达接收站的数量;uo=[x,y]T表示固定目标的真实位置,to=[xt,yt]T表示第三方外辐射源的真实位置;||*||表示计算欧式距离;
4.如权利要求2或3所述一种基于非合作外辐射源的多基雷达定位方法,其特征在于:考虑量测误差、雷达接收站站址误差和时钟偏差:
将公式(4)所示的误差代入公式(2)所示的非线性量测方程,并通过辅助变量将其转化为伪线性方程:
h1-C1ψ=B1n+D1ΔV
(6)
B1=diag(B1r,B1d);D1=diag(D1r,D1d),diag(*)表示对角矩阵,其对角项为*中的元素;ΔV=[Δv1 T,Δv2 T,…,ΔvM T]表示雷达接收站站址误差,接收站位置的协方差矩阵为Qs=E[ΔVΔVT]:
其中,0k表示长度为k的零向量。
6.如权利要求4所述一种基于非合作外辐射源的多基雷达定位方法,其特征在于:步骤4的具体过程为:
步骤4.2、整合线性估计方程,得到关联线性估计模型:
h2-C2θ=B2Aψ (12)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310214713.7A CN116148829A (zh) | 2023-03-08 | 2023-03-08 | 一种基于非合作外辐射源的多基雷达定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310214713.7A CN116148829A (zh) | 2023-03-08 | 2023-03-08 | 一种基于非合作外辐射源的多基雷达定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116148829A true CN116148829A (zh) | 2023-05-23 |
Family
ID=86339016
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310214713.7A Pending CN116148829A (zh) | 2023-03-08 | 2023-03-08 | 一种基于非合作外辐射源的多基雷达定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116148829A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116660855A (zh) * | 2023-07-31 | 2023-08-29 | 中国科学院空天信息创新研究院 | 基于多非协作式辐射源的动态三维空间目标探测定位方法 |
CN117092588A (zh) * | 2023-10-20 | 2023-11-21 | 中国科学院深海科学与工程研究所 | 一种估算水声定位系统时钟偏差方法 |
-
2023
- 2023-03-08 CN CN202310214713.7A patent/CN116148829A/zh active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116660855A (zh) * | 2023-07-31 | 2023-08-29 | 中国科学院空天信息创新研究院 | 基于多非协作式辐射源的动态三维空间目标探测定位方法 |
CN116660855B (zh) * | 2023-07-31 | 2023-09-29 | 中国科学院空天信息创新研究院 | 基于多非协作式辐射源的动态三维空间目标探测定位方法 |
CN117092588A (zh) * | 2023-10-20 | 2023-11-21 | 中国科学院深海科学与工程研究所 | 一种估算水声定位系统时钟偏差方法 |
CN117092588B (zh) * | 2023-10-20 | 2024-01-09 | 中国科学院深海科学与工程研究所 | 一种估算水声定位系统时钟偏差方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108717184B (zh) | 基于误差校正的联合doa与toa单站无源定位方法 | |
CN109633592B (zh) | 运动观测站误差下外辐射源雷达时差与频差协同定位方法 | |
Zhang et al. | Multistatic localization in the absence of transmitter position | |
CN107132505B (zh) | 直达与非直达混合场景中的多目标直接定位方法 | |
CN108387876B (zh) | 基于ctls的外辐射源雷达网双基距误差配准方法 | |
CN116148829A (zh) | 一种基于非合作外辐射源的多基雷达定位方法 | |
CN105607092B (zh) | 基于tdoa和功率测量值的gnss欺骗干扰定位方法 | |
CN104077498B (zh) | 一种结合目标角度的外辐射源雷达多目标跟踪方法 | |
CN107688179A (zh) | 基于多普勒信息辅助的综合概率数据互联方法 | |
CN112162244B (zh) | 一种相关噪声和随机丢包环境下的事件触发目标跟踪方法 | |
CN105137424A (zh) | 一种杂波背景下实波束扫描雷达角超分辨方法 | |
WO2016112758A1 (zh) | 终端的定位方法及装置 | |
CN105487072A (zh) | 一种基于t2/r的时差与多普勒频移联合定位方法及系统 | |
Atia et al. | A consistent zero-configuration GPS-Like indoor positioning system based on signal strength in IEEE 802.11 networks | |
CN116321410A (zh) | 一种基于toa的迭代加权最小二乘超宽带定位算法 | |
Ristic et al. | Recursive Bayesian state estimation from Doppler-shift measurements | |
Schmitz et al. | TDOA fingerprinting for localization in non-line-of-sight and multipath environments | |
Ren et al. | An improved indoor positioning algorithm based on RSSI filtering | |
CN116930963A (zh) | 一种基于无线通信系统的穿墙成像方法 | |
CN110536410B (zh) | 非视距环境下基于rss和tdoa测量的定位方法 | |
Hao et al. | Joint source localisation and sensor refinement using time differences of arrival and frequency differences of arrival | |
Huang et al. | 3D TDOA/AOA localization in MIMO passive radar with transmitter and receiver position errors | |
CN107219498B (zh) | 基于mml-pmht的多基站单频网无源协同定位方法 | |
CN114035182B (zh) | 一种基于电离层反射的多站时差多变量短波目标定位方法 | |
CN109991564A (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 |