CN113038364B - 一种基于uwb技术的tdoa与ds_twr结合的井下二维定位方法 - Google Patents
一种基于uwb技术的tdoa与ds_twr结合的井下二维定位方法 Download PDFInfo
- Publication number
- CN113038364B CN113038364B CN202110212885.1A CN202110212885A CN113038364B CN 113038364 B CN113038364 B CN 113038364B CN 202110212885 A CN202110212885 A CN 202110212885A CN 113038364 B CN113038364 B CN 113038364B
- Authority
- CN
- China
- Prior art keywords
- uwb
- clock
- locator
- time
- poll
- 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
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000005516 engineering process Methods 0.000 title claims abstract description 10
- 238000001914 filtration Methods 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 28
- 238000012937 correction Methods 0.000 claims description 13
- 238000005259 measurement Methods 0.000 claims description 5
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 239000003245 coal Substances 0.000 abstract description 5
- 230000008878 coupling Effects 0.000 abstract description 3
- 238000010168 coupling process Methods 0.000 abstract description 3
- 238000005859 coupling reaction Methods 0.000 abstract description 3
- 230000005540 biological transmission Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/023—Services making use of location information using mutual or relative location information between multiple location based services [LBS] targets or of distance thresholds
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/021—Services related to particular areas, e.g. point of interest [POI] services, venue services or geofences
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W56/00—Synchronisation arrangements
- H04W56/001—Synchronization between nodes
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W56/00—Synchronisation arrangements
- H04W56/0055—Synchronisation arrangements determining timing error of reception due to propagation delay
- H04W56/0065—Synchronisation arrangements determining timing error of reception due to propagation delay using measurement of signal travel time
- H04W56/007—Open loop measurement
- H04W56/0075—Open loop measurement based on arrival time vs. expected arrival time
-
- 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
Abstract
本发明公开了一种基于UWB技术的TDOA与DS_TWR结合的井下二维定位方法,该方法利用了识别卡到达多个UWB定位器之间的到达时间差值TDOA,和识别卡到达一个UWB定位器的TOA值,结合TDOA和TOA的测量计算值,采用Chan算法得到移动点二维坐标。其中,本发明采用了基于Kalman滤波的无线时钟同步方法实现TDOA的计算,实现了亚纳秒级的时钟精度,通过增加对一个基站的TOA测距值,有效提高了煤矿井下二维长距离定位的精度,解决了井下单一TDOA定位算法精度不稳定问题;同时,本发明方法还可以大大降低系统的耦合度,提高系统稳定性。
Description
技术领域
本发明属于煤矿井下UWB无线定位技术领域,涉及一种基于UWB技术的TDOA 与DS_TWR结合的井下二维定位方法。
背景技术
现有技术问题如下:
(1)一种基于ADS_TWR的井下UWB多基站定位区域判定方法(专利CN 107948920A),等间距布置基站,任意两个基站之间的距离小于定位覆盖半径,该专利实现二维定位问题在于:若获得识别卡的位置,需要识别卡与多个基站分别完成多次DS_TWR过程,占用空间时隙资源过多,系统容量和定位实时性差,可应用一维定位。
(2)一种基于UWB技术DS_TWR的井下精确定位方位判定方法及系统(专利CN111970760A),采用了单基站(双天线)测距结构,该专利实现二维定位问题在于:使用该方法只能实现一维定位,无法实现井下二维或多为定位。
(3)一种基于UWB的联合无线时钟同步和TOA的测距系统(专利 CN106842175A),采用了多基站测距结构,在实现TOA测距过程中,在通信的信息中插入无线时钟同步信息,其时钟同步的精度ms级别,其牺牲了时钟同步的精度,提升了系统的网络组建的易操作性。该专利实现二维定位问题在于:无线时钟同步的精度太差不能用于TDOA计算,其多维定位本质上还是使用DS_TWR 的改进形式,系统容量和定位实时性都较差。
发明内容
本发明解决的问题在于克服现有的煤矿井下精确人员定位技术中,二位或三维定位存在的问题,并提供一种适合井下环境的方法,可以提高系统定位精度,降低系统的耦合度,提高系统稳定性,提升系统容量和定位实时性。
本发明采用以下技术方案解决上述问题:
一种基于UWB技术的TDOA与DS_TWR结合的井下二维定位方法,其特征在于在井下巷道中间隔安装UWB定位基站单元,被定位目标佩戴识别卡,UWB定位基站单元由3个(或更多)UWB定位器组成,UWB定位器集中安装在集控器周围, UWB定位器之间通过无线时钟同步算法,实现UWB定位器之间的亚纳秒级别的时钟同步,各UWB定位器通过接收识别卡的Poll消息,系统可计算出定位器间的 TDOA差值;同时,使用DS_TWR方法完成UWB识别卡与UWB主定位器之间TOA测距,据此可以得到识别卡到达三个UWB定位器的TOA值,最终通过Chan算法得到识别卡的位置值。
进一步地,本发明实现流程如下:
(1)实现本发明方法的相关设备包括UWB定位基站单元(UWB定位器M、UWB 定位器S1、UWB定位器S2、集控器)、识别卡、UWB定位服务器等,如图1所示。
(2)实施本发明方法部署方案,如图2所示,沿着井下巷道两侧依次布设 UWB定位器M、集控器、UWB定位器S1、UWB定位器S2,其中三个UWB定位器不能同时共线,间距大于1米,小于20米。
(3)使用无线时钟同步法计算POLL消息到达UWB定位器之间的到达时间差值(TDOA)的主要步骤,如图4所示。
步骤3.1:UWB定位器M(主定位器)间隔150ms发送无线时钟校准包TX_CPP (以下简称CCP包),记录发送时刻本地时钟数据,并通过CAN总线把带有发送时间戳的CCP数据包传输到集控器。
步骤3.2:UWB定位器S1(从定位器1)和UWB定位器S2(从定位器2)接收RX_CCP包,并通过CAN总线把带有接收时间戳的RX_CCP1、RX_CCP2数据包传输到集控器。
步骤3.3:用基于Kalman滤波的时钟跟踪方法计算时钟偏移。
T为时钟周期,Uθ,K和Uγ,K分别是第K周期的时钟偏移和时钟偏移的修正值,θK和γK分别是kT时钟偏移和时钟漂移。在(K+1)T时刻,时钟迭代关系:
(θk+1)=θk-Uθ,k+(γk-Uγ,k)T+ωθ,k
(γk+1)=γk-Uγ,k+ωγ
ωθ,K是时钟偏差方差,ωγ,K是时钟漂移方差。
假设ωk=[ωθ,kωγ,k]T,协方差矩阵为Q,定义向量矩阵:
xk=[θkγk]T,uk=[Uθ,kUγ,k]T.
使用Kalman滤波迭代:
建立预测方程:(Xk+1)=Axk+Buk
其中:
Kalman增益矩阵:
Kk+1=Pk+1|k[Hk+1]T(Rk+1+Hk+1Pk+1|k[Hk+1]T)-1
RK+1是观测噪声的协方差矩阵,HK+1是测量矩阵,单位方阵;
修正:最小均方差误差矩阵:
Pk+1=(1-Kk+1)Pk+1|k
经过Kalman滤波,第K+1时刻时钟周期的修正值另是补偿时钟偏移和时钟漂移,则有UWB定位器S1和UWB定位器S2 的时钟偏移分别为,S1UK+1和S2UK+1。
步骤3.4:接收识别卡发送的POLL消息
UWB定位器M、UWB定位器S1和UWB定位器S2分别接收识别卡发送的POLL 消息,并获取POLL消息的接收时间戳,分别记作tRX_POLL_M、tRX_POLL_S1、tRX_POLL_S2;通过CAN总线发送这些时间戳到集控器。
步骤3.5:t时刻POLL消息时基校准
由于UWB定位器M、UWB定位器S1和UWB定位器S2设备之间的时钟系统是独立的,即非同的步时钟系统,它们所获得的tRX_POLL_M、tRX_POLL_S1、tRX_POLL_S2时间戳为非同步时间,根据步骤3.3所计算得到的时钟偏移S1UK+1和S2UK+1,把上述时间校准到以UWB定位器M为时基时间值,则有如下时间:
TRX_POLL_M=tRX_POLL_M
TRX_POLL_S1=tRX_POLL_S1+S1UK+1
TRX_POLL_S2=tRX_POLL_S1+S2UK+1
步骤3.6:计算t时刻识别卡到达UWB定位器M、UWB定位器S1和UWB定位器S2设备的距离差(TDOA)
dm_S1(t)是t时刻识别卡到UWB定位器M和UWB定位器S1距离的差值;
dm_S2(t)是t时刻识别卡到UWB定位器M和UWB定位器S2距离的差值;
dm_S1(t)=(TRX_POLL_S1-TRX_POLL_M)*C光速
dm_S2(t)=(TRX_POLL_S2-TRX_POLL_M)*C光速
(4)使用DS_TWR实现识别卡的TOA测距
应用DS_TWR方法UWB定位器M得到TOA结果。本实施例中,DS_TWR的时序方法,如图3所示。
计算识别卡与UWB定位器M的传播时间值,计算公式如下:
TTA:识别卡Tag到UWB定位器A的传播时间;
Tround1是从识别卡开始发送并收到返回信息的首轮时间,T是UWB定位器A的间隔时间,Tround2是从UWB定位器A开始发送并收到返回信息的次轮时间,TM是移动节点中的间隔时间。
计算t时刻识别卡到UWB定位器M的距离
DTM=TAT*C光速
(5)使用TDOA、DS_TWR结果计算t时刻识别卡坐标
使用Chan算法计算二维、三定位基站情况下的t时刻识别卡坐标,如图7所示。
UWB定位器M坐标(x1,y1),UWB定位器S1坐标(x2,y2),UWB定位器S2 坐标(x3,y3);r1、r2、r3分别是识别卡到UWB定位器M、S1、S2的距离;
设识别卡的位置坐标(x,y);
根据几何关系公式可知下列表达式:
根据公式(5-1)有如下关系:
ri 2=(xi-x)2+(yi-y)2=Ki-2xix-2yiy+x2+y2
令:Ki=xi 2+yi 2
其中ri 2=(ri,1+r1)2
则有如下矩阵:即Ax=C。可利用矩阵的逆即 x=A-1C。
最后,求得:
式(5-6)中三个基站的坐标都是已知,r2,1和r3,1是dm_S1(t)、dm_S12t)的TDOA 差值,r1是TOA距离DTM,因此,可直接求得识别卡t时刻二维坐标值x,y。
在识别卡远离三基站的时候,双曲线的渐近线会趋于平行,这产生强烈的误差放大,直接导致了TDOA求解结果数据的精度极具下降,因此,在加入一个 TOA测距结果可以有效降低误差。
本发明的有益效果
(1)本发明方法解决煤矿井下巷道环境的二维定位精度问题,采用TDOA与 DS_TWR结合方法,限制了TDOA的误差抖动问题,解决了识别卡远离三定位基站的定位误差问题,提高了系统定位精度,完全满足煤矿定位技术要求;
(2)本发明方法中TDOA使用了无线时钟同步,采用Kalman滤波跟踪主时钟,时钟校准精度达到亚纳秒级,有效抑制了测量过程中的噪声和干扰数据,;
(3)本发明方法通过集中布放三定位基站阵,可以得到更大的覆盖区域,从而减少井下了定位基站的数量50%,大大降低了系统建设成本以及施工的维护的工作量;
(4)本发明方法中由于采用了单独区域的时钟同步方法,无需实现跨区域全网级联时钟同步,大大降低了系统的耦合程度,杜绝了因本级时钟问题影响后面基站群的定位,提高定位系统稳定性。
附图说明
图1:UWB最小系统结构图;
图2:井下系统布置结构图;
图3:DS_TWR时序图;
图4:TDOA时序图;
图5:TDOA和DS_TWR通信流程图;
图6:TDOA与DS_TWR结合的井下二维定位方法流程图;
图7:TDOA和DS_TWR结合的Chan算法几何示意图。
具体实施方式
为使本发明实施例的目的和技术方案更加清楚,下面将附图,并结合本发明实施例对本发明的技术方案进行清楚、完整的描述。显然,所描述实施例是本发明的一部分实施例,而不是全部实施例。基于所描述的本发明的实施例,本领域普通技术人员在无需创造性劳动的前提下,所获得的所有其他实施例,都属于本发明的保护范围。
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义解释。
进一步地,本发明方法步骤:
步骤1:实现本发明方法的相关系统设备组成,如图1所示。
UWB定位基站单元:由UWB定位器M、UWB定位器S1、UWB定位器S2、集控器共同组成,UWB定位器与集控器采用RS-485或CAN总线连接,并实时上报接收到的定位数据和时钟同步数据,并在集控器内实现初步的位置信息解算。
UWB定位器M是主定位器,定时发送UWB无线时钟校准包TX_CCP,并完成与识别卡的DS_TWR测距,并上报这些消息。
UWB定位器S1是从定位器1,完成接收UWB无线时钟校准包RX_CCP,并上报。
UWB定位器S2是从定位器2,完成接收UWB无线时钟校准包RX_CCP,并上报。
UWB集控器是收集TX_CCP消息、RX_CCP消息、POLL消息,并计算TDOA,接收TOA消息。
UWB识别卡能够实现DS_TWR的UWB无线收发设备。
UWB定位服务器实现位置坐标解算。
步骤2:实施本发明方法的设备部署方案,如图2所示。
沿着井下巷道两侧依次布设UWB定位器M、集控器、UWB定位器S1、UWB定位器S2,其中三个UWB定位器不能同时共线,间距大于1米,小于20米(可依据CAN总线和RS-485传输速率、距离确定)该实际距离可通过系统自动测量(或人为实际测量)得到。
进一步,步骤3:使用基于Kalman滤波无线时钟同步法计算得到POLL消息的TDOA时间差值
步骤3.1:UWB定位器M(主定位器)间隔150ms发送无线时钟校准包TX_CPP (以下简称CCP包),记录发送时刻本地时钟数据,并通过CAN总线把带有发送时间戳的CCP数据包传输到集控器。
步骤3.2:UWB定位器S1(从定位器1)和UWB定位器S2(从定位器2)接收RX_CCP包,并通过CAN总线把带有接收时间戳的RX_CCP1、RX_CCP2数据包传输到集控器。
步骤3.3:用基于Kalman滤波的时钟跟踪方法计算时钟偏移
Kalman滤波采用最小均方差误差实现最佳的线性滤波方法。根据当前时刻的观测值、上一时刻的预测值以及预测误差来计算得到当前的最优值去预测下一时刻的值。第K+1同步周期主从时钟的时钟偏差和时钟漂移根据第K周期的信息预测。则根据Kalman标准方程可得到:
T为时钟周期,Uθ,K和Uγ,K分别是第K周期的时钟偏移和时钟偏移的修正值,θK和γK分别是kT时钟偏移和时钟漂移。在(K+1)T时刻,时钟迭代关系:
(θk+1)=θk-Uθ,k+(γk-Uγ,k)T+ωθ,k
(γk+1)=γk-Uγ,k+ωγ
ωθ,K是时钟偏差方差,ωγ,K是时钟漂移方差。假设ωk=[ωθ,kωγ,k]T,协方差矩阵为Q,
定义向量矩阵:xk=[θkγk]T,uk=[Uθ,kUγ,k]T.
使用Kalman滤波迭代:
建立预测方程:(Xk+1)=Axk+Buk
其中:Kalman增益矩阵:
Kk+1=Pk+1|k[Hk+1]T(Rk+1+Hk+1Pk+1|k[Hk+1]T)-1
RK+1是观测噪声的协方差矩阵,HK+1是测量矩阵,单位方阵;
修正:
最小均方差误差矩阵:
Pk+1=(1-Kk+1)Pk+1|k
经过Kalman滤波,第K+1时刻时钟周期的修正值
另是补偿时钟偏移和时钟漂移,则有UWB定位器S1和UWB定位器S2的时钟偏移分别为,S1UK+1和S2UK+1。
步骤3.4:接收识别卡发送的POLL消息
UWB定位器M、UWB定位器S1和UWB定位器S2分别接收识别卡发送的POLL 消息,并获取POLL消息的接收时间戳,分别记作tRX_POLL_M、tRX_POLL_S1、tRX_POLL_S2;通过CAN总线发送这些时间戳到集控器。
步骤3.5:t时刻POLL消息时基校准
由于UWB定位器M、UWB定位器S1和UWB定位器S2设备之间的时钟系统是独立的,即非同的步时钟系统,它们所获得的tRX_POLL_M、tRX_POLL_S1、tRX_POLL_S2时间戳为非同步时间,根据步骤3.3所计算得到的时钟偏移S1UK+1和S2UK+1,把上述时间校准到以UWB定位器M为时基时间值,则有如下时间:
TRX_POLL_M=tRX_POLL_M
TRX_POLL_S1=tRX_POLL_S1+S1UK+1
TRX_POLL_S2=tRX_POLL_S1+S2UK+1
步骤3.6:计算t时刻识别卡到达UWB定位器M、UWB定位器S1和UWB定位器S2设备的距离差(TDOA)
dm_S1(t)是t时刻识别卡到UWB定位器M和UWB定位器S1距离的差值;
dm_S2(t)是t时刻识别卡到UWB定位器M和UWB定位器S2距离的差值。
dm_S1(t)=(TRX_POLL_S1-TRX_POLL_M)*C光速
dm_S2(t)=(TRX_POLL_S2-TRX_POLL_M)*C光速
进一步,步骤4:使用DS_TWR方法得到识别卡的TOA测距值
步骤4.1识别卡向UWB定位器M发送POLL消息,并记录发送时刻时间TMtx(发送Poll消息的时间戳),UWB定位器M收到此消息并记录时间TArx(接收Poll消息的时间戳);
步骤4.2UWB定位器M向识别卡Respond消息,并记录发送时刻时间TAtx(发送Respon消息时间戳),识别卡s收到此消息,并记录时间戳TMrx(接收UWB定位器A的Respon消息时间戳);
步骤4.3识别卡向UWB定位器M发送Final消息,并记录发送时刻时间TMtx(发送Final消息的时间戳),UWB定位器M收到此消息并记录时间戳TBrx(接收Final消息的时间戳)。
步骤4.4UWB定位器M计算TOA结果
本实施例中,DS_TWR的时序方法,见图3所示。
(1)测量的识别卡时间戳
TMtx:发送Poll消息的时间戳;
TMrx:接收UWB定位器M的Respon消息时间戳;
TMtx:发送Final消息的时间戳。
(2)测量的UWB定位器M时间戳
TArx:接收Poll消息的时间戳;
TAtx:发送Respon消息时间戳;
TArx:接收Final消息的时间戳。
计算识别卡与UWB定位器M的传播时间值,计算公式如下:
TTA:识别卡Tag到UWB定位器A的传播时间;
Tround1是从识别卡开始发送并收到返回信息的首轮时间,T是UWB定位器A的间隔时间,Tround2是从UWB定位器A开始发送并收到返回信息的次轮时间,TM是移动节点中的间隔时间。
计算t时刻识别卡到UWB定位器M的距离
DTM=TAT*C光速
进一步,步骤5:根据TDOA与DS_TWR联合定位结果,使用Chan方法计算识别卡坐标值
步骤5.1计算t时刻识别卡与UWB定位器S之间的距离
根据步骤3.6中计算的t时刻识别卡发送POLL消息到达UWB定位器S的到达时间差(TDOA),以及步骤4.4的t时刻识别卡到UWB定位器M距离(TOA),有一下过程:
DTS1=DTM(t)+dm_S1(t)
DTS2=DTM(t)+dm_S2(t)
其中:DTS1、DTS2分别是t时刻识别卡到达UWB定位器S1和S2的距离值。
步骤5.2计算t时刻识别卡位置坐标
使用Chan算法计算二维、3定位基站情况下的t时刻识别卡坐标
(1)UWB定位器M坐标(x1,y1),UWB定位器S1坐标(x2,y2),UWB定位器 S2坐标(x3,y3);r1、r2、r3分别是识别卡到UWB定位器M、S1、S2的距离;
设识别卡坐标(x,y)位置;
据几何关系公式可知下列表达式:
根据公式(5-1)有如下关系:
ri 2=(xi-x)2+(yi-y)2=Ki-2xix-2yiy+x2+y2 (5-3)
令:Ki=xi 2+yi 2
根据公式(5-2)有如下关系:
ri 2=(ri,1+r1)2 (5-4)
将(5-3)带入(5-4)有如下关系:
则有如下矩阵:即Ax=C。可利用矩阵的逆即x=A-1C。最后,求得:
式(5-6)中三个基站的坐标都是已知,r2,1和r3,1分别是dm_S1(t)、dm_S12t)的 TDOA差值,r1是TOA距离DTM,因此,可直接求得t时刻识别卡二维坐标值x,y。
在识别卡远离三基站定位的时候,双曲线的渐近线会趋于平行,这产生强烈的误差放大,直接导致了TDOA求解结果数据的精度极具下降,因此,在加入一个TOA测距结果可以有效降低误差。
Claims (3)
1.一种基于UWB技术的TDOA与DS_TWR结合的井下二维定位方法,其特征在于,包括以下步骤:
步骤一:实现本发明方法的相关系统设备由UWB定位器M、UWB定位器S1、UWB定位器S2、集控器共同组成;
步骤二:实施本发明方法的设备部署方案为:沿着井下巷道两侧依次布设所述UWB定位器M、集控器、UWB定位器S1、UWB定位器S2,其中三个UWB定位器不能同时共线,间距大于1米,小于20米;
步骤三:使用基于Kalman滤波无线时钟同步法计算得到POLL消息TDOA时间差值;
其中,步骤三包括:
步骤3.1:UWB定位器M间隔150ms发送无线时钟校准包CPP,记录发送时刻本地时钟数据,并通过CAN总线把带有发送时间戳的CCP数据包传输到集控器;步骤3.2:UWB定位器S1和UWB定位器S2接收RX_CCP包,并通过CAN总线把带有接收时间戳的RX_CCP1、RX_CCP2数据包传输到集控器;
步骤3.3:用基于Kalman滤波的时钟跟踪方法计算时钟偏移,则有UWB定位器S1和UWB定位器S2的时钟偏移分别为,S1UK+1和S2UK+1;
步骤3.4:UWB定位器M、UWB定位器S1和UWB定位器S2分别接收识别卡发送的POLL消息,并获取POLL消息的接收时间戳,分别记作tRX_POLL_M、tRX_POLL_S1、tRX_POLL_S2;通过CAN总线发送这些时间戳到集控器;
步骤3.5:根据所述时钟偏移S1UK+1和S2UK+1,把tRX_POLL_M、tRX_POLL_S1、tRX_POLL_S2校准到以UWB定位器M为时基时间值,则有如下时间:
TRX_POLL_M=tRX_POLL_M
TRX_POLL_S1=tRX_POLL_S1+S1UK+1
TRX_POLL_S2=tRX_POLL_S1+S2UK+1;
步骤3.6:计算t时刻识别卡到达UWB定位器M和UWB定位器S1距离的差值和识别卡到UWB定位器M和UWB定位器S2设备的距离差(TDOA);
dm_S1(t)是t时刻识别卡到UWB定位器M和UWB定位器S1距离的差值;
dm_S2(t)是t时刻识别卡到UWB定位器M和UWB定位器S2距离的差值;
C光速是光速;
dm_S1(t)=(TRX_POLL_S1-TRX_POLL_M)*C光速
dm_S2(t)=(TRX_POLL_S2-TRX_POLL_M)*C光速;
步骤四:使用DS-TWR方法得到识别卡的TOA测距值;计算t时刻识别卡到UWB定位器M的距离;
TTA:识别卡Tag到UWB定位器A的传播时间;
Tround1是从识别卡开始发送并收到返回信息的首轮时间,T是UWB定位器A的间隔时间,Tround2是从UWB定位器A开始发送并收到返回信息的次轮时间,TM是移动节点中的间隔时间;
计算t时刻识别卡到UWB定位器M的距离:
DTM=TTA*C光速
步骤五:根据TDOA与DS-TWR联合定位结果,使用Chan方法计算t时刻识别卡坐标值,具体的:
UWB定位器M坐标(x1,y1),UWB定位器S1坐标(x2,y2),UWB定位器S2坐标(x3,y3);r1、r2、r3分别是识别卡到UWB定位器M、S1、S2的距离;
设识别卡的位置坐标(x,y);令:K=xi 2+yi 2;
得:
式中:y2,1是指y2,1=y2-y1;y3,1是指y3,1=y3-y1;x3,1是指x3,1=x3-x1;
x2,1是指x2,1=x2-x1;
r2,1和r3,1是dm_S1(t)、dm_S12(t)的TDOA差值,r1是TOA距离DTM,因此,可直接求得识别卡t时刻二维坐标值x,y。
2.根据权利要求1所述的一种基于UWB技术的TDOA与DS_TWR结合的井下二维定位方法,其特征在于:UWB定位器与集控器采用RS-485或CAN总线连接,并实时上报接收到的定位数据和时钟同步数据,并在集控器内实现初步的位置信息解算。
3.根据权利要求1所述的一种基于UWB技术的TDOA与DS_TWR结合的井下二维定位方法,其特征在于:使用基于Kalman滤波方法实现无线时钟同步法,并计算POLL消息到达各UWB定位器的到达时间差值(TDOA);
T为时钟周期,
Uθ,K是第K周期的时钟偏移修正值;
Uγ,K是第K周期的时钟漂移修正值;
θK是kT时刻时钟偏移;
γK是kT时刻时钟漂移;
在(K+1)T时刻,时钟迭代关系:
(θk+1)=θk-Uθ,k+(γk-Uγ,k)T+ωθ,k
(γk+1)=γk-Uγ,k+ωγ,k
ωθ,K是时钟偏差方差,ωγ,k是时钟漂移方差;
假设ωk=[ωθ,k ωγ,k]T,
ωk是时钟漂移方差和时钟偏移方差矩阵;
xk=[θkγk]T,uk=[Uθ,kUγ,k]T.
xk是kT时刻时钟偏移和时钟漂移矩阵;
uk是时钟偏移修正值和时钟漂移修正值矩阵;
使用Kalman滤波迭代:
建立预测方程:(Xk+1)=Axk+Buk
其中:
Kalman增益矩阵:
Kk+1=Pk+1|k[Hk+1]T(Rk+1+Hk+1Pk+1|k[Hk+1]T)-1
Kk+1是在(K+1)T时刻Kalman增益矩阵;
Pk+1|k是最小均方误差矩阵估计;
Rk+1是观测噪声的协方差矩阵,Hk+1是测量矩阵,单位方阵;修正:
Zk+1是第K+1时刻时钟周期误差;
是经过Kalman滤波,第K+1时刻时钟周期的修正值;
是第K时刻时钟周期对第K+1时刻的时钟预测值;
最小均方差误差矩阵:
Pk+1=(1-Kk+1)Pk+1|k
经过Kalman滤波,第K+1时刻时钟周期的修正值另是补偿时钟偏移和时钟漂移,则有UWB定位器S1和UWB定位器S2的时钟偏移分别为,S1Uk+1和S2Uk+1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110212885.1A CN113038364B (zh) | 2021-02-25 | 2021-02-25 | 一种基于uwb技术的tdoa与ds_twr结合的井下二维定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110212885.1A CN113038364B (zh) | 2021-02-25 | 2021-02-25 | 一种基于uwb技术的tdoa与ds_twr结合的井下二维定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113038364A CN113038364A (zh) | 2021-06-25 |
CN113038364B true CN113038364B (zh) | 2024-05-03 |
Family
ID=76462181
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110212885.1A Active CN113038364B (zh) | 2021-02-25 | 2021-02-25 | 一种基于uwb技术的tdoa与ds_twr结合的井下二维定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113038364B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113687301A (zh) * | 2021-08-20 | 2021-11-23 | 中煤科工集团重庆研究院有限公司 | 一种煤矿井下多支路口uwb联合定位方法 |
CN113706612B (zh) * | 2021-10-28 | 2022-02-11 | 天地(常州)自动化股份有限公司 | 融合uwb和单目视觉slam的煤矿井下车辆定位方法 |
CN114007261B (zh) * | 2022-01-05 | 2022-06-10 | 华清科盛(北京)信息技术有限公司 | 一种基于uwb定位系统的航吊吊钩三维定位方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20090002764A (ko) * | 2007-07-04 | 2009-01-09 | 한국전기연구원 | 무선측위에서의 클록 표류 영향을 고려한 시간정보추정방법 |
EP2117546A4 (en) * | 2007-02-05 | 2011-05-25 | Aeroscout Ltd | ARROW TIME DIFFERENCE SYSTEM (TDOA) WITH DOUBLE BAND WIDTH |
CN105376848A (zh) * | 2014-08-19 | 2016-03-02 | 丹东东方测控技术股份有限公司 | 一种实现井下矿人员及设备连续高精度定位的方法 |
CN106793060A (zh) * | 2017-03-08 | 2017-05-31 | 哈尔滨工程大学 | 一种超宽带室内定位方法 |
CN110099354A (zh) * | 2019-05-21 | 2019-08-06 | 合肥工业大学 | 一种结合tdoa与tof的超宽带通信二维定位方法 |
CN110174641A (zh) * | 2019-06-14 | 2019-08-27 | 中国科学院自动化研究所 | 基于时间校正的tdoa定位方法、系统、装置 |
CN111654908A (zh) * | 2020-07-03 | 2020-09-11 | 安徽理工大学 | 基于时间敏感网络的井下混合网络时间同步方法 |
CN111918207A (zh) * | 2020-08-10 | 2020-11-10 | 北京瑞华高科技术有限责任公司 | 定位系统及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011069552A1 (en) * | 2009-12-10 | 2011-06-16 | Nortel Networks Limited | Method, arrangement and computer program product for clocking |
-
2021
- 2021-02-25 CN CN202110212885.1A patent/CN113038364B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2117546A4 (en) * | 2007-02-05 | 2011-05-25 | Aeroscout Ltd | ARROW TIME DIFFERENCE SYSTEM (TDOA) WITH DOUBLE BAND WIDTH |
KR20090002764A (ko) * | 2007-07-04 | 2009-01-09 | 한국전기연구원 | 무선측위에서의 클록 표류 영향을 고려한 시간정보추정방법 |
CN105376848A (zh) * | 2014-08-19 | 2016-03-02 | 丹东东方测控技术股份有限公司 | 一种实现井下矿人员及设备连续高精度定位的方法 |
CN106793060A (zh) * | 2017-03-08 | 2017-05-31 | 哈尔滨工程大学 | 一种超宽带室内定位方法 |
CN110099354A (zh) * | 2019-05-21 | 2019-08-06 | 合肥工业大学 | 一种结合tdoa与tof的超宽带通信二维定位方法 |
CN110174641A (zh) * | 2019-06-14 | 2019-08-27 | 中国科学院自动化研究所 | 基于时间校正的tdoa定位方法、系统、装置 |
CN111654908A (zh) * | 2020-07-03 | 2020-09-11 | 安徽理工大学 | 基于时间敏感网络的井下混合网络时间同步方法 |
CN111918207A (zh) * | 2020-08-10 | 2020-11-10 | 北京瑞华高科技术有限责任公司 | 定位系统及方法 |
Non-Patent Citations (3)
Title |
---|
Performance analysis of Kalman filter-based clock synchronization in IEEE 1588 networks;G. Giorgi, C. Narduzzi;《2009 International Symposium on Precision Clock Synchronization for Measurement, Control and Communication》;第 III~IV章 * |
Performance analysis of Kalman filter-based clock synchronization in IEEE 1588 networks;G. Giorgi, C. Narduzzi;《2009 International Symposium on Precision Clock Synchronization for Measurement, d Communication》;第III~IV章 * |
联合TDOA改进算法和卡尔曼滤波的UWB室内定位研究;张桀;沈重;;现代电子技术(13);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113038364A (zh) | 2021-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113038364B (zh) | 一种基于uwb技术的tdoa与ds_twr结合的井下二维定位方法 | |
US20240069147A1 (en) | Method and system for multipath reduction for wireless synchronization and/or locating | |
CN1480873B (zh) | 用于检测终端在网络中位置的方法、系统和装置 | |
CN110856106B (zh) | 基于uwb和气压计的室内高精度三维定位方法 | |
US20120252503A1 (en) | System and method for passively determining own position listening to wireless time synchronization communications | |
CN103344942B (zh) | 控制节点、异步定位方法与系统 | |
CN110167135B (zh) | 一种免时钟同步的tdoa无线定位方法及系统 | |
CN108370551A (zh) | 基于到达时间差定位方法、用户设备及网络设备 | |
CN112073903A (zh) | 一种单基站高精度uwb室内定位系统及方法 | |
CN103675804A (zh) | 一种基于双星时间同步的星间测距方法 | |
CN108235427A (zh) | 一种测量Tof和Toda的方法 | |
CN109005584A (zh) | 基于tdoa技术的定位系统的无线时钟同步方案 | |
CN111193997A (zh) | 一种uwb定位系统的到达时间差测量与校准方法 | |
US20230110357A1 (en) | Wireless positioning | |
KR101260570B1 (ko) | 실시간 위치추적시스템에서의 시각 동기화 장치 및 그 방법 | |
CN115706592A (zh) | 利用基于到达角的干扰补偿的超宽带距离确定 | |
Chiasson et al. | Asynchronous hyperbolic UWB source-localization and self-localization for indoor tracking and navigation | |
KR101304849B1 (ko) | 클럭 오프셋 보정을 이용하는 무선 측위 서버 및 이를 이용한 무선 측위 방법 | |
US7359711B2 (en) | Methods and apparatus for improving accuracy of radio timing measurements | |
US20220385519A1 (en) | MOBILE DEVICE FREQUENCY OFFSET DETERMINATION AND TDoA LOCALIZATION | |
Yi et al. | Joint time synchronization and tracking for mobile underwater systems | |
CN110658491A (zh) | 测向系统、测向方法、定位系统及定位方法 | |
Hach et al. | Wireless synchronization in time difference of arrival based real time locating systems | |
CN113038365B (zh) | 一种标签无线定位方法和系统 | |
EP2420855A1 (en) | Location system and corresponding calibration method |
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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20240407 Address after: 2-4-2, 14 Huanteng street, Xigang District, Dalian, Liaoning, 116011 Applicant after: Wang Xuehui Country or region after: China Address before: 2-1, Xinghai green mountain city, 212 Shengli East Road, Xigang District, Dalian, Liaoning, 116011 Applicant before: Yang Yifei Country or region before: China |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant |