(海軍潛艇學院 青島 266071)
航空反潛憑借其機動快速、安全性高等特點在反潛作戰中發揮著重要作用。在搜索潛艇過程中,反潛機可以使用聲納浮標、吊放聲納等設備。主動聲納浮標雖具有探測精度高的特點,但是容易被敵潛艇偵察到進而進行規避,對抗反潛機偵測。被動聲納浮標通常只能進行定向,無法定位,但是憑借其隱蔽性高效性的特點在反潛對抗中發揮重要作用。在我國海軍奔向深藍的征途中,執行遠海反潛任務時,反潛機攜帶聲納浮標能迅速到達指定海區完成任務,可將機動性的優勢轉化為戰場優勢。
假設目標位置與聲納浮標位置關系如圖1所示。

圖1 陣元空間關系示意圖
M點為敵潛艇位置,其坐標為M(x,y,z),反潛機投放的被動聲納浮標坐標分別為A(x0,y0,z0),B(x1,y1,z1),C(x2,y2,z2),D(x3,y3,z3),在陸地上調節水聽器沒入海面深度,使四個水聽器具有相同的縱坐標,則聲納浮標在一定空間位置內可認為四個聲納浮標縱坐標均為0。目標到達水聽器距離為s0,s1,s2,s3,發出噪聲信號到達水聽器時間分別為t0,t1,t2,t3,則可近似可得[1~3]:

令τ10=t1-t0,τ20=t2-t0,τ30=t3-t0,則有:

已知時延差,取一點為坐標原點,其他三個傳感器空間坐標為

第i個傳感器到參考傳感器的時延差為τi,τ0為目標到參考傳感器的時間,目標坐標為xT=[xT,yT,zT]T,則有:

其中‖x0‖2表示x的L2范數。
由空間幾何關系知:

寫成矩陣形式則有:

RM×N表示M個元組的集合,每個元組有N個元素。對于M=3的情況,A是滿秩的,因此可以利用式(5)得:

代入式(3)得

令

則式(12)變為

上式為變量s的一個標量二次方程。此處的s為目標到參照傳感器之間的距離,且為正數。則由此可以得出目標的距離,將其代入式(11)中可以得到目標方位坐標[4~5]。
以三元直線陣為例,設目標距離參考陣元為r,陣元間距為d,目標與陣元直線夾角為θ,則遠場下,距離r為

估計量的方差:

其中l=dcosθ稱為等效基陣的半長,Carter曾證明過其接近克拉美羅下界[6~7],式(16)開方可以寫成:

該式等價于:

由誤差分析知,時延估計精度對目標定位精度影響較大。原始信號經濾波后進行數據采集,采用互相關法估計時延差。數據采集速度受硬件條件制約,因此可以采用插值法增加數據點數,提高時延估計精度。也可以通過提高陣元間距來提高時延估計精度[8~9]。
對時延估計的仿真如下:
截取頻率為8000Hz的wav音頻信號作為目標噪音信號,并由Matlab采集。構造0.6506s延遲信號,對兩路信號做互相關分析,并對兩路信號進行Lagrange差值處理后進行互相關時延估計,分析兩次時延估計值。仿真結果如下:

圖2 無插值下互相關時延估計

圖3 無插值下相關峰示意圖

圖4 Lagrange插值后互相關時延估計

圖5 Lagrange插值后相關峰示意圖
在數據采集過程中,設采樣頻率為Fs,則相鄰點的時間間隔為1/Fs,時延估計間隔最大誤差為1/2Fs。由Runge定理知,當插值點過多時誤差增大,因此取四點內插較為理想。由仿真結果可知,在無內插情況下,相關峰峰值位于0.6507s處,采用內插后,相關峰回歸到0.6506s位置,時延估計精度提高。
本文對四元被動聲納浮標陣的定位原理進行了闡述,分析了影響該方法的主要誤差因素,得出通過插值法可以提高時延估計精度進而提高定位精度的結論,并對插值法進行了仿真。本文所采用的仿真數據都是在理想條件下進行的仿真。在未來反潛作戰中,基于隱蔽性的反潛搜索將發揮越來越重要的作用。聲納浮標被動工作可以達到隱蔽性要求,但是由于其位于海表,受海表各種反射信號干擾較大,定位精度不高[10]。通過被動定位陣能夠抑制部分反射信號,減少干擾,提高定位精度,對反潛搜索過程有一定的幫助[11]。
[1]王新勇.互譜法被動測距改進研究[D].哈爾濱:哈爾濱工程大學,2004.
[2]劉聰鋒.無源定位與跟蹤[M].西安:西安電子科技大學出版社,2011:27-32.
[3]熊輝,雷斌.被動定位系統中時延估計算法的仿真分析[J].技術平臺,2011(3):47-49.
[4]王昭,李宏,趙俊渭,等.空氣聲被動定位的誤差分析[J].應用聲學,2000,19(2):39-43.
[5]毛衛寧.水下被動定位方法回顧與展望[J].東南大學學報(自然科學版),2001,31(6):129-132.
[6]王成,李少洪,王鑫全,等.測時差被動定位算法的研究[J].系統工程與電子技術,2001,23(11):9-12.
[7]馬雯,黃建國.廣義相關時延估計在被動定位系統中的應用研究[J].探測與控制學報,2000,22(3):51-54.
[8]殷世民,付小寧,劉上乾,等.對固定平臺紅外單站被動定位技術研究[J].光子學報,2004,33(2):237-239.
[9]崔瑋瑋,曹志剛,魏建強,等.聲源定位中的時延估計技術[J].數據采集與處理,2007,22(1):90-99.
[10]童峰,許肖梅,方世良,等.一種單頻水聲信號多徑時延估計算法[J].聲學學報,2008,33(1):62-68.
[11]行鴻彥,唐娟.時延估計方法的分析[J].聲學技術,2008,27(1):110-114.