CN1202688C - 利用信号达到时间和时间差估计移动台位置的方法及装置 - Google Patents
利用信号达到时间和时间差估计移动台位置的方法及装置 Download PDFInfo
- Publication number
- CN1202688C CN1202688C CN 01136462 CN01136462A CN1202688C CN 1202688 C CN1202688 C CN 1202688C CN 01136462 CN01136462 CN 01136462 CN 01136462 A CN01136462 A CN 01136462A CN 1202688 C CN1202688 C CN 1202688C
- Authority
- CN
- China
- Prior art keywords
- signal
- travelling carriage
- base station
- propagation time
- time difference
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Mobile Radio Communication Systems (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种利用信号达到时间估计移动台位置的估计方法,该方法通过对与移动台连接的所有基站的TOA以及相对主基站TDOA测量,获取与移动台连接的所有基站的TOA和相对主基站TDOA的有关统计信息;再根据基站位置信息以及所有TOA和TDOA的统计信息构造运算矩阵,两次使用加权线性最小二乘来估计移动台位置,最后利用最小距离差挑选最终位置估计值;该方法能够充分利用所有的TOA和TDOA信息,使移动台的位置估计精度,尤其是当移动台处于切换状态下的定位性能提高明显,以及GDOP性能均得到提高,同时该方法的位置估计过程无须迭代求解,运算量小,效率较高。本发明同时还提供了一种利用上述方法的移动台位置的估计装置。
Description
技术领域
本发明涉及移动通信系统中的无线定位技术,具体地说涉及移动通信系统中利用信号到达时间以及时间差进行移动台定位的方法及装置。
背景技术
通常,在对移动台进行定位时,必有基站与移动台保持连接,就信号测量而言,与移动台保持连接的基站的TOA(Time Of Arrival:单个信号到达时间)测量比较容易实现,由于不同基站相对某一参考基站的TDOA(Time Difference Of Arrival:两个信号到达时间差)测量也容易实现,使得与移动台保持连接的基站的TOA测量以及相对该基站的TDOA测量容易实现。因此,在移动通信系统中通常利用时间信息对移动台进行定位,即利用TOA进行定位和利用TDOA进行定位。
一般而言,直接利用TDOA进行进行移动台的位置估计,尽管信号的到达时间差TDOA能够消除一定程度的系统误差以及部分NLOS(非可视距离:Non_Line of Sight)误差,不同基站的TDOA测量也容易实现,但移动台的位置估计的性能对移动台与基站之间的相对位置比较敏感,即该类型方法的GDOP(几何精度因子:Geomtry Dilution of Precision)性能较差,特别是当移动台靠近基站时,该方法的位置估计性能明显下降。而单存利用TOA进行位置估计时,尽管GDOP性能较好,对移动台与基站位置敏感性较低,但移动台至不同基站的TOA直接测量较难实现,同时对系统固有误差以及信号的NLOS误差无抑制作用。因此,同时采用TOA和TDOA对移动台进行位置估计可以获得较高的精度。然而,现有方法只利用了一个与移动台保持连接的基站的TOA和其它相对于该基站的其它基站的TDOA,事实上,当移动台处于切换状态时,有多个基站与移动台保持连接,此时对移动台进行定位时,能得到多个TOA信息,同时相对某参考基站的一组TDOA信息也比较容易易得到。由于现有的利用TOA以及TDOA的方法并不能将所有的TOA以及TDOA测量信息利用上,而仅仅只是对部分信息的利用,使得移动台的定位精度并不能达到理想的要求。
发明内容
本发明的目的在于,提供一种利用信号传播时间和信号传播时间差、精度较高的估计移动台位置的方法及装置。
为达到上述目的,本发明提供的利用信号达到时间估计移动台位置的方法,包括:
(1)确定参与进行移动台定位的基站,发起所有与移动台保持连接的基站的信号到达时间的测量,获得有关基站信号的往返时间测量值、移动台的收发时间差测量值;
(2)将上述步骤(1)得到的一组信号到达时间的测量值转换为相应的信号从基站到移动台的传播时间,然后选取最小传播时间对应的基站为主基站;
(3)发起与所选主基站为参考基站的信号到达时间差测量,获得相对参考基站的一组信号到达时间差测量值,以及其它基站相对于主基站的信号发射时间差测量值,并将所得的上述两个测量值转换为信号从两个不同基站到同一个移动台的传播时间差;
(4)获取参与移动台定位的基站的坐标以及所有信号从基站到移动台的传播时间和信号从两个不同基站到同一个移动台的传播时间差的有关统计信息;
(5)利用上述步骤(4)提供的基站坐标以及所有信号从基站到移动台的传播时间、信号从两个不同基站到同一个移动台的传播时间差及其相关的统计信息进行移动台的位置估计,确定移动台的最终位置估计值。
上述步骤(5)进一步包括下述步骤:
(31)确定利用信号从基站到移动台的传播时间以及信号从两个不同基站到同一个移动台的传播时间差进行移动台位置估计的基本方程如下:
式中:i-1,2,.m,j-2,.m,而m>-3,i等于1的基站为移动台所属的主基站;ri=c×τi,τi为第i个基站的信号到移动台的传播时间,rj1=c×τj1,τj1为第j个基站相对主基站的信号传播时间差,c为信号传播的速度;xi、yi为第i个基站的几何位置坐标。x、y为待估计的移动台的几何位置坐标。
(32)根据上述步骤(31)的基本方程,使用加权线性最小二乘法估计移动台的粗略位置,得到移动台位置的初步估计解;
(33)对上述步骤(32)得到的移动台的位置估计值再次使用加权线性最小二乘法进行优化,降低所述位置估计值之间的相关性影响;
(34)根据上述步骤(33)得到的移动台的优化位置估计值,确定移动台最终的位置估计值。
其中所述步骤(32)中使用加权线性最小二乘法估计移动台的粗略位置进一步包括:
(41)对上述步骤(31)的基本方程进行变换,得到下述公式:
r1-d1=0
...
...
式中:i=2,.m,而m>-3,第1个基站为主基站;ri=c×τi,τi为第i个基站的信号到达移动台的传播时间,c为信号传播的速度;ri1=c×τi1,τi1为第i个基站信号相对第1个基站信号的传播时间差;xi、yi为第i个基站的几何位置坐标;
d2=(x-x1)2+(y-y1)2;x、y为待估计的移动台的几何位置坐标;
(42)构造矩阵h、Ga、Ksi,将上述步骤(41)中的公式表述为下述最小二乘的形式:Δ=h-GaZa;
式中:Δ为残差;Ksi为进行最小二乘估计时的加权系数矩阵,
Ksi=E(ΔΔT)=c2BQB
其中:B=diag{1,r1,r2,...,rm,2(r21+r1),2(r31+r1),...,2(rm1+r1)},diag表示对角阵,Q为基站信号到移动台的传播时间以及信号从两个不同基站到同一个移动台的传播时间差的误差的协方差矩阵,表示信号可信程度;
式中:Δ为残差;
(43)根据所述矩阵h、Ga、Ksi计算Za,得到移动台位置的初步估计解。在所述步骤(43)中,按下述方法计算Za:
Za=(GaTKsi-1Ga)-1GaTKsi-1h。
所述步骤(33)进一步包括:
(71)构造矩阵h′、Ga′、Ksi′,确定下述最小二乘形式:
Δ′=h′-Ga′Za′
(72)根据所述矩阵h′、Ga′、Ksi′计算Za′,得到移动台位置的精确位置解。
在所述步骤(72)中,按下述方法计算Za′:
Za′=(Ga′TKsi′-1Ga′)-1Ga′TKsi′-1h′。
所述步骤(34)进一步包括以下步骤:
(101)按照下述方法构造矩阵
得到移动台最终位置解;
(102)根据移动台位置的初步估计解与所述最终位置解的差值,从上述步骤(101)中的移动台的最终位置解中,挑选对应移动台位置的初步估计解与所述最终位置解的差值最小的解作为移动台的最终位置估计值。
本发明同时还提供了一种利用信号达到时间估计移动台位置的装置,包括:
信号传播时间和信号传播时间差生成器:用于接收所有与移动台保持连接的基站的信号传播时间与参与移动台定位的其它基站的相对于主基站的信号传播时间差相对应的测量值,所述主基站为最小传播时间对应的基站或信号质量最好的基站,所述测量值包括:所有与移动台保持连接的基站的往返时间测量值、移动台的信号收发时间差测量值,以及相对于主基站的其它基站的信号到达时间差测量值和信号发射时间差测量值,将上述测量值按下述方法转换为与移动台保持连接的基站到移动台的信号传播时间值和其它基站相对于主基站的信号传播时间差值:
基站到移动台的信号传播时间=(基站信号的往返时间测量值-移动台信号收发时间差测量值)/2;相对主基站的其它基站的信号传播时间差值=相应的信号到达时间差测量值-相应的信号发射时间差测量值;
信息数据库:用于存储参与移动台定位的所有基站的几何位置和有关的统计信息;
位置估计器:用于所述利用信号传播时间和信号传播时间差生成器输出的与移动台保持连接的一组信号传播时间和相对于主基站的一组信号传播时间差值和所述信息数据库提供的基站的几何位置以及有关的统计信息,进行移动台的位置估计,确定移动台的最终位置估计值。
所述位置估计器包括:
存储器:用于接收包括与移动台保持连接的所有基站信号到移动台的传播时间值、相对主基站的其它基站的信号到移动台的传播时间差值、基站坐标、所有信号传播时间以及信号传播时间差的统计信息,进行移动台位置估计时所需的数据,以及存储移动台位置估计时产生的中间数据和最终数据;
处理器:用于根据从所述存储器得到的进行移动台位置估计所需的一组信号传播时间数据和一组信号传播时间差数据,利用最小二乘法进行移动台位置的估计。
由于本发明能够综合利用与移动台保持连接的多个基站的一组信号传播时间和相对主基站的其它基站的一组信号传播时间差,利用最小二乘法进行移动台位置估计,更适合于移动台处于切换状态下的定位,这样,本发明在移动台的位置估计过程中充分利用了全部的信号传播时间和信号传播时间差信息,因此更能综合利用信号传播时间差进行移动台定位与利用信号传播时间进行移动台定位两种方法的优点,使移动台的位置估计精度以及GDOP性能进一步得到提高;同时本发明进行位置估计过程中无须迭代求解,运算量小,位置估计所需时间也较少,因此采用本发明可以使同一时间段内的系统处理的定位请求数量明显增加。
附图说明
图1是本发明所述方法进行移动台位置估计的几何原理图;
图2是本发明所述方法的实施例流程图;
图3是图2中步骤5的具体实施例流程图;
图4是本发明所述装置的实施例框图;
图5是图4所述的位置估计器的实施例框图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细的描述。
发明提出的方法和装置适用于3个或3个以上的基站参与对移动台进行位置估计。图1给出了综合利用与移动台保持连接的多个基站的一组信号传播时间以及相对于主基站的其它基站的信号传播时间差进行移动台位置估计的几何原理图。图1中,不同圆和双曲线的交点即为移动台的位置,当信号传播时间以及信号传播时间差存在误差时,曲线将不交于一点,而是大致区域。
图2是本发明所述方法的实施例流程图。按照图2实施本发明共具有5个步骤,步骤1确定参与进行移动台定位的基站,发起所有与移动台保持连接的基站的信号传播时间(单个信号到达时间)的测量,获得有关基站的RTT(往返时间:Round Trip Time)测量值、移动台的收发时间差测量值;在步骤2将上述步骤1得到的一组信号传播时间的测量值转换为相应的信号传播时间值,然后选取最小信号传播时间值对应的基站为主基站;在步骤3发起与所选主基站为参考基站的信号传播时间差测量,获得相对参考基站的一组信号到达时间差测量值,以及信号发射时间差测量值,并将所得的测量值转换为相应的信号传播时间差值;在步骤4获取参与移动台定位的基站的坐标以及所有信号传播时间和信号传播时间差的有关统计信息;最后在步骤5利用上述步骤4得到的基站坐标以及所有信号传播时间、信号传播时间差及其相关的统计信息进行移动台的位置估计,确定移动台的最终位置估计值。
上述步骤5的具体实施过程参考图3。图3中所述的流程主要包括三部分内容,第一部分是步骤41到步骤43,采用加权线性最小二乘法得到移动台位置的初步估计解。
首先在步骤41确定利用信号传播时间以及信号传播时间差进行移动台位置估计的基本方程如下:
式中:i-1,2,.m,j-2,.m,而m>-3,i等于1的基站为移动台所属的主基站;ri=c×τi,τi为第i个基站的信号到达移动台的传播时间,rj1=c×τj1,τj1为第j个基站相对主基站的信号传播时间差,c为信号传播的速度(光速);xi、yi为第i个基站的几何位置坐标。x、y为待估计的移动台的几何位置坐标。
在步骤42,首先对上述步骤41的基本方程进行变换,得到下述公式:
r1-d1=0
...
...
式中:1-2,.m,而m>-3,第1个基站为主基站;ri=c×τi,τi为第i个基站的信号到达移动台的传播时间,c为信号传播的速度;ri1=c×τi1,τi1为第i个基站信号相对第1个基站信号的传播时间差;xi、yi为第i个基站的几何位置坐标;
d2=(x-x1)2+(y-y1)2;x、y为待估计的移动台的几何位置坐标。实际中,若信号传播时间以及信号传播时间差的数量不是一一对应时,可以根据相应的信号传播时间与信号传播时间差的关系将不足的信号传播时间或信号传播时间差补足。例如,如果缺少第3个基站的τ3,可利用τ3=τ31+τ1的关系将τ3补齐,使之满足公式(2)的形式。另外也可在构造矩阵h、Ga、Ksi;时,将缺少的信号传播时间或信号传播时间差所对应的行去掉,其他不变化。
然后构造矩阵h、Ga、Ksi,将上述步骤41中的公式表述为下述最小二乘的形式:
Δ=h-GaZa (3);
最小二乘的目的是使‖Δ‖2最小。
式中:Δ即为残差;Ksi为进行最小二乘估计时的加权系数矩阵,
Ksi=E(ΔΔT)=c2BQB (4);
其中:B=diag{1,r1,r2,...,rm,2(r21+r1),2(r31+r1),...,2(rm1+r1)},diag表示对角阵,Q为信号传播时间以及信号传播时间差的误差的协方差矩阵,表示信号可信程度;
在步骤43,根据所述矩阵h、Ga、Ksi计算Za,假设x、y、d1、d2相互独立,利用加权线性最小二乘求解Za,
Za=(GaTKsi-1Ga)-1GaTKsi-1h (5)
进而得到移动台位置的初步估计解。
由于实际中信号传播时间、信号传播时间差测量误差的存在,以及系统有可能提供冗余的测量信息,使得从公式(5)中得到的移动台位置数据精度不高。因此图3中所述的流程的第二部分,是步骤44到步骤45,再次采用加权线性最小二乘法,降低x、y、d1、d2之间的相关性的影响。
在步骤44,构造矩阵h′、Ga、Ksi′,确定下述最小二乘形式:
Δ′=h′-Ga′Za′ (6)
最小二乘的目的是使‖Δ′‖2最小。
所述x0、y0实际中可以使用Za中求得的Za,1,Za,2来近似求解,r1 0用 来近似。
在步骤45,根据所述矩阵h′、Ga′、Ksi′计算Za′,也就是利用x,y,d1,d2间的联系以及Za的协方差矩阵,用最小二乘估计来求解 具体的计算公式如下:
Za′=(Ga′TKsi′-1Ga′)-1Ga′TKsi′-1h′ (8)
进而得到移动台位置的精确位置解。
图3中所述的流程的第三部分,是步骤46到步骤47,从上述式(8)得到的解中,挑选出移动台的最终位置估计值。
在步骤46,构造矩阵
然后按照下述方法计算Zp′:
进而得到移动台的最终位置解;
由于公式(9)得到的最终移动台位置的解的形式有四个,而其中只有一个是需要的位置估计值,所以,本实施例采用的策略为:根据移动台位置的初步估计解的距离差,从上述移动台的最终位置解中,挑选对应距离差最小的解作为移动台的最终位置估计值。
图4是本发明所述装置的实施例框图。图4描述的利用信号达到时间估计移动台位置的装置204,包括:
信号传播时间和信号传播时间差生成器201:用于接收所有与移动台保持连接的基站的信号传播时间与参与移动台定位的其它基站的相对于主基站的信号传播时间差相对应的测量值,所述测量值包括:所有与移动台保持连接的基站的RTT测量值、移动台的收发时间差测量值,以及相对于主基站的其它基站的信号到达时间差测量值和信号发射时间差测量值,将上述测量值按下述方法转换为与移动台保持连接的基站信号传播时间值和其它基站相对于主基站的信号传播时间差值:
与移动台保持连接的基站的信号传播时间值=(RTT测量值-移动台的收发时间差测量值)/2;
相对主基站的其它基站的信号传播时间差值=相应的信号到达时间差测量值-相应的信号发射时间差测量值;
信息数据库202:用于存储参与移动台定位的所有基站的几何位置和有关的统计信息;
位置估计器203:用于所述利用信号传播时间和信号传播时间差生成器输出的与移动台保持连接的一组信号传播时间和相对于主基站的一组信号传播时间差值和所述信息数据库提供的基站的几何位置以及有关的统计信息,进行移动台的位置估计,确定移动台的最终位置估计值。
首先信号传播时间和信号传播时间差生成器201选取与移动台保持连接的所有基站的信号传播时间的有关测量,信号传播时间和信号传播时间差生成器201将所得的信号传播时间有关测量值转换成相应的信号传播时间值,然后选取最小信号传播时间值所对应的基站作为主基站;接着发起以所选主基站为参考基站的信号传播时间差测量,信号传播时间和信号传播时间差生成器201将所得的信号传播时间差有关测量值转换成相应的信号传播时间差值;信息数据库202用于提供相应的基站坐标,以及信号传播时间和信号传播时间差的有关统计信息,由位置估计器203采用多个信号传播时间加多个信号传播时间差算法进行位置估计,最后位置估计器(203)给出移动台位置估计值,
所述位置估计器203包括:
存储器301:用于接收包括与移动台连接的所有基站的信号传播时间值、相对主基站的其它基站的信号传播时间差值、基站坐标、信号传播时间以及信号传播时间差的统计信息,进行移动台位置估计时所需的数据,以及存储移动台位置估计时产生的中间数据和最终数据;
处理器302:用于根据从所述存储器得到的进行移动台位置估计所需的数据,利用最小二乘法进行移动台位置的估计。
移动台位置估计的具体运算时,处理器302可以按照式(1)到式(9)的公式和要求运算得到。预算涉及到的中间数据和最终数据存储在存储器301中。
下面是本发明所述装置的一个具体应用的实例。
根据3GPP(第三代伙伴工程)的有关协议,WCDMA(宽带码分多址)系统中,基本的定位测量值为有关基站的RTT测量值、移动台的收发时间差测量值、相对主基站的一组信号到达时间差测量值,以及表示不同基站之间的信号发射时间差关系的RTD测量,其中RTT的测量必须是针对与移动台保持了连接的基站。此时采用发明的方法进行移动台位置估计的过程如以下描述:
首先信号传播时间和信号传播时间差生成器201将测量值转换成相应的信号传播时间和信号传播时间差值,具体的转换公式为:
保持连接的基站信号传播时间值=(基站的RTT测量值-对应UE的收发时间差测量值)/2;
相对主基站的信号传播时间差值=相应的信号到达时间差测量值-相应的信号发射时间差测量值;
其次将多个基站信号传播时间值以及多个信号传播时间差值输入位置估计器203。然后信息数据库202提供与信号传播时间、信号传播时间差相对应的基站坐标以及统计信息,最后位置估计器203根据得到的数据,利用多个信号传播时间和多个信号传播时间差等信息进行位置估计,进而得出移动台位置。
Claims (11)
1、一种利用信号达到时间和时间差估计移动台位置的估计方法,包括:
(1)确定参与进行移动台定位的基站,发起所有与移动台保持连接的基站的信号到达时间的测量,获得有关基站信号的往返时间测量值、移动台的信号收发时间差测量值;
(2)将上述步骤(1)得到的一组信号到达时间的测量值转换为信号从基站到移动台的传播时间,然后选取最小传播时间对应的基站或信号质量最好的基站为主基站;
(3)发起与所选主基站为参考基站的信号到达时间差测量,获得相对参考基站的一组信号到达时间差测量值,以及其它基站相对于主基站的信号发射时间差测量值,并将所得的上述两个测量值转换为信号从两个不同基站到同一个移动台的传播时间差;
(4)获取参与移动台定位的基站的坐标以及所有信号从基站到移动台的传播时间和信号从两个不同基站到同一个移动台的传播时间差的有关统计信息;
(5)利用上述步骤(4)提供的基站坐标以及所有信号从基站到移动台的传播时间、信号从两个不同基站到同一个移动台的传播时间差及其相关的统计信息进行移动台的位置估计,确定移动台的最终位置估计值。
2、根据权利要求1所述的移动台位置的估计方法,其特征在于所述步骤(2)和步骤(3)中的信号从基站到移动台的传播时间和信号从两个不同基站到同一个移动台的传播时间差的获得方法如下:
信号从基站到移动台的传播时间=(基站信号的往返时间测量值-移动台的信号收发时间差测量值)/2;
信号从两个不同基站到同一个移动台的传播时间差=相应的信号到达时间差测量值-其它基站相对于主基站的信号发射时间差测量值。
3、根据权利要求1所述的移动台位置的估计方法,其特征在于:所述步骤(5)进一步包括下述步骤:
(31)确定利用信号从基站到移动台的传播时间以及信号从两个不同基站到同一个移动台的传播时间差进行移动台位置估计的基本方程如下:
式中:i 1,2,.m,j 2,.m,而m>3,i等于1的基站为移动台所属的主基站;r1=c×τi,τi为第i个基站的信号到移动台的传播时间,rj1=c×τj1,τj1为第j个基站相对主基站的信号传播时间差,c为信号传播的速度;xi、yi为第i个基站的几何位置坐标。x、y为待估计的移动台的几何位置坐标。
(32)根据上述步骤(31)的基本方程,使用加权线性最小二乘法估计移动台的粗略位置,得到移动台位置的初步估计解;
(33)对上述步骤(32)得到的移动台的位置估计值再次使用加权线性最小二乘法进行优化,降低所述位置估计值之间的相关性影响;
(34)根据上述步骤(33)得到的移动台的优化位置估计值,确定移动台最终的位置估计值。
4、根据权利要求3所述的移动台位置的估计方法,其特征在于:所述步骤(32)中使用加权线性最小二乘法估计移动台的粗略位置进一步包括:
(41)对上述步骤(31)的基本方程进行变换,得到下述公式:
r1-d1=0
…
…
式中:i 2,.m,而m>3,第1个基站为主基站;ri=c×τi,τi为第i个基站的信号到移动台的传播时间,c为信号传播的速度;ri1=c×τi1,τi1为第i个基站信号相对第1个基站信号的传播时间差;xi、yi为第i个基站的几何位置坐标;
d2=(x-x1)2+(y-y1)2;x、y为待估计的移动台的几何位置坐标;
(42)构造矩阵h、Ga、Ksi,将上述步骤(41)中的公式表述为下述最小二乘的形式:Δ=h-GaZa;
式中:Δ为残差;Ksi为进行最小二乘估计时的加权系数矩阵,
Ksi=E(ΔΔT)=c2BQB
其中:B=diag{1,r1,r2,...,rm,2(r21+r1),2(r31+r1),...,2(rm1+r1)},diag表示对角阵,Q为信号到移动台的传播时间以及信号从两个不同基站到同一个移动台的传播时间差的误差的协方差矩阵,表示信号可信程度;
(43)根据所述矩阵h、Ga、Ksi计算Za,得到移动台位置的初步估计解。
5、根据权利要求4所述的移动台位置的估计方法,其特征在于:所述协方差矩阵Q配置为对角阵,对角线上元素为信号到移动台的传播时间以及信号从两个不同基站到同一个移动台的传播时间差的误差的方差。
6、根据权利要求4所述的移动台位置的估计方法,其特征在于:在所述步骤(43)按下述方法计算Za:
Za=(GaTKsi-1Ga)-1GaTKsi-1h。
7、根据权利要求6所述的移动台位置的估计方法,其特征在于:所述步骤(33)进一步包括:
(71)构造矩阵h′、Ga′、Ksi′,确定下述最小二乘形式:
Δ′=h′-Ga′Za′
式中Δ′为残差,Ksi′=E(Δ′Δ′T),并且,Ksi′=B′(GaTKsi-1Ga)-1B′,其中
(72)根据所述矩阵h′、Ga′、Ksi′计算Za′,得到移动台位置的精确位置解。
在所述步骤(71)中,x0,y0可以使用Za中求得的Za,1,Za,2来近似求解,r1 0用
来近似,所述Za为四行一列矩阵,所述Za,1为第一行第一列的元素,所述Za,2为第二行第一列的元素,所述Za,3为第三行第一列的元素,所述Za,4为第四行第一列的元素。
8、根据权利要求7所述的移动台位置的估计方法,其特征在于:在所述步骤(72)中按下述方法计算Za:
Za′=(Ga′TKsi′-1Ga′)-1Ga′TKsi′-1h′。
9、根据权利要求8所述的移动台位置的估计方法,其特征在于:所述步骤(34)进一步包括以下步骤:
(101)按照下述方法构造矩阵
得到移动台最终位置解;
(102)根据移动台位置的初步估计解与所述最终位置解的差值,从上述步骤(101)中的移动台的最终位置解中,挑选对应移动台位置的初步估计解与所述最终位置解的差值最小的解作为移动台的最终位置估计值。
10、一种利用信号达到时间和时间差估计移动台位置的装置,包括:
信号传播时间和信号传播时间差生成器:用于接收所有与移动台保持连接的基站的信号传播时间与参与移动台定位的其它基站的相对于主基站的信号到达时间差的测量值,所述主基站为最小传播时间对应的基站或信号质量最好的基站,所述测量值包括:所有与移动台保持连接的基站的往返时间测量值、移动台的信号收发时间差测量值,以及相对于主基站的其它基站的信号到达时间差测量值和信号发射时间差测量值,将上述测量值按下述方法转换为与移动台保持连接的基站到移动台的信号传播时间值和其它基站相对于主基站的信号传播时间差值:
基站到移动台的信号传播时间=(基站信号的往返时间测量值-移动台信号收发时间差测量值)/2;
相对主基站的其它基站的信号传播时间差值=相应的信号到达时间差测量值-信号发射时间差测量值;
信息数据库:用于存储参与移动台定位的所有基站的几何位置和有关的统计信息;
位置估计器:用于所述利用信号传播时间和信号传播时间差生成器输出的与移动台保持连接的一组信号传播时间和相对于主基站的一组信号传播时间差值和所述信息数据库提供的基站的几何位置以及有关的统计信息,进行移动台的位置估计,确定移动台的最终位置估计值。
11、根据权利要求10所述的估计移动台位置的装置,其特征在于:所述位置估计器包括:
存储器:用于接收包括与移动台保持连接的所有基站到移动台的信号传播时间、相对主基站的其它基站的信号传播时间差、基站坐标、所有信号传播时间以及信号传播时间差的统计信息,进行移动台位置估计时所需的数据,以及存储移动台位置估计时产生的中间数据和最终数据;
处理器:用于根据从所述存储器得到的进行移动台位置估计所需的一组信号传播时间数据和一组信号传播时间差数据,利用最小二乘法进行移动台位置的估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 01136462 CN1202688C (zh) | 2001-10-18 | 2001-10-18 | 利用信号达到时间和时间差估计移动台位置的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 01136462 CN1202688C (zh) | 2001-10-18 | 2001-10-18 | 利用信号达到时间和时间差估计移动台位置的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1413058A CN1413058A (zh) | 2003-04-23 |
CN1202688C true CN1202688C (zh) | 2005-05-18 |
Family
ID=4673679
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 01136462 Expired - Fee Related CN1202688C (zh) | 2001-10-18 | 2001-10-18 | 利用信号达到时间和时间差估计移动台位置的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1202688C (zh) |
Families Citing this family (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8483717B2 (en) | 2003-06-27 | 2013-07-09 | Qualcomm Incorporated | Local area network assisted positioning |
US8971913B2 (en) | 2003-06-27 | 2015-03-03 | Qualcomm Incorporated | Method and apparatus for wireless network hybrid positioning |
CN100428821C (zh) * | 2003-09-04 | 2008-10-22 | 华为技术有限公司 | 一种抑制到达时间附加时延误差的方法 |
US7319878B2 (en) | 2004-06-18 | 2008-01-15 | Qualcomm Incorporated | Method and apparatus for determining location of a base station using a plurality of mobile stations in a wireless mobile network |
CN100396021C (zh) * | 2004-09-30 | 2008-06-18 | 华为技术有限公司 | 避免级联组网往返时间测量误差累积的方法 |
CN100415039C (zh) * | 2004-10-26 | 2008-08-27 | 华为技术有限公司 | 基于gsm系统的wcdma系统中移动台定位的实现方法 |
CN100433897C (zh) * | 2004-12-27 | 2008-11-12 | 华为技术有限公司 | 往返时间的测量方法 |
CN100344207C (zh) * | 2005-04-29 | 2007-10-17 | 重庆邮电学院 | 一种单基站定位移动终端的方法 |
US7257413B2 (en) | 2005-08-24 | 2007-08-14 | Qualcomm Incorporated | Dynamic location almanac for wireless base stations |
RU2390791C2 (ru) | 2005-11-07 | 2010-05-27 | Квэлкомм Инкорпорейтед | Позиционирование для wlan и других беспроводных сетей |
CN101136678B (zh) * | 2006-08-29 | 2011-11-30 | 中国电信股份有限公司 | 一种用于移动台定位业务的定位数据处理方法及装置 |
US9226257B2 (en) | 2006-11-04 | 2015-12-29 | Qualcomm Incorporated | Positioning for WLANs and other wireless networks |
CN101075999B (zh) * | 2007-05-31 | 2011-06-29 | 西安电子科技大学 | 室内ofdm系统toa训练符号构建及toa估计方法与装置 |
CN101339238B (zh) * | 2007-07-06 | 2011-08-24 | 展讯通信(上海)有限公司 | 基于多系统联合测量的无线定位的位置计算方法与装置 |
CN101360322B (zh) * | 2007-07-30 | 2011-12-14 | 展讯通信(上海)有限公司 | 一种基于测量网络的无线定位系统与方法 |
CN101374249B (zh) * | 2007-08-21 | 2011-08-17 | 彭亮 | 一种基于移动终端ip地址获取该终端的实时精确地理位置的方法 |
CN101383855B (zh) * | 2007-09-04 | 2013-01-23 | 沈阳 | 一种基于ip地址信息获取计算机精确地理位置信息的方法 |
US7667649B2 (en) * | 2007-12-10 | 2010-02-23 | Trueposition, Inc. | Detection of time of arrival of CDMA signals in a wireless location system |
CN101848414B (zh) * | 2009-03-24 | 2013-04-17 | 电信科学技术研究院 | 确定信号传播时延的方法、对移动台定位的方法 |
CN101583078B (zh) * | 2009-06-15 | 2011-09-21 | 北京邮电大学 | 一种蜂窝定位方法 |
CN101931858B (zh) * | 2009-06-18 | 2013-01-30 | 华为技术有限公司 | 用户设备定位辅助数据配置的方法及装置 |
EP3650882A1 (en) * | 2009-08-13 | 2020-05-13 | InterDigital Patent Holdings, Inc. | Method and apparatus for supporting positioning measurements |
CN102026368A (zh) | 2009-09-29 | 2011-04-20 | 大唐移动通信设备有限公司 | 一种小区位置信息获取方法、系统及装置 |
US9699052B2 (en) * | 2013-05-30 | 2017-07-04 | Qualcomm Incorporated | Methods and systems for enhanced round trip time (RTT) exchange |
CN103983275B (zh) * | 2014-04-28 | 2016-11-09 | 上海融军科技有限公司 | 双基准信号源飞行器方向标定方法 |
US9907044B2 (en) * | 2014-09-15 | 2018-02-27 | Qualcomm Incorporated | IEEE 802.11 enhancements for high efficiency positioning |
CN105992335A (zh) * | 2015-01-30 | 2016-10-05 | 中兴通讯股份有限公司 | 一种定位方法、移动终端、定位装置及通信系统 |
EP3322998B1 (en) * | 2015-09-14 | 2020-02-26 | Red Point Positioning Corporation | Method to estimate and compensate for nlos bias in time difference of arrival estimate |
CN105548962B (zh) * | 2015-12-21 | 2018-05-22 | 李奇 | 非同步基站系统下的无线测距定位方法及系统 |
CN105738655A (zh) * | 2015-12-23 | 2016-07-06 | 深圳市金立通信设备有限公司 | 一种方向传感器校准方法及终端 |
CN105607039B (zh) * | 2016-01-27 | 2017-11-17 | 宁波大学 | 非视距环境下基于到达时间差的稳健最小二乘定位方法 |
CN106405496A (zh) * | 2016-08-31 | 2017-02-15 | 北斗时空信息技术(北京)有限公司 | 一种基于tdoa的室内定位方法 |
CN106932757A (zh) * | 2017-03-31 | 2017-07-07 | 南京理工大学 | 一种基于tdoa和多普勒的月球车联合定位方法 |
CN107450049A (zh) * | 2017-07-10 | 2017-12-08 | 武汉米风通信技术有限公司 | 到达定位系统时间差的精准计算方法 |
CN108093474B (zh) * | 2017-12-27 | 2020-09-29 | 重庆脉实智能制造有限公司 | 利用虚拟时间同步实现室内定位的方法及系统 |
CN111225439B (zh) * | 2018-11-23 | 2022-11-22 | 中兴通讯股份有限公司 | 确定终端位置的方法、装置及存储介质 |
CN111770568B (zh) * | 2019-04-02 | 2022-01-04 | 大唐移动通信设备有限公司 | 一种定位测量值的确定方法及装置 |
CN110988797A (zh) * | 2019-11-19 | 2020-04-10 | 北京智芯微电子科技有限公司 | 定位系统与定位方法 |
CN114222362B (zh) * | 2021-12-15 | 2024-03-15 | 武汉乾阳智能科技有限公司 | 定位方法及定位装置 |
-
2001
- 2001-10-18 CN CN 01136462 patent/CN1202688C/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN1413058A (zh) | 2003-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1202688C (zh) | 利用信号达到时间和时间差估计移动台位置的方法及装置 | |
CN1185897C (zh) | 一种移动台位置的估计方法及装置 | |
CN100347563C (zh) | 带有中继器检测和补偿的无线通信系统中的位置确定 | |
CN1852591A (zh) | 一种定位移动终端的方法 | |
CN1163091C (zh) | 无线电基站发射定时偏移的测定方法和偏移测定装置 | |
CN1130940C (zh) | 距离检测方法及其装置 | |
CN1150696C (zh) | 用于调整基站计时的系统和方法 | |
CN1225920C (zh) | 定位系统,定位服务器,基站和终端位置估算方法 | |
CN1259785C (zh) | 获取时分同步cdma(td-scdma)用户终端的载波频偏的方法 | |
CN1722905A (zh) | 位置信息提供系统、基站和随其使用的位置信息提供方法 | |
CN1248010C (zh) | 全球定位系统终端和定位测量系统 | |
CN1854752A (zh) | 位置估算方法与系统 | |
CN1616986A (zh) | Gps定位系统以及对gps终端进行定位的方法 | |
CN1751248A (zh) | 无线移动通信网络内使用移动站确定基站位置参数 | |
CN101049041A (zh) | 根据属于多个网络的基站的有关信息对无线装置进行定位 | |
CN1909732A (zh) | 使用混合网络的越区切换系统与方法 | |
CN1483148A (zh) | 适用于混合定位系统确定误差估计的方法和装置 | |
CN101076178A (zh) | 一种邻区规划方法和装置 | |
CN1752769A (zh) | 定位系统、定位装置、定位信息提供装置以及定位方法 | |
CN1679253A (zh) | 使用定位信息的移动通信系统 | |
CN1849007A (zh) | 一种移动台定位方法 | |
CN1628480A (zh) | 用于包括中继器的无线网络的信号路径检测 | |
CN1235429C (zh) | 一种位置估计方法 | |
CN1753550A (zh) | 利用波束成形增强定位信号发送的方法及设备 | |
CN1913703A (zh) | 多扇区小区对rtt测量值进行校正的方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
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: 20050518 Termination date: 20161018 |