沈家慶，姜維

（中國民航飛行學(xué)院航空工程學(xué)院，廣漢 618307）

0 引言

隨著中國經(jīng)濟的增長，航空運輸業(yè)也得到了快速的發(fā)展，為保證航空運輸?shù)陌踩?，工業(yè)界和學(xué)術(shù)界對新型監(jiān)控技術(shù)的研究力度也在不斷加深。多點定位技術(shù)[1]是通過在地面上布置多個地面基站來同時接收同一目標發(fā)出的信號，利用不同地面站接收到信號的時間差（TDOA）和接收角度（AOA）來對目標進行精確定位。根據(jù)目標到達兩個地面站的時間差可繪制一條以兩個地面站為焦點的雙曲線，目標位置就在這條雙曲線上。當有三個地面站時，就能繪制另外兩條雙曲線來確定目標在平面上的唯一位置。從理論上講，只要有四個地面站，就可以在空間上確定一點，在實際中就可以投入應(yīng)用。但在一般的實際應(yīng)用中都會布置四個以上的地方站以此來減少定位誤差，提高對監(jiān)控目標定位的精確度。多點定位系統(tǒng)具有成本低、精度高、易維護等一系列優(yōu)點，在國內(nèi)外已經(jīng)得到了廣泛的應(yīng)用。因此，對多點定位技術(shù)特別是對其定位算法的研究具有極其深刻的意義。

現(xiàn)有關(guān)于多點定位所運用的算法有很多且大多數(shù)都是基于TDOA的[2]，Taylor級數(shù)展開法[3]是最經(jīng)典的一種通過迭代來解決關(guān)于TDOA代數(shù)非線性問題的方法,迭代過程中每一步都是沿著當前點函數(shù)值下降方向。這種方法需要一個初始的預(yù)測值，然后通過求解局部線性最小二乘解來一步步的提高估計值的準確性。由于初始值的存在，這種方法也存在著一定的局限性，如果初始值與真實值距離過大，在操作中就不能保證迭代的收斂性，甚至可能會發(fā)散導(dǎo)致得不到所要求的結(jié)果。因為展開式每一步都要通過最小二乘法來求解，在實際計算過程中的計算量也是非常巨大的。而另外一種經(jīng)典的Chan算法[4]的計算量就要小很多，是較優(yōu)的求解雙曲線遞歸方程的方法。Chan算法只需要兩步迭代就可以得到最終結(jié)果，但如果其定位基站位置和TDOA定位精度存在較大誤差的情況下，其定位能力會受到較大影響。針對這個問題，本文利用誤差上限，利用文獻[5，6]提到的方法，解偽線性方程組，得到定位點的幾何位置以及它的斜距，再利用位置斜距的相關(guān)性，在定位點幾何位置的二范數(shù)、斜距相等的條件下，從而獲得一個相對誤差最小的定位解。再將此定位解作為初值代入Taylor級數(shù)展開，可得到一個精確度較高且穩(wěn)定的結(jié)果。

1 誤差上限法和Taylor級數(shù)展開法介紹

1.1 誤差上限法

假設(shè)估計目標物體的所在位置m=[x,y,z]T，地面接收站位置，i=1,2,3,4…,N（N≥6），地面接收站i與目標距離為ri，ri=‖ni-m ‖，地面接收站i與地面站1距離為為TDOA 的測量值，c為光速，由空間定位可得：

即：

又因為：

由式（2）式（3）可以得到：

將（4）改成矩陣形式：

其中：

A，b中都含有噪聲，在TDOA測量誤差或布站誤差（布站誤差是指第i個基站相對于原點的距離偏差）較大的情況下，以往定位方法的定位精度都有待提高?？梢愿鶕?jù)布站情況估計其布站誤差上限εA和測量誤差上限εb，令：

A0=A+δA，b0=b+δb

其中 A，b為測量值，δA，δb為誤差，A0，b0為真實值，我們可以得到

（‖x‖2為x的二范數(shù)）。

為了求出以上條件的解，首先需要求解：

根據(jù)文獻[5-6]，將（6）轉(zhuǎn)化為：

令A(yù)的SVD分解為：

求出定位解：

其中為方程的正根。

由于噪聲的存在，上述等式?jīng)]有嚴格成立?？梢赃x擇在噪聲干擾最小的條件下求得最精確的定位解?？傻孟率剑?/p>

由于α為加權(quán)系數(shù)，可以利用拉格朗日算子來進行計算（10）。

結(jié)果為：

其中，α不等于-1。

最終定位。

1.2 Taylor級數(shù)展開法

利用Taylor級數(shù)展開法計算，首先猜測目標位置為(xv,yv,zv) ，有 x=xv+ δx，y=yv+δy，z=zv+δz。 δx,δy，δz為估計誤差。用 Taylor級數(shù)令 fi等于（12）

在猜測位置( )xv,yv,zv線性化，得到：

再定義矩陣 M，δ，Z，e：

其中M為猜測誤差的權(quán)，δ為猜測誤差，Z為猜測位置觀測函數(shù)值與實際值的誤差，e為測量誤差。

則公式可近似表達為Mδ?Z-e,該公式的轉(zhuǎn)變主要是將對初始位置猜測的誤差值并入了測量誤差之中，目標位置未知，由此猜測xv,yv,zv，M,Z可以估算，δ，e無法估算 。因為有 x=xv+δx，y=yv+δy，z=zv+δz，設(shè)定一個閾值ε，當丨δx丨+丨δy丨+丨δz丨<ε時，可近似認為猜測位置為目標的真實位置。

2 協(xié)同算法

雖然Taylor算法的迭代求解計算量相對較小，但由于其計算初始值的存在，如果估測值與真實值誤差過大會導(dǎo)致整個算法發(fā)散從而得不到結(jié)果。本文提出一種新的協(xié)同算法能有效地規(guī)避這種情況，提高最終定位結(jié)果的精確度。首先根據(jù)式：

計算求得誤差上限法得到的初始定位點，然后將作為Taylor算法所需要的初始值代入式：

進行展開，通過Taylor算法迭代計算出最終的計算結(jié)果。圖1為本文提出的協(xié)同算法的計算流程圖。

圖1 協(xié)同算法計算流程圖

3 仿真分析

為了比較協(xié)同算法與其他算法的計算精度，假設(shè)六個站的地址為（0km,0km,0km），（9.08km,7.02km,3km）,（-9.08km,7.02km,3km）,（12.18km,16.06km,3km）,（12,18km,-16.06km,3km）,（10km，5km，10km）。測量目標的位置為（3km,6km,9km）。在使用MATLAB對協(xié)同算法進行仿真后，3種算法在標準相對誤差情況下定位精度誤差情況比較如圖2所示，從圖中可以看出，泰勒級數(shù)展開法的精度是相比較最低的，協(xié)同定位算法的精度比泰勒級數(shù)展開法和誤差上限法的精度都有了明顯的提高。

圖2 算法定位精度仿真圖

4 結(jié)語

以往的定位算法為了能夠簡單快速地計算出結(jié)果而忽略了許多誤差的存在，所以導(dǎo)致了最后計算出來的結(jié)果與真實值之間存在較大的誤差，本文提出的協(xié)同定位算法，在計算中引入了測量誤差和選址誤差并通過利用誤差上限對誤差大小的控制，減小了計算結(jié)果與真實值之間的差距，在通過與Taylor算法的優(yōu)劣互補，相互協(xié)同計算工作，從結(jié)果上看進一步提高其定位精度，具有一定的應(yīng)用價值。但文中方法的缺點也比較明顯，由于引進了兩項誤差并融合了Taylor算法協(xié)同計算，用到了矩陣的分解和代數(shù)計算，導(dǎo)致計算量相對巨大，在以后的應(yīng)用中可根據(jù)實際要求選擇所適合的算法進行計算研究。

參考文獻：

[1]Multilateration（MLAT）Concept of use[R].Bangkok:Internatonal Civil Aviation Organization Asia And Pacific Office,2007.

[2]郭華.TDOA定位技術(shù)的基本原理和算法[J].西安郵電學(xué)院學(xué)報，2007（1）:19-24.

[3]Foy W H.Position Location Solutions by Taylor Seriesestimation[J].IEEE Trans.on Aerospace and Electronic Systems,1976，12（2）:187-194.

[4]Chan,Y.T.Ho,K.C.A Simple and Efficient Estimator for Hyperbolic Location.Signal Processing,IEEE Transactions on,1994，42（8）:1905-1915.

[5]Chandrasekaran S,Golub G H,Gu M,Sayed A H.Parameterestimation in the Presence of Bounded Modeling Errors.IEEE Signal Processing,1997，4（7）:195-197.

[6]S.Chandrasekaran.Parameter Estimation in the Presence of Bounded Data Uncertainties.Siam Journal on Matrix Analysis&Applications,1998，19（1）:235-252.

[7]馬廣亮.一種用于多點定位系統(tǒng)的實時高精度位置解算方法.科技資訊,2016，14（20）：121-124.

[8]王洪,劉昌忠,汪學(xué)剛,吳宏剛.一種多點定位的目標位置精確解算方法.航空學(xué)報，2011，32（7）:1269-1274.

[9]宮峰勛.終端區(qū)及場面多點定位中TSOA與TDOA算法性能分析.南京航空航天大學(xué)學(xué)報,2015;47（6）:818-826.

[10]劉剛,趙國慶.時差定位與兩種測時差方法[J].電子對抗,2006,20（1）:21-25.

[11]汪波,薛磊.基于遺傳算法的TDOA定位系統(tǒng)的最優(yōu)布站算法[J].系統(tǒng)工程與電子技術(shù)，2009,31（9）:2125-2128.

[12]朱厚慶.到達時間差（TDOA）測向定位研究[J].電訊技術(shù),2007（1）:53-56.

[13]毛永毅,李明遠,張寶軍.基于RBF神經(jīng)網(wǎng)絡(luò)的TDOA/AOA定位算法[J].計算機工程，2008,34（2）:52-55.

[14]鄧平,李莉,范平志.一種TDOA/AOA混合定位算法及其性能分析[J].電波科學(xué)學(xué)報，2002,17（6）:633-636.

[15]俞麗娜,曾連蓀,金志華.基于DTOA的無線定位算法研究及精度分析[J].計算機測量與控制,2006（9）:1247-1249.