張 凱 趙建虎 陳 卓

1)武漢大學測繪學院，武漢 430079

2)中兵勘察設計研究院，北京100053

1 引言

水下地形匹配導航作為慣性導航(INS)的一種重要輔助手段，以其無積累、自主性和隱蔽性強等特點備受關注［1－7］。目前，水下地形匹配的常用方法為基于貝葉斯方法的連續(xù)濾波估計［8－10］，然而，這些算法各有其局限性。點群濾波對計算處理單元的要求較高，尤其是當參數(shù)估計的維度從二維增加到三維時，其計算負擔重。以此相反，粒子濾波由于其原理特性，適用于高維度的貝葉斯連續(xù)估計，但為了保證計算精度和防止樣本貧化，也產(chǎn)生了計算速度和濾波性能之間的矛盾。與此相反，基于等值線的互相關地形匹配則不存在上述問題，且具有不受深度方向系統(tǒng)誤差影響的優(yōu)點。但是該方法容易受到誤匹配結果的影響，從而影響連續(xù)定位的精度和穩(wěn)定性。為此，本文給出了一種互相關水下地形匹配導航算法，以期提高匹配導航解的可靠性。

2 水下地形導航原理

水下地形導航系統(tǒng)的狀態(tài)方程和量測方程可分別表示為:

其中水下載體的位置狀態(tài)向量通常表示為

式中ut代表由在不同水下地形匹配觀測歷元中INS所提供的狀態(tài)改變量，rt是由INS 漂移率和歷元間航行距離所決定的傳播過程噪聲，服從正態(tài)分布N(0，Qt)，其中Qt為傳播噪聲水平。與此相對，地形匹配導航的測量更新過程則由式(2)描述。在式(2)中，yt表示由聲吶系統(tǒng)在當前實際位置獲取的水下地形測量值，h(xt)則代表當前位置估值xt在地形模型h(·)中所對應的水深值，et表示yt的測量噪聲，服從正態(tài)分布N(0，σe)2。在導航過程中，et由聲吶系統(tǒng)中的測量噪聲和地形模型噪聲兩部分組成。

基于聲吶系統(tǒng)測量噪聲和地形模型噪聲相互獨立，以及聲吶系統(tǒng)各波束之間測量噪聲相互獨立的前提，載體基于當前位置所獲取的地形測量值的概率密度函數(shù)為:

式中N 代表聲吶系統(tǒng)的測量波束數(shù)。

通常，由于式(2)為非線性方程，因此需要使用基于貝葉斯估計的方法如點群濾波和粒子濾波等方法來進行連續(xù)導航。然而，當互相關匹配中波束觀測值的數(shù)量N 很大時，其匹配結果呈現(xiàn)出高斯分布的形態(tài)，收斂于高斯分布N(0，R)，

此時，將量測方程(2)重新定義為:

式中，zt為互相關匹配結果，wt為匹配誤差，由式(7)給出?；?/p>

的Kalman 濾波方程為:

3 基于等值線互相關水下地形匹配

該導航方法基于水深測量值中的測量噪聲為高斯白噪聲。然而，受潮位變化等因素影響，水深測量值中的噪聲常常出現(xiàn)系統(tǒng)誤差，其幅度隨時間而變化，并影響水下地形匹配導航解的精度，還常常導致濾波發(fā)散。為此使用基于等值線的地形匹配方法來消除深度方向上的系統(tǒng)誤差，即在地形匹配的過程中，同時將測量數(shù)據(jù)和每批待匹配區(qū)域數(shù)據(jù)同時減去各自的均值，然后進行互相關匹配。基于等值線的地形匹配方法可有效消除測深數(shù)據(jù)中系統(tǒng)偏差的影響，然而由于算法求差過程中消去了絕對水深信息，匹配結果中易產(chǎn)生多峰性，同時誤匹配的概率也隨之增大，嚴重影響匹配的精度和穩(wěn)定性。為此，Nygren［11］提出了使用假設檢驗對互相關匹配結果進行質(zhì)量控制，對潛在的誤匹配結果進行探測。

在假設檢驗中，通過構造

來判斷誤匹配的概率。式中，Δhi為相似地形之間的地形差異，σe為地形觀測噪聲水平。當誤匹配的概率大于指定闕值t 時，將不使用本次匹配結果進行量測更新。

4 利用Unscented 變換進行不同歷元測量數(shù)據(jù)的融合

誤匹配產(chǎn)生的原因是由于載體真實位置所對應的地形區(qū)域z1在匹配搜索范圍內(nèi)有與之相似的其他地形區(qū)域z2，在測量噪聲的影響下，z1與真實位置在海底地形圖中所對應的地形區(qū)域z1’兩者之間的差異大于z2與z1’之間的差異。當搜索范圍和測量噪聲固定，同時參與匹配的地形測量值增加時，產(chǎn)生誤匹配的幾率將以指數(shù)級的速度下降［11］。因此，可以通過增加單次匹配中地形測量值數(shù)量，即同時使用相鄰觀測歷元獲取的地形測量值進行匹配的方法來減小誤匹配出現(xiàn)的幾率。然而，由于不同地形測量歷元之間的INS 數(shù)據(jù)存在誤差，相應的地形匹配結果也將受到INS 誤差的影響。因此，在誤匹配診斷判斷可能產(chǎn)生誤匹配的時候，使用Unscented 變換［12］(UT)進行不同地形測量歷元間的地形數(shù)據(jù)融合，以減少INS 誤差轉換到地形測量值誤差過程中海底地形非線性的影響，從而通過增加地形匹配數(shù)據(jù)量的方法提高地形匹配結果的精度和穩(wěn)定性。

通過使用UT，sigma 點選擇如下:

λ=α2(n+κ)－n，其中α 和κ 分別為控制sigma 點分布距離的參數(shù)。每個sigma 點Xj對應兩個權值和，分別用于計算變換后的均值和方差。

通過使用sigma 點，新的地形測量值ynew，j～N(uy，i，Qy，i)為:

式中，i表示單次地形測量值的對應波束序數(shù)。

5 導航方法

本文使用的導航方法采用基于等值線的互相關方法進行地形匹配。匹配結束后，使用假設檢驗對匹配結果進行誤匹配診斷。對于通過檢驗的匹配結果，將其作為Kalman 濾波的量測更新觀測值進行濾波。而對于未能通過假設檢驗的匹配結果，算法將不使用其進行測量更新，并將其與下一次觀測獲得的地形測量數(shù)據(jù)進行融合，組成新的地形測量數(shù)據(jù)，然后使該數(shù)據(jù)重新進行地形匹配。只有當匹配結果通過假設檢驗時，算法才會使用其進行量測更新。

6 實驗結果及分析

為了驗證文中所述方法的有效性，以渤海灣3 km×7 km 水域多波束測深實測數(shù)據(jù)為背景場(圖1)進行仿真實驗。水下載體航在航行過程中，載體每隔110 m 使用多波束測深系統(tǒng)進行海底地形測量，多波束系統(tǒng)的波束數(shù)為121。仿真中，背景場地形圖噪聲誤差設為0.1 m。與此同時，INS 中的漂移誤差設為航行距離的1%，載體在水平方向上的初始搜索半徑為500 m。

實驗中，先將聲吶測量噪聲水平設為0.1 m，進行蒙特卡洛仿真實驗(圖2)。如圖2，在地形測量值的精度較好的情況下，基于等值線的互相關匹配的結果具有很高的精度，且不易受誤匹配的干擾。

圖1 實驗區(qū)域水下地形背景場及實驗航跡Fig.1 Background of underwater terrain and the test track in the test area

圖2 導航算法的定位偏差Fig.2 Positioning deviation of navigation algorithm

圖3 地形觀測噪聲水平為0.2 m 時的診斷結果Fig.3 Diagnosed results of terrain observation with noise level 0.2m

隨后將聲吶測量噪聲水平增大為0.2 m，并進行仿真實驗(圖3、圖4)。圖3 為匹配導航過程中進行誤匹配診斷的結果，從圖3 可以看出，隨著地形測量值中的噪聲增大，產(chǎn)生誤匹配的幾率也相應增大，同時誤匹配診斷的警報次數(shù)也隨之增加。而在誤匹配診斷發(fā)出警報之后，使用Unscented 變換融合相鄰歷元的地形觀測數(shù)據(jù)重新進行匹配，相應的誤匹配診斷結果見圖4，而未進行誤匹配診斷的導航定位結果由圖4(b)給出。從圖4 可以看出，在易于產(chǎn)生誤匹配的區(qū)域使用融合數(shù)據(jù)進行匹配，顯著地減小了誤匹配的幾率。在使用融合數(shù)據(jù)進行匹配之后，誤匹配診斷僅發(fā)出一次警報，而圖3 中的誤匹配結果也通過使用融合數(shù)據(jù)被消除。

圖4 導航算法的定位偏差及誤匹配的影響Fig.4 Influence of positioning deviation and false matching of navigation algorithm

7 結論

本文提出了一種基于假設檢驗和Unscented 變換數(shù)據(jù)融合的地形匹配導航方法，該方法通過使用假設檢驗進行誤匹配診斷，并在誤匹配被探測時，再使用Unscented 變換融合相鄰歷元的地形觀測數(shù)據(jù)重新進行匹配，以提高匹配定位的精度和穩(wěn)定性。實驗驗證了該方法的有效性。

1 劉勇.海底地形匹配導航技術研究［D］.哈爾濱工程大學，2009.(Liu Yong.Study on navigation technique with underwater terrain matching technology［D］.Harbin Engineering University.2009)

2 譚佳琳.基于粒子群優(yōu)化的地形匹配導航算法研究［D］.哈爾濱工程大學，2009.(Tan Jialin.Study on terrain matching navigation algorithm based on particle swarm optimization［D］.Harbin Engineering University.2009)

3 李恒.基于水下地形匹配的捷聯(lián)系統(tǒng)誤差估計方法［J］.海軍工程大學學報，2012，24(1):40－44.(Li Heng.SINS error estimation based on underwater terrain matching［J］.Journal of Naval University of Engineering.2012，24(1):40－44)

4 張濤.混沌優(yōu)化水下地形匹配算法研究［J］.中國慣性技術學報，2009，17(2):156－164.(Zhang Tao.Study on underwater terrain matching algorithm based on chaotic optimization［J］.Chinese Journal of inertial technology.2009，17(2):156－164)

5 龐永杰.基于改進貝葉斯估計的水下地形匹配輔助導航方法［J］.儀器與儀表學報，2012，33(10):2 161－2 167.(Pang Yongjie.Underwater terrain matching navigation TMNA method based on improved Bayesian estimation［J］.Chinese Journal of Scientific Instrument.2012，33(10):2 161－2 167)

6 諶劍.一種基于粒子濾波的水下地形匹配算法研究［J］.海軍工程大學學報，2009，12(6):107－112.(Shen Jian.Study on underwater terrain matching algorithm based on particle filter［J］.Journal of Naval University of Engineering.2009，12(6):107－112)

7 張凱，趙建虎.基于支持向量機的水下地形適配/誤配區(qū)劃分方法的研究［J］.大地測量與地球動力學，2013(待刊).(Zhang Kai and Zhao Jianhu.Study on recognition and classification of matching/unmatching underwater background area based on SVM［J］.Journal of Geodesy and Geodynamics，2013(to be published))

8 Anonsen K B and Hallingstad O.Terrain aided underwater navigation using point mass and particle filters［A］.Proceedings of the IEEE/ION Position Location and Navigation Symposium［C］.2006.

9 Meduna D K，Rock S M and McEwen R.Low-cost terrain relative navigation for long-range AUVs［R］.In OCEANS 2008.Oceans，Polesand Climate:Technological Challenges，QuebecCity，Quebec，September 2008.

10 Karlsson R，Gusfafsson F and Karlsson T.Particle filtering and Cramer-Rao lower bound for underwater navigation［R］.Proceedings of IEEE Int Conf Acoust Speech.，Signal Process.HongKong，Apr.2003.

11 Nygren I.Terrain navigation for underwater vehicles［D］，School of Electrical Engineering，Royal Institute of Technology(KTH)，Sweden，Stockholm，Sweden.2005.

12 Julier S J and Uhlmann J K.Unscented filtering and nonlinear estimation［A］.In Proceedings of the IEEE［C］.2004，92:401－422.doi:10.1109/JPROC.2003.823141.