基于分类匹配和分数阶扩散的心电信号去噪方法和系统技术方案

本发明专利技术提供一种基于分类匹配和分数阶扩散的心电信号去噪方法，包括以下步骤：S1.输入待处理的心电信号；S2.对心电信号进行滑动分类匹配，提取每个参考信号片段的相似数据矩阵；S3.通过变换系数的收缩对每个相似数据矩阵进行去噪；S4.将各去噪后相似数据矩阵中的信号片段返回原始位置，通过加权求和获得初步估计信号；S5.将获得的初步估计信号通过分数阶扩散进一步去噪，获得最终的去噪后的心电信号。本发明专利技术能够在去除心电信号高频噪声的同时，进一步去除心电信号中的同频段低频噪声，并保护心电信号中的细节特征信息。

ECG signal denoising method and system based on classification matching and fractional order diffusion

The invention provides a denoising method of ECG signal based on classification matching and fractional order diffusion, which includes the following steps: S1. input the ECG signal to be processed; S2. slide classification matching of ECG signal to extract the similar data matrix of each reference signal segment; S3. denoising each similar data matrix by contraction of transformation coefficient; S4. denoising the similar data moments after denoising. The signal fragments in the array return to the original position, and the preliminary estimated signal is obtained by weighted summation; S5. The preliminary estimated signal is further denoised by fractional diffusion to obtain the final denoised ECG signal. The invention can further remove the same frequency band low frequency noise in the ECG signal while removing the high frequency noise of the ECG signal, and protect the detailed feature information in the ECG signal.

【技术实现步骤摘要】
基于分类匹配和分数阶扩散的心电信号去噪方法和系统
本专利技术涉及信号噪声处理
，具体而言涉及一种基于分类匹配和分数阶扩散的心电信号去噪方法和系统。
技术介绍
随着临床医学的发展，心电信号已经成为了一个用于医疗诊断的重要工具。心电信号在提取过程中，往往会受到噪声的干扰，为了不影响进一步分析诊断，心电信号的去噪成为了必要的处理步骤之一。现有的心电信号去噪方法有很多种，包括小波变换方法，经验模态法等。它们在去除心电信号的高频噪声方面都取得了较好的效果，但是它们都有一个相同的缺点，它们无法有效地去除有心电信号同频率的低频段噪声。随着偏微分方程理论的发展，偏微分方程逐渐成为一种新型的信号处理工具，基于偏微分方程的信号处理方法将离散信号表示成连续的数学模型，利用偏微分方程完善的数值分析理论对信号进行处理。经典非线性扩散滤波就是基于偏微分方程的信号处理方法，但是它只考虑了上一时刻的扩散结果，并未将所有历史时刻作为参考，存在一定的局限性。
技术实现思路
本专利技术目的在于提供一种基于分类匹配和分数阶扩散的心电信号去噪方法和系统，在去除心电信号高频噪声的同时，进一步去除心电信号中的同频段低频噪声，并保护心电信号中的细节特征信息。为达成上述目的，本专利技术提出一种基于分类匹配和分数阶扩散的心电信号去噪方法，该方法包括以下步骤：S1.输入待处理的心电信号；S2.对心电信号进行滑动分类匹配，提取每个参考信号片段的相似数据矩阵；S3.通过变换系数的收缩对每个相似数据矩阵进行去噪；S4.将各去噪后相似数据矩阵中的信号片段返回原始位置，通过加权求和获得初步估计信号；S5.将获得的初步估计信号通过分数阶扩散进一步去噪，获得最终的去噪后心电信号。进一步的，在所述步骤S2中，对信号进行滑动分类匹配，提取每个参考片段的相似数据矩阵的步骤为：S21.选择待处理的心电信号的某一参考片段并在一大小固定为n的搜索窗口内搜索心电信号中与参考片段相似的信号片段，将所有相似的信号片段堆叠成一相似数据矩阵；将待处理的心电信号的某一参考片段设定成(空间域)，生成的相似数据矩阵设定为优选的，生成的相似数据矩阵为二维数据矩阵。步骤S21可以进一步表述为：选择待处理的心电信号的某一参考片段(空间域)并在大小固定为n的搜索窗口内搜索与参考片段相似的信号片段，在搜索结束后，将这些相似片段堆叠成二维数据矩阵S22.通过向后移位p来选择下一个参考片段，依旧在搜索窗口内搜索相似片段并堆叠形成下一个相似数据矩阵；S23.重复上述步骤，提取出每个参考片段对应的相似数据矩阵；所述n、p均大于零。进一步的，在所述步骤S3中，通过变换系数的收缩对每个相似数据矩阵进行去噪的步骤为：S31.对所述相似数据矩阵进行二维离散余弦变换；S32.将变换后的系数进行阈值收缩；S33.对阈值收缩后的系数进行离散余弦反变换得到去噪后的相似片段。以前述二维数据矩阵为例，代入具体的公式，通过变换系数的收缩对每个相似数据矩阵进行去噪的具体步骤为：S31.对得到的相似数据矩阵进行二维离散余弦变换S32.将变换后的系数进行阈值收缩S33.对阈值收缩后的系数进行离散余弦反变换得到去噪后的相似片段进一步的，在所述步骤S4中，将各去噪后相似数据矩阵中的信号片段返回原始位置，通过加权求和获得初步估计信号包括：处理所有参考信号段，获取每个参考信号片段位于心电信号不同位置坐标处的相似信号片段的估计集合。其中，如果心电信号的同一位置坐标上存在不同的局部片段估计值，为该位置坐标上参考信号不同的所有估计值设置与对应参考信号相关的权重，将所有局部估计的加权平均之和作为初步的估计信号。具体的，所述步骤S4包括：在处理完所有参考信号段后，可以得到以为参考信号片段且位于不同位置处的相似信号片段的估计集合需要注意一点，对于同一位置坐标，可能存在不同的局部片段估计值。例如，和都是xb处的估计值，但是它们是由不同的参考信号得到的。对位于同一位置但参考信号不同的所有估计值来说，它们分别有着与参考信号相关的权重初步的估计信号可以通过来计算获得，它是所有局部估计的加权平均之和。进一步的，在所述步骤S5中，将获得的初步估计信号通过分数阶扩散进一步去噪，获得最终的去噪后心电信号的步骤为：将初步的估计信号通过时间分数阶扩散模型进行扩散滤波，将时间分数阶扩散方程的数值解作为最终的滤波结果。优选的，利用有限差分法对所述时间分数阶扩散方程进行数值算法，所求的数值解作为最终的去噪后的心电信号具体的，将初步的估计信号通过时间分数阶扩散方程进行扩散滤波，方程的数值解就是最终的滤波结果。时间分数阶扩散方程设定如下：其中，α是时间分数阶导数的阶次，f(x)是在区间[0,L]上的初始信号，为扩散函数。利用有限差分法对时间分数阶扩散方程进行数值算法：对于α阶的时间分数阶Caputo导数定义如下：它的离散形式为:其中，并且则其隐式差分算法为：假定将上式写成矩阵格式为：U1＝[B0]-1U0其中，通过上述数值算法对时间分数阶扩散方程进行求解，所求的数值解即为扩散后的结果，也就是最终的去噪后的心电信号。基于前述方法，本专利技术还提及一种基于分类匹配和分数阶扩散的心电信号去噪系统，包括以下模块：1)用以输入待处理的心电信号的模块；2)用以对心电信号进行滑动分类匹配，提取每个参考信号片段的相似数据矩阵的模块；3)用以通过变换系数的收缩对每个相似数据矩阵进行去噪的模块；4)用以将各去噪后相似数据矩阵中的信号片段返回原始位置，通过加权求和获得初步估计信号的模块；5)用以将获得的初步估计信号通过分数阶扩散进一步去噪，获得最终的去噪后的心电信号的模块。由以上本专利技术的技术方案，与现有相比，其显著的有益效果在于，本专利技术提供的方法将分类匹配与分数阶扩散相结合，通过分类匹配去除了心电信号中的高频噪声；在利用分数阶扩散进一步去除心电信号中低频噪声的同时，更好的保护了心电信号的细节特征，有利于下一步的分析诊断。应当理解，前述构思以及在下面更加详细地描述的额外构思的所有组合只要在这样的构思不相互矛盾的情况下都可以被视为本公开的专利技术主题的一部分。另外，所要求保护的主题的所有组合都被视为本公开的专利技术主题的一部分。结合附图从下面的描述中可以更加全面地理解本专利技术教导的前述和其他方面、实施例和特征。本专利技术的其他附加方面例如示例性实施方式的特征和/或有益效果将在下面的描述中显见，或通过根据本专利技术教导的具体实施方式的实践中得知。附图说明附图不意在按比例绘制。在附图中，在各个图中示出的每个相同或近似相同的组成部分可以用相同的标号表示。为了清晰起见，在每个图中，并非每个组成部分均被标记。现在，将通过例子并参考附图来描述本专利技术的各个方面的实施例，其中：图1是本专利技术的基于分类匹配和分数阶扩散的心电信号去噪方法的流程图。具体实施方式为了更了解本专利技术的
技术实现思路
，特举具体实施例并配合所附图式说明如下。在本公开中参照附图来描述本专利技术的各方面，附图中示出了许多说明的实施例。本公开的实施例不必定意在包括本专利技术的所有方面。应当理解，上面介绍的多种构思和实施例，以及下面更加详细地描述的那些构思和实施方式可以以很多方式中任意一种来实施，这是因为本专利技术所公开的构思和实施例并不限于任何实施方式。另外，本专利技术公开的一些方面可以单独使用，或者与本专利技术本文档来自技高网...

【技术保护点】
1.一种基于分类匹配和分数阶扩散的心电信号去噪方法，其特征在于，所述方法包括以下步骤：S1.输入待处理的心电信号；S2.对心电信号进行滑动分类匹配，提取每个参考信号片段的相似数据矩阵；S3.通过变换系数的收缩对每个相似数据矩阵进行去噪；S4.将各去噪后相似数据矩阵中的信号片段返回原始位置，通过加权求和获得初步估计信号；S5.将获得的初步估计信号通过分数阶扩散进一步去噪，获得最终的去噪后的心电信号。

【技术特征摘要】
1.一种基于分类匹配和分数阶扩散的心电信号去噪方法，其特征在于，所述方法包括以下步骤：S1.输入待处理的心电信号；S2.对心电信号进行滑动分类匹配，提取每个参考信号片段的相似数据矩阵；S3.通过变换系数的收缩对每个相似数据矩阵进行去噪；S4.将各去噪后相似数据矩阵中的信号片段返回原始位置，通过加权求和获得初步估计信号；S5.将获得的初步估计信号通过分数阶扩散进一步去噪，获得最终的去噪后的心电信号。2.根据权利要求1所述的基于分类匹配和分数阶扩散的心电信号去噪方法，其特征在于，步骤S2中，对心电信号进行滑动分类匹配，提取每个参考信号片段的相似数据矩阵的步骤包括：S21.选择待处理的心电信号的某一参考片段并在一大小固定为n的搜索窗口内搜索心电信号中与参考片段相似的信号片段，将所有相似的信号片段堆叠成一相似数据矩阵；S22.通过向后移位p来选择下一个参考信号片段，在搜索窗口内搜索相似的信号片段并堆叠形成下一个相似数据矩阵；S23.重复上述步骤S22，提取每个参考信号片段对应的相似数据矩阵。3.根据权利要求2所述的基于分类匹配和分数阶扩散的心电信号去噪方法，其特征在于，所述相似数据矩阵为二维数据矩阵。4.根据权利要求1所述的基于分类匹配和分数阶扩散的心电信号去噪方法，其特征在于，在所述步骤S3中，通过变换系数的收缩对每个相似数据矩阵进行去噪的步骤包括：S31.对所述相似数据矩阵进行二维离散余弦变换；S32.将变换后的系数进行阈值收缩；S33.对阈值收缩后的系数进行离散余弦反变换得到去噪后的相似片段。5.根据权利要求1所述的基于分类匹...

【专利技术属性】
技术研发人员：李远禄，蒋民，赵伟静，李俊，周慧敏，孙双龙，
申请(专利权)人：南京信息工程大学，
类型：发明
国别省市：江苏,32