水下高速直線(xiàn)運(yùn)動(dòng)目標(biāo)的雙參考源航跡估計(jì)方法研究

王謙，侯宏，尹韶平，陳志菲（.西北工業(yè)大學(xué)航海學(xué)院，陜西西安7007；.中國(guó)船舶重工集團(tuán)公司第705研究所，陜西西安70075；.中國(guó)科學(xué)院聲學(xué)研究所，北京0090）

?

水下高速直線(xiàn)運(yùn)動(dòng)目標(biāo)的雙參考源航跡估計(jì)方法研究

王謙1，侯宏1，尹韶平2，陳志菲3
（1.西北工業(yè)大學(xué)航海學(xué)院，陜西西安710072；2.中國(guó)船舶重工集團(tuán)公司第705研究所，陜西西安710075；3.中國(guó)科學(xué)院聲學(xué)研究所，北京100190）

摘要：針對(duì)噪聲源識(shí)別中水下近場(chǎng)高速直線(xiàn)運(yùn)動(dòng)目標(biāo)的航跡估計(jì)問(wèn)題，將雙脈沖源方法擴(kuò)展到三維空間中。該方法利用公垂線(xiàn)最小化，將基線(xiàn)定位方法估計(jì)的兩個(gè)連續(xù)波脈沖參考源的定位結(jié)果投影到二維平面內(nèi)，并結(jié)合已知的雙參考源間距等信息構(gòu)造航跡參數(shù)的目標(biāo)函數(shù)，進(jìn)而利用遺傳算法估計(jì)三維空間內(nèi)直線(xiàn)運(yùn)動(dòng)目標(biāo)在各脈沖輻射時(shí)刻下的位置和航速。仿真結(jié)果表明，0 dB條件下當(dāng)航向角和俯仰角均在［-40°，40°］內(nèi)時(shí)，其估計(jì)偏差均在3. 8°以?xún)?nèi)，而正橫距離和深度差的估計(jì)偏差在5. 9 m以?xún)?nèi)。水庫(kù)實(shí)驗(yàn)進(jìn)一步驗(yàn)證了該方法的有效性。

關(guān)鍵詞：聲學(xué)；參數(shù)估計(jì)；直線(xiàn)航跡估計(jì)；雙脈沖源；基線(xiàn)定位

侯宏（1966—），男，教授，博士生導(dǎo)師。E-mail：houhong@ nwpu. edu. cn

0 引言

對(duì)水下航行器等高速水下運(yùn)動(dòng)目標(biāo)進(jìn)行噪聲源識(shí)別時(shí)，需要預(yù)估其航跡以便對(duì)輻射噪聲解多普勒［1 -3］，從而將動(dòng)態(tài)噪聲信號(hào)還原成靜止聲源進(jìn)行處理。精確的航跡估計(jì)不僅有利于解多普勒，也便于根據(jù)加裝在水下航行器上的已知參考源位置，解算噪聲源識(shí)別結(jié)果在水下航行器的對(duì)應(yīng)位置。

在運(yùn)動(dòng)目標(biāo)噪聲源識(shí)別時(shí)，陣列與航跡相距較近，此時(shí)多采用時(shí)頻分析方法或基線(xiàn)定位法估計(jì)航跡［4］。對(duì)于勻速直線(xiàn)運(yùn)動(dòng)目標(biāo)，其時(shí)頻曲線(xiàn)參數(shù)與航跡參數(shù)間有較明確的對(duì)應(yīng)關(guān)系，因此時(shí)頻分析方法，如利用Hough變換提取LOFAR譜圖中拋物線(xiàn)參數(shù)，或者估計(jì)瞬時(shí)頻率，可以反推航跡信息。時(shí)頻分析方法僅需要單個(gè)陣元，但是信噪比門(mén)限較高，非勻速直線(xiàn)運(yùn)動(dòng)情況下時(shí)頻參數(shù)的解讀較為困難。基線(xiàn)定位方法根據(jù)參考源輻射信號(hào)到達(dá)各陣元的時(shí)延差來(lái)定位目標(biāo)，它又分為同步脈沖和異步脈沖定位法。前者需要參考源和接收端的時(shí)鐘同步到世界標(biāo)準(zhǔn)時(shí)間，后者與三子陣法原理相似，通過(guò)構(gòu)造包含目標(biāo)坐標(biāo)和輻射時(shí)刻的超定方程進(jìn)行定位［5 -6］。同步脈沖法對(duì)時(shí)鐘要求較高，異步脈沖法降低了硬件要求，但是定位精度相對(duì)較差。整體上基線(xiàn)定位方法相對(duì)成熟，實(shí)際工程中應(yīng)用廣泛。

現(xiàn)有研究大多將水下航行器等運(yùn)動(dòng)目標(biāo)看成分布式點(diǎn)源，在二維平面內(nèi)作勻速直線(xiàn)運(yùn)動(dòng)。而實(shí)際中由于測(cè)噪船的漂移和水聽(tīng)器陣列的吊放深度等因素，水下航行器航跡與陣列的幾何關(guān)系需要在三維空間中描述，并且水下航行器實(shí)航時(shí)是一個(gè)加速、勻速到減速的運(yùn)動(dòng)過(guò)程。文獻(xiàn)［7］中雙脈沖參考源（DPS）方法通過(guò)在目標(biāo)上加裝兩個(gè)脈沖參考源，增加了雙源間距這個(gè)條件，從而提高了二維空間中異步脈沖法的定位精度，并且該方法適用于非勻速直線(xiàn)運(yùn)動(dòng)目標(biāo)。本文重新定義了DPS方法中的雙源間距（DSS）算子和聲程差（PD）算子，并增加公垂線(xiàn)長(zhǎng)度（CPL）最小化條件，從而將DPS方法擴(kuò)展到了三維空間。

1 算法的提出

1. 1 問(wèn)題的提出

如圖1所示，假設(shè)陣列位于Oxz平面內(nèi)，空間直線(xiàn)航跡PQ在Oxy平面的投影PU與x軸的夾角為航向角θ（取值范圍（-π/2，π/2）），PQ與PU的夾角為俯仰角φ（取值范圍（-π/2，π/2）），PU與y軸交點(diǎn)U的y坐標(biāo)為正橫距離R0（R0＞0），航跡與U點(diǎn)的垂直距離UQ為深度差H0，這樣航跡參數(shù)（θ，φ，R0，H0）就確定了空間直線(xiàn)航跡。本文假定目標(biāo)在水下作直線(xiàn)運(yùn)動(dòng)，平面陣從測(cè)噪船垂直吊放入水，由于測(cè)噪船起伏漂移等因素使得航跡與陣列出現(xiàn)圖1所示情形。

圖1 加裝雙參考源的直線(xiàn)運(yùn)動(dòng)目標(biāo)Fig. 1 Linear moving target with double reference source

假設(shè)目標(biāo)上有兩個(gè)間距為r0的同時(shí)輻射連續(xù)波脈沖信號(hào)的參考源A和B，輻射周期均為T(mén)c，載頻不同且都位于噪聲源識(shí)別分析頻帶以外。對(duì)接收信號(hào)濾波后可分離兩個(gè)連續(xù)波脈沖信號(hào)。對(duì)于某個(gè)連續(xù)波脈沖［8］有

式中：c為聲速；M為陣元數(shù)。Sm=（xm，ym，zm）和tm分別是第m個(gè)陣元的空間坐標(biāo)和其接收到的脈沖信號(hào)前沿到達(dá)時(shí)刻，（xs，ys，zs）和ts則分別是聲源輻射該脈沖時(shí)的空間坐標(biāo)和輻射時(shí)刻。（1）式中，xs、ys、zs和ts未知，陣元坐標(biāo)Sm已知，脈沖前沿到達(dá)時(shí)刻tm可由接收信號(hào)的時(shí)域處理得到。當(dāng)M＞4時(shí)，（1）式成為超定方程組，從而可以解得聲源坐標(biāo)（xs，ys，zs）和輻射時(shí)刻ts.與文獻(xiàn)［7］類(lèi)似，利用不同時(shí)刻下基線(xiàn)定位方法得到的參考源的定向結(jié)果和輻射時(shí)刻，可以構(gòu)造目標(biāo)函數(shù)來(lái)估計(jì)航跡參數(shù)（θ，φ，R0，H0）。

1. 2 DPS方法

脈沖基線(xiàn)定位方法的測(cè)向精度較高，而測(cè)距結(jié)果波動(dòng)較大。DPS方法僅利用其測(cè)向結(jié)果，使得基線(xiàn)定位結(jié)果分布于坐標(biāo)原點(diǎn)和真實(shí)航跡直線(xiàn)所確定的航跡平面內(nèi)，因此，該擬合平面包含的角度及深度等信息計(jì)算結(jié)果較為精確。由于基線(xiàn)定位算法的測(cè)距精度較差，導(dǎo)致水平定位結(jié)果精度較差。本文中DPS方法結(jié)合已知的雙參考源間距等約束條件來(lái)改善其水平定位結(jié)果，構(gòu)造了DSS和PD兩個(gè)目標(biāo)函數(shù)，并取它們的乘積作為DPS的目標(biāo)函數(shù)來(lái)估計(jì)航跡參數(shù)。

1. 2. 1 將DSS算子擴(kuò)展到三維空間中

圖2中假定Ek和Fk是第k個(gè)雙脈沖周期下兩個(gè)參考源的基線(xiàn)定位結(jié)果，二維空間內(nèi)OEk和OFk這兩條定向直線(xiàn)與航跡相交，兩個(gè)交點(diǎn)的間距可直接作為雙參考源間距r0的估計(jì)值^r0k.三維空間中定向直線(xiàn)OEk與航跡通常是異面直線(xiàn)，同樣，OFk和航跡直線(xiàn)也是如此；并且，多個(gè)脈沖時(shí)刻下估計(jì)結(jié)果在任意平面的投影點(diǎn)E1、E2、E3…通常不在一條直線(xiàn)上，同樣，F(xiàn)1、F2、F3…也如此。為此，圖2中取定向直線(xiàn)與航跡間的公垂線(xiàn)在航跡上的垂足Ak和Bk作為Ek和Fk在航跡上的投影點(diǎn)，進(jìn)而得到雙參考源間距估計(jì)值^r0k.

以Ak為例，要確定其坐標(biāo)需要航跡PQ和定向直線(xiàn)OEk的方向向量。由圖1的幾何關(guān)系可得航跡PQ 與Oyz平面交點(diǎn)Q =［0，R0，H0］，那么P點(diǎn)坐標(biāo)為

這樣航跡PQ的方向向量為

對(duì)于第k個(gè)雙脈沖周期，假設(shè)由（1）式估計(jì)的脈沖聲源坐標(biāo)為Ek=［xEk，yEk，zEk］，則定向直線(xiàn)OEk的方向向量為

由（3）式和（4）式可得Ak點(diǎn)坐標(biāo)為

將（4）式中Ek的坐標(biāo)替換為Fk，即得另一垂足Bk的坐標(biāo)，因此

式中：AkBk為Ak和Bk之間的距離，它是航跡參數(shù)（θ，φ，R0，H0）的函數(shù)；K為雙脈沖周期數(shù)。顯然，通過(guò)（6）式可將二維空間的DSS算子擴(kuò)展到了三維空間，即

1. 2. 2 將PD算子擴(kuò)展到三維空間中

對(duì)于某個(gè)頻率下的連續(xù)波脈沖聲源，其相鄰兩個(gè)周期的波形存在圖3所示的幾何關(guān)系。

圖3 相鄰連續(xù)波脈沖的聲程差示意圖Fig. 3 Schematic diagram of acoustic path difference of adjacent CW pulse

對(duì)于坐標(biāo)為Sm的第m個(gè)陣元的接收信號(hào)，A聲源的相鄰兩個(gè)脈沖波形前沿時(shí)刻的差值Tm可由（1）式估計(jì)的tm得到，即Tm= tm2- tm1.由此得到相鄰脈沖的聲程差為

式中：SmA1和SmA2分別是Sm到A1和A2的距離。這里Ak的定義與圖2中相同，即定向直線(xiàn)OEk與航跡PQ的公垂線(xiàn)在航跡上的垂足。不同陣元下對(duì)每組相鄰脈沖均可得到類(lèi)似（8）式的等式，這樣就把PD算子轉(zhuǎn)換到三維空間

式中：Tmk為第m個(gè)陣元接收信號(hào)中第k和k + 1個(gè)脈沖之間的時(shí)間間隔，即Tmk= tm（k +1）- tm（k）.

1. 2. 3 增加CPL最小化的限制條件

在真實(shí)航跡所在垂直平面，即圖1中的陰影平面內(nèi)，平行于真實(shí)航跡的一些直線(xiàn)對(duì)應(yīng)的DSS和PD算子值較為接近，從而難以確定航跡參數(shù)。因此公垂線(xiàn)最小化條件可以規(guī)避該情形，即最小化為

式中：dEk= EkAk，dFk= FkBk，即兩條公垂線(xiàn)的長(zhǎng)度。其中［9］

上述各變量的定義與（4）式和（5）式中對(duì)應(yīng)變量相同，類(lèi)似的可以得到dFk的表達(dá)式。

1. 2. 4 三維空間的DPS方法

將上述3個(gè)算子結(jié)合在一起就得到三維空間的DPS方法

式中：fPD1和fPD2是分別針對(duì)兩個(gè)連續(xù)波脈沖聲源的PD算子。采用遺傳算法使上述目標(biāo)函數(shù)最小化即可得到航跡估計(jì)結(jié)果。另外，（11）式采用相加而不是文獻(xiàn)［7］中相乘的方式得到DPS算子是因?yàn)镃PL算子的值有可能很小，容易使相乘方式構(gòu)造的DPS算子出現(xiàn)局部極小值。

為了準(zhǔn)確估計(jì)相鄰脈沖間的目標(biāo)航速，需要確定各脈沖輻射前沿時(shí)刻下雙參考源的坐標(biāo)。以第一參考源為例，其坐標(biāo)Ak可由航跡估計(jì)結(jié)果（θ，φ，R0，H0）代入（5）式得到，但是它對(duì)應(yīng)的輻射時(shí)刻tE（k）通常不等于（1）式中的ts，因?yàn)閠s對(duì)應(yīng)的是脈沖基線(xiàn)定位結(jié)果Ek的輻射時(shí)刻。為此，需要根據(jù)各陣元的脈沖前沿到達(dá)時(shí)刻tEm（k）（m =1，…，M），以及Ak來(lái)反推tE（k），即

依此方法，可以得到Bk對(duì)應(yīng)的輻射時(shí)刻tF（k）.這樣相鄰兩個(gè)脈沖信號(hào)間目標(biāo)的航速為

1. 3 理論分析

三維空間的DPS方法相當(dāng)于將脈沖定位結(jié)果對(duì)應(yīng)的定向直線(xiàn)OE和OF投影到平面OPQ上，進(jìn)而在該平面內(nèi)利用DSS和PD算子確定航跡參數(shù)，因此其核心仍是二維平面內(nèi)的航跡擬合。在二維投影平面內(nèi)估計(jì)正橫距離R＇0和航向角θ＇后，可根據(jù)平面OPQ在三維空間的位置得到航跡參數(shù)（θ，φ，R0，H0）.

圖4給出了雙參考源A和B的兩組脈沖定位結(jié)果投影到平面OPQ后的幾何關(guān)系。假定目標(biāo)從左向右作直線(xiàn)運(yùn)動(dòng)，r1和r3為第一參考源相鄰兩個(gè)脈沖對(duì)應(yīng)的定向直線(xiàn)，r2和r4為第二參考源相鄰兩個(gè)脈沖對(duì)應(yīng)的定向直線(xiàn)。由基線(xiàn)定向結(jié)果可得r1、r2和r3、r4所夾角度分別為β1和β2，它們正對(duì)的三角形邊長(zhǎng)均為雙源間距r0.

圖4 DPS方法原理示意圖Fig. 4 Schematic diagram of DPS method

由三角形三邊關(guān)系有

式中：r1、r2、r3、r4為未知量。若以三角形OA1B1為參考，OA2B2的對(duì)應(yīng)邊長(zhǎng)和夾角可改寫(xiě)為r3= r1+Δa，r4= r2+Δb，β2=β1+Δβ，其中Δa、Δb分別為同一參考源在不同位置的距離估計(jì)差，Δβ為兩次估計(jì)夾角的差，代入（14）式得

整理得

當(dāng)有多組雙脈沖定向結(jié)果時(shí)，即有多個(gè)三角形OAiBi，i =1，…，K，可利用（16）式構(gòu)建r1和r2的二元二次方程組

二維空間內(nèi)DPS方法就是要尋找一條航跡直線(xiàn)與rk所在的定向直線(xiàn)相交，并盡可能保證AkBk= r0.而對(duì)于任意給定的航跡直線(xiàn)，（17）式中獨(dú)立的未知量?jī)H有r1和r2.

當(dāng)K =4時(shí)，通過(guò)簡(jiǎn)單的變量代換可得

式中：Ui=Δai-Δbicos βi；Vi=Δbi-Δaicos βi；Wi=cos β0- cos βi；Xi=（Δa+Δb）/2.

（18）式、（19）式成立的條件是分母不為0，整理后得Wi/ Wj≠Ui/ Uj，Wi/ Wj≠Vi/ Vj，其中i≠j，這在實(shí)際中容易滿(mǎn)足。尤其是當(dāng)K＞4時(shí)，可采用非線(xiàn)性最小二乘求解（17）式得到r1和r2，從而得到二維空間航跡參數(shù)（θ＇，R＇0），并換算到三維空間中的（θ，φ，R0，H0）.

2 仿真

本文著重給出不同信噪比下，以及不同航向角和俯仰角下DPS方法在三維空間中的航跡參數(shù)估計(jì)性能。在圖1所示條件下，假定采用M = 9元均勻十字陣，其中垂直方向陣元數(shù)為2個(gè)，陣元間距d =0. 3 m.兩個(gè)間距r0= 0. 5 m的連續(xù)波脈沖聲源載頻f0=［6 000 Hz，7 000 Hz］，脈沖周期均為T(mén)c= 200 ms.水中聲速c = 1 500 m/ s，目標(biāo)作航速υ= 20 m/ s的勻速直線(xiàn)運(yùn)動(dòng)，航跡參數(shù)（θ，φ，R0，H0）= ［10°，15°，40 m，0 m］，從中截取K =7組雙脈沖進(jìn)行航跡參數(shù)估計(jì)。圖5給出了信噪比SNR =0 dB時(shí)的航跡估計(jì)結(jié)果，其中信噪比定義為雙脈沖源的能量和與背景噪聲能量的比值，基線(xiàn)定位方法是第一參考源定位結(jié)果則經(jīng)過(guò)了Kalman平滑處理。圖5 中DPS估計(jì)的航跡參數(shù)為（θ，φ，R0，H0）= ［10. 42°，15. 07°，39. 17 m，- 0. 01 m］，即估計(jì)偏差為［0. 42°，0. 07°，0. 83 m，0. 01 m］.

圖5 0 dB時(shí)DPS方法的航跡估計(jì)結(jié)果Fig. 5 Track estimation result of DPS method at 0 dB

圖6 不同信噪比下DPS方法的航跡參數(shù)估計(jì)性能Fig. 6 Track estimation performance of DPS method at different SNRs

相同仿真條件下，圖6給出了100次Monte Carlo實(shí)驗(yàn)后，不同信噪比下DPS方法的航跡參數(shù)估計(jì)結(jié)果的平均絕對(duì)偏差。隨著信噪比的提高，DPS方法的航跡估計(jì)性能不斷提高，其中航向角θ和正橫距離R0的估計(jì)偏差大于俯仰角φ和深度差H0的估計(jì)偏差。這與平面陣對(duì)垂直于其所在平面的縱深方向的分辨能力較弱有關(guān)，因而航向角θ和正橫距離R0也更易受到基線(xiàn)定位方法測(cè)向偏差，以及聲程差估計(jì)偏差的影響。圖7和圖8分別給出了0 dB下100次Monte Carlo實(shí)驗(yàn)后不同航向角和俯仰角下航跡參數(shù)的平均估計(jì)偏差，其中航向角θ和正橫距離R0的估計(jì)偏差仍大于相應(yīng)的俯仰角φ和深度差H0的估計(jì)偏差。另一方面，隨著航向角θ和俯仰角φ的增大，即航跡直線(xiàn)與平面陣所在平面夾角增大時(shí)，航跡估計(jì)性能逐漸下降。在0 dB下，當(dāng)航向角θ和俯仰角φ均在［-40°，40°］內(nèi)時(shí)，θ和φ的估計(jì)偏差在3. 8°以?xún)?nèi)，而R0和H0的估計(jì)偏差在5. 9 m以?xún)?nèi)。

圖7 不同航向角下DPS方法的航跡估計(jì)性能Fig. 7 Track estimation performance of DPS method at different course angles

3 水庫(kù)實(shí)驗(yàn)驗(yàn)證

受限于實(shí)驗(yàn)成本，下面以水庫(kù)中自由落體目標(biāo)的航跡估計(jì)為例來(lái)驗(yàn)證DPS方法。圖9給出了實(shí)驗(yàn)示意圖，其中水聽(tīng)器陣列與運(yùn)動(dòng)目標(biāo)均在垂直面內(nèi)，即（θ，R0）確定了航跡。水聽(tīng)器陣列為6 m長(zhǎng)的11元均勻十字陣，其中垂直方向陣元數(shù)為2個(gè)。運(yùn)動(dòng)目標(biāo)為從水面開(kāi)始作自由落體的金屬殼體，其上加裝兩個(gè)間距r0= 0. 4 m的脈沖參考源，連續(xù)波脈沖載頻為［7 000 Hz，10 000 Hz］，脈沖周期均為200 ms.

圖8 不同俯仰角下DPS方法的航跡估計(jì)性能Fig. 8 Track estimation performance of DPS method at different pitch angles

圖10給出了兩次獨(dú)立實(shí)驗(yàn)的航跡估計(jì)結(jié)果。兩次實(shí)驗(yàn)中僅改變了陣列與航跡的間距R0，其他實(shí)驗(yàn)參數(shù)和環(huán)境相同。兩次實(shí)驗(yàn)的真實(shí)航跡參數(shù)（θ，R0）分別為（0°，31. 3 m）和（0°，23 m）。DPS方法的航跡估計(jì)結(jié)果分別為（- 5. 7°，32. 5 m）和（-5. 1°，21. 2 m），即正橫距離R0估計(jì)偏差在2 m以?xún)?nèi)，航向角θ估計(jì)偏差在- 6°以?xún)?nèi)。航向角估計(jì)偏差過(guò)大可能是由于錨定的浮船仍在水面上緩慢旋轉(zhuǎn)漂移，使得垂直吊放的線(xiàn)陣有一定傾斜度，因此上述實(shí)驗(yàn)中真實(shí)航向角可能在-5°左右。

4 結(jié)論

本文在文獻(xiàn)［7］的雙參考源航跡估計(jì)方法基礎(chǔ)上，將原有的雙源間距和聲程差目標(biāo)函數(shù)擴(kuò)展到三維空間，并引入CPL最小化目標(biāo)函數(shù)，從而將二維空間中的DPS方法推廣到了三維空間中。DPS方法相當(dāng)于利用三維空間中基線(xiàn)定位方法的定向結(jié)果，以及已知的雙參考源間距等信息構(gòu)造航跡參數(shù)的最優(yōu)化目標(biāo)函數(shù)，進(jìn)而通過(guò)遺傳算法估計(jì)近場(chǎng)直線(xiàn)運(yùn)動(dòng)目標(biāo)航跡。所得到的航跡信息包含各脈沖輻射時(shí)刻下目標(biāo)的三維坐標(biāo)，因而也可以得到各時(shí)段下直線(xiàn)運(yùn)動(dòng)目標(biāo)的航速。仿真和水庫(kù)實(shí)驗(yàn)驗(yàn)證了DPS方法的有效性。

圖9 水庫(kù)航跡估計(jì)實(shí)驗(yàn)示意圖Fig. 9 Schematic diagram of track estimation test in reservoir

圖10 兩次水庫(kù)實(shí)驗(yàn)的航跡估計(jì)結(jié)果Fig. 10 Two result of track estimation test in reservoir

參考文獻(xiàn)（References）

［1］ Kim S M，Choi Y，Lim Y K. Nearfield noise identification with a passive array in the underwater anechoic basin at KRISO［C］∥OCEANS 2003. San Diego，CA，US：IEEE，2003：1389 -1394.

［2］ Yoon J H，Kim D Y，Bae S H，et al. Joint initialization and tracking of multiple moving objects using Doppler information［J］. IEEE Transactions on Signal Processing，2011，59（7）：3447 -3452.

［3］ Tegborg V，Pettersson M I，Claesson I. Experimental results of passive imaging of moving continuous broadband sound sources within a sensor field［J］. IEEE Journal of Oceanic Engineering，2011，36（1）：25 -36.

［4］ 董利英，陳志菲.近場(chǎng)條件下脈沖基線(xiàn)定位性能影響因素分析［J］.魚(yú)雷技術(shù)，2014，22（2）：100 -104. DONG Li-ying，CHEN Zhi-fei. Influencing factors of pulse baseline positioning performance in near field［J］. Torpedo Technology，2014，22（2）：100 -104.（in Chinese）

［5］ Hu B，Yang D，Shi S，et al. Underwater moving noise source identification based on hybrid near-field acoustical holography ［J］. Science China：Physics，Mechanics and Astronomy，2011，54（9）：1570 -1578.

［6］ Mei J D，Sheng X L，Guo L X，et al. Sea trial research on multisource underwater acoustic image measurement technology based on the sparse array［C］∥OCEANS 2010. Seattle，WA：IEEE，2010.

［7］ Chen Z，Hou H，Yang J，et al. Linear track estimation using double pulse sources for near-field underwater moving target［J］. Journal of Marine Science and Application，2013，12（2）：240 -244.

［8］ 陳曉忠，梁國(guó)龍，王逸林，等.非同步水聲定位技術(shù)及其性能評(píng)價(jià)［J］.聲學(xué)學(xué)報(bào)，2003，28（4）：357 -362. CHEN Xiao-zhong，LIANG Guo-long，WANG Yi-lin，et al. Nonsynchronous underwater acoustic location and it’s performance evaluation［J］. Acta Acustica，2003，28（4）：357 - 362.（in Chinese）

［9］ 董玉久，潘秀英，楊欣欣.兩異面直線(xiàn)公垂線(xiàn)垂足位置的計(jì)算方法［J］.哈爾濱科學(xué)技術(shù)大學(xué)學(xué)報(bào)，1990，14（1）：84 -88. DONG Yu-jiu，PAN Xiu-ying，YANG Xin-xin. A calculation method for finding the foot of common perpendicular of lines on different places［J］. Journal of Harbin University of Science & Technology，1990，14（1）：84 -88.（in Chinese）

Research on Double Reference Source Estimation Method for 3-dimension Space of Near-field Underwater Moving Target with Linear Track

WANG Qian1，HOU Hong1，YIN Shao-ping2，CHEN Zhi-fei3
（1. School of Marine Science and Technology，Northwestern Polytechnical University，Xi’an 710072，Shaanxi，China；2. The 705 Research Institute，China Shipbuilding Industry Corporation，Xi’an 710075，Shaanxi，China；3. Institute of Acoustics，Chinese Academy of Sciences，Beijing 100190，China）

Abstract：For the identification of noise source，the track information could be applied to Doppler compensation for the recorded noise emitted from moving vehicle. Linear track estimation of near-field underwater moving target by double pulse source（DPS）has been used successfully in 2-dimension. The use of DPS is extended to 3- dimension space in order to linearly track the target. The positioning result of double reference source obtained by the baseline positioning method is projected to a certain plane under the limitation of common perpendicular minimum length. The objective function of DPS is constructed with the known distance of double reference source. The track information in 3-dimension space is obtained by using genetic algorithm. The simulation results show that the estimated errors of azimuth and pitch angles are less than 3. 8°when their trace line lies under［-40°，40°］in case of 0 dB noise level. The estimated errors of abeam distance and depth difference are kept less than 5. 9 m at the same time.

Key words：acoustics；parameter estimation；linear track estimation；double reference source；baseline positioning

中圖分類(lèi)號(hào)：TB568

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1000-1093（2016）04-0677-07

DOI：10. 3969/ j. issn. 1000-1093. 2016. 04. 015

收稿日期：2015-05-20

基金項(xiàng)目：海軍裝備預(yù)先研究項(xiàng)目（4011003010201）

作者簡(jiǎn)介：王謙（1982—），男，高級(jí)工程師，博士研究生。E-mail：king2397@126. com；