CN106873013A - 一种基于约束凸加权的tdoa定位方法 - Google Patents
一种基于约束凸加权的tdoa定位方法 Download PDFInfo
- Publication number
- CN106873013A CN106873013A CN201710015046.4A CN201710015046A CN106873013A CN 106873013 A CN106873013 A CN 106873013A CN 201710015046 A CN201710015046 A CN 201710015046A CN 106873013 A CN106873013 A CN 106873013A
- Authority
- CN
- China
- Prior art keywords
- optimization problem
- satellite
- constraint
- tdoa
- target
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/421—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system
- G01S19/425—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system by combining or switching between signals derived from different satellite radio beacon positioning systems
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于约束凸加权的TDOA定位方法,TDOA定位是指利用目标辐射源信号到达不同观测站的时间差确定的多个双曲面相交进行定位。对于卫星定位系统,目标位于地球表面为一重要的先验信息。同时,利用TDOA测量信息进行定位,未知变量之间本身具有约束关系。为充分利用上述的约束信息以提高定位精度,本发明提供的方法首先建立关于目标位置的二次约束二次优化问题;而后针对多约束二次优化问题求解复杂度高,将原始的优化问题转化为两个单约束子优化问题进行求解,并对子优化问题的解以最优权值进行加权获得目标位置。本发明提供的一种基于约束凸加权的TDOA定位方法能以较低的计算复杂度获得较好的定位性能。
Description
技术领域
本发明涉及无线定位技术领域,特别涉及一种基于约束凸加权的TDOA定位方法。
背景技术
随着无线通信技术、互联网技术的发展,无线定位技术在紧急呼叫、船舶及车辆导航、声呐、智能交通等领域有着重要的应用。无线定位技术作为基于位置服务(LocationBased Service)的技术基础有着重要的研究价值。根据不同的测量参数,定位技术通常可分为:AOA定位技术、TOA定位技术、TDOA定位技术、FDOA定位技术以及利用两个测量参数的混合定位技术。在实际应用中,TOA定位技术和TDOA定位技术被证明往往能取得较好的定位精度。在卫星定位系统中,目标辐射源与卫星之间较小的时间偏差将会导致很大的定位误差,相较于TOA定位,TDOA定位只要求卫星之间同步而不要求目标辐射源与卫星之间同步,因而TDOA定位更具实用性。
对于星载TDOA定位,利用目标位于地球表面的先验信息,能够提高定位精度。由于TDOA定位方程的非线性度高,在考虑地球约束的地球面方程后,目标位置的求解更加困难。因此,研究定位精度较高且求解复杂度较低的TDOA定位算法具有重要意义。
发明内容
发明目的:在卫星定位系统中,为了充分利用未知变量之间的约束关系,以及目标位于地球表面的先验信息以提高定位精度,本发明提供的一种基于约束凸加权的TDOA定位方法能以较低的计算复杂度实现较好的定位精度。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种基于约束凸加权的TDOA定位方法,包括如下步骤:
(1)TDOA测量信息-距离差的转化:卫星定位系统有M个卫星,位置为ri=(xi,yi,zi)T,不失一般性,选择卫星1为参考卫星,目标辐射源信号到达卫星1与卫星i的时间差为τi1和Δτi1分别为信号到达卫星1与卫星i时间差的真实值与测量误差。假定信号的传播速度为光速c,由此得到相关的距离差为其中di1=||ri-r||-||r1-r||,r=(x,y,z)T为目标的位置;
(2)二次约束二次优化问题:将平方并整理得到
考虑到d1未知,定义变量w=(x,y,z,d1)T,将上式转化为矩阵形式b=Aw+n,其中A为(M-1)×4的矩阵,为(M-1)×1的向量;为(M-1)×1的向量,k=1,2,…,M-1。
目标位置r=(x,y,z)T与d1的约束关系,即变量约束为:以及目标位于地球表面的先验信息,地球约束为:rTr=re,re为地球半径。为得到表达形式简单的约束表达式,定义此时其中变量约束及地球约束为其中P=diag{1,1,1,-1},Q=diag{1,1,1,0},q=Qr1。二次约束二次优化问题表述为:
(3)子优化问题:上述的二次优化问题含有两个等式约束,求解复杂度高且难以获得全局最优解。为降低计算复杂度,将原优化问题转化为两个单约束子优化问题求解,得到和从而获得和
(4)目标位置估计值:充分利用两个约束条件以提高定位精度,通过最优加权因子α1和α2,对步骤(3)中的单约束子优化问题的解和进行加权,得到目标位置的最优估计值为其中α1+α2=1,且α1>0,α2>0。
有益效果:本发明提供的一种基于约束凸加权的TDOA定位方法,具有如下优点:1、本方法能充分利用变量之间的约束关系,及目标位于地球表面的先验信息;2、本方法求解复杂度低;3、本方法能获得较好的定位精度。
本发明附加的方面和优点将在下面的描述中部分给出,这些将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
图1示出了一种基于约束凸加权的TDOA定位方法流程图。
图2示出了采用五星定位系统,不同目标位置的定位精度随TDOA测量误差标准差的变化曲线。
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
考虑卫星定位系统有M个卫星,位置为ri=(xi,yi,zi)T,目标位置为r=(x,y,z)T,不失一般性,选择卫星1为参考卫星。为使本发明中的技术方案更加清楚,下面对本方案进行具体描述,算法流程图见图1:
S1.TDOA测量信息-距离差的转化:
目标辐射源信号到达卫星1与卫星i的时间差测量值为和Δτi1分别为信号到达卫星1与卫星i时间差的真实值与测量误差。假定信号的传播速度为光速c,由此得到相关的距离差为其中di1=||ri-r||-||r1-r||。
S2.二次约束二次优化问题:
对平方并整理得到
d1为目标到卫星1的距离。考虑到d1未知,定义变量w=(x,y,z,d1)T,将(1)式转化为矩阵形式b=Aw+n,其中A为(M-1)×4的矩阵, 为(M-1)×1的向量;的向量,k=1,2,…,M-1。
目标位置r=(x,y,z)T与d1的约束关系,即变量约束为:以及目标位于地球表面的先验信息,地球约束为:rTr=re,re为地球半径。为得到表达形式简单的约束表达式,定义此时其中变量约束及地球约束为其中P=diag{1,1,1,-1},Q=diag{1,1,1,0},q=Qr1。二次约束二次优化问题表述为:
S3.子优化问题:
上述(2)式的二次优化问题含有两个等式约束,求解复杂度高且难以获得全局最优解。为降低计算复杂度,将原优化问题转化为包含变量约束的子优化问题:
以及包含地球约束的子优化问题:
采用拉格朗日乘子法,通过(3)式子优化问题获得通过(4)式子优化问题获得将和分别代入(3)式和(4)式的约束条件得到关于λ1和λ2的多项式分别为求解多项式f(λ1)=0,得到λ11,...,λ1p,对于实值优化问题,去除其中的复数根,由剩余的λ1实数解得到由于采用如下准则选择最优的
同理得到最优的
S4.目标位置估计值:
充分利用两个约束条件以提高定位精度,通过最优加权因子α1和α2,对步骤S3中的单约束子优化问题的最优解和进行加权,得到目标位置的最优估计值为其中α1+α2=1,且α1>0,α2>0。
Claims (2)
1.一种基于约束凸加权的TDOA定位方法,其特征在于具体步骤包括:
(1)TDOA测量信息-距离差的转化:卫星定位系统有M个卫星,位置为ri=(xi,yi,zi)T,不失一般性,选择卫星1为参考卫星,目标辐射源信号到达卫星1与卫星i的时间差为τi1和Δτi1分别为信号到达卫星1与卫星i时间差的真实值与测量误差。假定信号的传播速度为光速c,由此得到相关的距离差为其中di1=||ri-r||-||r1-r||,r=(x,y,z)T为目标的位置;
(2)二次约束二次优化问题:将平方并整理得到
考虑到d1未知,定义变量w=(x,y,z,d1)T,将上式转化为矩阵形式b=Aw+n,其中A为(M-1)×4的矩阵,为(M-1)×1的向量;为(M-1)×1的向量,k=1,2,…,M-1。
目标位置r=(x,y,z)T与d1的约束关系,即变量约束为:以及目标位于地球表面的先验信息,地球约束为:rTr=re,re为地球半径。为得到表达形式简单的约束表达式,定义此时其中变量约束及地球约束为其中P=diag{1,1,1,-1},Q=diag{1,1,1,0},q=Qr1。二次约束二次优化问题表述为:
(3)子优化问题:上述的二次优化问题含有两个等式约束,求解复杂度高且难以获得全局最优解。为降低计算复杂度,将原优化问题转化为两个单约束的子优化问题求解,得到和从而获得和
(4)目标位置估计值:充分利用两个约束条件以提高定位精度,通过最优加权因子α1和α2,对步骤(3)中的单约束子优化问题的解和进行加权,得到目标位置的最优估计值为其中α1+α2=1,且α1>0,α2>0。
2.根据权利要求1所述的一种基于约束凸加权的TDOA定位方法,其特征在于:所述方法将原始的多约束二次优化问题转化为单约束子优化问题,并通过最优加权因子对子优化问题的解加权获得目标位置,既能充分利用约束条件,又能降低求解复杂度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710015046.4A CN106873013B (zh) | 2017-01-09 | 2017-01-09 | 一种基于约束凸加权的tdoa定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710015046.4A CN106873013B (zh) | 2017-01-09 | 2017-01-09 | 一种基于约束凸加权的tdoa定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106873013A true CN106873013A (zh) | 2017-06-20 |
CN106873013B CN106873013B (zh) | 2019-04-12 |
Family
ID=59165627
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710015046.4A Expired - Fee Related CN106873013B (zh) | 2017-01-09 | 2017-01-09 | 一种基于约束凸加权的tdoa定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106873013B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080231511A1 (en) * | 2007-03-22 | 2008-09-25 | Nortel Networks Limited | Beacon-Assisted Precision Location of Untethered Client in Packet Networks |
CN104283597A (zh) * | 2013-11-25 | 2015-01-14 | 北京邮电大学 | 一种波束赋形方法及设备 |
CN104600706A (zh) * | 2014-12-26 | 2015-05-06 | 国家电网公司 | 一种基于二次约束二次优化的含海上风电场电压无功优化控制方法 |
CN105306123A (zh) * | 2015-10-15 | 2016-02-03 | 哈尔滨工程大学 | 一种抗阵列系统误差的稳健波束形成方法 |
-
2017
- 2017-01-09 CN CN201710015046.4A patent/CN106873013B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080231511A1 (en) * | 2007-03-22 | 2008-09-25 | Nortel Networks Limited | Beacon-Assisted Precision Location of Untethered Client in Packet Networks |
CN104283597A (zh) * | 2013-11-25 | 2015-01-14 | 北京邮电大学 | 一种波束赋形方法及设备 |
CN104600706A (zh) * | 2014-12-26 | 2015-05-06 | 国家电网公司 | 一种基于二次约束二次优化的含海上风电场电压无功优化控制方法 |
CN105306123A (zh) * | 2015-10-15 | 2016-02-03 | 哈尔滨工程大学 | 一种抗阵列系统误差的稳健波束形成方法 |
Non-Patent Citations (1)
Title |
---|
朱伟强等: "基于TDOA_FDOA不相关估计的选择融合定位方法", 《宇航学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN106873013B (zh) | 2019-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105911518A (zh) | 机器人定位方法 | |
CN102928860B (zh) | 基于局部定位信息提高gps定位精度的方法 | |
CN102231912A (zh) | 一种基于rssi测距的室内无线传感器网络定位方法 | |
CN111157943B (zh) | 在异步网络中基于toa的传感器位置误差抑制方法 | |
Xu et al. | Deep learning based channel covariance matrix estimation with user location and scene images | |
CN102395197A (zh) | 一种基于残差加权的tdoa蜂窝定位方法 | |
CN105635963A (zh) | 多智能体分布式协同定位方法 | |
CN105044667A (zh) | 一种运动目标的双星跟踪方法、装置和系统 | |
CN108519578A (zh) | 一种基于群智感知的室内定位指纹库构建方法 | |
Wang et al. | An investigation and solution of angle based rigid body localization | |
Destino | Positioning in wireless networks: non-cooperative and cooperative algorithms | |
CN113239952A (zh) | 一种基于空间尺度注意力机制和矢量地图的航空图像地理定位方法 | |
CN103096465B (zh) | 一种环境自适应的多目标直接定位方法 | |
CN113313176A (zh) | 一种基于动态图卷积神经网络的点云分析方法 | |
Zhu et al. | Intelligent fingerprint-based localization scheme using CSI images for internet of things | |
Zhang et al. | Tips: Transformer based indoor positioning system using both csi and doa of wifi signal | |
CN105158730A (zh) | 基于mds子空间第4和第5特征向量的tdoa定位方法 | |
CN106873013A (zh) | 一种基于约束凸加权的tdoa定位方法 | |
CN105704676A (zh) | 利用信号时间相关性提高指纹室内定位精度的方法 | |
Shan et al. | Cooperative localization based on GPS correction and EKF in urban environment | |
Wisanmongkol et al. | An ensemble approach to deep‐learning‐based wireless indoor localization | |
CN113194401B (zh) | 一种基于生成式对抗网络的毫米波室内定位方法及系统 | |
Chen et al. | Applying geometric dilution of precision approximation to adaptive neural network learning for precise mobile station positioning | |
Nguyen et al. | Convolutional Neural Network-based UWB System Localization | |
CN113727433A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190412 Termination date: 20200109 |