CN112526449B - 一种利用运动目标对接收站位置信息进行校准的方法 - Google Patents
一种利用运动目标对接收站位置信息进行校准的方法 Download PDFInfo
- Publication number
- CN112526449B CN112526449B CN202011365091.0A CN202011365091A CN112526449B CN 112526449 B CN112526449 B CN 112526449B CN 202011365091 A CN202011365091 A CN 202011365091A CN 112526449 B CN112526449 B CN 112526449B
- Authority
- CN
- China
- Prior art keywords
- receiving station
- positioning
- self
- target
- information
- 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
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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/10—Position of receiver fixed by co-ordinating a plurality of position lines defined by path-difference measurements, e.g. omega or decca systems
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供一种利用运动目标对接收站位置信息进行校准的方法,包括:步骤一、根据目标定位测量值、目标定位估计值和接收站自身定位信息,构造一个关于接收站定位误差的第一方程组;步骤二、将接收站定位误差统计信息纳入所述第一方程组,形成第二方程组;步骤三、求解所述第二方程组,得到WLS解,由WLS解得接收站自身定位估计误差,利用所述接收站自身定位估计误差对接收站自身定位信息进行校准。本发明利用运动目标来对接收站位置信息进行校准时,将目标定位估计值及其估计误差同时纳入考虑范围,通过综合优化,提升了算法的鲁棒性,在较高的噪声环境下,接收站自身的定位精度也能达到CRLB。
Description
技术领域
本发明涉及通信辐射源无源定位跟踪领域,具体是一种利用运动目标对接收站位置信息进行校准的方法。
背景技术
在现代的无源定位系统中,无源定位接收站通常布设在移动平台上,比如无人机、飞机和汽车。由于目标和接收站之间存在相对运动,移动平台上的接收站可以同时获得目标的到达时间差(Time Difference of Arrival,TDOA)[1]和到达频率差(FrequancyDifference of Arrival,FDOA)[2]测量值,对目标的位置和速度进行联合估计[3]。但是,在某些环境下,移动平台自身的位置信息并不精确,造成接收站自身定位信息存在误差。即使接收站自身定位误差比较小,但其通常会导致目标定位精度急剧恶化[4]。为了保持目标定位精度,需要对接收站的位置信息进行校准。
一些研究人员采用引入校准站的方法[5-10],来对接收站的位置信息进行校准。通过仿真发现,该方法能够有效地提高接收站的定位精度,但是,校准站在实际的环境中并不总是存在的。近年来,有提出利用运动目标来对接收站位置信息进行校准,通过先对目标进行定位,再利用目标定位的估计值校准接收站位置信息,该方法仅在低噪声的环境下,才能取得较好的效果[11,12]。
参考文献:
[1]Ho K C.Bias reduction for an explicit solution of sourcelocalization using TDOA.IEEE Transactions on Signal Processing 2012;60(5):2101-2114.
[2]Yu H G,Huang G M,Gao J.Constrained total least-squareslocalisation algorithm using time difference of arrival and frequencydifference of arrival measurements with sensor location uncertainties.IETRadar,Sensor and Navigation 2012;6(9):891-899.
[3]Ho K C,Xu W.An accurate algebraic solution for moving sourcelocation using TDOA and FDOA measurements.IEEE Transactions on SignalProcessing 2004;52(9):2453-2463.
[4]Ho K C,Lu X,Kovavisaruch L.Source localization using TDOA and FDOAmeasurements in the presence of receiver location errors:Analysis andsolution.IEEE Transactions on Signal Processing 2007;55(2):684-696.
[5]Ho K C,Yang L.On the use of a calibration emitter for sourcelocalization in the presence of sensor position uncertainty.IEEE Transactionson Signal Processing 2008;56(12):5758-5772.
[6]Yang L,Ho K C.Alleviating sensor position error in sourcelocalization using calibration emitters at inaccurate locations.IEEETransactions on Signal Processing 2010;58(1):67-83.
[7]Yang L,Ho K C.On using multiple calibration emitters and theirgeometric effects for removing sensor position errors in TDOAlocalization.Proceedings IEEE International Conference on Acoustics Speechand Signal Processing.2010;2702-2705.
[8]Li J Z,Guo F C,Jiang W L.Source localization and calibration usingTDOA and FDOA measurements in the presence of sensor locationuncertainty.Science China Information Sciences 2014;57(4):1-12.
[9]Li J Z,Guo F C,Le Y,et al.On the use of calibration sensors insource localization using TDOA and FDOA measurements.Digital SignalProcessing 2014;27:33-43.
[10]Li J Z,Guo F C,Jiang W L.Source localization and calibrationusing TDOA and FDOA measurements in the presence of sensor locationuncertainty.Science China Information Sciences 2014;57(4):1-12.
[11]Sun M,Ho K C.Refining inaccurate sensor positions using target atunknown location.Signal Processing 2012;92(9):2097-2104.
[12]B Hao,Z Li,J Si,et al.Joint source localisation and sensorrefinement using time differences of arrival and frequency differences ofarrival.IET Signal Processing 2014;8(6):588-600.
发明内容
针对现有技术存在的问题,本发明提供一种利用运动目标对接收站位置信息进行校准的方法,利用运动目标来对接收站位置信息进行校准时,将目标定位估计值及其估计误差同时纳入考虑范围,通过综合优化,提升了算法的鲁棒性,在较高的噪声环境下,接收站自身的定位精度也能达到克拉美罗下界(Cramer-Rao Lower Bound,CRLB)。
一种利用运动目标对接收站位置信息进行校准的方法,包括如下步骤:
步骤一、根据目标定位测量值、目标定位估计值和接收站自身定位信息,构造一个关于接收站定位误差的第一方程组;
步骤二、将接收站定位误差统计信息纳入所述第一方程组,形成第二方程组;
步骤三、求解所述第二方程组,得到WLS解,由WLS解得接收站自身定位估计误差,利用所述接收站自身定位估计误差对接收站自身定位信息进行校准。
进一步的,假设有M个接收站,它们位置和速度的真实值未知,含噪值为其中,i=1,2,...,M,和分别是位置和速度的真实值,Δsi和分别是位置和速度估计误差,将si、构造成一个接收站自身定位向量其中, 是接收站自身定位误差向量,假设其服从均值为零,协方差为Qβ的高斯分布,假设待定位目标的位置和速度分别为uo和
将接收站1设置为参考站,目标到第i个接收站与参考站的距离差和距离差对时间的微分真实值分别为:
利用接收站自身定位值β和测量值α可实现目标定位,目标定位估计值为
θ=θo+Δθ (16)
所述步骤一具体包括:
通过利用TDOA测量值信息、目标定位估计值和接收站的位置信息,得到:
为充分利用FDOA信息,将公式(30)对时间进行微分,得到:
其中
其中
进一步的,所述步骤二具体包括:将所述第一方程组,即公式(32),与-Δβ=-Δβ进行联立,得到所述第二方程组:
公式(34)的WLS解为:
其中
进一步的,所述步骤三具体为:
本发明先对目标进行定位,然后利用目标的定位估计值及其定位误差,结合目标的TDOA/FDOA测量值,改善接收站自身定位精度。该方法仅需要采用一次WLS算法。通过计算机模拟仿真,发现在较高的噪声环境下,该算法对接收站自身定位的估计精度仍然能够达到CRLB,并且,利用提升精度后的接收站位置信息,可以有效提升后续目标的估计精度。
本专利申请受国家自然科学基金项目(61501484)资助支持。
附图说明
图1为本发明接收站与目标分布图;
图2为本发明与现有技术利用运动目标对接收站位置进行校准的精度比较;
图3为本发明与现有技术利用运动目标对接收站速度进行校准的精度比较;
图4为本发明与现有技术利用校准之后的接收站对后续目标定位的位置估计精度比较;
图5为本发明与现有技术利用校准之后的接收站对后续目标定位的速度估计精度比较。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
本申请提出了一种利用运动目标对接收站位置信息进行校准的方法,该方法将待定位目标作为接收站的校准站,通过利用目标的定位估计值和TDOA/FDOA测量值,并充分考虑目标定位误差对接收站校准的影响,实现了对接收站的校准。所提方法具有闭合解的形式,经计算机仿真,发现该算法在接收站定位误差较大时,仍能达到CRLB。
假设有M个接收站,它们位置和速度的真实值是未知的,仅知道它们的含噪值为其中,i=1,2,...,M,和分别是位置和速度的真实值,Δsi和分别是位置和速度估计误差。将si、构造成一个接收站自身定位向量其中,此处,是接收站自身定位误差向量,假设其服从均值为零,协方差为Qβ的高斯分布。假设待定位目标的位置和速度分别为uo和
将接收站1设置为参考站。目标到第i个接收站与参考站的距离差和距离差对时间的微分真实值分别为:
利用接收站自身定位值β和测量值α可以实现目标定位,目标定位估计值为
θ=θo+Δθ (29)
为了实现对接收站自身定位精度进行校准,最直接的方法是构建一个包含接收站自身定位误差的方程组,通过求解该方程组,可以得到接收站自身定位误差。为了进一步提升校准精度,本发明在构建该方程组后,充分考虑了接收站自身定位误差的统计信息,将它纳入方程组,求解出接收站自身定位误差,并进行校准。本发明包括如下步骤:
步骤一,根据目标定位测量值、目标定位估计值和接收站自身定位信息,构造一个关于接收站定位误差的第一方程组。
通过利用TDOA测量值信息、目标定位估计值和接收站的位置信息,可以得到:
为了充分利用FDOA信息,将公式(30)对时间进行微分,可以得到:
其中
其中
步骤二,将接收站定位误差统计信息纳入所述第一方程组,形成第二方程组。
由于公式(32)未利用到接收站定位误差Δβ统计特性,将公式(32)与-Δβ=-Δβ进行联立,得:
公式(34)的WLS解为:
其中
步骤三,求解所述第二方程组,得到WLS解,由WLS解得接收站自身定位估计误差,利用所述接收站自身定位估计误差对接收站自身定位信息进行校准。
因此,可以对接收站定位进行校准,得:
计算机模拟:
该部分将进行两组计算机仿真,分别对本申请所提算法的接收站自身位置校准精度和对后续目标的定位能力进行评估。接收站的布置如图1所示。有6个接收站、1个待定位目标(用于接收站位置信息校准)和1个后续待定位目标(用于检验校准后的精度提升),其真实的位置和速度如表1所示。
表1接收站和目标的位置和速度
通过引入零均值,协方差分别为Qβ和Qα的高斯白噪声,生成含噪的接收站自身定位信息β和TDOA/FDOA测量值α。其中,,Qα=diag{Jt,Jf},Jt=10-3×J,Jf=10-5×J,J是一个对角线元素为1,其余元素为0.5的(M-1)×(M-1)维矩阵,为噪声方差。接收站自身定位校准能力和后续目标的定位精度通过计算均方误差来进行评估。
仿真过程如下:首先,利用含噪的接收站位置信息和TDOA/FDOA测量值对待定位目标进行位置估计;然后,利用目标定位结果结合TDOA/FDOA测量值来改善接收站自身的定位精度;最后,利用自身定位精度改善后的接收站对后续目标进行定位。
模拟仿真一:接收站校准
图2和图3比较了本申请所提算法与文献[11,12]所提算法对接收站自身定位精度的改善情况。从图中可以看出:本申请所提算法相较于接收站未校准前,定位精度提升了大概1dB;3种接收站校准算法在噪声方差小于-8dB时,都能够达到CRLB,但随着噪声增大,其余两种算法分别在噪声方差等于-4dB和12dB处,偏离了CRLB;通过比较,本申请所提算法在噪声方差等于20dB时,仅稍稍偏离CRLB,且仍小于或等于未经过校准的接收站定位误差。
模拟仿真二:后续目标定位
对接收站定位误差校准的目的就是为了提高对后续目标的定位精度。图4和图5分别比较了使用校准之后的接收站和未校准的接收站对目标定位的位置和速度估计误差。从图4和图5可以看出:由于进行了接收站定位误差校准,目标位置和速度的估计精度相比未校准前分别提高了4dB和5dB;相较于其他两种算法,本申请所提算法在较高接收站定位误差的条件下,具有相对较小的定位误差值。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何属于本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (1)
1.一种利用运动目标对接收站位置信息进行校准的方法,其特征在于:包括如下步骤:
步骤一、根据目标定位测量值、目标定位估计值和接收站自身定位信息,构造一个关于接收站定位误差的第一方程组;
步骤二、将接收站定位误差统计信息纳入所述第一方程组,形成第二方程组;
步骤三、求解所述第二方程组,得到WLS解,由WLS解得接收站自身定位估计误差,利用所述接收站自身定位估计误差对接收站自身定位信息进行校准;
假设有M个接收站,它们位置和速度的真实值未知,含噪值为其中,i=1,2,...,M,和分别是位置和速度的真实值,Δsi和分别是位置和速度估计误差,将si、构造成一个接收站自身定位向量其中, 是接收站自身定位误差向量,假设其服从均值为零,协方差为Qβ的高斯分布,假设待定位目标的位置和速度分别为uo和
将接收站1设置为参考站,目标到第i个接收站与参考站的距离差和距离差对时间的微分真实值分别为:
利用接收站自身定位值β和测量值α可实现目标定位,目标定位估计值为
θ=θo+Δθ (3)
所述步骤一具体包括:
通过利用TDOA测量值信息、目标定位估计值和接收站的位置信息,得到:
为充分利用FDOA信息,将公式(4)对时间进行微分,得到:
其中
其中
所述步骤二具体包括:将所述第一方程组,即公式(6),与-Δβ=-Δβ进行联立,得到所述第二方程组:
公式(8)的WLS解为:
其中
所述步骤三具体为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011365091.0A CN112526449B (zh) | 2020-11-27 | 2020-11-27 | 一种利用运动目标对接收站位置信息进行校准的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011365091.0A CN112526449B (zh) | 2020-11-27 | 2020-11-27 | 一种利用运动目标对接收站位置信息进行校准的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112526449A CN112526449A (zh) | 2021-03-19 |
CN112526449B true CN112526449B (zh) | 2022-12-27 |
Family
ID=74994792
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011365091.0A Active CN112526449B (zh) | 2020-11-27 | 2020-11-27 | 一种利用运动目标对接收站位置信息进行校准的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112526449B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000304854A (ja) * | 1999-04-26 | 2000-11-02 | Mitsubishi Electric Corp | 目標追尾方法及び目標追尾装置 |
CN110389327A (zh) * | 2019-07-29 | 2019-10-29 | 杭州电子科技大学 | 接收站位置误差下多站多外辐射源雷达双基距定位方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102017209002A1 (de) * | 2017-05-29 | 2018-11-29 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e. V. | Vorrichtung und Verfahren zum Ermitteln von relativen Positionen von Empfängerstationen eines Positionsbestimmungssystems |
CN109633581B (zh) * | 2018-12-26 | 2020-11-06 | 杭州电子科技大学 | 基于外辐射源tdoa/fdoa误差校正下的定位方法 |
CN109633592B (zh) * | 2019-01-18 | 2020-11-17 | 杭州电子科技大学 | 运动观测站误差下外辐射源雷达时差与频差协同定位方法 |
-
2020
- 2020-11-27 CN CN202011365091.0A patent/CN112526449B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000304854A (ja) * | 1999-04-26 | 2000-11-02 | Mitsubishi Electric Corp | 目標追尾方法及び目標追尾装置 |
CN110389327A (zh) * | 2019-07-29 | 2019-10-29 | 杭州电子科技大学 | 接收站位置误差下多站多外辐射源雷达双基距定位方法 |
Non-Patent Citations (3)
Title |
---|
一种TDOA_FDOA联合校准定位算法;周成 等;《海军工程大学学报》;20150831;第27卷(第4期);第45-49页 * |
传感器参数误差下的运动目标TDOA/FDOA无源定位算法;孙霆 等;《航空学报》;20200225;第41卷(第2期);第323317-1--323317-10页 * |
考虑站址误差的稳健TDOA定位算法;高向颖等;《信号处理》;20200825(第08期);第1196-1204页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112526449A (zh) | 2021-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108490473B (zh) | 一种融合gnss和uwb的无人机增强定位方法与系统 | |
CN108307301B (zh) | 基于rssi测距和轨迹相似性的室内定位方法 | |
KR101914550B1 (ko) | 레이더의 표적 위치 추적 방법 | |
CN109917333B (zh) | 融合aoa观测量与tdoa观测量的无源定位方法 | |
CN113074739A (zh) | 基于动态鲁棒容积卡尔曼的uwb/ins融合定位方法 | |
CN107708204B (zh) | 基于卡尔曼滤波的uwb定位系统基站自标定方法 | |
CN110673089B (zh) | 未知视距和非视距分布情况下基于到达时间的定位方法 | |
WO2006059155A1 (en) | Emitter geolocation | |
CN110231638A (zh) | 一种无线网络下的载波相位定位方法及装置 | |
CN110132281B (zh) | 一种基于询问应答模式的水下高速目标高精度自主声学导航方法 | |
CN110673181B (zh) | 一种基于网格能量遍历搜索的gnss干扰源定位方法 | |
CN107371133B (zh) | 一种提高基站定位精度的方法 | |
CN110636436A (zh) | 基于改进chan算法的三维uwb室内定位方法 | |
CN110657806A (zh) | 一种基于CKF、chan解算和Savitzky-Golay平滑滤波的位置解算方法 | |
Lee et al. | TDOA-/FDOA-based adaptive active target localization using iterated dual-EKF algorithm | |
CN108445446B (zh) | 一种无源测速定位方法及装置 | |
CN112526449B (zh) | 一种利用运动目标对接收站位置信息进行校准的方法 | |
CN111624549B (zh) | 非共视条件下无源滤波跟踪方法 | |
CN113093241A (zh) | 一种考虑高度角的单测站对流层斜延迟计算方法 | |
CN112835020B (zh) | 面向非视距参数估计的刚体定位方法 | |
CN110673088B (zh) | 混合视距和非视距环境中基于到达时间的目标定位方法 | |
Liang et al. | Application of Taylor-Chan algorithm based on TDOA in sound source location | |
CN110996248B (zh) | 一种基于rss的凸优化无线传感器网络定位方法 | |
CN111103573B (zh) | 一种基于卫星外辐射源雷达量测tdoa的折射修正方法 | |
Sung et al. | TDoA based UGV localization using adaptive Kalman filter algorithm |
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 |