技术领域
[0001] 本
发明涉及超声探测领域的
信号参数估计技术，更具体地说，尤其涉及一种基于卷积神经网络的窄带超声回波个数的估计方法。
背景技术
[0002] 超声探测成像技术是当今社会
无损检测技术领域中一种非常重要的手段和方法。它是在不损坏材料
工件或成品的条件下，利用材料内部结构异常或
缺陷存在所引起的对声、磁、光、电、热场等反应的变化，来探测各种工程材料、零部件、结构件等内部和表面缺陷，并对缺陷的类型、性质、数量、形状、
位置、尺寸、分布及其变化做出判断和评价。无损检测诊断的目的在于定量掌握缺陷与强度的关系，评价构件的允许
载荷、寿命或剩余寿命，检测设备(构件)在制造或使用过程中产生的结构不完整性，以及缺陷的情况，从而方便改进制造工艺，提高产品
质量，及时发现故障，保证设备安全、高效、可靠地运行。它在机械制造、石油化工、舰艇
船舶、
汽车、
铁道、建筑、
冶金、航空航天和医学等众多领域均有广泛的应用。
[0003] 在实际中，由于材料结构的复杂性和噪声对目标回波信号的影响，实测中所接收的回波信号总会发生一定的变化，这些变化主要表现为
频率消散、
相位漂移或者是相对参考回波发生偏移等。而且，随着超声射入目标的深度增加，回波衰减将进一步加重。另外，对于超声回波的多层材料检测技术，主要是应用层与层之间的时延来估计层的厚度，但随着层厚度减小，不同层之间的反射回波会发生重叠，造成时间测量困难和误差，这种情况还会表现在材料内部两个距离很近或小尺寸缺陷上，由于缺陷回波的严重重叠，无法有效成像，进而影响判断缺陷的个数、大小及准确位置。目前，回波参数的估计基本都是通过信号分离技术、参数化方法或者时频分析技术去实现，这些方法在复杂环境中或者
信噪比低的情况下，估计误差比较大。因此，亟待发明一种能对在复杂的环境下获取到的窄带超声回波信号进行有效的回波个数估计的方法。
发明内容
[0004] 本发明的目的在于提供一种基于卷积神经网络的窄带超声回波个数的估计方法，利用该方法能够有效的估算出窄带超声回波信号的回波个数，具有鲁棒性好、
精度高的特点，还具备学习能
力。
[0005] 本发明采用的技术方案如下：
[0006] 一种基于卷积神经网络的窄带超声回波个数的估计方法，包括以下步骤：
[0007] (1)设定窄带超声回波信号模型参数，在不同噪声条件下生成多个回波信号集模型；
[0008] (2)在步骤(1)所生成的多个回波信号集模型随机提取一部分回波信号集模型进行短时傅里叶变换并生成时
频谱图作为训练样本集，将剩下的回波信号集模型进行短时傅里叶变换并生成时频谱图作为测试样本集；
[0009] (3)将步骤(2)中获得的训练样本集输入卷积神经网络进行训练，并输出训练好的卷积神经网络；
[0010] (4)将步骤(2)中获得的测试样本集输入步骤(3)中训练好的卷积神经网络，即可估计出对应回波信号的回波个数。
[0011] 进一步的，在所述的步骤(1)中，具体包括以下步骤：
[0012] (1a)根据窄带超声回波的物理特性设定窄带超声回波信号集初始模型参数，窄带超声回波信号集初始模型参数的计算公式为：
[0013]
[0014] 其中，θ＝[α,τ,fc,φ,β]，
[0015] α为窄带超声回波的带宽，τ为窄带超声回波的到达时间，fc为窄带超声回波的中心频率，φ为窄带超声回波的相位，β为窄带超声回波的幅度系数，t为窄带超声回波的传递时间；
[0016] 再将窄带超声回波信号集初始模型参数中的窄带超声回波信号逐个做加性白噪声处理后得回波信号集模型参数，窄带超声回波信号做加性白噪声处理的计算公式为：
[0017]
[0018] 其中，w(t)为加性高斯白噪声，M是窄带超声回波信号集初始模型参数中的窄带超声回波个数。
[0019] (1b)设定噪声信噪比取值范围、步长值以及仿真次数，根据蒙特卡洛法，生成多个回波信号集模型。
[0020] 进一步的，在所述的步骤(1a)中，α的取值范围为2-6MHz，τ的取值范围为0.5-10s，2
fc的取值范围为2-10MHz ，φ的取值范围为0-2rad，β的取值范围为0-2，M的取值范围为1-10个。
[0021] 进一步的，在所述的步骤(1b)中，噪声信噪比的取值范围为-10dB-30dB，步长为1dB，仿真次数为100次。
[0022] 进一步的，在所述的步骤(2)中，具体包括以下步骤：
[0023] (2a)在步骤(1)所生成的多个回波信号模型中提取一部分回波信号集模型作短时傅里叶变换获取回波信号x(t)的时频谱图，作为训练样本集，所述的短时傅里叶变换公式为：
[0024]
[0025] 其中，t为窄带超声回波的传递时间，γ(t)为的分析窗函数,γ*(t)为γ(t)的复共轭函数，取STFT(τ,f)的模|STFT(τ,f)|为x(t)的时频谱图；
[0026] (2b)将步骤(1)所生产的多个回波信号模型中剩下的回波信号模型作短时傅里叶变换获取回波信号x(t)的时频谱图，作为测试样本集，所述的短时傅里叶变换公式为：
[0027]
[0028] 其中，t为窄带超声回波的传递时间，γ(t)为的分析窗函数,γ*(t)为γ(t)的复共轭函数，取STFT(τ,f)的模|STFT(τ,f)|为x(t)的时频谱图。
[0029] 进一步的，在步骤(2)中，在步骤(1)所生成的多个回波信号集模型中随机提取80％的回波信号集模型进行短时傅里叶变换并生成时频谱图作为训练样本集，将步骤(1)所生产的多个回波信号模型中剩下的回波信号模型作短时傅里叶变换并生成时频谱图作为测试样本集。
[0030] 进一步的，在所述的步骤(3)中，具体包括以下步骤：
[0031] (3a)所述的卷积神经网络依次设为
输入层、隐层和全连接层三层结构；
[0032] (3b)将训练样本集中所有时频谱图进行预处理为的灰度图像，并对应生成回波个数标签集，将训练样本集中所有时频谱图作为标签输入层的输入，对应生成的回波个数标签集作为处理输出；
[0033] (3c)所述的隐层依次设为第一卷积层、第一下
采样层、第二卷积层以及第二
下采样层这四层结构，所述的第一卷积层以灰度图像作为输入，采用4个5
像素×5像素的卷积核进行处理输出，所述的第一下采样层以第一卷积层的输出作为输入，采用窗函数大小为2像素×2像素作最大
池化抽样处理输出，所述的第二卷积层以第一下采样层的输出作为输入，采用8个卷积核进行处理输出，所述的第二下采样层以第二卷积层的输出作为输入，采用窗函数大小为2像素×2像素作最大池化抽样处理输出；
[0034] (3d)全连接层采用反向传播神经网络模型，以第二下采样层的输出作为输入，最后输出回波个数的估计结果。
[0035] 进一步的，在所述的步骤(3c)中，各个卷积层的输出计算公式为：
[0036]
[0037] 其中，yl为第l卷积层的输出， 为第l卷积层的第j个卷积核， 的计算公式为：
[0038]
[0039] randj(5)是产生第j个5×5卷积核矩阵的随机函数， 为卷积符号，sigmoid()为神经元的激活函数，
[0040] 进一步的，在所述的步骤(3c)中，各个下采样层的输出计算公式为：(3e)下采样层，按以下公式计算：
[0041] xl＝maxyl-1u(n,n)
[0042] 其中，u(n,n)为窗函数，且n＝2。
[0043] 进一步的，在所述的步骤(3b)中，将训练样本集中所有时频谱图进行预处理时生成41像素×50像素的灰度图像。
[0044] 与
现有技术相比，本发明具有的有益效果为：
[0045] 本发明的一种基于卷积神经网络的窄带超声回波个数的估计方法，先设定窄带超声回波信号模型参数，在不同噪声条件下生成多个回波信号集模型；在所生成的多个回波信号集模型随机提取一部分回波信号集模型进行短时傅里叶变换并生成时频谱图作为训练样本集，将剩下的回波信号集模型进行短时傅里叶变换并生成时频谱图作为测试样本集；将获得的训练样本集输入卷积神经网络进行训练，并输出训练好的卷积神经网络；将获得的测试样本集输入训练好的卷积神经网络，即可估计出对应回波信号的回波个数。本发明的一种基于卷积神经网络的窄带超声回波个数的估计方法能够在任何复杂环境下收集到的窄带超声回波信号都可以作为训练学习的样本，经过多个对比分析后，能够更精确的估算出窄带超声回波的个数，鲁棒性也更好。
附图说明
[0046] 图1是本发明的估算方法的处理
流程图；
[0047] 图2是本发明的估算方法中第三步的卷积神经网络结构图；
[0048] 图3是本发明的估算方法中第一步的随机生成的部分回波信号；
[0049] 图4是图3中的回波信号按本发明的估算方法进行短时傅里叶变换后生成的时频谱图；
[0050] 图5是本发明的估算方法中第三步中输入训练样本学习得到的
迭代次数与分类误差率关系图；
[0051] 图6是本发明的估算方法中第三步中8000个测试样本的真实值与估计值的误差图；
[0052] 图7是本发明的估算方法中第三步中8000个测试样本前200个样本的真实值与估计值的误差图。
具体实施方式
[0053] 下面结合具体实施方式，对本发明的技术方案作进一步的详细说明，但不构成对本发明的任何限制。
[0054] 参照图1所示，一种基于卷积神经网络的窄带超声回波个数的估计方法，包括以下步骤：
[0055] (1)设定窄带超声回波信号模型参数，在不同噪声条件下生成多个回波信号集模型。其中，具体包括以下步骤：
[0056] (1a)根据窄带超声回波的物理特性设定窄带超声回波信号集初始模型参数，窄带超声回波信号集初始模型参数的计算公式为：
[0057]
[0058] 其中，θ＝[α,τ,fc,φ,β]，
[0059] α为窄带超声回波的带宽，τ为窄带超声回波的到达时间，fc为窄带超声回波的中心频率，φ为窄带超声回波的相位，β为窄带超声回波的幅度系数，t为窄带超声回波的传递时间。
[0060] 随机设定的各参数取值为：α的取值范围为2-6MHz，τ的取值范围为0.5-10s，fc的取值范围为2-10MHz2，φ的取值范围为0-2rad，β的取值范围为0-2，M的取值范围为1-10个。
[0061] 再将窄带超声回波信号集初始模型参数中的窄带超声回波信号逐个做加性白噪声处理后得回波信号集模型参数，窄带超声回波信号做加性白噪声处理的计算公式为：
[0062]
[0063] 其中，w(t)为加性高斯白噪声，M是窄带超声回波信号集初始模型参数中的窄带超声回波个数。
[0064] (1b)设定噪声信噪比取值范围、步长值以及仿真次数，根据蒙特卡洛法，生成多个回波信号集模型。噪声信噪比的取值范围为-10dB-30dB，步长为1dB，则信噪比个数为41个，再设定仿真次数为100次，则总共生产的回波信号集模型为10×41×100＝41000个。如图3所示，图中给出了随机生成的部分回波信号。
[0065] (2)在步骤(1)所生成的多个回波信号集模型随机提取一部分回波信号集模型进行短时傅里叶变换并生成时频谱图作为训练样本集，一般提取所生成的多个回波信号集模型的80％的回波信号集模型，即将41000个回波信号集模型中的80％进行短时傅里叶变换并生成时频谱图作为训练样本集。将剩下的回波信号集模型进行短时傅里叶变换并生成时频谱图作为测试样本集。
[0066] 其中，具体包括以下步骤：
[0067] (2a)在步骤(1)所生成的多个回波信号模型中提取一部分回波信号集模型作短时傅里叶变换获取回波信号x(t)的时频谱图，作为训练样本集，所述的短时傅里叶变换公式为：
[0068]
[0069] 其中，t为窄带超声回波的传递时间，γ(t)为的分析窗函数,γ*(t)为γ(t)的复共轭函数，取STFT(τ,f)的模|STFT(τ,f)|为x(t)的时频谱图。
[0070] (2b)将步骤(1)所生产的多个回波信号模型中剩下的回波信号模型作短时傅里叶变换获取回波信号x(t)的时频谱图，作为测试样本集，所述的短时傅里叶变换公式为：
[0071]
[0072] 其中，t为窄带超声回波的传递时间，γ(t)为的分析窗函数,γ*(t)为γ(t)的复共轭函数，取STFT(τ,f)的模|STFT(τ,f)|为x(t)的时频谱图。如图4所示，图中根据图3所给出的部分回波信号对应生成的时频谱图。
[0073] (3)将步骤(2)中获得的训练样本集输入卷积神经网络进行训练，并输出训练好的卷积神经网络，卷积神经网络结构如图2所示。具体包括以下步骤：
[0074] (3a)所述的卷积神经网络依次设为输入层、隐层和全连接层三层结构。
[0075] (3b)将训练样本集中所有时频谱图进行预处理为的灰度图像，并对应生成回波个数标签集，将训练样本集中所有时频谱图作为标签输入层的输入，对应生成的回波个数标签集作为处理输出。其中，将训练样本集中所有时频谱图进行预处理时生成41像素×50像素的灰度图像。
[0076] (3c)所述的隐层依次设为第一卷积层、第一下采样层、第二卷积层以及第二下采样层这四层结构，所述的第一卷积层以灰度图像作为输入，采用4个5像素×5像素的卷积核进行处理输出，所述的第一下采样层以第一卷积层的输出作为输入，采用窗函数大小为2像素×2像素作最大池化抽样处理输出，所述的第二卷积层以第一下采样层的输出作为输入，采用8个卷积核进行处理输出，所述的第二下采样层以第二卷积层的输出作为输入，采用窗函数大小为2像素×2像素作最大池化抽样处理输出；
[0077] (3d)全连接层采用反向传播神经网络模型，以第二下采样层的输出作为输入，最后输出回波个数的估计结果。
[0078] 各个卷积层的输出计算公式为：
[0079]
[0080] 其中，yl为第l卷积层的输出， 为第l卷积层的第j个卷积核， 的计算公式为：
[0081]
[0082] randj(5)是产生第j个5×5卷积核矩阵的随机函数， 为卷积符号，sigmoid()为神经元的激活函数，
[0083] 各个下采样层的输出计算公式为：(3e)下采样层，按以下公式计算：
[0084] xl＝maxyl-1u(n,n)
[0085] 其中，u(n,n)为窗函数，且n＝2。
[0086] (4)将步骤(2)中获得的测试样本集输入步骤(3)中训练好的卷积神经网络，即可估计出对应回波信号的回波个数。
[0087] 图5为在本发明的一种基于卷积神经网络的窄带超声回波个数的估计方法中步骤(3)中训练卷积神经网络时，迭代次数与分类误差率的关系情况，图6为在步骤(3)中训练卷积神经网络时，所输入的8000个训练样本的真实值与估计值的误差情况，图7为在步骤(3)中训练卷积神经网络时，所输入的8000个训练样本集的前200个训练样本的真实值与估计值的误差情况，本仿真实验测试误差率为16.55％，即8000个测试样本中，正确估计出结果的有6676个，能够比较准确的估算出窄带超声回波个数。
[0088] 本发明的一种基于卷积神经网络的窄带超声回波个数的估计方法能够在任何复杂环境下收集到的窄带超声回波信号都可以作为训练学习的样本，经过多个对比分析后，能够更精确的估算出窄带超声回波的个数，鲁棒性也更好。
[0089] 以上所述仅为本发明的较佳
实施例，凡在本发明的精神和原则范围内所作的任何
修改、等同替换和改进等，均应包含在本发明的保护范围之内。
