單廣超 趙漢波

(1.海軍陸戰(zhàn)學(xué)院 廣州 510430)(2.國(guó)防科技大學(xué)電子科學(xué)與工程學(xué)院 長(zhǎng)沙 410073)

艦船輻射噪聲線譜檢測(cè)與分析*

單廣超1趙漢波2

(1.海軍陸戰(zhàn)學(xué)院 廣州 510430)(2.國(guó)防科技大學(xué)電子科學(xué)與工程學(xué)院 長(zhǎng)沙 410073)

結(jié)合經(jīng)典功率譜和DEMON譜分析各自的特點(diǎn),對(duì)艦船輻射噪聲線譜進(jìn)行綜合分析?；谥芷趫D譜法進(jìn)行艦船輻射噪聲功率譜估計(jì),通過(guò)線譜與連續(xù)譜分離、取除虛警及歸并線譜,有效地對(duì)輻射噪聲功率譜中的特征線譜進(jìn)行了提取,并結(jié)合DEMON譜分析了艦船輻射噪聲的調(diào)制效應(yīng)和調(diào)制周期,獲得諸如艦船螺旋槳轉(zhuǎn)速、螺旋槳葉片數(shù)等不變的艦船物理特性,從而得到艦船輻射噪聲線譜比較全面的特性,為艦船的識(shí)別和線譜輻射噪聲的控制具有重要提供參考價(jià)值。

周期圖譜法; DEMON譜分析; 船舶輻射噪聲線譜

Class Number TP391.9; TN911.72

1 引言

艦船輻射噪聲的功率譜是由能夠反映隨機(jī)噪聲部分能量分布的寬帶連續(xù)譜和離散的窄帶線譜分量迭加而成的。其中,線譜作為艦船輻射噪聲中的重要成分,不僅能量集中且頻率穩(wěn)定,而且通常分布在低頻段,往往與艦船上的某些機(jī)械振動(dòng)特性相關(guān)聯(lián),從而攜帶有重要的艦船特征信息。一方面,這些強(qiáng)線譜可以用來(lái)估計(jì)艦船的運(yùn)動(dòng)參數(shù),通過(guò)高質(zhì)量的線譜估計(jì)對(duì)目標(biāo)進(jìn)行識(shí)別分類,利用窄帶檢測(cè)系統(tǒng)提高被動(dòng)聲納的作用距離,另一方面,譜特征必然直接與產(chǎn)生它們的現(xiàn)象有關(guān),分析到線譜的準(zhǔn)確頻率,就意味著可能找到該線譜產(chǎn)生源的振動(dòng)周期,這無(wú)疑對(duì)艦船輻射線譜進(jìn)行有源或無(wú)源控制提供了重要依據(jù)。國(guó)外海軍試驗(yàn)表明,低頻線譜噪聲信號(hào)在海中可以傳播相當(dāng)遠(yuǎn)的距離[1]?，F(xiàn)代艦船輻射噪聲中低頻線譜能量比較豐富和低頻信號(hào)遠(yuǎn)距離傳播的事實(shí),為低頻遠(yuǎn)程探測(cè)系統(tǒng)的研究與發(fā)展提供了線譜檢測(cè)的依據(jù)。因此,檢測(cè)和控制艦船輻射噪聲中強(qiáng)線譜分量是一個(gè)十分重要的課題,無(wú)論是對(duì)保持艦船自身的隱蔽性還是遠(yuǎn)距離探測(cè)目標(biāo)都具有重要的意義,同時(shí)也可能為艦船的減震降噪帶新的進(jìn)展。

目前,基于功率譜估計(jì)的特征提取和DEMON譜分析分別是水下目標(biāo)特征提取的兩個(gè)重要方面。信號(hào)的功率譜反映了信號(hào)的許多重要特征,利用功率譜的連續(xù)譜特征和線譜特征進(jìn)行目標(biāo)識(shí)別和分類,是聲納、雷達(dá)、語(yǔ)音識(shí)別和噪聲分析等領(lǐng)域信號(hào)處理的重要內(nèi)容[2～3]。然而,由于螺旋槳噪聲對(duì)空化噪聲等連續(xù)分量存在幅度調(diào)制現(xiàn)象,直接從艦船輻射噪聲功率譜無(wú)法得出螺旋槳的軸頻、葉頻等調(diào)制特征,因此,為了對(duì)艦船輻射噪聲線譜進(jìn)行完整分析還需要通過(guò)DEMON譜分析對(duì)艦船螺旋槳軸頻等不變的物理特征進(jìn)行分析。本文結(jié)合經(jīng)典功率譜和DEMON譜分析各自的特點(diǎn),對(duì)艦船輻射噪聲線譜進(jìn)行綜合分析。通過(guò)線譜提取方法對(duì)輻射噪聲功率譜中的特征線譜進(jìn)行提取,通過(guò)DEMON譜分析提取艦船輻射噪聲的周期調(diào)制譜,獲得諸如艦船螺旋槳轉(zhuǎn)速、螺旋槳葉片數(shù)等不變的艦船物理特性。

2 艦船輻射噪聲線譜理論模型

艦船輻射噪聲功率譜有兩種不同的類型:一種是在噪聲級(jí)上表現(xiàn)為連續(xù)譜的寬帶噪聲,連續(xù)譜可以采用寬帶平穩(wěn)隨機(jī)過(guò)程擬合;另一種是具有非連續(xù)單頻噪聲,這種噪聲由出現(xiàn)在離散頻率上的線譜組成[4],機(jī)械振動(dòng)及螺旋槳轉(zhuǎn)動(dòng)都是周期性的振動(dòng)源,所以用周期信號(hào)作為線譜模型是合適的[5]。艦船的噪聲可以表示為

(1)

式中,{x(t)}表示寬帶平穩(wěn)隨機(jī)過(guò)程;{li(t)}表示初相位隨機(jī)的周期信號(hào),i=1,2,3,…,n。

艦船噪聲的功率譜可以用下式計(jì)算:

(2)

式中:T為進(jìn)行Fourier變換時(shí)截取的每一段信號(hào)長(zhǎng)度,E表示集合平均,K表示信號(hào)段的編號(hào)。

式(2)是艦船噪聲譜在數(shù)學(xué)模型下的嚴(yán)格定義,即要有無(wú)窮多個(gè)信號(hào)段且每個(gè)信號(hào)段的時(shí)間長(zhǎng)度趨于無(wú)窮,但實(shí)際上只能是有限長(zhǎng)和有限數(shù)的集合平均。此時(shí)它的單邊功率譜估計(jì)為

(3)

式中,n為平均次數(shù)。

3 輻射噪聲特征線譜提取

利用現(xiàn)代功率譜和DEMON譜分析各自的特點(diǎn),通過(guò)線譜提取方法對(duì)輻射噪聲功率譜中的特征線譜進(jìn)行提取,通過(guò)DEMON譜分析提取艦船輻射噪聲的周期調(diào)制譜,確定螺旋槳的軸頻和葉頻,進(jìn)而確定螺旋槳葉片數(shù)。

3.1 基于功率譜的線譜分析

功率譜估計(jì)方法一般分為兩大類[6]:經(jīng)典譜估計(jì)和現(xiàn)代譜估計(jì)。經(jīng)典譜估計(jì)是基于Fourier變換的功率譜估計(jì)方法,又可分為“間接法”和“直接法”。相對(duì)于間接法而言,直接法(也稱周期圖法),由信號(hào)序列的Fourier變換直接得到的,具有物理感念清晰、方法簡(jiǎn)便、計(jì)算效率高及分辨率高等優(yōu)點(diǎn),但該方法估計(jì)出的線譜有偽峰,難以辨別特征線譜的真?zhèn)?這里通過(guò)選取合適的線譜分析參數(shù),如信號(hào)段長(zhǎng)度、平均次數(shù)、線譜檢測(cè)門(mén)限等參數(shù),來(lái)克服這些不利因素的影響,進(jìn)行線譜提取。

1) 連續(xù)譜與線譜分離

由于艦船輻射噪聲功率譜中不僅有線譜成分,也有連續(xù)譜噪聲和其它噪聲成分,并且部分線譜疊加在連續(xù)譜上,提取線譜特征時(shí),若直接在含有連續(xù)譜的譜中提取,可能由于連續(xù)譜的趨勢(shì)走向引起誤判和漏判,因此需要先從整個(gè)譜中提取趨勢(shì)項(xiàng),并減去該慢變化的趨勢(shì)項(xiàng),再在剩下部分中提取線譜。

將慢變化分量看成是連續(xù)譜,需要采用背景噪聲平滑技術(shù)估計(jì)出連續(xù)譜和其他噪聲成分的變化趨勢(shì)曲線,通常背景噪聲平滑技術(shù)主要有:線性相位濾波器方法、雙通道分離窗算法、排序截短算法[7]、局部分段擬合方法[8]等。這里采用在空間域和頻率域都十分有效的高斯平滑濾波器提取連續(xù)譜分量的趨勢(shì)項(xiàng),同時(shí)為了克服在信號(hào)邊界由于數(shù)據(jù)量小而導(dǎo)致的平滑失真的現(xiàn)象,因此針對(duì)連續(xù)譜兩段的局部區(qū)域,可采用最小二乘高階擬合進(jìn)行對(duì)應(yīng)趨勢(shì)項(xiàng)提取,從而得到精確的整個(gè)趨勢(shì)項(xiàng)。

2) 譜線提取

在得到連續(xù)譜的變化趨勢(shì)后,從原功率譜中減去這一變化趨勢(shì),即得到拉平后的線譜圖。由于在該線譜圖中存在許多偽峰,對(duì)后續(xù)線譜特征的提取是非常不利的,所以還需要采用進(jìn)一步的邏輯算法對(duì)“偽”線譜進(jìn)行剔除,以凈化線譜。

首先,譜峰往往為所在局部區(qū)域?yàn)樽畲簏c(diǎn),不可能出現(xiàn)在中間點(diǎn)上,由此剔除連續(xù)上升或連續(xù)下降中間點(diǎn),只留下轉(zhuǎn)折點(diǎn)以得到初步簡(jiǎn)化的線譜圖。可采用簡(jiǎn)化方法為[9]:假設(shè)k-1、k、k+1三點(diǎn)的譜值分別為yk-1、yk、yk+1,求一階差分:Δyk=yk+1-yk,Δyk-1=yk-yk-1。若Δyk×Δyk-1<0,則可判定yk點(diǎn)為轉(zhuǎn)折點(diǎn),保留該頻率位置。反之,若Δyk×Δyk-1>0,則yk為中間點(diǎn),將其剔除。

其次,對(duì)于經(jīng)過(guò)初步簡(jiǎn)化的線譜圖進(jìn)一步處理時(shí),若單純依據(jù)卡門(mén)限的方法,會(huì)把一根線判成兩根線,所以需把小變化部分去掉。具體方法為:給定以門(mén)限值Δgate,對(duì)于任意一點(diǎn)yk,如果Δyk-Δyk-1<Δgate或Δyk-Δyk+1<Δgate,則yk不是線譜,應(yīng)進(jìn)行剔除,反之yk為線譜。

最后,對(duì)所剩下的局部最大點(diǎn)進(jìn)行卡門(mén)限處理,將各局部最大點(diǎn)的均值乘一比例因子設(shè)為門(mén)限Δgate′,若yk-2-yk>Δgate′或yk+2-yk>Δgate′,且其頻率間隔小于3Hz,則yk不是線譜,將其剔除,反之保留。由此便可得到了艦船噪聲反而線譜分布。

3.2 DEMON譜分析

由于艦船螺旋槳噪聲存在“螺旋槳拍”現(xiàn)象,是周期調(diào)制譜,直接從艦船輻射噪聲功率譜中提取的線譜無(wú)法得知螺旋槳的軸頻、葉頻等調(diào)制特征,因此需要通過(guò)對(duì)高頻寬帶噪聲信號(hào)進(jìn)行解調(diào)分析獲得線譜信息,即DEMON譜分析(Detection of Envelope Modulation on Noise)。艦船噪聲用周期性局部平穩(wěn)過(guò)程作模型可表示為[10]

s(t)=(1+m(t))*x(t)

(4)

式中,x(t)為窄帶平穩(wěn)高斯隨機(jī)過(guò)程,m(t)是調(diào)制信號(hào),是慢變化的周期函數(shù),它所在的頻域比x(t)所在的頻域要低的多。

通常利用DEMON譜對(duì)螺旋槳的軸頻信息進(jìn)行分析的整個(gè)流程圖如圖1[11]所示。

圖1 DEMON譜分析流程圖

在上述流程圖中,由于通常在低信噪比背景噪聲下,輻射噪聲中連續(xù)譜分量可能將線譜成分淹沒(méi),為此需要對(duì)接收到的艦船輻射噪聲進(jìn)行帶通濾波以得到高頻輻射噪聲中的線譜成分。線性檢波有兩種方法,絕對(duì)值低通檢波和平方低通檢波。雖然兩種方法就極大多數(shù)樣本來(lái)說(shuō),沒(méi)有顯著的差別,但平方低通解調(diào)方法得到的DEMON譜具有更顯著的諧波特性。

4 數(shù)據(jù)處理與分析

為了結(jié)合經(jīng)典功率譜和DEMON譜分析各自的特點(diǎn),對(duì)實(shí)測(cè)的艦船輻射噪聲線譜進(jìn)行綜合分析。下面采用某低速商船實(shí)測(cè)數(shù)據(jù)進(jìn)行分析,為了去除不確定隨機(jī)誤差的影響,選取中間較穩(wěn)定測(cè)量的信號(hào),其長(zhǎng)度為10.8s,采樣頻率為22kHz,采用周期圖譜法得到該商船的功率譜估計(jì)如圖2所示。

為了在線譜提取時(shí),減小由于連續(xù)譜的趨勢(shì)走向引起的誤判和漏判,通過(guò)高斯平滑濾波得到的連續(xù)譜趨勢(shì)項(xiàng)如圖3所示。

圖2 商船輻射噪聲功率譜

圖3 連續(xù)譜趨勢(shì)項(xiàng)

從整個(gè)譜中減去連續(xù)譜趨勢(shì)項(xiàng),將線譜與連續(xù)譜分離,得到拉平的線譜圖,如圖4所示。

圖4 拉平后的線譜

為了剔除噪聲中的虛假“線譜”,根據(jù)文中第3.1節(jié)的方法對(duì)拉平后的線譜進(jìn)一步凈化線譜,選擇集中且能量穩(wěn)定的線譜,根據(jù)文獻(xiàn)[9]的結(jié)論,平均均次數(shù)選擇3,門(mén)限值Δgate取6dB較合適。這里對(duì)所剩下的局部最大點(diǎn)進(jìn)行卡門(mén)限處理時(shí),取門(mén)限Δgate′為4dB。由此可得到該船線譜分布如圖5所示。

圖5 凈化后的線譜分布圖

前面通過(guò)基于周期圖的功率譜方法對(duì)艦船輻射噪聲中線譜成分進(jìn)行提取與分析,有效剔除了過(guò)渡點(diǎn),并對(duì)可能的偽線譜進(jìn)行了剔除,得到了艦船噪聲譜中與艦船物理特性相關(guān)且穩(wěn)定存在線譜分布,這可對(duì)艦船目標(biāo)特性的檢測(cè)與輻射噪聲線譜的控制提供重要參考。同時(shí),由于應(yīng)用功率譜方法并不能從譜圖上看出艦船的軸頻、葉頻等調(diào)制特性,因此,為了對(duì)艦船輻射噪聲線譜進(jìn)行完整分析,下進(jìn)一步通過(guò)DEMON譜分析對(duì)艦船螺旋槳軸頻等不變的物理特征進(jìn)行分析。相對(duì)于直接寬帶解調(diào)而言,采用窄帶解調(diào)平均后的DEMON譜效果更好,這里選取窄帶信號(hào)寬度為250Hz,從1kHz～4kHz共取8個(gè)窄帶,低通濾波截止頻率為100Hz,經(jīng)過(guò)平方解調(diào)后的DEMON譜如圖6所示。

圖6 DEMON譜分析圖

在DEMON譜圖中可以清晰看到對(duì)應(yīng)軸頻基頻和葉頻基頻及其諧波數(shù)。同時(shí),由于該船航速較低,圖5中線譜數(shù)豐富且葉頻基頻諧波數(shù)較多,其中某些線譜也可能是機(jī)械噪聲成份。這樣通過(guò)對(duì)艦船輻射噪聲的DEMON譜分析,可以得出反應(yīng)艦船有關(guān)物理屬性的節(jié)奏信息,即螺旋槳轉(zhuǎn)速和槳葉片數(shù)。

5 結(jié)語(yǔ)

本文結(jié)合經(jīng)典功率譜和DEMON譜分析各自的特點(diǎn),對(duì)艦船輻射噪聲線譜進(jìn)行綜合分析。基于周期圖譜法進(jìn)行艦船輻射噪聲功率譜估計(jì),通過(guò)線譜與連續(xù)譜分離、取除虛警及歸并線譜,有效地對(duì)輻射噪聲功率譜中的特征線譜進(jìn)行了提取,并結(jié)合DEMON譜分析了艦船輻射噪聲的調(diào)制效應(yīng)和調(diào)制周期,獲得諸如艦船螺旋槳轉(zhuǎn)速、螺旋槳葉片數(shù)等不變的艦船物理特性,從而得到艦船輻射噪聲線譜比較全面的特性,為艦船的識(shí)別和線譜輻射噪聲的控制具有重要提供參考價(jià)值。

[1] 胡偉文,肖鋒,苑秉成.目標(biāo)輻射噪聲線譜的特征建模與檢測(cè)[J].海軍工程大學(xué)學(xué)報(bào),2006,18(6):45-48.

[2] Khotanzad A, Lu JH, Serinath M D. TargetDetection U-sing a NeuralNetwork Based on Passive Sonar[J]. Procof IJCNN,1989(1):335-340.

[3] 湯和,劉勝利.噪聲譜的自動(dòng)分析與識(shí)別[J].振動(dòng)工程學(xué)報(bào),1993,6(3):255-263.

[4] 王盼盼,張國(guó)軍,關(guān)凌綱,等.艦船輻射噪聲特征線譜提取方法研究[J].海洋技術(shù),2010,29(3):47-50.

[5] 馬徐坤.淺析水下高速航行體對(duì)螺旋槳輻射噪聲線譜建模[J].聲學(xué)學(xué)報(bào),2002,27(6):502-506.

[6] 黃甫堪,陳建文,樓生強(qiáng).現(xiàn)代數(shù)字信號(hào)處理[M].北京:電子工業(yè)出版社,2003,6.

[7] 劉銘.水下目標(biāo)特征提取與線譜跟蹤技術(shù)研究[D].哈爾濱:哈爾濱工程大學(xué),2011.

[8] 楊向鋒,張效民,龐學(xué)亮,等.基于分段擬合的艦船噪聲線譜特征提取[J].聲學(xué)學(xué)報(bào),2006(4):394-400.

[9] 吳國(guó)清,李靖.艦船噪聲識(shí)別(1)—總體框架、線譜分析和提取[J].聲學(xué)學(xué)報(bào),1998,23(5):88-90.

[10] 熊紫英,朱錫清.基于LOFAR譜和DEMON譜特征的艦船輻射噪聲研究[J].船舶力學(xué),2007,11(2):300-306.

[11] 孫炎.被動(dòng)目標(biāo)分段調(diào)制特征提取方法研究[D].哈爾濱:哈爾濱工程大學(xué),2005.

Detection and Analysis of Ship Radiated-noise Line Spectrum

SHAN Guangchao1ZHAO Hanbo2

(1. Navy Marine-corps Academy, Guangzhou 510430) (2. Depantment of Electromic Science and Engineering, National University of Defense Technology, Changsha 410073)

According to the respective characteristics of classical power spectrum and DEMON spectrum, ship radiated-noise line spectrum is comprehensively analyzed. Based on periodogram algorithm, the spectral estimation of ship radiated-noise is made. After separating line spectrum and continuum spectrum, eliminating fasle alarm and censoring line spectrumm, characteristic line spectrum from the ship radiated-noise is effectivly extracted. Then, the periodic modulation of ship radiated-noise is analyzed in order to obtain the ship invariant physical speciality, such as ship propeller rotation, propeller blade-number, etc. In the end, the comprehensive line spectrum characteristics of ship radiated-noise is obtained, which has important reference value for ship recognition and controlling line spectrum radiated-noise.

periodogram algorithm, DEMON analysis, ship radiated-noise line spectrum

2014年4月7日,

2014年3月12日

單廣超,男,碩士,研究方向:潛艇戰(zhàn)術(shù)。趙漢波,男,博士研究生,研究方向:星基導(dǎo)航與定位技術(shù)。

TP391.9; TN911.72

10.3969/j.issn1672-9730.2014.10.034