一种超声信号稀疏分解方法及其信号降噪与缺陷检测方法技术

本发明专利技术公开了一种超声信号稀疏分解方法及其降噪与缺陷检测方法，涉及超声无损检测领域，步骤一、采用Gabor模型对超声信号进行近似拟合，选取Gabor字典参数集合，组合形成字典矩阵，利用稀疏贝叶斯学习算法求解系数向量在稀疏约束条件下的最优近似解。利用原子向量参数来排除系数向量中对应于噪声的项并通过设置能量阈值进一步去噪，缺陷的数量即为系数向量中非零项的个数。本发明专利技术可以实现低信噪比的、包含复杂噪声的超声无损检测信号的降噪、缺陷反射信号的恢复和缺陷检测与定位。

【技术实现步骤摘要】
一种超声信号稀疏分解方法及其信号降噪与缺陷检测方法
本专利技术涉及超声无损检测领域，一种超声无损检测信号稀疏分解方法和降噪与缺陷检测方法。
技术介绍
在土木工程领域，超声无损检测在建筑结构、桥梁、管道等结构损伤检测中已经有着广泛的应用。借助于超声无损检测技术，可以有效检测结构中是否存在缺陷或损伤、进行缺陷定位并测量缺陷大小，以便对结构的安全状况进行及时评估。与结构健康监测的整体监测和少量局部监测技术相比，超声无损检测技术可以更好地检测结构局部微小损伤，是对结构健康监测技术的有益补充。将无损检测技术与结构健康监测技术相融合于一体，可以实现结构局部损伤至整体损伤的全尺度监测，将有效提高结构安全评估的准确性。土木工程结构复杂，构件材料与服役环境具备多样性，导致超声无损检测信号非常容易受到噪声等因素干扰，使得信号的有效解读变得异常困难，难以准确地获取损伤大小和部位，限制了超声无损检测技术在结构检测方面的有效性和实用性。近些年来，基于信号稀疏表示概念的数字信号处理技术已经成为超声信号处理方面的研究热点。利用冗余表示和稀疏性作为驱动力的信号处理方法在过去的数年中已经在图像去噪研究的一大焦点。在信号稀疏表示框架下，信号在某一过完备字典下可以直接进行稀疏分解。尽管此时信号的分解过程是欠静定的，过去的研究显示这种处理方法在某些应用中具有独特的优势，比如信号恢复和数据压缩。一方面，在捕捉数据特征方面，信号稀疏表示方法更具灵活性，对于特定一种信号，只需要从过完备字典中选取较少数量的原子来进行表征。由于每一个原子都描述这一种数据特征，因此，这样的表征结果会非常紧致。另外，信号稀疏表示相比常规方法具有超分辨率特性，同时还提升了信号存在微小扰动时表示结果的稳定性。目前，信号稀疏表示在超声信号处理方面已经得到较为广泛的研究，例如超声无损检测信号去噪和缺陷反射信号的恢复等，显示出优于诸如小波降噪等常规方法的性能。而目前的信号稀疏表示在超声无损检测信号处理方面的研究及应用，在噪声特征、字典设计和信号稀疏分解算法方面均存在着若干不足之处，使得信号稀疏表示的优势没有得到充分发掘，方法的抗噪声性能也有进一步提升的空间。
技术实现思路
本专利技术的目的是为了解决现有信号稀疏表示在超声无损检测信号处理中存在的字典设计方法和信号稀疏分解算法不够完善等问题，进而导致超声信号分析得到的结果不够精确，而提出的一种基于稀疏贝叶斯学习(SBL)的信号稀疏分解、恢复以及缺陷检测方法，适用于接触式脉冲反射法超声检测。首先，本专利技术提供一种超声信号稀疏分解方法，包括以下步骤：步骤一、采用Gabor模型分别对超声换能器发射的信号和构件底面预定最大探测深度处反射脉冲信号进行近似拟合；步骤二、选取Gabor字典的时间参数集合；步骤三、选取Gabor字典的频率参数集合；步骤四、选取Gabor字典的带宽参数集合；步骤五、基于步骤二-四选取的Gabor模型参数集合，组合形成字典矩阵并对原子向量(字典矩阵中的列向量称为原子向量)，同时对原子向量做能量归一化；步骤六、利用步骤五处理后的字典矩阵，依据信号稀疏表示模型，利用稀疏贝叶斯学习算法求解系数向量在稀疏约束条件下的最优近似解。2.根据权利要求1所述的方法，其特征在于：所述步骤一具体为：采用Gabor模型分别对超声换能器发射的信号和构件底面预定最大探测深度处反射脉冲信号进行近似拟合，获得超声换能器激发信号对应的Gabor模型频率参数f0、带宽参数s0，和试件底面或预定最大探测深度处反射脉冲信号的Gabor模型频率参数fb、带宽参数sb，Gabor模型表达式为：g(t)＝C0exp[-(t-u)2/(2s2)]exp[i2πf(t-u)]，其中，u表示Gabor脉冲信号的中心在时间轴上的位置、f表示脉冲信号中心频率、s表示脉冲信号带宽、t表示时间轴上的时间序列、i为虚数、C0是用于信号能量归一化的参数；其中u、f、s为基本参数。所述步骤二具体为：选取Gabor字典的时间参数集合：对于包含N个采样点的待分析信号y，参数ui的选取方法如下：ui＝aT0Δu,a＝0,1,…,(N-1)/Δu((N-1)/Δu是整数)；T0为信号的采样时间间隔，取值为采样频率的倒数，Δu为相邻时间参数之间的间隔(指相邻的两个u之间做差)，取值为T0的整数倍；参数ui的个数Nu满足N/4≤Nu≤N/2。所述步骤三具体为：选取Gabor字典的频率参数集合：对于包含N个采样点的待分析信号y，频率参数fi的选取方法如下：用竖线将信号y的频谱与频率轴之间的面积均分为Nf+1部分，则Nf条竖线与频率轴的交点对应的频率值即为fi的取值，参数fi的个数Nf满足5≤Nf≤15。所述步骤四具体为：选取Gabor字典的带宽参数集合：频率参数si在(0.5s0,1.5sb)范围内服从均匀分布，参数si的个数Nf满足5≤Nf≤10。所述步骤五具体为：基于步骤二-四选取的Gabor模型参数集合，组合形成字典矩阵Φ并对字典中的原子向量做能量归一化，能量归一化参数与步骤一相同，字典Φ的维数为N×(Nu×Nf×Ns)。所述步骤六具体为：利用步骤五处理后的字典矩阵Φ，依据信号稀疏表示模型：y＝Φc+ε，其中，ε为稀疏模型中的残余项，包含信号中的噪声和稀疏模型的误差；利用稀疏贝叶斯学习算法求解系数向量c在稀疏约束条件下的最优近似解：其中，γ为约束系数向量c的参数，通过稀疏贝叶斯学习算法可以计算得到。本专利技术还提供一种基于上述信号稀疏分解方法的信号降噪与缺陷检测方法，包括以下步骤：步骤一、依据权利要求1-7任一项所述方法得到的系数向量提取所有非零系数对应原子向量的三个参数(uj,fj,sj)j＝1,…L，L为系数向量中非零项的个数；步骤二、利用原子向量参数来排除系数向量中对应于噪声的项：步骤三、通过设置能量阈值进一步去噪：步骤四、经过去噪得到的信号：步骤五、缺陷检测：缺陷的数量Nd即为系数向量中非零项的个数；缺陷信号的时间延迟即为系数向量中非零项对应原子向量的时间参数各缺陷距离构件表面的距离为其中，v为超声波在材料中的纵波波速。本专利技术的有益效果为：1、由于噪声与缺陷信号的频谱特征比较相似而且有相当部分的重叠，以往方法在去噪过程中不可避免的也丢失了部分缺陷信号成分，而本方法可以实现无损去噪；2、本方法在设计字典矩阵时充分考虑了信号的能量分布，相比之前的方法在信号表示效率上更具有优势；3、本方法采用基于稀疏贝叶斯学习的稀疏优化求解算法，在处理信号稀疏表示这类病态问题方面具有独特的优势，对噪声的鲁棒性也更好。附图说明图1为本专利技术所涉及的接触式超声无损检测系统示意图；图2为本专利技术的信号截取示意图：(A).换能器接收的全部信号；(B).截取用于分析的信号y；图3为本专利技术中超声换能器发射信号中心频率f0和带宽参数s0的估测方法示意图；图4为本专利技术中测试构件底面反射信号中心频率fb和带宽参数sb的估测方法示意图；图5为本专利技术中的频率参数选取示意图；图6为本专利技术实施例的降噪与缺陷检测方法的关键步骤结果展示。具体实施方式实施例1一种超声信号稀疏分解方法和降噪与缺陷检测方法：本实施方式基于金属构件内部缺陷的接触式超声检测方法，如图1所示，包括：将超声无损检测采集到的超声信号进行降采样处理，使得处理之后的信号有效采本文档来自技高网...

【技术保护点】
1.一种超声信号稀疏分解方法，其特征在于：包括以下步骤：步骤一、采用Gabor模型分别对超声换能器发射的信号和构件底面预定最大探测深度处反射脉冲信号进行近似拟合；步骤二、选取Gabor字典的时间参数集合；步骤三、选取Gabor字典的频率参数集合；步骤四、选取Gabor字典的带宽参数集合；步骤五、基于步骤二‑四选取的Gabor模型参数集合，组合形成字典矩阵并对字典矩阵的列向量即原子向量做能量归一化；步骤六、利用步骤五处理后的字典矩阵，依据信号稀疏表示模型，利用稀疏贝叶斯学习算法求解系数向量在稀疏约束条件下的最优近似解。

【技术特征摘要】
1.一种超声信号稀疏分解方法，其特征在于：包括以下步骤：步骤一、采用Gabor模型分别对超声换能器发射的信号和构件底面预定最大探测深度处反射脉冲信号进行近似拟合；步骤二、选取Gabor字典的时间参数集合；步骤三、选取Gabor字典的频率参数集合；步骤四、选取Gabor字典的带宽参数集合；步骤五、基于步骤二-四选取的Gabor模型参数集合，组合形成字典矩阵并对字典矩阵的列向量即原子向量做能量归一化；步骤六、利用步骤五处理后的字典矩阵，依据信号稀疏表示模型，利用稀疏贝叶斯学习算法求解系数向量在稀疏约束条件下的最优近似解。2.根据权利要求1所述的方法，其特征在于：所述步骤一具体为：采用Gabor模型分别对超声换能器发射的信号和构件底面预定最大探测深度处反射脉冲信号进行近似拟合，获得超声换能器激发信号对应的Gabor模型频率参数f0、带宽参数s0，和试件底面或预定最大探测深度处反射脉冲信号的Gabor模型频率参数fb、带宽参数sb，Gabor模型表达式为：g(t)＝C0exp[-(t-u)2/(2s2)]exp[i2πf(t-u)]，其中，u表示Gabor脉冲信号在时间轴上的位置、f表示脉冲信号频率、s表示脉冲信号带宽、t表示时间轴上的时间序列、i为虚数、C0是用于信号能量归一化的参数；其中u、f、s为基本参数。3.根据权利要求2所述的方法，其特征在于：所述步骤二具体为：选取Gabor字典的时间参数集合：对于包含N个采样点的待分析信号y，参数ui的选取方法如下：ui＝aT0Δu,a＝0,1,…,(N-1)/Δu；T0为信号的采样时间间隔，取值为采样频率的倒数，Δu为相邻时间参数之间的间隔，取值为T0的整数倍；参数ui的个数Nu满足N/4≤Nu≤N/2。4.根据权利要求3所述的方法，其特征在于：所述步骤三...

【专利技术属性】
技术研发人员：吴彪，黄永，李惠，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江,23