留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

PCA特征提取和弹性BP神经网络的水下气泡识别

张银波 李思宁 姜鹏 孙剑峰

张银波, 李思宁, 姜鹏, 孙剑峰. PCA特征提取和弹性BP神经网络的水下气泡识别[J]. 红外与激光工程, 2021, 50(6): 20200352. doi: 10.3788/IRLA20200352
引用本文: 张银波, 李思宁, 姜鹏, 孙剑峰. PCA特征提取和弹性BP神经网络的水下气泡识别[J]. 红外与激光工程, 2021, 50(6): 20200352. doi: 10.3788/IRLA20200352
Zhang Yinbo, Li Sining, Jiang Peng, Sun Jianfeng. Underwater bubbles recognition based on PCA feature extraction and elastic BP neural network[J]. Infrared and Laser Engineering, 2021, 50(6): 20200352. doi: 10.3788/IRLA20200352
Citation: Zhang Yinbo, Li Sining, Jiang Peng, Sun Jianfeng. Underwater bubbles recognition based on PCA feature extraction and elastic BP neural network[J]. Infrared and Laser Engineering, 2021, 50(6): 20200352. doi: 10.3788/IRLA20200352

PCA特征提取和弹性BP神经网络的水下气泡识别

doi: 10.3788/IRLA20200352
详细信息
    作者简介:

    张银波,男,博士生,主要从事激光雷达水下信号处理方面的研究

    李思宁,男,教授,博士生导师,博士,主要从事激光成像雷达技术及图像处理方面的研究

  • 中图分类号: TN959

Underwater bubbles recognition based on PCA feature extraction and elastic BP neural network

  • 摘要: 针对水下激光雷达探测得到的尾流回波信号由于非稳态造成特征提取困难、不易识别的问题,提出了基于PCA特征提取与弹性BP神经网络结合的水下气泡识别算法。首先对连续采集的回波信号进行切片预处理,然后采用PCA算法对拼接的高维样本进行主要特征提取,确定特征值个数,其次对弹性BP神经网络进行参数的选择,确定能实现最优分类的隐含层节点数、特征个数等,最后根据室内搭建的尾流探测模拟平台,实现对气泡群和干扰目标的识别。实验结果表明:在隐含节点为12,增量因子为1.15,减量因子为0.55时,选取两个特征值能对有气泡、无气泡及干扰物进行有效分类;识别率随着气泡群密度的增大提升13.4%,在低密度下的识别率随激光能量的增加平均提升6.3%,识别率随距离的增加先增大后减小,气泡群在2.2 m时的目标峰特征明显,平均识别率提升3.5%。通过与自适应附加动量BP对比,该方法在减少识别时间的同时准确率达到99.1%,证明该算法可有效运用于激光雷达舰船尾流气泡的识别。
  • 图  1  (a) 水下气泡数据采集实验装置图;(b) 水下激光雷达系统(集成在金属盒内);(c) 没有背景光时激光穿过气泡群;(d) 气泵(电压:220 V,频率:50 Hz,功率:60 W,排气量:50 L/min);(e) 气泡板(气泡的直径约为10~200 μm);(f) 阀门:控制气流

    Figure  1.  Experimental device for data acquisition of underwater bubbles; (b) Photo of the underwater lidar system (inside the metal box); (c) The laser runs through the bubbles without background light; (d) Air pump (voltage: 220 V, frequency: 50 Hz, power: 60 W, transmission volume: 50 L/min); (e) Bubbles plate (the diameter of the air bubbles is about 10-200 μm); (f) Valve: control airflow

    图  2  气泡群回波信号

    Figure  2.  Echo signals of bubbles

    图  3  前2个主成分得分分布图

    Figure  3.  Score distribution of the first two principal components

    图  4  不同算法的迭代收敛情况。(a) 弹性BP算法;(b) 自适应与附加动量BP算法

    Figure  4.  Iterative convergence of different algorithms. (a) Elastic BP algorithm; (b) Adaptive and additional momentum BP algorithm

    图  5  PCA+弹性BP神经网络的水下气泡识别流程图

    Figure  5.  Underwater bubbles recognition process based on PCA and elastic BP neural network

    图  6  不同条件下的回波信号及识别情况。(a) 不同目标回波曲线;(b) 低密度的气泡群识别率;(c) 高密度气泡群识别率

    Figure  6.  Echo signals and recognition results under different conditions. (a) Different targets echo curves; (b) Low density bubbles recognition rate; (c) High density bubbles recognition rate

    表  1  系统指标参数

    Table  1.   System parameters

    ParameterValueParameterValue
    Wavelength/nm532Receiving diameter/mm80
    Pulse width/ns10Emission diameter/mm20
    Pulse energy/mJ10Emission angle/mrad1.7
    Center distance/mm100Receiving angle/mrad2.5
    下载: 导出CSV

    表  2  不同隐含节点数的训练情况

    Table  2.   Training of different hidden nodes

    Number of hidden nodesTraining timesConvergence time/s
    31451.713.43
    61079.432.29
    9685.861.43
    12653.861.43
    下载: 导出CSV

    表  3  不同累计贡献率下的识别结果

    Table  3.   Recognition results under different cumulative contribution rates

    Number of eigenvalues21128
    Cumulative contribution rate39.7%71.2%90.6%
    No bubbles(correct)200200200199200198198195199
    Bubbles(correct)200199198199200200198198198
    Glass plate(correct)196197194198196199198195193
    Recognition rate99.3%99.3%98.7%99.3%99.3%99.5%99.0%98.0%99.0%
    Average recognition rate99.1%99.4%98.7%
    下载: 导出CSV

    表  4  不同算法识别对比结果

    Table  4.   Recognition and contrast results of different algorithms

    MethodRecognition
    rate
    Time(training+
    recognition)/s
    Adaptive and additional
    momentum BP
    98.3%48.94
    PCA+Adaptive and
    additional momentum BP
    96.6%10.15
    Elastic BP98.8%1.60
    PCA+elastic BP99.1%1.36
    下载: 导出CSV
  • [1] Rao Jionghui, Ge Weilong, Zhu Donghua, et al. Ship wake detecting system based on backscattering light of bubbles [J]. Journal of Naval University of Engineering, 2008, 20(4): 62-65. (in Chinese)
    [2] Wan Jun, Zhang Xiaohui, Rao Jionghui, et al. Processing of backscattering signal of warship wake flow based on indep- endent component analysis [J]. Infrared and Laser Enginee- ring, 2013, 42(1): 244-250. (in Chinese) doi:  10.3969/j.issn.1007-2276.2013.01.044
    [3] Peng Xiaolei, Ma Aoling, Liu Yimin. The study on backward optical detection method for ship wake bubbles [J]. Ship Science and Technology, 2016, 38(11A): 133-135. (in Chinese)
    [4] Liang Shanyong, Wang Jiang’an, Zong Siguang, et al. Laser detection method of ship wake bubbles based on multiple scattering intensity and polarization characteristics [J]. Acta Physica Sinica, 2013, 62(6): 87-97. (in Chinese)
    [5] Han Biao, Liu Jifang, Liu Kunlun, et al. Study of backward optical detection method for ship wake bubbles [J]. Acta Optica Sinica, 2012, 32(1): 0101001. (in Chinese) doi:  10.3788/AOS201232.0101001
    [6] Liu Tao, Wang Jiang’an, Zong Siguan, et al. Experimental study of laser-generated cavitation bubble motion near a free liquid surface [J]. Acta Optica Sinica, 2012, 32(7): 0714003. (in Chinese) doi:  10.3788/AOS201232.0714003
    [7] Asante-Okyere S, Shen C, Ziggah Y Y. Principal component analysis (PCA) based hybrid models for the accurate estimation of reservoir water saturation [J]. Computers & Geosciences, 2020, 145: 104555.
    [8] Wang S Q, Zhou G A, Song B J, et al. Research on radar emitter signal feature extraction method based on fuzzy entropy [J]. Procedia Computer Science, 2019, 154: 508-513. doi:  10.1016/j.procs.2019.06.081
    [9] Xie Bing, Duan Zhemin. UAV target recognition algorithm based on fusion of SAE and bottom visual feature [J]. Infrared and Laser Engineering, 2018, 47(S1): S126004. (in Chinese) doi:  10.3788/IRLA201847.S126004
    [10] Zhao Dongbo, Li Hui. Radar target recognition based on central moment feature and GA-BP neural network [J]. Infrared and Laser Engineering, 2018, 47(8): 0826005. (in Chinese)
    [11] Jin Baoming, Lu Guangyi, Wang Wei, et al. Research on BP neural network rainfall runoff forecasting model based on elastic gradient descent algorithm [J]. Journal of Shandong University(Engineering Science), 2020, 50(3): 117-124. (in Chinese)
    [12] Wang Lili, Liu Hongbo, Chen Deyun, et al. Identification of flow regimes based on adaptive learning and additional momentum BP neural network for electrical capacitance tomography [J]. Journal of Harbin University of Science and Technology, 2018, 23(1): 105. (in Chinese)
    [13] Guo Wenzhong, Huang Yuansheng. Fault diagnosis of rotor based on improved models of resilient back -propagation neural network [J]. Gas Turbine Technology, 2009, 22(1): 61-65. (in Chinese) doi:  10.3969/j.issn.1009-2889.2009.01.015
  • [1] 黄宜帆, 贺岩, 胡善江, 侯春鹤, 朱小磊, 李凯鹏, 刘芳华, 陈勇强, 郭守川.  海洋激光雷达图像处理提取海水深度的方法 . 红外与激光工程, 2021, 50(6): 20211034-1-20211034-8. doi: 10.3788/IRLA20211034
    [2] 杨程, 鄢秋荣, 祝志太, 王逸凡, 王明, 戴伟辉.  基于深度学习的压缩光子计数激光雷达 . 红外与激光工程, 2020, 49(S2): 20200380-20200380. doi: 10.3788/IRLA20200380
    [3] 沈振民, 赵彤, 王云才, 郑永超, 尚卫东, 王冰洁, 李静霞.  混沌脉冲激光雷达水下目标探测 . 红外与激光工程, 2019, 48(4): 406004-0406004(7). doi: 10.3788/IRLA201948.0406004
    [4] 张浩, 李向春, 杨倩, 吴承璇, 雷卓.  水下气泡光学图像识别方法 . 红外与激光工程, 2019, 48(3): 326001-0326001(7). doi: 10.3788/IRLA201948.0326001
    [5] 高健, 周安然, 孙东松, 郑俊, 李梓霂, 韩於利.  多普勒激光雷达中种子激光注入的一种鉴别方法 . 红外与激光工程, 2018, 47(2): 230001-0230001(7). doi: 10.3788/IRLA201847.0230001
    [6] 王强, 郝利丽, 唐红霞, 李贤丽, 牟海维, 韩连福, 赵远.  实际环境对量子激光雷达性能的影响 . 红外与激光工程, 2018, 47(S1): 29-35. doi: 10.3788/IRLA201847.S106006
    [7] 任文核, 宋小全, 王芳涵.  星载激光雷达观测东中国海气溶胶时空分布的统计特征 . 红外与激光工程, 2018, 47(12): 1230002-1230002(6). doi: 10.3788/IRLA201847.1230002
    [8] 杜玉红, 王鹏, 史屹君, 王璐瑶, 赵地.  环境特征自适应激光雷达数据分割方法 . 红外与激光工程, 2018, 47(8): 830001-0830001(8). doi: 10.3788/IRLA201847.0830001
    [9] 夏文泽, 韩绍坤, 曹京亚, 王亮, 翟倩.  激光雷达距离估计技术 . 红外与激光工程, 2016, 45(9): 906005-0906005(6). doi: 10.3788/IRLA201645.0906005
    [10] 吴超, 刘春波, 韩香娥.  光波导相控阵激光雷达接收系统设计 . 红外与激光工程, 2016, 45(10): 1030003-1030003(6). doi: 10.3788/IRLA201645.1030003
    [11] 单坤玲, 刘新波, 卜令兵, 郜海阳, 黄兴友.  激光雷达和毫米波雷达的卷云微物理特性的联合反演方法 . 红外与激光工程, 2015, 44(9): 2742-2746.
    [12] 史风栋, 刘文皓, 汪鑫, 丁娟, 史屹君, 修春波.  室内激光雷达导航系统设计 . 红外与激光工程, 2015, 44(12): 3570-3575.
    [13] 曹开法, 黄见, 胡顺星.  边界层臭氧差分吸收激光雷达 . 红外与激光工程, 2015, 44(10): 2912-2917.
    [14] 杨成华, 张勇, 赵远.  基于点阵照明的APD阵列外差激光雷达外差效率的提高 . 红外与激光工程, 2014, 43(10): 3269-3275.
    [15] 靳辰飞, 胡鹏, 张思琦, 刘丽萍, 赵远, 谭立英.  增益调制激光雷达弱信号的双阈值处理分析 . 红外与激光工程, 2014, 43(8): 2447-2452.
    [16] 赵培娥, 罗雄, 曹文勇, 赵彬, 冯立天, 李晓锋, 谭锦, 周鼎富.  应用Zoom FFT方法提高相干测风激光雷达频谱分辨率 . 红外与激光工程, 2014, 43(1): 98-102.
    [17] 姜成昊, 杨进华, 张丽娟, 李祥.  新型多普勒成像激光雷达原理设计与仿真 . 红外与激光工程, 2014, 43(2): 411-416.
    [18] 杨文秀, 付文兴, 周志伟, 余巍, 马杰.  基于投影降维的激光雷达快速目标识别 . 红外与激光工程, 2014, 43(S1): 1-7.
    [19] 孙崇利, 苏伟, 武红敢, 刘睿, 刘婷, 黄健熙, 朱德海, 张晓东, 刘峻明.  改进的多级移动曲面拟合激光雷达数据滤波方法 . 红外与激光工程, 2013, 42(2): 349-354.
    [20] 梁善勇, 张伟, 王江安, 吴荣华, 马治国.  舰船尾流共轴激光雷达近场抑制比的计算方法 . 红外与激光工程, 2012, 41(8): 2063-2067.
  • 加载中
图(6) / 表(4)
计量
  • 文章访问数:  21
  • HTML全文浏览量:  17
  • PDF下载量:  3
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-11-20
  • 修回日期:  2021-01-25
  • 刊出日期:  2021-06-30

PCA特征提取和弹性BP神经网络的水下气泡识别

doi: 10.3788/IRLA20200352
    作者简介:

    张银波,男,博士生,主要从事激光雷达水下信号处理方面的研究

    李思宁,男,教授,博士生导师,博士,主要从事激光成像雷达技术及图像处理方面的研究

  • 中图分类号: TN959

摘要: 针对水下激光雷达探测得到的尾流回波信号由于非稳态造成特征提取困难、不易识别的问题,提出了基于PCA特征提取与弹性BP神经网络结合的水下气泡识别算法。首先对连续采集的回波信号进行切片预处理,然后采用PCA算法对拼接的高维样本进行主要特征提取,确定特征值个数,其次对弹性BP神经网络进行参数的选择,确定能实现最优分类的隐含层节点数、特征个数等,最后根据室内搭建的尾流探测模拟平台,实现对气泡群和干扰目标的识别。实验结果表明:在隐含节点为12,增量因子为1.15,减量因子为0.55时,选取两个特征值能对有气泡、无气泡及干扰物进行有效分类;识别率随着气泡群密度的增大提升13.4%,在低密度下的识别率随激光能量的增加平均提升6.3%,识别率随距离的增加先增大后减小,气泡群在2.2 m时的目标峰特征明显,平均识别率提升3.5%。通过与自适应附加动量BP对比,该方法在减少识别时间的同时准确率达到99.1%,证明该算法可有效运用于激光雷达舰船尾流气泡的识别。

English Abstract

    • 随着世界各国对海上军事装备的加强,对于舰船、航母、核潜艇的侦查显得尤为重要,舰船尾流的识别受到了国内外学者的广泛关注。海军工程大学[1]设计研制并验证了基于后向散射的光尾流探测方法的可行性;万俊等[2]采用盲源分离和独立成分分析的方法对水体与气泡尾流后向散射光信号进行分离,但得到的信号噪声对判决结果有很大影响;彭晓雷等[3]通过对气泡回波信号的后沿位置、峰值大小和位置等脉冲特征进行统计分析,实现了对气泡的识别,但由于采集数据多造成脉冲特征统计过于繁琐;梁善勇等[4]利用回波强度、偏振信息等实现对舰船尾流的探测识别。综上可以看出,目前采用激光雷达对舰船尾流的识别主要依靠回波信号特征和偏振信息[2-4]等,这些方法是对探测结果较为明显的单帧回波信号进行特征提取,并未对连续采集的多帧动态回波信号进行分析。

      水中气泡群对激光有很强的后向散射效应[5],入射到气泡表面与透过气泡群的回波信号有很大差别。气泡直径、密度、运动速度时刻变化,使得回波信号处于一个非线性、非稳态状态[6],导致气泡群回波信号特征不稳定。因此,利用传统方法对非稳态的尾流回波信号进行特征提取和识别,存在很大的局限性。

      针对以上问题,文中提出了基于主分量分析法 (Principal Component Analysis, PCA)[7]与弹性BP (back propagation)神经网络结合的水下气泡识别算法。通过自主设计的水下激光雷达探测系统进行室内模拟实验,针对非稳态的回波信号进行数据预处理得到满足条件的高维特征矩阵,利用主分量分析法对预处理后的拼接样本进行特征提取,采用弹性BP神经网络进行水下气泡识别,可有效解决连续采集的非稳态回波信号特征提取困难和识别率低的问题。实验结果表明,该算法在减少识别时间的同时,对气泡群有较好的识别效果。

    • 实验数据经笔者所在实验室自主研发的水下单元APD激光雷达探测系统得到,如图1(a)所示为搭建的激光雷达室内气泡群探测平台,主要分为3个部分:lidar系统结构、水池环境、数据采集结构。数据采集结构包括采样频率为500 MHz的AD采集卡和计算机。水池环境包括内壁贴满黑色卡纸的透明玻璃水池(0.48 m×0.67 m×3.70 m),水深0.36 m,激光探测高度0.26 m。利用电磁式空气压缩机与增氧气盘产生气泡,气泡直径约为10~200 μm,并用气阀对气流量进行调节。雷达系统指标参数见表1

      图  1  (a) 水下气泡数据采集实验装置图;(b) 水下激光雷达系统(集成在金属盒内);(c) 没有背景光时激光穿过气泡群;(d) 气泵(电压:220 V,频率:50 Hz,功率:60 W,排气量:50 L/min);(e) 气泡板(气泡的直径约为10~200 μm);(f) 阀门:控制气流

      Figure 1.  Experimental device for data acquisition of underwater bubbles; (b) Photo of the underwater lidar system (inside the metal box); (c) The laser runs through the bubbles without background light; (d) Air pump (voltage: 220 V, frequency: 50 Hz, power: 60 W, transmission volume: 50 L/min); (e) Bubbles plate (the diameter of the air bubbles is about 10-200 μm); (f) Valve: control airflow

      表 1  系统指标参数

      Table 1.  System parameters

      ParameterValueParameterValue
      Wavelength/nm532Receiving diameter/mm80
      Pulse width/ns10Emission diameter/mm20
      Pulse energy/mJ10Emission angle/mrad1.7
      Center distance/mm100Receiving angle/mrad2.5

      在测量过程中,信号的起始时间由激光器内触发信号决定,回波信号是以512个采样点为一个周期的连续一维曲线。系统对位于2.8 m处的气泡群进行连续探测,将前5个周期的有效回波信号截取出来,如图2所示。曲线1、2、3呈双峰分布,表明此刻激光系统探测到气泡群回波,由于激光器、AD采集卡具有一定延时以及气泡群的运动,导致回波信号产生时间起点上的相对漂移,连续采集的回波信号在幅值和峰值位置具有一定差异;曲线4呈单峰分布,即只有水体回波信号,由于系统的低采样率和气泡群的非稳态特性,造成此刻未接收到气泡群回波信号,这类数据的存在降低了气泡群回波信号的识别率;曲线5没有回波峰,为无效数据,需在预处理中去除。因此,数据中存在有气泡群回波峰的信号、无气泡群回波峰的信号和空采的无效数据等。与传统的经过人为选择方式得到满足条件的单帧回波信号相比,文中得到的原始数据与实际舰船尾流探测得到的回波信号更加接近,对该情况下的气泡群回波进行探测识别更具有实用性。

      对目标回波信号的特征提取一直是进行目标识别与分类的主要技术[8-9],回波漂移、采样率低和无效采集等问题对于气泡群的识别有极大影响,因此必须对测量结果做相应的预处理,具体步骤如下:

      (1)原始数据分段切片。根据水池长度及回波信号时间漂移量可知,回波信号主要集中在前80个采样点,为实现后续的波形平移,因此,截取每个分段的前100个采样点作为单次采集的有效数据,共采集N组数据。

      (2)去除无效数据。首先将各组数据中小于零值的强度置零,求得峰值求平均值$ E\left( x \right)$,确定峰值有效范围,若超出有效范围,则去除该组数据。经多次实验,有效范围为$ \left( {1 \pm 0.5} \right) \times E\left( x \right)$时,可对无效数据进行有效去除。

      (3)波形平移。设第i次采集的信号为$ x\left( i \right)$,峰值位置为$ a\left( i \right)$,找到N组数据中峰值位置最大值MAX和最小值MIN,设基准点为$ n$,满足公式(1)的条件。每次以步长1增加,寻找最优基准点使得各$ a\left( i \right)$与基准点$ n$的平方差$ E\left( n \right)$达到最小,最后将$ x\left( i \right)$$ a\left( i \right)$$ n$的相对距离在时域上进行左右平移,数据空位按零值补位。

      $$ {\rm MIN} \leqslant n \leqslant {\rm MAX} $$ (1)

      (4)阈值处理。在不失真的前提下选择强度阈值,将小于阈值的置零,减少噪声的影响,文中的阈值为最大峰值的5%。

      (5)将定长的矢量数据按采集时间顺序进行拼接,形成一个高维特征矩阵。

      图  2  气泡群回波信号

      Figure 2.  Echo signals of bubbles

    • PCA是将多个指标化为少数几个综合指标来代替原始数据中较多的变量的统计方法,而且所含信息互不重复。为验证该算法的识别能力,在实验中增加一类干扰目标(透明玻璃板),将其置于气泡群所处的位置,得到其回波信号与气泡接近。分别对距离激光雷达2.2 m的无气泡、气泡、透明玻璃板3类目标进行连续1 min测量,测得的回波信号组数分别为589、593、578;然后对各组数据进行预处理,得到有效数据各为545、560、569;随机选取其中各300组数据,每组数据按行排列,以100个采样点作为每组数据的特征维数,最后得到由维数100的900组数据拼接组成的高维样本矩阵X900×100

      为选取能有效识别的特征,利用PCA算法对样本矩阵进行降维处理,特征值以权重大小依次往后排列。一般认定当前n个特征值累计权重值足够大,可取前n个主分量作为提取后的样本特征向量。

      为最大化减少特征维数,选择前2个特征值的得分,第1个作为横坐标,第2个作为纵坐标画离散图,如图3所示。气泡主要分布于PC1轴正半轴区域,无气泡和玻璃板主要分布于PC1轴负半轴区域,除个别没有完全分离以外,3类目标基本分为3个部分,且气泡与非气泡的区别比较明显。通过前2个主成分可以将气泡识别出来,说明将预处理后的数据拼接为高维样本矩阵,利用PCA进行特征提取是有效的。

      图  3  前2个主成分得分分布图

      Figure 3.  Score distribution of the first two principal components

    • BP神经网络[10]以其能够处理知识背景模糊、环境信息复杂等系列问题, 被广泛应用于信号处理、模式识别等领域。

      目前,弹性梯度下降法能有效避免梯度大导致收敛速度慢的问题[11]。为消除梯度幅度的影响,权值修正时仅与梯度的符号有关,2次迭代的梯度方向决定权值大小。梯度方向相同时,权值的修正值乘以一个增量因子;相反时,乘以一个减量因子;梯度等于零时,保持不变。权值修正的更新公式为:

      $$ \begin{split} \begin{array}{*{20}{c}} {\Delta w(i + 1) = \Delta w(i + 1) \cdot sign(gra(i))}=\\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{\!\!\! \left\{ {\begin{array}{*{20}{c}} {\!\!\!\!\! \Delta w(i) \times \alpha \times sign(gra(i))}\\ {\!\!\!\!\! \Delta w(i) \times \beta \times sign(gra(i))}\\ {\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\! \Delta w(i)} \end{array}} \right.} \end{array} \end{split} $$ (2)

      式中:gra(i)为第i次迭代的梯度;$ \alpha $为增量因子;$ \;\beta$为减量因子,$ 0 < \beta < 1 < \alpha $

      在相同条件下,设定目标函数为0.001,迭代5000次,得到弹性BP算法和自适应学习与附加动量BP算法[12]随迭代次数的增加的收敛情况,如图4所示,横坐标为训练次数,纵坐标为归一化后数据的均方误差。弹性BP算法在迭代415次时就达到目标,而自适应与附加动量BP在迭代4996次时还未达到目标函数,弹性BP算法收敛速度快且收敛性更好。

      图  4  不同算法的迭代收敛情况。(a) 弹性BP算法;(b) 自适应与附加动量BP算法

      Figure 4.  Iterative convergence of different algorithms. (a) Elastic BP algorithm; (b) Adaptive and additional momentum BP algorithm

      因此,利用PCA算法与弹性BP神经网络结合的气泡群识别过程如图5所示。

      图  5  PCA+弹性BP神经网络的水下气泡识别流程图

      Figure 5.  Underwater bubbles recognition process based on PCA and elastic BP neural network

    • 首先对隐含节点数的确定,采用经验公式[13]为:

      $$ i = \sqrt{ {j + l} + a} $$ (3)

      式中:$ j$为输入节点;$ l$为输出节点;$ a$为1~10之间的调节常数。文中输入神经元是经PCA提取的特征个数,为2,输出神经元个数为1,得到隐含节点范围在(3,12)之间,学习率为0.05,训练5000次,目标函数为0.001。选择不同的节点数逐个计算,每组进行7次实验取平均值,选出最优节点数,如表2所示,随着隐含节点的增加,训练次数减少,收敛用时也减少。综合考虑,确定隐含节点数为12。

      表 2  不同隐含节点数的训练情况

      Table 2.  Training of different hidden nodes

      Number of hidden nodesTraining timesConvergence time/s
      31451.713.43
      61079.432.29
      9685.861.43
      12653.861.43

      将输入层和输出层的激发函数设置为线性函数,隐含层的激发函数设置Sigmoid函数,增量因子、减量因子由实验结果所决定,其中在$ \alpha =1.15$$ \;\beta=0.55 $时有最优输出。无气泡、气泡和透明玻璃板分别用0、1、−1表示,输出偏差阈值设定为0.2,即输出偏差大于0.2时表示类别不确定,反之表示识别正确。

    • 选择对距离激光雷达2.2 m的无气泡、气泡和透明玻璃板的回波信号共3类目标进行识别,随机选取每类各300组,共计900组数据作为训练集,另随机选取各200组共计600组作为测试集。选择特征值的不同累计贡献率,得到3组不同测试集的识别率如表3所示,选择特征值个数为2时,其识别率可以达到99.1%,累计贡献率提升到71.2%时,识别率只提高0.3%。随着特征值个数的增加,由于在识别过程中进行划分的特征维数的增加,造成特征区分不明显,最终导致识别率低。因此,在进行PCA算法和弹性BP神经网络识别的过程中,特征个数对识别结果有很大的影响。综合考虑,特征值个数为2时为最优情况。

      表 3  不同累计贡献率下的识别结果

      Table 3.  Recognition results under different cumulative contribution rates

      Number of eigenvalues21128
      Cumulative contribution rate39.7%71.2%90.6%
      No bubbles(correct)200200200199200198198195199
      Bubbles(correct)200199198199200200198198198
      Glass plate(correct)196197194198196199198195193
      Recognition rate99.3%99.3%98.7%99.3%99.3%99.5%99.0%98.0%99.0%
      Average recognition rate99.1%99.4%98.7%

      为验证PCA算法和弹性BP神经网络结合的优越性,采用自适应与附加动量BP和无PCA进行降维的网络作为对比,对比结果如表4所示。利用神经网络对预处理的数据进行识别,弹性BP的用时明显少于自适应与附加动量BP,PCA降维的用时少于不降维时间,其中PCA+弹性BP识别率高达99.1%。对于大量水下气泡数据的进行处理时,该算法能体现出更明显的优势。

      表 4  不同算法识别对比结果

      Table 4.  Recognition and contrast results of different algorithms

      MethodRecognition
      rate
      Time(training+
      recognition)/s
      Adaptive and additional
      momentum BP
      98.3%48.94
      PCA+Adaptive and
      additional momentum BP
      96.6%10.15
      Elastic BP98.8%1.60
      PCA+elastic BP99.1%1.36

      为证明该算法的泛化能力,改变激光器电流实现激光发射能量的调节,控制气阀进气量实现气泡群密度的变化,选择对距离激光雷达1.6 m、2.2 m和2.8 m的气泡群进行探测,另外选用水体和透明玻璃板的回波信号作为干扰信号。当激光器电流为72 mA,气泡群和玻璃板距离为2.2 m时,经预处理后得到的回波信号如图6(a)所示,目标的存在使得回波具有双峰现象,并且回波强度明显大于水体回波;气泡密度较小时的目标峰与玻璃板回波一致,强度小于水体回波峰;随着气泡密度的增加,目标峰值强度明显增大。

      图  6  不同条件下的回波信号及识别情况。(a) 不同目标回波曲线;(b) 低密度的气泡群识别率;(c) 高密度气泡群识别率

      Figure 6.  Echo signals and recognition results under different conditions. (a) Different targets echo curves; (b) Low density bubbles recognition rate; (c) High density bubbles recognition rate

      将预处理后的气泡群回波信号与水体和玻璃板回波信号混合,随机抽取部分数据进行训练和识别。经算法处理后,得到气泡群在不同位置和激光能量下的识别情况,结果如图6(b)(c)所示。当气泡密度较低时,由于玻璃板的回波信号与其相似,导致识别率小于高密度时的识别率约13.4%;激光器电流从72 mA增加至74 mA时,在低密度下的识别率平均提升6.3%,高密度下变化不明显。由于水体回波峰一直存在,当气泡群距离为1.6 m时,气泡回波峰与水体回波峰重叠,单峰现象明显增多,导致识别率降低;当距离增大至2.2 m时,回波信号为双峰,此时目标峰特征明显,平均识别率提升约3.5%,识别结果较为稳定;随着距离增加至2.8 m时,由于水体对激光具有强散射性质,导致回波能量较低,对气泡群回波峰的有效探测次数降低,识别率减小。

    • 对于尾流气泡的探测,利用水下激光雷达进行连续采集,得到多帧动态回波信号,可以实现对非稳态气泡群的强度的检测,将原始一维信号切片预处理能减小气泡群非稳态造成的回波漂移、采样率低和无效采集等影响,利用主分量分析法对拼接的高维样本进行降维,减少了神经网络输入的维数,不但有99.1%的识别精度,还有更快的学习速度,对3类目标回波信号进行识别,结果表明PCA+弹性BP算法对水下气泡识别的有效性。改变激光发射能量和气泡群密度,实现对不同距离的气泡群探测,在有干扰信号的情况下,识别率随着气泡群密度的增大提升13.4%,在低密度下的识别率随激光能量的增加平均提升6.3%,识别率随距离的增加先增大后减小,气泡群在2.2 m时的目标峰特征明显,平均识别率提升3.5%。此次实验的样本维数太少,没有完全体现出PCA特征提取和弹性BP神经网络的优越性,后续工作将会展开大量多维样本的识别。

参考文献 (13)

目录

    /

    返回文章
    返回