专利名称：一种脉搏波信号的处理方法及处理系统的制作方法
一种脉搏波信号的处理方法及处理系统技术领域
本发明属于信号处理领域，尤其涉及一种脉搏波信号的处理方法及处理系统。
背景技术：
脉搏波信号分析，是针对脉搏波图的分析和处理。目前大多数分析方法主要从时域或频域中提取特征信息。时域分析法是在时间方向上分析脉搏波信号的动态特征，而频域分析方法是通过傅立叶变换将时间信号以频谱的形式表示出来，从频域上去研究一个复杂信号的特征。主要的脉搏分析方法包括如下:
1、高斯函数法，脉搏波在形式上可看成由主波、重搏波和重搏前波等三者叠加而成，并且这三波的位置、形状的变化和生理病理状态密切相关。一个周期的脉搏波可三个高斯函数来合成，分别称这三个高斯函数为钟形主波，钟形重搏波和钟形重搏前波。求出高斯函数的参数来反映脉搏波的特征。
缺点:计算繁琐，拟和求出的特征参数与真实的生理过程存在差异，和生理意义关系不明确。
2、脉图面积法，由脉搏波波图面积提取的特征量K值，从宏观上描述出脉搏波平均特征。
缺点:反映出脉搏曲线比较粗糙易造成误差。
3、频域法，采用频域方法研究脉搏波，就是把时域的脉搏波信号通过离散快速傅立叶变换到频域，从脉搏波频谱中提取出与人体心血管生理病理相应的信息。把一个很复杂的由许多重叠波构成的脉搏波分解成不同的谐波，这样就能把其中所包含的丰富信息和能量提取出来，因此较时域分析法能更清楚地反映脉搏波。
缺点:计算复杂，结果抽象。
4、特征点法，脉搏波波图的特征点(或称标志点)实质上就是脉搏波压力曲线的拐点，它是心动周期中从一个力学过程转变为另一个力学过程的转变点，因而这些拐点都有明确的生理意义。通过对不同脉图特征点参数的差异来寻求辨别脉图的客观标准。
缺点:特征点不明显、难以辨认带有误差。发明内容
本发明的目的在于提供一种脉搏波信号的处理方法，旨在解决计算繁琐复杂、结构抽象、反映出脉搏曲线比较粗糙易造成误差、特征点不明显、难以辨认带有误差的问题。
本发明是这样 实现的，一种脉搏波信号的处理方法，所述方法包括以下步骤:
A、采集并显示脉搏波信号；
B、采用数字滤波、周期识别、基线调整、特征点识别的方法将采集的脉搏波信号进行处理；
C、将处理过的脉搏波信号进行存储；
所述数字滤波，通过一定的算法，对原始数据进行平滑、滤波预处理，去除波形中毛刺，最大限度的保证脉搏波在在幅度及相位上的不失真；
所述周期识别，即一段多周期的脉搏信号中辨识出每个周期的起始点坐标，将脉搏波形序列分割成一个一个的单周期波形数据；
所述基线调整，即有效的去除脉搏波形的基线漂移，消除呼吸运动、肌肉运动或身体位移的因素影响；
所述特征点识别，即标识出单周期内脉搏波形的主波峰值点、降中峡和重搏波。
本发明的进一步技术方案是:所述B步骤中包括以下步骤:
B1、利用三点均值滤波算法对原始数据进行平滑、滤波预处理，去除波形中毛刺，最大限度的保证脉搏波在幅度及相位上的不失真；
B2、对一段多周期脉搏波信号进行求导寻找最大值，设置阈值找各个周期峰值点，寻找各个压力上升最快之间的最小值，此最小值点就是一个脉搏波周期的起始点；
B3、利用插值拟合法拟合基线漂移曲线，用原始信号减去漂移曲线得到消除基线漂移的脉搏信号；
B4、标识出单周期内脉搏波形的主波峰值点、降中峡和重搏波的位置。
本发明的进一步技术方案是:所述B2步骤中还包括以下步骤:
B21、对一段·多周期的脉搏波信号进行一阶求导；
B22、寻找求导结果中的最大值点，设置阈值；
B23、求导结果与阈值比较寻找各周期的峰值点；
B24、寻找脉搏周期的起始点。
本发明的进一步技术方案是:所述步骤B3还包括以下步骤:
B31、确定单个脉搏波的起始点和结束点的位置和幅值；
B32、利用插值拟合法拟合基线漂移曲线；
B33、用原始信号减去漂移曲线得到消除基线漂移的脉搏信号。
本发明的进一步技术方案是:所述步骤B4还包括以下步骤:
B41、求脉搏波信号在起始点与结束点之间的主波峰值点；
B42、计算单周期内脉搏波形数据的二阶导数及单周期内1/2心动周期位置；
B43、对二阶导数结果进行峰值检测，得到波峰波谷数组；
B44、取波峰数组中靠近1/2心动周期的波峰点为降中峡位置；
B45、取波谷数组中靠近降中峡之后的第一个波谷点位置为重搏波位置。
本发明还提供一种脉搏波信号的处理系统，该处理系统包括:
采集显示模块，用于采集并显示脉搏波信号；
信号处理模块，用于采用数字滤波、周期识别、基线调整、特征点识别的方法将采集的脉搏波信号进行处理；
信号存储模块，用于将处理过的脉搏波信号进行存储；
所述数字滤波，通过一定的算法，对原始数据进行平滑、滤波预处理，去除波形中毛刺，最大限度的保证脉搏波在在幅度及相位上的不失真；
所述周期识别，即一段多周期的脉搏信号中辨识出每个周期的起始点坐标，将脉搏波形序列分割成一个一个的单周期波形数据；
所述基线调整，即有效的去除脉搏波形的基线漂移，消除呼吸运动、肌肉运动或身体位移的因素影响；
所述特征点识别，即标识出单周期内脉搏波形的主波峰值点、降中峡和重搏波。
本发明的进一步技术方案是:所述信号处理模包括:
数字滤波单元，用于利用三点均值滤波算法对原始数据进行平滑、滤波预处理，去除波形中毛刺，最大限度的保证脉搏波在幅度及相位上的不失真；
周期识别单元，用于对一段多周期脉搏波信号进行求导寻找最大值，设置阈值找各个周期峰值点，寻找各个压力上升最快之间的最小值，此最小值点就是一个脉搏波周期的起始点；
基线调整单元，用于利用插值拟合法拟合基线漂移曲线，用原始信号减去漂移曲线得到消除基线漂移的脉搏信号；
特征点识别单元，用于标识出单周期内脉搏波形的主波峰值点、降中峡和重搏波的位置。
本发明的进一步技术方案是:所述周期识别单元包括:
信号一阶求导 模块，用于对一段多周期的脉搏波信号进行一阶求导；
阈值设置模块，用于寻找求导结果中的最大值点，设置阈值；
对比寻峰值点模块，用于求导结果与阈值比较寻找各周期的峰值点；
最小值点寻找模块，用于寻找脉搏周期的起始点。
本发明的进一步技术方案是:所述基线调整单元包括:
幅值与位置确定模块，用于确定单个脉搏波的起始点和结束点的位置和幅值；
拟合曲线模块，用于利用插值拟合法拟合基线漂移曲线；
消除基线漂移模块，用于用原始信号减去漂移曲线得到消除基线漂移的脉搏信号。
本发明的进一步技术方案是:所述特征点识别单元包括:
主波峰值点计算模块，用于求脉搏波信号在起始点与结束点之间的主波峰值点；
计算模块，用于计算单周期内脉搏波形数据的二阶导数及单周期内1/2心动周期位置；
峰值检测模块，对二阶导数结果进行峰值检测，得到波峰波谷数组；
降中峡模块，用于取波峰数组中靠近1/2心动周期的波峰点为降中峡位置；
重搏波模块，用于取波谷数组中靠近降中峡之后的第一个波谷点位置为重搏波位置。
本发明的有益效果是:在单周期识别时，对原始信号波形数据进行截取，截取了完整周期部分，剔除了原始记录数据中起始部分和结束部分不是完整周期的数据。减小了识别错误，也降低了分析算法的复杂程度。在利用原有数据一阶导数寻找降中峡和重博波位置方法上增加了数据二阶导数处理，提高了降中峡和重博波的识别准确性。克服了特征点不明显时产生的误差。通过基线调整，反映出脉搏曲线比较平滑，降低误差率。
图1是本发明实施例提供的脉搏波信号的处理方法流程图2是本发明实施例提供的脉搏波信号处理方法中周期识别流程图3是本发明实施例提供的脉搏波信号处理方法中基线调整流程图4是本发明实施例提供的脉搏波信号处理方法中特征点识别流程图5是本发明实施例提供的脉搏波信号处理系统结构图6是本发明实施例提供的脉搏波信号处理系统中周期识别单元结构图7是本发明实施例提供的脉搏波信号处理系统中基线调整单元结构图8是本发明实施例提供的脉搏波信号处理系统中特征点识别单元结构图。
具体实施方式
附图标记:10_采集显示模块20-信号处理模块30-信号存储模块201-数字滤波单元202-周期识别单元203-基线调整单元204-特征点识别单元2021-信号一阶求导模块2022-阈值设置模块2023对比寻峰值点模块2024-最小值点寻找模块2031-幅值与位置确定模块2032-拟合曲线模块2033-消除基线漂移模块2041-主波峰值点计算模块2042-计算模块2043-峰值检测模块2044-识别降中峡模块2045-识别重搏波模块
图1示出了本发明提供的脉搏波信号的处理方法流程图，详述如下:
步骤SlO中，通过采集装置自人体采集来脉搏波信号，并且将脉搏波信号通过波形的方式显示出来。
步骤S20中，利用三点均值滤波算法对脉搏波信号的原始进行平滑、滤波预处理，去除脉搏波信号波形中毛刺，最大限度的保证脉搏波在幅度及相位上的不失真。三点均值滤波法的具体算法是每一点的值等于和它相邻的前一点和后一点以及该点平均值，即公示
权利要求
1.一种脉搏波信号的处理方法，其特征在于，所述方法包括以下步骤: A、采集脉搏波信号； B、依次采用数字滤波、周期识别、基线调整、特征点识别的方法将采集的脉搏波信号进行处理； C、将处理过的脉搏波信号进行存储； 所述数字滤波，通过一定的算法，对原始数据进行平滑、滤波预处理，去除波形中毛刺，最大限度的保证脉搏波在在幅度及相位上的不失真； 所述周期识别，即一段多周期的脉搏信号中辨识出每个周期的起始点坐标，将脉搏波形序列分割成一个一个的单周期波形数据； 所述基线调整，即有效的去除脉搏波形的基线漂移，消除呼吸运动、肌肉运动或身体位移的因素影响； 所述特征点识别，即标识出单周期内脉搏波形的主波峰值点、降中峡和重搏波。
2.根据权利要求1所述的处理方法，其特征在于，所述B步骤中包括以下步骤: B1、利用三点均值滤波算法对原始数据进行平滑、滤波预处理，去除波形中毛刺，最大限度的保证脉搏波 在幅度及相位上的不失真； B2、对一段多周期脉搏波信号进行求导寻找最大值，设置阈值找各个周期峰值点，寻找各个压力上升最快之间的最小值，此最小值点就是一个脉搏波周期的起始点； B3、利用插值拟合法拟合基线漂移曲线，用原始信号减去漂移曲线得到消除基线漂移的脉搏信号； B4、标识出单周期内脉搏波形的主波峰值点、降中峡和重搏波的位置。
3.根据权利要求2所述的处理方法，其特征在于，所述B2步骤中还包括以下步骤: B21、对一段多周期的脉搏波信号进行一阶求导； B22、寻找求导结果中的最大值点，设置阈值； B23、求导结果与阈值比较寻找各周期的峰值点； B24、寻找脉搏周期的起始点。
4.根据权利要求3所述的处理方法，其特征在于，所述步骤B3还包括以下步骤: B31、确定单个脉搏波的起始点和结束点的位置和幅值； B32、利用插值拟合法拟合基线漂移曲线； B33、用原始信号减去漂移曲线得到消除基线漂移的脉搏信号。
5.根据权利要求4所述的处理方法，其特征在于，所述步骤B4还包括以下步骤: B41、求脉搏波信号在起始点与结束点之间的主波峰值点； B42、计算单周期内脉搏波形数据的二阶导数及单周期内二分之一心动周期位置； B43、对二阶导数结果进行峰值检测，得到波峰波谷数组； B44、取波峰数组中靠近二分之一心动周期的波峰点为降中峡位置； B45、取波谷数组中靠近降中峡之后的第一个波谷点位置为重搏波位置。
6.一种脉搏波信号的处理系统，其特征在于:该处理系统包括: 采集显示模块，用于采集并显示脉搏波信号； 信号处理模块，用于采用数字滤波、周期识别、基线调整、特征点识别的方法将采集的脉搏波信号进行处理；信号存储模块，用于将处理过的脉搏波信号进行存储； 所述数字滤波，通过一定的算法，对原始数据进行平滑、滤波预处理，去除波形中毛刺，最大限度的保证脉搏波在在幅度及相位上的不失真； 所述周期识别，即一段多周期的脉搏信号中辨识出每个周期的起始点坐标，将脉搏波形序列分割成一个一个的单周期波形数据； 所述基线调整，即有效的去除脉搏波形的基线漂移，消除呼吸运动、肌肉运动或身体位移的因素影响； 所述特征点识别，即标识出单周期内脉搏波形的主波峰值点、降中峡和重搏波。
7.根据权利要求6所述的处理系统，其特征在于，所述信号处理模包括: 数字滤波单元，用于利用三点均值滤波算法对原始数据进行平滑、滤波预处理，去除波形中毛刺，最大限度的保证脉搏波在幅度及相位上的不失真； 周期识别单元，用于对一段多周期脉搏波信号进行求导寻找最大值，设置阈值找各个周期峰值点，寻找各个压力上升最快之间的最小值，此最小值点就是一个脉搏波周期的起始点；基线调整单元，用于利用插值拟合法拟合基线漂移曲线，用原始信号减去漂移曲线得到消除基线漂移的脉搏信号；特征点识别单元，用于标识出单周期内脉搏波形的主波峰值点、降中峡和重搏波的位置。
8.根据权利要求7所述的处理系统，其特征在于，所述周期识别单元包括:信号一阶求导模块，用于对一段多周期的脉搏波信号进行一阶求导； 阈值设置模块，用于寻找求导结果中的最大值点，设置阈值；对比寻峰值点模块，用于求导结果与阈值比较寻找各周期的峰值点； 最小值点寻找模块，用于寻找脉搏周期的起始点。
9.根据权利要求8所述 的处理系统，其特征在于，所述基线调整单元包括:幅值与位置确定模块，用于确定单个脉搏波的起始点和结束点的位置和幅值； 拟合曲线模块，用于利用插值拟合法拟合基线漂移曲线； 消除基线漂移模块，用于用原始信号减去漂移曲线得到消除基线漂移的脉搏信号。
10.根据权利要求9所述的处理系统，其特征在于，所述特征点识别单元包括: 主波峰值点计算模块，用于求脉搏波信号在起始点与结束点之间的主波峰值点； 计算模块，用于计算单周期内脉搏波形数据的二阶导数及单周期内1/2心动周期位置； 峰值检测模块，对二阶导数结果进行峰值检测，得到波峰波谷数组； 降中峡模块，用于取波峰数组中靠近1/2心动周期的波峰点为降中峡位置； 重搏波模块，用于取波谷数组中靠近降中峡之后的第一个波谷点位置为重搏波位置。
全文摘要
本发明适用于信号处理领域，提供了一种脉搏波信号的处理方法，所述方法包括以下步骤A、采集脉搏波信号；B、依次采用数字滤波、周期识别、基线调整、特征点识别的方法将采集的脉搏波信号进行处理；C、将处理过的脉搏波信号进行存储。在单周期识别时，对原始信号波形数据进行截取，截取了完整周期部分，剔除了原始记录数据中起始部分和结束部分不是完整周期的数据。减小了识别错误，也降低了分析算法的复杂程度。在利用原有数据一阶导数寻找降中峡和重博波位置方法上增加了数据二阶导数处理，提高了降中峡和重博波的识别准确性。克服了特征点不明显时产生的误差。
文档编号A61B5/02GK103169456SQ201310110100
公开日2013年6月26日 申请日期2013年3月29日 优先权日2013年3月29日
发明者韩君 申请人:深圳职业技术学院