CN106353720A - 基于tdoa/groa的多站连续定位模型 - Google Patents
基于tdoa/groa的多站连续定位模型 Download PDFInfo
- Publication number
- CN106353720A CN106353720A CN201610800756.3A CN201610800756A CN106353720A CN 106353720 A CN106353720 A CN 106353720A CN 201610800756 A CN201610800756 A CN 201610800756A CN 106353720 A CN106353720 A CN 106353720A
- Authority
- CN
- China
- Prior art keywords
- overbar
- rsqb
- lsqb
- beta
- alpha
- 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/0257—Hybrid positioning
- G01S5/0268—Hybrid positioning by deriving positions from different combinations of signals or of estimated positions in a single positioning system
-
- 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/16—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using electromagnetic waves other than radio waves
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明属于多传感器无源定位技术,针对固定辐射源目标,提供一种基于TDOA/GROA的多站连续定位模型。考虑利用中间变量的传统定位模型对定位场景的局限,从TDOA和GROA的量测方程出发,得到了只与目标状态向量有关的定位模型,从而避免了引入中间变量。然后根据量测模型推导了定位模型的误差项,最后利用约束加权最小二乘算法计算目标状态估计。该定位模型适用于多运动站场景下的固定辐射源目标连续定位。
Description
技术领域
本发明属于多传感器无源定位技术,涉及利用到达增益比等量测信息对辐射源目标的定位问题,提供了一种基于到达增益比的多站连续定位模型。
背景技术
无源定位技术一直是信号处理领域的重要研究课题,具有隐蔽性好、探测距离较远、适用性强等优点,在雷达,无线传感器网络等领域受到国内外学者的广泛关注。目前主要使用的量测信息有方位角、到达时差(TDOA)、到达频差以及到达增益比(GROA)等,或以上量测组合得到的多种联合定位系统。
传统的基于方位角的辐射源定位方法,需要在不同的角度对同一信号进行检测,并且容易因为坐标系不同和观测站姿态引入新的误差项。而基于到达时差和到达增益比量测的无源定位技术无需测量方位角、俯仰角等信息,因此不会受到坐标系转换和传感器载体姿态误差的影响,具有稳定性强、精度高等优点。目前的时差定位算法中,多通过引入额外的中间变量(辐射源到参考观测站的径向距离)建立线性时差方程。在此定位模型基础上,利用两步加权最小二乘(TS-WLS)算法、约束总体最小二乘(CTLS)算法、束加权最小二乘(CWLS)算法等进一步提高了定位精度。另一方面,由于辐射源信号会随着传播而产生能量衰减,将达到增益比量测信息加入到TDOA体制中,构造TDOA/GROA联合定位体制,并利用GROA量测对中间变量的修正,改善对辐射源的定位精度。
以上都是使用的基于中间变量的定位模型的定位算法,对于多运动站的情况下,由于观测站不断运动,中间变量即辐射源到参考观测站的径向距离也在不断变化,使得上述定位模型及其算法无法直接推广到多运动站连续定位,需要建立新的定位模型。
发明内容
本发明的目的在于针对固定辐射源目标,提供一种基于TDOA/GROA的多站连续定位模型,并利用CWLS算法求解辐射源位置信息。考虑到中间变量造成的局限性,从TDOA和GROA的量测方程出发,得到了只与目标状态向量有关的定位模型,从而避免了引入中间变量。然后根据量测模型推导了定位模型的误差项,得到了一种基于CWLS的TDOA/GROA联合定位算法。该算法适用于多运动站场景下的固定辐射源目标连续定位,其方法流程如图1所示,包括以下步骤:
1问题描述
假设定位场景中有M个运动观测站接收固定辐射源的信号,通过测量两个不同站之间的达到时间差和到达增益比来估计辐射源位置。N维空间时,由文献(15]可知得到闭式解的条件是M≥N+2。本文以2维空间为研究背景,所采用的方法可推广到3维空间中。记k时刻未知辐射源的位置向量为s0=[x0(k)y0(k)]T,第i个观测站位置向量为si=[xi(k)yi(k)]T(i=1,2,…,M)。以第1个观测站为参考站,则真实时差量测为
式中:为真实距离差量测,ri=|si-s0|(i=1,2,…,M)表示第i个观测站到辐射源目标的径向距离,c为电磁波传播速度。典型定位场景示意图如图2所示。
根据信号传播理论,信号的传播损耗与信号源和传感器之间距离的n次方成正比,在自由空间中,不考虑多径效应时,可将损耗因子n设定为常数2。同样以第1个观测站为参考站,则真实增益比量测为
2 TDOA/GROA定位模型
(1)TODA定位模型
基于TDOA量测方程,为表示方便,省略时间标签,该模型的向量形式为
对式(3)做变换,可得无测量误差时的TDOA线性方程为
式中:
实际到达距离差(RDOA)测量值中含有测量噪声,则实际RDOA量测模型为
式中:δgi1为噪声引起的RDOA测量误差,服从零均值的高斯分布。记到RDOA误差向量为δr(k)=[δr21(k),δr31(k),…,δrM1(k)]T,其量测噪声的协方差矩阵为
式中:T为观测时间,ω为信号带宽,ρ=Psignal/Pnoise为信号的信噪比,Psignal和Pnoise分别为信号和噪声的功率谱密度函数,IM-1为M-1阶单位矩阵,1(M-1)×(M-1)为M-1阶全1矩阵。
(2)GRDA定位模型
注意到对式(2)进行等式变换,省略时间标签可得
令i=2,上式可写作
式(9)和式(10)中均含有未知项-|s0|2,两两相减消去,整理后得到基于GROA定位模型的向量形式为
这样我们就得到了无测量误差时关于辐射源目标状态的GROA线性方程
式中:
实际到达增益比测量值中含有测量噪声δgi1(k),则实际GROA量测模型为
式中:δgi1为噪声引起的到达增益比测量误差,服从零均值的高斯分布。记到达增益比误差向量为δg(k)=[δg21(k),δg31(k),…,δgM1(k)]T,其协方差矩阵为
原有需要中间变量的GROA定位模型通过联合TDOA和GROA量测得到ri1=(gi1-1)r1,可知原定位模型通过加入GROA量测信息提高对中间变量r1的估计精度,进而提高算法性能,并没有通过GROA量测直接改善目标定位精度。由式(3)和式(11)可知所提无需中间变量的GROA定位模型能够直接估计目标状态,并避免引入中间变量,适用于多运动站的无源定位。
综上所述,通过引入GROA定位模型,构建TDOA/GROA联合定位模型,可以充分利用GROA中包含的目标位置信息,提高算法性能,并且保持了无需中间变量的优点,适用于多运动站对固定目标的连续定位。其联合定位模型为
式中:
3基于CWLS的定位算法
将式(6)和式(14)带入到和中,忽略高阶误差项后,整理可得
式中:A(k)和b(k)分别是带入实际测量值的观测矩阵和观测向量,且△A(k)=[-B1(k)v(k)-B2(k)v(k)],△b(k)=B3(k)v(k),忽略时间标签,误差项推导结果如式(19)所示。
式中:blkdaig(·)为分块对角矩阵。
将式(15)带入式(13),整理后可得
b(k)-A(k)x=△b(k)-△A(k)x=C(k)v(k) (20)
式中:C(k)=B1(k)x0+B2(k)y0+B3(k)
定义加权矩阵
W(k)=E(C(k)v(k)vT(k)CT(k))=C(k)RCT(k) (21)
式中协方差矩阵R=E(v(k)vT(k))=blkdiag(Rrr,Rgg)。
定义增广矩阵
Au=[-A b] (22)
和增广解向量
X=[x 1]T (23)
式中:A=[AT(1)AT(2)…AT(k)]T,b的定义与A类似。
易得
从而得到CWLS算法的最小化代价函数为
式中:ε的定义与A类似;W=blkdiag(W(1),W(2),…W(k))。
带入式(24)并对ε求期望,可得平均代价函数为
等式右边两项均为非负且与X相关。第一项为理想代价函数,当时,可取得最小值0。正是第二项的存在,当最小化E(ε)时,使解偏离真值从而产生偏差。因此,当约束第二项为常数时,E(ε)将在X取得真值时达到最小值。于是可构造CWLS问题为
式中:const为任意非负常数。利用拉格朗日乘子法可以解上述约束最小化问题。记λmin为 最小的广义特征值,λmin对应的广义特征值向量为,则辐射源的状态估计为
下面求解约束矩阵Ω,带入△Au=[B1v B2v B3v]整理后可得
式中:
注意到随着时间的推进,矩阵维数快速增大,导致计算量剧增,可利用式(31)和式(32)避免高维矩阵直接相乘,从而减小计算量,也可以对该算法加窗处理,提高算法的实时性。
注意到加权矩阵W与未知的辐射源的状态有关,计算时首先将W设为单位矩阵,得到辐射源的状态估计后,重新计算W,进而得到最终的目标状态估计。计算时,用实际测量值代替W和约束矩阵中的真实距离差和真实达到增益比。
附图说明
图1:技术方案流程图;
图2:典型定位场景。
具体实施方式
以下结合说明书附图对本发明作进一步详细描述。参照说明书附图1,本发明中基于TDOA/GROA的多站连续定位方法分为以下几个步骤:
步骤1:令k=1,初始化协方差矩阵R
R=blkdiag(Rrr,Rgg)
式中Rrr和Rgg分别为RDOA量测和GROA量测的协方差矩阵。
步骤2:根据TDOA和GORA量测信息以及站址信息计算定位模型Au(k)
Au(k)=[-A(k)b(k)]
其中ri1(k)和gi1(k)分别为k时刻的RODA和GROA的实际测量值且i=2,3,…,M,xi和yi为k时刻观测站位置且i=1,2,…,M;
省略时间标签k,误差项B1、B2和B3为
其中:blkdaig(·)为分块对角矩阵,IM-2为M-2阶单位矩阵,且i=1,2,…,M。
步骤3:若已经得到目标状态估计则加权矩阵W(k)为
若目标状态估计未知,令加权矩阵W(k)为单位阵;
步骤4:计算和Ω
并通过对其广义特征值分解,得到最小广义特征值对应的广义特征值向量;
步骤5:更新目标状态估计
令k=k+1,并从步骤2继续运行。
Claims (2)
1.基于TDOA/GROA的多站连续定位方法,其特征在于包括以下步骤:
步骤1:令k=1,初始化变量协方差矩阵R;
步骤2:根据TDOA和GORA量测信息以及站址信息计算定位模型Au(k)和误差项B1(k)、B2(k)和B3(k);
步骤3:若已经得到目标状态估计则计算加权矩阵W(k),若目标状态估计未知,令加权矩阵W(k)为单位阵;
步骤4:计算和Ω,并通过对其广义特征值分解,得到最小广义特征值对应的广义特征值向量
步骤5:更新目标状态估计令k=k+1,并从步骤2继续运行。
2.根据权利要求1所述的基于TDOA/GROA的多站连续定位方法,其特征在于步骤2采取以下技术措施:
初始化变量协方差矩阵R后,步骤2:根据TDOA和GORA量测信息以及站址信息计算定位模型Au(k)
Au(k)=[-A(k) b(k)]
其中:ri1(k)和gi1(k)分别为k时刻的RODA和GROA的实际测量值且i=2,3,…,M,xi和yi为k时刻观测站位置且i=1,2,…,M;
省略时间标签k,误差项B1、B2和B3为
其中:blkdaig(·)为分块对角矩阵,IM-2为M-2阶单位矩阵,且i=1,2,…,M。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610800756.3A CN106353720B (zh) | 2016-09-04 | 2016-09-04 | 基于tdoa/groa的多站连续定位模型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610800756.3A CN106353720B (zh) | 2016-09-04 | 2016-09-04 | 基于tdoa/groa的多站连续定位模型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106353720A true CN106353720A (zh) | 2017-01-25 |
CN106353720B CN106353720B (zh) | 2019-09-24 |
Family
ID=57858793
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610800756.3A Active CN106353720B (zh) | 2016-09-04 | 2016-09-04 | 基于tdoa/groa的多站连续定位模型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106353720B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108445445A (zh) * | 2018-02-06 | 2018-08-24 | 中国人民解放军战略支援部队信息工程大学 | 一种分布式无源定位方法及装置 |
CN108761442A (zh) * | 2018-06-01 | 2018-11-06 | 中国人民解放军战略支援部队信息工程大学 | 一种无源雷达目标定位方法及装置 |
CN108761387A (zh) * | 2018-06-01 | 2018-11-06 | 西安邮电大学 | 一种固定辐射源的双站时差频差联合定位方法 |
CN109212475A (zh) * | 2018-10-17 | 2019-01-15 | 哈尔滨工业大学 | 基于方位角和俯仰角信息的双机无源定位方法 |
CN109375159A (zh) * | 2018-10-10 | 2019-02-22 | 中国人民解放军空军预警学院雷达士官学校 | 纯方位加权约束总体最小二乘定位方法 |
CN113945888A (zh) * | 2021-10-19 | 2022-01-18 | 江南大学 | 基于tdoa的区间化无源定位方法及系统 |
CN115372902A (zh) * | 2022-08-05 | 2022-11-22 | 中国人民解放军战略支援部队信息工程大学 | 一种基于水下多基地声纳的tdoa偏置削减定位方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102026370A (zh) * | 2010-12-21 | 2011-04-20 | 西安电子科技大学 | 基于监测节点呈圆周分布的tdoa定位方法 |
CN103135094A (zh) * | 2013-01-31 | 2013-06-05 | 西安电子科技大学 | 基于bfgs拟牛顿法的信号源定位方法 |
CN105607092A (zh) * | 2016-01-27 | 2016-05-25 | 中国人民解放军国防科学技术大学 | 基于tdoa和功率测量值的gnss欺骗干扰定位方法 |
CN105807253A (zh) * | 2016-03-16 | 2016-07-27 | 中国人民解放军总参谋部工程兵科研三所 | 一种无须先验信息的多基站冗余信息法定位技术 |
-
2016
- 2016-09-04 CN CN201610800756.3A patent/CN106353720B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102026370A (zh) * | 2010-12-21 | 2011-04-20 | 西安电子科技大学 | 基于监测节点呈圆周分布的tdoa定位方法 |
CN103135094A (zh) * | 2013-01-31 | 2013-06-05 | 西安电子科技大学 | 基于bfgs拟牛顿法的信号源定位方法 |
CN105607092A (zh) * | 2016-01-27 | 2016-05-25 | 中国人民解放军国防科学技术大学 | 基于tdoa和功率测量值的gnss欺骗干扰定位方法 |
CN105807253A (zh) * | 2016-03-16 | 2016-07-27 | 中国人民解放军总参谋部工程兵科研三所 | 一种无须先验信息的多基站冗余信息法定位技术 |
Non-Patent Citations (3)
Title |
---|
徐征 等: "无需中间变量的多运动站时差定位新算法", 《航空学报》 * |
杜彦伸 等: "一种新的基于TDOA与GROA的无源定位算法", 《航空学报》 * |
郝本建 等: "基于TDOA与GROA的信号源被动定位偏差消除技术", 《电子学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108445445A (zh) * | 2018-02-06 | 2018-08-24 | 中国人民解放军战略支援部队信息工程大学 | 一种分布式无源定位方法及装置 |
CN108761442A (zh) * | 2018-06-01 | 2018-11-06 | 中国人民解放军战略支援部队信息工程大学 | 一种无源雷达目标定位方法及装置 |
CN108761387A (zh) * | 2018-06-01 | 2018-11-06 | 西安邮电大学 | 一种固定辐射源的双站时差频差联合定位方法 |
CN108761442B (zh) * | 2018-06-01 | 2021-06-11 | 中国人民解放军战略支援部队信息工程大学 | 一种无源雷达目标定位方法及装置 |
CN108761387B (zh) * | 2018-06-01 | 2022-03-08 | 西安邮电大学 | 一种固定辐射源的双站时差频差联合定位方法 |
CN109375159A (zh) * | 2018-10-10 | 2019-02-22 | 中国人民解放军空军预警学院雷达士官学校 | 纯方位加权约束总体最小二乘定位方法 |
CN109375159B (zh) * | 2018-10-10 | 2023-02-28 | 中国人民解放军空军预警学院雷达士官学校 | 纯方位加权约束总体最小二乘定位方法 |
CN109212475A (zh) * | 2018-10-17 | 2019-01-15 | 哈尔滨工业大学 | 基于方位角和俯仰角信息的双机无源定位方法 |
CN113945888A (zh) * | 2021-10-19 | 2022-01-18 | 江南大学 | 基于tdoa的区间化无源定位方法及系统 |
CN115372902A (zh) * | 2022-08-05 | 2022-11-22 | 中国人民解放军战略支援部队信息工程大学 | 一种基于水下多基地声纳的tdoa偏置削减定位方法 |
CN115372902B (zh) * | 2022-08-05 | 2023-12-01 | 中国人民解放军战略支援部队信息工程大学 | 一种基于水下多基地声纳的tdoa偏置削减定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106353720B (zh) | 2019-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106353720A (zh) | 基于tdoa/groa的多站连续定位模型 | |
CN103969622B (zh) | 一种基于多运动接收站的时差定位方法 | |
Poisel | Electronic warfare target location methods | |
CN106405533B (zh) | 基于约束加权最小二乘的雷达目标联合同步与定位方法 | |
CN110174643B (zh) | 一种无需噪声功率信息的基于到达时间差的定位方法 | |
CN107770859A (zh) | 一种考虑基站位置误差的tdoa‑aoa定位方法 | |
CN108872932B (zh) | 基于神经网络的超视距目标直接定位结果纠偏方法 | |
CN106597363A (zh) | 一种室内wlan环境下的行人定位方法 | |
CN102231912A (zh) | 一种基于rssi测距的室内无线传感器网络定位方法 | |
CN101561499B (zh) | 一种单站多普勒测距定位方法 | |
CN102981144A (zh) | 空中运动平台对目标的三维无源定位方法 | |
CN102186242A (zh) | 固定区域无线传感器网络移动节点定位方法 | |
CN107193023B (zh) | 一种具有闭式解的高精度北斗卫星系统单点定位方法 | |
CN106842118A (zh) | 用于分布式多点定位监视系统的时差定位方法 | |
CN111199280B (zh) | 短波信道模型误差存在下联合信号复包络和载波相位信息的多站目标源地理坐标估计方法 | |
CN109917333A (zh) | 融合aoa观测量与tdoa观测量的无源定位方法 | |
CN111551895A (zh) | 基于加权多维标度和拉格朗日乘子技术的运动辐射源tdoa和fdoa定位方法 | |
CN106842203B (zh) | 一种三维天线阵列综合孔径辐射计的图像反演方法 | |
CN110954865A (zh) | 一种基于电离层信息的短波时差定位方法 | |
CN107991676A (zh) | 星载单航过InSAR系统对流层误差校正方法 | |
CN104181513A (zh) | 一种雷达天线阵元位置的校正方法 | |
CN108279411B (zh) | 一种基于mds的被动mimo时差定位方法 | |
Cui et al. | Approximate closed-form TDOA-based estimator for acoustic direction finding via constrained optimization | |
CN109239686A (zh) | 一种针对分布式mimo雷达目标定位的发射机和接收机布局方法 | |
CN108445446B (zh) | 一种无源测速定位方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20190823 Address after: 264001 Research and Academic Department, 188 Erma Road, Zhifu District, Yantai City, Shandong Province Applicant after: Naval Aviation University of PLA Address before: 264001 Yantai City, Zhifu Province, No. two road, No. 188, Department of research, Applicant before: Naval Aeronautical Engineering Institute PLA |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |