吳 芳，吳 銘，林義杰

（1.海軍航空大學(xué)，山東煙臺(tái)264001；2.92485部隊(duì)，遼寧大連116113）

潛艇聲隱身性能的提高，使得傳統(tǒng)的聲探潛手段面臨著諸多挑戰(zhàn)，而磁異探測(cè)設(shè)備因可以清晰地探明海洋深處磁場(chǎng)的微弱變化，準(zhǔn)確定位出磁體與探測(cè)設(shè)備之間的相對(duì)位置，及時(shí)有效地發(fā)現(xiàn)敵方潛艇，受到各國(guó)海軍的重視[1-4]。

本文基于潛艇磁場(chǎng)時(shí)域信號(hào)的信息冗余，對(duì)潛艇磁場(chǎng)信號(hào)進(jìn)行了頻域分析，得到更為直觀的特征數(shù)據(jù)，有助于降低目標(biāo)磁異常信號(hào)的誤判概率[5-14]。本文主要分析了我周邊海域的水下目標(biāo)磁異常分布、磁場(chǎng)信號(hào)的頻域特征，仿真研究了潛艇不同下潛深度、不同運(yùn)動(dòng)速度下的磁場(chǎng)頻域信號(hào)特點(diǎn)，為我國(guó)附近海域水下潛艇的發(fā)現(xiàn)與識(shí)別提供了一定的判別依據(jù)。

1 潛艇磁異常分布

對(duì)于潛艇目標(biāo)，在無電流區(qū)域，存在方程：

式（1）中，H 為潛艇磁場(chǎng)強(qiáng)度。

標(biāo)量磁勢(shì)Vm就可以根據(jù)以下關(guān)系式定義：

磁場(chǎng)強(qiáng)度H 與磁感應(yīng)強(qiáng)度B 存在以下關(guān)系：

式（3）中：μ0=4π×10-7H/m，為真空磁導(dǎo)率；μr為相對(duì)磁導(dǎo)率。

由式（1）～（3）可得標(biāo)量磁勢(shì)Vm的方程：

在這個(gè)模型中使用了約化場(chǎng)公式，所以只需要求解對(duì)應(yīng)于擾動(dòng)場(chǎng)的標(biāo)量磁勢(shì)Vm，其求解公式為：

式中，Be為背景地磁場(chǎng)。

對(duì)于潛艇艇殼，設(shè)置磁屏蔽，則式（5）為：

式（6）中：n 為法向量；B1為艇殼內(nèi)磁感應(yīng)強(qiáng)度；B2為艇殼外磁感應(yīng)強(qiáng)度；ds為艇殼的厚度。

基于上述模型，可仿真不同地磁背景條件下的潛艇磁異常分布，如圖1～4所示。

圖1 地磁場(chǎng)方向豎直向下Fig.1 The direction of the geomagnetic field is vertical downward

圖2 地磁場(chǎng)方向豎直向下時(shí)，潛艇磁場(chǎng)分布Fig.2 The submarine’s magnetic field distribution when the direction of the geomagnetic field is vertical downward

圖3 地磁場(chǎng)方向從南向北沿45°向下Fig.3 The direction of the geomagnetic field is 45°downward from south to north

圖4 地磁場(chǎng)方向從南向北沿45°向下時(shí)，潛艇磁場(chǎng)分布Fig.4 The submarine’s magnetic field distribution when the direction of the geomagnetic field is 45° downward from south to north

圖1～4中，潛艇在水下引起磁場(chǎng)發(fā)生了明顯變化，顏色條為總磁感應(yīng)強(qiáng)度值，白色箭頭顯示了船體中切向磁場(chǎng)的方向和強(qiáng)度，為潛艇在背景場(chǎng)中產(chǎn)生的擾動(dòng)。當(dāng)背景磁場(chǎng)方向不同時(shí)，潛艇引起的磁異常磁場(chǎng)也不同。當(dāng)背景磁場(chǎng)方向豎直向下時(shí)，潛艇引起的磁異常大致上是橢圓狀，從中心向四周逐漸減小一段距離后稍微增加回歸50 000 nT，并且是大致均勻且關(guān)于潛艇的艏艉軸線對(duì)稱的；而當(dāng)背景磁場(chǎng)方向從南向北沿45°向下時(shí)，潛艇艏艉附近引起的異常磁場(chǎng)出現(xiàn)了2個(gè)前后順序排列的峰值。

當(dāng)背景磁場(chǎng)方向從南向北沿45°向下時(shí)，潛艇縱軸方向正上方50 m、100 m 的磁感應(yīng)強(qiáng)度分布如圖5、6所示。

圖5 潛艇下潛50 m 處磁場(chǎng)分布Fig.5 Magnetic field distribution of submarine at 50 m

圖6 潛艇下潛100 m 處磁場(chǎng)分布Fig.6 Magnetic field distribution of submarine at 100 m

如圖5、6所示，紅色區(qū)域指示較高的磁場(chǎng)強(qiáng)度，藍(lán)色區(qū)域表示較低的磁場(chǎng)強(qiáng)度。背景磁場(chǎng)方向斜向下時(shí)，潛艇引起磁場(chǎng)變化，可以發(fā)現(xiàn)總磁場(chǎng)存在高低2個(gè)大的峰值，且峰值出現(xiàn)在潛艇艏艉附近，2個(gè)峰都均勻變化，并且是關(guān)于潛艇的艏艉軸線對(duì)稱的。

2 潛艇磁信號(hào)的頻域分析

2.1 平均周期圖（Welch）法功率譜估計(jì)原理

信號(hào)的功率譜反映了信號(hào)功率隨頻率的變化。平穩(wěn)信號(hào)的功率譜S(f) 可以通過對(duì)其自相關(guān)函數(shù)Rs(τ)進(jìn)行傅里葉變換得到[15]，即：

傅里葉變換周期圖方法對(duì)信號(hào)序列XN(n)執(zhí)行離散傅立葉變換。因?yàn)殡x散傅里葉變換的周期是N ，且在傅里葉變換后的值XN(ejω)的平方再除以序列長(zhǎng)度N，所以其功率譜S(ejω)也以N 為周期[16]，可表示為：

化簡(jiǎn)后得到：

即功率譜S(k)可以經(jīng)由XN(n)做快速傅里葉變換很便利地得到。但是，這樣得到的功率譜估計(jì)值頻譜特性不是很好，屬于有偏估計(jì)，適用于長(zhǎng)信號(hào)序列的情況。周期圖法得到的功率譜曲線，當(dāng)N 太大時(shí)，功率譜變化太大；當(dāng)N 太小時(shí)，曲線的峰值又不好辨識(shí)。為了對(duì)其進(jìn)行改進(jìn)，誕生了平均周期圖（Welch）法。

平均周期圖法[17]是把長(zhǎng)度為N 的數(shù)據(jù)XN(n)分成K 個(gè)小段，每小段L 個(gè)點(diǎn)，每相鄰小段之間重疊P 個(gè)點(diǎn)（P ＜L），一般采用P=L/2，則K 可以表示為：

對(duì)各個(gè)小段加入窗函數(shù)wd(n)進(jìn)行窗處理，得到第i 段的數(shù)值：

再對(duì)其進(jìn)行離散傅里葉變換，得到：

式（13）中，M 為離散傅里葉變換周期。

最后可以得到各段平均功率譜為：

窗函數(shù)wd(n)可以是矩形窗、漢明窗等，選擇一個(gè)好的窗函數(shù)對(duì)于結(jié)果是很重要的，一般要求窗函數(shù)的主瓣不能太寬，但主瓣太窄會(huì)導(dǎo)致旁瓣下的面積過大。因此，須根據(jù)處理的要求進(jìn)行選擇，平均周期圖法功率譜估計(jì)使得曲線方差特性得到了較大的改善，但是當(dāng)方差較小時(shí)會(huì)導(dǎo)致頻率分辨率較低，二者須權(quán)衡考慮。

2.2 仿真分析

2.2.1 不同下潛深度的潛艇磁異常頻域信號(hào)特征

以我國(guó)周邊海域地磁為背景，令采樣點(diǎn)為256，采樣頻率為10 Hz ，對(duì)南北航向、深度分別為100 m 、200 m、300 m、400 m 的潛艇磁場(chǎng)信號(hào)進(jìn)行平均周期圖法功率譜密度估計(jì)，如圖7所示。

圖7 不同深度處潛艇磁異常信號(hào)功率譜Fig.7 Power spectrum of submarine magnetic anomaly signal at different depths

圖7為潛艇磁場(chǎng)信號(hào)功率譜密度函數(shù)隨深度變化的關(guān)系。在淺水時(shí)，潛艇磁場(chǎng)頻率分量較多、帶寬較寬，隨著潛艇潛航深度的增加，潛艇磁場(chǎng)信號(hào)的高頻分量逐漸減少，低頻分量一直占有信號(hào)的主要能量，曲線趨于平緩，且強(qiáng)度隨著潛艇深度的增加而減弱。我國(guó)附近海域潛艇磁場(chǎng)頻域信號(hào)的主要能量集中也在0～0.25 Hz 內(nèi)，頻率越高，所含有的能量越少，最終趨于平緩并穩(wěn)定在某一個(gè)較低的值，屬于極低頻信號(hào)。

2.2.2 不同運(yùn)動(dòng)速度的潛艇磁異常頻域信號(hào)特征

為研究潛艇航行速度對(duì)潛艇磁異常信號(hào)功率譜的影響，以我國(guó)周邊海域地磁為背景，在下潛100 m 深度處，對(duì)高速（12 n mile）與低速（6 n mile）潛艇的磁異信號(hào)功率譜曲線進(jìn)行比較，如圖8所示。

圖8 不同航速下潛艇磁場(chǎng)信號(hào)功率譜比較Fig.8 Comparison of power spectrum of submarine magnetic anomaly signals at different speeds

從圖8 中可以看出，在水深不變，低速時(shí)，潛艇磁場(chǎng)的頻率成分較少，帶寬較窄，譜峰位于較低頻位置；高速時(shí)，帶寬變寬，出現(xiàn)多個(gè)較小的峰值，且頻率成分向高頻方向拓展。

3 結(jié)論

本文主要對(duì)在我國(guó)周邊海域活動(dòng)的潛艇的磁異常分布、磁場(chǎng)信號(hào)的頻域特征進(jìn)行分析。通過對(duì)潛艇引起的磁異常信號(hào)采用Welch 法進(jìn)行功率譜分析，得到了潛艇磁異常信號(hào)能量與頻率的關(guān)系，仿真分析了潛艇不同下潛深度、不同運(yùn)動(dòng)速度下的磁場(chǎng)頻域信號(hào)特點(diǎn)。因此，利用潛艇磁異信號(hào)的頻域特征能夠更好地對(duì)潛艇進(jìn)行探測(cè)與識(shí)別。