唐 兵，朱 艷，陳衛(wèi)東

(中國(guó)船舶重工集團(tuán)公司第七一○研究所，湖北 宜昌 443003)

0 引言

艦船輻射噪聲重構(gòu)技術(shù)是水中兵器仿真研究中的重要分支，一直受到各國(guó)的關(guān)注?，F(xiàn)代人類已知在水中最有效的場(chǎng)信號(hào)是聲信號(hào)，所以艦船在水中的輻射噪聲是對(duì)艦船目標(biāo)進(jìn)行探測(cè)、識(shí)別、定向和跟蹤等目的最常用的目標(biāo)信息。

艦船的輻射噪聲主要由機(jī)械噪聲、螺旋槳噪聲及水動(dòng)力噪聲組成，艦船在不同的航行狀態(tài)下，這些噪聲對(duì)輻射噪聲有不同的貢獻(xiàn)。輻射噪聲包含寬帶連續(xù)譜分量和窄帶線譜分量。連續(xù)譜主要來(lái)自螺旋槳轉(zhuǎn)動(dòng)引起的空化噪聲和艦船機(jī)械噪聲，當(dāng)航速增加時(shí)，螺旋槳空化噪聲級(jí)增加，使得低頻線譜幅度相對(duì)減小甚至被連續(xù)譜淹沒(méi)。線譜主要的噪聲源為艦船的機(jī)械噪聲，其分布情況取決于艦船的類型和航行狀態(tài)。

文獻(xiàn)［1－2］都提出了利用特定頻率響應(yīng)FIR濾波器得到寬帶連續(xù)譜噪聲的方法，該方法能夠產(chǎn)生任意長(zhǎng)度的給定功率譜形狀的噪聲信號(hào)。文獻(xiàn)［3－4］通過(guò)把艦船輻射噪聲的連續(xù)譜與線譜分開(kāi)重構(gòu)再疊加以獲得完整的艦船輻射噪聲重構(gòu)信號(hào)。上述方法中線譜與連續(xù)譜特性的獲得都是通過(guò)對(duì)不同艦船輻射噪聲的共性研究，經(jīng)簡(jiǎn)化擬合而得到的數(shù)理模型。而現(xiàn)代艦船經(jīng)過(guò)一系列的消聲減噪處理，其輻射噪聲差異很大，并處于不斷變化中。所以由上述方法得到的仿真信號(hào)與現(xiàn)代艦船真實(shí)輻射噪聲存在較大差異，不能做到很好地將目標(biāo)分類。同時(shí)將線譜與連續(xù)譜分開(kāi)重構(gòu)也增加了艦船輻射噪聲重構(gòu)的復(fù)雜程度。

隨著被動(dòng)目標(biāo)識(shí)別技術(shù)的發(fā)展，以被動(dòng)聲信號(hào)作為目標(biāo)對(duì)象的水中兵器及探測(cè)裝置的仿真需要越來(lái)越逼真的目標(biāo)噪聲信號(hào)，以便為其研制和驗(yàn)證提供更可靠更有效的仿真手段。本文介紹了利用實(shí)測(cè)艦船輻射噪聲信號(hào)，實(shí)現(xiàn)具有高逼真度的艦船輻射噪聲重構(gòu)。通過(guò)對(duì)實(shí)測(cè)艦船輻射噪聲的分析，探討結(jié)合其功率譜，設(shè)計(jì)特定頻率響應(yīng)FIR濾波器來(lái)實(shí)現(xiàn)艦船輻射噪聲功率譜與連續(xù)譜同時(shí)重構(gòu)的方法與可行性，在一定條件下得到與實(shí)測(cè)艦船輻射噪聲信號(hào)頻域分布相一致的重構(gòu)信號(hào)。

1 實(shí)測(cè)艦船輻射噪聲譜估計(jì)

要實(shí)現(xiàn)基于實(shí)測(cè)艦船輻射噪聲功率譜的艦船輻射噪聲重構(gòu)，首先要得到實(shí)測(cè)艦船的高分辨率功率譜。一般實(shí)測(cè)得到的艦船輻射噪聲都具有通過(guò)特性，并非是一個(gè)平穩(wěn)隨機(jī)信號(hào)，但通過(guò)分段的方式，在時(shí)間足夠短的情況下，可將其近似為平穩(wěn)隨機(jī)信號(hào)，不影響一般性。利用經(jīng)典譜估計(jì)算法中的Welch法對(duì)平穩(wěn)隨機(jī)信號(hào)進(jìn)行譜估計(jì)是目前較為常用的方法［5］。

圖1為某艦船實(shí)測(cè)輻射噪聲信號(hào)，圖1(a)是實(shí)測(cè)具有通過(guò)特性的艦船輻射噪聲的時(shí)域信號(hào)，圖1(b)為對(duì)該噪聲信號(hào)進(jìn)行平穩(wěn)化處理后，即對(duì)噪聲信號(hào)在時(shí)域上進(jìn)行截短處理，再利用Welch法做譜估計(jì)得到的該艦船實(shí)測(cè)輻射噪聲的功率譜圖。從圖中可以看出，該方法所得到的功率譜圖能夠在全頻段上很好地表現(xiàn)出艦船輻射噪聲的功率譜特性——低頻線譜加上連續(xù)譜，并為后面的艦船輻射噪聲重構(gòu)提供有效的依據(jù)。

Welch法所得到的信號(hào)功率譜分辨率正比于2π/N，功率譜數(shù)據(jù)長(zhǎng)度為N/2，其中N為對(duì)信號(hào)作FFT變換的長(zhǎng)度。為使得其功率譜在低頻段的線譜能夠被識(shí)別出來(lái)，要求N取一個(gè)較大的值。N值的增大，增加了所得到的功率譜長(zhǎng)度，同時(shí)增加了以此為期望頻率響應(yīng)的FIR濾波器設(shè)計(jì)中的頻率約束點(diǎn)，并增加了FIR濾波器的復(fù)雜度。

圖1 某艦船實(shí)測(cè)噪聲信號(hào)Fig.1 Measure signal of ship noise

為了減少設(shè)計(jì)FIR濾波器的復(fù)雜程度即減少計(jì)算所需時(shí)間，有必要在不影響功率譜特性的情況下減少功率譜的數(shù)據(jù)長(zhǎng)度。由于艦船輻射噪聲的線譜多集中在低頻段(通常為0～1 kHz)，高頻段主要是連續(xù)譜，所以功率譜在高頻段的分辨率要求并不高。這樣可以通過(guò)對(duì)信號(hào)功率譜高頻段進(jìn)行降采以使得到的實(shí)測(cè)艦船輻射噪聲功率譜數(shù)據(jù)變短同時(shí)能保持原有的特征。

圖2為重新采樣后的艦船輻射噪聲功率譜曲線與原始功率譜曲線對(duì)比圖。從對(duì)比圖可以看出，降采后的功率譜曲線與原曲線基本一致，在低頻段能夠很好地表現(xiàn)線譜，在高頻段同時(shí)能夠表現(xiàn)其趨勢(shì)。

圖2 降采后功率譜與噪聲原功率譜對(duì)比Fig.2 Contrast picture of sampling and original power spectrum

2 艦船輻射噪聲重構(gòu)

艦船輻射噪聲可用寬帶平穩(wěn)過(guò)程來(lái)擬合。線性系統(tǒng)理論中，當(dāng)一個(gè)平穩(wěn)隨機(jī)信號(hào)的取樣序列通過(guò)一個(gè)線性系統(tǒng)后，所得到的系統(tǒng)輸出仍是平穩(wěn)隨機(jī)信號(hào)。根據(jù)以上理論，設(shè)計(jì)一個(gè)具有與實(shí)測(cè)艦船輻射噪聲頻譜特性相一致的頻率響應(yīng)的線性系統(tǒng)，輸入高斯白噪聲，以得到平穩(wěn)隨機(jī)信號(hào)的取樣序列，其頻譜特性與線性系統(tǒng)的頻率響應(yīng)相一致。

根據(jù)上述思路，文中構(gòu)建了一個(gè)有限沖擊響應(yīng)(FIR)濾波器，以實(shí)測(cè)艦船輻射噪聲頻譜特性作為其期望的頻率響應(yīng)，并利用該濾波器對(duì)高斯白噪聲濾波以獲得艦船輻射噪聲的重構(gòu)信號(hào)。圖3為艦船輻射噪聲重構(gòu)原理圖。

圖3 艦船輻射噪聲重構(gòu)原理圖Fig.3 Principle drawing of ship radiated noise reconstruction

根據(jù)圖3所示，首先需獲得服從高斯分布且具有純白功率譜的高斯白噪聲隨機(jī)序列。高斯白噪聲采用Box-Muller方法:由相互獨(dú)立的均服從N(0，1)分布的隨機(jī)數(shù)β1和β2，可以得到服從N(0，1)分布、幅值服從高斯分布、具有純白功率譜的隨機(jī)數(shù):

圖4為利用上述方法得到的高斯白噪聲時(shí)間序列。其幅值在0附近波動(dòng)。

圖4 高斯白噪聲時(shí)域信號(hào)Fig.4 Time domain signal of Gaussian white noise

特定頻率響應(yīng)濾波器的設(shè)計(jì)是艦船輻射噪聲重構(gòu)的關(guān)鍵，其頻率響應(yīng)直接決定了重構(gòu)信號(hào)的頻譜特性。為使重構(gòu)信號(hào)與實(shí)測(cè)艦船噪聲信號(hào)在頻域上分布一致，文中以實(shí)測(cè)艦船輻射噪聲的頻譜特性作為設(shè)計(jì)FIR濾波器的期望頻率響應(yīng)。前面利用Welch法已經(jīng)得到信號(hào)的功率譜，由數(shù)字信號(hào)處理可知功率譜與頻率響應(yīng)有如下關(guān)系［5］:

式中:P(k)為信號(hào)的功率譜;N為FFT變換的長(zhǎng)度;XN(k)為信號(hào)傅立葉變換，即信號(hào)的頻譜。

從而得到FIR濾波器的期望頻率響應(yīng)Hd(ejω)。在已知頻率響應(yīng)的情況下，本文利用窗函數(shù)法設(shè)計(jì)FIR濾波器。由頻率特性知，由傅立葉逆變換得到該濾波器的單位抽樣響應(yīng)為

這樣，h(n)是因果的，且有限長(zhǎng)，由此可以計(jì)算出所設(shè)計(jì)的濾波器的轉(zhuǎn)移函數(shù)為

當(dāng)濾波器的階數(shù)足夠長(zhǎng)，即截取的h(n)增加時(shí)，H(z)逐漸逼近Hd(ejω)。通過(guò)式(3)～式(5)，就能夠設(shè)計(jì)出與期望頻譜特性具有一致頻率響應(yīng)的FIR濾波器。在Matlab中可以利用函數(shù)FIR2來(lái)實(shí)現(xiàn)上述算法。

3 仿真結(jié)果

以第2節(jié)中使用的某艦船實(shí)測(cè)輻射噪聲信號(hào)為實(shí)例，在計(jì)算機(jī)上利用Matlab對(duì)上述方法進(jìn)行數(shù)學(xué)仿真。圖5為設(shè)計(jì)出的FIR濾波器頻率響應(yīng)與艦船噪聲頻譜特性對(duì)比圖。從圖中可以看出，設(shè)計(jì)的濾波器頻率響應(yīng)與艦船頻譜特性基本一致，在低頻段能夠比較好地反應(yīng)出實(shí)測(cè)噪聲的線譜特征，高頻段與實(shí)測(cè)艦船噪聲頻譜特性相吻合，滿足設(shè)計(jì)要求。

圖5 FIR濾波器頻率響應(yīng)與艦船噪聲頻譜特性對(duì)比Fig.5 Frequency response of FIR filter and frequency spectrum of ship noise contrast picture

圖6 重構(gòu)信號(hào)與實(shí)測(cè)信號(hào)的功率譜對(duì)比Fig.6 Reconstruction signal and measure signal power spectrum contrast picture

利用設(shè)計(jì)得到的濾波器對(duì)高斯白噪聲濾波得到重構(gòu)的艦船輻射噪聲。然后利用Welch法對(duì)重構(gòu)信號(hào)進(jìn)行譜估計(jì)并與實(shí)測(cè)艦船輻射噪聲的功率譜相對(duì)比。圖6是重構(gòu)信號(hào)與實(shí)測(cè)信號(hào)的功率譜對(duì)比圖?？梢?jiàn)，該方法能夠很好地實(shí)現(xiàn)基于實(shí)測(cè)噪聲數(shù)據(jù)的艦船輻射噪聲重構(gòu)，重構(gòu)信號(hào)的功率譜不僅能夠很好地?cái)M合艦船輻射噪聲的連續(xù)譜，同時(shí)也能夠很好地表現(xiàn)出低頻線譜。

4 結(jié)語(yǔ)

本文通過(guò)對(duì)實(shí)測(cè)艦船輻射噪聲的分析，利用設(shè)計(jì)特定頻率響應(yīng)FIR濾波器的方法，實(shí)現(xiàn)了艦船連續(xù)譜與線譜的同時(shí)重構(gòu)，仿真結(jié)果表明，重構(gòu)信號(hào)在頻域上能夠很好地與實(shí)測(cè)噪聲信號(hào)相吻合。該方法使得重構(gòu)信號(hào)能夠更加逼真地反應(yīng)艦船輻射噪聲頻譜特性，能夠?yàn)橐耘灤晥?chǎng)作為目標(biāo)識(shí)別手段的探測(cè)裝置的研制、試驗(yàn)和驗(yàn)證提供新的更加逼真和可靠的仿真信號(hào)源。

［1］向輝平，羅建，傅瑞錦.艦船寬帶連續(xù)譜噪聲的模擬［J］.聲學(xué)技術(shù)，2005，9(24):140 －143.

［2］羅建，程越位，相敬林.艦船輻射噪聲波形重構(gòu)研究［J］.魚(yú)雷技術(shù)，2010，12(18):431 －435.

［3］李琴，苑秉成，張文娟.艦船輻射噪聲建模及仿真模擬器的實(shí)現(xiàn)［J］.艦船科學(xué)技術(shù)，2010，32(4):121 －124.

［4］陳曦，林懷清.基于時(shí)變調(diào)制的艦船輻射噪聲信號(hào)仿真算法研究［J］.系統(tǒng)仿真學(xué)報(bào)，2012，24(2):297－300.

［5］胡廣書(shū).數(shù)字信號(hào)處理［M］.北京，清華大學(xué)出版社，2003.