CN112861731B - 一种基于参数寻优的超声信号去噪方法 - Google Patents

一种基于参数寻优的超声信号去噪方法 Download PDF

Info

Publication number
CN112861731B
CN112861731B CN202110181618.2A CN202110181618A CN112861731B CN 112861731 B CN112861731 B CN 112861731B CN 202110181618 A CN202110181618 A CN 202110181618A CN 112861731 B CN112861731 B CN 112861731B
Authority
CN
China
Prior art keywords
fruit fly
parameter
population
ith
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110181618.2A
Other languages
English (en)
Other versions
CN112861731A (zh
Inventor
董明
李敬
张广明
陈渊
万翔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian University of Science and Technology
Original Assignee
Xian University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian University of Science and Technology filed Critical Xian University of Science and Technology
Priority to CN202110181618.2A priority Critical patent/CN112861731B/zh
Publication of CN112861731A publication Critical patent/CN112861731A/zh
Application granted granted Critical
Publication of CN112861731B publication Critical patent/CN112861731B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02—Preprocessing
    • G06F2218/04—Denoising
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Data Mining & Analysis (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Signal Processing (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于参数寻优的超声信号去噪方法，包括步骤：一、超声回波信号的获取和存储；二、超声回波信号去噪：201、基于参数寻优的稀疏分解，202、寻找最优时频参数，203、迭代停止条件判断。本发明采用基于参数寻优的稀疏分解方法寻找最佳匹配原子，同时结合自适应迭代停止条件和自适应搜索步长，能有效提高去噪后信号的信噪比，加快去噪的速度，确保重构信号误差小。

Description

一种基于参数寻优的超声信号去噪方法
技术领域
本发明属于超声信号去噪技术领域，具体涉及一种基于参数寻优的超声信号去噪方法。
背景技术
在超声检测中，超声探头可以将电信号转换为超声信号，通过耦合剂传入被检测材料，当超声波在传播过程中遇到缺陷，会发生反射、折射等现象。通过对回波分析，可以估计缺陷的位置、大小。但是在实际检测中，超声波在材料内部传播时会与缺陷、界面、边界等相互作用，产生透射、散射等现象导致能量衰减、模式转换等问题，导致回波中包含大量噪声，很难将缺陷回波分离出来，影响对缺陷的准确评价。
对超声信号去噪常用的方法是将信号在特定的基函数上进行分解，例如小波变换和傅里叶变换等，基函数的特征决定了信号的特征。超声信号具有非平稳性，波形随时间发生变化，在单一的、有限数量的基函数上分解信号，得到的信号仅仅依赖于基函数的特征，信号自身的特征则被忽略，容易失真。稀疏分解是另一种被广泛应用的信号处理方法，方法提出了过完备字典的概念，将信号在字典中的原子上分解，基函数是过完备的，可以根据信号特征灵活选取，提高去噪效果。近几年提出的稀疏分解算法有LASSO、稀疏贝叶斯学习、贪婪基追踪、非凸正则化、匹配追踪等。其中匹配追踪算法算法被广泛用于超声信号去噪，但是其仍存在不足，这是因为匹配追踪算法在匹配原子时是贪婪搜索的，这会导致最初选择的原子不是最稀疏的，后续原子的选择都是在给最初的原子进行补偿。为此提出了正交匹配追踪，对每一次匹配到的最佳原子进行正交化化处理，则残差与选择过的原子都是正交的，同一个原子不会被选择两次。与匹配追踪算法相比，正交匹配追踪算法不会重复选择同一个原子，迭代次数少，但是其引入了正交化，每一次迭代的计算量较大，计算复杂度升高。
针对算法运行时间长的问题，可以将字典的有限离散参数空间变为无限连续参数空间，在连续空间进行稀疏分解，可以适应超声信号的非平稳性，加快匹配追踪的速度，提高重构回波精度。这可以通过群搜索算法实现，其中包括人工鱼群算法、人工蜂群算法、粒子群算法等。但是这些算法复杂度高，计算量大，运行时间长。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种基于参数寻优的超声信号去噪方法，采用基于参数寻优的稀疏分解方法寻找最佳匹配原子，同时结合自适应迭代停止条件和自适应搜索步长，能有效提高去噪后信号的信噪比，加快去噪的速度，确保重构信号误差小，便于推广使用。
为解决上述技术问题，本发明采用的技术方案是：一种基于参数寻优的超声信号去噪方法，其特征在于，该方法包括以下步骤：
步骤一、超声回波信号的获取和存储：采用信号采集设备获取超声回波信号f(t)，并将超声回波信号f(t)存储在计算机中，信号采集设备获取超声回波信号f(t)为待处理信号且f(t)＝[f(t1),f(t2),...,f(tn),...,f(tN)]，其中，n为采样时刻编号且n＝1,2,...,N，N为信号的长度，f(tn)为f(t)在tn时刻的采样点；
步骤二、超声回波信号去噪：采用计算机对超声回波信号f(t)进行去噪，过程如下：
步骤201、基于参数寻优的稀疏分解：调用正交匹配追踪算法模块对步骤一中的超声回波信号f(t)进行分解，即
Figure BDA0002941614260000021
得到K个最佳匹配原子，
Figure BDA0002941614260000022
为K个最佳匹配原子中第k个最佳匹配原子，k是正整数且k＝1，2，...，K，
Figure BDA0002941614260000023
为超声回波信号f(t)第k次迭代后的残差量，ak为第k次迭代产生的残差在对应的最佳匹配原子上的投影；
其中，
Figure BDA0002941614260000031
为Gabor原子且
Figure BDA0002941614260000032
为高斯窗函数且
Figure BDA0002941614260000033
rk为
Figure BDA0002941614260000034
的时频参数且rk＝(sk,uk,vk,wk)，sk为第k次迭代中的尺度参数，uk为第k次迭代中的平移参数，vk为第k次迭代中的频率参数，wk为第k次迭代中的相位参数；
步骤202、寻找最优时频参数：调用参数寻优算法模块寻找最优时频参数，在时频参数范围内找出最大味道浓度smellbest(rk)的时频参数rk，再根据
Figure BDA0002941614260000035
得出最佳匹配原子；
其中，
Figure BDA0002941614260000036
为第k-1次迭代所得残差和第k次迭代得到的最佳匹配原子的内积，未开始迭代时，初始残差
Figure BDA0002941614260000037
和超声回波信号f(t)的初始信号相等；
所述参数寻优算法模块为果蝇优化算法模块，参数寻优的过程如下：
步骤2021、参数初始化：采用计算机设定果蝇种群体的最大迭代次数maxgen和种群大小sizepop；
步骤2022、生成果蝇位置：采用计算机随机生成两个果蝇种群的初始位置，两个果蝇种群的初始位置均为四维向量，两个果蝇种群分别用X和Y表示，并赋予果蝇种群X的初始位置
Figure BDA0002941614260000038
果蝇种群Y的初始位置
Figure BDA0002941614260000039
为随机生成的果蝇种群X的初始位置时频参数的初始值，
Figure BDA00029416142600000310
为随机生成的果蝇种群Y的初始位置时频参数的初始值；
步骤2023、给定果蝇搜索的方向和步长：每个果蝇从初始位置离开，开始搜索最佳味道浓度蜜源的位置，调用自适应步长模块，根据公式
Figure BDA00029416142600000311
确定果蝇种群X中第i个果蝇此时的位置Xi和果蝇种群Y中第i个果蝇此时的位置Yi，i为果蝇编号且i＝1、2、3、...、sizepop，Lg为果蝇的搜索步长，g是果蝇优化算法当前迭代次数，g＝1、2、3、...、maxgen；
其中，
Figure BDA0002941614260000041
为果蝇种群X中第i个果蝇的尺度参数且
Figure BDA0002941614260000042
为果蝇种群X中第i个果蝇的平移参数且
Figure BDA0002941614260000043
为果蝇种群X中第i个果蝇的频率参数且
Figure BDA0002941614260000044
为果蝇种群X中第i个果蝇的相位参数且
Figure BDA0002941614260000045
为果蝇种群Y中第i个果蝇的尺度参数且
Figure BDA0002941614260000046
为果蝇种群Y中第i个果蝇的平移参数且
Figure BDA0002941614260000047
为果蝇种群Y中第i个果蝇的频率参数且
Figure BDA0002941614260000048
为果蝇种群Y中第i个果蝇的相位参数且
Figure BDA0002941614260000049
步骤2024、确定最佳蜜源的大小和位置：根据公式
Figure BDA00029416142600000410
计算果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中时频参数到原点的有效距离，其中，
Figure BDA00029416142600000411
为果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中尺度参数到原点的有效距离，
Figure BDA00029416142600000412
为果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中平移参数到原点的有效距离，
Figure BDA00029416142600000413
为果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中频率参数到原点的有效距离，
Figure BDA00029416142600000414
为果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中相位参数到原点的有效距离；
根据公式
Figure BDA00029416142600000415
计算果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇的尺度参数判别式
Figure BDA00029416142600000416
平移参数判别式
Figure BDA00029416142600000417
频率参数判别式
Figure BDA00029416142600000418
相位参数判别式
Figure BDA0002941614260000051
将果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇的尺度参数判别式
Figure BDA0002941614260000052
平移参数判别式
Figure BDA0002941614260000053
频率参数判别式
Figure BDA0002941614260000054
相位参数判别式
Figure BDA0002941614260000055
代入
Figure BDA0002941614260000056
并与第k-1次迭代所得残差
Figure BDA0002941614260000057
内积运算，此时smelli作为蜜源的味道浓度值并保存在蜜源位置味道浓度集合中；
步骤2025、获取果蝇的位置和味道浓度：在蜜源位置味道浓度集合中获取时频参数范围内最大味道浓度smellbest和果蝇的位置bestindex，群体中的其他果蝇朝着味道浓度最高的蜜源飞去；
步骤2026、多次循环2023至步骤2025，直至果蝇种群体的最大迭代次数maxgen迭代完毕或时频参数范围内最大味道浓度smellbest不变，获取最优时频参数和最佳匹配原子；
步骤203、迭代停止条件判断：经过k次迭代，此时得到重构信号fk(t)和重构信号fk(t)的信噪比SNRk，再进行第k+1次迭代得到重构信号fk+1(t)和重构信号fk+1(t)的信噪比SNRk+1，调用迭代停止条件模块判断是否可以停止迭代，若满足条件则停止迭代，输出超声回波信号f(t)去噪信号；若不满足条件则返回步骤201继续分解。
上述的一种基于参数寻优的超声信号去噪方法，其特征在于：所述时频参数范围包括尺度参数范围、平移参数范围、频率参数范围、相位参数范围，所述尺度参数范围为
Figure BDA0002941614260000058
平移参数范围为[1,N]，频率参数范围为
Figure BDA0002941614260000059
相位参数范围为[0,π]，其中，f0为超声回波信号f(t)的超声探头的中心频率。
上述的一种基于参数寻优的超声信号去噪方法，其特征在于：所述迭代停止条件为q＜ε，其中，q为误差且
Figure BDA00029416142600000510
ε为给定阈值且ε取0.01～0.1，＜fk+1(t),fk(t)＞为fk+1(t)和fk(t)的内积。
上述的一种基于参数寻优的超声信号去噪方法，其特征在于：所述果蝇的搜索步长
Figure BDA0002941614260000061
其中，L0为初始步长，T为步长收敛系数且
Figure BDA0002941614260000062
α为正交匹配追踪算法的迭代次数，β为步长收敛因子且β取0.1～0.25，γ为非均匀变异因子且γ取0.6，r为0-1之间的随机数。
本发明与现有技术相比具有以下优点：
1、本发明采用基于参数寻优的稀疏分解方法寻找最佳匹配原子，同时结合自适应迭代停止条件和自适应搜索步长，能有效提高去噪后信号的信噪比，信号去噪的速度大幅加快，相比于传统稀疏分解，该方法能在更短时间内完成信号去噪，实现信号实时去噪，便于推广使用。
2、本发明能够根据重构信号和其信噪比自适应停止迭代，该方法具有很强的鲁棒性，能够克服传统迭代停止条件无法选择迭代停止阈值的问题，可靠稳定，使用效果好。
3、本发明方法步骤简单，改进了果蝇优化算法的搜索步长，提出的自适应步长相比于随机步长，在搜索过程中搜索精度更高，不容易漏掉味道浓度高的蜜源，克服了随机步长容易陷入局部最优的问题，便于推广使用。
综上所述，本发明采用基于参数寻优的稀疏分解方法寻找最佳匹配原子，同时结合自适应迭代停止条件和自适应搜索步长，能有效提高去噪后信号的信噪比，加快去噪的速度，确保重构信号误差小，便于推广使用。
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明的方法流程框图。
图2为本发明采集含缺陷试块的缺陷和底面回波所用***的连接关系示意图。
图3为本发明计算机接收到的含缺陷试块的缺陷和底面回波时域波形图。
图4为本发明计算机接收到的含缺陷试块的缺陷和底面回波频谱图。
图5为本发明对原始波形处理后的信号的时域波形图。
图6为本发明对原始波形处理后的信号的频谱图。
附图标记说明:
1—信号采集设备；     2—计算机；     3—含缺陷试块；
4—超声探头。
具体实施方式
如图1至图6所示，本发明的一种基于参数寻优的超声信号去噪方法，包括以下步骤：
步骤一、超声回波信号的获取和存储：采用信号采集设备1获取超声回波信号f(t)，并将超声回波信号f(t)存储在计算机2中，信号采集设备1获取超声回波信号f(t)为待处理信号且f(t)＝[f(t1),f(t2),...,f(tn),...,f(tN)]，其中，n为采样时刻编号且n＝1,2,...,N，N为信号的长度，f(tn)为f(t)在tn时刻的采样点；
实际实施时，制作标准的含缺陷试块3，标准的含缺陷试块3的形状为长方体，其长和宽均为5D，标准的含缺陷试块3的高度范围为
Figure BDA0002941614260000071
λ为超声波在含缺陷试块3中的波长，当超声探头4为直探头时，D是直探头的直径；当超声探头4是斜探头时，D是斜探头的较长边的边长，缺陷是一个直径2mm的平底孔，孔的高度为30mm。
实际使用时，超声探头4是型号为2.5Z20N的圆形直探头，直探头的直径D为20mm，中心频率是2.5MHz；标准的含缺陷试块的材料是45钢，超声波在标准的含缺陷试块中的传播速度为5920m/s，标准是含缺陷试块的高度为65mm，长和宽均为100mm。
如图2所示，将超声探头3放在含缺陷试块的顶部中心处，用信号采集设备1获取缺陷和底面的超声回波信号f(t)，并将超声回波信号f(t)存储在计算机2中。
本实施例中，信号采集设备1优选的采用CTS-4020超声探伤仪。存储在计算机2中的缺陷和底面回波如图3所示，其频率如图4所示。
步骤二、超声回波信号去噪：采用计算机2对超声回波信号f(t)进行去噪，过程如下：
步骤201、基于参数寻优的稀疏分解：调用正交匹配追踪算法模块对步骤一中的超声回波信号f(t)进行分解，即
Figure BDA0002941614260000081
得到K个最佳匹配原子，
Figure BDA0002941614260000082
为K个最佳匹配原子中第k个最佳匹配原子，k是正整数且k＝1，2，...，K，
Figure BDA0002941614260000083
为超声回波信号f(t)第k次迭代后的残差量，ak为第k次迭代产生的残差在对应的最佳匹配原子上的投影；
其中，
Figure BDA0002941614260000084
为Gabor原子且
Figure BDA0002941614260000085
为高斯窗函数且
Figure BDA0002941614260000086
rk为
Figure BDA0002941614260000087
的时频参数且rk＝(sk,uk,vk,wk)，sk为第k次迭代中的尺度参数，uk为第k次迭代中的平移参数，vk为第k次迭代中的频率参数，wk为第k次迭代中的相位参数；
步骤202、寻找最优时频参数：调用参数寻优算法模块寻找最优时频参数，在时频参数范围内找出最大味道浓度smellbest(rk)的时频参数rk，再根据
Figure BDA0002941614260000088
得出最佳匹配原子；
其中，
Figure BDA0002941614260000089
为第k-1次迭代所得残差和第k次迭代得到的最佳匹配原子的内积，未开始迭代时，初始残差
Figure BDA00029416142600000810
和超声回波信号f(t)的初始信号相等；
本实施例中，所述时频参数范围包括尺度参数范围、平移参数范围、频率参数范围、相位参数范围，所述尺度参数范围为
Figure BDA00029416142600000811
平移参数范围为[1,N]，频率参数范围为
Figure BDA00029416142600000812
相位参数范围为[0,π]，其中，f0为超声回波信号f(t)的超声探头的中心频率。
所述参数寻优算法模块为果蝇优化算法模块，参数寻优的过程如下：
步骤2021、参数初始化：采用计算机2设定果蝇种群体的最大迭代次数maxgen和种群大小sizepop；
本实施例中，采用越界回折的方法保证果蝇在参数范围内寻优，对果蝇搜索到的新蜜源的进行边界检测，若超出了参数范围的上下边界，则对新蜜源的四个元素均进行越界回折操作，越界回折后获得蜜源作为新蜜源，保证搜索的效率和精度。
步骤2022、生成果蝇位置：采用计算机2随机生成两个果蝇种群的初始位置，两个果蝇种群的初始位置均为四维向量，两个果蝇种群分别用X和Y表示，并赋予果蝇种群X的初始位置
Figure BDA0002941614260000091
果蝇种群Y的初始位置
Figure BDA0002941614260000092
为随机生成的果蝇种群X的初始位置时频参数的初始值，
Figure BDA0002941614260000093
为随机生成的果蝇种群Y的初始位置时频参数的初始值；
步骤2023、给定果蝇搜索的方向和步长：每个果蝇从初始位置离开，开始搜索最佳味道浓度蜜源的位置，调用自适应步长模块，根据公式
Figure BDA0002941614260000094
确定果蝇种群X中第i个果蝇此时的位置Xi和果蝇种群Y中第i个果蝇此时的位置Yi，i为果蝇编号且i＝1、2、3、...、sizepop，Lg为果蝇的搜索步长，g是果蝇优化算法当前迭代次数，g＝1、2、3、...、maxgen；
其中，
Figure BDA0002941614260000095
为果蝇种群X中第i个果蝇的尺度参数且
Figure BDA0002941614260000096
为果蝇种群X中第i个果蝇的平移参数且
Figure BDA0002941614260000097
为果蝇种群X中第i个果蝇的频率参数且
Figure BDA0002941614260000098
为果蝇种群X中第i个果蝇的相位参数且
Figure BDA0002941614260000099
为果蝇种群Y中第i个果蝇的尺度参数且
Figure BDA00029416142600000910
为果蝇种群Y中第i个果蝇的平移参数且
Figure BDA00029416142600000911
为果蝇种群Y中第i个果蝇的频率参数且
Figure BDA00029416142600000912
为果蝇种群Y中第i个果蝇的相位参数且
Figure BDA00029416142600000913
本实施例中，所述果蝇的搜索步长
Figure BDA00029416142600000914
其中，L0为初始步长，T为步长收敛系数且
Figure BDA00029416142600000915
α为正交匹配追踪算法的迭代次数，β为步长收敛因子且β取0.1～0.25，γ为非均匀变异因子且γ取0.6，r为0-1之间的随机数。
本实施例中，随机数r是由随机数rand产生的。
步骤2024、确定最佳蜜源的大小和位置：根据公式
Figure BDA0002941614260000101
计算果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中时频参数到原点的有效距离，其中，
Figure BDA0002941614260000102
为果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中尺度参数到原点的有效距离，
Figure BDA0002941614260000103
为果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中平移参数到原点的有效距离，
Figure BDA0002941614260000104
为果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中频率参数到原点的有效距离，
Figure BDA0002941614260000105
为果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中相位参数到原点的有效距离；
根据公式
Figure BDA0002941614260000106
计算果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇的尺度参数判别式
Figure BDA0002941614260000107
平移参数判别式
Figure BDA0002941614260000108
频率参数判别式
Figure BDA0002941614260000109
相位参数判别式
Figure BDA00029416142600001010
将果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇的尺度参数判别式
Figure BDA00029416142600001011
平移参数判别式
Figure BDA00029416142600001012
频率参数判别式
Figure BDA00029416142600001013
相位参数判别式
Figure BDA00029416142600001014
代入
Figure BDA00029416142600001015
并与第k-1次迭代所得残差
Figure BDA00029416142600001016
内积运算，此时smelli作为蜜源的味道浓度值并保存在蜜源位置味道浓度集合中；
步骤2025、获取果蝇的位置和味道浓度：在蜜源位置味道浓度集合中获取时频参数范围内最大味道浓度smellbest和果蝇的位置bestindex，群体中的其他果蝇朝着味道浓度最高的蜜源飞去；
步骤2026、多次循环2023至步骤2025，直至果蝇种群体的最大迭代次数maxgen迭代完毕或时频参数范围内最大味道浓度smellbest不变，获取最优时频参数和最佳匹配原子；
步骤203、迭代停止条件判断：经过k次迭代，此时得到重构信号fk(t)和重构信号fk(t)的信噪比SNRk，再进行第k+1次迭代得到重构信号fk+1(t)和重构信号fk+1(t)的信噪比SNRk+1，调用迭代停止条件模块判断是否可以停止迭代，若满足条件则停止迭代，输出超声回波信号f(t)去噪信号；若不满足条件则返回步骤201继续分解。
本实施例中，所述迭代停止条件为q＜ε，其中，q为误差且
Figure BDA0002941614260000111
ε为给定阈值且ε取0.01～0.1，＜fk+1(t),fk(t)＞为fk+1(t)和fk(t)的内积。
本实施例中，给定阈值ε的大小决定了超声信号稀疏分解的程度，若选取的过大，则会导致信号过匹配，若选取的过小，则会导致信号中的噪声还没有完全去除就停止迭代。
本实施例中，超声信号去噪后的时频如图5所示，频率如图6所示；
从时域波形上来看，去噪后的信号和原始信号相似度较高。图4是原始信号的频谱图，由于频率混叠等因素的影响，在原始信号中产生了杂乱的低频信号，影响对正常信号的判断，去噪后的信号低频部分明显被抑制，高频部分与原始信号基本重合。
本发明使用时，采用基于参数寻优的稀疏分解方法寻找最佳匹配原子，同时结合自适应迭代停止条件和自适应搜索步长，能有效提高去噪后信号的信噪比，信号去噪的速度大幅加快，相比于传统稀疏分解，该方法能在更短时间内完成信号去噪，实现信号实时去噪；能够根据重构信号和其信噪比自适应停止迭代，该方法具有很强的鲁棒性，能够克服传统迭代停止条件无法选择迭代停止阈值的问题；改进了果蝇优化算法的搜索步长，提出的自适应步长相比于随机步长，在搜索过程中搜索精度更高，不容易漏掉味道浓度高的蜜源，克服了随机步长容易陷入局部最优的问题。
以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化，均仍属于本发明技术方案的保护范围内。

Claims (2)

1.一种基于参数寻优的超声信号去噪方法，其特征在于，该方法包括以下步骤：
步骤一、超声回波信号的获取和存储：采用信号采集设备(1)获取超声回波信号f(t)，并将超声回波信号f(t)存储在计算机(2)中，信号采集设备(1)获取超声回波信号f(t)为待处理信号且f(t)＝[f(t1),f(t2),...,f(tn),...,f(tN)]，其中，n为采样时刻编号且n＝1,2,...,N，N为信号的长度，f(tn)为f(t)在tn时刻的采样点；
步骤二、超声回波信号去噪：采用计算机(2)对超声回波信号f(t)进行去噪，过程如下：
步骤201、基于参数寻优的稀疏分解：调用正交匹配追踪算法模块对步骤一中的超声回波信号f(t)进行分解，即
Figure FDA0004019956480000011
得到K个最佳匹配原子，
Figure FDA0004019956480000012
为K个最佳匹配原子中第k个最佳匹配原子，k是正整数且k＝1，2，...，K，
Figure FDA0004019956480000013
为超声回波信号f(t)第k次迭代后的残差量，ak为第k次迭代产生的残差在对应的最佳匹配原子上的投影；
其中，
Figure FDA0004019956480000014
为Gabor原子且
Figure FDA0004019956480000015
Figure FDA0004019956480000016
为高斯窗函数且
Figure FDA0004019956480000017
rk为
Figure FDA0004019956480000018
的时频参数且rk＝(sk,uk,vk,wk)，sk为第k次迭代中的尺度参数，uk为第k次迭代中的平移参数，vk为第k次迭代中的频率参数，wk为第k次迭代中的相位参数；
步骤202、寻找最优时频参数：调用参数寻优算法模块寻找最优时频参数，在时频参数范围内找出最大味道浓度smellbest(rk)的时频参数rk，再根据
Figure FDA0004019956480000019
得出最佳匹配原子；
其中，
Figure FDA00040199564800000110
Figure FDA00040199564800000111
为第k-1次迭代所得残差和第k次迭代得到的最佳匹配原子的内积，未开始迭代时，初始残差
Figure FDA00040199564800000112
和超声回波信号f(t)的初始信号相等；
所述参数寻优算法模块为果蝇优化算法模块，参数寻优的过程如下：
步骤2021、参数初始化：采用计算机(2)设定果蝇种群体的最大迭代次数maxgen和种群大小sizepop；
步骤2022、生成果蝇位置：采用计算机(2)随机生成两个果蝇种群的初始位置，两个果蝇种群的初始位置均为四维向量，两个果蝇种群分别用X和Y表示，并赋予果蝇种群X的初始位置
Figure FDA0004019956480000021
果蝇种群Y的初始位置
Figure FDA0004019956480000022
Figure FDA0004019956480000023
为随机生成的果蝇种群X的初始位置时频参数的初始值，
Figure FDA0004019956480000024
为随机生成的果蝇种群Y的初始位置时频参数的初始值；
步骤2023、给定果蝇搜索的方向和步长：每个果蝇从初始位置离开，开始搜索最佳味道浓度蜜源的位置，调用自适应步长模块，根据公式
Figure FDA0004019956480000025
确定果蝇种群X中第i个果蝇此时的位置Xi和果蝇种群Y中第i个果蝇此时的位置Yi，i为果蝇编号且i＝1、2、3、...、sizepop，Lg为果蝇的搜索步长，g是果蝇优化算法当前迭代次数，g＝1、2、3、...、maxgen；
其中，
Figure FDA0004019956480000026
Figure FDA0004019956480000027
为果蝇种群X中第i个果蝇的尺度参数且
Figure FDA0004019956480000028
Figure FDA0004019956480000029
为果蝇种群X中第i个果蝇的平移参数且
Figure FDA00040199564800000210
Figure FDA00040199564800000211
为果蝇种群X中第i个果蝇的频率参数且
Figure FDA00040199564800000212
Figure FDA00040199564800000213
为果蝇种群X中第i个果蝇的相位参数且
Figure FDA00040199564800000214
Figure FDA00040199564800000215
为果蝇种群Y中第i个果蝇的尺度参数且
Figure FDA00040199564800000216
Figure FDA00040199564800000217
为果蝇种群Y中第i个果蝇的平移参数且
Figure FDA00040199564800000218
Figure FDA00040199564800000219
为果蝇种群Y中第i个果蝇的频率参数且
Figure FDA00040199564800000220
Figure FDA00040199564800000221
为果蝇种群Y中第i个果蝇的相位参数且
Figure FDA00040199564800000222
步骤2024、确定最佳蜜源的大小和位置：根据公式
Figure FDA00040199564800000223
计算果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中时频参数到原点的有效距离，其中，
Figure FDA0004019956480000031
为果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中尺度参数到原点的有效距离，
Figure FDA0004019956480000032
为果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中平移参数到原点的有效距离，
Figure FDA0004019956480000033
为果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中频率参数到原点的有效距离，
Figure FDA0004019956480000034
为果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇搜索到的蜜源中相位参数到原点的有效距离；
根据公式
Figure FDA0004019956480000035
计算果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇的尺度参数判别式
Figure FDA0004019956480000036
平移参数判别式
Figure FDA0004019956480000037
频率参数判别式
Figure FDA0004019956480000038
相位参数判别式
Figure FDA0004019956480000039
将果蝇种群X中第i个果蝇和果蝇种群Y中第i个果蝇的尺度参数判别式
Figure FDA00040199564800000310
平移参数判别式
Figure FDA00040199564800000311
频率参数判别式
Figure FDA00040199564800000312
相位参数判别式
Figure FDA00040199564800000313
代入
Figure FDA00040199564800000314
并与第k-1次迭代所得残差
Figure FDA00040199564800000315
内积运算，此时smelli作为蜜源的味道浓度值并保存在蜜源位置味道浓度集合中；
步骤2025、获取果蝇的位置和味道浓度：在蜜源位置味道浓度集合中获取时频参数范围内最大味道浓度smellbest和果蝇的位置bestindex，群体中的其他果蝇朝着味道浓度最高的蜜源飞去；
步骤2026、多次循环2023至步骤2025，直至果蝇种群体的最大迭代次数maxgen迭代完毕或时频参数范围内最大味道浓度smellbest不变，获取最优时频参数和最佳匹配原子；
步骤203、迭代停止条件判断：经过k次迭代，此时得到重构信号fk(t)和重构信号fk(t)的信噪比SNRk，再进行第k+1次迭代得到重构信号fk+1(t)和重构信号fk+1(t)的信噪比SNRk+1，调用迭代停止条件模块判断是否可以停止迭代，若满足条件则停止迭代，输出超声回波信号f(t)去噪信号；若不满足条件则返回步骤201继续分解；
所述迭代停止条件为q＜ε，其中，q为误差且
Figure FDA0004019956480000041
ε为给定阈值且ε取0.01～0.1，＜fk+1(t),fk(t)＞为fk+1(t)和fk(t)的内积；
所述果蝇的搜索步长
Figure FDA0004019956480000042
其中，L0为初始步长，T为步长收敛系数且
Figure FDA0004019956480000043
α为正交匹配追踪算法的迭代次数，β为步长收敛因子且β取0.1～0.25，γ为非均匀变异因子且γ取0.6，r为0-1之间的随机数。
2.按照权利要求1所述的一种基于参数寻优的超声信号去噪方法，其特征在于：所述时频参数范围包括尺度参数范围、平移参数范围、频率参数范围、相位参数范围，所述尺度参数范围为
Figure FDA0004019956480000044
平移参数范围为[1,N]，频率参数范围为
Figure FDA0004019956480000045
相位参数范围为[0,π]，其中，f0为超声回波信号f(t)的超声探头的中心频率。
CN202110181618.2A 2021-02-09 2021-02-09 一种基于参数寻优的超声信号去噪方法 Active CN112861731B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110181618.2A CN112861731B (zh) 2021-02-09 2021-02-09 一种基于参数寻优的超声信号去噪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110181618.2A CN112861731B (zh) 2021-02-09 2021-02-09 一种基于参数寻优的超声信号去噪方法

Publications (2)

Publication Number Publication Date
CN112861731A CN112861731A (zh) 2021-05-28
CN112861731B true CN112861731B (zh) 2023-04-18

Family

ID=75989596

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110181618.2A Active CN112861731B (zh) 2021-02-09 2021-02-09 一种基于参数寻优的超声信号去噪方法

Country Status (1)

Country Link
CN (1) CN112861731B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114235970B (zh) * 2021-12-20 2024-04-23 西安科技大学 一种自适应超声重叠回波分离方法
CN115236206B (zh) * 2022-06-30 2023-06-20 江南大学 基于多路径匹配追踪算法改进的高频超声去噪方法及***

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110208869A (zh) * 2019-07-08 2019-09-06 湖南师范大学 一种基于稀疏分解阈值设定的大地电磁信号去噪方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2993362B1 (fr) * 2012-07-12 2016-07-01 Commissariat Energie Atomique Procede de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage a ultrasons correspondants
CN106951956A (zh) * 2017-03-14 2017-07-14 中国矿业大学 一种基于果蝇优化算法的多机器人气味源定位方法
CN107070254B (zh) * 2017-04-13 2019-03-19 湖南科技大学 一种Buck-Boost矩阵变换器参数优化方法及装置
CN109507292B (zh) * 2018-12-26 2021-08-06 西安科技大学 一种信号提取方法
CN109632973B (zh) * 2018-12-26 2021-07-13 西安科技大学 一种基于多尺度匹配追踪的超声回波信号提取方法
CN110333054B (zh) * 2019-06-13 2021-07-06 电子科技大学 一种针对白车身焊接设备的缓变微小故障检测方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110208869A (zh) * 2019-07-08 2019-09-06 湖南师范大学 一种基于稀疏分解阈值设定的大地电磁信号去噪方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Application of basis pursuit in signal denoising of ultrasonic testing flaw;Ailing Qi等;《2017 First International Conference on Electronics Instrumentation & Information Systems (EIIS)》;第1-5页 *

Also Published As

Publication number Publication date
CN112861731A (zh) 2021-05-28

Similar Documents

Publication Publication Date Title
CN112861731B (zh) 一种基于参数寻优的超声信号去噪方法
CN108872962B (zh) 基于分数阶傅里叶变换的激光雷达微弱信号提取和分解方法
CN110926386B (zh) 基于改进小波阈值去噪的变压器绕组超声检测成像方法
CN109871733B (zh) 一种自适应海杂波信号去噪方法
CN107037136B (zh) 一种管道换热污垢阵列式超声检测方法
CN113449683B (zh) 基于k-svd训练局部字典的高频超声稀疏去噪方法及***
CN110133643B (zh) 植物根系探测方法及装置
CN110554049B (zh) 基于太赫兹波的复合绝缘子缺陷检测装置、方法及介质
CN112255607B (zh) 一种海杂波的抑制方法
CN109632973A (zh) 一种基于多尺度匹配追踪的超声回波信号提取方法
CN115236206B (zh) 基于多路径匹配追踪算法改进的高频超声去噪方法及***
CN112858183A (zh) 一种基于波形分离的增材制造激光超声信号缺陷成像方法
CN114117912A (zh) 一种数据模型双驱动下的海杂波建模与抑制方法
CN110515063A (zh) 基于迭代平稳离散小波变换的水声信号处理方法和装置
CN109410178A (zh) 一种工件裂纹检测方法及***
Bharadwaj et al. Target identification with wave-based matched pursuits and hidden Markov models
CN109584256A (zh) 一种基于霍夫直线检测的脉冲星dm估计算法
CN109507292A (zh) 一种信号提取方法
CN109682892A (zh) 一种基于时频分析的信号去燥方法
Shpigler et al. Detection of overlapping ultrasonic echoes with deep neural networks
CN112784811A (zh) 一种全局最优超声信号去噪方法
CN116203634A (zh) 一种基于低秩约束的鬼波去除方法
Tran-Duc et al. Application of l 1 Regularization for High-Quality Reconstruction of Ultrasound Tomography
Tsui et al. An adaptive threshold filter for ultrasound signal rejection
CN113030861B (zh) 一种时间反转和分数阶小波变换的水声信号检测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant