CN113727433A - 一种未知信号传播速度下远近场到达时间差测量定位方法 - Google Patents
一种未知信号传播速度下远近场到达时间差测量定位方法 Download PDFInfo
- Publication number
- CN113727433A CN113727433A CN202111288597.0A CN202111288597A CN113727433A CN 113727433 A CN113727433 A CN 113727433A CN 202111288597 A CN202111288597 A CN 202111288597A CN 113727433 A CN113727433 A CN 113727433A
- Authority
- CN
- China
- Prior art keywords
- time difference
- far
- parameter
- arrival
- field
- 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.)
- Pending
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
- H04W64/006—Locating users or terminals or network equipment for network management purposes, e.g. mobility management with additional information processing, e.g. for direction or speed determination
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提出了一种未知信号传播速度下远近场到达时间差测量定位方法,通过设计近场信号时间达到时间差测量模型和远场信号达到时间差测量模型,然后设计估计参数g,通过估计参数g对近场信号时间达到时间差测量模型和远场信号达到时间差测量模型进行关系变换,得到近场信号时间达到时间差测量和远场信号达到时间差测量通用的远近场定位估计模型。本发明通过以上设置设计了一种将近场和远场统一解形式的定位方法,解决了在未知信号的传播速度的情况下的远场和近场的方法不可互相通用的问题。
Description
技术领域
本发明属于无线通讯技术领域,具体地说,涉及一种未知信号传播速度下远近场到达时间差测量定位方法。
背景技术
源目标定位问题在雷达、声纳、导航搜索与营救、无线网络节点的位置获取等领域中有着广泛的应用。为实现源目标定位,通常的做法是利用各种距离测量相关技术,包括到达时间、到达时间差、到达角度与信号接收强度等,采用已知位置坐标的信标节点去推算源节点目标的位置坐标。到达时间差测量方法通过时间差测量方法消除了发送端的时钟误差,在近场和远场情况下的源目标位置参数获取中都有大量的应用。
TDOA定位是一种利用时间差进行定位的方法。通过测量信号到达监测站的时间,可以确定信号源的距离。利用信号源到各个监测站的距离(以监测站为中心,距离为半径作圆),就能确定信号的位置。但是绝对时间一般比较难测量,通过比较信号到达各个监测站的绝对时间差,就能作出以监测站为焦点,距离差为长轴的双曲线,双曲线的交点就是信号的位置。
TDOA是一种无线定位技术。不同于TOA,TDOA(到达时间差)是通过检测信号到达两个基站的绝对时间差,而不是到达的飞行时间来确定移动台的位置,降低了信号源与各个监测站的时间同步要求,但提高了各个监测站的时间同步要求。采用三个不同的基站可以测到两个TDOA,移动站位于两个TDOA决定的双曲线的交点上。
该定位技术可应用于各种移动通信系统,尤其适用于CDMA系统,CDMA系统用扩频方式将信号频谱扩展到很宽的范围,使系统具有较强的抗多径能力。CDMA属"非功率敏感系统",信号衰减对时间测量的精度影响较小。TDOA法与TOA法比较优点之一是:当计算TDOA值时,计算误差对所有的基站是相同的且其和为零,这些误差包括公共的多径时延和同步误差。但由于功率控制造成离服务基站近的移动台发射功率小,使得与服务基站相邻的参与定位的另一基站接收到的功率非常小(即相邻基站的SNR太小),造成比较大的测量误差。
TDOA算法是对TOA算法的改进,他不是直接利用信号到达时间,而是用多个基站接收到信号的时间差来确定移动台位置,与TOA算法相比他不需要加入专门的时间戳,定位精度也有所提高。TDOA值的获取一般有2种形式:
第1种形式是利用移动台到达2个基站的时间TOA,取其差值来获得,这时仍需要基站时间的严格同步,但是当两基站间移动信道传输特性相似时,可减少由多径效应带来的误差。
第2种形式是将一个移动台接收到的信号与另一个移动台接收到的信号进行相关运算,从而得到TDOA的值,这种算法可以在基站和移动台不同步时,估计出TDOA的值,由于实际应用中,往往很难做到基站与移动台的同步,所以利用相关估计得到TDOA值,再进行定位计算能获得较高精度。对于蜂窝网中的移动台定位而言,TDOA更具有实际意义,这种方法对网络的要求相对较低,并且定位精度较高,已经成为研究的热点。
通过建立TDOA测量方程的近场优化模型,并对模型进行近似等效、线性化处理,并采用各种算法包括极大似然、凸优化算法及最小二乘代数法等进行计算是近场TDOA定位常采用的方法。采用数值计算方法的ML估计方法依赖于初始解,有可能陷入局部最优。为此,常采用凸优化松弛法及最小二乘代数法进行求解。常用的凸优化方法有半正定规划( SDP,Semidefinite Programming)、二次锥规划等方法,但凸优化算法的计算复杂度较高。为降低计算复杂度,也可将近场TDOA测量模型作近似线性化处理,设计源目标位置估计的线性最小二乘代数法。当源目标位置远离信标节点位置坐标时,采用近场定位方法难以准确估计。为此在远场情况下,方向角( DOA,Direction of Angle) 估计目前远场估计中常采用的方法[12-13]。在阵列信号处理领域中,采用相位差进行 DOA 估计的原理类似于TDOA。经典的DOA估计方法包括多元信号分类算法( MUSIC,Multiple Signal Classification)和基于旋转不变性原理的信号参数估算技术(ESPRIT,Esti-mation Signal Parametervia Rotational Invariant Tech- nique)对信号参数进行估计。远场DOA估计与近场定位算法分别实现了远近场两种不同情况下的参数估计问题,需要根据源目标与信标节点的距离远近来选择合适的估计算法。并且TDOA的到达时间差测量方法难以直接确定目前是处于近场或远场,而设计远近场统一解形式的定位算法可以有效解决以上问题。
发明内容
本发明针对现有技术的上述需求,提出了一种未知信号传播速度下远近场到达时间差测量定位方法,提供设计了一种将近场和远场统一解形式的定位方法来解决远场和近场的方法不可互相通用的问题。
本发明具体实现内容如下:
本发明提出了一种未知信号传播速度下远近场到达时间差测量定位方法,用于对源目标进行定位,包括以下步骤:
步骤1:建立三维和二维的下的位置坐标模型;
步骤2:建立近场信号到达时间差测量模型;
步骤3:建立远场信号到达时间差测量模型;
步骤4:设计待估计参数g,基于估计参数g根据得到的近场信号时间达到时间差测量模型和远场信号达到时间差测量模型计算得到三维和二维下的远近场定位估计模型;
步骤5:计算待估计参数g、信号速度c和源节点位置u的位置坐标参数;
步骤6:使用计算得到待估计参数g、信号速度c和源节点位置u的位置坐标参数,结合得到的远近场定位估计模型进行到达时间差测量定位计算。
为了更好地实现本发明,进一步地,所述步骤1的具体操作为:
步骤1.1:定义已知位置坐标的M+1个信标节点为Si,i=0,1,2,3,...,M;并定义待确定坐标的源节点位置为u;
步骤1.2:定义二维下信标节点Si的位置坐标Si=[xi , yi]T,源节点位置坐标u=[x , y]T;定义三维下信标节点Si的位置坐标Si=[xi , yi ,zi]T,源节点位置坐标u=[x ,y , z]T;
将位置坐标分别转换为极坐标形式:
二维空间下:Si=[ρi , αi]T,u=[ρ , θ]T;
三维空间下:Si=[ρi , αi ,βi]T,u=[ρ , θ , Φ]T;
为了更好地实现本发明,进一步地,所述步骤2中的近场信号到达时间差测量模型中,对于到达时间差ti采用距离差di和信号速度c的比值进行计算获得。
为了更好地实现本发明,进一步地,所述步骤3的具体操作为:
在三维空间下,需要先计算信标节点Si的位置参数以及源节点u的位置参数;
所述位置参数为信标节点Si的位置坐标除以信标节点Si的极径,所述位置参数为源节
点u的位置坐标除以源节点的极径;位置参数和位置参数中都包括三个元素,分别为:
各自对应的方向角的余弦值和仰角的余弦值的和,各自对应的方向角的正弦值和仰角的正
弦值的和,各自的仰角的正弦值;然后通过位置参数和位置参数进行反余弦处理得到
三维空间下的角度差参数。
为了更好地实现本发明,进一步地,所述步骤4的具体操作为:
为了更好地实现本发明,进一步地,所述步骤4具体操作为:
步骤4.1:将近场信号到达时间差测量模型的等式两边进行平方处理,得到如下关系式:
式中,c为信号传播速度,ti为源节点u分别与信标节点si之间的信号到达时间差,c为信号传播速度,ni为时间测量噪声;
步骤4.3:根据估计参数方法模型,结合近场信号时间达到时间差测量模型和远场信号达到时间差测量模型得到如下远近场定位估计模型:
采用远近场定位估计模型对未知信号传播速度下的远场和近场信号进行定位。
为了更好地实现本发明,进一步地,在所述步骤5中对信号传播速度c、仰角β、方向角α和待估计参数g进行计算,具体步骤为:
步骤5.2:设置如下最小化问题约束条件方程:
步骤5.4:对步骤5.2的表达式的等式两边取平方,忽略二次高阶误差项,得到如下表达式:
式中,v=1,...,p;
然后将v=p代入真实值的表达式中,并结合步骤5.3中的误差公式,得到:
采用微分方法展开后,忽略二次高阶误差项,得到如下表达式:
式中,1表示全1矩阵,下标表示全1矩阵的阶数;0表示全0矩阵,下标表示全0矩阵的阶数;
步骤5.8:计算得到在二维平面的情况下:
在三维空间的情况下:
本发明与现有技术相比具有如下优点及有益效果:
本发明设计的方法可以实现远场和近场的定位统一,且通过约束条件使得算法精度保持一定高度。
附图说明
图1为近场示意图;
图2为远场示意图。
具体实施方式
为了更清楚地说明本发明实施例的技术方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,应当理解,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例,因此不应被看作是对保护范围的限定。基于本发明中的实施例,本领域普通技术工作人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“设置”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;也可以是直接相连,也可以是通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
实施例1:
本实施例提出了一种未知信号传播速度下远近场到达时间差测量定位方法,用于对源目标进行定位,包括以下步骤:
步骤1:建立三维和二维的下的位置坐标模型;
步骤2:建立近场信号到达时间差测量模型;
步骤3:建立远场信号到达时间差测量模型;
步骤4:设计待估计参数g,基于估计参数g根据得到的近场信号时间达到时间差测量模型和远场信号达到时间差测量模型计算得到三维和二维下的远近场定位估计模型;
步骤5:计算待估计参数g、信号速度c和源节点位置u的位置坐标参数;
步骤6:使用计算得到待估计参数g、信号速度c和源节点位置u的位置坐标参数,结合得到的远近场定位估计模型进行到达时间差测量定位计算。
工作原理:通过将远场和近场进行统一,实现了获取到共同的参数即可进行快速定位计算,实现了快速简便且精准的定位计算。
实施例2:
本实施例在上述实施例1的基础上,如图1、图2所示,为了更好地实现本发明,进一步地,所述步骤1的具体操作为:
步骤1.1:定义已知位置坐标的M+1个信标节点为Si,i=0,1,2,3,...,M;并定义待确定坐标的源节点位置为u;
步骤1.2:定义二维下信标节点Si的位置坐标Si=[xi , yi]T,源节点位置坐标u=[x , y]T;定义三维下信标节点Si的位置坐标Si=[xi , yi ,zi]T,源节点位置坐标u=[x ,y , z]T;
将位置坐标分别转换为极坐标形式:
二维空间下:Si=[ρi , αi]T,u=[ρ , θ]T;
三维空间下:Si=[ρi , αi ,βi]T,u=[ρ , θ , Φ]T;
所述步骤2中的近场信号到达时间差测量模型中,对于到达时间差ti采用距离差di和信号速度c的比值进行计算获得,具体如下:
式中:di为距离差,c为信号传播速度,ti为源节点u分别与信标节点si之间的信号到达时间差,ni为时间测量噪声。
所述步骤3的具体操作为:
所述步骤4的具体操作为:
步骤4.1:定义待估计参数,通过数学关系换算,将得到的近场信号时间
达到时间差测量模型和远场信号达到时间差测量模型进行结合,得到以位置参数进行表
达的远近场定位估计模型,通过位置参数换算得到到达时间差信息。
工作原理:TDOA 定位问题模型是通过源节点与信标节点间的 TDOA 测量及已知位置坐标的信标节点si,以此估计源节点位置坐标u。在步骤2中,可将TDOA测量方程表示为:
式中: ti 表示了源节点分别与信标节点si、s0之间的信号到达时间差( TDOA),c为信号传播速度,则距离差di = cti,ni 为时间测量噪声,可设ni服从均值为零,方差为δ2i的高斯噪声。不失一般性,可假设s0=[0 0]T,故式子可重写为:
显然,‖u‖ = ρ。
在步骤3中,建立远场信号到达时间差测量模型,具体模型建立如下:
式中:di为距离差,
工作原理:由于源节点距离信标节点si 与s0 都相当远,若不考虑噪声,通过推导有远场测量近似方程:
式中:di为距离差,
所述步骤4具体操作为:
步骤4.1:将近场信号到达时间差测量模型的等式两边进行平方处理,得到如下关系式:
式中,c为信号传播速度,ti为源节点u分别与信标节点si之间的信号到达时间差,c为信号传播速度,ni为时间测量噪声;
步骤4.3:根据估计参数方法模型,结合近场信号时间达到时间差测量模型和远场信号达到时间差测量模型得到如下远近场定位估计模型:
采用远近场定位估计模型对未知信号传播速度下的远场和近场信号进行定位。
当源节点处于远场情况时,g接近于零,不难看出,不考虑噪声时所有的公式都是一致的,故本方法既适用于近场 TDOA 定位问题,也适用于远场估计问题。
本实施例的其他部分与上述实施例1相同,故不再赘述。
实施例3:
在三维空间下,需要先计算信标节点Si的位置参数以及源节点u的位置参数;
所述位置参数为信标节点Si的位置坐标除以信标节点Si的极径,所述位置参数为源节
点u的位置坐标除以源节点的极径;位置参数和位置参数中都包括三个元素,分别为:
各自对应的方向角的余弦值和仰角的余弦值的和,各自对应的方向角的正弦值和仰角的正
弦值的和,各自的仰角的正弦值;然后通过位置参数和位置参数进行反余弦处理得到
三维空间下的角度差参数。
本实施例的其他部分与上述实施例1-2任一项相同,故不再赘述。
实施例4:
本实施例在上述实施例1-3任一项的基础上,为了更好地实现本发明,进一步地,在所述步骤5中对信号传播速度c、仰角β、方向角α和待估计参数g进行计算,具体步骤为:
步骤5.2:设置如下最小化问题约束条件方程:
步骤5.4:对步骤5.2的表达式的等式两边取平方,忽略二次高阶误差项,得到如下表达式:
式中,v=1,...,p;
然后将v=p代入真实值的表达式中,并结合步骤5.3中的误差公式,得到:
采用微分方法展开后,忽略二次高阶误差项,得到如下表达式:
式中,1表示全1矩阵,下标表示全1矩阵的阶数;0表示全0矩阵,下标表示全0矩阵的阶数;
步骤5.8:计算得到在二维平面的情况下:
在三维空间的情况下:
本实施例的其他部分与上述实施例1-3任一项相同,故不再赘述。
以上所述,仅是本发明的较佳实施例,并非对本发明做任何形式上的限制,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化,均落入本发明的保护范围之内。
Claims (6)
1.一种未知信号传播速度下远近场到达时间差测量定位方法,用于对源目标进行定位,其特征在于,包括以下步骤:
步骤1:建立三维和二维的下的位置坐标模型;
步骤2:建立近场信号到达时间差测量模型;
步骤3:建立远场信号到达时间差测量模型;
步骤4:设计待估计参数g,基于估计参数g根据得到的近场信号时间达到时间差测量模型和远场信号达到时间差测量模型计算得到三维和二维下的远近场定位估计模型;
步骤5:计算待估计参数g、信号速度c和源节点位置u的位置坐标参数;
步骤6:使用计算得到待估计参数g、信号速度c和源节点位置u的位置坐标参数,结合得到的远近场定位估计模型进行到达时间差测量定位计算。
2.如权利要求1所述的一种未知信号传播速度下远近场到达时间差测量定位方法,其特征在于,所述步骤1的具体操作为:
步骤1.1:定义已知位置坐标的M+1个信标节点为Si,i=0,1,2,3,...,M;并定义待确定坐标的源节点位置为u;
步骤1.2:定义二维下信标节点Si的位置坐标Si=[xi , yi]T,源节点位置坐标u=[x ,y]T;定义三维下信标节点Si的位置坐标Si=[xi , yi ,zi]T,源节点位置坐标u=[x , y ,z]T;
将位置坐标分别转换为极坐标形式:
二维空间下:Si=[ρi , αi]T,u=[ρ , θ]T;
三维空间下:Si=[ρi , αi ,βi]T,u=[ρ , θ , Φ]T;
3.如权利要求2所述的一种未知信号传播速度下远近场到达时间差测量定位方法,其特征在于,所述步骤2中的近场信号到达时间差测量模型中,对于到达时间差ti采用距离差di和信号速度c的比值进行计算获得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111288597.0A CN113727433A (zh) | 2021-11-02 | 2021-11-02 | 一种未知信号传播速度下远近场到达时间差测量定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111288597.0A CN113727433A (zh) | 2021-11-02 | 2021-11-02 | 一种未知信号传播速度下远近场到达时间差测量定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113727433A true CN113727433A (zh) | 2021-11-30 |
Family
ID=78686462
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111288597.0A Pending CN113727433A (zh) | 2021-11-02 | 2021-11-02 | 一种未知信号传播速度下远近场到达时间差测量定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113727433A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114358094A (zh) * | 2022-03-18 | 2022-04-15 | 成都迅翼卫通科技有限公司 | 基于雷达通信系统的信号去噪方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103731796A (zh) * | 2013-10-10 | 2014-04-16 | 华南理工大学 | 用于远场和近场头相关传输函数的多声源自动测量系统 |
-
2021
- 2021-11-02 CN CN202111288597.0A patent/CN113727433A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103731796A (zh) * | 2013-10-10 | 2014-04-16 | 华南理工大学 | 用于远场和近场头相关传输函数的多声源自动测量系统 |
Non-Patent Citations (1)
Title |
---|
吴晓平: "《未知信号传播速度的TDOA测量远近场统一定位方法》", 《传感技术学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114358094A (zh) * | 2022-03-18 | 2022-04-15 | 成都迅翼卫通科技有限公司 | 基于雷达通信系统的信号去噪方法及系统 |
CN114358094B (zh) * | 2022-03-18 | 2022-06-03 | 成都迅翼卫通科技有限公司 | 基于雷达通信系统的信号去噪方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Unified near-field and far-field localization for AOA and hybrid AOA-TDOA positionings | |
CN102231912A (zh) | 一种基于rssi测距的室内无线传感器网络定位方法 | |
Li et al. | Improved two-step constrained total least-squares TDOA localization algorithm based on the alternating direction method of multipliers | |
KR20130020295A (ko) | 무선센서네트워크 상에서 효율적인 다변측량을 이용한 무선측위방법과 이를 실시하기 위한 프로그램이 기록된 기록매체 | |
Sun et al. | Wireless sensor network-based localization method using TDOA measurements in MPR | |
CN109581281B (zh) | 基于到达时间差和到达频率差的移动目标定位方法 | |
CN105353351A (zh) | 一种基于多信标到达时间差改进型定位方法 | |
Hu et al. | A moving source localization method using TDOA, FDOA and Doppler rate measurements | |
Li et al. | A new close form location algorithm with AOA and TDOA for mobile user | |
Sun et al. | Array geometry calibration for underwater compact arrays | |
CN113727433A (zh) | 一种未知信号传播速度下远近场到达时间差测量定位方法 | |
CN104735779A (zh) | 一种基于troa的nlos传输环境无线定位方法 | |
Meng et al. | Communication aware optimal sensor motion coordination for source localization | |
Zeng et al. | Radio frequency based direction sensing using massive MIMO | |
Nakamura et al. | An accurate technique for simultaneous measurement of 3D position and velocity of a moving object using a single ultrasonic receiver unit | |
Yang et al. | Robust Direct position determination against sensor gain and phase errors with the use of calibration sources | |
Zhou et al. | Optimal location method of spontaneous data fusion based on TDOA/AOA | |
Wang et al. | Robust direct position determination methods in the presence of array model errors | |
Gu et al. | The effect of ground truth accuracy on the evaluation of localization systems | |
CN114051209B (zh) | 一种基于智能反射面和场景几何模型的指纹定位方法 | |
CN112954591B (zh) | 一种协同分布式定位方法及系统 | |
Manikas et al. | Positioning in wireless sensor networks using array processing | |
Qing et al. | Wireless positioning method based on dynamic objective function under mixed LOS/NLOS conditions | |
KR101459915B1 (ko) | 위치인식 방법 | |
Peng et al. | Comparison of wireless location algorithms in mobile communication 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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211130 |
|
RJ01 | Rejection of invention patent application after publication |