曲付勇 孟祥偉

?

基于約束總體最小二乘方法的到達時差到達頻差無源定位算法

曲付勇 孟祥偉*

(海軍航空工程學院電子信息工程系 煙臺 264001)

兩步加權最小二乘方法(two-stage WLS)是求解TDOA/FDOA無源定位問題的經(jīng)典線性方法，但也存在著定位偏差和均方誤差對測量噪聲的適應能力較差的缺點。該文根據(jù)TDOA/FDOA的偽線性定位方程組特點，將其建立為一種帶約束條件的約束總體最小二乘(CTLS)模型，并采用拉格朗日乘子法求解帶約束條件的CTLS問題，建立了幾種最小二乘類定位方法的統(tǒng)一解，從而將約束加權最小二乘(CWLS)定位解和約束最小二乘(CLS)定位解變?yōu)樵撐腃TLS定位解的特例。仿真表明，該文方法比兩步加權最小二乘方法具有更低的均方誤差，并能夠有效減小定位偏差，因而具有更好的測量噪聲適應能力。

無源定位；到達時差；到達頻差；最小二乘；偏差；均方誤差

1 引言

2 時差頻差定位方程的描述

將式(1)兩邊求導，得到距離差變化率的公式為

這里采用和文獻[7]中一致的表達式。

3 帶約束條件的CTLS定位求解

利用Lagrange乘子法，將式(9)轉(zhuǎn)化為對式(10)的最小化：

整理式(14)可得

本文基于CTLS的定位算法計算過程如下：

4 討論

(1)本文定位解與two-stage WLS的關系。two- stage WLS中給出的偽線性方程組[2]是

(3)文獻[7]中CTLS定位解是通過極小化

5 仿真分析

圖1給出了CTLS, CWLS與two-stage WLS 3種定位算法對近場源的定位位置偏差和速度偏差隨接收站位置誤差變化的性能曲線。由圖1可以看出，CTLS方法的位置偏差和速度偏差相比two-stage WLS和CWLS都得到了顯著的降低，表明CTLS方法能夠有效減弱偏差的影響。

圖2給出了3種定位算法對近場源的定位位置均方誤差和速度均方誤差隨接收站位置誤差變化的性能曲線。圖2中，CTLS方法的位置均方誤差的“門限效應”相比two-stage WLS有了明顯的改善。當接收站位置誤差分別大于7.5 dB, 12.5 dB時，two- stage WLS和CWLS的定位位置均方誤差開始偏離CRLB，性能曲線開始發(fā)散，而CTLS方法直到大于15 dB時其位置均方誤差才大于CRLB，并且隨著接收站位置誤差的增大，CTLS方法的位置均方誤差要低于CWLS與two-stage WLS定位位置均方誤差。從圖2可以看出，相比位置均方誤差，速度均方誤差對時頻差測量誤差的變化更為敏感。當接收站位置誤差分別大于5 dB, 7.5 dB時，two-stage WLS和CWLS的定位位置均方誤差開始偏離CRLB，而CTLS方法大于10 dB時其位置均方誤差才大于CRLB，并且隨著接收站位置誤差的增大，CTLS方法的速度均方誤差要低于CWLS與two-stage WLS定位位置均方誤差。

圖3給出了3種定位算法對遠場輻射源的位置偏差和速度偏差隨接收站位置誤差變化的性能曲線。顯然，隨著位置誤差的增大，3種方法的偏差不斷增大，這表明3種方法只能隨接收站位置誤差的減小獲得漸進無偏的性能。由圖3可以看出，CTLS定位方法的位置偏差和速度偏差要低于CWLS和two-stage WLS的位置偏差和速度偏差。

圖1 3種算法得到的近場輻射源位置和速度偏差隨接收機位置誤差的變化曲線

圖2 3種算法得到的近場輻射源位置和速度均方誤差隨接收機位置誤差的變化曲線

圖3 3種算法得到的遠場輻射源位置和速度偏差隨接收機位置誤差的變化曲線

圖4給出了3種定位算法對遠場輻射源的位置均方誤差和速度均方誤差隨接收站位置誤差變化的性能曲線，以及位置和速度的CRLB。從圖4可以看出，當接收站位置誤差分別大于-7.5 dB, -2.5 dB時，two-stage WLS和CWLS的定位位置均方誤差開始偏離CRLB，性能曲線開始發(fā)散，而CTLS方法直到大于0 dB時其位置均方誤差才大于CRLB，并且隨著接收站位置誤差的增大，CTLS方法的位置均方誤差要低于CWLS與two-stage WLS定位位置均方誤差。相比位置均方誤差，速度均方誤差對時頻差測量誤差的變化更為敏感。當接收站位置誤差分別大于-20 dB, -12.5 dB時，two-stage WLS和CWLS的定位速度均方誤差開始偏離CRLB，而CTLS方法直到大于-10 dB時其速度均方誤差才大于CRLB，并且隨著接收站位置誤差的增大，CTLS方法的速度均方誤差要低于CWLS與two-stage WLS定位速度均方誤差。

仿真2 考察遠場源中3種算法定位偏差和均方誤差隨時頻差測量誤差變化的性能。

圖6給出了3種定位算法的定位位置和速度均方誤差隨時頻差測量誤差的變化曲線，以及位置和速度的CRLB。從圖6可以看出，當測量誤差分別大于7.5 dB, 12.5 dB時，two-stage WLS和CWLS的定位位置均方誤差開始偏離CRLB，而CTLS方法直到大于15 dB時其位置均方誤差才大于CRLB。相比位置均方誤差，速度均方誤差對時頻差測量誤差的變化更為敏感。當測量誤差分別大于-2.5 dB, 5 dB時，two-stageWLS和CWLS的定位位置均方誤差開始偏離CRLB，而CTLS方法直到大于7.5 dB時，其位置均方誤差才大于CRLB。

綜合分析3種方法的定位偏差和均方誤差，當均方誤差開始偏離CRLB時，此時估計已為有偏估計，由于CWLS方法考慮定位方程的約束關系，能夠減弱其對定位偏差的影響，而本文中CTLS方法不僅考慮了約束關系，而且降低了偽線性方程中由于觀測矩陣和測量向量相關導致的偏差影響，因而具有更低的偏差，且均方誤差相比CWLS與two-stage WLS也得到了有效改善。

仿真3 在仿真1條件下以近場源為例，仿真本文算法的收斂性。

圖4 3種算法得到的遠場輻射源位置和速度均方誤差隨接收機位置誤差的變化曲線

圖5 3種算法得到的遠場輻射源位置和速度偏差隨時頻差誤差的變化曲線

圖6 3種算法得到的遠場輻射源位置和速度均方誤差隨時頻差誤差的變化曲線

圖7 本文算法得到的近場輻射源位置迭代曲線

6 結(jié)論

本文根據(jù)TDOA/FDOA偽線性方程組的特點，研究了基于約束總體最小二乘的無源定位算法。基于拉格朗日乘子法建立了幾種最小二乘類定位方法的統(tǒng)一解，并給出了具體的求解方法。同時討論了two-stage WLS, CWLS, CLS等定位方法解與本文定位解的關系，說明CWLS和CLS為本文定位解的特例。通過仿真表明，在時頻差誤差和接收機位置誤差適度的條件下本文方法比two-stage WLS和CWLS方法具有更低的均方誤差，并能夠有效減小定位偏差，因而具有更好的測量噪聲適應能力。

[1] Ho K C and Xu W. An accurate algebraic solution for moving source location using TDOA and FDOA measurements[J]., 2004, 52(9): 2453-2463.

[2] Ho K C, Lu X, and Kovavisaruch L. Source localization using TDOA and FDOA measurements in the presence of receiver location errors: analysis and solution[J]., 2007, 55(2): 684-696.

[3] Sun M and Ho K C. An asymptotically efficient estimator for TDOA and FDOA positioning of multiple disjoint sources in the presence of sensor location uncertainties[J]., 2011, 59(7): 3434-3440.

[4] Sun X Y, Li J D, Huang P Y,..Total least-squares solution of active target localization using TDOA and FDOA measurements in WSN[C]. The 22th International Conference on Advanced Information Networking Applications Workshops, Okinawa, Japan, 2008: 995-999.

[5] Guo F C and Ho K C. A quadratic constraint solution method for TDOA and FDOA localization[C]. International Conference on Acoustics Speech and Signal Processing, Prague, Czech, 2011: 2588-2591.

[6] Qu F Y, Guo F C,.. Constrained location algorithms based on total least squares method using TDOA and FDOA measurements[C]. IET International Conference on Automatic Control and Artificial Intelligence, Xiamen, China, 2012: 2587-2590.

[7] Yu H, Huang G, and Gao J. Constrained total least-squares localisation algorithm using time difference of arrival and frequency difference of arrival measurements with sensor location uncertainties[J].,&, 2012, 6(9): 891-899.

[8] Ho K C. Bias reduction for an explicit solution of source localization using TDOA[J]., 2012, 60(5): 2101-2114.

[9] Xu B, Qi W D, Wei L,.. Turbo-TSWLS: enhanced two-step weighted least squares estimator for TDOA-based Localisation[J]., 2012, 48(25): 1597-1598.

[10] Amar A, Leus G, and Friedlander B. Emitter localization given time delay and frequency shift measurements[J]., 2012, 48(2): 1826-1837.

[11] Cheung K W, So H C,.. A constrained least squares approach to mobile positioning: algorithms and optimality[J]., 2006: 1-23.

[12] Yu H, Huang G, Gao J,..Practical constrained least-square algorithm for moving source location using TDOA and FDOA measurements[J]., 2012, 23(4): 488-494.

[13] Stoica P and Li J. Source localization from range-difference measurements[J]., 2006, 23(6): 63-65.

[14] Abatzoglou T J, Mendel J M, and Harada G A. The constrained total least squares technique and its applications to harmonic superresolution[J]., 1991, 39(5): 1070-1087.

[15] Schaffrin B. A note on constrained total least-squares estimation[J]., 2006, 417: 245-258.

[16] Foy W H. Position-location solutions by Taylor-series estimation[J]., 1976, 12(2): 187-194.

曲付勇： 男，1984年生，博士生，研究方向為無源定位、數(shù)據(jù)壓縮.

孟祥偉： 男，1966年生，教授，博士生導師，研究方向為雷達信號檢測、信號理論、定位等.

Source Localization Using TDOA and FDOA MeasurementsBased on Constrained Total Least Squares Algorithm

Qu Fu-yong Meng Xiang-wei

(,,264001,)

The two-stage Weighted Least Squares (WLS) method is a well-known linear approach in Time- Difference-Of-Arrival (TDOA) and Frequency-Difference-Of-Arrival (FDOA) passive localization. But this method can only attain the CRLB in a modest noise environment and the bias of the localization result is significant for strong noise. This paper discusses a Constrained Total Least Square (CTLS) solution to the pseudo linear equations with two constrains for TDOA/FDOA localization. A unified expression for several LS solutions is derived based on Lagrange multiplier. The Constrained Weighted Least Square (CWLS) method and Constrained Least Square (CLS) localization method reduce to the special cases of the localization solution. The simulation results show that the proposed method has lower Mean Square Error (MSE) and lower bias compared with the two-stage WLS method, and it is more robust to noise.

Passive localization; Time-Difference-Of-Arrival (TDOA); Frequency-Difference-Of-Arrival (FDOA); Least Square (LS); Bias; Mean Square Error (MSE)

TN971

A

1009-5896(2014)05-1075-07

10.3724/SP.J.1146.2013.01019

孟祥偉 mengxw163@163.com

2013-07-11收到，2013-11-26改回

國家自然科學基金(61179016)資助課題