Chan定位算法与TDOA估计精度的关系

（12）
其中：
∑{ } σmax
= ΔD ⋅
λ
/
2
⋅
⎧ ⎨(M
c⎩
M
−1) −
m=2
cos ⎡⎣ϕ1,0
−ϕm,0 ⎤⎦
⎫−1/ 2 ⎬
。(13)
⎭
2 物理意义和算法流程
式（12）具有的物理含义非常重要。通过分析，本文得到 以下几个结论：
结论 1 在给定定位估计误差ΔD 下，TDOA 允许的最大估计 误差和定位系统的拓扑结构存在关联。在式（13）中，本文采 用了余弦的形式表达了系统的拓扑结构，使其更加直观和清晰。 显然，拓扑结构是由定位系统各基站的布站结构和目标移动台 的位置来共同决定的。它们共同影响着 TDOA 允许的最大误 差值。
【摘 要】在基于 TDOA（到达时间差）的定位算法中，Chan 方法由于其独特的优势得到广泛应用，但如何估计 Chan 定
位算法所需的时延估计精度暂时没有明确的函数表达式。针对 Chan 定位算法，推导了该算法定位估计误差与定位参数(TDOA)
估计精度之间的函数关系。给出了推导的具体步骤和函数关系的表达式，分析了该表达式所具有的物理意义。在移动台无线
结论 2 当所有的基站与目标手机呈直线排列时，σmax 趋向 于无穷大。这个结论实际上说明，在此种情况下，无论 TDOA 估计的精度有多高，最后的定位误差都会为无穷大。
结论 3 在给定定位估计误差ΔD 下，TDOA 估计精度的要 求与基站数有关。由式（13）可以看出，当基站数目不同时， 即使在相同的ΔD 下，σmax 会可能不同。具体的关系可以从数值 仿真中看出。
阵的逆。式（3）中，[.](M-1)×2 中的下标表示矩阵的二维维数。 式（2）的结果同时也表示以第一个基站作为系统的参考
基站。
假设 M-1 个 TDOA 的估计值都具有相同的方差值大小，
且各 TDOA 估计值之间是相互独立的，则有： Q = σ 2 I(M −1)×(M −1) ，
（4）
其中，I(M-1)× (M-1)为(M-1)× (M-1)维的单位矩阵；σ2 为 TDOA 估计精度的方差。
【Abstract】Chan method, for its unique advantages, is widely applied in TDOA (Time Difference of Arrival) positioning algorithm. However, there is now no explicit function expression on how to estimate the time-delay estimation accuracy required by Chan location algorithm. In this paper, the relation function between estimation errors of Chan location algorithm and the estimation accuracy of location parameter (TDOA) is deduced. The detailed deduction and the expression of the relation function are given and the physical meanings of the relation function also analyzed. In the design of MS (Mobile Station) radio location system, the least time delay estimation accuracy could be directly computed by this relation function according to the location accuracy required by the location system.
【Key words】radio location; location accuracy; TDOA; time delay estimation
0 引言
1996 年，美国 FCC 明确要求各移动通信运营商为移动 用户提供 E-911 服务[1]，因此多个基站利用移动台无线信号 到达时刻 TOA（Time-of-Arrival）[2]和到达时间差 TDOA (Time-Difference-of-Arrival) 进 行 定 位 的 方 式 在 近 年 来 受到了广泛关注[3-4]。
rm,0 = (xm − x0 )2 + ( ym − y0 )2 ，
（1）
由文献[6]，估计误差的协方差矩阵为： φ = c2 (Gt,0TQ−1Gt,0 )−1 ，
（2）
其中，c = 3*108(m/s)，且
⎡(x1 − x0 ) / r1,0 − (x2 − x0 ) / r2,0 ⎤
⎢ ⎢( y1
−
y0
)
/
r1,0
−
( y2
−
y0
)
/
r2,0
⎥ ⎥
Gt , 0
=
⎢ ⎢
#
#
⎥ ⎥
，
⎢( ⎢
x1
−
x0 ) /
r1,0
− (xM
−
x0 )
/
rM ,0
⎥ ⎥
⎢⎣( y1 − y0 ) / r1,0 − ( yM − y0 ) / rM ,0 ⎥⎦(M −1)×2
（3）
式（2）中，矩阵 Q 表示 TDOA 估计误差的协方差矩阵，其 维数为(M-1)× (M-1)；[.]T 表示取矩阵的转置；[.]-1 表示求矩
在众多基于 TDOA 的定位算法[5]中，Chan 算法[6]得到了 广泛的应用。这主要是因为该算法具备三大优势：①算法不 需要初值；②仅进行两次迭代就可求得最终结果；③算法的
收稿日期：2009-01-04。 作者简介：丁宏毅（1974-），男，工程师，主要从事卫星导航定位的研
究工作；柳其许（1966-），男，高级工程师，主要从事 GPS 伪码测量、卫星频率、卫星导航定位的研究工作；王 巍 （1976-），男，工程师，主要从事移动通信的研究工作。
将式（3）和式（4）代入式（2）中，可得：
∑ ⎡
⎢ ⎢ ⎢
M m=
2
⎡ ⎢ ⎢⎣
(
x1
−x r1,0
0
)
−
( xm
−
x0
)
⎤2 ⎥
rm,0 ⎥⎦
⎤ −1 ⎥ ⎥ ⎥
φ
=
c2σ
2
⎢ E (( xˆ ⎢⎢E(( yˆ
− −
x0 )( yˆ − y0 )(xˆ −
y0 x0
)) ))
⎥ ⎥
，
⎥
∑ ⎢
⎢M ⎢ ⎣m=2
=
c2σ λ
2
⋅
∑M
⎪⎧⎡ ⎨⎢
(x1
−
x0
)
−
( xm
−
x0
)
⎤2 ⎥
+
⎡ ⎢
(
y1
−
y0
)
−
(
ym
−
y0
)
⎤2 ⎥
⎪⎫ ⎬
,（8）
m=2 ⎪⎩⎢⎣ r1,0
rm,0 ⎥⎦ ⎢⎣ r1,0
rm,0 ⎥⎦ ⎪⎭
其中，ΔD 表示 Chan 算法定位估计的误差且大于等于 0。令：
tan(ϕk )
=
yk xk
定位系统的设计过程中，根据所需的定位精度，通过该表达式可以直接计算出系统所必须的最小时延估计精度。
【关键词】无线定位；定位精度；信号到达时间差；时延估计
【中图分类号】TN929.53
【文献标识码】A
【文章编号】1002-0802(2010)03-0134-03
Relation Between Chan Location Algorithm and TDOA Estimation Accuracy
DING Hong-yi①， LIU Qi-xu①， WANG Wei②
(①China National Technical Service Corporation for Surveying and Mapping, Beijing 100088, China; ②Southwest Inst Of Electron & Telecom Techn, Chengdu Sichuan 610041, China))⎤ ⎥，（6）m=2 ⎢⎣ r1,0
rm,0 ⎥⎦ ⎢⎣ r1,0
rm,0 ⎥⎦
现需对式（5）所表示的 2×2 阶矩阵求逆。
假设λ为式（5）中[.]的行列式的值，则有：
∑ ⎡
⎢
M
⎡ ⎢
(
y1
−
y0
)
−
(
ym
−
y0
)
φ
=
c2σ
2
1
⎢m=2 ⎢
⎢⎣
r1,0
rm,0
−E((xˆ
−
x0
)(
yˆ
−
y0
⎤2 ))⎥
基于式（12），当给定定位估计误差ΔD 时，可以求出 TDOA 允许的最大估计误差，算法的具体流程如下：
① 确定定位系统的基站坐标和目标手机的坐标； ② 计算式（5）中[.]的行列式值，得到λ； ③ 由式（13），求出 TDOA 允许的最大估计误差。
3 仿真和分析
考虑几种具有对称结构的定位网络，坐标的单位都为 米。A 类（三角形）基站的坐标分别为：(200,200)、(-200,200) 和(0,-282.8)，此时 M=3；B 类（四边形）基站的坐标分别为： (200,200)、(-200,200)、(-200,-200)和(200,-200)，此时 M = 4； C 类（六边形）基站的坐标分别为：(200,200)、(-200,200)、 (-200,-200)、(200,-200)、(-282.8,0)和(282.8,0)，此时 M=6。