（91388部隊(duì) 湛江 524022）

1 引言

長基線水聲定位導(dǎo)航系統(tǒng)［1］是水下目標(biāo)試驗(yàn)必備的測量設(shè)備。但是試驗(yàn)過程中在回收應(yīng)答器時(shí)，應(yīng)答器是否在原校準(zhǔn)點(diǎn)位或應(yīng)答器上浮后未能及時(shí)發(fā)現(xiàn)，都將會(huì)發(fā)生應(yīng)答器丟失的安全隱患。開發(fā)本軟件，根據(jù)長基線潛艇導(dǎo)航系統(tǒng)測陣的原理，實(shí)時(shí)解算出待回收應(yīng)答器的點(diǎn)位或運(yùn)動(dòng)軌跡，并給出包括工作母船位置、航跡、航速、應(yīng)答器預(yù)定回收點(diǎn)、再校準(zhǔn)后位置等回收態(tài)勢圖，引導(dǎo)測量船完成應(yīng)答器的回收工作，將大大降低應(yīng)答器丟失的風(fēng)險(xiǎn)。

2 測陣原理

帶DGPS接收機(jī)的工作母船在海底應(yīng)答器斜上方海面按預(yù)定航路勻速航行并記錄軌跡，同時(shí)周期性的向海底應(yīng)答器基陣發(fā)射詢問聲脈沖，通過測量聲信號(hào)在母船和應(yīng)答器之間的雙程傳播時(shí)延Ti，采用橢球面交匯算法即可確定應(yīng)答器的位置［2］。設(shè)有效測陣脈沖個(gè)數(shù)為N，母船以航速（vx，vy）沿航路航行（如圖1），得定位方程組［3］：

式（1）中：i為有效測陣脈沖序號(hào)；（xi，yi，zi）為第i個(gè)應(yīng)答器的空間位置；Ti為目標(biāo)到應(yīng)答器i之間的聲雙程傳播時(shí)延；C為平均聲速；Tdi＝T0i＋T1＋Ti，Tdi是水聲收發(fā)機(jī)發(fā)出詢問聲脈沖到接收到第i個(gè)應(yīng)答器的應(yīng)答聲脈沖為止的總時(shí)延，T0i為第i號(hào)應(yīng)答器的電路延時(shí)，T1為水聲收發(fā)機(jī)的電路延時(shí)（T0i、T1試驗(yàn)前可通過計(jì)量測得）。

圖1 長基線導(dǎo)航定位原理示意圖

母船在多個(gè)位置測量，并記錄該位置坐標(biāo)和水聲收發(fā)機(jī)所測距離，即對(duì)于任何一個(gè)應(yīng)答器都可如式（1）建立相應(yīng)的定位方程，一般來說，聯(lián)立三個(gè)方程就可以解得應(yīng)答器（xi，yi，zi）的精確位置。

3 算法模型

3.1 最速下降法迭代模型［4～6］

水聲收發(fā)機(jī)發(fā)出詢問聲脈沖時(shí)刻母船到應(yīng)答器的距離Rik，即得到下列式（2）所示的等效球面交匯方程組：

（x，y，z）為發(fā)出測距詢問聲脈沖時(shí)刻母船的位置。

先選取一組初值（x0，y0，z0），按最快收斂的方向進(jìn)行迭代計(jì)算，以“失敗—成功”法求最佳步長，從而直指最優(yōu)解，得到滿足一定精度要求的目標(biāo)坐標(biāo)。迭代計(jì)算數(shù)學(xué)模型為

式中，（x，y，z）為本次迭代計(jì)算定位點(diǎn)的坐標(biāo)，（x0，y0，z0）為前次迭代計(jì)算點(diǎn)的坐標(biāo)，α為探索步長。

N為有效應(yīng)答器數(shù)目。

在式（4），（5）中

其中（Xi，Yi，Zi）為應(yīng)答器i坐標(biāo)，D（i，1）＝Rik，母船到應(yīng)答器i的距離。

均方誤差為

在程序中，式（3）的迭代計(jì)算不是無限的迭代計(jì)算下去，一般選用幾次迭代就可滿足測量定位要求［8］。

3.2 脈沖挑選［7～8］

為有效的從帶有脈沖干擾的脈沖序列中挑選出測距信號(hào)的直達(dá)聲，需對(duì)脈沖序列結(jié)構(gòu)和特征的分析，建立一系列判據(jù)，根據(jù)這些判據(jù)為每一個(gè)有效脈沖評(píng)定品質(zhì)因數(shù)，然后依據(jù)品質(zhì)因數(shù)選出最優(yōu)脈沖參與解算。

設(shè)每個(gè)陣元記錄了最多m個(gè)脈沖的信息傳給計(jì)算機(jī)，設(shè)Qai為直達(dá)脈沖品質(zhì)因數(shù)，記f0為信號(hào)載頻，tki為第k幀收到的第i個(gè)脈沖信號(hào)的時(shí)延值，tk-1表示上個(gè)周期選出的最優(yōu)脈沖，Ti表示脈沖的脈寬，則評(píng)定品質(zhì)因數(shù)的判據(jù)如下：

判據(jù)①：制保留每個(gè)脈沖群中的首脈沖，若某脈沖與前一脈沖的時(shí)延差大于200ms，即可認(rèn)為屬于另一脈沖群。

對(duì)于同一周期，當(dāng)

對(duì)于相鄰周期，當(dāng)

時(shí)，令Qai＝Qai＋5。

判據(jù)②：由于信道的影響，聲波的能量將逐漸衰減，一般情況下，直達(dá)聲先于反射聲到達(dá)，因而能量較大。

判據(jù)③：信道串漏脈沖和尖脈沖干擾的瞬時(shí)頻率方差較大，增加方差小的脈沖的品種因數(shù)。

當(dāng)Vki＞Vmax時(shí)，令Qai＝0。

判據(jù)④：以最高航速估算，接收到的脈沖信號(hào)多普勒Δf在一定范圍內(nèi)。當(dāng)時(shí)，令Qai＝Qai-4。

判據(jù)⑤：直達(dá)聲的脈寬不會(huì)過于小。

當(dāng)Tki＜Tmin時(shí)，令Qai＝Qai-4。

判據(jù)⑥：比較當(dāng)前脈沖與其它通道所有脈沖的時(shí)延，若其差值接近于0，則幅度小的脈沖的品質(zhì)因數(shù)置0。如果幅度接近，則脈寬小的脈沖的品質(zhì)因數(shù)置0。

當(dāng)目標(biāo)位于兩陣元中點(diǎn)時(shí)，理想情況下，其時(shí)延也應(yīng)該相等，應(yīng)用判據(jù)⑥必然會(huì)產(chǎn)生誤判［9］。在實(shí)際應(yīng)用中，可以忽略這種誤判，原因是目標(biāo)位于這種特殊位置點(diǎn)的概率很小，而且即使發(fā)生這種情況，可利用冗余陣元參加定位解算。

判據(jù)⑦：由于目標(biāo)是運(yùn)動(dòng)的，如果是直達(dá)脈沖，則在相鄰周期一定的時(shí)延差內(nèi)必定存在與之對(duì)應(yīng)的脈沖，也就是時(shí)延的連續(xù)性。為此，我們將本幀的每個(gè)脈沖與上一幀的選出的脈沖相比，判斷其差值是否在一周期目標(biāo)所能航行的最大距離范圍內(nèi)。

4 誤差分析

對(duì)軟件算法模型及實(shí)際環(huán)境分析可知，軟件誤差主要受以下因素的影響［10］：

1）系統(tǒng)誤差。系統(tǒng)本身的電路時(shí)延造成的，通過計(jì)量得到其時(shí)延值，在解算中消減。

2）GPS定位誤差。由于DGPS本身的定位精度造成的偏差，這里主要表現(xiàn)在監(jiān)測浮標(biāo)的GPS位置數(shù)據(jù)上。實(shí)際應(yīng)用中，我們將其當(dāng)作真實(shí)值來參與解算，但其誤差存在必然性對(duì)定位的準(zhǔn)確度造成的影響是不可忽略的，因此選用高精度的DGPS成為降低該誤差的一種途徑。

3）測時(shí)誤差。測距儀通過測量“詢問-應(yīng)答”時(shí)間與平均聲速的乘積得到母船與應(yīng)答器的距離，測時(shí)誤差的影響取決于應(yīng)答器到兩測量點(diǎn)的距離差與測量點(diǎn)間距之比。

4）聲速裝訂誤差。由于水文條件的影響，且海洋環(huán)境的不確定性，參與計(jì)算的聲速值是一個(gè)折算后的平均值，由于聲速貫穿整個(gè)測量和解算的過程，故其精度對(duì)整個(gè)導(dǎo)引軟件定位精度都有影響。

5 實(shí)驗(yàn)室功能性測試和驗(yàn)證

在實(shí)驗(yàn)室讀取海上試驗(yàn)數(shù)據(jù)，裝訂當(dāng)時(shí)試驗(yàn)時(shí)聲速、應(yīng)答器延時(shí)、聲頭深度、測陣后待回收應(yīng)答器經(jīng)緯度及深度值，對(duì)軟件進(jìn)行公正性測試，運(yùn)行圖如圖2，圖中顯示了待回收應(yīng)答器實(shí)時(shí)解算位置及運(yùn)動(dòng)軌跡，工作母船當(dāng)前位置及軌跡，并以圖表形式實(shí)時(shí)顯示斜距值，在列表中顯示工作母船和待回收應(yīng)答器當(dāng)前及歷史經(jīng)緯度值、斜距值及航速。

圖2 應(yīng)答器回收導(dǎo)引軟件運(yùn)行圖

6 結(jié)語

通過實(shí)驗(yàn)室功能性測試和驗(yàn)證，本軟件能夠?qū)崿F(xiàn)實(shí)時(shí)解算并顯示出當(dāng)前應(yīng)答器位置及斜距值，同時(shí)在界面上還顯示工作母船當(dāng)前位置及運(yùn)動(dòng)軌跡、航速應(yīng)答器運(yùn)動(dòng)軌跡，用列表的形式顯示母船、應(yīng)答器經(jīng)緯度坐標(biāo)、母船航速及斜距值并保存文件以備后期處理，可以完成應(yīng)答器回收導(dǎo)引的功能。

［1］吳永亭，周興華，楊龍.水下聲學(xué)定位系統(tǒng)及其應(yīng)用［J］.海洋測繪，2003：18-2.

［2］王先華.長基線水下導(dǎo)航定位系統(tǒng)測陣校陣及系統(tǒng)集成［D］.哈爾濱：哈爾濱工程大學(xué)碩士學(xué)位論文，2004：35-38.

［3］李莉.長基線陣測陣校陣技術(shù)研究［D］.哈爾濱：哈爾濱工程大學(xué)碩士學(xué)位論文，2007：44-48.

［4］陳曉忠，梁國龍，等.非同步水聲定位技術(shù)及其性能評(píng)價(jià)［J］.聲學(xué)學(xué)報(bào)，2003，28（4）：358-360.

［5］封金星，丁士圻，惠俊英.水下運(yùn)動(dòng)目標(biāo)長基線定位解算研究［J］.聲學(xué)學(xué)報(bào)，1996，18（5）：22-23.

［6］于平.兩種水聲導(dǎo)航定位算法的應(yīng)用與分析［J］.靶場試驗(yàn)與管理，2004，1：15-18.

［7］于平，吳波，謝勝.水聲定位信號(hào)脈沖挑選方法［J］.艦船電子工程，2012，32（2）：123.

［8］王正明，易東云，等.測量數(shù)據(jù)建模與參數(shù)估計(jì)［M］.國防科技大學(xué)出版社，1996：158-159.

［9］惠娟.水下多目標(biāo)跟蹤定位算法研究［D］.哈爾濱：哈爾濱工程大學(xué)碩士學(xué)位論文，2004：8-18.

［10］劉百峰，楊志權(quán)，趙珩.一種適用長基線系統(tǒng)絕對(duì)測陣的新航路［J］.艦船電子工程，2011，31（6）：173.