周家新， 陳建勇， 單志超， 陳長(zhǎng)康

(1.海軍海洋測(cè)繪研究所， 天津 300061； 2.海軍航空工程學(xué)院 電子信息工程系, 山東 煙臺(tái) 264001)

0 引言

鐵磁性物體受到地磁場(chǎng)磁化產(chǎn)生的磁場(chǎng)，與環(huán)境原有的背景磁場(chǎng)疊加在一起，使得鐵磁性目標(biāo)所在區(qū)域的空間磁場(chǎng)分布發(fā)生改變，從而形成磁異常。通過(guò)識(shí)別這種磁異常信號(hào)進(jìn)行目標(biāo)檢測(cè)與定位的技術(shù)稱(chēng)為磁異常檢測(cè)[1-2]。由于航空磁探儀受空氣、海水、泥沙土壤等介質(zhì)的影響小，且具有可以連續(xù)搜索、使用簡(jiǎn)單可靠、定位精度高、反應(yīng)迅速、搜索效率高等優(yōu)點(diǎn)[3]，作為被動(dòng)檢測(cè)技術(shù)的一種，已成為水下潛艇目標(biāo)的重要探測(cè)手段，被廣泛應(yīng)用于地磁勘探、管線(xiàn)搜索以及未爆彈藥探測(cè)等[4-6]領(lǐng)域。

目前，磁異常檢測(cè)技術(shù)主要分為兩類(lèi)：一是基于目標(biāo)特征，將鐵磁性目標(biāo)視為磁偶極子模型，從而給出目標(biāo)磁異常的數(shù)學(xué)表達(dá)式，進(jìn)而使用統(tǒng)計(jì)分析理論進(jìn)行目標(biāo)信號(hào)檢測(cè)，主要有正交基函數(shù)(OBF)磁異常檢測(cè)器與主成分分析(PCA)檢測(cè)器[7-10]；二是基于磁背景噪聲的統(tǒng)計(jì)特性，通過(guò)噪聲與目標(biāo)信號(hào)之間的差異，借助統(tǒng)計(jì)學(xué)中的假設(shè)檢驗(yàn)進(jìn)行判決，主要有最小熵檢測(cè)器(MED)與高階過(guò)零(HOC)檢測(cè)器[11-14]。

傳統(tǒng)基于目標(biāo)特征的OBF磁異常檢測(cè)器[7]和PCA檢測(cè)器[10]通常基于Anderson函數(shù)將目標(biāo)信號(hào)進(jìn)行正交基分解，得到能量信號(hào)進(jìn)行閾值檢測(cè)。該方法需要對(duì)檢測(cè)平臺(tái)磁補(bǔ)償后剩余磁噪聲的特性進(jìn)行分析和處理，使其滿(mǎn)足能量正交分解最優(yōu)檢測(cè)的相關(guān)要求?；谠肼曁匦缘腗ED磁異常檢測(cè)通常將磁背景噪聲視為高斯分布，但在復(fù)雜磁環(huán)境中，磁噪聲分布具有非高斯性，主要的磁噪聲源為地磁噪聲[13]，其功率譜密度為1/fα，其中f為地磁噪聲頻率，指數(shù)α的范圍為0～2. OBF檢測(cè)算法被證明是高斯白噪聲情況下的一個(gè)最優(yōu)檢測(cè)方法，但該方法依賴(lài)于背景噪聲特性。而在實(shí)際磁場(chǎng)測(cè)量環(huán)境中，噪聲功率譜密度為1/fα更常見(jiàn)。另外，MED檢測(cè)器[11]與HOC檢測(cè)器[12]需要大量的先驗(yàn)噪聲樣本，在低信噪比的實(shí)際應(yīng)用中往往難以滿(mǎn)足檢測(cè)需求。

文獻(xiàn)[15]提出了一種適用于非高斯噪聲背景的帶通濾波結(jié)合OBF分解的磁異常信號(hào)檢測(cè)算法，將非高斯噪聲進(jìn)行白化處理，基于OBF分解得到能量信號(hào)并用于檢測(cè)鐵磁性目標(biāo)。但由于航空磁探潛的背景噪聲干擾強(qiáng)，該檢測(cè)算法的性能有限。文獻(xiàn)[16]提出了一種航空磁探中水下目標(biāo)的自適應(yīng)探測(cè)方法，基于學(xué)習(xí)飛行補(bǔ)償背景磁場(chǎng)，并對(duì)噪聲進(jìn)行自適應(yīng)白化處理，從而對(duì)目標(biāo)信號(hào)進(jìn)行檢測(cè)。但該方法的學(xué)習(xí)飛行不穩(wěn)定，對(duì)背景磁場(chǎng)補(bǔ)償效果有限，難以將目標(biāo)信號(hào)從強(qiáng)磁背景干擾中準(zhǔn)確識(shí)別。

為提高航空磁探儀在低信噪比條件下對(duì)水下鐵磁性潛艇目標(biāo)的檢測(cè)概率，降低航空磁探潛過(guò)程中的虛警概率，本文提出了一種聯(lián)合估計(jì)檢測(cè)方法?；诖排紭O子建立潛艇目標(biāo)磁場(chǎng)信號(hào)模型，使用航空磁力儀觀測(cè)樣本，通過(guò)遺傳算法(GA)在潛艇目標(biāo)先驗(yàn)信息范圍內(nèi)搜索目標(biāo)的磁矩矢量、位置矢量和航向等參數(shù)，初步檢測(cè)目標(biāo)是否存在并得到估計(jì)信號(hào)。使用目標(biāo)的估計(jì)信號(hào)與測(cè)量平臺(tái)的實(shí)測(cè)信號(hào)進(jìn)行平均相關(guān)檢測(cè)，從而更進(jìn)一步確認(rèn)潛艇目標(biāo)的存在，得到目標(biāo)的位置、磁矩等參數(shù)，同時(shí)實(shí)現(xiàn)目標(biāo)的檢測(cè)、定位和分類(lèi)。

1 潛艇目標(biāo)磁場(chǎng)信號(hào)模型

航空磁異常探潛中，航空磁探儀和水下潛艇目標(biāo)的距離一般都在數(shù)百米之遙，甚至更遠(yuǎn)，根據(jù)文獻(xiàn)[17]的研究結(jié)論，在這種情況下，目標(biāo)磁場(chǎng)可以等效為偶極子磁場(chǎng)[18-19]。

對(duì)于一般的鐵磁性潛艇目標(biāo)，其在空間點(diǎn)P(x,y,z)處產(chǎn)生的磁場(chǎng)強(qiáng)度H可以根據(jù)磁偶極子模型計(jì)算得到，以水下潛艇目標(biāo)中心為坐標(biāo)原點(diǎn)(見(jiàn)圖1)，則可以得到P點(diǎn)處的磁場(chǎng)強(qiáng)度為

(1)

式中：M為潛艇目標(biāo)的磁矩矢量；r0為目標(biāo)源點(diǎn)O指向空間測(cè)量點(diǎn)P的單位矢量；r為目標(biāo)源點(diǎn)O與空間測(cè)量點(diǎn)P之間的距離；o(1/r5)為高階次極子在鐵磁性潛艇目標(biāo)磁場(chǎng)中的貢獻(xiàn)。隨著測(cè)量點(diǎn)和目標(biāo)之間距離的增大，高階次極子的作用逐漸變小。在忽略高階次極子的作用之后，水下鐵磁性潛艇目標(biāo)磁場(chǎng)信號(hào)模型可表示為

(2)

圖1中潛艇目標(biāo)磁矩M可以表示為3個(gè)分量的形式，分別為Mx、My、Mz，將磁偶極子模型轉(zhuǎn)換到直角坐標(biāo)系下，磁偶極子模型在空間產(chǎn)生的磁場(chǎng)矢量可以表示為H=ixHx+iyHy+izHz，ix、iy、iz為3個(gè)坐標(biāo)軸的單位矢量，則三分量Hx、Hy、Hz為

(3)

在航空磁異常探潛中，將磁場(chǎng)強(qiáng)度H在真空(空氣)中引起的磁感應(yīng)強(qiáng)度記為B，并有如下關(guān)系：

(4)

式中：μ0為真空磁導(dǎo)率，μ0=4π×10-7H/m. 在海水中磁導(dǎo)率近似等于真空磁導(dǎo)率，而且由于磁感應(yīng)強(qiáng)度的基本單位T非常大，航空磁探測(cè)中多使用nT作為單位。

將(3)式代入(4)式，得到水下潛艇目標(biāo)在空間產(chǎn)生的磁感應(yīng)強(qiáng)度三分量Bx、By、Bz為

(5)

使用(5)式計(jì)算時(shí)，右邊仍采用國(guó)際單位制，左邊所得結(jié)果單位為nT.

將潛艇目標(biāo)空間三分量磁場(chǎng)轉(zhuǎn)換為航空光泵磁探儀測(cè)量得到的空間磁場(chǎng)總場(chǎng)信號(hào)，設(shè)θ為磁北方向到潛艇目標(biāo)航向的逆時(shí)針角度(航向角)，將目標(biāo)的三分量磁感應(yīng)強(qiáng)度投影到地磁場(chǎng)方向，則目標(biāo)磁感應(yīng)強(qiáng)度的水平分量Bh和垂直分量Bv的形式為

(6)

沿地磁場(chǎng)方向的單位磁場(chǎng)向量為(cosγ,sinγ)，其中γ為測(cè)量點(diǎn)處的地磁傾角，則鐵磁性潛艇目標(biāo)的光泵測(cè)量值即為目標(biāo)磁場(chǎng)向量沿地磁場(chǎng)方向的投影：

B=Bhcosγ+Bvsinγ.

(7)

將(5)式、(6)式代入(7)式中，整理得到鐵磁性潛艇目標(biāo)光泵測(cè)量信號(hào)模型為

(8)

目標(biāo)磁感應(yīng)強(qiáng)度三分量Bx、By、Bz及其投影到地磁場(chǎng)方向的光泵總場(chǎng)信號(hào)如圖2所示，其中該信號(hào)以磁偶極子為中心，測(cè)線(xiàn)中心位于目標(biāo)中心正上方，高度為200 m，測(cè)線(xiàn)的范圍x∈[-450 m,450 m]，磁矩為Mx=80 000 A·m2，My=-20 000 A·m2，Mz=40 000 A·m2，目標(biāo)航向角為0°，地磁傾角為45°.

由(8)式可知，鐵磁性目標(biāo)信號(hào)有8個(gè)參數(shù)，若其中地磁傾角γ可以認(rèn)為已知，則仍有位置x、y、z，磁矩Mx、My、Mz，航向角θ等7個(gè)未知參數(shù)。

航空磁探儀通過(guò)在任務(wù)區(qū)域的飛行來(lái)實(shí)施航空磁異常探潛任務(wù)，建立如圖3所示的探測(cè)坐標(biāo)系。

經(jīng)典探測(cè)寬度模型[20]可以表示為

(9)

式中：Rm表示航空磁探儀的作用距離；Hp表示磁探儀所在高度；Hs表示目標(biāo)所在深度。設(shè)航空磁探儀直線(xiàn)飛行，從磁探儀采樣信號(hào)開(kāi)始到結(jié)束的飛行距離為L(zhǎng)，為簡(jiǎn)化探測(cè)過(guò)程，將探測(cè)區(qū)域近似為一個(gè)長(zhǎng)方體，從而可以得到目標(biāo)所在區(qū)域長(zhǎng)為L(zhǎng)+W，單位均為m. 基于磁探儀的探測(cè)寬度模型和磁偶極子信號(hào)模型，能夠?qū)崿F(xiàn)磁性目標(biāo)的遠(yuǎn)程定位，有利于快速縮小目標(biāo)的搜索范圍。

2 基于GA和Fréchet距離的搜索估計(jì)量

航空磁探潛過(guò)程中，磁力儀采樣一段信號(hào)，信號(hào)長(zhǎng)度為N，觀測(cè)信號(hào)分別記為B1，B2，…，Bn，…，BN，即

(10)

由(8)式的信號(hào)模型可知，待估計(jì)參量矢量為

a=[x,y,z,Mx,My,Mz,θ]T.

(11)

根據(jù)文獻(xiàn)[21]設(shè)置探測(cè)區(qū)域以及待估計(jì)參量矢量的范圍。基于目標(biāo)先驗(yàn)信息，在探測(cè)區(qū)域內(nèi)，使用GA對(duì)待估計(jì)參量進(jìn)行搜索[22-24]，基于Fréchet距離構(gòu)造適應(yīng)度函數(shù)?；贕A的參量估計(jì)使得適應(yīng)度函數(shù)在搜索過(guò)程中總能收斂到一個(gè)極小值點(diǎn)，能以很大概率找到全局最優(yōu)解。

Fréchet距離應(yīng)用廣泛[25]，可描述兩曲線(xiàn)之間的相似性程度[26]。在航空磁探潛水下目標(biāo)的聯(lián)合估計(jì)檢測(cè)中，可以用于量化估計(jì)信號(hào)與實(shí)測(cè)信號(hào)的相似性。鐵磁性目標(biāo)的估計(jì)信號(hào)為B()，實(shí)測(cè)信號(hào)B(a)與B()為兩個(gè)離散的點(diǎn)序列串，其中B(a)為〈B1(a),B2(a),…，BN(a)〉，估計(jì)信號(hào)B()為〈B1()，B2()，…，BM()〉，M為估計(jì)信號(hào)的長(zhǎng)度。則Fréchet距離為

Fd(B(a),B())=max (Ed(B(a),B()),
min (Fd(〈B1(a),B2(a),…,BN-1(a)〉,
〈B1(),B2(),…,BM()〉),?N≠1,
Fd(〈B1(a),B2(a),…,BN(a)〉,
〈B1(),B2(),…,BM-1()〉),?M≠1,
Fd(〈B1(a),B2(a),…,BN-1(a)〉,〈B1(),
B2(),…,BM-1()〉),?M≠1,N≠1))，

(12)

式中：Ed(B(a),B())為B(a)與B()之間的歐氏距離。根據(jù)文獻(xiàn)[21]中的數(shù)據(jù)設(shè)定一組目標(biāo)位置、磁矩及航向等參數(shù)，生成目標(biāo)空間磁場(chǎng)分布，模擬磁力儀采樣目標(biāo)信號(hào)，得到仿真的實(shí)測(cè)信號(hào)B(a). 然后基于GA在區(qū)域范圍內(nèi)搜索待估計(jì)參量，搜索過(guò)程如圖4所示，在規(guī)定的遺傳代數(shù)(3 000)內(nèi)搜索，若適應(yīng)度未達(dá)到門(mén)限值則判定目標(biāo)不存在。針對(duì)搜索過(guò)程，在空間直角坐標(biāo)系中標(biāo)記出估計(jì)參數(shù)和設(shè)定參數(shù)位置。圖5為基于GA的潛艇目標(biāo)在航空磁探儀搜索區(qū)域范圍內(nèi)的待估計(jì)參數(shù)搜索定位圖，目標(biāo)位置參數(shù)為x=500 m，y=100 m，z=150 m，磁矩參數(shù)為Mx=80 000 A·m2，My=-20 000 A·m2，Mz=40 000 A·m2，航向角參數(shù)為θ=190°. 在圖5(a)和圖5(b)中，藍(lán)色方塊為目標(biāo)的真實(shí)位置以及真實(shí)磁矩，紅色十字點(diǎn)為每次搜索定位的結(jié)果；圖5(c)中，藍(lán)色直線(xiàn)為θ=190°. 由搜索定位結(jié)果可以得到，基于GA的搜索目標(biāo)參量，能夠快速收斂到實(shí)際真值。

3 使用目標(biāo)估計(jì)信號(hào)的平均相關(guān)檢測(cè)

航空磁探潛中的磁異常檢測(cè)屬于雙擇檢測(cè)，水下潛艇目標(biāo)存在狀態(tài)對(duì)應(yīng)假設(shè)H0和假設(shè)H1兩種情況。根據(jù)磁探儀采樣的觀測(cè)信號(hào)，水下潛艇目標(biāo)雙擇檢測(cè)問(wèn)題表示如下：

(13)

式中：H0和H1對(duì)應(yīng)兩種假設(shè)情況，H0表示目標(biāo)不存在，H1表示目標(biāo)存在；x(n)為觀測(cè)信號(hào)；ω(n)為噪聲信號(hào)；n=1，2，…，N.

通過(guò)估計(jì)目標(biāo)的磁矩矢量、位置矢量以及航向等參數(shù)，得到目標(biāo)磁感應(yīng)強(qiáng)度三分量的估計(jì)，進(jìn)一步得到目標(biāo)的估計(jì)信號(hào)。基于目標(biāo)的估計(jì)信號(hào)與測(cè)量平臺(tái)的觀測(cè)信號(hào)進(jìn)行平均相關(guān)檢測(cè)，從而判定目標(biāo)是否存在。構(gòu)造如下檢驗(yàn)統(tǒng)計(jì)量：

(14)

式中：輸入x為觀測(cè)信號(hào)矢量；RxB(m)為觀測(cè)信號(hào)與估計(jì)信號(hào)的互相關(guān)函數(shù)。根據(jù)檢驗(yàn)統(tǒng)計(jì)量構(gòu)造的聯(lián)合估計(jì)檢測(cè)器實(shí)驗(yàn)框圖如圖6所示。根據(jù)(8)式的目標(biāo)光泵測(cè)量信號(hào)模型生成目標(biāo)的估計(jì)信號(hào)，通過(guò)相關(guān)器對(duì)估計(jì)信號(hào)與觀測(cè)信號(hào)做平均相關(guān)，將得到的檢驗(yàn)統(tǒng)計(jì)量與門(mén)限β進(jìn)行比較，判定目標(biāo)是否存在。檢測(cè)門(mén)限β根據(jù)奈曼- 皮爾遜準(zhǔn)則得到，由大量的觀測(cè)噪聲樣本進(jìn)行統(tǒng)計(jì)分析，計(jì)算得到最終的門(mén)限值。

4 檢測(cè)器性能分析

采用光泵磁力儀實(shí)測(cè)的磁噪聲數(shù)據(jù)和磁偶極子仿真目標(biāo)信號(hào)數(shù)據(jù)進(jìn)行聯(lián)合估計(jì)檢測(cè)器的性能分析，實(shí)測(cè)噪聲數(shù)據(jù)片段如圖7(a)所示。航空磁探儀飛行高度為150 m，飛行速度為100 m/s，飛行航向?yàn)?°，測(cè)量地區(qū)的背景磁場(chǎng)為49 699 nT，地磁傾角為52°，其中磁力儀的分辨率為0.01 nT，采樣頻率為100 Hz. 潛艇目標(biāo)信號(hào)使用(8)式的磁偶極子模型仿真產(chǎn)生，目標(biāo)信號(hào)加實(shí)測(cè)磁噪聲經(jīng)過(guò)去趨勢(shì)處理后如圖7(b)所示，目標(biāo)出現(xiàn)在采樣點(diǎn)2 500～2 800段。目標(biāo)信號(hào)寬度為300，檢測(cè)器輸出中目標(biāo)信號(hào)的寬度接近300. 為與文獻(xiàn)[7]中的OBF檢測(cè)器性能進(jìn)行比較，將目標(biāo)信號(hào)與實(shí)測(cè)磁噪聲之間的幅度信噪比設(shè)置為0.3. 基于目標(biāo)估計(jì)信號(hào)的聯(lián)合估計(jì)檢測(cè)器(ES-COD)歸一化輸出如圖7(c)所示，其輸出的幅度信噪比為10，能夠大大提高檢測(cè)器的輸出信噪比。文獻(xiàn)[7]中OBF檢測(cè)器的輸出幅度信噪比為5，相比之下，ES-COD檢測(cè)器的輸出幅度信噪比提升了1倍，檢測(cè)性能有所提高。文獻(xiàn)[15]提出的帶通濾波結(jié)合OBF檢測(cè)器分解的磁異常信號(hào)檢測(cè)算法適用于非高斯噪聲環(huán)境，但其輸出幅度信噪比僅為2. 在壓制噪聲能力方面，與該算法相比，本文提出的ES-COD檢測(cè)器性能更優(yōu)，能夠顯著提高信噪比。

選取實(shí)際鐵磁性潛艇模型代替潛艇目標(biāo)對(duì)檢測(cè)器性能進(jìn)行Monte Carlo仿真分析，艇模三視圖如圖8所示，噪聲信號(hào)采用實(shí)測(cè)噪聲數(shù)據(jù)。實(shí)驗(yàn)中，在潛艇模型中心上方10 m的高度面上沿x軸方向每0.1 m取1個(gè)測(cè)量點(diǎn)，共計(jì)151個(gè)考核點(diǎn)，x∈[-7.5 m, 7.5 m]。在考核點(diǎn)處使用G858光泵磁力儀進(jìn)行測(cè)量，磁力儀性能參數(shù)如表1所示。橫向測(cè)量范圍為15 m，潛艇模型長(zhǎng)度為2 m，測(cè)量范圍達(dá)到15倍艇長(zhǎng)。測(cè)量高度為10 m，是艇長(zhǎng)的5倍。測(cè)量結(jié)果如圖9所示?；谀温? 皮爾遜準(zhǔn)則，在虛警概率為0.001的條件下，各信噪比情況仿真10 000次，得到如圖10所示的檢測(cè)概率曲線(xiàn)。

銫光泵參數(shù)銫光泵數(shù)值測(cè)量范圍/nT18000～95000分辨率/nT0.05溫度漂移/(nT·℃-1)<0.05梯度容限/(nT·in-1) 500測(cè)量精度/%<1頻域噪聲/(pT·Hz-12)<4

分別將文獻(xiàn)[3]、文獻(xiàn)[7]和文獻(xiàn)[12]中的最小均方估計(jì)- 正交基函數(shù)(LMS-OBF)檢測(cè)器、OBF檢測(cè)器和HOC檢測(cè)器性能曲線(xiàn)繪于圖10中，與ES-COD檢測(cè)器性能進(jìn)行對(duì)比分析。

由圖10可以看出：在4種檢測(cè)器中，ES-COD檢測(cè)器的檢測(cè)概率最大；HOC檢測(cè)器相較于OBF檢測(cè)器性能有所提高；文獻(xiàn)[3]提出的LMS-OBF檢測(cè)器能夠克服海浪磁噪聲對(duì)航空磁異常檢測(cè)器的性能影響，相對(duì)于HOC檢測(cè)器，其性能大大提高；在低信噪比條件下，ES-COD檢測(cè)性能較OBF檢測(cè)器、HOC檢測(cè)器和LMS-OBF檢測(cè)器更優(yōu)。當(dāng)信噪比≥-5 dB時(shí)，檢測(cè)概率大于90%，能夠準(zhǔn)確發(fā)現(xiàn)目標(biāo)，可以滿(mǎn)足實(shí)際需求；當(dāng)信噪比≥-7 dB時(shí)，檢測(cè)概率仍保持在50%，相對(duì)于文獻(xiàn)[3]中的LMS-OBF檢測(cè)器，ES-COD檢測(cè)器的弱信號(hào)檢測(cè)能力更突出?；趯?shí)測(cè)數(shù)據(jù)分析可知，地磁異常信號(hào)與潛艇目標(biāo)信號(hào)相比信號(hào)寬度通常較大，飛行器平臺(tái)干擾信號(hào)寬度較小，二者與潛艇目標(biāo)磁場(chǎng)信號(hào)通常呈弱相關(guān)性。綜上所述，使用聯(lián)合估計(jì)檢測(cè)方法，在低信噪比條件下能夠降低地磁異常和飛行器平臺(tái)干擾信號(hào)引起的虛警概率。

5 結(jié)論

本文針對(duì)航空磁探測(cè)中潛艇目標(biāo)信號(hào)被強(qiáng)磁背景噪聲和干擾淹沒(méi)導(dǎo)致檢測(cè)概率低、虛警概率過(guò)高的問(wèn)題，提出了一種航空磁探中的聯(lián)合估計(jì)檢測(cè)方法，并給出了具體的實(shí)現(xiàn)方案。該方法的主要思想是通過(guò)磁偶極子建立潛艇目標(biāo)磁場(chǎng)信號(hào)模型，根據(jù)觀測(cè)信號(hào)，基于目標(biāo)先驗(yàn)信息，使用GA和Fréchet距離對(duì)目標(biāo)參量進(jìn)行估計(jì)?；谀繕?biāo)估計(jì)信號(hào)構(gòu)造了平均相關(guān)檢驗(yàn)統(tǒng)計(jì)量，設(shè)計(jì)了使用潛艇目標(biāo)估計(jì)信號(hào)的ES-COD聯(lián)合估計(jì)檢測(cè)器。得到如下結(jié)論：

1)使用GA搜索目標(biāo)參量，能夠快速收斂到真實(shí)值，得到目標(biāo)的估計(jì)參量。

2)由估計(jì)參量得到的潛艇目標(biāo)磁場(chǎng)估計(jì)信號(hào)與實(shí)測(cè)信號(hào)相近。

3)實(shí)驗(yàn)中ES-COD檢測(cè)器輸出的幅度信噪比約為10，該檢測(cè)器對(duì)航空磁探中的背景磁噪聲具有強(qiáng)壓制能力，能夠顯著提高信噪比。

4)當(dāng)信噪比≥-5 dB時(shí)，ES-COD檢測(cè)器檢測(cè)概率大于90%，表明該檢測(cè)器在低信噪比條件下仍具有良好的檢測(cè)能力。

5)基于奈曼- 皮爾遜準(zhǔn)則，在相同虛警概率條件下，ES-COD檢測(cè)器的檢測(cè)性能優(yōu)于OBF檢測(cè)器、HOC檢測(cè)器和LMS-OBF檢測(cè)器，并且可以同時(shí)實(shí)現(xiàn)目標(biāo)的檢測(cè)、定位以及分類(lèi)，具有一定的實(shí)際應(yīng)用價(jià)值。

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

[1] Liu Y, Zhang Y, Yi H. The new magnetic survey method for underwater pipeline detection [J]. Applied Mechanics and Materials, 2013, 239(2): 338-343.

[2] Dames P M, Schwager M, Schwager D. Active magnetic anomaly detection using multiple micro aerial vehicles [J]. IEEE Robotics and Automation Letters, 2016, 1(1): 153-160.

[3] 熊雄, 楊日杰, 王鴻吉. 海浪磁噪聲背景中動(dòng)目標(biāo)航空磁異常檢測(cè)算法[J].華中科技大學(xué)學(xué)報(bào)：自然科學(xué)版, 2015, 43(5): 100-105.

XIONG Xiong, YANG Ri-jie, WANG Hong-ji. Airborne magnetic anomaly detection algorithm for moving target under ocean wave generated magnetic noise [J]. Journal of Huazhong University of Science and Technology: Natural Science Edition, 2015, 43(5): 100-105. (in Chinese)

[4] Song L, Billings S, Pasion L, et al. Transient electromagnetic scattering of a metallic object buried in underwater sediments [J]. IEEE Transactions on Geoscience and Remote Sensing, 2016, 54(2): 1091-1102.

[5] Liu Z, Pang H, Pan M. Calibration and compensation of geomagnetic vector measurement system and improvement of magnetic anomaly detection [J]. IEEE Transactions on Geoscience and Remote Sensing Letters, 2016, 13(3): 447-451.

[6] Robert F, Saiful H, Mojtaba A, et al. Magnetic signature attenuation of an unmanned aircraft system for aeromagnetic survey [J]. IEEE/ASME Transactions on Mechatronics, 2014, 19(4): 1436-1446.

[7] Ginzburg B, Frumkis L, Kaplan B Z. Processing of magnetic scalar gradiometer signals using orthonormalized functions [J]. Sensors and Actuators, 2002, 102(1): 67-75.

[8] Ginzburg B, Frumkisb L, Kaplan B Z. An efficient method for processing scalar magnetic gradiometer signals [J]. Sensors and Actuators, 2004, 114(1): 73-79.

[9] Frumkisa L, Ginzburg B, Salomonski N. Optimization of scalar magnetic gradiometer signal processing [J]. Sensors and Actuators, 2005, 121(1): 88-94.

[10] Sheinker A, Moldwin M B. Magnetic anomaly detection (MAD) of ferromagnetic pipelines using principal component analysis (PCA) [J]. Measurement Science and Technology, 2016, 27(4): 1-7.

[11] Sheinker A, Salomonskil N, Ginzburg B. Magnetic anomaly detection using entropy filter [J]. Measurement Science and Technology, 2008, 19(5): 1-5.

[12] Sheinker A, Ginzburg B, Salomonski N. Magnetic anomaly detection using high-order crossing method [J]. IEEE Transactions on Geoscience and Remote Sensing, 2012, 50(4): 1095-1103.

[13] Sheinker A, Shkalim A, Salomonski N. Processing of a scalar magnetometer signal contaminated by 1/fαnoise [J]. Sensors and Actuators, 2007, 138(1): 105-111.

[14] Sheinker A, Frumkis L, Ginzburg B. Magnetic anomaly detection using a three-axis magnetometer [J]. IEEE Transactions on Magnetics, 2009, 45(1): 160-167.

[15] 張堅(jiān), 林春生, 鄧鵬, 等. 非高斯背景噪聲下的微弱磁異常信號(hào)檢測(cè)算法[J]. 海軍工程大學(xué)學(xué)報(bào), 2011, 23(4): 22-26.

ZHANG Jian, LIN Chun-sheng，DENG Peng, et al. Detection of weak magnetic anomaly signal in non-Gaussian noise background [J]. Journal of Naval University of Engineering, 2011, 23(4): 22-26. (in Chinese)

[16] 鄧鵬, 張堅(jiān), 林春生. 航空磁探中水下目標(biāo)的自適應(yīng)探測(cè)方法[J]. 艦船科學(xué)技術(shù), 2012, 34(3): 76-79.

DENG Peng, ZHANG Jian, LIN Chun-sheng. Adaptive detection method of magnetic anomaly detection in aeromagnetic survey [J]. Ship Science and Technology, 2012, 34(3): 76-79. (in Chinese)

[17] 張朝陽(yáng), 肖昌漢, 高俊吉, 等. 磁性物體磁偶極子模型適用性的試驗(yàn)研究[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào), 2010, 18(5): 862-868.

ZHANG Zhao-yang, XIAO Chang-han, GAO Jun-ji, et al. Experiment research of magnetic dipole model applicability for a magnetic object [J]. Journal of Basic Science and Engineering, 2010, 18(5): 862-868. (in Chinese)

[18] 郭志馗, 陳超, 陶春輝, 等. 有限長(zhǎng)圓柱體磁異常場(chǎng)全空間正演方法[J]. 地球物理學(xué)報(bào), 2017, 60(4): 1557-1570.

GUO Zhi-kui, CHEN Chao, TAO Chun-hui, et al. Forward modeling of magnetic anomalies of finite-length cylinder in 3D space [J]. Chinese Journal of Geophysics, 2017, 60(4): 1557-1570. (in Chinese)

[19] 周家新, 陳建勇, 單志超, 等. 航空磁探中潛艇磁場(chǎng)建模方法分析[J]. 海軍航空工程學(xué)院學(xué)報(bào), 2017, 32(1): 143-148.

ZHOU Jia-xin, CHEN Jian-yong, SHAN Zhi-chao, et al. Analysis of submarine magnetic field modeling method for aeromagnetic detection [J]. Journal of Naval Aeronautical and Astronautical University, 2017, 32(1): 143-148. (in Chinese)

[20] 楊日杰, 熊雄, 郭新奇, 等. 基于潛艇磁偶極子模型的航空磁探潛探測(cè)寬度模型與仿真[J]. 兵工學(xué)報(bào), 2014, 35(9): 1458-1465.

YANG Ri-jie, XIONG Xiong, GUO Xin-qi, et al. Research on model and simulation of airborne magnetic anomaly detection sweep width based on magnetic dipole model [J]. Acta Armamentarii, 2014, 35(9): 1458-1465. (in Chinese)

[21] Sheinker A, Salomonski N, Ginzburg B, et al. Aeromagnetic search using genetic algorithm [C]∥Progress in Electromagnetics Research Symposium (PIERS 2005). Hangzhou: The Electromagnetics Academy, 2005: 492-495.

[22] 劉淵, 楊永輝, 張春瑞, 等. 一種基于遺傳算法的Fuzzing測(cè)試用例生成新方法[J]. 電子學(xué)報(bào), 2017, 45(3): 552-556.

LIU Yuan, YANG Yong-hui, ZHANG Chun-rui, et al. A novel method for fuzzing test cases generating based on genetic algorithm [J]. Acta Electronica Sinica, 2017, 45(3): 552-556. (in Chinese)

[23] 孫躍, 夏金鳳, 唐春森, 等. 采用元胞遺傳算法的無(wú)線(xiàn)電能傳輸網(wǎng)路徑尋優(yōu)[J]. 西安交通大學(xué)學(xué)報(bào), 2017, 51(4): 30-36.

SUN Yue, XIA Jin-feng, TANG Chun-sen, et al. Path optimization for wireless power transfer networks with cellular genetic algorithm [J]. Journal of Xi’an Jiaotong University, 2017, 51(4): 30-36. (in Chinese)

[24] 高翔, 嚴(yán)勝剛, 李斌. 一種混合算法下單個(gè)磁通門(mén)定位運(yùn)動(dòng)磁性目標(biāo)研究[J].大連理工大學(xué)學(xué)報(bào), 2016, 56(3): 292-298.

GAO Xiang, YAN Sheng-gang, LI Bin. Study of a hybrid algorithm for localization of mobile magnetic target by a single fluxgate [J]. Journal of Dalian University of Technology, 2016, 56(3): 292-298. (in Chinese)

[25] 高孝杰, 簡(jiǎn)季, 戴曉愛(ài), 等. 基于Fréchet距離的光譜曲線(xiàn)匹配應(yīng)用分析[J]. 武漢大學(xué)學(xué)報(bào)：信息科學(xué)版, 2016, 41(3): 408-414.

GAO Xiao-jie, JIAN Ji, DAI Xiao-ai, et al. Spectral curve matching application analysis based on Fréchet distance [J]. Geoma-tics and Information Science of Wuhan University, 2016, 41(3): 408-414. (in Chinese)

[26] 安曉亞, 劉平芝, 楊云, 等. 一種線(xiàn)狀要素兒何相似性度量方法及其應(yīng)用[J]. 武漢大學(xué)學(xué)報(bào)：信息科學(xué)版, 2015, 40(9): 1225-1229.

AN Xiao-ya, LIU Ping-zhi, YANG Yun, et al. A geometric similarity measurement method and applications to linear feature [J]. Geomatics and Information Science of Wuhan University, 2015, 40(9): 1225-1229. (in Chinese)