CN1413058A - 利用信号达到时间和时间差估计移动台位置的方法及装置 - Google Patents

利用信号达到时间和时间差估计移动台位置的方法及装置 Download PDF

Info

Publication number
CN1413058A
CN1413058A CN 01136462 CN01136462A CN1413058A CN 1413058 A CN1413058 A CN 1413058A CN 01136462 CN01136462 CN 01136462 CN 01136462 A CN01136462 A CN 01136462A CN 1413058 A CN1413058 A CN 1413058A
Authority
CN
China
Prior art keywords
tdoa
base station
travelling carriage
toa
location
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
CN 01136462
Other languages
English (en)
Other versions
CN1202688C (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.)
Huawei Technologies Co Ltd
Original Assignee
Huawei Technologies Co Ltd
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 Huawei Technologies Co Ltd filed Critical Huawei Technologies Co Ltd
Priority to CN 01136462 priority Critical patent/CN1202688C/zh
Publication of CN1413058A publication Critical patent/CN1413058A/zh
Application granted granted Critical
Publication of CN1202688C publication Critical patent/CN1202688C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

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测量信息利用上,而仅仅只是对部分信息的利用,使得移动台的定位精度并不能达到理想的要求。
发明内容
本发明的目的在于,提供一种利用TOA和TDOA、精度较高的估计移动台位置的方法及装置。
为达到上述目的,本发明提供的利用信号达到时间估计移动台位置的方法,包括:
(1)确定参与进行移动台定位的基站,发起所有与移动台保持连接的基站的TOA(单个信号到达时间)的测量,获得有关基站的RTT(往返时间:Round Trip Time)测量值、移动台的Rx(接收时间)-Tx(发送时间)测量值;
(2)将上述步骤(1)得到的一组TOA的测量值转换为相应的TOA值,然后选取最小TOA值对应的基站为主基站;
(3)发起与所选主基站为参考基站的TDOA测量,获得相对参考基站的一组TDOA测量值,以及RTD(相对时间差)测量值,并将所得的TDOA测量值转换为相应的TDOA值;
(4)获取参与移动台定位的基站的坐标以及所有TOA和TDOA的有关统计信息;
(5)利用上述步骤(4)提供的基站坐标以及所有TOA、TDOA及其相关的统计信息进行移动台的位置估计,确定移动台的最终位置估计值。上述步骤(5)进一步包括下述步骤:
(31)确定利用TOA以及TDOA进行移动台位置估计的基本方程如下:
             (x-xi)2+(y-yi)2=(ri)2 ( ( x - x j ) 2 + ( y - y j 2 ) - ( x - x 1 ) 2 + ( y - y 1 2 ) ) 2 = r j 1 2
式中:i=1,2,.m,j=2,.m,而m>=3,i等于1的基站为移动台所属的主基站;ri=c×τi,τi为第i个基站的信号到达时间TOA,ri1=c×τi1,τi1为第j个基站相对主基站的信号到达时间差TDOA,c为信号传播的速度(光速);xi、yi为第i个基站的几何位置坐标。x、y为待估计的移动台的几何位置坐标。
(32)根据上述步骤(31)的基本方程,使用加权线性最小二乘法估计移动台的粗略位置,得到移动台位置的初步估计解;
(33)对上述步骤(32)得到的移动台的位置估计值再次使用加权线性最小二乘法进行优化,降低所述位置估计值之间的相关性影响;
(34)根据上述步骤(33)得到的移动台的优化位置估计值,确定移动台最终的位置估计值。
其中所述步骤(32)中使用加权线性最小二乘法估计移动台的粗略位置进一步包括:
(41)对上述步骤(31)的基本方程进行变换,得到下述公式:
 r1-d1=0 r 1 2 - d 2 = 0 r 2 2 + x 1 2 - x 2 2 + y 1 2 - y 2 2 - 2 ( x 1 - x 2 ) x - 2 ( y 1 - y 2 ) y - d 2 = 0 r 3 2 + x 1 2 - x 3 2 + y 1 2 - y 3 2 - 2 ( x 1 - x 3 ) x - 2 ( y 1 - y 3 ) y - d 2 = 0 r m 2 + x 1 2 - x m 2 + y 1 2 - y m 2 - 2 ( x 1 - x m ) x - 2 ( y 1 - y m ) y - d 2 = 0 r 21 2 + x 1 2 - x 2 2 + y 1 2 - y 2 2 - 2 ( x 1 - x 2 ) x - 2 ( y 1 - y 2 ) y + 2 r 21 d 1 = 0 r 31 2 + x 1 2 - x 3 2 + y 1 2 - y 3 2 - 2 ( x 1 - x 3 ) x - 2 ( y 1 - y 3 ) y + 2 r 31 d 1 = 0 r m 1 2 + x 1 2 - x m 2 + y 1 2 - y m 2 - 2 ( x 1 - x m ) x - 2 ( y 1 - y m ) y + 2 r m 1 d 1 = 0
式中:i=2,.m,而m>=3,第1个基站为主基站;ri=c×τi,τi为第i个基站的信号到达时间TOA,c为信号传播的速度;ri1=c×τi1,τi1为第i个基站信号相对第1个基站信号的到达时间差TDOA;xi、yi为第i个基站的几何位置坐标; d 1 = ( x - x 1 ) 2 + ( y - y 1 ) 2 , 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为TOA以及TDOA的误差的协方差矩阵,表示信号可信程度;
式中:Δ为残差;
Figure A0113646200121
(43)根据所述矩阵h、Ga、Ksi计算Za,得到移动台位置的初步估计解。
在所述步骤(43)中,按下述方法计算Za:
           Za=(GaTKsi-1Ga)-1GaTKsi-1h。
所述步骤(33)进一步包括:
(71)构造矩阵h′、Ga′、Ksi′,确定下述最小二乘形式:
              Δ′=h′-Ga′Za′
式中Δ′为残差,Ksi′=E(Δ′Δ′T),并且,
Figure A0113646200123
其中 B ′ = diag { 2 ( x 0 - x 1 ) , 2 ( y 0 - y 1 ) , 2 r 1 0 , 1 } ; Ga ′ = 1 0 0 1 1 1 1 1 , h ′ = ( Z a , 1 - x 1 ) 2 ( Z a , 2 - y 1 ) 2 Z a , 3 2 Z a , 4 , Za ′ = ( x - x 1 ) 2 ( y - y 1 ) 2 ;
(72)根据所述矩阵h′、Ga′、Ksi′计算Za′,得到移动台位置的精确位置解。
在所述步骤(72)中,按下述方法计算Za′:
所述步骤(34)进一步包括以下步骤:
(101)按照下述方法计算Zp′: Zp ′ = ± Za ′ + x 1 y 1 ,
其中:
Figure A0113646200133
得到移动台最终位置解;
(102)根据移动台位置的初步估计解的距离差,从上述步骤(101)中的移动台的最终位置解中,挑选对应距离差最小的解作为移动台的最终位置估计值。
本发明同时还提供了一种利用信号达到时间估计移动台位置的装置,包括:
TOA和TDOA生成器:用于接收所有与移动台保持连接的基站的TOA与参与移动台定位的其它基站的相对于主基站的TDOA相对应的测量值,即所有与移动台保持连接的基站的RTT测量值、移动台的Rx和Tx测量值,以及相对于主基站的其它基站的TDOA测量值和RTD测量值,将上述测量值按下述方法转换为与移动台保持连接的基站TOA值和其它基站相对于主基站的TDOA值:
与移动台保持连接的基站的TOA值=(RTT测量值-移动台的Rx-Tx测量值)/2;
相对主基站的其它基站的TDOA值=相应的TDOA测量值-RTD测量值;
信息数据库:用于存储参与移动台定位的所有基站的几何位置和有关的统计信息;
位置估计器:用于所述利用TOA和TDOA生成器输出的与移动台保持连接的一组TOA和相对于主基站的一组TDOA值和所述信息数据库提供的基站的几何位置以及有关的统计信息,进行移动台的位置估计,确定移动台的最终位置估计值。
所述位置估计器包括:
存储器:用于接收包括与移动台保持连接的所有基站TOA值、相对主基站的其它基站的TDOA值、基站坐标、所有TOA以及TDOA的统计信息,进行移动台位置估计时所需的数据,以及存储移动台位置估计时产生的中间数据和最终数据;
处理器:用于根据从所述存储器得到的进行移动台位置估计所需的一组TOA数据和一组TDOA数据,利用最小二乘法进行移动台位置的估计。
由于本发明能够综合利用与移动台保持连接的多个基站的一组TOA和相对主基站的其它基站的一组TDOA,利用最小二乘法进行移动台位置估计,更适合于移动台处于切换状态下的定位,这样,本发明在移动台的位置估计过程中充分利用了全部的TOA和TDOA信息,因此更能综合TDOA与TOA方法的优点,使移动台的位置估计精度以及GDOP性能进一步得到提高;同时本发明进行位置估计过程中无须迭代求解,运算量小,位置估计所需时间也较少,因此采用本发明可以使同一时间段内的系统处理的定位请求数量明显增加。
附图说明
图1是本发明所述方法进行移动台位置估计的几何原理图;
图2是本发明所述方法的实施例流程图;
图3是图2中步骤5的具体实施例流程图;
图4是本发明所述装置的实施例框图;
图5是图4所述的位置估计器的实施例框图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细的描述。
发明提出的方法和装置适用于3个或3个以上的基站参与对移动台进行位置估计。图1给出了综合利用与移动台保持连接的多个基站的一组TOA以及相对于主基站的其它基站的TDOA进行移动台位置估计的几何原理图。图1中,不同圆和双曲线的交点即为移动台的位置,当TOA以及TDOA存在误差时,曲线将不交于一点,而是大致区域。
图2是本发明所述方法的实施例流程图。按照图2实施本发明共具有5个步骤,步骤1确定参与进行移动台定位的基站,发起所有与移动台保持连接的基站的TOA(单个信号到达时间)的测量,获得有关基站的RTT(往返时间:Round Trip Time)测量值、移动台的Rx(接收时间)-Tx(发送时间)测量值;在步骤2将上述步骤1得到的一组TOA的测量值转换为相应的TOA值,然后选取最小TOA值对应的基站为主基站;在步骤3发起与所选主基站为参考基站的TDOA测量,获得相对参考基站的一组TDOA测量值,以及RTD(相对时间差)测量值,并将所得的TDOA测量值转换为相应的TDOA值;在步骤4获取参与移动台定位的基站的坐标以及所有TOA和TDOA的有关统计信息;最后在步骤5利用上述步骤4得到的基站坐标以及所有TOA、TDOA及其相关的统计信息进行移动台的位置估计,确定移动台的最终位置估计值。
上述步骤5的具体实施过程参考图3。图3中所述的流程主要包括三部分内容,第一部分是步骤41到步骤43,采用加权线性最小二乘法得到移动台位置的初步估计解。
首先在步骤41确定利用TOA以及TDOA进行移动台位置估计的基本方程如下:
            (x-xi)2+(y-yi)2=(ri)2             (1) ( ( x - x j ) 2 + ( y - y j 2 ) - ( x - x 1 ) 2 + ( y - y 1 2 ) ) 2 = r j 1 2
式中:i=1,2,.m,j=2,.m,而m>=3,i等于1的基站为移动台所属的主基站;ri=c×τi,τi为第i个基站的信号到达时间TOA,rj1=c×τj1,τi1为第j个基站相对主基站的信号到达时间差TDOA,c为信号传播的速度(光速);xi、yi为第i个基站的几何位置坐标。x、y为待估计的移动台的几何位置坐标。
在步骤42,首先对上述步骤41的基本方程进行变换,得到下述公式:
r1-d1=0 r 1 2 - d 2 = 0 r 2 2 + x 1 2 - x 2 2 + y 1 2 - y 2 2 - 2 ( x 1 - x 2 ) x - 2 ( y 1 - y 2 ) y - d 2 = 0 r 3 2 + x 1 2 - x 3 2 + y 1 2 - y 3 2 - 2 ( x 1 - x 3 ) x - 2 ( y 1 - y 3 ) y - d 2 = 0 r m 2 + x 1 2 - x m 2 + y 1 2 - y m 2 - 2 ( x 1 - x m ) x - 2 ( y 1 - y m ) y - d 2 = 0 r 21 2 + x 1 2 - x 2 2 + y 1 2 - y 2 2 - 2 ( x 1 - x 2 ) x - 2 ( y 1 - y 2 ) y + 2 r 21 d 1 = 0 r 31 2 + x 1 2 - x 3 2 + y 1 2 - y 3 2 - 2 ( x 1 - x 3 ) x - 2 ( y 1 - y 3 ) y + 2 r 31 d 1 = 0 r m 1 2 + x 1 2 - x m 2 + y 1 2 - y m 2 - 2 ( x 1 - x m ) x - 2 ( y 1 - y m ) y + 2 r m 1 d 1 = 0 - - - - ( 2 )
式中:i=2,.m,而m>=3,第1个基站为主基站;ri=c×τi,τi为第i个基站的信号到达时间TOA,c为信号传播的速度;ri1=c×τi1,τi1为第i个基站信号相对第1个基站信号的到达时间差TDOA;xi、yi为第i个基站的几何位置坐标; d 1 = ( x - x 1 ) 2 + ( y - y 1 ) 2 , d2=(x-x1)2+(y-y1)2;x、y为待估计的移动台的几何位置坐标。实际中,若TOA以及TDOA的数量不是一一对应时,可以根据相应的TOA与TDOA的关系将不足的TOA或TDOA补足。例如,如果缺少第3个基站的τ3,可利用τ3=τ311的关系将τ3补齐,使之满足公式(2)的形式。另外也可在构造矩阵h、Ga、Ksi时,将缺少的TOA或TDOA所对应的行去掉,其他不变化。
然后构造矩阵h、Ga、Ksi,将上述步骤41中的公式表述为下述最小二乘的形式:
         Δ=h-GaZa                          (3):
最小二乘的目的是使
Figure A0113646200173
最小。
式中:Δ即为残差;Ksi为进行最小二乘估计时的加权系数矩阵,
           Ksi=E(ΔΔT)=c2BQB              (4);
其中:B=diag{1,r1,r2,…,rm,2(r21+r1),2(r31+r1),…,2(rm1+r1)}diag表示对角阵,Q为TOA以及TDOA的误差的协方差矩阵,表示信号可信程度;
Figure A0113646200182
在步骤43,根据所述矩阵h、Ga、Ksi计算Za,假设x、y、d1、d2相互独立,利用加权线性最小二乘求解Za,
      Za=(GaTKsi-1Ga)-1GaTKsi-1h          (5)
进而得到移动台位置的初步估计解。
由于实际中TOA、TDOA测量误差的存在,以及系统有可能提供冗余的测量信息,使得从公式(5)中得到的移动台位置数据精度不高。因此图3中所述的流程的第二部分,是步骤44到步骤45,再次采用加权线性最小二乘法,降低x、y、d1、d2之间的相关性的影响。
在步骤44,构造矩阵h′、Ga′、 Ksi′,确定下述最小二乘形式:
             Δ′=h′-Ga′Za′             (6)
最小二乘的目的是使
Figure A0113646200183
最小。
式中Δ′为残差,Ksi′=E(Δ′Δ′T),并且, Ksi ′ = B ′ ( Ga T Ksi - 1 Ga ) - 1 B ′ , 其中 B ′ = diag { 2 ( x 0 - x 1 ) , 2 ( y 0 - y 1 ) , 2 r 1 0 , 1 } ; Ga ′ = 1 0 0 1 1 1 1 1 , h ′ = ( Z a , 1 - x 1 ) 2 ( Z a , 2 - y 1 ) 2 Z a , 3 2 Z a , 4 , Za ′ = ( x - x 1 ) 2 ( y - y 1 ) 2 ;
所述x0、y0实际中可以使用Za中求得的Za,1,Za,2来近似求解,r1 0 r 1 0 = ( Z a , 3 + Z a , 4 ) / 2 来近似。
在步骤45,根据所述矩阵h′、Ga′、 Ksi′计算Za′,也就是利用x,y,d1,d2间的联系以及Za的协方差矩阵,用最小二乘估计来求解 Za ′ = ( x - x 1 ) 2 ( y - y 1 ) 2 . 具体的计算公式如下:
       Za′=(Ga′TKsi′-1Ga′)-1Ga′TKsi′-1h′    (8)
进而得到移动台位置的精确位置解。
图3中所述的流程的第三部分,是步骤46到步骤47,从上述式(8)得到的解中,挑选出移动台的最终位置估计值。
在步骤46,构造矩阵 Zp ′ = x y , 然后按照下述方法计算Zp′: Zp ′ = ± Za ′ + x 1 y 1 - - - ( 9 )
进而得到移动台的最终位置解;
由于公式(9)得到的最终移动台位置的解的形式有四个,而其中只有一个是需要的位置估计值,所以,本实施例采用的策略为:根据移动台位置的初步估计解的距离差,从上述移动台的最终位置解中,挑选对应距离差最小的解作为移动台的最终位置估计值。
图4是本发明所述装置的实施例框图。图4描述的利用信号达到时间估计移动台位置的装置204,包括:
TOA和TDOA生成器201:用于接收所有与移动台保持连接的基站的TOA与参与移动台定位的其它基站的相对于主基站的TDOA相对应的测量值,即所有与移动台保持连接的基站的RTT测量值、移动台的Rx和Tx测量值,以及相对于主基站的其它基站的TDOA测量值和RTD测量值,将上述测量值按下述方法转换为与移动台保持连接的基站TOA值和其它基站相对于主基站的TDOA值:
与移动台保持连接的基站的TOA值=(RTT测量值-移动台的Rx-Tx测量值)/2;
相对主基站的其它基站的TDOA值=相应的TDOA测量值-RTD测量值;
信息数据库202:用于存储参与移动台定位的所有基站的几何位置和有关的统计信息;
位置估计器203:用于所述利用TOA和TDOA生成器输出的与移动台保持连接的一组TOA和相对于主基站的一组TDOA值和所述信息数据库提供的基站的几何位置以及有关的统计信息,进行移动台的位置估计,确定移动台的最终位置估计值。
首先TOA和TDOA生成器201选取与移动台保持连接的所有基站的TOA的有关测量,TOA和TDOA生成器201将所得的TOA有关测量值转换成相应的TOA值,然后选取最小TOA值所对应的基站作为主基站;接着发起以所选主基站为参考基站的TDOA测量,TOA和TDOA生成器201将所得的TDOA有关测量值转换成相应的TDOA值;信息数据库202用于提供相应的基站坐标,以及TOA和TDOA的有关统计信息,由位置估计器203采用多个TOA加多个TDOA算法进行位置估计,最后位置估计器(203)给出移动台位置估计值,
所述位置估计器203包括:
存储器301:用于接收包括与移动台连接的所有基站的TOA值、相对主基站的其它基站的TDOA值、基站坐标、TOA以及TDOA的统计信息,进行移动台位置估计时所需的数据,以及存储移动台位置估计时产生的中间数据和最终数据;
处理器302:用于根据从所述存储器得到的进行移动台位置估计所需的数据,利用最小二乘法进行移动台位置的估计。
移动台位置估计的具体运算时,处理器302可以按照式(1)到式(9)的公式和要求运算得到。预算涉及到的中间数据和最终数据存储在存储器301中。
下面是本发明所述装置的一个具体应用的实例。
根据3GPP2(第三代伙伴工程2)的有关协议,WCDMA(宽带码分多址)系统中,基本的定位测量值为有关基站的RTT测量值、移动台的Rx、Tx测量值、相对主基站的一组TDOA测量值,以及表示不同基站之间的时差关系的RTD测量,其中RTT的测量必须是针对与移动台保持了连接的基站。此时采用发明的方法进行移动台位置估计的过程如以下描述:
首先TOA和TDOA生成器201将测量值转换成相应的TOA和TDOA值,具体的转换公式为:
保持连接的基站TOA值=(基站的RTT测量值-对应UE的Rx-Tx测量值)/2;
相对主基站的TDOA值=相应的TDOA测量值-RTD测量;
其次将多个基站TOA值以及多个TDOA值输入位置估计器203。然后信息数据库202提供与TOA、TDOA相对应的基站坐标以及统计信息,最后位置估计器203根据得到的数据,利用多个TOA和多个TDOA等信息进行位置估计,进而得出移动台位置。

Claims (12)

1、一种利用信号达到时间和时间差估计移动台位置的估计方法,包括:
(1)确定参与进行移动台定位的基站,发起所有与移动台保持连接的基站的TOA(单个信号到达时间)的测量,获得有关基站的RTT(往返时间:Round Trip Time)测量值、移动台的Rx(接收时间)-Tx(发送时间)测量值;
(2)将上述步骤(1)得到的一组TOA的测量值转换为相应的TOA值,然后选取最小TOA值对应的基站或信号质量最好的基站为主基站;
(3)发起与所选主基站为参考基站的TDOA测量,获得相对参考基站的一组TDOA测量值,以及RTD(相对时间差)测量值,并将所得的TDOA测量值转换为相应的TDOA值;
(4)获取参与移动台定位的基站的坐标以及所有TOA和TDOA的有关统计信息;
(5)利用上述步骤(4)提供的基站坐标以及所有TOA、TDOA及其相关的统计信息进行移动台的位置估计,确定移动台的最终位置估计值。
2、根据权利要求1所述的移动台位置的估计方法,其特征在于所述步骤(2)和步骤(3)中的TOA值和TDOA值的获得方法如下:
与移动台保持连接的基站TOA值=(RTT测量值-移动台的Rx-Tx测量值)/2;
相对主基站的TDOA值=相应的TDOA测量值-RTD测量值。
3、根据权利要求1所述的移动台位置的估计方法,其特征在于:所述步骤(5)进一步包括下述步骤:
(31)确定利用TOA以及TDOA进行移动台位置估计的基本方程如下:
          (x-xi)2+(y-yi)2=(ri)2 ( ( x - x j ) 2 + ( y - y j 2 ) - ( x - x 1 ) 2 + ( y - y 1 2 ) ) 2 = r j 1 2
式中:i=1,2,.m,j=2,.m,而m>=3,i等于1的基站为移动台所属的主基站;ri=c×τi,τi为第i个基站的信号到达时间TOA,ri1=c×τi1,τi1为第j个基站相对主基站的信号到达时间差TDOA,c为信号传播的速度(光速);xi、yi为第i个基站的几何位置坐标。x、y为待估计的移动台的几何位置坐标。
(32)根据上述步骤(31)的基本方程,使用加权线性最小二乘法估计移动台的粗略位置,得到移动台位置的初步估计解;
(33)对上述步骤(32)得到的移动台的位置估计值再次使用加权线性最小二乘法进行优化,降低所述位置估计值之间的相关性影响;
(34)根据上述步骤(33)得到的移动台的优化位置估计值,确定移动台最终的位置估计值。
4、根据权利要求3所述的移动台位置的估计方法,其特征在于:所述步骤(32)中使用加权线性最小二乘法估计移动台的粗略位置进一步包括:
(41)对上述步骤(31)的基本方程进行变换,得到下述公式:
r1-d1=0 r 1 2 - d 2 = 0 r 2 2 + x 1 2 - x 2 2 + y 1 2 - y 2 2 - 2 ( x 1 - x 2 ) x - 2 ( y 1 - y 2 ) y - d 2 = 0 r 3 2 + x 1 2 - x 3 2 + y 1 2 - y 3 2 - 2 ( x 1 - x 3 ) x - 2 ( y 1 - y 3 ) y - d 2 = 0 r m 2 + x 1 2 - x m 2 + y 1 2 - y m 2 - 2 ( x 1 - x m ) x - 2 ( y 1 - y m ) y - d 2 = 0 r 21 2 + x 1 2 - x 2 2 + y 1 2 - y 2 2 - 2 ( x 1 - x 2 ) x - 2 ( y 1 - y 2 ) y + 2 r 21 d 1 = 0 r 31 2 + x 1 2 - x 3 2 + y 1 2 - y 3 2 - 2 ( x 1 - x 3 ) x - 2 ( y 1 - y 3 ) y + 2 r 31 d 1 = 0 r m 1 2 + x 1 2 - x m 2 + y 1 2 - y m 2 - 2 ( x 1 - x m ) x - 2 ( y 1 - y m ) y + 2 r m 1 d 1 = 0
式中:i=2,.m,而m>=3,第1个基站为主基站;ri=c×τi,τi为第i个基站的信号到达时间TOA,c为信号传播的速度;ri1=c×τi1,τi1为第i个基站信号相对第1个基站信号的到达时间差TDOA;xi、yi为第i个基站的几何位置坐标; d 1 = ( x - x 1 ) 2 + ( y - y 1 ) 2 , 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为TOA以及TDOA的误差的协方差矩阵,表示信号可信程度;
式中:Δ为残差;
Figure A0113646200046
Figure A0113646200051
(43)根据所述矩阵h、Ga、Ksi计算Za,得到移动台位置的初步估计解。
5、根据权利要求4所述的移动台位置的估计方法,其特征在于:所述协方差矩阵Q配置为对角阵,对角线上元素为TOA以及TDOA误差的方差。
6、根据权利要求4所述的移动台位置的估计方法,其特征在于:在所述步骤(43)按下述方法计算Za:
          Za=(GaTKsi-1Ga)-1GaTKsi-1h。
7、根据权利要求3所述的移动台位置的估计方法,其特征在于:所述步骤(33)进一步包括:
(71)构造矩阵h′、Ga′、Ksi′,确定下述最小二乘形式:
          Δ′=h′-Ga′Za′
式中Δ′为残差,Ksi′=E(Δ′Δ′T),并且,
Figure A0113646200052
,其中 B ′ = diag { 2 ( x 0 - x 1 ) , 2 ( y 0 - y 1 ) , 2 r 1 0 , 1 } ; Ga ′ = 1 0 0 1 1 1 1 1 , h ′ = ( Z a , 1 - x 1 ) 2 ( Z a , 2 - y 1 ) 2 Z a , 3 2 Z a , 4 , = Za ′ = ( x - x 1 ) 2 ( y - y 1 ) 2 ;
(72)根据所述矩阵h′、Ga′、Ksi′计算Za′,得到移动台位置的精确位置解。
8、根据权利要求7所述的移动台位置的估计方法,其特征在于:在所述步骤(71)中,x0,y0可以使用Za中求得的Za,1,Za,2来近似求解,r1 0 r 1 0 = ( Z a , 3 + Z a , 4 ) / 2 来近似。
9、根据权利要求7所述的移动台位置的估计方法,其特征在于:在所述步骤(72)中按下述方法计算Za′:
       Za′=(Ga′TKsi′-1Ga′)-1Ga′TKsi′-1h′。
10、根据权利要求3所述的移动台位置的估计方法,其特征在于:所述步骤(34)进一步包括以下步骤:
(101)按照下述方法计算Zp′: Zp ′ = ± Za ′ + x 1 y 1 ,
其中: Zp ′ = x y , 得到移动台最终位置解;
(102)根据移动台位置的初步估计解的距离差,从上述步骤(101)中的移动台的最终位置解中,挑选对应距离差最小的解作为移动台的最终位置估计值。
11、一种利用信号达到时间和时间差估计移动台位置的装置,包括:
TOA和TDOA生成器:用于接收所有与移动台保持连接的基站的TOA与参与移动台定位的其它基站的相对于主基站的TDOA相对应的测量值,即所有与移动台保持连接的基站的RTT测量值、移动台的Rx和Tx测量值,以及相对于主基站的其它基站的TDOA测量值和RTD测量值,将上述测量值按下述方法转换为与移动台保持连接的基站TOA值和其它基站相对于主基站的TDOA值:
与移动台保持连接的基站的TOA值=(RTT测量值-移动台的Rx-Tx测量值)/2;
相对主基站的其它基站的TDOA值=相应的TDOA测量值-RTD测量值;
信息数据库:用于存储参与移动台定位的所有基站的几何位置和有关的统计信息;
位置估计器:用于所述利用TOA和TDOA生成器输出的与移动台保持连接的一组TOA和相对于主基站的一组TDOA值和所述信息数据库提供的基站的几何位置以及有关的统计信息,进行移动台的位置估计,确定移动台的最终位置估计值。
12、根据权利要求1所述的估计移动台位置的装置,其特征在于:所述位置估计器包括:
存储器:用于接收包括与移动台保持连接的所有基站TOA值、相对主基站的其它基站的TDOA值、基站坐标、所有TOA以及TDOA的统计信息,进行移动台位置估计时所需的数据,以及存储移动台位置估计时产生的中间数据和最终数据;
处理器:用于根据从所述存储器得到的进行移动台位置估计所需的一组TOA数据和一组TDOA数据,利用最小二乘法进行移动台位置的估计。
CN 01136462 2001-10-18 2001-10-18 利用信号达到时间和时间差估计移动台位置的方法及装置 Expired - Fee Related CN1202688C (zh)

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 true CN1413058A (zh) 2003-04-23
CN1202688C 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)

Cited By (39)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100344207C (zh) * 2005-04-29 2007-10-17 重庆邮电学院 一种单基站定位移动终端的方法
CN100396021C (zh) * 2004-09-30 2008-06-18 华为技术有限公司 避免级联组网往返时间测量误差累积的方法
CN100415039C (zh) * 2004-10-26 2008-08-27 华为技术有限公司 基于gsm系统的wcdma系统中移动台定位的实现方法
CN100428821C (zh) * 2003-09-04 2008-10-22 华为技术有限公司 一种抑制到达时间附加时延误差的方法
CN100433897C (zh) * 2004-12-27 2008-11-12 华为技术有限公司 往返时间的测量方法
CN101848414A (zh) * 2009-03-24 2010-09-29 大唐移动通信设备有限公司 确定信号传播时延的方法、对移动台定位的方法
WO2010145559A1 (zh) * 2009-06-18 2010-12-23 华为技术有限公司 用户设备定位辅助数据配置的方法及装置
CN101075999B (zh) * 2007-05-31 2011-06-29 西安电子科技大学 室内ofdm系统toa训练符号构建及toa估计方法与装置
CN101374249B (zh) * 2007-08-21 2011-08-17 彭亮 一种基于移动终端ip地址获取该终端的实时精确地理位置的方法
CN101339238B (zh) * 2007-07-06 2011-08-24 展讯通信(上海)有限公司 基于多系统联合测量的无线定位的位置计算方法与装置
CN101583078B (zh) * 2009-06-15 2011-09-21 北京邮电大学 一种蜂窝定位方法
CN101136678B (zh) * 2006-08-29 2011-11-30 中国电信股份有限公司 一种用于移动台定位业务的定位数据处理方法及装置
CN101360322B (zh) * 2007-07-30 2011-12-14 展讯通信(上海)有限公司 一种基于测量网络的无线定位系统与方法
US8320934B2 (en) 2005-08-24 2012-11-27 Qualcomm Incorporated Dynamic location almanac for wireless base stations
CN101383855B (zh) * 2007-09-04 2013-01-23 沈阳 一种基于ip地址信息获取计算机精确地理位置信息的方法
CN101821645B (zh) * 2007-12-10 2013-06-19 真实定位公司 无线定位系统中的cdma信号的到达时间检测
US8483717B2 (en) 2003-06-27 2013-07-09 Qualcomm Incorporated Local area network assisted positioning
CN103983275A (zh) * 2014-04-28 2014-08-13 上海融军科技有限公司 双基站飞行器方向标定方法
US8971913B2 (en) 2003-06-27 2015-03-03 Qualcomm Incorporated Method and apparatus for wireless network hybrid positioning
US9042917B2 (en) 2005-11-07 2015-05-26 Qualcomm Incorporated Positioning for WLANS and other wireless networks
CN104869590A (zh) * 2009-08-13 2015-08-26 交互数字专利控股公司 支持定位测量的方法、本地服务器、wtru及其执行的方法
USRE45808E1 (en) 2004-06-18 2015-11-17 Qualcomm Incorporated Method and apparatus for determining location of a base station using a plurality of mobile stations in a wireless mobile network
US9226257B2 (en) 2006-11-04 2015-12-29 Qualcomm Incorporated Positioning for WLANs and other wireless networks
CN105548962A (zh) * 2015-12-21 2016-05-04 李奇 非同步基站系统下的无线测距定位方法及系统
CN105607039A (zh) * 2016-01-27 2016-05-25 宁波大学 非视距环境下基于到达时间差的稳健最小二乘定位方法
CN105738655A (zh) * 2015-12-23 2016-07-06 深圳市金立通信设备有限公司 一种方向传感器校准方法及终端
WO2016119462A1 (zh) * 2015-01-30 2016-08-04 中兴通讯股份有限公司 一种定位方法、移动终端、定位装置及通信系统
CN106405496A (zh) * 2016-08-31 2017-02-15 北斗时空信息技术(北京)有限公司 一种基于tdoa的室内定位方法
CN106932757A (zh) * 2017-03-31 2017-07-07 南京理工大学 一种基于tdoa和多普勒的月球车联合定位方法
CN107079399A (zh) * 2014-09-15 2017-08-18 高通股份有限公司 用于高效定位的ieee 802.11增强技术
US9743377B2 (en) 2009-09-29 2017-08-22 China Academy Of Telecommunications Technology Method, system and apparatus for acquiring cell location information
CN107450049A (zh) * 2017-07-10 2017-12-08 武汉米风通信技术有限公司 到达定位系统时间差的精准计算方法
CN107923964A (zh) * 2015-09-14 2018-04-17 红点定位公司 用于估计和补偿到达时间差估计中的nlos偏差的方法
CN108093474A (zh) * 2017-12-27 2018-05-29 重庆脉实智能制造有限公司 利用虚拟时间同步实现室内定位的方法及系统
CN110087186A (zh) * 2013-05-30 2019-08-02 高通股份有限公司 用于增强型往返时间(rtt)交换的方法和系统
CN110988797A (zh) * 2019-11-19 2020-04-10 北京智芯微电子科技有限公司 定位系统与定位方法
WO2020103908A1 (zh) * 2018-11-23 2020-05-28 中兴通讯股份有限公司 确定终端位置的方法、装置及存储介质
TWI745909B (zh) * 2019-04-02 2021-11-11 大陸商大唐移動通信設備有限公司 一種定位測量值的確定方法、裝置、電子設備及電腦存儲介質
CN114222362A (zh) * 2021-12-15 2022-03-22 武汉乾阳智能科技有限公司 定位方法及定位装置

Cited By (57)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8483717B2 (en) 2003-06-27 2013-07-09 Qualcomm Incorporated Local area network assisted positioning
US9335419B2 (en) 2003-06-27 2016-05-10 Qualcomm Incorporated Wireless network hybrid positioning
US9778372B2 (en) 2003-06-27 2017-10-03 Qualcomm Incorporated Wireless network hybrid positioning
US8971913B2 (en) 2003-06-27 2015-03-03 Qualcomm Incorporated Method and apparatus for wireless network hybrid positioning
US9810761B2 (en) 2003-06-27 2017-11-07 Qualcomm Incorporated Local area network assisted positioning
US10849092B2 (en) 2003-06-27 2020-11-24 Qualcomm Incorporated Local area network assisted positioning
US9749876B2 (en) 2003-06-27 2017-08-29 Qualcomm Incorporated Local area network assisted positioning
US9814016B2 (en) 2003-06-27 2017-11-07 Qualcomm Incorporated Local area network assisted positioning
US10841892B2 (en) 2003-06-27 2020-11-17 Qualcomm Incorporated Local area network assisted positioning
US10895648B2 (en) 2003-06-27 2021-01-19 Qualcomm Incorporated Method and apparatus for wireless network hybrid positioning
CN100428821C (zh) * 2003-09-04 2008-10-22 华为技术有限公司 一种抑制到达时间附加时延误差的方法
USRE45808E1 (en) 2004-06-18 2015-11-17 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 重庆邮电学院 一种单基站定位移动终端的方法
US8320934B2 (en) 2005-08-24 2012-11-27 Qualcomm Incorporated Dynamic location almanac for wireless base stations
US9042917B2 (en) 2005-11-07 2015-05-26 Qualcomm Incorporated Positioning for WLANS and other wireless networks
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
US10568062B2 (en) 2006-11-04 2020-02-18 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地址信息获取计算机精确地理位置信息的方法
CN101821645B (zh) * 2007-12-10 2013-06-19 真实定位公司 无线定位系统中的cdma信号的到达时间检测
CN101848414B (zh) * 2009-03-24 2013-04-17 电信科学技术研究院 确定信号传播时延的方法、对移动台定位的方法
CN101848414A (zh) * 2009-03-24 2010-09-29 大唐移动通信设备有限公司 确定信号传播时延的方法、对移动台定位的方法
CN101583078B (zh) * 2009-06-15 2011-09-21 北京邮电大学 一种蜂窝定位方法
CN101931858B (zh) * 2009-06-18 2013-01-30 华为技术有限公司 用户设备定位辅助数据配置的方法及装置
WO2010145559A1 (zh) * 2009-06-18 2010-12-23 华为技术有限公司 用户设备定位辅助数据配置的方法及装置
CN104869590A (zh) * 2009-08-13 2015-08-26 交互数字专利控股公司 支持定位测量的方法、本地服务器、wtru及其执行的方法
CN104869590B (zh) * 2009-08-13 2019-04-19 交互数字专利控股公司 支持定位测量的方法、位置服务器、wtru及其执行的方法
US10064156B2 (en) 2009-08-13 2018-08-28 Interdigital Patent Holdings, Inc. Method and apparatus for supporting positioning measurements
US9743377B2 (en) 2009-09-29 2017-08-22 China Academy Of Telecommunications Technology Method, system and apparatus for acquiring cell location information
CN110087186B (zh) * 2013-05-30 2021-02-12 高通股份有限公司 用于增强型往返时间(rtt)交换的方法和系统
CN110087186A (zh) * 2013-05-30 2019-08-02 高通股份有限公司 用于增强型往返时间(rtt)交换的方法和系统
CN103983275B (zh) * 2014-04-28 2016-11-09 上海融军科技有限公司 双基准信号源飞行器方向标定方法
CN103983275A (zh) * 2014-04-28 2014-08-13 上海融军科技有限公司 双基站飞行器方向标定方法
CN107079399A (zh) * 2014-09-15 2017-08-18 高通股份有限公司 用于高效定位的ieee 802.11增强技术
WO2016119462A1 (zh) * 2015-01-30 2016-08-04 中兴通讯股份有限公司 一种定位方法、移动终端、定位装置及通信系统
CN107923964B (zh) * 2015-09-14 2022-02-25 红点定位公司 用于估计和补偿到达时间差估计中的nlos偏差的方法
CN107923964A (zh) * 2015-09-14 2018-04-17 红点定位公司 用于估计和补偿到达时间差估计中的nlos偏差的方法
CN105548962B (zh) * 2015-12-21 2018-05-22 李奇 非同步基站系统下的无线测距定位方法及系统
CN105548962A (zh) * 2015-12-21 2016-05-04 李奇 非同步基站系统下的无线测距定位方法及系统
CN105738655A (zh) * 2015-12-23 2016-07-06 深圳市金立通信设备有限公司 一种方向传感器校准方法及终端
CN105607039A (zh) * 2016-01-27 2016-05-25 宁波大学 非视距环境下基于到达时间差的稳健最小二乘定位方法
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 武汉米风通信技术有限公司 到达定位系统时间差的精准计算方法
CN108093474A (zh) * 2017-12-27 2018-05-29 重庆脉实智能制造有限公司 利用虚拟时间同步实现室内定位的方法及系统
WO2020103908A1 (zh) * 2018-11-23 2020-05-28 中兴通讯股份有限公司 确定终端位置的方法、装置及存储介质
TWI745909B (zh) * 2019-04-02 2021-11-11 大陸商大唐移動通信設備有限公司 一種定位測量值的確定方法、裝置、電子設備及電腦存儲介質
CN110988797A (zh) * 2019-11-19 2020-04-10 北京智芯微电子科技有限公司 定位系统与定位方法
CN114222362A (zh) * 2021-12-15 2022-03-22 武汉乾阳智能科技有限公司 定位方法及定位装置
CN114222362B (zh) * 2021-12-15 2024-03-15 武汉乾阳智能科技有限公司 定位方法及定位装置

Also Published As

Publication number Publication date
CN1202688C (zh) 2005-05-18

Similar Documents

Publication Publication Date Title
CN1413058A (zh) 利用信号达到时间和时间差估计移动台位置的方法及装置
CN1185897C (zh) 一种移动台位置的估计方法及装置
JP5450081B2 (ja) ユーザ装置の位置を判定する方法、ソフトウエア及びシステム
CN1267746C (zh) 获取卫星定位系统信号的方法和装置
CN1607858A (zh) 用于检测网络中终端位置的系统和方法
CN101049041A (zh) 根据属于多个网络的基站的有关信息对无线装置进行定位
CN1468023A (zh) 位置测定方法,在该方法中使用的终端以及服务器
CN1852591A (zh) 一种定位移动终端的方法
CN1483148A (zh) 适用于混合定位系统确定误差估计的方法和装置
CN1524188A (zh) 根据初始粗位置估计提供经改善的位置估计的方法和装置
CN1416530A (zh) 定位方法和装置
CN1722905A (zh) 位置信息提供系统、基站和随其使用的位置信息提供方法
CN1480873A (zh) 用于检测终端在网络中位置的方法、系统和装置
CN1841084A (zh) 混合测距方法
CN1898975A (zh) 到达时间差(tdoa)/全球定位系统(gps)混合无线定位系统
CN1500357A (zh) 验证移动站位置方位的方法和系统
CN1465015A (zh) 利用传播天文历的广域参考网定位移动接收机的方法和装置
CN1776448A (zh) 基于数字广播电视信号的无线电组合定位方法
JPWO2009054069A1 (ja) 測位用受信装置及び移動端末の測位方法
TWI444073B (zh) Location and Tracking of Low Complexity Decentralized Message Transmission with Multiple Heterogeneous Data Fusion
CN1859721A (zh) 一种到达时间差定位方法
CN1235429C (zh) 一种位置估计方法
CN109889971B (zh) 应用于大型室内环境的基站三维协作定位方法
WO2019001175A1 (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

Granted publication date: 20050518

Termination date: 20161018

CF01 Termination of patent right due to non-payment of annual fee