CN109375165A - 一种tdoa多测量站的定位方法 - Google Patents

一种tdoa多测量站的定位方法 Download PDF

Info

Publication number
CN109375165A
CN109375165A CN201811399645.1A CN201811399645A CN109375165A CN 109375165 A CN109375165 A CN 109375165A CN 201811399645 A CN201811399645 A CN 201811399645A CN 109375165 A CN109375165 A CN 109375165A
Authority
CN
China
Prior art keywords
solution
algorithm
measuring
combination
fuzzy
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
Application number
CN201811399645.1A
Other languages
English (en)
Other versions
CN109375165B (zh
Inventor
胡正
郭利强
杨青
刘军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Electronics Technology Instruments Co Ltd CETI
Original Assignee
China Electronics Technology Instruments Co Ltd CETI
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Electronics Technology Instruments Co Ltd CETI filed Critical China Electronics Technology Instruments Co Ltd CETI
Priority to CN201811399645.1A priority Critical patent/CN109375165B/zh
Publication of CN109375165A publication Critical patent/CN109375165A/zh
Application granted granted Critical
Publication of CN109375165B publication Critical patent/CN109375165B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-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/08Position of single direction-finder fixed by determining direction of a plurality of spaced sources of known location

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)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提出了一种TDOA多测量站的定位方法,包括以下步骤:首先,选用三个测量站进行解算,为Taylor算法提供初值;当三个测量站解算有解时,判断该解是否位于测量站周围或者延长线周围,如果是,直接用LCLS算法进行计算;否则,将三个测量站解算结果作为Taylor算法的初始迭代值。本发明从初值的角度出发,设计了一种低复杂度的初值计算方法,不需要测量值的先验统计信息,可以为Taylor算法提供有效初值;选取一个好的初值,即初值与真实值之间的距离很近,不仅可以保证迭代的收敛,而且大大降低了迭代过程中对矩阵进行修正的次数,降低了计算复杂度和迭代次数。

Description

一种TDOA多测量站的定位方法
技术领域
本发明涉及通信技术领域,特别涉及一种TDOA多测量站的定位方法。
背景技术
TDOA定位是一种无源定位技术,是信号处理领域中的一个研究重点,广泛应用于雷达、声呐、无线通信等领域。在空间布局多台接收机,利用信号到达不同接收机的时间差,形成多个双曲线方程,通过对双曲线方程组的求解实现对目标信号的定位。对于该非线性方程组的求解存在非迭代和迭代算法。典型的非迭代算法有CHAN算法,LCLS(linear-correction least square)算法,迭代算法有Taylor算法,CTLS(Constrained TotalLeast Square)算法。
Taylor算法是一种基于泰勒级数展开的最小二乘估计迭代算法,基本适用于所有的定位系统,并利用测量参量来改善定位精度。首先在目标位置的初始估计点利用泰勒级数展开,忽略二次以上项,将非线性方程变为线性方程,采用最小二乘法进行估计;其次,利用估计的偏移量修正估计的目标位置,不断迭代,最终得到对目标位置的最优估计。
虽然Taylor算法在一定条件下能够得到方程组的精确解,但是,Taylor算法需要一个比较好的初始值,初始值必须靠近理想目标位置才能够保证算法的收敛。如果初始值没有设定好,那么Taylor级数方法可能不收敛,双曲线方程组无解。多数情况下,对双曲线方程组进行线性化不会给位置估计带来更多的误差。但是,GDOP(Geometric DilutionOfPrecision)较差或者噪声比较大条件下进行定位时,线性化操作会带来较大误差,可能导致最小二乘算法求伪逆值偏差极大,最终导致算法不收敛。
已有的算法采用CHAN算法为Taylor算法的提供初值,但是CHAN算法需要测量值的先验统计信息,计算量大,并且噪声测量误差较大时,该算法的性能会显著下降,这样会导致Taylor迭代不收敛。也有的算法对Taylor迭代过程中的迭代矩阵进行修正,如正则化Taylor修正算法,这种算法虽然可以改进Taylor算法收敛情况,但是计算量大,复杂度高。
发明内容
为解决上述现有技术中的不足,本发明提出了一种TDOA多测量站的定位方法,实现了低复杂度初始值的计算,同时,为了进一步保证收敛,采用正则化修正Taylor算法,保证算法的稳健收敛,提高算法的定位精度。
本发明的技术方案是这样实现的:
一种TDOA多测量站的定位方法,包括以下步骤:
首先,选用三个测量站进行解算,为Taylor算法提供初值;当三个测量站解算有解时,判断该解是否位于测量站周围或者延长线周围,如果是,直接用LCLS算法进行计算;否则,将三个测量站解算结果作为Taylor算法的初始迭代值。
可选地,所述方法还包括:如果迭代次数超过门限值或者计算结果无穷大,那么判定为迭代不收敛,采用正则化Taylor修正算法进行处理。
可选地,所述选用三个测量站进行解算,为Taylor算法提供初值的步骤,包括:
首先,选取三个测量站作为第一组合进行解算,获得三种解:1表示单解,2表示模糊解,0表示无解;
然后,再选取三个测量站作为第二组合进行解算,也会获得三种解:1表示单解,2表示模糊解,0表示无解;
两组测量站组合会出现9种解组合,一个单解和一个模糊解的解组合,或者两个模糊解的解组合找到距离最近两个解,取平均确定最后的解。
可选地,如果出现一个无解和一个模糊解的解组合,或者两个无解的解组合,则无法判定最后的解,则再取三个测量站作为第三组合进行解算,三组测量站的解分别是模糊解、无解、单解的解组合,或者是无解、模糊解、单解的解组合,或者是模糊解、模糊解、单解的解组合,能够找到距离最近两个解,取平均确定最后的解。
本发明的有益效果是:
(1)从初值的角度出发,设计了一种低复杂度的初值计算方法,不需要测量值的先验统计信息,可以为Taylor算法提供有效初值。
(2)选取一个好的初值,即初值与真实值之间的距离很近,不仅可以保证迭代的收敛,而且大大降低了迭代过程中对矩阵进行修正的次数,降低了计算复杂度和迭代次数。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明一种TDOA多测量站的定位方法的流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
通常Taylor算法的初值可以采用CNAN算法给出,但是CHAN算法需要测量值的先验统计信息,计算量大,并且噪声测量误差较大时,该算法的性能会显著下降,这样会导致Taylor迭代不收敛。本发明提出了一种TDOA多测量站的定位方法,不需要测量值的先验统计信息,可以为Taylor算法提供有效初值。
本发明提出了一种TDOA多测量站的定位方法,包括以下步骤:首先,直接选用三个测量站进行解算,为Taylor算法提供初值;当三个测量站解算有解时,判断该解是否位于测量站周围或者延长线周围,如果是,直接用LCLS算法进行计算;否则,三个测量站解算结果可以作为Taylor算法的初始迭代值。
当目标位于测量站周围或者测量站的连线的外延长线时,虽然三站组合可能给出相应的解,但是在这种情况下,三站组合给出的解一般误差很大。因此当三站解算有解时,判断该解是否位于测量站周围或者延长线周围,如果是这种情况,直接用LCLS算法进行计算;否则,三站解算结果可以作为Taylor算法的初始迭代值。
可选地,所述方法还包括:如果迭代次数超过门限值或者计算结果无穷大,那么判定为迭代不收敛,采用已有的正则化Taylor修正算法进行处理。
采用上述给出的初始值进行迭代,在测量噪声特别大的极端情况下,可能会出现不收敛的情况,在这种极端情况下,如果迭代次数超过门限值或者计算结果无穷大,那么判定为迭代不收敛,在这种情况下,采用已有的正则化修正Taylor方法进行处理,进一步保证迭代的收敛。
图1给出了TDOA多测量站的定位方法的一个可选实施例。
该可选实施例中,假设有N个测量站同时接收辐射源发射的信号,首先选取三个测量站作为第一组合进行解算,例如选取测量站1、测量站2、测量站3进行解算,获得三种解:1单解,2模糊解(双解),0无解。
然后,再选取三个测量站作为第二组合进行解算,例如选取测量站1、测量站2、测量站4进行解算,也会获得三种解:1单解,2模糊解(双解),0无解。
这样,两个组合会出现9种解的组合,例如,11组合表示第一组合与第二组合的解算结果都是单解的情况,12解组合表示第一组合是单解、第二组合是模糊解的组合,以此类推;一个单解和一个模糊解的组合,或者两个模糊解的组合可以找到距离最近两个解,例如12解组合、21解组合、22解组合可以通过找到距离最近两个解,取平均确定最后的解。
如果出现02、20、00这三种解组合,那么无法判定最后的解,则再取三个测量站作为第三组合进行解算,例如测量站1、测量站4、测量站5的组合。
第一组合与第二组合的解的20解组合与第三组合的单解的组合为201解组合,即三组三站的解分别是模糊解、无解、单解,可以通过找到距离最近两个解,取平均确定最后的解。
同理,如果当出现200,020,000这三种解组合时,无法为Taylor算法提供初始迭代值。
大部分情况下,两组三站组合基本可以给出初始迭代值,计算量低。出现计算三组组合给出初值的情况极少,如果三组组合依然不能给出解,那么直接选用LCLS算法计算。
本发明从初值的角度出发,设计了一种低复杂度的初值计算方法,不需要测量值的先验统计信息,可以为Taylor算法提供有效初值。选取一个好的初值,即初值与真实值之间的距离很近,不仅可以保证迭代的收敛,而且大大降低了迭代过程中对矩阵进行修正的次数,降低了计算复杂度和迭代次数。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种TDOA多测量站的定位方法,其特征在于,包括以下步骤:
首先,选用三个测量站进行解算,为Taylor算法提供初值;当三个测量站解算有解时,判断该解是否位于测量站周围或者延长线周围,如果是,直接用LCLS算法进行计算;否则,将三个测量站解算结果作为Taylor算法的初始迭代值。
2.如权利要求1所述的TDOA多测量站的定位方法,其特征在于,还包括:如果迭代次数超过门限值或者计算结果无穷大,那么判定为迭代不收敛,采用正则化Taylor修正算法进行处理。
3.如权利要求1所述的TDOA多测量站的定位方法,其特征在于,所述选用三个测量站进行解算,为Taylor算法提供初值的步骤,包括:
首先,选取三个测量站作为第一组合进行解算,获得三种解:1表示单解,2表示模糊解,0表示无解;
然后,再选取三个测量站作为第二组合进行解算,也会获得三种解:1表示单解,2表示模糊解,0表示无解;
两组测量站组合会出现9种解组合,一个单解和一个模糊解的解组合,或者两个模糊解的解组合找到距离最近两个解,取平均确定最后的解。
4.如权利要求3所述的TDOA多测量站的定位方法,其特征在于,
如果出现一个无解和一个模糊解的解组合,或者两个无解的解组合,则无法判定最后的解,则再取三个测量站作为第三组合进行解算,三组测量站的解分别是模糊解、无解、单解的解组合,或者是无解、模糊解、单解的解组合,或者是模糊解、模糊解、单解的解组合,能够找到距离最近两个解,取平均确定最后的解。
CN201811399645.1A 2018-11-22 2018-11-22 一种tdoa多测量站的定位方法 Active CN109375165B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811399645.1A CN109375165B (zh) 2018-11-22 2018-11-22 一种tdoa多测量站的定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811399645.1A CN109375165B (zh) 2018-11-22 2018-11-22 一种tdoa多测量站的定位方法

Publications (2)

Publication Number Publication Date
CN109375165A true CN109375165A (zh) 2019-02-22
CN109375165B CN109375165B (zh) 2021-12-14

Family

ID=65376996

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811399645.1A Active CN109375165B (zh) 2018-11-22 2018-11-22 一种tdoa多测量站的定位方法

Country Status (1)

Country Link
CN (1) CN109375165B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111967197A (zh) * 2020-08-27 2020-11-20 中电科仪器仪表有限公司 一种基于圆概率误差的多测量站布局方法及系统
JP2023506803A (ja) * 2020-04-10 2023-02-20 中興通訊股▲ふん▼有限公司 協調測位方法、装置、機器及び記憶媒体

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7667640B2 (en) * 2007-04-13 2010-02-23 Glowlink Communications Technology, Inc. Determining a geolocation solution of an emitter on earth using satellite signals
CN101813763A (zh) * 2009-12-22 2010-08-25 民航数据通信有限责任公司 多点定位方法、系统及中心站
CN103033825A (zh) * 2012-12-26 2013-04-10 江苏科技大学 一种gnss接收机的定位解算方法
CN104113911A (zh) * 2014-05-14 2014-10-22 浙江工业大学 基于mle和ukf组合的wsn节点定位方法
CN106443655A (zh) * 2016-09-21 2017-02-22 河海大学 一种多输入多输出雷达近场定位算法
CN106483496A (zh) * 2016-10-19 2017-03-08 河南城建学院 基于chan算法与改进牛顿迭代的联合时差定位方法
CN107229033A (zh) * 2017-05-16 2017-10-03 电子科技大学 基于高度维分段搜索的多目标到达时间差定位方法
CN107300687A (zh) * 2017-03-22 2017-10-27 哈尔滨工程大学 一种基于运动多站的高精度无源时差定位方法
CN108738130A (zh) * 2017-04-24 2018-11-02 温州市鹿城区中津先进科技研究院 一种基于泰勒级数展开的蜂窝tdoa定位算法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7667640B2 (en) * 2007-04-13 2010-02-23 Glowlink Communications Technology, Inc. Determining a geolocation solution of an emitter on earth using satellite signals
CN101813763A (zh) * 2009-12-22 2010-08-25 民航数据通信有限责任公司 多点定位方法、系统及中心站
CN103033825A (zh) * 2012-12-26 2013-04-10 江苏科技大学 一种gnss接收机的定位解算方法
CN104113911A (zh) * 2014-05-14 2014-10-22 浙江工业大学 基于mle和ukf组合的wsn节点定位方法
CN106443655A (zh) * 2016-09-21 2017-02-22 河海大学 一种多输入多输出雷达近场定位算法
CN106483496A (zh) * 2016-10-19 2017-03-08 河南城建学院 基于chan算法与改进牛顿迭代的联合时差定位方法
CN107300687A (zh) * 2017-03-22 2017-10-27 哈尔滨工程大学 一种基于运动多站的高精度无源时差定位方法
CN108738130A (zh) * 2017-04-24 2018-11-02 温州市鹿城区中津先进科技研究院 一种基于泰勒级数展开的蜂窝tdoa定位算法
CN107229033A (zh) * 2017-05-16 2017-10-03 电子科技大学 基于高度维分段搜索的多目标到达时间差定位方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
房嘉奇 等: "TDOA中的修正牛顿及泰勒级数方法", 《西安电子科技大学学报(自然科学版)》 *
谢红 等: "基于TDOA的一种简化的非视距误差抑制算法", 《哈尔滨工程大学学报》 *
邓勇 等: "平面三站时差定位中的模糊及无解研究", 《系统工程与电子技术》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2023506803A (ja) * 2020-04-10 2023-02-20 中興通訊股▲ふん▼有限公司 協調測位方法、装置、機器及び記憶媒体
JP7344389B2 (ja) 2020-04-10 2023-09-13 中興通訊股▲ふん▼有限公司 協調測位方法、装置、機器及び記憶媒体
CN111967197A (zh) * 2020-08-27 2020-11-20 中电科仪器仪表有限公司 一种基于圆概率误差的多测量站布局方法及系统

Also Published As

Publication number Publication date
CN109375165B (zh) 2021-12-14

Similar Documents

Publication Publication Date Title
US8331953B2 (en) System and method for estimating the location of a mobile device
US8280384B2 (en) System and method for predicting location accuracy of a mobile unit
Yedavalli et al. Sequence-based localization in wireless sensor networks
CN100449330C (zh) 位置计算方法
CN104507050A (zh) 一种WiFi室内定位中概率型指纹匹配方法
US9660740B2 (en) Signal strength distribution establishing method and wireless positioning system
CN105792129A (zh) 一种非视距环境下的定位方法及装置
CN103376447A (zh) 一种实现非合作双多基地雷达目标三维定位的方法
CN109031190B (zh) 一种高重频脉冲信号无源时差定位方法
CN109375165A (zh) 一种tdoa多测量站的定位方法
CN110636436A (zh) 基于改进chan算法的三维uwb室内定位方法
CN111273266A (zh) 多站雷达定位系统
CN115598593A (zh) 一种等长短基线高精度测向定位方法、系统、设备及终端
CN107371133B (zh) 一种提高基站定位精度的方法
CN107613458A (zh) 一种tdoa条件下最优联合时间同步与定位的定位方法
CN107592654B (zh) 一种基于压缩感知的同频多辐射源场强定位方法
CN105807253A (zh) 一种无须先验信息的多基站冗余信息法定位技术
CN111487586A (zh) 基于分布式无源定位技术的定位精度提高方法
KR101509569B1 (ko) 지구곡률과 신호방향정보를 이용한 신호원 위치추정방법
CN109709555B (zh) 一种识别天气雷达相邻体扫数据差异性的方法及系统
Zhang et al. DOA‐Based Localization Method with Multiple Screening K‐Means Clustering for Multiple Sources
Li et al. Mobile station location estimation for MIMO communication systems
Yan et al. Robust Target Localization for Multistatic Passive Radar Networks
CN107786939B (zh) 基于蒙特卡洛最小二乘法室内定位模型及构建方法和应用
CN109991571B (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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province

Applicant after: CLP kesiyi Technology Co.,Ltd.

Address before: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province

Applicant before: CHINA ELECTRONICS TECHNOLOGY INSTRUMENTS Co.,Ltd.

GR01 Patent grant
GR01 Patent grant