CN109541540A - 一种运动单站伪距速度联合定位方法 - Google Patents
一种运动单站伪距速度联合定位方法 Download PDFInfo
- Publication number
- CN109541540A CN109541540A CN201811572118.6A CN201811572118A CN109541540A CN 109541540 A CN109541540 A CN 109541540A CN 201811572118 A CN201811572118 A CN 201811572118A CN 109541540 A CN109541540 A CN 109541540A
- Authority
- CN
- China
- Prior art keywords
- toa
- pseudorange
- vector
- velocity
- change rate
- 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.)
- Granted
Links
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/06—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
-
- 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/0009—Transmission of position information to remote stations
- G01S5/0045—Transmission from base station to mobile station
- G01S5/0054—Transmission from base station to mobile station of actual mobile position, i.e. position calculation on base station
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种运动单站伪距速度联合定位方法,只需单个接收天线和处理通道,只需要进行TOA测量,通过分段数据处理,分别提取伪距和速度进行定位。该方法利用运动单站相对目标的位置和速度变化信息,而目前常见的方法通常仅利用到相对位置信息,不利用相对速度信息,或利用速度信息时需要对全部数据进行批处理,计算量较大。该方法基于定位模型分析,将位置信息分解为伪距和速度信息,同时解决了提高定位精度和降低计算量的问题,定位精度优于其它方法,等价于已有的虚拟时差方法,且复杂度大大降低。本发明的方法具有系统简单、设备量小、易于实现等优点,对无源定位具有重要的参考价值。
Description
技术领域
本发明属于无源定位领域,具体涉及一种运动单站伪距速度联合定位方法。
背景技术
无源定位技术通过接收目标辐射源信号确定其位置,具有独特的技术优势,具有十分重要的作用。单站无源定位由于不存在多站协同、数据传输、时频同步等问题,得到了广泛的研究和应用。常见的单站定位技术以干涉仪体制为主,包括测向定位、角度相位差联合定位、测相位差变化率定位等,由于需要多个接收天线和通道,并形成较长的基线,因此系统体积大、代价高,且对接收站的姿态测量精度要求较高。为此利用单个接收通道进行定位的方法近年来得到了越来越多的研究。
TOA是最为容易获得的单通道观测量,且估计方法成熟。文献(周龙健,罗景青,孔辉.基于虚拟时差的运动阵列空间无源定位算法[J].电子与信息学报,2017,39(7):1759-1763.)建立了包含相对运动速度信息的TOA定位模型,并提出了虚拟时差定位方法(下面称之为STOA-DPD方法)和最大似然定位算法。该方法的定位模型建立一一对应的脉冲TOA序列与观测站位置序列的基础上,而实际场景下,观测站的导航数据率通常远低于TOA序列的数据率,无法满足上述模型的求解条件。
本发明采用与文献(周龙健,罗景青,孔辉.基于虚拟时差的运动阵列空间无源定位算法[J].电子与信息学报,2017,39(7):1759-1763.)相同的假设,即在较短的观测时间内,观测站与目标的相对速度保持不变。在此基础上,利用短时观测得到的TOA序列估计得到伪距和速度进行定位。该发明实现了观测数据的压缩,同时避免了导航数据率与TOA数据率不匹配的问题。
发明内容
本发明的目的在于提供一种运动单站伪距速度联合定位方法,只需要单个接收天线和通道进行TOA测量,算法简单,易于实现,定位精度优于其它方法,等价于已有的虚拟时差方法且复杂度大大降低,具有一定的理论和工程应用价值。
实现本发明目的的技术解决方案为:一种运动单站伪距速度联合定位方法,包括以下步骤:
步骤一:利用单个运动观测站,持续接收目标辐射源脉冲串信号,并估计TOA,得到TOA序列,将得到的序列分为连续的M个时间片,每个时间片内包含N次对TOA的连续观测,收到第m个时间片内第n个脉冲的时刻记为tm,n,记各个时间片内观测站的位置为、速度矢量分别为Pm、vm,m=1,…,M,n=1,…,N,粗估计得到脉冲重复周期Tr0=t1,2-t1,1;
步骤二:根据接收到的TOA序列,估计出各个时间片内的TOA变化率
步骤三:利用估计得到的TOA变化率采用BFGS算法,得到目标位置的估计
本发明与现有技术相比,其显著优点在于:
(1)只需单个接收天线和处理通道,系统简单,设备量小;
(2)只需要进行TOA测量,易于实现;
(3)定位精度优于其它方法,等价于已有的虚拟时差方法且复杂度大大降低。
附图说明
图1为TOA变化率定位原理图。
图2为斜率和截距包含了完整的TOA序列信息图。
图3为不同TOA测量精度下TOA/TOA变化率联合方法与虚拟时差定位方法的比较图。
图4为观测子段分段时长对TOA/TOA变化率联合方法的影响图。
图5为本发明运动单站伪距速度联合定位方法的流程图。
具体实施方式
下面结合附图对本发明作进一步详细描述。
1、定位观测模型
结合图1和图5所示,假设静止目标以恒定的重复周期(pulse repetitioninterval,PRI)辐射脉冲信号,观测站采用单个通道接收目标信号并估计出脉冲到达时间(time of arrival,TOA)序列。一段时间内,观测站沿一定的航迹运动,在空间上形成虚拟的观测孔径。由于多普勒效应,在A、B、C等不同时刻,观测站接收到的脉冲间隔不等于脉冲发出的间隔,而是与观测站的速度和观测站与目标的相对位置有关,表现为当观测站接近目标时,脉冲间隔小于脉冲发出的重复周期,反之则大于脉冲发出的重复周期。因而脉冲到达时间以及脉冲的间隔包含了相对位置信息,如果持续观测目标信号,尽管单个脉冲无法确定脉冲的位置,但通过观测不同时刻脉冲到达时间和脉冲间隔的变化,可以实现定位。
1.1、TOA矢量模型
记未知的脉冲重复周期为Tr,电磁波传播速度记为c。目标位置为接收到第n个脉冲时的观测站的位置为
TOA基本模型中包含了观测站位置,但却没有包含观测站的速度。根据已有的研究可以知道,如果利用观测站的速度,可以提高定位的精度。下面将从速度矢量出发建立TOA观测模型。以第一个脉冲的发出时刻为时间零点,假定初始时刻观测站位置为观测站的瞬时速度矢量为假定观测站速度保持不变,则第n个脉冲时观测站与目标的距离记为rn
如果以第一个脉冲的发出时刻为时间零点,则有
将式(1)代入式(2)得到从而求解可以得到
这一模型可以刻画观测站以特定速度运动时,观测得到的TOA序列与相对位置和相对速度的关系,在本文的仿真分析中,均采用上述矢量模型生成TOA序列。
1.2、TOA标量模型
由于切向速度分量不影响TOA的变化,因此只可以考虑径向速度分量的影响。记片内第一个脉冲的发出时刻观测站与目标的距离、瞬时相对径向速度分别为r0、在一个时间片内,认为瞬时相对径向加速度不变。分别满足
式中,分别为观测站相对目标的位置矢量和距离,为观测站的位置。则有
将式(5)代入式(2)构成关于tn的二次方程,可解得
同时考虑到观测站无法获知脉冲发出的时间零点,即观测站与目标存在时钟偏差,记为Δt,得到
其中,t0为与时间无关的未知常量。
矢量模型和标量模型分别建立在速度不变和瞬时径向不变的基础上,但速度不变与瞬时径向速度不变是矛盾的。标量模型建立在较短的观测时间范围内,对于匀速直线运动的观测站,在无限长的观测时间段后,不能认为相对目标的瞬时径向速度保持不变。
1.3、联合定位模型
文献(徐义,郭福成,冯道旺.一种单星仅测TOA无源定位方法[J].宇航学报.2010,31(2):502-508.)采取了相同的分析思路,建立了存在相对运动加速度条件下的定位模型,通过估计TOA的二次变化率进行定位的方法。类似的,如图2所示,由于短时间内TOA呈一次变化的规律,TOA序列的截距和斜率即代表了该段TOA序列的全部信息,因而可以估计出各段观测中TOA截距和一次变化率进行定位。TOA截距即为初始时刻的TOA,也可以称为伪距,TOA变化率的本质则为相对速度,因此本文也称这种方法为伪距/速度联合定位方法。
截距本质上是初始的TOA,和传统TOA定位相比,截距定位是利用TOA序列估计出一个最为准确的TOA代表完整的TOA序列参与定位计算。因此,从一段TOA序列中估计出TOA变化率以及初始的TOA,即保留了完整的TOA序列信息,而无需直接处理数百个TOA。TOA序列截距和TOA变化率联合定位单时刻定位方程组为
其中, 分别为TOA序列的变化率和截距,仅与组序号m有关,与组内脉冲计数n无关。
2、TOA变化率估计及误差分析
根据最小二乘准则,可估计得到
式中,化简得到如果各次观测独立同零均值高斯分布,记TOA估计方差为σt,则可推知TOA变化率估计误差式为
继而得到估计误差的协方差矩阵
可以看出,和的估计误差并不是独立的,即为非对角矩阵,
对式(8)进行误差分析求得可以把TOA变化率估计误差等效为相对径向速度估计误差(下面简称“等效速度估计误差”)。结合的估计误差分析可以知道
因此:1)等效速度估计误差与TOA估计误差成正比;2)相对径向速度越大,等效速度估计误差越小,但由于远小于光速,径向速度大小的影响较小;3)在观测时间TrN一定的前提下,Tr越小、N越大,等效速度估计误差越小。
3、定位误差分析
未知的参数包括T、Tr、Δt,对式(10)微分得到误差分析式
其中,整理为矩阵形式
其中,从而得到
TOA/TOA变化率联合定位与文献0方法均是基于同样的观测模型,均利用了观测站的位置和速度信息,两种方法的参数估计和定位求解均是统计最优的,那么两种方法的定位精度必然也是相同的。由式(11)可以看出,TOA/TOA变化率定位方程组实质上是文献(周龙健,罗景青,孔辉.基于虚拟时差的运动阵列空间无源定位算法[J].电子与信息学报,2017,39(7):1759-1763.)方法定位方程组左乘了该矩阵与自身的转置相乘为单位矩阵时,定位误差保持不变。本质上来看,STOA-DPD方法可以看成是TOA/TOA变化率定位方法基于直接定位技术的实现。STOA-DPD方法更适合采用批处理的定位算法,而TOA/TOA变化率定位方法更适合采用传统的滤波跟踪算法。新方法实现了数据的压缩,而STOA-DPD方法需要保存全部的TOA序列。脉冲序列越长、PRI越小或相对速度越小、短时模型持续时间越长,越适合采用本文方法。
4、最大似然迭代定位算法
经过M次观测,整理为如式(18)所示的向量形式
式中,式(18)所示的方程组非线性程度较高,难以寻找解析方法,这里采用搜索算法求解,令 建立位置代价函数
式中, 为TOA截距及TOA变化率估计的协方差矩阵,,其元素由式(13)得到,采用BFGS算法进行迭代计算,对Tr、Δt同时进行求解。根据BFGS算法,需要计算梯度矩阵,梯度矩阵为其中,
迭代估计的步骤即首先选取粗估计值,根据粗估计值进行迭代,当定位精度满足要求时,停止迭代。
5、仿真试验
1)实施例一:不同TOA测量精度下的定位精度对比
图3对不同TOA估计精度下两种定位方法的性能进行了仿真对比分析,仿真条件为:观测站起始位置为[0,0],目标位置为[50,50]km,观测站速度为[250,0]m/s,观测站加速度为[5,0]m/s2,导航位置和速度数据率为0.5s,对应的观测站位置误差为80m,对应的观测站速度误差为0.2m。脉冲重复周期为1ms,持续观测80s进行定位,进行500次蒙塔卡洛仿真。由图4可以看出,TOA/TOA变化率联合定位方法定位精度等价于虚拟时差定位方法;提出的基于BFGS的最大似然定位算法可以逼近理论定位误差限。
2)实施例二:不同分段时长下的定位精度对比
下面仿真分析不同的观测子段长度对TOA/TOA变化率联合定位精度的影响。分段时长为0.25~2.5s,其它仿真条件与图4对应的条件相同。导航位置和速度数据率为0.5s,在0.25s的分段长度下,对导航位置和速度进行相邻点一阶插值。当0.75s、1.25s等分段长度下子段长度作取整近似。如4给出了观测子段分段时长对TOA/TOA变化率联合方法的影响。
由图4可以看出,在观测子段时长小于1s时,定位误差可以逼近理论定位精度。而当观测子段时长越长时,定位误差越大,这是因为此时短时观测模型不再成立,从而引入了模型误差。而在观测子段时长为0.25s并对导航数据进行插值时,定位误差也大于理论值,其主要误差为插值引入的误差,整体上小于观测子段时长过长时引入的误差。由此可以看出,选择合适的观测子段长度尤其重要,在进行定位算法设计时,需要考虑短时观测模型的适用范围。
Claims (5)
1.一种运动单站伪距速度联合定位方法,其特征在于,包括以下步骤:
步骤一:利用单个运动观测站,持续接收目标辐射源脉冲串信号,并估计TOA,得到TOA序列,将得到的序列分为连续的M个时间片,每个时间片内包含N次对TOA的连续观测,收到第m个时间片内第n个脉冲的时刻记为tm,n,记各个时间片内观测站的位置为、速度矢量分别为Pm、vm,m=1,…,M,n=1,…,N,粗估计得到脉冲重复周期Tr0=t1,2-t1,1;
步骤二:根据接收到的TOA序列,估计出各个时间片内的TOA变化率
步骤三:利用估计得到的TOA变化率采用BFGS算法,得到目标位置的估计
2.根据权利要求1所述的运动单站伪距速度联合定位方法,其特征在于,在步骤二中,估计各个时间片内的TOA变化率的公式如下:
3.根据权利要求1所述的运动单站伪距速度联合定位方法,其特征在于,在步骤三中,BFGS算法如下:
步骤3-1:在观测范围内任意选择初始位置T0,选定时间偏差的粗估计值Δt0,转入步骤3-2;
步骤3-2:根据粗估计值Δt0和梯度矩阵迭代更新目标位置估计,得到位置估计值转入步骤3-3;
步骤3-3:如果ξ为所需的定位精度,则完成定位,目标位置最终估计值即为否则令返回步骤3-2。
4.根据权利要求3所述的运动单站伪距速度联合定位方法,其特征在于:步骤3-2中,梯度矩阵为W为TOA截距及TOA变化率估计的协方差矩阵,F(α)为定位估计的代价函数,α为待估计矢量,α=[TT Tr Δt]T,T为待估计的目标位置,Tr为待估计的脉冲重复周期,Δt为待估计的时间偏差,为TOA变化率关于位置的梯度向量, 为TOA变化率关于位置的梯度,rm、rm、vm、分别为观测站相对目标的位置矢量、距离、速度矢量、速度,c为光速,为TOA截距关于位置的梯度向量, 为TOA截距关于位置的梯度,km为时间压缩因子向量, 为TOA变化率关于脉冲重复周期的梯度向量,km为时间压缩因子向量,km=[k1 … kM],为TOA截距关于脉冲重复周期的梯度向量, 为TOA截距关于脉冲重复周期的梯度, 为TOA截距关于时间偏差的梯度向量,
5.根据权利要求4所述的运动单站伪距速度联合定位方法,其特征在于:TOA截距及TOA变化率估计的协方差矩阵W公式如下:
N为时间片内脉冲总个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811572118.6A CN109541540B (zh) | 2018-12-21 | 2018-12-21 | 一种运动单站伪距速度联合定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811572118.6A CN109541540B (zh) | 2018-12-21 | 2018-12-21 | 一种运动单站伪距速度联合定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109541540A true CN109541540A (zh) | 2019-03-29 |
CN109541540B CN109541540B (zh) | 2020-10-20 |
Family
ID=65856375
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811572118.6A Active CN109541540B (zh) | 2018-12-21 | 2018-12-21 | 一种运动单站伪距速度联合定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109541540B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109991570A (zh) * | 2019-04-05 | 2019-07-09 | 电子科技大学 | 一种基于拟牛顿法的运动目标直接定位方法 |
CN116184449A (zh) * | 2022-11-09 | 2023-05-30 | 湖南六九零六信息科技股份有限公司 | 一种基于单星分时的时差定位方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102495392A (zh) * | 2011-12-05 | 2012-06-13 | 中国电子科技集团公司第三十八研究所 | 高精度测量脉冲到达时间的方法和装置 |
CN106597364A (zh) * | 2016-11-18 | 2017-04-26 | 烟台职业学院 | 一种单天线单站无源定位的目标辐射源初始位置估计方法 |
US20180093156A1 (en) * | 2013-06-14 | 2018-04-05 | Vg Buyer, Llc | Enhanced golf simulation system |
US20180113160A1 (en) * | 2016-10-24 | 2018-04-26 | Xia Llc | Ratio-Reference Measurement of the Arrival Time and/or Amplitude of a Digitized Electronic Pulse |
CN108120954A (zh) * | 2017-12-21 | 2018-06-05 | 中国航天科工集团八五研究所 | 一种高精度toa变化率定位方法 |
CN108205121A (zh) * | 2016-12-20 | 2018-06-26 | 中国航天科工集团八五研究所 | 一种基于toa变化率的单通道无源定位方法 |
-
2018
- 2018-12-21 CN CN201811572118.6A patent/CN109541540B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102495392A (zh) * | 2011-12-05 | 2012-06-13 | 中国电子科技集团公司第三十八研究所 | 高精度测量脉冲到达时间的方法和装置 |
US20180093156A1 (en) * | 2013-06-14 | 2018-04-05 | Vg Buyer, Llc | Enhanced golf simulation system |
US20180113160A1 (en) * | 2016-10-24 | 2018-04-26 | Xia Llc | Ratio-Reference Measurement of the Arrival Time and/or Amplitude of a Digitized Electronic Pulse |
CN106597364A (zh) * | 2016-11-18 | 2017-04-26 | 烟台职业学院 | 一种单天线单站无源定位的目标辐射源初始位置估计方法 |
CN108205121A (zh) * | 2016-12-20 | 2018-06-26 | 中国航天科工集团八五研究所 | 一种基于toa变化率的单通道无源定位方法 |
CN108120954A (zh) * | 2017-12-21 | 2018-06-05 | 中国航天科工集团八五研究所 | 一种高精度toa变化率定位方法 |
Non-Patent Citations (1)
Title |
---|
周龙健 等: "基于虚拟时差的运动阵列空间无源定位算法", 《电子与信息学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109991570A (zh) * | 2019-04-05 | 2019-07-09 | 电子科技大学 | 一种基于拟牛顿法的运动目标直接定位方法 |
CN109991570B (zh) * | 2019-04-05 | 2023-02-03 | 电子科技大学 | 一种基于拟牛顿法的运动目标直接定位方法 |
CN116184449A (zh) * | 2022-11-09 | 2023-05-30 | 湖南六九零六信息科技股份有限公司 | 一种基于单星分时的时差定位方法 |
CN116184449B (zh) * | 2022-11-09 | 2023-11-10 | 湖南六九零六信息科技股份有限公司 | 一种基于单星分时的时差定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109541540B (zh) | 2020-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Musicki et al. | Geolocation using TDOA and FDOA measurements | |
CN108872971B (zh) | 一种基于运动单阵列的目标定位方法与装置 | |
EP2817652B1 (en) | Method and system for simultaneous receiver calibration and object localisation for multilateration | |
CN108872932B (zh) | 基于神经网络的超视距目标直接定位结果纠偏方法 | |
CN109407126A (zh) | 一种多模接收机联合定位解算的方法 | |
CN105044667A (zh) | 一种运动目标的双星跟踪方法、装置和系统 | |
US5610609A (en) | Passive doppler emitter classification and ranging system utilizing a time-frequency constant | |
CN103363994A (zh) | 仅基于无线电载波相位观测的卫星精密定轨技术 | |
Zhu et al. | Direct localization based on motion analysis of single-station using TOA | |
CN109541540A (zh) | 一种运动单站伪距速度联合定位方法 | |
CN104820221B (zh) | 多基合成孔径雷达的目标三维定位方法 | |
Zhang et al. | Emitter localization algorithm based on passive synthetic aperture | |
Grabbe et al. | Geo-location using direction finding angles | |
CN110471029A (zh) | 一种基于扩展卡尔曼滤波的单站无源定位方法及装置 | |
Davey et al. | Detection and tracking of multipath targets in over-the-horizon radar | |
CN106597421A (zh) | 异地天线阵中基于预测模型的时延和时延率快速搜索方法 | |
CN105548959A (zh) | 一种基于稀疏重建的多传感器多目标的定位方法 | |
Zhao et al. | Sequential doppler-shift-based optimal localization and synchronization with TOA | |
Stefanski | Asynchronous wide area multilateration system | |
CN108205121B (zh) | 一种基于toa变化率的单通道无源定位方法 | |
Kumar et al. | Performance analysis of TDOA and FDOA for missile tracking application using extended Kalman filter | |
Bolelov et al. | A Study of Aircraft Positioning Precision in a MLAT Surveillance System with Different Flight Paths and Ground Station Layouts | |
Li et al. | Performance analysis for space-air-ground integrated passive localization using TDOA measurements | |
Mazidi | Introducing new localization and positioning system for aerial vehicles | |
Komaragiri et al. | A sag monitoring device based on a cluster of code based GPS receivers |
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 |