一种改进的小波阈值膝关节摆动信号去噪算法制造技术

本发明专利技术涉及生物医学信号领域，由于采集的膝关节摆动信号存在一些噪声的干扰，为了得到高质量的VAG信号，我们需要对信号做一些预处理，有利于后续的特征分析，这对计算机的辅助诊断和治疗是必要的。本发明专利技术根据已有的经验模态分解方法，提出一种改进的小波阈值膝关节摆动信号去噪算法，并对小波阈值去噪函数做了改进。实验仿真结果证明，本发明专利技术提出的算法和其他已有的算法相比，得到了较高的信噪和较小的均方根误差，使膝关节摆动信号的质量得到改善。

An improved wavelet threshold denoising algorithm for knee joint swing signals

This invention is involved in the field of biomedical signal. Because there are some noises in the swing signal of the knee joint, in order to get high quality VAG signal, we need to do some pretreatment to the signal, which is beneficial to the subsequent feature analysis, which is necessary for the computer aided diagnosis and treatment. Based on the existing empirical mode decomposition method, an improved wavelet threshold denoising algorithm for the knee joint swing signal is proposed, and the wavelet threshold de-noising function is improved. The experimental simulation results show that the proposed algorithm has higher signal to noise and smaller root mean square error compared with other existing algorithms, and improves the quality of the swing signal of the knee joint.

【技术实现步骤摘要】
一种改进的小波阈值膝关节摆动信号去噪算法
本专利技术涉及膝关节摆动信号领域，涉及到集成经验模态分解算法和小波阈值去噪算法，特别涉及到一种改进的小波阈值膝关节摆动信号去噪算法。
技术介绍
膝关节损伤和疾病的早期诊断和治疗，对于保障人口健康具有重要意义。目前临床上主要采用有创或微创方法来检测和评估膝关节的损伤与疾病情况，其中，MRI和CT等无创检测方法能够较为有效地实现对膝关节损伤疾病的检测，但不能提供早期诊断，且这些大型仪器设备成本高，不便于日常检查。作为微创手段的关节镜检查能提供诊断的详细信息，但其对膝关节的反复切口，会对患者造成痛苦，且不利于损伤疾病的恢复。另一方面，这些现行的检测方法在展现膝关节动态特征方面均存在一定的缺陷。因此，基于膝关节摆动(VAG)信号分析方法对膝关节损伤病变进行检测和辅助诊断是近年来兴起的一种新型无创检测方法。通过对膝关节摆动信号的特征进行分析，采用机器学习等分类方法，能够有效区分与病理学有关的正常与异常VAG信号，从而可使医生清楚地了解膝关节内部的情况，实现对膝关节损伤和病变的诊断与监测。由加速度传感器记录的膝关节摆动信号会受一些干扰，主要包括由低频移动产生的噪声，由医疗器械的热效应和半导体引起的随机噪声和膝关节运动时的肌肉收缩干扰。为了得到高质量的VAG信号，我们需要对信号做一些预处理，有利于后续的特征分析，这对计算节的辅助诊断和治疗是必要的。本专利技术根据已有的经验模态分解减少VAG信号随机噪声的方法，提出了一种改进的集成经验模态分解和小波阈值联合的去噪方法，并对小波阈值去噪函数做了改进，得到了较高的信噪比并减小了均方误差，使VAG信号的质量得到改善。
技术实现思路
由于采集的膝关节摆动信号存在一些噪声的干扰，为了得到高质量的VAG信号，需要对信号做一些预处理，有利于后续的特征分析，这对计算节的辅助诊断和治疗是必要的。本专利技术的主要目的就是为了得到质量更高的膝关节摆动信号，提出了一种改进的集成经验模态分解和小波阈值联合的去噪方法，并对小波阈值去噪函数做了改进，得到了较高的信噪比并减小了均方误差，使VAG信号的质量得到改善。为了达到上述目的，本专利技术的技术方案为：一种改进的小波阈值膝关节摆动信号去噪算法，包括以下步骤：第一步，对膝关节摆动信号进行集成经验模态分解(EEMD)第二步，对分解后得到的固有模态分量做去趋势波动分析(DFA)2.1)用DFA方法确定各阶模态分量的标度指数；2.2)通过相关分析确定k值。第三步，提出改进的小波阈值函数3.1)对各阶固有模态分量做小波变换，得到相应的小波变换系数；3.2)通过改进的小波阈值函数对小波系数做阈值计算，得到去除噪声后估计的小波系数；第四步，对第三步得到的小波阈值处理后的前k阶模态分量，和EEMD分解的后k阶模态分量进行重构，得到最终去噪后的膝关节摆动信号。与现有技术相比，本专利技术的有益效果为：本专利技术提出的改进的小波阈值函数和其他几种经典的小波阈值函数相比，有较好的性能，对仿真信号用不同的阈值函数对其去噪，本专利技术改进的小波阈值函数由最高的信噪比和最小的均方误差。本专利技术提出的去噪方法用在膝关节摆动信号上，可以增加膝关节摆动信号的信噪比和减小均方误差，并改善膝关节摆动信号的质量。附图说明图1原始信号和集经验模态分解的前3阶的固有模态分量(IMF)；图2是集成经验模态EMD分解后的IMF4-IMF8分量；图3是集成经验模态EMD分解后的IMF9-IMF12及剩余分量；图4是不同小波阈值函数的波形图；图5使用本专利技术提出的方法和其他方法对膝关节摆动信号去噪后的图；(a)原始的含噪声VAG信号，(b)使用本文方法对VAG信号去噪，(c)EMD对VAG信号去噪。具体实施方式为使本专利技术实施的目的、技术方案和优点更加清晰，以下结合本专利技术的技术方案和附图作进一步详细描述：一种改进的小波阈值膝关节摆动信号去噪算法，该方法主要可分为四个环节，分别为对VAG信号做EEMD分解；对分解后得到的固有模态分量做去趋势波动分析(DFA)；提出改进的小波阈值函数，采用这个阈值函数对各阶模态分量做小波阈值去噪；对第三步得到的小波阈值处理后的前k阶模态分量和EEMD分解的后k阶模态分量进行重构，得到最终去噪后的膝关节摆动信号。其具体步骤如下：步骤A.对原始的膝关节摆动信号进行集成经验模态分解(EEMD)，得到各阶固有模态分量，包括如下子步骤：A1.在第n次试验中，信号x(t)加入一个白噪声系列un(t)，获得一个新的时间序列：yn(t)＝x(t)+un(t),n＝1,2,…,N式中，N表示集成的数量。A2.通过三次样条函数获得yn(t)的极大值包络线un(t)和极小值包络线vn(t)，其均值mn(t)表示为：A3.令hn(t)＝yn(t)-mn(t)，如果hn(t)不满足IMF条件(a)和(b)，令hn(t)＝yn(t)，重复第二个步骤，直到hn(t)为IMF函数；(a)信号的零点和极值点个数相同或相差一个；(b)信号的上下包络线关于时间轴对称；A4.将信号作为yn(t)重复以上两个步骤，直到rn(t)不能继续分解，原信号分解为：式中，i表示每次分解过程中IMFs的个数，表示第j个IMF，rin表示yn(t)在第n次分解过程的剩余量；A5.重复上述A1～A4四个步骤n次，每次给信号加入不同的白噪声；A6.对n次分解过程中的第j个IMF取平均值，得到的值就是EEMD最终的固有模态分量IMF；步骤B.对分解后得到的固有模态分量做去趋势波动分析(DFA)采用去趋势波动分析方法确定各阶固有模态分量的标度指数，通过相关分析确定k值，使得C1-Ck为噪声主导分量，k之后的IMF为信号主导分量。DFA的具体算法如下：B1.加入给定长度为L分解后的IMF其平均值为ωj，完整的IMF时间序列s(m)定义为：B2.完整的IMF时间序列s(m)被分成几个等长且长度为k1的窗口段，最小二乘拟合得到的趋势定义为sk1(m)，并适应所选择的窗口样本，局部线性趋势波动的计算是在每一个窗口段中通过从完整的时序信号s(m)减去局部线性趋势波动sk1(m)，平均波动F(k1)是用局部线性趋势波动的均方根来表示：采用分形标度指数α定义log10F(k1)和log10k1的线性关系，可以表示为F(k1)～k1α。当0.5<α<1时，时间序列具有长程相关性；当0<α<0.5时，时间序列具有反相关性，此时的时间序列被认为是噪声。因为随机噪声不是远程相关的，膝关节摆动信号可以用远程相关的IMFs(α>0.5)来构造重建，而具有反相关的IMFs(α>0.5)和最终的单调余量(基线漂移)被认为是噪声。所以分解点k确定为6。步骤C.提出改进的小波阈值函数，采用这个阈值函数对各阶模态分量做小波阈值去噪，具体步骤如下：C1.对前k阶模态分量做小波变换，得到相应的小波变换系数；C2.通过改进的小波阈值函数对小波系数做阈值计算，得到去除噪声后估计的小波系数；假设dj,k2为小波变换后的系数，d'j,k2为做阈值处理后的小波系数，j为小波分解的尺度，k2∈Z，t1为阈值，t1取为σ为噪声的标准差，可以由每一层小波系数估计，其中，表示第二层小波系数本文档来自技高网...

【技术保护点】
1.一种改进的小波阈值膝关节摆动信号去噪算法，其特征在于以下步骤：步骤A.对原始的膝关节摆动信号进行集成经验模态分解EEMD，得到各阶固有模态分量IMF，如下式所示：

【技术特征摘要】
1.一种改进的小波阈值膝关节摆动信号去噪算法，其特征在于以下步骤：步骤A.对原始的膝关节摆动信号进行集成经验模态分解EEMD，得到各阶固有模态分量IMF，如下式所示：其中：N表示集成的数量；n表示信号分解过程的次数；表示第j个固有模态分量；步骤B.对分解后得到的固有模态分量做去趋势波动分析DFA采用去趋势波动分析方法确定各阶固有模态分量的标度指数，通过相关分析确定k值，使得C1-Ck为噪声主导分量，k之后的IMF为信号主导分量；步骤C.提出改进的小波阈值函数，采用这个阈值函数对各阶模态分量做小波阈值去噪，具体步骤如下：C1.对前k阶模态分量做小波变换，得到相应的小波变换系数；C2.通过改进的小波阈值函数对小波系数做阈值计算，得到去除噪声后估计的小波系数；假设dj,k2为小波变换后的系数，d'j,k2为做阈值处理后的小波系数，k2∈Z，t1为阈值，小波阈值函数如下：其中，u＝1-exp[-m(|dj,k2|-t1)2]，u的范围是0≤u≤1；β为可调参数，令β的范围为[0,1]；C3.得到前k阶模态分量估计的小波系数，将这些估计的小波系数相加，得到一个小波阈值处理后的前k阶模态分量；步骤D.对第三步得到的小波阈值处理后的前k阶模态分量，和EEMD分解的后k阶模态分量进行重构，得到最终去噪后的膝关节摆动信号。2.根据权利要求1所述的一种改进的小波阈值膝关节摆动信号去噪算法，其特征在于，所述的步骤A包括如下子步骤：A1.在第n次试验中，信号x(t)加入一个白噪声系列un(t)，获得一个新的时间序列：yn(t)＝x(t)+un(t),n＝1,2,…,N式中，N表示集成的数量；A2.通过三次样条函数获得yn(t)的极大值包络线un(t)和极小值包络线vn(t)，其均值mn(t)表示为：A3.令hn(t)＝yn(t)-mn(t)，如果hn(t)不满足IMF条件(a)和(b)，令hn(t)＝yn(t)，重复第二个步骤，直到hn(t)为IMF函数；(a)信号的零点和极值点个数相同或相差一个；(b)信号的上下包络线...

【专利技术属性】
技术研发人员：邱天爽，徐一平，张耘歌，汪靖华，于泓颖，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：辽宁,21