崔劍鋒，楊寶山，杜敬林

（中船重工大連測控技術(shù)研究所，遼寧 大連 116013）

多站式時差定位方法廣泛應(yīng)用于雷達(dá)和聲納等定位系統(tǒng)中，探測網(wǎng)絡(luò)由若干個測量站構(gòu)成，時差定位系統(tǒng)是通過處理多個測量站采集到的信號到達(dá)時間數(shù)據(jù)對輻射源進(jìn)行定位.在三維空間中，輻射源信號到達(dá)兩站的時間差確定了一對以兩站為焦點(diǎn)的雙曲面，三維定位則需要四個測量站產(chǎn)生3對雙曲面，交匯解出輻射源的位置.在輻射源某一坐標(biāo)值（如深度信息）可以事先確定的情況下，采用3站式測量可以得到非線性解，采用4站測量可以求出線性解.

實(shí)際工作中發(fā)現(xiàn)，在人為誤差、設(shè)備誤差均較小的情況下，定位系統(tǒng)仍可能在某些區(qū)域出現(xiàn)較大的定位誤差（本文把定位誤差較大的區(qū)域表述為低精度區(qū)）.例如在某導(dǎo)彈落點(diǎn)（落點(diǎn)位于海上）測量過程中，采用聲、光等不同定位系統(tǒng)得到的測量結(jié)果在某些情況下會出現(xiàn)較大差異，本文旨在基于數(shù)值模擬對時差定位誤差的空間分布規(guī)律進(jìn)行探究，并分析其規(guī)律的形成原因.

時差定位系統(tǒng)由若干個位置已知的接收站（或接收陣元）組成，設(shè)輻射源的位置為未知量M（x，y，z），輻射源信號到達(dá)站點(diǎn)i的距離為di，則有距離方程

任意兩站點(diǎn)的距離方程做差，以1作為參考站點(diǎn)，整理可得距離差方程

式中：di1為第i號站點(diǎn)與1號站點(diǎn)到輻射源的距離差；ri為i號站點(diǎn)到坐標(biāo)原點(diǎn)的距離.

由式（2）構(gòu)成的方程組就可以解算輻射源目標(biāo)位置M（x，y，z），本文著重分析已知目標(biāo)深度信息時的誤差空間分布特性，即z坐標(biāo)作為已知量.

1 三站式定位分析

若有三個站點(diǎn)接收到信號，根據(jù)式（2）可以得到二維方程組，改寫為矩陣形式

解方程式可以求出坐標(biāo)用d1表達(dá)，帶入式（1）中，得到［3］

根據(jù)距離為正可以排除負(fù)解，當(dāng)兩解都為正時出現(xiàn)模糊解.

設(shè)定聲速測量誤差小于0.1%，時差測量誤差小于0.04ms，站點(diǎn)位置測量誤差小于3m，目標(biāo)深度誤差小于3m，設(shè)定原始誤差服從均勻分布，利用計算機(jī)在確定空間區(qū)域內(nèi)進(jìn)行數(shù)值模擬，得到三站式定位誤差的空間分布如圖1，圖2所示，誤差值根據(jù)圖示色棒比對，深色表示高精度區(qū)，淺色反之，站點(diǎn)坐標(biāo)在圖中給出.分析圖示，可見三站式定位方案中的高精度區(qū)域具有指向性，在站點(diǎn)相對的方向高精度區(qū)覆蓋范圍較大，在站點(diǎn)兩側(cè)沿著三角形邊長延長方向出現(xiàn)低精度區(qū)；且站點(diǎn)對邊相對越寬，高精度區(qū)越擴(kuò)展；對邊相對越窄，高精度區(qū)域越壓縮.因此，三站式定位誤差空間分布受到站點(diǎn)位置的影響.

圖1 三陣元定位Fig.1 Three element location

圖2 三陣元定位Fig.2 Three element location

2 四站式定位分析

2.1 解算方法與數(shù)值模擬

采用四站式定位方法，可以得到三維線性方程組，求得唯一解，方程組形式如下

設(shè)定誤差來源與大小同前面三站式定位誤差分析一致，對不同陣型下四站式定位誤差進(jìn)行數(shù)值模擬，結(jié)果如圖3～圖8所示.

圖3 正方形陣Fig.3 Square arry

圖4 菱形陣Fig.4 Diamond array

圖5 菱形陣Fig.5 Diamond array

圖6 菱形陣Fig.6 Diamond array

圖7 平行四邊形陣Fig.7 Parallelogram array

圖8 任意四邊形陣Fig.8 Quadrilateral array

由圖可見，對正方形陣，低精度區(qū)分布在陣型的兩條中線附近，從陣中心到陣外圍，低精度區(qū)呈擴(kuò)大趨勢，從高精度區(qū)到低精度區(qū)的過渡時漸變的；對菱形陣，低精度區(qū)近似分布在某雙曲線附近，菱形越扁，低精度區(qū)越壓縮，并逐漸覆蓋陣元后的高精度區(qū)；進(jìn)而分析平行四邊形陣及任意四邊形陣，均有相同的規(guī)律.

2.2 分布規(guī)律原因分析

定位結(jié)果的誤差來源主要是人為誤差和設(shè)備誤差，但特定的分布規(guī)律應(yīng)該有其內(nèi)在原因.

對正方形陣，設(shè)陣元坐標(biāo)如圖3上所標(biāo)注，對系數(shù)矩陣做初等變換，當(dāng)矩陣不滿秩時，推知

對菱形陣，設(shè)陣元坐標(biāo)如圖6所示，推知系數(shù)矩陣的不滿秩條件與式（6）一樣；對平行四邊形陣，設(shè)陣元坐標(biāo)如圖7所示，推知系數(shù)矩陣的不滿秩條件也是式（6）；對任意四邊形，設(shè)陣元坐標(biāo)如圖8所示，推知矩陣不滿秩的條件為

依據(jù)奇異矩陣導(dǎo)出的方程可以繪出滿足條件的聲源位置，構(gòu)成一族曲線，如圖9～圖12中虛線所示.

圖9 正方形陣Fig.9 Square array

圖10 菱形陣Fig.10 Diamond array

圖11 平行四邊形陣Fig.11 Parallelogram array

圖12 任意四邊形陣Fig.12 Quadrilateral array

在任意陣型的四陣元定位中，對系數(shù)矩陣做初等變換，令矩陣不滿秩，可得到d1，d2，d3，d4的關(guān)系滿足

式（8）為非標(biāo)準(zhǔn)的雙曲方程，其中，A，B，C，D為常數(shù)參量，完全由陣元位置決定.特殊的，當(dāng)陣元分布為正方形時，雙曲線退化為正交的兩條直線.

解出的曲線與定位系統(tǒng)數(shù)值試驗(yàn)中的低精度區(qū)分布規(guī)律和位置一致，可互為驗(yàn)證.

由于存在測量誤差，當(dāng)聲源的空間位置逐漸逼近奇異矩陣導(dǎo)出的雙曲線所在位置時，空間定位誤差將逐漸增大.

對比三站式與四站式定位誤差的空間分布，可以發(fā)現(xiàn)三站式定位的高精度區(qū)域明顯小于四站式定位，但三站式定位方案得到的高精度區(qū)可以基本覆蓋四站式定位方案中的低精度區(qū)域，因此在實(shí)際定位過程中，采用兩種解算方案進(jìn)行互補(bǔ)，可能實(shí)現(xiàn)大部分區(qū)域的精確定位.

3 結(jié) 論

實(shí)際上，在深水情況下，聲速沿深度方向不是一個常數(shù)，聲線由聲源到接收換能器的傳播聲線不是直線，需要利用水平傳播距離與傳播時間的關(guān)系來進(jìn)行聲速補(bǔ)償.若忽略聲線彎曲，或在對聲速進(jìn)行補(bǔ)償?shù)那疤嵯?，基于文中涉及的多站式定位解算方法，可以得到如下結(jié)論：① 多站式時差定位方法都存在低精度區(qū)的問題；相同條件下，四站式時差定位方法的低精度區(qū)域遠(yuǎn)小于三站式時時差定位方法；② 三站式定位低精度區(qū)的產(chǎn)生原因是非線性解出現(xiàn)無解與模糊解；高精度定位區(qū)域的空間分布具有指向性；③ 四站式定位低精度區(qū)的原因是解算矩陣不滿秩；低精度區(qū)域與奇異矩陣導(dǎo)出的雙曲線重合.

多站式時差定位方法在解算方法上不可避免地存在低精度區(qū)的問題，如文中分析：三站式、四站式解算方法中的低精度區(qū)域和高精度區(qū)域能夠相互覆蓋，若結(jié)合多種定位方法對定位算法進(jìn)行互補(bǔ)與整合，可能減小甚至消除定位低精度區(qū).本文暫未對此進(jìn)行展開分析.

［1］ 俞志強(qiáng)，葉朝謀.四站三維時差定位模糊分析［J］.空軍雷達(dá)學(xué)院學(xué)報，2009，23（5）：370－372.

Yu Zhiqiang，Ye Chaomei.Study on ambiguity and non－solution of 3－D TDOA location systems［J］.Journal of Air Force Radar Academy，2009，23（5）：370－372.（in Chinese）

［2］ 顧曉東，袁志勇，付天舒，等.多聲納基陣時差定位精度分析［J］.測控技術(shù).2005，24（4）：69－71.

Gu Xiaodong，Yuan Zhiyong，F(xiàn)u Tianshu，et al.Analysis of the time division of arrival of multiple sonar array［J］.Measurement and control technology，2005，24（4）：69－71.（in Chinese）

［3］ 陳曉忠，梁國龍，王逸林，等.非同步水聲定位技術(shù)及其性能評價［J］.聲學(xué)學(xué)報，2003，28（4）：357－362.

Chen Xiaozhong，Liang Guolong，Wang Yilin，et al.Non－synchronous underwater acoustic location and it’s performance evaluation［J］.Acta Acustica，2003，28（4）：357－362.（in Chinese）

［4］ 鄧勇，周宇.無源時差定位系統(tǒng)的模糊無解問題［J］.國防科技大學(xué)學(xué)報，1999，21（3）：84－87.

DengYong，Zhou Yu.An study of ambiguity and non－solution of passive TDOA location systems［J］.Journal of national university of defense technology，1999，21（3）：84－87.（in Chinese）

［5］ 西北工業(yè)大學(xué)應(yīng)用數(shù)學(xué)系線性代數(shù)教學(xué)組.線性代數(shù)［M］.西北工業(yè)大學(xué)出版社，1999.

［6］ 劉伯勝.分層介質(zhì)中聲波傳播時間和傳播距離間的一個近似公式［J］.船工科技，1977（1）：29－32.

Liu Bosheng.An approximate equation between sound spread time and spread distance in layered medium［J］.Ship industry science and technology，1977（1）：29－32.（in Chinese）

［7］ 王燕，梁國龍.一種適用于長基線水聲定位系統(tǒng)的聲速修正方法［J］.哈爾濱工程大學(xué)學(xué)報，2002，23（5）：32－34.

Wang Yan，Liang Guolong.Correction of sound velocity in long baseline acoustic positioning system［J］.Jounal of Harbin Engineering University，2002，23（5）：32－34.（in Chinese）

［8］ 徐衛(wèi)良，張啟先.機(jī)器人誤差分析的蒙特卡洛方法［J］.機(jī)器人，1988，2（4）：1－6.

Xu Weiliang，Zhang Qixian.Monte－carlo method used on robot error analysis［J］.Robot，1988，2（4）：1－6.（in Chinese）

［9］ 趙永濤，王田苗，孫磊.基于蒙特卡洛方法的導(dǎo)航機(jī)器人的誤差分析［J］.航空制造技術(shù)，2004（3）：52－55.

Zhao Yongtao，Wang Tianmiao，Sun Lei.Error analysis of navigating robot based on monte－carlo method［J］.Aeronautical Manufacturing Technology，2004（3）：52－55.（in Chinese）

［10］ 葉朝謀，俞志強(qiáng).三維時差定位系統(tǒng)的模糊無解分析［J］.現(xiàn)代電子技術(shù)，2005（21）：22－27.

Ye Chaomei，Yu Zhiqiang.Study on ambiguity and non－solution of 3－D TDOA location systems［J］.Modern electronic technology，2005（21）：22－27.（in Chinese）