一种超声波探伤用回波信号分离方法与流程

本发明属于超声波探伤技术领域，尤其是涉及一种超声波探伤用回波信号分离方法。

背景技术：

超声波探伤(也称超声检测或超声波检测)是利用材料及其缺陷的声学性能差异对超声波传播波形反射情况和穿透时间的能量变化来检验材料内部缺陷的无损检测方法。超声波探伤的原理是超声波在固体中的传输损失很小，探测深度大，由于超声波在异质界面上会发生反射、折射等现象，尤其是不能通过气体固体界面，如果金属中有气孔、裂纹、分层等缺陷(缺陷中有气体)或夹杂，超声波(也称超声波信号或超声信号)传播到金属与缺陷的界面处时，就会全部或部分反射，反射回来的超声波信号(也称超声回波信号，简称回波信号)被探头接收，根据所接收超声回波信号的波形变化特征便能判断缺陷在工件中的深度、位置和形状。

由于金属材料的特性和机械加工工艺固有缺陷等因素的影响，工业机械设备的生产过程中不可避免地存在不同程度的内部缺陷，常见的缺陷包括裂纹、气孔、缩孔、夹杂、夹砂冷隔等。例如，采煤机箱体等煤矿机械设备，生产过程中不可避免地存在上述缺陷，由于上述煤矿机械设备的制造工艺复杂，并且长期处于重载和交变载荷的状态下工作，工作环境恶劣，设备受力状况差，造成设备在使用过程中重要部件从小损伤逐渐发展成大事故，产生各种安全隐患，尤其是由于这些故障的产生，导致生产停产，直接影响着煤炭的产量和煤矿的经济效益。

现如今，在采煤机箱体缺陷检测中广泛使用的是超声检测方法。超声检测是一种重要的无损检测方法，超声信号是一个被探头中心频率调制的宽带脉冲信号，其回波信号中包含了大量与缺陷有关的信息，但是超声波信号(即超声回波信号)常常被测量系统和测试工件的随机噪声和相关噪声污染，特别是粗晶材料中的晶粒噪声，这些噪声会使超声信号的缺陷识别变得困难，严重制约缺陷检测的精度和可靠性。因此，对超声检测回波信号进行去噪处理，保证所得缺陷信号的真实性十分重要。这对企业提高产品缺陷检出率、保证产品质量和产品使用寿命具有重要意义。

由上述内容可知，对采煤机等煤矿机械的缺陷进行超声检测过程中，噪声严重影响了检测结果的可靠性和质量。从混有各种干扰噪声的信号(即超声检测信号)中提取出缺陷信号是确保回波信号准确性的关键，用超声波对材料内部缺陷进行检测时，缺陷信息将通过接收到的超声回波信号的幅值、频率或相位来表现。上述缺陷信号指的是检测到的超声回波信号，该超声回波信号中含有缺陷信息。但是由于仪器噪声和测试环境因素的影响，检测信号会伴有各种电噪声、结构噪声以及脉冲噪声等各种干扰噪声，特别是当缺陷信号的微结构噪声很大，或者材料本身晶粒粗大时，缺陷与噪声的信号强度相比会显得微弱，这种强噪声背景下缺陷信号的提取是超声信号处理领域研究的难题。如何从这些信号中提取所需信息是一项困难而又重要的课题。因此对超声检测回波信号进行去噪处理，保证所得缺陷信号的真实性十分重要。

目前，已有许多超声信号提取方法，如非线性滤波、傅里叶变换和小波变换等，这些方法对提高一般超声信号的信噪比具有较好效果，但是对于小缺陷或者强噪声背景下缺陷的提取具有局限性，检测结果不准确，可靠性不高。稀疏分解是一种新的信号分析理论，可以根据待提取信号的特征自适应的选择合适的展开函数，用很少的函数可以表示出所提取信号的基本特征，它可以在低信噪比情况下，更好地提取出弱小信号，最大化的逼近原始信号。稀疏分解算法首先是由mallat提出的，它是众所周知的匹配追踪算法。因此，在产品内部缺陷检测中已逐渐推广。但是该算法有两个缺陷，一是稀疏分解算法的计算量很大，计算时间在目前现有计算条件下十分巨大，无法进行实时检测；二是稀疏分解算法是连续条件下求的最优解，对于弱小缺陷的检测精度仍然有局限。

在超声无损检测中，超声波探伤装置中的超声波传感器发送脉冲通过被测对象不连续界面，则接受的反射回波(即超声回波信号)中包含与缺陷特征有关的位置信息和缺陷尺寸信息，因此准确的检测缺陷的位置和大小以对缺陷进行精确估计是超声无损评价的重要内容。并且，超声回波信号表示为沿时间轴在不同深度位置不同界面反射的缺陷波(也称为缺陷信号)的叠加组合。在实际应用中，经常出现一个超声回波信号中包含多个缺陷的情形，多个缺陷所处位置存在重叠，相应地该超声回波信号中对应多个缺陷的缺陷信号之间同样存在重叠，此时需要将该超声回波信号中对应多个缺陷的缺陷信号分离出来，否则便无法进行缺陷信号提取，并无从进行缺陷分类识别。待该超声回波信号中的多个缺陷信号均分离出来后，再对每个分离出的缺陷信号分别进行信号提取、特征提取与缺陷分类识别。

另外，目前采用稀疏分解方法进行信号去燥(也称为信号提取)时，其本质均是采用匹配追踪算法进行信号提取，由于在匹配追踪算法中，首先gabor原子的选择没有物理意义，其次如果在某一次迭代中选择了一个错误的gabor原子，则后续的迭代过程则不得不修正这个错误，导致误跟踪，以至算法性能差。并且，根据稀疏信号表示理论，字典原子的大小小于一个阈值时，由于累计一致性约束，则稀疏表示是不稳定不可靠的。因而，现有的稀疏分解方法存在诸多缺陷。

技术实现要素：

本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种超声波探伤用回波信号分离方法，其方法步骤简单、设计合理且实现方便、使用效果好，利用信号稀疏分解实现同一个超声回波信号中多个缺陷所处位置处的叠加信号的简便、快速分离。

为解决上述技术问题，本发明采用的技术方案是：一种超声波探伤用回波信号分离方法，其特征在于，该方法包括以下步骤：

步骤一、超声回波信号获取及同步上传：采用超声波探伤装置对被测对象进行超声波检测，获得被测对象的超声回波信号f(t)，并将所获得的超声回波信号f(t)同步传送至数据处理设备；

其中，f(t)＝[f(t1),f(t2),...,f(tn')]^t，t表示时间参数，ti为超声波探伤装置的第i个采样时刻，f(ti)为超声波探伤装置的第i个采样时刻采样得到的信号值，i为正整数且i＝1、2、3、…、n'，n'为正整数且其为超声回波信号f(t)的信号长度；

所述超声回波信号f(t)中包括被测对象中k个缺陷所处位置处的超声回波信号；

步骤二、信号稀疏分解：采用数据处理设备且调用稀疏分解模块对步骤一中所述超声回波信号f(t)进行迭代分解处理，将超声回波信号f(t)转换为并获得此时的迭代分解最佳原子集合；此时的所述迭代分解最佳原子集合中包含k个最佳匹配原子，为所述迭代分解最佳原子集合中的第z个所述最佳匹配原子；

式中r^k(t)为超声回波信号f(t)经过k次迭代分解后的残差量，其中k为预先设定的迭代分解总次数且k为正整数，z为正整数且z＝1、2、…、k；az为第z次迭代分解后的最佳匹配原子与上一次迭代分解后残差量的展开系数；

为第z次迭代分解时找出的最佳匹配原子，为gabor原子且式中函数ψ(t)为高斯窗函数且rz为的时频参数，rz＝(sz,uz,vz,wz)，sz为尺度参数，uz为位移参数，vz为频率参数，wz为相位参数；

步骤三、回波信号分离：采用数据处理设备且根据此时所述迭代分解最佳原子集合中的k个所述最佳匹配原子，对步骤一中所述超声回波信号f(t)进行分离，获得k个分离信号；

k个所述分离信号中第k'个所述分离信号记作fk'(t)，其中，为此时所述迭代分解最佳原子集合中的第k'个所述最佳匹配原子，k'为正整数且k'＝1、2、…、k；表示r^k'^-1(t)与的内积；r^k'^-1(t)为调用匹配追踪算法模块且利用此时所述迭代分解最佳原子集合中的前k'-1个所述最佳匹配原子对超声回波信号f(t)进行k'-1次迭代分解后的残差量，r⁰(t)＝f(t)。

上述一种超声波探伤用回波信号分离方法，其特征是：步骤一中所述超声波探伤装置与数据处理设备通过有线或无线方式进行连接。

上述一种超声波探伤用回波信号分离方法，其特征是：步骤四中每个所述分离信号均为被测对象中一个缺陷所处位置处的超声回波信号；

步骤一中被测对象中k个缺陷所处位置处的超声回波信号相互重叠；

步骤一中所述超声回波信号f(t)为被测对象中k个缺陷所处位置处的超声回波信号叠加后的信号，步骤四中每个所述分离信号均为一个叠加信号。

上述一种超声波探伤用回波信号分离方法，其特征是：步骤二中进行信号稀疏分解之前，先对所述超声回波信号f(t)中所述叠加信号的数量k进行确定；

对所述超声回波信号f(t)中所述叠加信号的数量k进行确定时，采用数据处理设备进行确定，过程如下：

步骤b1、波峰与波谷确定：采用数据处理设备对步骤一中所述超声回波信号f(t)的所有波峰与所有波谷分别进行确定，并对所确定的每个波峰与每个波谷的采样时刻与信号值分别进行同步记录；

本步骤中，所确定的每个波峰与每个波谷均为所述超声回波信号f(t)的一个极值点；

步骤b2、极值点剔除：采用数据处理设备且调用时域极值点剔除模块或频域极值点剔除模块进行极值点剔除，获得剔除后的m'个极值点，并按照各极值点的采样时间先后顺序对m'个所述极值点由前至后进行排列；其中，m'为正整数且其为本步骤中进行极值点剔除后获得的极值点的总数量；

其中，采用数据处理设备且调用时域极值点剔除模块进行极值点剔除时，将步骤b1中所确定的所有极值点中信号值的绝对值小于β'的极值点剔除，获得剔除后的m'个极值点；其中，β'＝α'×max|f(t)|，α'为常数且其取值范围为0.1～0.35，max|f(t)|为超声回波信号f(t)中信号值的绝对值最大值；

采用数据处理设备且调用频域极值点剔除模块进行极值点剔除时，将步骤b1中所确定的所有极值点中信号值的绝对值小于β的极值点剔除，获得剔除后的m'个极值点；其中，β为预先设定的剔除阈值且β＝α×max|y(f)|，α为常数且其取值范围为0.25～0.35，y(f)为超声回波信号f(t)的频谱，max|y(f)|为超声回波信号f(t)的频谱中幅值的绝对值最大值；

步骤b3、叠加判断：采用数据处理设备且按照采样时间先后顺序，由前至后对步骤b2中m'个所述极值点中相邻两个所述极值点分别进行叠加判断，并根据叠加判断结果对m'个所述极值点中叠加极值点对的总数量进行确定；

m'个所述极值点中相邻两个所述极值点的叠加判断方法均相同；对m'个所述极值点中第m'个所述极值点与第m'+1个所述极值点进行叠加判断时，采用数据处理设备且调用波峰判断模块或波谷判断模块进行叠加判断；其中，m'为正整数且m'＝1、2、…、m'-2、m'-1；

采用数据处理设备且调用波峰判断模块进行叠加判断时，当第m'个所述极值点与第m'+1个所述极值点均为波峰时，判断为第m'个所述极值点与第m'+1个所述极值点为一个所述叠加极值点对；

采用数据处理设备且调用波谷判断模块进行叠加判断时，当第m'个所述极值点与第m'+1个所述极值点均为波谷时，判断为第m'个所述极值点与第m'+1个所述极值点为一个所述叠加极值点对；

待m'个所述极值点中所有相邻两个所述极值点均完成叠加判断后，获得m'个所述极值点中的所有叠加极值点对，并获得m'个所述极值点中所述叠加极值点对的总数量，所获得的m'个所述极值点中所述叠加极值点对的总数量记作k'；再根据m'个所述极值点中所述叠加极值点对的总数量k'，并按照公式k＝k'+1，对所述超声回波信号f(t)中所述叠加信号的数量k进行确定。

上述一种超声波探伤用回波信号分离方法，其特征是：步骤b3中进行叠加判断之前，先采用数据处理设备对m'个所述极值点中所述叠加极值点对的总数量k'进行初始化，此时k'＝0；

对m'个所述极值点中第m'个所述极值点与第m'+1个所述极值点进行叠加判断时，当第m'个所述极值点与第m'+1个所述极值点均为波峰或均为波谷时，判断为第m'个所述极值点与第m'+1个所述极值点为一个所述叠加极值点对，并将此时的k'加1；否则，此时的k'保持不变；

待m'个所述极值点中所有相邻两个所述极值点均完成叠加判断后，此时的k'为m'个所述极值点中所述叠加极值点对的总数量。

上述一种超声波探伤用回波信号分离方法，其特征是：步骤二中进行信号稀疏分解时，采用数据处理设备由先至后对所述迭代分解最佳原子集合中k个所述最佳匹配原子分别进行查找；

所述的为采用数据处理设备查找出的最佳匹配原子，的时频参数rz为最佳时频参数；

采用数据处理设备且对进行查找时，对的时频参数rz进行查找，过程如下：

步骤c1、时频参数寻优：调用寻优算法模块且根据预先设定的时频参数中尺度参数、位移参数、频率参数和相位参数的取值范围，找出适应度值最大的最优时频参数，并找出适应度值＞thrdz的多个较优时频参数；

本步骤中，所找出的所述较优时频参数的数量记作kz，kz为整数且kz≥0；

其中，c'为常数且2≤c'≤3；rz-1为此时采用数据处理设备已查找出的前z-1个所述最佳匹配原子的残差量，rz-1＝f(t)-ψz-1(t)，ψz-1(t)为此时采用数据处理设备已确定的前z-1个所述最佳匹配原子之和，||rz-1||2为rz-1的2-范数；当z＝1时，r0＝f(t)；

所述最优时频参数为所找出的适应度值最大的时频参数，所述较优时频参数为所找出的适应度值＞thrdz的时频参数，每个所述时频参数的适应度值均为该时频参数所对应的gabor原子与rz-1的内积；

步骤c2、最佳时频参数确定：首先对步骤c1中所述较优时频参数的数量kz进行判断：当kz＝0或kz＝1时，将步骤c1中找出的最优时频参数作为最佳时频参数rz；当kz≥2时，对步骤c1中找出的kz个所述较优时频参数的稀疏度分别进行计算，并将稀疏度最小的所述较优时频参数作为最佳时频参数rz；

对kz个所述较优时频参数中第d个所述较优时频参数的稀疏度xsd进行计算时，xsd＝||r^d||ξ，r^d为ψd(t)的残差量，r^d＝f(t)-ψd(t)，ψd(t)为kz个所述较优时频参数中第d个所述较优时频参数对应的gabor原子；||r^d||ξ表示r^d的ξ-范数，ξ为常数且0≤ξ≤1；d为正整数且d＝1、2、…、kz。

上述一种超声波探伤用回波信号分离方法，其特征是：步骤三中进行回波信号分离之前，还需进行迭代分解最佳原子集合更新判断；

实际进行迭代分解最佳原子集合更新判断时，判断||r^k(t)||2是否小于ε：当||r^k(t)||2＜ε时，判断为无需对步骤二中所述迭代分解最佳原子集合进行更新，并进入步骤三；否则，当||r^k(t)||2≥ε时，判断为需对步骤二中所述迭代分解最佳原子集合进行更新，并进入步骤d；

其中，||r^k(t)||2为步骤二中所述r^k(t)的2-范数，ε为预先设定的残差量判断阈值；

步骤d、迭代分解最佳原子集合更新：采用数据处理设备对步骤二中所述迭代分解最佳原子集合进行更新，包括以下步骤：

步骤d1、最佳匹配原子优化，过程如下：

步骤d11、原子随机选取：采用数据处理设备从此时的所述迭代分解最佳原子集合中随机取出一个最佳匹配原子作为待优化原子，所述待优化原子记作其中j1为正整数且1≤j1≤k；

此时的所述迭代分解最佳原子集合中除所述待优化原子之外的k-1个最佳匹配原子均为待处理原子，k-1个待处理原子组成此时的待处理原子集合；

步骤d12、最佳匹配原子查找：所找出的最佳匹配原子记作的时频参数记作时频参数rj1'，时频参数rj1'＝(sj1',uj1',vj1',wj1')；

对最佳匹配原子进行查找时，根据预先设定的sj1'、uj1'、vj1'和wj1'的取值范围，采用数据处理设备且调用所述寻优算法模块，找出使适应度值fitness(rj1')最大的最佳时频参数，所找出的最佳时频参数为时频参数rj1'；再根据公式求解出最佳匹配原子

其中，表示与的内积；ψ'0(t)为步骤d11中k-1个所述待处理原子之和；

步骤d13、原子替换判断、原子替换及迭代分解最佳原子集合更新：采用数据处理设备且调用残值判断模块、适应度值判断模块或稀疏性判断模块，判断是否需对步骤d11中所述待优化原子进行替换，并根据判断结果对所述待优化原子进行替换；

采用数据处理设备且调用残值判断模块判断是否需对步骤d11中所述待优化原子进行替换时，根据替换后残值||rj1'^k(t)||ξ是否小于替换前残值||rj1^k(t)||ξ进行判断：当||rj1'^k(t)||ξ＜||rj1^k(t)||ξ时，判断为需对步骤d11中所述待优化原子进行替换，将步骤d11中所述待优化原子替换为步骤d12中所述最佳匹配原子获得更新后的所述迭代分解最佳原子集合；否则，判断为无需对步骤d11中所述待优化原子进行替换，并进入步骤d15；

其中，rj1'^k(t)＝f(t)-ψj1'(t)，rj1^k(t)＝f(t)-ψj1(t)，ψj1(t)为本步骤中进行原子替换判断之前所述迭代分解最佳原子集合中k个最佳匹配原子之和；||rj1'^k(t)||ξ表示rj1'^k(t)的ξ-范数，||rj1^k(t)||ξ表示rj1^k(t)的ξ-范数，ξ为常数且0≤ξ≤1；

采用数据处理设备且调用适应度值判断模块判断是否需对步骤d11中所述待优化原子进行替换时，根据替换后适应度值fitness(rj1')是否大于替换前适应度值fitness(rj1)进行判断：当fitness(rj1')＞fitness(rj1)时，判断为需对步骤d11中所述待优化原子进行替换，将步骤d11中所述待优化原子替换为步骤d12中所述最佳匹配原子获得更新后的所述迭代分解最佳原子集合；否则，判断为无需对步骤d11中所述待优化原子进行替换，并进入步骤d15；

其中，表示r^j1-1(t)与的内积，r^j1-1(t)＝f(t)-ψj1-1(t)，ψj1-1(t)为此时所述迭代分解最佳原子集合中前j1-1个最佳匹配原子之和；表示r^j1-1(t)与的内积；

采用数据处理设备且调用稀疏性判断模块判断是否需对步骤d11中所述待优化原子进行替换时，根据||r^j1'||ξ是否小于||r^j1||ξ进行判断：当||r^j1'||ξ＜||r^j1||ξ时，判断为需对步骤d11中所述待优化原子进行替换，将步骤d11中所述待优化原子替换为步骤d12中所述最佳匹配原子获得更新后的所述迭代分解最佳原子集合；否则，判断为无需对步骤d11中所述待优化原子进行替换，并进入步骤d15；

其中，r^j1'为的残差量且r^j1为的残差量且||r^j1'||ξ表示r^j1'的ξ-范数，||r^j1||ξ表示r^j1的ξ-范数；

本步骤中，完成原子替换判断、原子替换及迭代分解最佳原子集合更新后，完成步骤d11中选取的一个所述最佳匹配原子的优化过程；

步骤d14、残差量判断：对步骤d13中所述最佳匹配原子优化后的残差量进行判断：当||r'j1^k(t)||2＜ε时，完成迭代分解最佳原子集合更新过程，进入步骤三；否则，当||r'j1^k(t)||2≥ε时，进入步骤d15；

其中，||r'j1^k(t)||2为r'j1^k(t)的2-范数；r'j1^k(t)为根据此时所述迭代分解最佳原子集合中的k个所述最佳匹配原子对f(t)进行k次迭代分解后的残差量；

步骤d15、下一个最佳匹配原子优化：按照步骤d11中至步骤d13中所述的方法，对此时所述迭代分解最佳原子集合中未进行优化的一个所述最佳匹配原子进行优化；

步骤d16、残差量判断：对步骤d15中所述最佳匹配原子优化后的残差量进行判断：当||r”j1^k(t)||2＜ε时，完成迭代分解最佳原子集合更新过程，进入步骤三；否则，当||r”j1^k(t)||2≥ε时，返回步骤d15；

其中，||r”j1^k(t)||2为r”j1^k(t)的2-范数；r”j1^k(t)为根据此时所述迭代分解最佳原子集合中的k个所述最佳匹配原子对f(t)进行k次迭代分解后的残差量。

上述一种超声波探伤用回波信号分离方法，其特征是：步骤d1中进行信号稀疏分解后，采用数据处理设备将所述迭代分解最佳原子集合同步存储至数据存储器内，所述数据存储器与数据处理设备连接；

步骤d13中进行原子替换判断及原子替换后，采用数据处理设备对更新后的所述迭代分解最佳原子集合进行同步存储；

步骤d1中进行信号稀疏分解后，采用数据处理设备将所述迭代分解最佳原子集合同步存储至数据存储器内时，按照迭代分解先后顺序对所述迭代分解最佳原子集合中的k个所述最佳匹配原子分别进行存储；其中，为步骤d1中对所述超声回波信号f(t)进行第z次迭代分解时找出的最佳匹配原子。

上述一种超声波探伤用回波信号分离方法，其特征是：步骤d1中最佳匹配原子优化时，按照存储先后顺序对所述迭代分解最佳原子集合中的所述最佳匹配原子进行优化；

步骤d1中最佳匹配原子优化时，最先进行优化的所述最佳匹配原子为步骤d1中所述迭代分解最佳原子集合中的第1个所述最佳匹配原子。

上述一种超声波探伤用回波信号分离方法，其特征是：步骤三中回波信号分离完成后，还需对k个所述分离信号分别进行信号提取；k个所述分离信号的信号提取方法均相同；

对k个所述分离信号中任一个所述分离信号进行信号提取时，采用数据处理设备对该分离信号进行信号提取，该分离信号为待处理信号且其记作信号f(t)；

采用数据处理设备对信号f(t)进行信号提取时，过程如下：

步骤a1、基于寻优算法的稀疏分解：采用数据处理设备且调用稀疏分解模块对信号f(t)进行迭代分解处理，将信号f(t)转换为并获得此时的迭代分解最佳原子集合；此时的所述迭代分解最佳原子集合中包含m个最佳匹配原子，为所述迭代分解最佳原子集合中的第n个所述最佳匹配原子；

式中r^m(t)为信号f(t)经过m次迭代分解后的残差量，其中m为预先设定的迭代分解总次数且m为正整数，n为正整数且n＝1、2、…、m；an为第n次迭代分解后的最佳匹配原子与上一次迭代分解后残差量的展开系数；

为第n次迭代分解时采用数据处理设备且调用寻优算法模块找出的最佳匹配原子；为gabor原子且式中函数ψ(t)为高斯窗函数且rn为的时频参数，rn＝(sn,un,vn,wn)，sn为尺度参数，un为位移参数，vn为频率参数，wn为相位参数；

本步骤中，寻找最佳匹配原子时，根据预先设定的sn、un、vn和wn的取值范围，采用数据处理设备且调用寻优算法模块，找出使适应度值fitness(rn)最大的最佳时频参数，所找出的最佳时频参数为时频参数rn；

其中，fitness(rn)为时频参数rn的适应度值，表示rⁿ-¹(t)与的内积；r^n-1(t)为信号f(t)经过n-1次迭代分解后的残差量，r⁰(t)＝f(t)；

步骤a2、残差量判断：判断||r^m(t)||2是否小于ε：当||r^m(t)||2＜ε时，进入步骤a4；否则，当||r^m(t)||2≥ε时，进入步骤a3；

其中，||r^m(t)||2为步骤a1中所述r^m(t)的2-范数，ε为预先设定的残差量判断阈值；

步骤a3、最佳匹配原子优化，过程如下：

步骤a31、原子随机选取：采用数据处理设备从此时的所述迭代分解最佳原子集合中随机取出一个最佳匹配原子作为待优化原子，所述待优化原子记作其中j为正整数且1≤j≤m；

此时的所述迭代分解最佳原子集合中除所述待优化原子之外的m-1个最佳匹配原子均为待处理原子，m-1个待处理原子组成此时的待处理原子集合；

步骤a32、最佳匹配原子查找：所找出的最佳匹配原子记作的时频参数记作时频参数rj'，时频参数rj'＝(sj',uj',vj',wj')；

对最佳匹配原子进行查找时，根据预先设定的sj'、uj'、vj'和wj'的取值范围，采用数据处理设备且调用所述寻优算法模块，找出使适应度值fitness(rj')最大的最佳时频参数，所找出的最佳时频参数为时频参数rj'；再根据公式求解出最佳匹配原子

其中，表示与的内积；ψ0(t)为步骤a31中m-1个所述待处理原子之和；

步骤a33、原子替换判断及原子替换：采用数据处理设备且调用残值判断模块、适应度值判断模块或稀疏性判断模块，判断是否需对步骤a31中所述待优化原子进行替换，并根据判断结果对所述待优化原子进行替换；

采用数据处理设备且调用残值判断模块判断是否需对步骤a31中所述待优化原子进行替换时，根据替换后残值||rj'^m(t)||ξ是否小于替换前残值||rj^m(t)||ξ进行判断：当||rj'^m(t)||ξ＜||rj^m(t)||ξ时，判断为需对步骤a31中所述待优化原子进行替换，将步骤a31中所述待优化原子替换为步骤a32中所述最佳匹配原子获得更新后的所述迭代分解最佳原子集合；否则，判断为无需对步骤a31中所述待优化原子进行替换，并进入步骤a35；

其中，rj'^m(t)＝f(t)-ψj'(t)，rj^m(t)＝f(t)-ψj(t)，ψj(t)为本步骤中进行原子替换判断之前所述迭代分解最佳原子集合中m个最佳匹配原子之和；||rj'^m(t)||ξ表示rj'^m(t)的ξ-范数，||rj^m(t)||ξ表示rj^m(t)的ξ-范数，ξ为常数且0≤ξ≤1；

采用数据处理设备且调用适应度值判断模块判断是否需对步骤a31中所述待优化原子进行替换时，根据替换后适应度值fitness(rj')是否大于替换前适应度值fitness(rj)进行判断：当fitness(rj')＞fitness(rj)时，判断为需对步骤a31中所述待优化原子进行替换，将步骤a31中所述待优化原子替换为步骤a32中所述最佳匹配原子获得更新后的所述迭代分解最佳原子集合；否则，判断为无需对步骤a31中所述待优化原子进行替换，并进入步骤a35；

其中，表示r^j-1(t)与的内积，r^j-1(t)＝f(t)-ψj-1(t)，ψj-1(t)为此时所述迭代分解最佳原子集合中前j-1个最佳匹配原子之和；表示r^j-1(t)与的内积；

采用数据处理设备且调用稀疏性判断模块判断是否需对步骤a31中所述待优化原子进行替换时，根据||r^j'||ξ是否小于||r^j||ξ进行判断：当||r^j'||ξ＜||r^j||ξ时，判断为需对步骤a31中所述待优化原子进行替换，将步骤a31中所述待优化原子替换为步骤a32中所述最佳匹配原子获得更新后的所述迭代分解最佳原子集合；否则，判断为无需对步骤a31中所述待优化原子进行替换，并进入步骤a35；

其中，r^j'为的残差量且r^j为的残差量且r^j'||ξ表示r^j'的ξ-范数，||r^j||ξ表示r^j的ξ-范数；

本步骤中，完成原子替换判断及原子替换后，完成步骤a31中选取的一个所述最佳匹配原子的优化过程；

步骤a34、残差量判断：对步骤a33中所述最佳匹配原子优化后的残差量进行判断：当||r'j^m(t)||2＜ε时，进入步骤a4；否则，当||r'j^m(t)||2≥ε时，进入步骤a35；

其中，||r'j^m(t)||2为r'j^m(t)的2-范数；r'j^m(t)为根据此时所述迭代分解最佳原子集合中的m个所述最佳匹配原子对f(t)进行m次迭代分解后的残差量；

步骤a35、下一个最佳匹配原子优化：按照步骤a31中至步骤a33中所述的方法，对此时所述迭代分解最佳原子集合中未进行优化的一个所述最佳匹配原子进行优化；

步骤a36、残差量判断：对步骤a35中所述最佳匹配原子优化后的残差量进行判断：当||r”j^m(t)||2＜ε时，进入步骤a4；否则，当||r”j^m(t)||2≥ε时，返回步骤a35；

其中，||r”j^m(t)||2为r”j^m(t)的2-范数；r”j^m(t)为根据此时所述迭代分解最佳原子集合中的m个所述最佳匹配原子对f(t)进行m次迭代分解后的残差量；

步骤a4、信号重构：根据此时的所述迭代分解最佳原子集合，采用数据处理设备得出信号f(t)的近似信号f″(t)；其中，近似信号f″(t)为从信号f(t)提取出的信号，

其中为此时所述迭代分解最佳原子集合中的第n'个所述最佳匹配原子，n'为正整数且n'＝1、2、…、m；an'为与根据此时所述迭代分解最佳原子集合中的前n'-1个所述最佳匹配原子对f(t)进行n'-1次迭代分解后残差量的展开系数。

本发明与现有技术相比具有以下优点：

1、方法步骤简单、设计合理且实现方便，投入成本较低。

2、信号分离速度快，采用数据处理器自动完成信号分离过程，并且能在几分钟甚至更短的时间内完成信号分离过程，实现信号实时分离。

3、所采用的基于寻优算法的信号稀疏分解方法简单、设计合理且实现方便、使用效果好，不仅提高了信号提取的速度，而且经信号提取后能有效提高原信号的质量和性能指标，尤其在超声无损探伤中具有重要的作用。同时，将频率参数v的取值范围限定为并且fo的单位为mhz，一方面能有效减少稀疏分解算法的计算量，实现进行实时检测；另一方面，有效提高mp算法(即匹配追踪算法)的性能，使得稀疏表示的信号能有效满足弱小缺陷的检测精度，达到更简易、更准确地获得信号中所蕴含的有效信息。通过将频率参数v的取值范围限定为能进一步凸显信号中所蕴含有效的信息，使稀疏表示的信号更加侧重表达有效信息，弱化冗余信息，从而能更准确地表达信号本征特征，使分离信号的精度能得到有效保证。

4、基于寻优算法查找最佳匹配原子时，除考虑适应度值这一指标外，还需考虑信号稀疏度，以便选择与信号f(t)最优匹配的原子，从而使得去燥后的信号(或所提取信号)能更准确地表达信号本征特征，使信号提取精度能得到有效保证。

5、增加分离信号精度验证步骤且验证方法简单、易于实现，信号稀疏分解后还需通过残差量判断判断此时的迭代分解最佳原子集合是否满足预先设定的信号提取精度要求，并根据判断结果进行最佳匹配原子优化，从而能进一步提高信号提取的准确性，进一步提高信号提取精度，使提取后信号进一步逼近原始信号，实现与原始信号的最佳匹配，提高信号的提取准确性及提取速度，确保信号分离效果。

6、所采用的最佳匹配原子优化方法设计合理、实现方便且使用效果好，从此时的迭代分解最佳原子集合中随机选取一个最佳匹配原子进行优化，当优化完成一个所述最佳匹配原子后，便通过残差量判断对此时迭代分解最佳原子集合是否满足信号提取精度要求进行判断，并根据判断结果决定是否需要继续对其余的最佳匹配原子进行优化。因而，实现简便，并能实现快速优化与优化结果实时判断相结合，能有效简化最佳匹配原子优化过程，并能快速达到最佳匹配原子优化的目的，进一步有效提高信号提取精度。同时，所采用的原子替换判断方法设计合理、实现简便且使用效果好，采用残值判断、适应度值判断或稀疏性判断任一种方法进行原子替换判断，可选取任一种方法进行原子替换判断，使用方式灵活，并且每一种原子替换判断方法均能实现有效的原子替换判断。

7、采用改进的人工蜂群算法进行寻优实现最佳匹配原子查找，信号稀疏分解之前不需要生成过完备字典中的全部原子，只需要产生蜜源的位置来替代原子库中的gabor原子，大大节省了存储空间。此外，人工蜂群算法是在连续的空间中进行最佳匹配原子的搜寻，而匹配追踪算法是在离散的搜索空间中进行原子的搜寻，因此人工蜂群算法搜索的范围更广，提取出的原子能更好的反映原始信号的特点，不但提高了计算速度，由于其在连续解空间范围内寻优，因此，也提高了参数提取的准确性。相比离散空间范围，能更准确的提取出于信号匹配的最佳匹配原子，从而提高了信号提取的精度，能有效提取强噪声背景下的有用信号。

8、所采用的超声回波信号f(t)中叠加信号的数量k的确定方法设计合理、实现简便且使用效果好，先进行波峰与波谷确定，并相应获得超声回波信号f(t)中的所有极值点；再根据按照特定的极值点剔除方法进行有效地极值点剔除，不仅有效减小了步骤b3中叠加判断的计算量，并且能更加凸显叠加信号，通过b3中叠加判断识别出极值点剔除后的m'个所述极值点中所述叠加极值点对的总数量k'，并按照公式k＝k'+1，便可对所述超声回波信号f(t)中所述叠加信号的数量k进行确定。只有步骤b2中进行有效地极值点剔除与步骤b3中叠加判断有效结合，才能将m'个所述极值点中所述叠加极值点对与两个叠加信号的叠加之处一一对应，从而实现所述超声回波信号f(t)中所述叠加信号数量k的确定，使得所述超声回波信号f(t)中所述叠加信号的数量k既不会过大，也不会过小。步骤b2中进行有效地极值点剔除后，可将一个叠加信号的没有实际分析意义的极值点剔除，仅保留该叠加信号中具有真正价值的极值点，这样能有效避免所确定的所述叠加信号的数量k过大；同时，步骤b3中通过对m'个所述极值点中相邻两个所述极值点分别进行叠加判断，不漏过任何一个叠加之处，从而能有效避免所确定的所述叠加信号的数量k过小，因而所确定的所述叠加信号的数量k的准确性非常高。

9、信号分离后，还需对每个分离信号进行提取，并且所采用的信号提取方法步骤简单、设计合理且实现方便，投入成本较低，信号提取速度快，采用数据处理器自动完成信号提取过程，并且能在几分钟甚至更短的时间内完成信号提取过程，实现信号实时提取。实际进行信号提取时，采用的基于寻优算法的信号稀疏分解方法简单、设计合理且实现方便、使用效果好，不仅提高了信号提取的速度，而且经信号提取后能有效提高原信号的质量和性能指标，尤其在超声无损探伤中具有重要的作用。同时，将频率参数v的取值范围限定为并且fo的单位为mhz，一方面能有效减少稀疏分解算法的计算量，实现进行实时检测；另一方面，有效提高mp算法(即匹配追踪算法)的性能，使得稀疏表示的信号能有效满足弱小缺陷的检测精度，达到更简易、更准确地获得信号中所蕴含的有效信息。通过将频率参数v的取值范围限定为能进一步凸显信号中所蕴含有效的信息，使稀疏表示的信号更加侧重表达有效信息，弱化冗余信息，从而能更准确地表达信号本征特征，使信号提取精度能得到有效保证。

同时，增加信号提取精度验证步骤且验证方法简单、易于实现，信号稀疏分解后还需通过残差量判断判断此时的迭代分解最佳原子集合是否满足预先设定的信号提取精度要求，并根据判断结果进行最佳匹配原子优化，从而能进一步提高信号提取的准确性，进一步提高信号提取精度，使提取后信号进一步逼近原始信号，实现与原始信号的最佳匹配，提高信号的提取准确性及提取速度。

能采用改进的人工蜂群算法进行寻优实现对分离后的信号进行提取，尤其对于超声无损探伤而言，为缺陷准确检测提供可靠的依据，同时为缺陷的定性和定量分析提供理论基础，能有效解决强噪声背景下弱缺陷提取难的问题，能解决弱缺陷的提取速度和精度问题，能准确地提取出强噪声背景下的缺陷信息，提高了超声信号提取的速度，为实时自动化检测提供了技术支持，因而能有效解决现有匹配追踪算法存在的算法复杂度大，过匹配等问题。因此，本发明利用人工蜂群算法从连续字典库中选择与超声信号最优匹配的原子，从而恢复出待处理信号。

10、使用效果好且实用价值高，方法步骤简单、设计合理且实现方便、使用效果好，利用信号稀疏分解实现同一个超声回波信号中多个缺陷所处位置处的叠加信号的简便、快速分离。

下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。

附图说明

图1为本发明的方法流程框图。

图2为本发明所采用信号分离处理系统的电路原理框图。

附图标记说明:

1—超声波探伤装置；2—数据处理设备；3—数据存储器。

具体实施方式

如图1所示的一种超声波探伤用回波信号分离方法，包括以下步骤：

步骤一、超声回波信号获取及同步上传：采用超声波探伤装置1对被测对象进行超声波检测，获得被测对象的超声回波信号f(t)，并将所获得的超声回波信号f(t)同步传送至数据处理设备2；

其中，f(t)＝[f(t1),f(t2),...,f(tn')]^t，t表示时间参数，ti为超声波探伤装置1的第i个采样时刻，f(ti)为超声波探伤装置1的第i个采样时刻采样得到的信号值，i为正整数且i＝1、2、3、…、n'，n'为正整数且其为超声回波信号f(t)的信号长度；

其中，[]^t表示矩阵的转置；[f(t1),f(t2),...,f(tn')]^t为矩阵[f(t1),f(t2),...,f(tn')]的转置；

所述超声回波信号f(t)中包括被测对象中k个缺陷所处位置处的超声回波信号；

步骤二、信号稀疏分解：采用数据处理设备2且调用稀疏分解模块对步骤一中所述超声回波信号f(t)进行迭代分解处理，将超声回波信号f(t)转换为并获得此时的迭代分解最佳原子集合；此时的所述迭代分解最佳原子集合中包含k个最佳匹配原子，为所述迭代分解最佳原子集合中的第z个所述最佳匹配原子；

式中r^k(t)为超声回波信号f(t)经过k次迭代分解后的残差量，其中k为预先设定的迭代分解总次数且k为正整数，z为正整数且z＝1、2、…、k；az为第z次迭代分解后的最佳匹配原子与上一次迭代分解后残差量的展开系数；

为第z次迭代分解时找出的最佳匹配原子，为gabor原子且式中函数ψ(t)为高斯窗函数且rz为的时频参数，rz＝(sz,uz,vz,wz)，sz为尺度参数，uz为位移参数，vz为频率参数，wz为相位参数；

步骤三、回波信号分离：采用数据处理设备2且根据此时所述迭代分解最佳原子集合中的k个所述最佳匹配原子，对步骤一中所述超声回波信号f(t)进行分离，获得k个分离信号；

k个所述分离信号中第k'个所述分离信号记作fk'(t)，其中，为此时所述迭代分解最佳原子集合中的第k'个所述最佳匹配原子，k'为正整数且k'＝1、2、…、k；表示r^k'^-1(t)与的内积；r^k'^-1(t)为调用匹配追踪算法模块且利用此时所述迭代分解最佳原子集合中的前k'-1个所述最佳匹配原子对超声回波信号f(t)进行k'-1次迭代分解后的残差量，r⁰(t)＝f(t)。

稀疏分解算法(也称为mp算法)是通过多次迭代进行进行原子的逼近，每次迭代中添加一个单个的原子与超声回波信号进行匹配，使得选择的原子与信号残差的内积最大，从而使得残差的能量最大程度地减少，这种原子的选择方法使得原子的解释缺乏物理含义。因而，mp算法是利用原子的最小子集来匹配信号特征，恢复信号，不仅计算量大，更重要的是缺乏物理解释。如果能给出所选择的原子一个准确的物理意义，让原子与超声回波信号一一对应，通过从过完备原子库中依次选择占主导地位的原子，添加到信号逼近表示中，来恢复多个缺陷信号，这种具有物理解释的原子的选择使得多个缺陷信号的提取有了可靠依据。

本发明中将所述迭代分解最佳原子集合中的k个最佳匹配原子分别与超声回波信号f(t)中包含的k个分离信号(即被测对象中k个缺陷所处位置处的超声回波信号，也称为k个缺陷信号)一一对应，通过信号稀疏分解实现k个缺陷信号的分离，以便对k个缺陷信号分别进行准确地识别。

其中，步骤四中每个所述分离信号均为被测对象中一个缺陷所处位置处的超声回波信号。

步骤一中被测对象中k个缺陷所处位置处的超声回波信号相互重叠。

步骤一中所述超声回波信号f(t)为被测对象中k个缺陷所处位置处的超声回波信号叠加后的信号，步骤四中每个所述分离信号均为一个叠加信号。

本实施例中，步骤一中所述超声波探伤装置1与数据处理设备2通过有线或无线方式进行连接。

步骤二中所述的

实际使用时，当超声回波信号f(t)所包含叠加信号的数量k已知时，按照步骤一至步骤三中所述的方法，便能简便、快速完成超声回波信号f(t)中k个缺陷所处位置处的超声回波信号(即k个所述叠加信号)的分离过程；但当超声回波信号f(t)所包含叠加信号的数量k未知时，还需对超声回波信号f(t)中所述叠加信号的数量k进行确定。

本实施例中，步骤二中进行信号稀疏分解之前，先对所述超声回波信号f(t)中所述叠加信号的数量k进行确定；

对所述超声回波信号f(t)中所述叠加信号的数量k进行确定时，采用数据处理设备2进行确定，过程如下：

步骤b1、波峰与波谷确定：采用数据处理设备2对步骤一中所述超声回波信号f(t)的所有波峰与所有波谷分别进行确定，并对所确定的每个波峰与每个波谷的采样时刻与信号值分别进行同步记录；

本步骤中，所确定的每个波峰与每个波谷均为所述超声回波信号f(t)的一个极值点；

步骤b2、极值点剔除：采用数据处理设备2且调用时域极值点剔除模块或频域极值点剔除模块进行极值点剔除，获得剔除后的m'个极值点，并按照各极值点的采样时间先后顺序对m'个所述极值点由前至后进行排列；其中，m'为正整数且其为本步骤中进行极值点剔除后获得的极值点的总数量；

其中，采用数据处理设备2且调用时域极值点剔除模块进行极值点剔除时，将步骤b1中所确定的所有极值点中信号值的绝对值小于β'的极值点剔除，获得剔除后的m'个极值点；其中，β'＝α'×max|f(t)|，α'为常数且其取值范围为0.1～0.35，max|f(t)|为超声回波信号f(t)中信号值的绝对值最大值；

采用数据处理设备2且调用频域极值点剔除模块进行极值点剔除时，将步骤b1中所确定的所有极值点中信号值的绝对值小于β的极值点剔除，获得剔除后的m'个极值点；其中，β为预先设定的剔除阈值且β＝α×max|y(f)|，α为常数且其取值范围为0.25～0.35，y(f)为超声回波信号f(t)的频谱，max|y(f)|为超声回波信号f(t)的频谱中幅值的绝对值最大值；

步骤b3、叠加判断：采用数据处理设备2且按照采样时间先后顺序，由前至后对步骤b2中m'个所述极值点中相邻两个所述极值点分别进行叠加判断，并根据叠加判断结果对m'个所述极值点中叠加极值点对的总数量进行确定；

m'个所述极值点中相邻两个所述极值点的叠加判断方法均相同；对m'个所述极值点中第m'个所述极值点与第m'+1个所述极值点进行叠加判断时，采用数据处理设备2且调用波峰判断模块或波谷判断模块进行叠加判断；其中，m'为正整数且m'＝1、2、…、m'-2、m'-1；

采用数据处理设备2且调用波峰判断模块进行叠加判断时，当第m'个所述极值点与第m'+1个所述极值点均为波峰时，判断为第m'个所述极值点与第m'+1个所述极值点为一个所述叠加极值点对；

采用数据处理设备2且调用波谷判断模块进行叠加判断时，当第m'个所述极值点与第m'+1个所述极值点均为波谷时，判断为第m'个所述极值点与第m'+1个所述极值点为一个所述叠加极值点对；

待m'个所述极值点中所有相邻两个所述极值点均完成叠加判断后，获得m'个所述极值点中的所有叠加极值点对，并获得m'个所述极值点中所述叠加极值点对的总数量，所获得的m'个所述极值点中所述叠加极值点对的总数量记作k'；再根据m'个所述极值点中所述叠加极值点对的总数量k'，并按照公式k＝k'+1，对所述超声回波信号f(t)中所述叠加信号的数量k进行确定。

本实施例中，步骤b3中采用数据处理设备2且调用波峰判断模块进行叠加判断，当第m'个所述极值点与第m'+1个所述极值点均为波峰时，判断为第m'个所述极值点与第m'+1个所述极值点为一个所述叠加极值点对。

实际使用时，步骤b3中也可以采用数据处理设备2且调用波谷判断模块进行叠加判断，当第m'个所述极值点与第m'+1个所述极值点均为波谷时，判断为第m'个所述极值点与第m'+1个所述极值点为一个所述叠加极值点对。

本实施例中，步骤b2中采用数据处理设备2且调用时域极值点剔除模块进行极值点剔除，并且可根据实际需要，对α的取值大小进行相应调整。

实际使用时，步骤b2中也可以采用数据处理设备2且调用频域极值点剔除模块进行极值点剔除，并且可根据实际需要，对α'的取值大小进行相应调整。

由于信号采样过程中，受环境因素、采样系统等多方面因素影响，所述超声回波信号f(t)中不可避免会存在多个非真正的极值点(即干扰极值点)，上述极值点并非所述超声回波信号f(t)中的缺陷信号的极值点，因而需要将上述非真正的极值点剔除，这样不仅能有效减少计算量，并且能剔除干扰极值点，以便对所存在的叠加信号数量k进行准确确定。

本发明所采用时域剔除与频域剔除两种剔除方法，均能简便、快速且有效剔除干扰极值点，上述两种剔除方法均采用数据处理设备2便能自动完成，因而实现方便，处理速度快，并且剔除结果能有效保证。

采用时域剔除方法或频域剔除方法进行极值点剔除时，所采用的剔除阈值β'和剔除阈值β均设定合理。剔除阈值β'和剔除阈值β的确定均跟超声回波信号f(t)直接相关，针对不同的待处理信号并非采用一个固定值，而是由当前所处理的超声回波信号f(t)进行确定的，具有一定的适应性，能有效确保剔除效果，确保不会过多剔除真正的极值点。并且，结合实际的极值点剔除效果与超声缺陷信号识别目的，上述极值点剔除方法的实际使用效果非常好。

其中，β'＝α'×max|f(t)|，β'由超声回波信号f(t)中信号值的绝对值最大值max|f(t)|决定，根据max|f(t)|能对当前所处理的超声回波信号f(t)中所包含叠加信号的信号值进行了解，而信号值的绝对值小于β'的极值点均与超声回波信号f(t)中所包含叠加信号基本上无关联，实际分析的意义基本不存在，因而将信号值的绝对值小于β'的极值点剔除。

相应地，β＝α×max|y(f)|，β由超声回波信号f(t)中幅值的绝对值最大值max|y(f)|决定，根据max|y(f)|能对当前所处理的超声回波信号f(t)中所包含叠加信号的幅值进行了解，信号值的绝对值小于β的极值点均与超声回波信号f(t)中所包含叠加信号基本上无关联，实际分析的意义基本不存在，因而将信号值的绝对值小于β的极值点剔除。

步骤b3中进行叠加判断时，所采用的叠加判断方法设计合理、实现方便且使用效果好，通过对剔除后的m'个所述极值点分别进行叠加判断，识别出m'个所述极值点中所述叠加极值点对的总数量k'，根据m'个所述极值点中所述叠加极值点对的总数量k'便能直接确定所述超声回波信号f(t)中所述叠加信号的数量k。

上述步骤b1至步骤b3中采用的超声回波信号f(t)中所述叠加信号的数量k的确定方法设计合理、实现简便且使用效果好，先进行波峰与波谷确定，并相应获得超声回波信号f(t)中的所有极值点；再根据按照特定的极值点剔除方法进行有效地极值点剔除，不仅有效减小了步骤b3中叠加判断的计算量，并且能更加凸显叠加信号，通过b3中叠加判断识别出极值点剔除后的m'个所述极值点中所述叠加极值点对的总数量k'，并根据公式k＝k'+1，便可对所述超声回波信号f(t)中所述叠加信号的数量k进行确定。只有步骤b2中进行有效地极值点剔除与步骤b3中叠加判断有效结合，才能将m'个所述极值点中所述叠加极值点对与两个叠加信号的叠加之处一一对应，从而实现所述超声回波信号f(t)中所述叠加信号数量k的确定，使得所述超声回波信号f(t)中所述叠加信号的数量k既不会过大，也不会过小。步骤b2中进行有效地极值点剔除后，可将一个叠加信号的没有实际分析意义的极值点剔除，仅保留该叠加信号中具有真正价值的极值点，这样能有效避免所确定的所述叠加信号的数量k过大；同时，步骤b3中通过对m'个所述极值点中相邻两个所述极值点分别进行叠加判断，不漏过任何一个叠加之处，从而能有效避免所确定的所述叠加信号的数量k过小，因而所确定的所述叠加信号的数量k的准确性非常高。

本实施例中，步骤b3中进行叠加判断之前，先采用数据处理设备2对m'个所述极值点中所述叠加极值点对的总数量k'进行初始化，此时k'＝0；

对m'个所述极值点中第m'个所述极值点与第m'+1个所述极值点进行叠加判断时，当第m'个所述极值点与第m'+1个所述极值点均为波峰或均为波谷时，判断为第m'个所述极值点与第m'+1个所述极值点为一个所述叠加极值点对，并将此时的k'加1；否则，此时的k'保持不变；

待m'个所述极值点中所有相邻两个所述极值点均完成叠加判断后，此时的k'为m'个所述极值点中所述叠加极值点对的总数量。

本实施例中，步骤二中进行信号稀疏分解时，采用数据处理设备2由先至后对所述迭代分解最佳原子集合中k个所述最佳匹配原子分别进行查找；

所述的为采用数据处理设备2查找出的最佳匹配原子，的时频参数rz为最佳时频参数；

采用数据处理设备2且对进行查找时，对的时频参数rz进行查找，过程如下：

步骤c1、时频参数寻优：调用寻优算法模块且根据预先设定的时频参数中尺度参数、位移参数、频率参数和相位参数的取值范围，找出适应度值最大的最优时频参数，并找出适应度值＞thrdz的多个较优时频参数；

本步骤中，所找出的所述较优时频参数的数量记作kz，kz为整数且kz≥0；

其中，c'为常数且2≤c'≤3；rz-1为此时采用数据处理设备2已查找出的前z-1个所述最佳匹配原子的残差量，rz-1＝f(t)-ψz-1(t)，ψz-1(t)为此时采用数据处理设备2已确定的前z-1个所述最佳匹配原子之和，||rz-1||2为rz-1的2-范数；当z＝1时，r0＝f(t)；

所述最优时频参数为所找出的适应度值最大的时频参数，所述较优时频参数为所找出的适应度值＞thrdz的时频参数，每个所述时频参数的适应度值均为该时频参数所对应的gabor原子与rz-1的内积；

步骤c2、最佳时频参数确定：首先对步骤c1中所述较优时频参数的数量kz进行判断：当kz＝0或kz＝1时，将步骤c1中找出的最优时频参数作为最佳时频参数rz；当kz≥2时，对步骤c1中找出的kz个所述较优时频参数的稀疏度分别进行计算，并将稀疏度最小的所述较优时频参数作为最佳时频参数rz；

对kz个所述较优时频参数中第d个所述较优时频参数的稀疏度xsd进行计算时，xsd＝||r^d||ξ，r^d为ψd(t)的残差量，r^d＝f(t)-ψd(t)，ψd(t)为kz个所述较优时频参数中第d个所述较优时频参数对应的gabor原子；||r^d||ξ表示r^d的ξ-范数，ξ为常数且0≤ξ≤1；d为正整数且d＝1、2、…、kz。

根据本领域公知常识，2-范数是指向量各元素的平方和然后求平方根(即l2范数)。

步骤c1中所述的rz-1为n'×1维向量，||rz-1||2为rz-1中n'个元素绝对值的2次方和的1/2次幂。

步骤c2中所述的r^d为n'×1维向量，||r^d||ξ为r^d中n'个元素绝对值的ξ次方和的1/ξ次幂。

本实施例中，所述超声波探伤装置1为a型数字超声探伤仪。实际使用时，也可以采用其它类型的超声波探伤设备。

步骤二中进行信号稀疏分解之前，按照常规信号稀疏分解时时频参数中尺度参数、位移参数、频率参数和相位参数的取值范围确定方法，对时频参数中尺度参数、位移参数、频率参数和相位参数的取值范围分别进行确定。步骤二中所述的为对待处理信号f(t)进行第z次迭代分解时的最佳匹配原子。

每个gabor原子均与其时频参数相对应，并且每个gabor原子均与其时频参数一一对应。

2006年4月《电子与信息学报》(第28卷第4期)公开的《利用fft实现基于mp的信号稀疏分解》(作者：尹忠科)一文中公开：“…，r＝(s,u,v,w)，时频参数可以按以下方法离散化：r＝(α^j,pα^jδu,kα^-jδv,iδw)，其中α＝2，δu＝1/2，δv＝π，δw＝π/6，0＜j＜log2n,0≤p≤n2^-j+1，0≤k≤n2^j+1，0≤i≤12。上面的描述就给出了一个具体的过完备原子库”。由上述内容可知，频率参数v根据kα^-jδv进行离散化，由于0≤k≤n2^j+1、0＜j＜log2n、α＝2且δv＝π，此时频率参数v的取值范围非常大，频率参数v的最小值为0且其最大值为因而频率参数v的取值范围为即使离散化，但频率参数v的取值范围仍非常大。

本实施例中，步骤c1中尺度参数的取值范围为[1,n']，位移参数的取值范围为[0,n']，频率参数的取值范围为相位参数的取值范围为[0,2π]。其中，fo为超声波探伤装置1的采样频率，fo的单位为mhz。

根据本领域公知常识，稀疏分解算法(也称为mp算法)存在两个缺陷，一是稀疏分解算法的计算量很大，计算时间在目前现有计算条件下十分巨大，无法进行实时检测；二是稀疏分解算法是连续条件下求的最优解，对于弱小缺陷的检测精度仍然有局限。

信号稀疏表示的目的就是在给定的超完备字典中用尽可能少的原子来表示信号，可以获得信号更为简洁的表示方式，从而使我们更容易地获取信号中所蕴含的信息，更方便进一步对信号进行加工处理，如压缩、编码等。信号稀疏表示方向的研究热点主要集中在稀疏分解算法、超完备原子字典(也称为原子库，gabor字典)和稀疏表示的应用等方面。信号稀疏表示的两大主要任务就是字典的生成和信号的稀疏分解。但现有的研究已证实在不增加原子库大小的情况下，从一个粗尺度到细尺度上在尺度和频率上搜索原子能显著提高mp算法(即匹配追踪算法)的性能。因而，对频率参数v的取值范围进行进一步的缩小化能有效提高mp算法(即匹配追踪算法)的性能。尤其对于频率参数而言，其中取值范围对mp算法(即匹配追踪算法)的性能影响更大。

由于频率参数(即频率参数v)的取值范围与信号的实际采样频率有关，在多年稀疏分解的研究经验基础上，同时借助时频参数的取值范围对提高mp算法(即匹配追踪算法)性能的影响进行充分、长期的研究与验证后，得出将频率参数的取值范围与被处理信号的实际采样频率(即超声波探伤装置1的采样频率fo)密切相关，并且并不是完全地一一对应关系，从简化稀疏分解算法的计算量与细化时频参数的取值范围并提高匹配追踪算法性能这一综合角度出发，将频率参数的取值范围限定为并且fo的单位为mhz，一方面能有效减少稀疏分解算法的计算量，实现进行实时检测；另一方面，有效提高mp算法(即匹配追踪算法)的性能，使得稀疏表示的信号能有效满足弱小缺陷的检测精度，达到更简易、更准确地获得信号中所蕴含的有效信息。通过将频率参数的取值范围限定为能进一步凸显信号中所蕴含有效的信息，使稀疏表示的信号更加侧重表达有效信息，弱化冗余信息，从而能更准确地表达信号本征特征，使信号提取精度能得到有效保证。

根据本领域公知常识，同时结合2006年4月《电子与信息学报》(第28卷第4期)公开的《利用fft实现基于mp的信号稀疏分解》(作者：尹忠科)一文可知，目前进行稀疏分解之前，通常均需要对时频参数的四个参数分别进行离散化，并生成过完备原子库，但所生成过完备原子库中的原子数量通常均非常大，所占用的存储空间非常大，并且计算量大，计算工程复杂，需对过完备原子库中的所有原子分别进行分析判断，并找出最佳匹配原子；同时，参数的取值范围与离散化方法对所生成的过完备原子库也会产生极大影响，不可避免会造成所生成过完备原子库(也称为过完备字典，gabor字典)的准确性较差，从而不能准确地表达信号本征特征，使信号提取精度不能得到保证。

而本发明中步骤二中进行信号稀疏分解之前，不需要生成过完备字典中的全部原子，只需采用数据处理设备2且调用寻优算法模块进行寻优后，便可简便、快速逐一找出最佳匹配原子，从而大大节省了存储空间。此外，所述寻优算法模块是在各参数的取值范围内(具体是在连续的空间中)进行最佳匹配原子的搜寻，与传统的匹配追踪算法在离散的搜索空间中(即经离散化获得的过完备字典或过完备原子库)进行最佳匹配原子的搜寻，因此寻优算法模块搜索的范围更广，搜索出的最佳匹配原子能更好地反映原始信号(即超声回波信号f(t))的特点，从而能进一步确保信号提取的精度，相应地能确保k个叠加信号的分离效果。

步骤c1中所述寻优算法模块为遗传算法模块、人工鱼群算法模块或人工蜂群算法模块。实际使用时，所述寻优算法模块也可以为其它类型的寻优算法模块。其中，调用遗传算法模块进行寻优时，采用常规的遗传算法即可；调用人工鱼群算法模块进行寻优时，采用常规的人工鱼群算法即可；调用遗传算法模块人工蜂群算法模块进行寻优时，采用常规的人工蜂群算法即可。

采用数据处理设备2且调用寻优算法模块进行寻优确定最佳匹配原子的方法，具有以下优点：第一、摆脱了傅里叶变换与小波变换等传统方法只能在正交基上进行分解的缺点，能够用更准确地表达信号本征特征，从而提高信号提取的精度；第二、能有效避免局部最优值的产生，并且可以进行连续空间的寻优查找，与原匹配追踪算法进行离散空间的寻优相比，扩大了搜索范围，从而进一步有效提高了信号提取的准确率；第三、通过寻优算法模块寻优找出最佳匹配原子，实现简便且提取速度快，能有效解决原始匹配算法复杂度大的问题，极大提高了降噪处理的收敛速度和信号提取的速度，提高了信号提取的实时性；第四、能有效提高信号提取的精度，解决了强噪声背景下的信号提取以及弱小信号的提取问题；第五、使用效果好，能解决超声无损探伤领域弱小缺陷等检测问题，提高了生产企业的产品质量，避免了安全隐患；第六、适用范围广，能有效适用于多种信号的提取过程，特别是能对非平稳难检声信号进行有效提取。因而，调用寻优算法模块进行寻优确定最佳匹配原子的方法，设计合理、效果好且实用价值高，不仅提高了信号提取的速度，而且经信号提取后能有效提高原信号的质量和性能指标，尤其在超声无损探伤中具有重要的作用。

实际使用过程中，不论是常规匹配追踪算法中建立过完备原子库后进行稀疏分解时，还是利用寻优算法模块寻优找出最佳匹配原子完成信号稀疏分解时，对最佳匹配原子进行确定，均以原子的适应度值(即gabor原子与对应残差量的内积)作为唯一的匹配判断标准，因而具有片面性。相应地，所确定的最佳匹配原子仅是适应度值最佳的gabor原子，然而对信号进行进行提取时，为选择与待处理信号f(t)最优匹配的原子，除考虑适应度值这一指标外，还需考虑信号稀疏度，从而能更准确地表达信号本征特征，使信号提取精度能得到有效保证。

而本发明步骤二中采用数据处理设备2对的时频参数rz进行查找进行查找时，先按照步骤c1进行时频参数寻优，找出适应度值最大的最优时频参数，并找出适应度值＞thrdz的多个较优时频参数，thrdz为适应度值判断阈值；再按照步骤c2，结合各较优时频参数的稀疏度进行最佳时频参数确定。因而，对最佳匹配原子进行确定时，不仅考虑到适应度值这一匹配判断标准，同时兼顾各原子的稀疏度(也称为鲁棒支撑)判断，具体是根据公式xsd＝||r^d||ξ计算得出第d个所述较优时频参数的稀疏度xsd，并根据计算结果(即各原子的鲁棒支撑)，选择稀疏度最小的所述较优时频参数作为最佳时频参数rz。因而，所确定的最佳匹配原子更准确，能更准确地表达信号本征特征。

为进一步确保信号提取的精度，确保分离出的每个叠加信号的准确性，需更加准确地选择最佳匹配原子。本实施例中，步骤三中进行回波信号分离之前，还需进行迭代分解最佳原子集合更新判断，以便选取更加准确的最佳匹配原子。

实际进行迭代分解最佳原子集合更新判断时，判断||r^k(t)||2是否小于ε：当||r^k(t)||2＜ε时，判断为无需对步骤二中所述迭代分解最佳原子集合进行更新，并进入步骤三；否则，当||r^k(t)||2≥ε时，判断为需对步骤二中所述迭代分解最佳原子集合进行更新，并进入步骤d；

其中，||r^k(t)||2为步骤二中所述r^k(t)的2-范数，ε为预先设定的残差量判断阈值；

步骤d、迭代分解最佳原子集合更新：采用数据处理设备2对步骤二中所述迭代分解最佳原子集合进行更新，包括以下步骤：

步骤d1、最佳匹配原子优化，过程如下：

步骤d11、原子随机选取：采用数据处理设备2从此时的所述迭代分解最佳原子集合中随机取出一个最佳匹配原子作为待优化原子，所述待优化原子记作其中j1为正整数且1≤j1≤k；

此时的所述迭代分解最佳原子集合中除所述待优化原子之外的k-1个最佳匹配原子均为待处理原子，k-1个待处理原子组成此时的待处理原子集合；

步骤d12、最佳匹配原子查找：所找出的最佳匹配原子记作的时频参数记作时频参数rj1'，时频参数rj1'＝(sj1',uj1',vj1',wj1')；

对最佳匹配原子进行查找时，根据预先设定的sj1'、uj1'、vj1'和wj1'的取值范围，采用数据处理设备2且调用所述寻优算法模块，找出使适应度值fitness(rj1')最大的最佳时频参数，所找出的最佳时频参数为时频参数rj1'；再根据公式求解出最佳匹配原子

其中，表示与的内积；ψ'0(t)为步骤d11中k-1个所述待处理原子之和；

步骤d13、原子替换判断、原子替换及迭代分解最佳原子集合更新：采用数据处理设备2且调用残值判断模块、适应度值判断模块或稀疏性判断模块，判断是否需对步骤d11中所述待优化原子进行替换，并根据判断结果对所述待优化原子进行替换；

采用数据处理设备2且调用残值判断模块判断是否需对步骤d11中所述待优化原子进行替换时，根据替换后残值||rj1'^k(t)||ξ是否小于替换前残值||rj1^k(t)||ξ进行判断：当||rj1'^k(t)||ξ＜||rj1^k(t)||ξ时，判断为需对步骤d11中所述待优化原子进行替换，将步骤d11中所述待优化原子替换为步骤d12中所述最佳匹配原子获得更新后的所述迭代分解最佳原子集合；否则，判断为无需对步骤d11中所述待优化原子进行替换，并进入步骤d15；

其中，rj1'^k(t)＝f(t)-ψj1'(t)，rj1^k(t)＝f(t)-ψj1(t)，ψj1(t)为本步骤中进行原子替换判断之前所述迭代分解最佳原子集合中k个最佳匹配原子之和；||rj1'^k(t)||ξ表示rj1'^k(t)的ξ-范数，||rj1^k(t)||ξ表示rj1^k(t)的ξ-范数，ξ为常数且0≤ξ≤1；

采用数据处理设备2且调用适应度值判断模块判断是否需对步骤d11中所述待优化原子进行替换时，根据替换后适应度值fitness(rj1')是否大于替换前适应度值fitness(rj1)进行判断：当fitness(rj1')＞fitness(rj1)时，判断为需对步骤d11中所述待优化原子进行替换，将步骤d11中所述待优化原子替换为步骤d12中所述最佳匹配原子获得更新后的所述迭代分解最佳原子集合；否则，判断为无需对步骤d11中所述待优化原子进行替换，并进入步骤d15；

其中，表示r^j1-1(t)与的内积，r^j1-1(t)＝f(t)-ψj1-1(t)，ψj1-1(t)为此时所述迭代分解最佳原子集合中前j1-1个最佳匹配原子之和；表示r^j1-1(t)与的内积；

采用数据处理设备2且调用稀疏性判断模块判断是否需对步骤d11中所述待优化原子进行替换时，根据||r^j1'||ξ是否小于||r^j1||ξ进行判断：当||r^j1'||ξ＜||r^j1||ξ时，判断为需对步骤d11中所述待优化原子进行替换，将步骤d11中所述待优化原子替换为步骤d12中所述最佳匹配原子获得更新后的所述迭代分解最佳原子集合；否则，判断为无需对步骤d11中所述待优化原子进行替换，并进入步骤d15；

其中，r^j1'为的残差量且r^j1为的残差量且||r^j1'||ξ表示r^j1'的ξ-范数，||r^j1||ξ表示r^j1的ξ-范数；

本步骤中，完成原子替换判断、原子替换及迭代分解最佳原子集合更新后，完成步骤d11中选取的一个所述最佳匹配原子的优化过程；

步骤d14、残差量判断：对步骤d13中所述最佳匹配原子优化后的残差量进行判断：当||r'j1^k(t)||2＜ε时，完成迭代分解最佳原子集合更新过程，进入步骤三；否则，当||r'j1^k(t)||2≥ε时，进入步骤d15；

其中，||r'j1^k(t)||2为r'j1^k(t)的2-范数；r'j1^k(t)为根据此时所述迭代分解最佳原子集合中的k个所述最佳匹配原子对f(t)进行k次迭代分解后的残差量；

步骤d15、下一个最佳匹配原子优化：按照步骤d11中至步骤d13中所述的方法，对此时所述迭代分解最佳原子集合中未进行优化的一个所述最佳匹配原子进行优化；

步骤d16、残差量判断：对步骤d15中所述最佳匹配原子优化后的残差量进行判断：当||r”j1^k(t)||2＜ε时，完成迭代分解最佳原子集合更新过程，进入步骤三；否则，当||r”j1^k(t)||2≥ε时，返回步骤d15；

其中，||r”j1^k(t)||2为r”j1^k(t)的2-范数；r”j1^k(t)为根据此时所述迭代分解最佳原子集合中的k个所述最佳匹配原子对f(t)进行k次迭代分解后的残差量。

其中，所述的r^k(t)为n'×1维向量，||r^k(t)||2为r^k(t)中n'个元素绝对值的2次方和的1/2次幂。

步骤d14中r'j1^k(t)为n'×1维向量，||r'j1^k(t)||2为r'j1^k(t)r'j^m(t)中n个元素绝对值的2次方和的1/2次幂。

步骤d16中r”j1^k(t)为n'×1维向量，||r”j1^k(t)||2为r”j1^k(t)中n'个元素绝对值的2次方和的1/2次幂。

步骤d13中所述的rj1'^k(t)为n'×1维向量，||rj1'^k(t)ξ为rj1'^k(t)中n'个元素绝对值的ξ次方和的1/ξ次幂。

所述的rj1^k(t)为n'×1维向量，||rj1^k(t)ξ为rj1^k(t)中n'个元素绝对值的ξ次方和的1/ξ次幂。

所述的r^j1'为n'×1维向量，||r^j1'||ξ为r^j1'中n'个元素绝对值的ξ次方和的1/ξ次幂。所述的r^j1为n'×1维向量，||r^j1||ξ为r^j1中n'个元素绝对值的ξ次方和的1/ξ次幂。

步骤d12中进行最佳匹配原子查找之前，先对sj1'、uj1'、vj1'和wj1'的取值范围分别进行设定，并且所设定的sj1'的取值范围与步骤c1中所设定尺度参数的取值范围相同，所设定的uj1'的取值范围与步骤c1中所设定位移参数的取值范围相同，所设定的vj1'的取值范围与步骤c1中所设定的频率参数的取值范围相同，所设定的wj1'的取值范围与步骤c1中所设定相位参数的取值范围相同。因此，sj1'的取值范围为[1,n']且sj1'∈[1,n']，uj1'的取值范围为[0,n']且uj1'∈[0,n']，vj1'的取值范围为且wj1'的取值范围为[0,2π]且wj1'∈[0,2π]。

本实施例中，sz的取值范围为[1,n']且sz∈[1,n']，uz的取值范围为[0,n']且uz∈[0,n']，vz的取值范围为且wz的取值范围为[0,2π]且wz∈[0,2π]。

实际使用过程中，不论是常规匹配追踪算法中建立过完备原子库后进行稀疏分解，还是利用寻优算法模块寻优找出最佳匹配原子完成信号稀疏分解，均具有一定的局限性，均是在一定限制条件下获取的最佳匹配原子，因而采用上述两种方法进行信号提取时，仅能说信号提取的准确性相对较高。其中，采用过完备原子库进行稀疏分解时，由于时频参数中各参数的取值范围与离散化方法均对所生成的过完备原子库也会产生极大影响，而最终确定的过完备原子库中不可能包括所有原子，不可避免会遗漏一个或多个最佳匹配原子，从而影响信号提取的准确性。而利用寻优算法模块寻优找出最佳匹配原子时，虽能提高信号提取速度，并能实现连续区间上的搜索，受寻优算法模块中算法本身的优劣性和性能影响，如搜索步长、搜索策略、搜索终止条件等，所找出的最佳匹配原子也仅是在一定程度上或一定范围内的最佳匹配原子，因而在一定程度上也会影响信号提取的准确性，相应地会影响最佳匹配原子的准确性，从而影响分离信号的准确度。

由上述内容可知，步骤二中完成信号稀疏分解后，还需进行残差量判断，判断此时的迭代分解最佳原子集合是否满足预先设定的信号提取精度要求，如不符合，需进入进行步骤d1进行最佳匹配原子优化，以便进一步提高最佳匹配原子的准确性。因而，步骤二中完成信号稀疏分解后，根据残差量判断结果，对步骤二中信号稀疏分解后迭代分解最佳原子集合是否满足预先设定的要求进行判断，增设这一信号提取精度(也称为原子选取精度)的验证环节，从而能进一步提高信号提取精度，从而确保分离信号的准确性。

步骤d1中进行最佳匹配原子优化时，所采用的最佳匹配原子优化方法设计合理、实现方便且使用效果好，从此时的迭代分解最佳原子集合中随机选取一个最佳匹配原子进行优化，当优化完成一个所述最佳匹配原子后，便通过残差量判断对此时迭代分解最佳原子集合是否满足信号提取精度要求进行判断，并根据判断结果决定是否需要继续对其余的最佳匹配原子进行优化。因而，实现简便，可随机选取一个最佳匹配原子进行优化，原子优化顺序不限，并且每完成一个最佳匹配原子的优化过程，便进行一次残差量判断，并能实现快速优化与优化结果实时判断相结合，能有效简化最佳匹配原子优化过程，并能快速达到最佳匹配原子优化的目的，有效提高信号分离精度。

对所述待优化原子进行优化时，所采用的与该待优化原子对应的最佳匹配原子查找方法(即步骤d12中进行最佳匹配原子查找的方法)设计合理，并且能简便、快速找出比该待优化原子更佳的最佳匹配原子。

所找出的时频参数rj1'为使适应度值fitness(rj1')最大的最佳时频参数；

由于且ψ'0(t)为步骤d11中k-1个所述待处理原子之和，因而为超声回波信号f(t)(即原始信号)减去除该待优化原子之外的k-1个所述待处理原子后的残差，因而是与该待优化原子直接相关的残差信号，因而利用作为评价所找出时频参数rj1'的指标更具有针对性，除去此时迭代分解最佳原子集合中k-1个所述待处理原子之外的残差信号与该待优化原子直接相关，利用寻优算法模块找出时频参数rj1'的过程不受其它原子(即k-1个所述待处理原子)的影响，并且找出比该待优化原子更佳的最佳匹配原子的概率更高，同时获得的最佳匹配原子对此时迭代分解最佳原子集合中k-1个所述待处理原子不会造成影响，无需重新进行信号稀疏分解，只需按照步骤d13中所述的方法完成该待优化原子的原子替换，使用效果非常好，能简便、快速提高信号提取精度，确保所选择原子的准确性。

步骤d13中进行原子替换判断及原子替换时，采用残值判断、适应度值判断或稀疏性判断任一种方法进行原子替换判断，可选取任一种方法进行原子替换判断，使用方式灵活，并且每一种原子替换判断方法均能实现有效的原子替换判断。

其中，调用残值判断模块判断是否需对步骤d11中所述待优化原子进行替换时，根据替换后残值||rj1'^k(t)||ξ是否小于替换前残值||rj1^k(t)||ξ的判断结果判断是否对待优化原子进行替换，选取残值较小的原子，使信号残差量更小，从而有效提高信号提取精度，使提取后信号更进一步逼近原始信号。

调用适应度值判断模块判断是否需对步骤d11中所述待优化原子进行替换时，根据替换后适应度值fitness(rj1')是否大于替换前适应度值fitness(rj1)判断是否对待优化原子进行替换，选取适应度值较大的原子，使信号残差量更小，从而有效提高信号提取精度，使提取后信号更进一步逼近原始信号。

调用稀疏性判断模块判断是否需对步骤d11中所述待优化原子进行替换时，根据最小鲁棒支撑来判断是否进行原子替换，选取更低的鲁棒支撑的原子，从而能更好的匹配信号特征，使得信号的表示更稀疏，从而达到有效提高信号提取精度的目的，使提取后信号更进一步逼近原始信号。

其中，r^j1'(ti)为r^j1'中第i个采样时刻的信号值，即r^j1'中的第i个信号值。

本实施例中，步骤d1中进行信号稀疏分解后，采用数据处理设备2将所述迭代分解最佳原子集合同步存储至数据存储器3内，所述数据存储器3与数据处理设备2连接。所述超声波探伤装置1、数据处理设备2和数据存储器3组成信号分离处理系统，详见图2。

步骤d13中进行原子替换判断及原子替换后，采用数据处理设备2对更新后的所述迭代分解最佳原子集合进行同步存储。

步骤d1中进行信号稀疏分解后，采用数据处理设备2将所述迭代分解最佳原子集合同步存储至数据存储器3内时，按照迭代分解先后顺序对所述迭代分解最佳原子集合中的k个所述最佳匹配原子分别进行存储；其中，为步骤d1中对所述超声回波信号f(t)进行第z次迭代分解时找出的最佳匹配原子。

本实施例中，步骤d1中最佳匹配原子优化时，按照存储先后顺序对所述迭代分解最佳原子集合中的所述最佳匹配原子进行优化；

步骤d1中最佳匹配原子优化时，最先进行优化的所述最佳匹配原子为步骤d1中所述迭代分解最佳原子集合中的第1个所述最佳匹配原子。

实际使用时，步骤d1中最佳匹配原子优化时，也可以不按存储先后顺序对所述迭代分解最佳原子集合中的所述最佳匹配原子进行优化。

所述的ε为预先设定的残差量判断阈值，实际使用时，可根据具体需要，对ε的取值大小进行限定。

本实施例中，所述的ε＝e^-5。

实际使用时，可根据具体需要，对ε的取值大小进行相应调整。

本实施例中，步骤c2中所述的ξ＝1。

实际使用时，可根据具体需要，对ξ的取值大小进行相应调整。

步骤二中进行信号稀疏分解时，采用数据处理设备2且按照步骤c1至步骤c2中所述的方法由先至后找出所述迭代分解最佳原子集合中的k个所述最佳匹配原子；

步骤c1中所述寻优算法模块为人工蜂群算法模块。

实际使用时，所述寻优算法模块也可以为其它寻优算法模块，如遗传算法模块、人工鱼群算法模块等。

本实施例中，步骤二中所述的为对所述超声回波信号f(t)进行第z次迭代分解时找出的最佳匹配原子。

步骤c1中所述寻优算法模块为人工蜂群算法模块；

步骤c1中进行时频参数寻优时，过程如下：

步骤2011、参数初始化：采用数据处理设备2设定所述人工蜂群算法模块的最大迭代次数mc、蜜源的数量sn、雇佣蜂的数量、观察蜂的数量和蜜源最大开采次数limit；同时，采用数据处理设备2随机生成sn个不同的蜜源，sn个蜜源均为待开采蜜源，所生成的sn个蜜源中第p个所述蜜源记作一个4维向量xp＝(x1p,x2p,x3p,x4p)，每个所述蜜源均为一个时频参数；所述雇佣蜂的数量和观察蜂的数量均为sn，所生成的每个蜜源均分配给一个雇佣蜂；

其中，p为正整数且p＝1、2、…、sn；x1p的取值范围与步骤201中预先设定的sn的取值范围相同，x2p的取值范围与步骤201中预先设定的un的取值范围相同，x3p的取值范围与步骤201中预先设定的vn的取值范围相同，x4p的取值范围与步骤201中预先设定的wn的取值范围相同；

本步骤中参数初始化完成后，采用数据处理设备2对sn个蜜源的适应度值分别进行计算，将sn个蜜源中所有适应度值＞thrdz的蜜源进行同步记录，并采用数据处理设备2将每个适应度值＞thrdz的蜜源均标记为所述较优时频参数；

步骤2012、雇佣蜂邻域搜索：每个雇佣蜂均对其所分配的蜜源进行邻域搜索，如果搜索到的新蜜源的适应度值大于原有蜜源的适应度值，将新蜜源作为雇佣蜂搜索到的待开采蜜源，并将已开采次数置0；否则，将原有蜜源的已开采次数加1；

本步骤中，每个雇佣蜂对其所分配的蜜源进行邻域搜索过程中，均采用数据处理设备2将搜索到的所有适应度值＞thrdz的新蜜源进行同步记录，并采用数据处理设备2将所记录的每个适应度值＞thrdz的蜜源均标记为所述较优时频参数；

步骤2013、观察蜂邻域搜索：根据步骤2012中雇佣蜂搜索到的所有蜜源的适应度值，计算得出雇佣蜂搜索到的各蜜源的被选择概率；观察蜂再按照计算得出的各蜜源的被选择概率，从雇佣蜂搜索到的所有蜜源中选择采蜜的蜜源作为新的蜜源；

观察蜂对所选择的蜜源进行邻域搜索，如果搜索到的新蜜源的适应度值大于原有蜜源的适应度值，观察蜂转变为雇佣蜂，将新的蜜源作为搜索到的待开采蜜源，并将已开采次数置0；否则，否则蜜源和雇佣蜂保持不变，将原有蜜源的已开采次数加1；

本步骤中，观察蜂对所选择的蜜源进行邻域搜索过程中，均采用数据处理设备2将搜索到的所有适应度值＞thrdz的新蜜源进行同步记录，并采用数据处理设备2将所记录的每个适应度值＞thrdz的蜜源均标记为所述较优时频参数；

步骤2014、最优蜜源实时记录：待雇佣蜂邻域搜索与观察蜂邻域搜索结束后，获得此时的最优蜜源并同步记录，所述人工蜂群算法模块的迭代次数加1；

雇佣蜂邻域搜索与观察蜂邻域搜索过程中，若蜜源的已开采次数达到蜜源最大开采次数limit，则观察蜂转变为侦查蜂，再通过侦查蜂生成新的蜜源，并将已开采次数置0；

步骤2015、多次重复步骤2012至步骤2014，直至所述人工蜂群算法模块的迭代次数达到最大迭代次数mc，此时获得的最优蜜源为步骤c1中所述最优时频参数；

步骤2011中对sn个蜜源的适应度值分别进行计算时、步骤2012中进行雇佣蜂邻域搜索时和步骤2013中进行观察蜂邻域搜索时，任一个蜜源的适应度值均为该蜜源所对应的gabor原子与r^z-1(t)的内积；

步骤2015中多次重复步骤2012至步骤2014后，获得步骤c1中所述的kz个所述较优时频参数；

步骤2011中、步骤2012中、步骤2013中和步骤2015中，采用数据处理设备2标记为所述较优时频参数的所有蜜源的总数量为步骤c1中所述的kz，每个采用数据处理设备2标记为所述较优时频参数的蜜源均为一个所述较优时频参数。

其中，r^z-1(t)为步骤二中对超声回波信号f(t)经过z-1次迭代分解后的残差量。

所述时频参数rz对应的gabor原子为

步骤2014中所获得的最优蜜源为一次迭代过程中获得的最优蜜源，步骤2015中获得的最优蜜源为mc次迭代过程中所获得的最优蜜源中适应度值最大的最优蜜源。

本实施例中，所述原有蜜源的适应度值表示r^z-1(t)与的内积；

本步骤中，雇佣蜂搜索到的待开采蜜源的数量为多个，雇佣蜂搜索到的所有待开采蜜源均为雇佣蜂搜索到的蜜源。

任一个搜索到的新蜜源的适应度值均为该蜜源所对应的gabor原子与r^z-1(t)的内积。

本实施例中，步骤d12中对最佳匹配原子进行查找时，采用数据处理设备2且调用所述寻优算法模块查找的时频参数rj1'，所述寻优算法模块为人工蜂群算法模块，过程如下：

步骤d21、参数初始化：采用数据处理设备2设定所述人工蜂群算法模块的最大迭代次数mc’、蜜源的数量sn’、雇佣蜂的数量、观察蜂的数量和蜜源最大开采次数limit’；同时，采用数据处理设备2随机生成sn’个不同的蜜源，sn’个蜜源均为待开采蜜源，所生成的sn’个蜜源中第p'个所述蜜源记作一个4维向量xp'＝(x1p',x2p',x3p',x4p')，每个所述蜜源均为一个时频参数；所述雇佣蜂的数量和观察蜂的数量均为sn’，所生成的每个蜜源均分配给一个雇佣蜂；

其中，p'为正整数且p'＝1、2、…、sn’；x1p'的取值范围与步骤二中预先设定的sz的取值范围相同，x2p'的取值范围与步骤二中预先设定的uz的取值范围相同，x3p'的取值范围与步骤二中预先设定的vz的取值范围相同，x4p'的取值范围与步骤二中预先设定的wz的取值范围相同；

步骤d22、雇佣蜂邻域搜索：每个雇佣蜂均对其所分配的蜜源进行邻域搜索，如果搜索到的新蜜源的适应度值大于原有蜜源的适应度值，将新蜜源作为雇佣蜂搜索到的待开采蜜源，并将已开采次数置0；否则，将原有蜜源的已开采次数加1；

步骤d23、观察蜂邻域搜索：根据步骤d22中雇佣蜂搜索到的所有蜜源的适应度值，计算得出雇佣蜂搜索到的各蜜源的被选择概率；观察蜂再按照计算得出的各蜜源的被选择概率，从雇佣蜂搜索到的所有蜜源中选择采蜜的蜜源作为新的蜜源；

观察蜂对所选择的蜜源进行邻域搜索，如果搜索到的新蜜源的适应度值大于原有蜜源的适应度值，观察蜂转变为雇佣蜂，将新的蜜源作为搜索到的待开采蜜源，并将已开采次数置0；否则，否则蜜源和雇佣蜂保持不变，将原有蜜源的已开采次数加1；

步骤d24、最优蜜源实时记录：待雇佣蜂邻域搜索与观察蜂邻域搜索结束后，获得此时的最优蜜源并同步记录，所述人工蜂群算法模块的迭代次数加1；

雇佣蜂邻域搜索与观察蜂邻域搜索过程中，若蜜源的已开采次数达到蜜源最大开采次数limit，则观察蜂转变为侦查蜂，再通过侦查蜂生成新的蜜源，并将已开采次数置0；

步骤d25、多次重复步骤d22至步骤d23，直至所述人工蜂群算法模块的迭代次数达到最大迭代次数mc，此时获得的最优蜜源为的时频参数rj1'，rj1'＝(sj1',uj1',vj1',wj1')；

步骤d22中进行雇佣蜂邻域搜索时和步骤d23中进行观察蜂邻域搜索时，任一个蜜源的适应度值均为该蜜源所对应的gabor原子与r^z-1(t)的内积。

本实施例中，x1p'的取值范围为[1,n']且x1p∈[1,n']，x2p'的取值范围为[0,n']且x2p∈[0,n']，x3p'的取值范围为且x4p的取值范围为[0,2π]且x4p'∈[0,2π]。

所述时频参数rj1'对应的gabor原子为

步骤d24中所获得的最优蜜源为一次迭代过程中获得的最优蜜源，步骤d25中获得的最优蜜源为mc’次迭代过程中所获得的最优蜜源中适应度值最大的最优蜜源。

并且，本发明中采用越界折回处理，雇佣蜂和观察蜂进行邻域搜索，生成新的蜜源后，对新蜜源进行边界检测，若超出上下界，则要对新蜜源进行越界折回操作。对新蜜源进行越界折回操作时，根据蜜源的四个元素的最大值和最小值，对新蜜源的4个元素分别进行越界折回操作。对新蜜源进行边界检测，根据蜜源的四个元素的最大值和最小值，对新蜜源的4个元素分别进行超界判断；并根据超界判断结果对新蜜源的4个元素分别进行越界折回操作，并获得越界折回操作后的蜜源，避免了误搜索现象。

其中，新蜜源对新蜜源的第q个元素进行超界判断时，当时，判断为未超界，无需对进行越界折回操作；当时，判断为超出下界，根据公式获得越界折回操作处理后的当时，判断为超出上界，根据公式获得越界折回操作处理后的

步骤2013中根据步骤2012中雇佣蜂搜索到的所有蜜源的适应度值，计算得出雇佣蜂搜索到的各蜜源的被选择概率时，根据轮盘赌方式计算各蜜源的被选择概率。其中，任一个蜜源的被选择概率均为该蜜源的适应度值与雇佣蜂搜索到的所有蜜源的适应度值之和的比值。步骤2013中观察蜂再按照计算得出的各蜜源的被选择概率，从雇佣蜂搜索到的所有蜜源中选择采蜜的蜜源作为新的蜜源时，选取被选择概率最大的蜜源作为新的蜜源。

相应地，步骤d23中根据步骤d22中雇佣蜂搜索到的所有蜜源的适应度值，计算得出雇佣蜂搜索到的各蜜源的被选择概率时，根据轮盘赌方式计算各蜜源的被选择概率。其中，任一个蜜源的被选择概率均为该蜜源的适应度值与雇佣蜂搜索到的所有蜜源的适应度值之和的比值。步骤d23中观察蜂再按照计算得出的各蜜源的被选择概率，从雇佣蜂搜索到的所有蜜源中选择采蜜的蜜源作为新的蜜源时，选取被选择概率最大的蜜源作为新的蜜源。

步骤2013中进行观察蜂邻域搜索时和步骤d23进行观察蜂邻域搜索时，为了加快搜索速度，搜索方式由随机搜索改为以下搜索方式：判断下一个随机搜索的蜜源的适应度值是否大于此时蜂群中心位置的蜜源的适应度值，当下一个随机搜索的蜜源的适应度值大于此时蜂群中心位置的蜜源的适应度值时，则将下一个随机搜索的蜜源作为新的蜜源；否则，将此时蜂群中心位置的蜜源作为新的蜜源，以提高算法的搜索速度。其中，此时蜂群中心位置的蜜源为此时所搜索到的所有蜜源之和的平均值。

由于随着蜂群搜索次数的增多，距离最优原子也越来越近，为了加快寻优速度并且避免陷入局部最优，则在观察蜂进行搜索时，比较下一搜索位置和蜜蜂中心位置的蜜源浓度(即适应度值)，并根据比较结果确定新的蜜源，这样加大了搜索的步长，加快蜜蜂向着最优原子方向而去的速度。

本实施例中，步骤2011中和步骤d21中进行参数初始化时，初始蜂群采用均匀分布方法生成。

原有人工蜂群算法中初始蜂群分布的随机性会造成搜索空间的不确定性，若初始蜂群搜索空间不包含全局最优解，且又不能在有限次的搜索内覆盖到全局最优解的区域，则会造成过早收敛现象。而本发明中初始蜂群采用均匀分布方法生成，因而能有效避免上述过早收敛现象。

本实施例中，步骤2011中进行参数初始化过程中，对sn个蜜源进行生成时，根据公式计算得出sn个蜜源中第p个所述蜜源的第q个元素xqp，其中q为正整数且q＝1、2、3或4；xqup为蜜源的第q个元素的最大值，xqlow为蜜源的第q个元素的最小值。

其中，蜜源的第1个元素的最大值为n'且其最小值为1，因而x1up＝n'且x1low＝1。蜜源的第2个元素的最大值为n'且其最小值为0，因而x2up＝n'且x2low＝0。蜜源的第3个元素的最大值为且其最小值为0，因而且x3low＝0。蜜源的第4个元素的最大值为2π且其最小值为0，因而x4up＝2π且x4low＝0。

本实施例中，步骤2012中雇佣蜂邻域搜索时，通过雇佣蜂在当前蜜源位置附近进行邻域搜索并产生一个新蜜源，新蜜源位置根据公式xp*＝xp+φp(xp-xl)进行确定，其中xp为当前所搜索的原蜜源，φp为[-1,1]范围内的一个随机数，xl为一个随机蜜源，xp*为新蜜源，通过φp限定了新蜜源的范围。

本实施例中，步骤d21中进行参数初始化过程中，对sn’(即sn')个蜜源进行生成时，根据公式计算得出sn’个蜜源中第p'个所述蜜源的第q个元素xqp'。

本实施例中，步骤d22中雇佣蜂邻域搜索时，通过雇佣蜂在当前蜜源位置附近进行邻域搜索并产生一个新蜜源，新蜜源位置根据公式xp'*＝xp'+φp(xp'*-xl)进行确定，其中xp'为当前所搜索的原蜜源，φp为[-1,1]范围内的一个随机数，xl为一个随机蜜源，xp'*为新蜜源，通过φp限定了新蜜源的范围。

本实施例中，步骤d14中所述的r'j1^k(t)根据公式进行计算，其中为此时所述迭代分解最佳原子集合中的第n1'个所述最佳匹配原子，n1'为正整数且n1'＝1、2、…、k；an1'为与根据此时所述迭代分解最佳原子集合中的前n1'-1个所述最佳匹配原子对f(t)进行n1'-1次迭代分解后残差量的展开系数；

步骤d16中所述的r”j1^k(t)根据公式进行计算，其中为此时所述迭代分解最佳原子集合中的第n2'个所述最佳匹配原子，n2'为正整数且n2'＝1、2、…、k；an2'为与根据此时所述迭代分解最佳原子集合中的前n2'-1个所述最佳匹配原子对f(t)进行n2'-1次迭代分解后残差量的展开系数。

本实施例中，步骤d13中所述的r^j1-1(t)为根据本步骤中进行原子替换判断之前所述迭代分解最佳原子集合中的前j1-1个所述最佳匹配原子对f(t)进行j1-1次迭代分解后的残差量。

步骤d13中对r^j1-1(t)进行计算时，根据本步骤中进行原子替换判断之前所述迭代分解最佳原子集合与进行计算，其中k'为正整数且k'＝1、2、…、j1-1，k'＜j1；为本步骤中进行原子替换判断之前所述迭代分解最佳原子集合中的第k'个所述最佳匹配原子，ak'为与根据本步骤中进行原子替换判断之前所述迭代分解最佳原子集合中的前k'-1个所述最佳匹配原子对f(t)进行k'-1次迭代分解后残差量的展开系数。

本实施例中，步骤三中回波信号分离完成后，还需对k个所述分离信号分别进行信号提取；k个所述分离信号的信号提取方法均相同。

对k个所述分离信号中任一个所述分离信号进行信号提取时，采用数据处理设备2对该分离信号进行信号提取，该分离信号为待处理信号且其记作信号f(t)；

采用数据处理设备2对信号f(t)进行信号提取时，过程如下：

步骤a1、基于寻优算法的稀疏分解：采用数据处理设备2且调用稀疏分解模块对信号f(t)进行迭代分解处理，将信号f(t)转换为并获得此时的迭代分解最佳原子集合；此时的所述迭代分解最佳原子集合中包含m个最佳匹配原子，为所述迭代分解最佳原子集合中的第n个所述最佳匹配原子；

式中r^m(t)为信号f(t)经过m次迭代分解后的残差量，其中m为预先设定的迭代分解总次数且m为正整数，n为正整数且n＝1、2、…、m；an为第n次迭代分解后的最佳匹配原子与上一次迭代分解后残差量的展开系数；

为第n次迭代分解时采用数据处理设备2且调用寻优算法模块找出的最佳匹配原子；为gabor原子且式中函数ψ(t)为高斯窗函数且rn为的时频参数，rn＝(sn,un,vn,wn)，sn为尺度参数，un为位移参数，vn为频率参数，wn为相位参数；

本步骤中，寻找最佳匹配原子时，根据预先设定的sn、un、vn和wn的取值范围，采用数据处理设备2且调用寻优算法模块，找出使适应度值fitness(rn)最大的最佳时频参数，所找出的最佳时频参数为时频参数rn；

其中，fitness(rn)为时频参数rn的适应度值，表示rⁿ-¹(t)与的内积；r^n-1(t)为信号f(t)经过n-1次迭代分解后的残差量，r⁰(t)＝f(t)；

步骤a2、残差量判断：判断||r^m(t)||2是否小于ε：当||r^m(t)||2＜ε时，进入步骤a4；否则，当||r^m(t)||2≥ε时，进入步骤a3；

其中，||rm(t)||2为步骤a1中所述r^m(t)的2-范数，ε为预先设定的残差量判断阈值；

步骤a3、最佳匹配原子优化，过程如下：

步骤a31、原子随机选取：采用数据处理设备2从此时的所述迭代分解最佳原子集合中随机取出一个最佳匹配原子作为待优化原子，所述待优化原子记作其中j为正整数且1≤j≤m；

此时的所述迭代分解最佳原子集合中除所述待优化原子之外的m-1个最佳匹配原子均为待处理原子，m-1个待处理原子组成此时的待处理原子集合；

步骤a32、最佳匹配原子查找：所找出的最佳匹配原子记作的时频参数记作时频参数rj'，时频参数rj'＝(sj',uj',vj',wj')；

对最佳匹配原子进行查找时，根据预先设定的sj'、uj'、vj'和wj'的取值范围，采用数据处理设备2且调用所述寻优算法模块，找出使适应度值fitness(rj')最大的最佳时频参数，所找出的最佳时频参数为时频参数rj'；再根据公式求解出最佳匹配原子

其中，表示与的内积；ψ0(t)为步骤a31中m-1个所述待处理原子之和；

步骤a33、原子替换判断及原子替换：采用数据处理设备2且调用残值判断模块、适应度值判断模块或稀疏性判断模块，判断是否需对步骤a31中所述待优化原子进行替换，并根据判断结果对所述待优化原子进行替换；

采用数据处理设备2且调用残值判断模块判断是否需对步骤a31中所述待优化原子进行替换时，根据替换后残值||rj'^m(t)||ξ是否小于替换前残值||rj^m(t)||ξ进行判断：当||rj'^m(t)||ξ＜||rj^m(t)||ξ时，判断为需对步骤a31中所述待优化原子进行替换，将步骤a31中所述待优化原子替换为步骤a32中所述最佳匹配原子获得更新后的所述迭代分解最佳原子集合；否则，判断为无需对步骤a31中所述待优化原子进行替换，并进入步骤a35；

其中，rj'^m(t)＝f(t)-ψj'(t)，rj^m(t)＝f(t)-ψj(t)，ψj(t)为本步骤中进行原子替换判断之前所述迭代分解最佳原子集合中m个最佳匹配原子之和；||rj'^m(t)||ξ表示rj'^m(t)的ξ-范数，||rj^m(t)||ξ表示rj^m(t)的ξ-范数，ξ为常数且0≤ξ≤1；

采用数据处理设备2且调用适应度值判断模块判断是否需对步骤a31中所述待优化原子进行替换时，根据替换后适应度值fitness(rj')是否大于替换前适应度值fitness(rj)进行判断：当fitness(rj')＞fitness(rj)时，判断为需对步骤a31中所述待优化原子进行替换，将步骤a31中所述待优化原子替换为步骤a32中所述最佳匹配原子获得更新后的所述迭代分解最佳原子集合；否则，判断为无需对步骤a31中所述待优化原子进行替换，并进入步骤a35；

其中，表示r^j-1(t)与的内积，r^j-1(t)＝f(t)-ψj-1(t)，ψj-1(t)为此时所述迭代分解最佳原子集合中前j-1个最佳匹配原子之和；表示r^j-1(t)与的内积；

采用数据处理设备2且调用稀疏性判断模块判断是否需对步骤a31中所述待优化原子进行替换时，根据||r^j'||ξ是否小于||r^j||ξ进行判断：当||r^j'||ξ＜||r^j||ξ时，判断为需对步骤a31中所述待优化原子进行替换，将步骤a31中所述待优化原子替换为步骤a32中所述最佳匹配原子获得更新后的所述迭代分解最佳原子集合；否则，判断为无需对步骤a31中所述待优化原子进行替换，并进入步骤a35；

其中，r^j'为的残差量且r^j为的残差量且r^j'||ξ表示r^j'的ξ-范数，||r^j||ξ表示r^j的ξ-范数；

本步骤中，完成原子替换判断及原子替换后，完成步骤a31中选取的一个所述最佳匹配原子的优化过程；

步骤a34、残差量判断：对步骤a33中所述最佳匹配原子优化后的残差量进行判断：当||r'j^m(t)||2＜ε时，进入步骤a4；否则，当||r'j^m(t)||2≥ε时，进入步骤a35；

其中，||r'j^m(t)||2为r'j^m(t)的2-范数；r'j^m(t)为根据此时所述迭代分解最佳原子集合中的m个所述最佳匹配原子对f(t)进行m次迭代分解后的残差量；

步骤a35、下一个最佳匹配原子优化：按照步骤a31中至步骤a33中所述的方法，对此时所述迭代分解最佳原子集合中未进行优化的一个所述最佳匹配原子进行优化；

步骤a36、残差量判断：对步骤a35中所述最佳匹配原子优化后的残差量进行判断：当||r”j^m(t)||2＜ε时，进入步骤a4；否则，当||r”j^m(t)||2≥ε时，返回步骤a35；

其中，||r”j^m(t)||2为r”j^m(t)的2-范数；r”j^m(t)为根据此时所述迭代分解最佳原子集合中的m个所述最佳匹配原子对f(t)进行m次迭代分解后的残差量；

步骤a4、信号重构：根据此时的所述迭代分解最佳原子集合，采用数据处理设备2得出信号f(t)的近似信号f″(t)；其中，近似信号f″(t)为从信号f(t)提取出的信号，

其中为此时所述迭代分解最佳原子集合中的第n'个所述最佳匹配原子，n'为正整数且n'＝1、2、…、m；an'为与根据此时所述迭代分解最佳原子集合中的前n'-1个所述最佳匹配原子对f(t)进行n'-1次迭代分解后残差量的展开系数。

步骤a2中所述的r^m(t)为n×1维向量，||r^m(t)||2为r^m(t)中n个元素绝对值的2次方和的1/2次幂。其中，n为正整数且其为信号f(t)的信号长度。其中，n＝n'。

步骤a34中r'j^m(t)为n×1维向量，||r'j^m(t)||2为r'j^m(t)中n个元素绝对值的2次方和的1/2次幂。

步骤a36中r”j^m(t)为n×1维向量，||r”j^m(t)||2为r”j^m(t)中n个元素绝对值的2次方和的1/2次幂。

步骤a33中所述的rj'^m(t)为n×1维向量，||rj'^m(t)||ξ为rj'^m(t)中n个元素绝对值的ξ次方和的1/ξ次幂。

所述的rj^m(t)为n×1维向量，||rj^m(t)||ξ为rj^m(t)中n个元素绝对值的ξ次方和的1/ξ次幂。

所述的r^j'为n×1维向量，||r^j'||ξ为r^j'中n个元素绝对值的ξ次方和的1/ξ次幂。所述的r^j为n×1维向量，||r^j||ξ为r^j中n个元素绝对值的ξ次方和的1/ξ次幂。

其中，步骤a1中

步骤a1中所述的为对待处理信号f(t)进行第n次迭代分解时的最佳匹配原子。

步骤a32中进行最佳匹配原子查找之前，先对sj'、uj'、vj'和wj'的取值范围分别进行设定，并且所设定的sj'的取值范围与步骤a1中所设定的sn的取值范围相同，所设定的uj'的取值范围与步骤a1中所设定的un的取值范围相同，所设定的vj'的取值范围与步骤a1中所设定的vn的取值范围相同，所设定的wj'的取值范围与步骤a1中所设定的wn的取值范围相同。

每个gabor原子均与其时频参数相对应，并且每个gabor原子均与其时频参数一一对应。

本实施例中，步骤a32中sj'的取值范围与sn的取值范围相同，uj'的取值范围与un的取值范围相同，vj'的取值范围与vn的取值范围相同，wj'的取值范围与wn的取值范围相同。因此，sj'的取值范围为[1,n]且sj'∈[1,n]，uj'的取值范围为[0,n]且uj'∈[0,n]，vj'的取值范围为且wj'的取值范围为[0,2π]且wj'∈[0,2π]。

步骤a1中完成信号稀疏分解后，还需进入步骤a2进行残差量判断，判断此时的迭代分解最佳原子集合是否满足预先设定的信号提取精度要求，如不符合，需进入进行步骤a3进行最佳匹配原子优化，以便进一步提高信号提取的准确性。因而，步骤a1中完成信号稀疏分解后，根据步骤a2中残差量判断结果，对步骤a1中信号稀疏分解后迭代分解最佳原子集合是否满足预先设定的信号提取精度要求进行判断，增设这一信号提取精度的验证环节，从而能进一步提高信号提取精度，使提取后信号进一步逼近原始信号。

步骤a3中进行最佳匹配原子优化时，所采用的最佳匹配原子优化方法设计合理、实现方便且使用效果好，从此时的迭代分解最佳原子集合中随机选取一个最佳匹配原子进行优化，当优化完成一个所述最佳匹配原子后，便通过残差量判断对此时迭代分解最佳原子集合是否满足信号提取精度要求进行判断，并根据判断结果决定是否需要继续对其余的最佳匹配原子进行优化。因而，实现简便，可随机选取一个最佳匹配原子进行优化，原子优化顺序不限，并且每完成一个最佳匹配原子的优化过程，便进行一次残差量判断，并能实现快速优化与优化结果实时判断相结合，能有效简化最佳匹配原子优化过程，并能快速达到最佳匹配原子优化的目的，有效提高信号提取精度。

对所述待优化原子进行优化时，所采用的与该待优化原子对应的最佳匹配原子查找方法(即步骤a32中进行最佳匹配原子查找的方法)设计合理，并且能简便、快速找出比该待优化原子更佳的最佳匹配原子。

所找出的时频参数rj'为使适应度值fitness(rj')最大的最佳时频参数；

由于且ψ0(t)为步骤a31中m-1个所述待处理原子之和，因而为待处理信号f(t)(即原始信号)减去除该待优化原子之外的m-1个所述待处理原子后的残差，因而是与该待优化原子直接相关的残差信号，因而利用作为评价所找出时频参数rj'的指标更具有针对性，除去此时迭代分解最佳原子集合中m-1个所述待处理原子之外的残差信号与该待优化原子直接相关，利用寻优算法模块找出时频参数rj'的过程不受其它原子(即m-1个所述待处理原子)的影响，并且找出比该待优化原子更佳的最佳匹配原子的概率更高，同时获得的最佳匹配原子对此时迭代分解最佳原子集合中m-1个所述待处理原子不会造成影响，无需重新进行信号稀疏分解，只需按照步骤a33中所述的方法完成该待优化原子的原子替换，最后直接进入步骤a4进行信号重构即可，因而使用效果非常好，能简便、快速提高信号提取精度。

步骤a33中进行原子替换判断及原子替换时，采用残值判断、适应度值判断或稀疏性判断任一种方法进行原子替换判断，可选取任一种方法进行原子替换判断，使用方式灵活，并且每一种原子替换判断方法均能实现有效的原子替换判断。

其中，调用残值判断模块判断是否需对步骤a31中所述待优化原子进行替换时，根据替换后残值||rj'^m(t)||ξ是否小于替换前残值||rj^m(t)||ξ的判断结果判断是否对待优化原子进行替换，选取残值较小的原子，使信号残差量更小，从而有效提高信号提取精度，使提取后信号更进一步逼近原始信号。

调用适应度值判断模块判断是否需对步骤a31中所述待优化原子进行替换时，根据替换后适应度值fitness(rj')是否大于替换前适应度值fitness(rj)判断是否对待优化原子进行替换，选取适应度值较大的原子，使信号残差量更小，从而有效提高信号提取精度，使提取后信号更进一步逼近原始信号。

调用稀疏性判断模块判断是否需对步骤a31中所述待优化原子进行替换时，根据最小鲁棒支撑来判断是否进行原子替换，选取更低的鲁棒支撑的原子，从而能更好的匹配信号特征，使得信号的表示更稀疏，从而达到有效提高信号提取精度的目的，使提取后信号更进一步逼近原始信号。

其中，r^j'(ti)为r^j'中第i个采样时刻的信号值，即r^j'中的第i个信号值。

本实施例中，步骤a1中进行信号稀疏分解后，采用数据处理设备2将所述迭代分解最佳原子集合同步存储至数据存储器3内，所述数据存储器3与数据处理设备2连接；

步骤a33中进行原子替换判断及原子替换后，采用数据处理设备2对更新后的所述迭代分解最佳原子集合进行同步存储。

其中，所述信号采样系统1、数据处理设备2与数据存储器3组成信号采样及提取系统。

步骤a35中所优化的所述最佳匹配原子为步骤a1中所述迭代分解最佳原子集合中的一个所述最佳匹配原子。已经完成优化的所述最佳匹配原子不能再次进行优化。

本实施例中，步骤a33中完成一个所述最佳匹配原子的优化过程后，将步骤a31中选取的所述最佳匹配原子标注为已优化原子。因而，步骤a35中所优化的最佳匹配原子为此时所述迭代分解最佳原子集合中除所述已优化原子之外的一个所述最佳匹配原子。其中，此时所述迭代分解最佳原子集合中未优化的一个所述最佳匹配原子为此时所述迭代分解最佳原子集合中除所述已优化原子之外的一个所述最佳匹配原子。

本实施例中，步骤a1中进行信号稀疏分解后，采用数据处理设备2将所述迭代分解最佳原子集合同步存储至数据存储器3内时，按照迭代分解先后顺序对所述迭代分解最佳原子集合中的m个所述最佳匹配原子分别进行存储；其中，为步骤a1中对所述待处理信号f(t)进行第n次迭代分解时找出的最佳匹配原子。

本实施例中，步骤a3中最佳匹配原子优化时，按照存储先后顺序对所述迭代分解最佳原子集合中的所述最佳匹配原子进行优化；

步骤a3中最佳匹配原子优化时，最先进行优化的所述最佳匹配原子为步骤a1中所述迭代分解最佳原子集合中的第1个所述最佳匹配原子。

实际使用时，步骤a3中最佳匹配原子优化时，也可以不按存储先后顺序对所述迭代分解最佳原子集合中的所述最佳匹配原子进行优化。

步骤a2中所述的ε为预先设定的残差量判断阈值，实际使用时，可根据具体需要，对ε的取值大小进行限定。

本实施例中，步骤a2中所述的ε＝e^-5。

实际使用时，可根据具体需要，对ε的取值大小进行相应调整。

本实施例中，步骤a33中所述的ξ＝1。

实际使用时，可根据具体需要，对ξ的取值大小进行相应调整。

本实施例中，步骤a1中所述的为步骤a1中对所述待处理信号f(t)进行第n次迭代分解时找出的最佳匹配原子；

步骤a1中进行信号稀疏分解时，采用数据处理设备(2)由先至后找出步骤a1中所述迭代分解最佳原子集合中的m个所述最佳匹配原子；

步骤a1中所述寻优算法模块为人工蜂群算法模块。实际使用时，所述寻优算法模块也可以为其它寻优算法模块，如遗传算法模块、人工鱼群算法模块等。

本实施例中，步骤a34中所述的r'j^m(t)根据公式进行计算，其中为此时所述迭代分解最佳原子集合中的第n1个所述最佳匹配原子，n1为正整数且n1＝1、2、…、m；an1为与根据此时所述迭代分解最佳原子集合中的前n1-1个所述最佳匹配原子对f(t)进行n1-1次迭代分解后残差量的展开系数；

步骤a36中所述的r”j^m(t)根据公式进行计算，其中为此时所述迭代分解最佳原子集合中的第n2个所述最佳匹配原子，n2为正整数且n2＝1、2、…、m；an2为与根据此时所述迭代分解最佳原子集合中的前n2-1个所述最佳匹配原子对f(t)进行n2-1次迭代分解后残差量的展开系数。

本实施例中，步骤a33中所述的r^j-1(t)为根据本步骤中进行原子替换判断之前所述迭代分解最佳原子集合中的前j-1个所述最佳匹配原子对f(t)进行j-1次迭代分解后的残差量。

步骤a33中对r^j-1(t)进行计算时，根据本步骤中进行原子替换判断之前所述迭代分解最佳原子集合与进行计算，其中k为正整数且k＝1、2、…、j-1，k＜j；为本步骤中进行原子替换判断之前所述迭代分解最佳原子集合中的第k个所述最佳匹配原子，ak为与根据本步骤中进行原子替换判断之前所述迭代分解最佳原子集合中的前k-1个所述最佳匹配原子对f(t)进行k-1次迭代分解后残差量的展开系数。

以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化，均仍属于本发明技术方案的保护范围内。