CN109375165B - 一种tdoa多测量站的定位方法 - Google Patents
一种tdoa多测量站的定位方法 Download PDFInfo
- Publication number
- CN109375165B CN109375165B CN201811399645.1A CN201811399645A CN109375165B CN 109375165 B CN109375165 B CN 109375165B CN 201811399645 A CN201811399645 A CN 201811399645A CN 109375165 B CN109375165 B CN 109375165B
- Authority
- CN
- China
- Prior art keywords
- solution
- solutions
- combination
- algorithm
- 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.)
- 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/08—Position 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)
- Radar Systems Or Details Thereof (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提出了一种TDOA多测量站的定位方法,包括以下步骤:首先,选用三个测量站进行解算,为Taylor算法提供初值;当三个测量站解算有解时,判断该解是否位于测量站周围或者延长线周围,如果是,直接用LCLS算法进行计算;否则,将三个测量站解算结果作为Taylor算法的初始迭代值。本发明从初值的角度出发,设计了一种低复杂度的初值计算方法,不需要测量值的先验统计信息,可以为Taylor算法提供有效初值;选取一个好的初值,即初值与真实值之间的距离很近,不仅可以保证迭代的收敛,而且大大降低了迭代过程中对矩阵进行修正的次数,降低了计算复杂度和迭代次数。
Description
技术领域
本发明涉及通信技术领域,特别涉及一种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 (1)
1.一种TDOA多测量站的定位方法,其特征在于,包括以下步骤:
首先,选用三个测量站进行解算,为Taylor算法提供初值;当三个测量站解算有解时,判断该解是否位于测量站周围或者延长线周围,如果是,直接用LCLS算法进行计算;否则,将三个测量站解算结果作为Taylor算法的初始迭代值;
所述选用三个测量站进行解算,为Taylor算法提供初值的步骤,包括:
首先,选取三个测量站作为第一组合进行解算,获得三种解:1表示单解, 2表示模糊解,0表示无解;
然后,再选取三个测量站作为第二组合进行解算,也会获得三种解:1表示单解, 2表示模糊解,0表示无解;
两组测量站组合会出现9种解组合,一个单解和一个模糊解的解组合,或者两个模糊解的解组合找到距离最近两个解,取平均确定最后的解;
如果出现一个无解和一个模糊解的解组合,或者两个模糊解的解组合,则无法判定最后的解,则再取三个测量站作为第三组合进行解算,三组测量站的解分别是模糊解、无解、单解的解组合,或者是无解、模糊解、单解的解组合,或者是模糊解、模糊解、单解的解组合,能够找到距离最近两个解,取平均确定最后的解;还包括:如果迭代次数超过门限值或者计算结果无穷大,那么判定为迭代不收敛,采用正则化Taylor修正算法进行处理。
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 CN109375165A (zh) | 2019-02-22 |
CN109375165B true 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) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111896914A (zh) * | 2020-04-10 | 2020-11-06 | 中兴通讯股份有限公司 | 一种协同定位方法、装置、设备和存储介质 |
CN111967197A (zh) * | 2020-08-27 | 2020-11-20 | 中电科仪器仪表有限公司 | 一种基于圆概率误差的多测量站布局方法及系统 |
Citations (6)
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 |
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算法与改进牛顿迭代的联合时差定位方法 |
CN108738130A (zh) * | 2017-04-24 | 2018-11-02 | 温州市鹿城区中津先进科技研究院 | 一种基于泰勒级数展开的蜂窝tdoa定位算法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813763B (zh) * | 2009-12-22 | 2012-03-07 | 民航数据通信有限责任公司 | 多点定位方法、系统及中心站 |
CN107300687B (zh) * | 2017-03-22 | 2020-06-16 | 哈尔滨工程大学 | 一种基于运动多站的高精度无源时差定位方法 |
CN107229033B (zh) * | 2017-05-16 | 2020-04-10 | 电子科技大学 | 基于高度维分段搜索的多目标到达时间差定位方法 |
-
2018
- 2018-11-22 CN CN201811399645.1A patent/CN109375165B/zh active Active
Patent Citations (6)
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 |
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算法与改进牛顿迭代的联合时差定位方法 |
CN108738130A (zh) * | 2017-04-24 | 2018-11-02 | 温州市鹿城区中津先进科技研究院 | 一种基于泰勒级数展开的蜂窝tdoa定位算法 |
Non-Patent Citations (3)
Title |
---|
TDOA中的修正牛顿及泰勒级数方法;房嘉奇 等;《西安电子科技大学学报(自然科学版)》;20161231;第43卷(第6期);第28页 * |
基于TDOA的一种简化的非视距误差抑制算法;谢红 等;《哈尔滨工程大学学报》;20050228;第26卷(第1期);文章第一部分一种NLOS抑制算法 * |
平面三站时差定位中的模糊及无解研究;邓勇 等;《系统工程与电子技术》;20001231;第22卷(第3期);文章第三部分模糊研究、第四部分无解研究 * |
Also Published As
Publication number | Publication date |
---|---|
CN109375165A (zh) | 2019-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101260647B1 (ko) | 무선센서네트워크 상에서 효율적인 다변측량을 이용한 무선측위방법과 이를 실시하기 위한 프로그램이 기록된 기록매체 | |
Yedavalli et al. | Sequence-based localization in wireless sensor networks | |
JP2004361412A (ja) | 信頼領域フィルタリング法を用いた位置推定方法および装置 | |
CN109375165B (zh) | 一种tdoa多测量站的定位方法 | |
CN107820206B (zh) | 基于信号强度的非视距定位方法 | |
CN108882149B (zh) | 距离相关概率的nlos补偿定位方法 | |
JP2010539483A (ja) | Los条件及びnlos条件の下での地図情報を使用した移動端末装置の位置の向上した線形最小二乗推定のための方法 | |
CN111123197A (zh) | 基于tdoa的目标辐射源定位方法 | |
CN110806561B (zh) | 一种多基站的自标定方法 | |
CN108414974B (zh) | 一种基于测距误差矫正的室内定位方法 | |
KR20080078576A (ko) | 순방향 링크 aoa를 이용한 위치 추정 장치 및 방법 | |
KR101865120B1 (ko) | 테스트 노드 기반의 무선 측위 방법 및 그 장치 | |
CN106575825B (zh) | 用于使用模态天线寻找信号方向的方法 | |
Wu et al. | Effective error control of iterative localization for wireless sensor networks | |
CN111487586A (zh) | 基于分布式无源定位技术的定位精度提高方法 | |
CN112083375A (zh) | 一种基于位置指纹和Taylor的协同定位算法 | |
US11979848B2 (en) | Access point based location system for high density WiFi deployments | |
CN108848447B (zh) | 一种采用未知节点修正的差分DV_Distance节点定位方法 | |
CN110856101B (zh) | 一种基于曲线拟合的无线传感器网络节点定位方法 | |
CN114205749B (zh) | 一种适用于穿墙场景的超宽带迭代定位算法及装置 | |
KR101509569B1 (ko) | 지구곡률과 신호방향정보를 이용한 신호원 위치추정방법 | |
CN109282816B (zh) | 一种利用三角形约束的测距定位误差抑制方法 | |
KR101750502B1 (ko) | 센서노드의 자세 보정 및 측위 방법 | |
CN107786939B (zh) | 基于蒙特卡洛最小二乘法室内定位模型及构建方法和应用 | |
Leca et al. | Determining optimum base stations configuration for TOA localization inside celullar networks |
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 |
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. |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |