基于辐射噪声能量分布的舰船目标分类算法
2021-03-17
来源:飒榕旅游知识分享网
第32卷第10期 2010年lO月 舰船科学技术 Vo1.32.No.10 Oct..2010 SHIP SCIENCE AND TECHNOLOGY 基于辐射噪声能量分布的舰船目标分类算法 鞠熠昊,耿伯英,刘 健,王文涛 (海军3-程大学电子工程学院,湖北武汉430033) 摘 要: 有效的特征提取技术是水中目标识别的基础。为提高基于舰船辐射噪声的水中目标识别准确率,选 用小波变换完成信号预处理和滤波,并在信号变换后的多尺度子空间上提取信号特征参数,归一化处理后构建分类 特征向量,最后用支持向量机算法进行训练和测试。仿真结果表明,利用小波变换的多分辨率分析方法和支持向量 机算法对舰船辐射噪声信号进行分类识别,特征提取算法有效,分类速度较快。 关键词: 小波变换;特征提取;支持向量机;目标识别 中图分类号:TB182;TP319 文献标识码:A 文章编号: 1672—7649(2010)10—0045—05 DOI:10.3404/j.issn.1672—7649.2010.10.012 Algorithm of ship target recognition based on energy distribution of the ship radiated noise JU Yi—hao,GENG Bo—ying,LIU Jian,WANG Wen—tao (Electrical Engineering College,Naval University of Engineering,Wuhan 430033,China) Abstract:The effeetive feature extraction is the basis for underwater target recognition.To enhance the accuracy of underwater target passive recognition based on the ship radiated signals,the signals of ship radiated are preprocessed and filtered by WT and the energy feature parameters are extracted upon multi— scale subspace after WT.At last,the characteristic vector is established with normalization process and the signals are classified using SVM.The simulation experiment results demonstrate that,the ship radiated signals can be classified or recognized availably based on WT and SVM.The classification speed rapidly. Key words: wavelet transform;feature extraction;suppo ̄vector machine;target recognition 0 引 言 水中目标识别 是通过水听器(如被动声呐)接 目标信号 (含噪声) 类别判定 收信号(包括目标的辐射噪声以及与其所处的水介 质共同作用后产生的环境噪声)的处理,从而对辐射 源进行类别判定。在民用和军用领域(如海洋环境 图l 水中目标识别的总体框架 Fig.1 The framework of the underwater target recognition 开发、海洋生物研究、海底资源开发、武器开发研制等 方面)都有广泛的应用价值和重要意义。水中目标 特征提取与识别是水中目标识别的关键技术之~,也 是国内外一直公认的难题。作为一类特殊的模式识 别问题,水中目标识别的总体框架如图1所示。 键环节。由于海洋环境的复杂性和水声信号通道的 特殊性,要从舰船辐射噪声信号中提取1种既能反映 目标本质特征,又能满足水中探测要求的有效特征表 示,一直是水中目标识别领域的难题和热点 。小 波变换(Wavelet Transform,wT)不仅具有良好的时频 舰船辐射噪声信号是水中目标分类识别的主要 信息来源,其特征信号的提取是水中目标识别中的关 定位特性,通过它提取的特征可以更好地反映信号的 时间与尺度特性,而且小波变换后信号相当于在不同 收稿日期:2010—04—07;修回日期:2010—05—04 基金项目:海军工程大学自然科学基金(HGDJJ20O8029) 作者简介:鞠熠吴(1980一),男,硕士研究生,讲师,主要研究方向为军用目标识别技术。 ・46・ 舰船科学技术 第32卷 频带上进行滤波处理,更容易准确 地提取物理特 征。 待分析信号,(t)的小波变换定义是把基本小波 (t)作位移r后,再在不同尺度a下与待分析的信 号.厂(t)作内积: 统计学习理论(Statistical Learning Theory, SLT)Es]是一种小样本统计理论,为研究有限样本情 况下的统计模式识别和更广泛的机器学习(Machine Learning,ML)问题建立了1个较好的理论框架,同时 (n,'.)= f,4(£) f 11d£,口>0。(3) a 一 、 0 在实际应用中,往往需要把连续小波及其变换离 散化,对连续小波的参数o,b进行采样;口=口 ,b: 也发展了1种新的基于支持向量机(Suppo ̄Vector Machine,SVM)的模式识别方法。SVM是Vapnik于 1995年提出的1种新的机器学习方法 ,由于其具 有良好的泛化能力,在解决小样本、非线形及高维模 式识别问题中表现出许多特有优势,被推广应用到函 数拟合等其他机器学习问题中。目前,SVM已成为 机器学习域新的研究热点并已被应用于人脸、文本、 手写体识别等多个研究工作中" 。 为提高基于舰船辐射噪声的水中目标识别准确 率,本文首先利用小波变换对舰船辐射噪声信号进行 预处理和滤波,然后在小波变换后的各层子空间上提 取能量特征参数,再进行归一化处理,构建出信号的 特征向量,最后将原始测量信号分为训练部分和测试 部分,运用C—SVM的分类算法对舰船辐射噪声信号 进行了分类识别。 1 算法原理 1.1小波变换理论 小波分析(Wavelet Analysis)是在20世纪80年 代,由Morlet在分析研究地理物理信号时提出的,是 1种分析非平稳信号的非常有效方法。它克服了短 时傅立叶变换固定分辨率的缺点,既能分析信号的整 个轮廓,又可以进行信号细节的分析,实现了对信号 不同分辨率下的分析。 小波函数定义:设 (t)为一平方可积函数,即 (f)∈L ( ),若其傅立叶变换 (to)满足条件: f—j一 2 d <∞, (1)1 J R 则称 (t)为1个基本小波或小波母函数,并称式 (1)为小波的可容性条件。原则上,任何满足可容性 条件的 ( )空间的函数都可作为小波母函数。 将小波母函数 (t)进行伸缩和平移,设其伸缩 因子(也称尺度因子)为a,平移因子为 ,令其伸缩 平移后的函数为: …(£)=a一丁 ( —J_),a>0,7-∈R, (2) Ⅱ 称 …(t)为依赖于参数a和.r的小波基函数。 n60Ⅱ ,从而: . (z)=口一T (ao t—n6o), (4) .+∞ (m, )=fJ 一∞ £) …(t)dt= t), (‘)}。(5) 因此,小波变换可以描述为函数 t)∈ ( )。 通过带通滤波器的滤波输出响应。研究证明 ,人 耳之所以能有效区分声信号,在于耳蜗处理复杂声信 号的能力及内在机理,即在低频时使用较窄的频域 窗,提高谱分辨率;在高频时使用较宽的频域窗,降低 谱分辨率,其实质就是1个小波变换。因此,利用小 波变换进行声信号特征描述,并用于计算机识别,具 有可靠的理论依据和实践基础。 1.2多分辨率分析 ’ 多分辨率分析(Multi—Resolution Analysis,MRA) 由S.Mallat和Y.Meyer于1986年提出,它将此前所 有的正交小波基的构造统一起来,并在计算机上实现 了正交小波变换的快速算法(Mallat算法),使小波理 论研究产生突破性的进展。 MRA的主要思想是用不同的时间分辨率来逐级 逼近待分析的信号_厂( )。利用多分辨率分析,可以将 1个已知函数 t)投影到相互正交的函数空间中。 尺度相同的小波函数或尺度函数构成了1个函数空 间的基,而不同函数空间其基的尺度是不相同的,从 而表征了不同的频率特性。通过不断地投影,即可逐 步得到原函数在其任意子空间中的分量,由粗及细地 实现信号的分析。这里把信号分为近似部分和细节 部分,所谓近似部分,即指信号中可以投影到大尺度 (即尺度值较大)空间中的低频分量,而细节部分则 指信号中和小尺度空间相对应的高频分量。 1.3支持向量机 支持向量机用于分类问题可描述为:给定1组训 练数据( ‘,Y1),i=1,…,z; ∈R ,YiE{一1,+I}, 要寻找1个分类规则,( ),使它能对未知类别的新 样本(新样本与训练样本独立同分布)作尽可能正确、 的划分。其实质就是寻找1个最优分类超平面作为 分类决策面;同时还通过引进核函数巧妙地解决了低 第10期 鞠熠昊,等:基于辐射噪声能量分布的舰船目标分类算法 ・47・ 维空间向量映射到高维空问向量时带来的“维数灾 难”问题。 在训练集线性可分情形下,SVM就是要构造1 个最优超平面: (09・ )+b=0。 (6) 这个超平面既要满足下面的约束条件: Y [( ・ )+b]≥1,i=1,…,f; (7) 同时还要使下面的函数取得最小值: 1 1 (to)=÷llot lI =÷(to・ )。 (8) 通过求解最优化问题可得最优超平面的形式如 下: ∑YiOL?( ・ )+b。=0。 (9) 其中:SV为支持向量,Ol 为Lagrange乘子。 在训练集线性不可分时,引进松弛因子 >0及惩 罚参数C。此时需在约束Y [( ・ )+b]≥1一 , f 1 i:1,…,z下求最小化函数 ( )=÷ 『一 l+c’∑ ,I 类似可得最优超平面。有了最优超平面,分类规则或 分类函数只要取,( ):sgn i∑Yi ?( ‘Xi)+b。I V 即可。 支持向量机特点之一在于核函数的引入。众所 周知,低维空间向量集往往难于划分,若把它们映射 到高维空间,随之而来的是计算复杂度的大大增加, 核函数巧妙地解决了这个问题。 若函数 ( ・Y)满足Mercer条件 ,则K( ・Y) = ( )・ (Y)。其中西()表示某个映射(未必知其 具体表达式)。这样,只要适当选取核函数就可以得 到对应高维空间的分类函数: ,( )=sgn l∑YSV i ?K( ~ )+b。 。(10) 其中: ( ),(b(Y)为比 ,Y更高维的向量(注意:不 必知道其具体形式),由于K( ・Y)只涉及( ・Y), 因此计算没有涉及高维运算。 为此,我们将分类函数(决策函数)类型为式 (10)的学习机称为支持向量机。 由上可知,支持向量机的基本思想为:在进行分 类时,对于线性不可分的样本首先通过一个非线性映 射将原空间的样本映射到1个高维的特征空间(也 称核空问),使得在核空问中变得线性可分或近似线 性可分,然后在核空间中进行线性分类(或其他线性 算法),从而实现相对于原空间的非线性算法。 2 舰船辐射噪声的特征提取 2.1 小波变换提取多尺度空间能量特征基本原理 按照小波变换理论,基于小波变换提取多尺度空 间能量特征的基本原理可描述为:小波变换可在“放 大”了的不同频带内分析信号,使本不易察觉的信号 特征在不同分辨率的若干子空间中显露出来,具备数 字“显微镜”功能。因此,选取不同的小波函数,分解 结果也不同。但无论小波函数如何选取,每一分解尺 度所用的滤波器中心频率与带宽成固定的比例,即具 有所谓“恒Q”特性。因此,各尺度空间内的平滑信 号和细节信号提供了原始信号的时频局域信息,特别 是能提供不同频段上信号的构成信息。由于能量是 反映信号特征的一个重要物理量,求解出不同分解尺 度上信号能量,将这些能量值按尺度顺序排列成特征 向量以备分类识别使用。 2.2基于舰船辐射噪声的水中目标小波分解能量特 征提取 对于舰船辐射噪声来说,信号所包含的能量频谱 分布与目标类型、形状和大小密切相关。因此,小波 分解后尺度空间上的能量分布如同舰船的低频线谱 一样,反映出目标的本质特征,可用于分类识别 。 因此可选取各尺度上细节信号的能量作为特征向量, 主要原因在于:①不同目标辐射信号在各个尺度上 的细节能量分布不同;②小波变换后的各个尺度上 的细节信号是正交的,因此它们是不相关的,符合目 标特征向量最小维数的要求;③用归一化能量作为 特征向量,既可以消除原始信号幅值的影响,也易于 使特征向量归一化。故采用小波变换方法进行目标 特征提取,从多分辨率分析出发,对信号进行Mallat 多级小波分解,可得出信号在各级近似空间和细节空 问上的信号表示。 鉴于此,基于舰船辐射噪声的水中目标小波分解 能量特征提取步骤如下: 1)利用信号的多尺度分解及其快速算法对舰船 辐射噪声信号进行小波变换; 2)重构小波变换后的各层细节信号; 3)在第1级到第n级细节子空问上,分别计算 这n个子空间中的信号能量; 4)将这n个子空间中的信号能量进行归一化处 理,作为分类特征向量。 这样,经过/1,级小波分解后,便可提取出特征向 量,并且其特征向量具备以下2个优点:一是特征维 ・48・ 舰船科学技号名术 0.5 0。0.5第32卷 i 数低,因为小波分解的级数Ⅳ一般不会太高;二是选 用此方法提取的特征包含了信号在不同频带k的信 息,并未造成信息的大量损失。 2.3特征提取实例 一 ——— ——— }!!!!!蔷 。。 重豳一 0.O.50 5 以2种类型舰船辐射噪声信号为例进行多尺度小 写0‘5波分解,用22 kHz频率对原始信号采样,采样时间为 1 s,信号长度为22 050个点,取db5小波,进行5级小 写一0.5 ———菇 0—-— —■■●—■ 波分解。图2和图3分别表示噪声的原始波形及分解 结构。图4和图5分别表示各级小波细节信号。 图2类型1原始信号及小波分解结构 Fig.2 Type 1 structure of the original signal and wavelet decomposition 1 O —O -1 (S)第2类原始信号 (b)低频系数和各层高频系数 图3类型2原始信号及小波分解结构 Fig.3 Type 2 structure of the original signal and wavelet decomposition 计算各层细节信号的能量并进行归一化处理,构 成5维小波变换后的尺度一能量特征向量,如图6所 示。 2种类型舰船辐射噪声经小波变换的多分辨率 分析后提取的典型特征向量见表1,以db5小波的5 层分解为例。 2.4结果分析 从图表中可以看出:2种不同舰船辐射噪声信号 小波变换后的不同尺度细节信号能量值的分布明显 0.2 亳_().2 O.O5 O -0.o5 0 0.5 1_0 1.5 2.0 xlo4 时间 s 图4类型1小波细节信号 Fig.4 Type 1 wavelet detail signal 0.5 . .....1一 号一o-9}!!!!! 0 _口 0 时间t/s 图5类型2小波细节信号 Fig.5 Type 2 wavelet detail signal 归一化能量 归一化能量 图6 归一化后的尺度一能量特征向量 Fig.6 Normalized scale—the energy feature vector 第l0期 鞠熠昊,等:基于辐射噪声能量分布的舰船目标分类算法 表1分类特征向量 ・49・ Tab.1 Classificati0n of the vector 舰船辐射噪声 、 —二]信号 / 归l 一一 化能量 小波变换 ) / 构建尺度能 量特征向量 一多分辩挈分辑 ) 计算各尺度空 、—二二[二 ( 一SVM分类器 一 一 、 一旦 ~一一 丝垦— ~ ’—]—:一一 图7计算总能量 ) ( _ 分类结果 _二 水中目标分类流程不意图 Fig.7 Flow diagram of underwater target classiicatifon 不同。表明了不同目标的物理特性,而这些特征具有 相对的稳定性,可用于舰船目标的分类识别。 利用LIBSVM¨ “ 进行分类,采用C—SVC分类类型, RBF核函数,测试结果见表2。其中,最优参数是采 用LIBSVM软件自带选优工具所得的结果;KNN— SVM是采用文献[11—12]提出的算法。 分类实验结果表明: 此外,通过分析可知:①选取不同的小波基函数 有着不同的特征向量,因此小波基函数的选取将直接 影响分类结果。②小波变换后的小波系数表明了小 波与被处理信号之间的相似程度。如果小波变换后 的小波系数较大,表明小波和信号的波形相似程度较 大;反之则较小。如果小波变换要反映信号整体的近 似特性,往往选用较大的尺度;要反映细节上的变换 则选用小尺度。由于小波函数家族成员多,进行小波 变换的目的各异,目前如何选取小波基函数还没有一 个通用的标准,不过在实际应用中已取得了一些经 验,具体可参考文献[9]。 1)信号的能量结构能真实反映舰船辐射噪声信 号的本质特征。 2)SVM分类算法能进行有效的分类,平均识别 率在70%以上,分类速度较快。 3)若进行核函数参数优选或采用KNN—SVM算 法优化,可进一步提高分类的准确率。 4 结 语 本文探讨了基于舰船辐射噪声的水中目标小波 特征提取及SVM分类方法。针对水中目标识别的特 3 基于舰船辐射噪声的水中目标分类实验 根据小波变换和支持向量机的基本原理,在实际 应用中,具体多尺度子空问能量分布特征的提取及目 标分类识别过程的流程如图7所示。 殊性,选用目标能量分布特征进行特征提取,直接提 取所有分解后的尺度空问能量,对目标先验信息的依 赖较弱。提取的能量特征形式简单且特征维数低,算 法简便、存储及运算量小、速度快。实测信号的仿真 测试结果表明,基于小波特征提取和支持向量机分类 选用实测舰船噪声辐射信号,以22 kHz频率对 原始信号采样,采样时间为1 s,信号长度为22 050个 点。待分类的信号为2类,各有120个样本,其中36 个样本作为训练样本,其他84个样本作为测试样本。 算法能有效识别舰船辐射噪声类型。 表2分类实验结果 Tab.2 Classi cati0n results 特征提取方法 db5小波4层分解 。 70.987490 0.005472 0.005824 0.069692 sV数目 72 72 66 15 25 28 257 271 237 0.2857 1.3095 1.7143 scc 0.7500 0 2086 0.0769 92.857l 67.26l9 57 1429 95.2381 84.5238 79.7619 73.21 70.24 db5小波5层分解 db5小波6层分解 —69.462890 35.6O1947 dmey小波4层分解 dmey小波5层分解 —7O.955385 0.006905 0.005459 72 72 44 34 193 256 0.3333 0.3571 0 7143 0.6970 91 6667 89 8810 93.4524 91.0714 71.127677 32.23l379 dmey小波6层分解 0.074806 60 37 l94 1.7143 0.0769 57.1429 88.0952 64.88 f下转第83页) 第10期 李 响,等:舰船地震波场检测系统及实验研究 如 ∞ ∞ ∞ m 5 O ・83・ 1.5 耗和大容量存储的难题,并取得了丰富的海上试验数 1.4 据。数据分析表明:舰船地震波场信号有着明显的频 h I1.3 -_^上‘ 乜IJ 谱特性,有利于舰船目标的探测识别。传感器性能的 之1 2 提高,测量方法的改进以及信号处理方法的丰富,为 ‘”P ’ 再’ rr 1.1 Ir, 以后深入研究舰船地震波场打下了良好的基础。 lI 1.0 参考文献: 0.9 [1] 傅金祝.国外水雷战装备参考资料[M].宜昌:710所, O.8 0 5 10 15 20 25 30 35 l994.1—80. 时间/s [2] 张胜业,潘玉玲.应用地球物理学原理[M].武汉:中国 地质大学出版社,2004,299—370. [3] 肖汉光,蔡从中,廖克俊.利用声波和地震波识别军事车 辆类型[J].系统工程理论与实践,2006,4(4):108— J脚 113. [4]李响,颜冰.舰船地震波场分析[J].噪声与振动控制, 2007,27(4):l20—122. [5] 李响,颜冰.海底沉积物特性对舰船地震波传播的影响 [J].舰船科学技术,2008,30(4):43—46. [6]HECHT R J.Investigation of the Potentialities of Using 图7实验数据及其时频图 Seismic Sensors for the Detection of Ships and Other Naval Fig.7 Time and ̄equency analysis of the signal Platforms[R].AD repo ̄,AO63218,1976. [7] 蔡素雄.MSP430F149在电力测控保护产品中的应用[J]. 3 结 语 单片机与嵌入式系统应用,2004,25(8):58—59,74. 『8]Texas In,trument¥.MSP430XIXX User Guide[M].USA: 设计的舰船地震波场检测系统成功解决了微功 Texas Instruments,2002. (上接第49页) [8] 张静远,张冰,蒋兴舟.基于小波变换的特征提取方法分 参考文献: 析[J].信号处理,2000,16(2):156—162. [9] 葛哲学,沙威.小波分析理论与MATLAB R2007实现 [1] 章新华.水中目标自动识别的核心技术[J].兵工学报, [M].北京:电子工业出版社,2007. 1998,19(3):275—280. [10] CHANG C C。LIN C J.LIBSVM:a library for suppo ̄ [2] 郭丽华,王大成,丁士圻.水中目标特征提取方法研究 vector machines[EB/OL].http://www.csie.ntu.edu.tw/ [J].声学技术,2005,9(3):148—152. ~cjlin/libsvm. [3] 丁玉薇.被动声纳目标识别技术的现状与发展[J].声 SHEN Xiao—qiao,LIN Ya—ping.Gene Expression Data 学技术,2004,23(4):253—257. Classiifcation Using SVM—KNN Classifier[A].Proceedings [4] 章新华,王骥程,林良骥.基于小波变换的舰船辐射噪声 of 2004 International Symposium on Intelligent Multimedia, 特征提取[J].声学学报,1997,22(2):139—144. Video and Speech Processing[c].Hong Kong,2004,149一 [5] CRISTIANINI N.支持向量机导论[M].李国正,译.北 l52. 京:机械工业出版社,2004. [12] 刘健,薛蒙.基于K最近邻决策的支持向量机分类算法 [6] 张学工.关于统计学习理论与支持向量机[J].自动化 及仿真[J].舰船电子工程,2009,29(3):88—91. 学报,2000,26(1):32—42. [13] LIBSVM软件,http://www. csie. ntu. edu.tw/~ [7] YANG X W,et a1.Auditory representation of acoustic cjlin/libsvm/acknowledgements. signals[J].IEEE trans.on IT,1992,38:824—839.