CN1935084A - 一种基于提升小波的脉搏信号处理方法 - Google Patents
一种基于提升小波的脉搏信号处理方法 Download PDFInfo
- Publication number
- CN1935084A CN1935084A CNA2006101132925A CN200610113292A CN1935084A CN 1935084 A CN1935084 A CN 1935084A CN A2006101132925 A CNA2006101132925 A CN A2006101132925A CN 200610113292 A CN200610113292 A CN 200610113292A CN 1935084 A CN1935084 A CN 1935084A
- Authority
- CN
- China
- Prior art keywords
- sample
- waveform
- pulse
- threshold value
- pulse wave
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 12
- 238000000034 method Methods 0.000 claims abstract description 51
- 230000008569 process Effects 0.000 claims abstract description 15
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 13
- 238000001914 filtration Methods 0.000 claims abstract description 11
- 230000003044 adaptive effect Effects 0.000 claims abstract description 6
- 230000001427 coherent effect Effects 0.000 claims abstract description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 22
- 238000012545 processing Methods 0.000 claims description 15
- 238000000605 extraction Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 10
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 7
- 238000011002 quantification Methods 0.000 claims description 6
- 230000003321 amplification Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 claims description 4
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 3
- 230000005055 memory storage Effects 0.000 claims description 3
- 230000000737 periodic effect Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 238000013500 data storage Methods 0.000 claims 1
- 230000010365 information processing Effects 0.000 abstract description 4
- 238000007621 cluster analysis Methods 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 abstract description 2
- 230000003750 conditioning effect Effects 0.000 abstract 1
- 238000013139 quantization Methods 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 6
- 238000003491 array Methods 0.000 description 5
- 239000000284 extract Substances 0.000 description 4
- 230000001575 pathological effect Effects 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明涉及信息处理以及医学信号处理,特别是一种基于提升小波的脉搏信号处理方法。本方法实现步骤:脉搏传感器摄取人体的脉搏波形,脉搏波形经调理电路、AD转换器变成数字信号后输入微处理器,在微处理器中首先进行波形预处理,即采用自适应相干模板法进行滤波。滤波后的脉搏信号采用阈值搜索法获取连续稳定波形,再用时域聚类法剔除稳定波形中在周期上不合格者,得到其中至少包含5个以上波形样本的最优波形。对最优波形做多层小波分解,再对小波分解的高频、低频系数的阈值量化处理,最终得到脉搏波形的特征向量W。本发明通过上述处理方法能够获得脉搏波形的过零点、极值点、突变点等信息组成的特征向量,为进一步的分析诊断打下良好基础。
Description
技术领域
本发明涉及信息处理以及医学信号处理,特别是一种基于提升小波的脉搏信号处理方法。
背景技术
非侵入性测量脉搏的方法,主要是采用压力传感器摄取人体显著动脉部位的脉搏,经过信号预处理,模拟数字转换后,直接利用脉搏波数据进行分析处理及提取特征,计算相关参数。
目前,在对脉搏信息进行分析处理时大多是只对单个周期的脉搏波进行分析,或者采用多个脉搏波进行简单平均处理,并且在选取波形时随意选取经过简单预处理后的波形。这种方法存在一个问题:由于硬件采样频率和各种干扰的问题,使得传感器采得的脉搏波形数据具有不稳定性,如可能存在某些段时间内对脉搏重要信息的漏采和非脉搏信息的误采。这就导致对脉搏信息的后续分析处理失去较好的准确度和正确性,在一定情况下甚至导致误诊。另一方面,目前,在对脉搏波进行分析特征提取时,所用的方法大多比较繁杂,对硬件要求较高,也不适合进行实时处理,使用简单算法的则大多停留在简单的时域特征提取,医用价值不大,无法满足广大家庭用市场的实用型需求。
发明内容
本发明的目的在于克服了现有的处理方法是基于时域或频域分析,并只对单个准周期完整脉搏波进行处理,且算法复杂、对硬件要求高,不利实时处理的缺陷,提出了一种基于提升小波的脉搏信号处理方法。本方法通过脉搏波预处理、波形聚类优选得到优选波形,并从优选波形中提取出特征向量,从而获得信息丰富准确的良好中间参数,为进一步的分析诊断打下良好基础。
为达此目的,本发明提出的基于提升小波的脉搏信号处理方法的实现步骤包括:1)脉搏波形预处理;2)多个完整的典型脉搏波形的优选获取;3)脉搏波形特征的提取。
脉搏波形预处理主要是滤波,利用自适应相干模板法进行滤波处理,该方法能有效虑除工频干扰以及由传感器温差及人体呼吸等引起的基线漂移。并且,该方法运算简单,可以在单片机上实现。
多个完整的典型脉搏波形的优选提取主要利用类似于距离阈值法的聚类分析算法,从而有效避免只利用单个脉搏波形带来的信息误取或漏取,找到一段连续稳定的含有多个完整脉搏波的波形数据后,对该段数据信号进行提升小波变换,获取突变点、极值点、过零点等的提升小波系数构成特征向量。
在进行脉搏波形特征的提取时采用一种改进的整数小波变换算法,该算法不需乘法运算,因此运算量小,对硬件要求低,能够实现实时处理,适合单片机或DSP等硬件实现。通过这种小波变换对小波系数进行软阈值处理,实现特征压缩降低脉搏波特征向量维数。因此,所需存储空间较小。从而能够将大量脉搏病例信息保存在FLASH或外部存储器中。
本发明所采用了如下的技术方案,主要是按以下步骤实现的:
1)由脉搏传感器拾取人体的脉搏波形,脉搏波形经调理电路、AD转换器变成数字信号后输入微处理器,在微处理器中首先经过如下的波形预处理,即采用自适应相干模板法滤波:
采样频率取为工频的整数倍即k×50HZ,设X(n)为原始信号,S(n)为其中的有用信号,N(n)为工频干扰信号,则:
X(n)=S(n)+N(n)
定义模板信号为:
k为采样频率与工频之比,取为自然数;
由于N(n)为周期信号,若S(n)为零均值信号,要进行处理的脉搏波信号可以看成是零均值信号,当M足够大即M≥100时,有:
所以,只要从原始输入信号中减去模板信号就能达到滤除工频干扰的目的,即:
其中参数M的选择:
对于M的取值一般为大于等于100的整数就行,若想达最优的效果与运算速度,则应综合以下几方面来考虑:
①由于(1)式对应的滤波器的截至频率为
为达虑除基线漂移目的,需选择合适的M值,例如为消除频率为0.5Hz的基线漂移则取M=100即可。
②M以2的整数次幂为佳,这样在中低档的嵌入式微处理器中应用时,可用移位操作代替乘除法运算,从而有效地提高运算速度,并节省内存空间。特别是当M取值为256时,甩掉低位字节就可得到误差模板,而无需做任何运算。
③综合考虑运算速度及数据存储量的问题。M值也不宜太大,M值越大,滤波器的响应越慢;而且随着M值的增大,存储原始数据所需内存也相应加大,例如选择M为256时,采样率以200Hz计,仅存储原始数据就需1k字节的空间。
2)多个完整的典型脉搏波形的优选:对于脉搏信号特征向量的提取及分析诊断来说是非常重要的。因为,由压力传感器获取的脉搏信号是非常微弱的生理信号,且极易受各种干扰的影响。本发明算法主要通过以下几个步骤来确保获取的脉搏波形的信息的准确性、丰富性:
①采用阈值搜索法获取连续稳定波形
i)先把前L秒的脉搏波形平均分成N段,其中,15≤L≤60,2≤N≤10,再用差分法分别搜索出每段中的极小值点,即搜索这样的点,差分值由负值到正值的转折点。并将每段中的极小值点的幅值与对应的采样点数,即时间点保存在一个二维数组中,由于把脉搏波形分成了N段,所以分别保存在N个二维数组中。由脉搏波理论及实测波形可知,若为正常波形——指测量过程中的正常,未受严重干扰变形,非指健康。则在每个脉搏周期内的波形有3个极小值点,2个为最低点,1个为重博前波极小值点。
ii)在N个数组中比较各自数组中首个幅值与第二个幅值的大小,找出值较小的点,以此点为起点,做隔点差分计算,并把差分值的绝对值与一设定阈值比较,该阈值具体与硬件电路的放大倍数有关,可根据实验经验获得,一般至小于0.01,找到差分绝对值小于阈值的连续5个以上的点,其对应原波形中的4个以上连续稳定波形。将从这N个数组中搜索出的符合条件的连续点的幅值及其对应的采样时间点保存在另外N个数组中;如果有数组中不包含差分绝对值小于阈值的连续5个以上的点,则继续按步骤2)-①-i)中方法搜索接下来L/N秒中的点,直到搜索出N段均含有4个以上稳定波形为止。
②用时域聚类法剔除在周期上不合格者
i)由步骤2)-①-ii)中得到的信息,回到原波形,就得到多于4×N个准周期的样本波形,对这些准周期样本进行分类;首先设第一个脉搏波的准周期,称之为样本1,为第一聚类中心U1;如果第二个准周期脉搏波即样本U2与U1的距离小于阈值,则将样本2归入第一聚类,U1修正为两个样本的平均值,否则设样本2为第二聚类中心U2。样本U2与U1的距离采用q=1时的闵可夫斯基距离,即周期长度之差,该阈值可根据实验经验具体设定,一般至少小于50ms。如果已经产生U2,再计算样本U3与U1、U2的距离,如果没有产生U2,就计算样本U3与U1的距离,如其中之一小于阈值则可将其归类,并依照上述方法修正聚类中心;否则设样本U3为新的聚类中心,如此可将所有样本分类一次;
ii)如果分类后,没有任何一类包含5个以上的样本,则将阈值放大一倍,重复步骤2)-②-i重新进行分类,如果阈值放大到大于50ms,仍然没有任何一类包含5个以上的样本,则跳到脉搏波形的第二个L秒的信号,按照步骤1)、2)中所述步骤重新开始寻找4×N个以上的稳定周期脉搏波并重复上述聚类过程;
iii)得到其中至少包含5个以上波形样本的聚类后,则此类样本即为所需最优波形;将其中包含样本数最多的一类的聚类中心Umax定为标准脉搏周期,其倒数即为脉率。则此类样本即为所需最优波形。
3)脉搏波形特征的提取:本步骤的目的是为了获取脉搏波形中隐含着丰富的生理病理信息的特征参数。本发明算法以提升小波来提取脉搏波形中的奇异点等信息作为特征向量。其要点如下:
①提升算法
i)使用Donoho小波取其合成尺度函数和分解尺度函数分别有N和M阶消失矩,进行提升;
ii)其提升格式为:
分解:
预测:
更新:
其中sl,dl为分解小波的逼近系数与细节系数;
iii)为将上述提升格式进行分解为只需加法与移位运算,定义下面的运算格式:设P,U为两个位移加法处理单元:
P(a,b,c,d,)=(-a+b+c-d+b<<3+c<<3+8)>>4
U(a,b,c,d,)=(-a+b+c-d+b<<3+c<<3+16)>>5
其中,<<是左位移符号,>>是右位移符号,则
sl=x2l,dl=x2l+1;
dl=dl-p(sl-1,sl,sl+1,sl+2);
sl=sl+(dl-2,dl-1,dl,dl+1).
②对步骤2)中得到的5个以上周期的优选波形按照3)-①中所述的提升算法进行提升小波分解,并作如下处理:
i、对所得脉搏波形做多层小波分解;
ii、小波分解的高频系数的阈值量化,对深层的高频系数,选择阈值进行软阈值量化处理,表达式如下:
式中w是高频小波系数,sgn是符号函数,σ是一个由经验所得阈值,一般在0.5~1之间。由此得到多个一维向量A,B,C,D,E,对这些一维向量对应元素a0,b0,c0,d0,e0,……排序,取中值作为最终的小波系数w′0,同理,依次对a1,b1,c1,d1,e1,……直至an,bn,cn,dn,en,……作相同处理,可得w′1……w′n,如此可得一个小波高频系数一维向量w′。
iii、小波分解的低频系数的阈值量化,对第1层的低频系数,选择阈值进行软阈值量化处理,阈值选择及处理方法同3)-②-ii,得到另一个一维向量w″。
iv、记W=[w′,w″],最终获得脉搏波形的特征向量W。
本发明中所得特征信息是脉搏信号处理的中间结果,可用做进一步分析诊断之用。比如可以直接根据提取的特征向量运用聚类分析等模式识别方法进行疾病的诊断分析,也可以利用过零点、极值点、突变点等的小波系数进行信号重建然后在微机上运用谱分析等各种时频域方法进行处理,以供进一步研究之用。因此,本发明对于硬件没有特别的要求,计算简单,易于在实际简易***中实现,并且能够实现实时处理。
本发明提出的是一种脉搏处理及特征提取方案,而现有的处理方法多数是基于时域或频域分析的,而且都是单个准周期完整脉搏波进行处理。本发明采用一种基于提升小波变换的奇异点信息分析的特征提取方法。与当前其他方法相比,本发明的新颖性主要体现在:
1)在波形预处理上采取一种自适应相干模板法,能够有效消除工频干扰及由温漂、呼吸等引起的基线漂移等干扰。且该算法运算简单,能够在单片机等硬件上实时实现。
2)在脉搏波形优选上先采取简单的类似距离阈值法的聚类算法,同时获取多个优质波形,再采用q=1时的闵可夫斯基距离,即周期长度之差阈值法进行聚类处理,进一步剔除非优波形。同时可获得波形准周期,即得脉率。以上算法简单有效,能较好地避免信息的漏采、误采,同时获得了较为准确的脉率。
3)在特征提取上,先对获取的多个优选脉波形按准周期分段进行提升小波变换,对得到的小波系数进行阈值处理,从而得到一维特征向量,简洁有效。本方法中的小波提升算法是一种改进的整数到整数的提升小波变换方法,不需进行乘法运算,因而对硬件要求低,运算速度快,甚至能够由单片机实现实时处理。
附图说明
图1***处理流程框图
图2具体算法处理流程图
图3a是采集的原始脉搏波形
图3b是经滤波处理后的脉搏波形
图4脉搏信号波形图
图5a为第六层高频系数
图5b为第一层低频小波系数
图6***程序模块数据处理流程图
具体实施方式
下面结合图1~图6说明本实施例。
本实施例包括三个主要步骤:信号预处理、多个完整的典型脉搏波形的优选和脉搏波形特征的提取。
本发明的技术实现方案如图1~图2所示。由脉搏传感器获的波形经过由一级放大、滤波、二级放大等电路组成的调理电路通过AD转换器变成数字流进入微处理器***进行一系列的处理最终提取出小波特征系数向量,将这些向量数据保存在单片机或dsp等仪器的存储器或外部存储器FLASH中,以便后续使用。
本实施例提出的基于提升小波及聚类分析的脉搏信号特征提取方法的实现包括以下三个步骤:
1.脉搏波形预处理
1)自适应相干模板法滤波
取采样频率为工频的整数倍即50×k=200HZ,设X(n)为原始信号,S(n)为其中的有用信号,N(n)为工频干扰信号,则:
X(n)=S(n)+N(n)
定义模板信号为:
k为采样频率与工频之比,取为自然数。
由于N(n)为周期信号,若S(n)为零均值信号,当M足够大(M>=100)时,有:
所以,只要从原始输入信号中减去模板信号就能达到滤除工频干扰的目的,即:
2)参数选择:本实施例中选择M=256,k=4。
2、多个完整的典型脉搏波形的优选
1)阈值搜索法获取连续稳定波形。包括以下几步:
①先把前30秒分成3段,则每段一般会包含10个以上波形,用差分法搜索出每段中的极小值点,即搜索这样的点,差分值由负值到正值的转折点。并将三段中的这些点的幅值与对应的采样点数,即时间点,分别保存在3个二维数组中,即每段中的这些点的幅值与对应的采样点数都保存在1个二维数组中,3段共保存在三个二维数组中。由脉搏波理论及实测波形可知,若为正常波形(指测量过程中的正常,未受严重干扰变形,非指健康),则在每个脉搏周期内的波形有3个极小值点,2个为最低点,1个为重博前波极小值点,见图3a。
②在3个数组中比较各自数组中首个幅值与第二个幅值大小,找出值较小的点,以此点为起点,做隔点差分计算,并把差分值的绝对值与一设定阈值比较,该阈值取为0.005。找到差分绝对值小于阈值的连续5个以上的点,其对应原波形中的4个以上连续稳定波形。将从这3个数组中搜索出的符合条件的连续点的幅值及其对应的采样时间点保存在另外3个数组中。如果有数组中不包含这样的点,则继续按步骤①中方法搜索接下来10秒中的点,直到搜索出3段均含有4个以上稳定波形为止;
2)用时域聚类法剔除在周期上不合格者。
①由步骤1)-②中得到的信息,回到原波形,就得到多于12个准周期的样本波形。对这些准周期样本进行分类。首先设第一个脉搏波的准周期,称之为样本1,为第一聚类中心U1。如果第二个准周期脉搏波即样本2与U1的距离(此处采用q=1时的闵可夫斯基距离,即周期长度之差)小于阈值,该阈值取为40ms,则将样本2归入第一聚类,U1修正为两个样本的平均值;否则设样本2为第二聚类中心U2。而后计算样本3与U1,U2(如果已经产生)的距离,如其中之一小于阈值则可将其归类,修正聚类中心;否则设样本3为新的聚类中心。如此可将所有样本分类一次;
②如果分类后,没有任何一类包含5个以上的样本,则将阈值放大一倍,重复步骤①,重新进行分类。如果阈值放大到大于50ms,仍然没有任何一类包含5个以上的样本,则跳到下一段信号按照步骤1)中所述重新开始寻找12个以上的稳定周期脉搏波并重复上述聚类过程;
③得到其中至少包含5个以上波形样本的聚类后,将其中包含样本数最多的一类的聚类中心Umax定为标准脉搏周期,其倒数即为脉率。则此类样本即为所需最优波形。在本实例中含最多样本的聚类中样本个数为5个。
3.脉搏波形特征的提取的目的是为了获取脉搏波形中隐含着丰富的生理病理信息的特征参数。本发明算法以提升小波来提取脉搏波形中的奇异点等信息作为特征向量。其要点如下:
1)提升算法
①使用Donoho小波取其合成尺度函数和分解尺度函数分别有N和M阶消失矩,取N=M=4,进行提升;
②其提升格式为:
分解:
预测:
更新:
其中sl,dl为分解小波的逼近系数与细节系数。
③为将上述提升格式进行分解为只需加法与移位运算,定义下面的运算格式:设P,U为两个位移加法处理单元:
P(a,b,c,d)=(-a+b+c-d+b<<3+c<<3+8)>>4
U(a,b,c,d,)=(-a+b+c-d+b<<3+c<<3+16)>>5
其中,<<是左位移符号,>>是右位移符号。则
sl=x2l,dl,=x2l+1;
dl=dl-p(sl-1,sl,sl+1,sl+2);
sl=sl+(dl-2,dl-1,dl,dl+1).
经过上述处理,小波提升运算过程就完全能由加法及移位实现。
2)对2中所得5个周期的优选波形按照1)中所述算法进行提升小波分解,并作如下处理:
①对所得脉搏波形做6层小波分解;
②小波分解的高频系数的阈值量化。对深层的高频系数,选择阈值进行软阈值量化处理。表达式如下:
式中w是高频小波系数,sgn是符号函数,σ是一个由经验所得阈值,一般在0.5到1之间。由此得到多个一维向量A,B,C,D,E。对这些一维向量对应元素a0,b0,c0,d0,e0排序,取中值做为最终的小波系数w′0,同理,依次对a1,b1,c1,d1,e1,直至an,bn,cn,dn,en作相同处理,可得w′1……w′n。如此可得一个小波高频系数一维向量w′。
③小波分解的低频系数的阈值量化。对第1层的低频系数,选择阈值进行软阈值量化处理。处理方法同②,得到另一个一维向量w″。
④记W=[w′,w″],最终获得脉搏波形的特征向量W。
Claims (1)
1、一种基于提升小波的脉搏信号处理方法,其特征在于,该方法是按以下步骤实现的:
1)由脉搏传感器摄取人体的脉搏波形,脉搏波形经调理电路、AD转换器变成数字信号后输入微处理器中的数据存储器RAM,微处理器调用其程序存储器ROM中的各处理模块进行如下处理:
首先调用波形预处理模块进行运算处理,该模块采用自适应相干模板法进行滤波,其算法如下:
取采样频率为工频的整数倍即k×50HZ,设X(n)为原始信号,S(n)为其中的有用信号,N(n)为工频干扰信号,则:
X(n)=S(n)+N(n)
定义模板信号为:
k为采样频率与工频之比,取为自然数;
由于N(n)为周期信号,若S(n)为零均值信号,要进行处理的脉搏波信号可看成是零均值信号,当M足够大即M≥100时,有:
从原始输入信号中减去模板信号就能达到滤除工频干扰的目的,即:
经过上述程序模块处理后将数据保存RAM中,微处理器再调用波形优选模块进行处理,该模块算法实现如下:
2)多个完整的典型脉搏波形的优选,
①采用阈值搜索法获取连续稳定波形
i)先把前L秒的脉搏波形平均分成N段,其中,15≤L≤60,2≤N≤10,再用差分法分别搜索出每段中的极小值点,并将N段中的极小值点的幅值与对应的采样点数,即时间点分别保存在N个二维数组中;
ii)在N个数组中比较各自数组中首个幅值与第二个幅值的大小,找出值较小的点,以此点为起点,做隔点差分计算,并把差分值的绝对值与一设定阈值比较,该阈值具体与硬件电路的放大倍数有关,可根据实验经验获得,一般至少小于0.01,找到差分绝对值小于阈值的连续5个以上的点,其对应原波形中的4个以上连续稳定波形;将从这N个数组中搜索出的符合条件的连续点的幅值及其对应的采样时间点分别保存在另外N个数组中;如果有数组中不包含差分绝对值小于阈值的连续5个以上的点,则继续按步骤2)-①-i)中方法搜索接下来L/N秒中的点,直到搜索出N段均含有4个以上稳定波形为止;
②用时域聚类法剔除在周期上不合格者
i)由步骤2)-①-ii)中得到的信息,回到原波形,就得到多于4×N个准周期的样本波形,对这些准周期样本进行分类;首先设第一个脉搏波的准周期,称之为样本1,为第一聚类中心U1;如果第二个准周期脉搏波即样本U2与U1的距离小于阈值,则将样本2归入第一聚类,U1修正为两个样本的平均值,否则设样本2为第二聚类中心U2;样本U2与U1的距离采用q=1时的闵可夫斯基距离,即周期长度之差,阈值可根据实验经验具体设定,一般至少小于50ms;如果已经产生U2,再计算样本U3与U1、U2的距离,如果没有产生U2,就计算样本U3与U1的距离,如其中之一小于阈值则可将其归类,并依照上述方法修正聚类中心;否则设样本U3为新的聚类中心,如此可将所有样本分类一次;
ii)如果分类后,没有任何一类包含5个以上的样本,则将阈值放大一倍,重复步骤2)-②-i重新进行分类,如果阈值放大到大于50ms,仍然没有任何一类包含5个以上的样本,则跳到脉搏波形的第二个L秒的信号,按照步骤1)、2)中所述步骤重新开始寻找4×N个以上的稳定周期脉搏波并重复上述聚类过程;
iii)得到其中至少包含5个以上波形样本的聚类后,则此类样本即为所需最优波形;将其中包含样本数最多的一类的聚类中心Umax定为标准脉搏周期,其倒数即为脉率。
3)脉搏波形特征的提取
得到优选波形后,可以只保存下该波形数据在RAM或者外部存储器中,同时微处理器调用最后一个模块即特征提取模块,该程序模块算法采用如下提升小波算法:
①提升算法
i)使用Donoho小波取其合成尺度函数和分解尺度函数分别有N和M阶消失矩,进行提升;
ii)其提升格式为:
分解:
预测:
更新:
其中Sl,dl为分解小波的逼近系数与细节系数;
iii)为将上述提升格式进行分解为只需加法与移位运算,定义下面的运算格式:设P,U为两个位移加法处理单元:
P(a,b,c,d)=(-a+b+c-d+b<<3+c<<3+8)>>4
U(a,b,c,d)=(-a+b+c-d+b<<3+c<<3+16)>>5
其中,<<是左位移符号,>>是右位移符号,则
Sl=x2l,dl=x2l+1;
dl=dl-p(Sl-1,Sl,Sl+1,Sl+2);
Sl=Sl+(dl-2,dl-1,dl,dl+1)。
②对步骤2)中得到的5个以上周期的优选波形按照3)-①中所述的提升算法进行提升小波分解,并作如下处理:
i、对所得脉搏波形做多层小波分解;
ii、小波分解的高频系数的阈值量化,对深层的高频系数,选择阈值进行软阈值量化处理,表达式如下:
式中w是高频小波系数,sgn是符号函数,σ是一个由经验所得阈值,一般在0.5到1之间;由此得到多个一维向量A,B,C,D,E,……对这些一维向量对应元素a0,b0,c0,d0,e0,……排序,取中值做为最终的小波系数w′0,同理,依次对a1,b1,c1,d1,e1,……直至an,bn,cn,dn,en,……作相同处理,可得w′1……w′n,如此可得一个小波高频系数一维向量w′;
iii、小波分解的低频系数的阈值量化,对第1层的低频系数,选择阈值进行软阈值量化处理,处理方法同3)-②-ii,得到另一个一维向量w″;
iv、记W=[w′,w″],最终获得脉搏波形的特征向量W。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006101132925A CN100457023C (zh) | 2006-09-22 | 2006-09-22 | 一种基于提升小波的脉搏信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006101132925A CN100457023C (zh) | 2006-09-22 | 2006-09-22 | 一种基于提升小波的脉搏信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1935084A true CN1935084A (zh) | 2007-03-28 |
CN100457023C CN100457023C (zh) | 2009-02-04 |
Family
ID=37952961
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006101132925A Expired - Fee Related CN100457023C (zh) | 2006-09-22 | 2006-09-22 | 一种基于提升小波的脉搏信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100457023C (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101332084B (zh) * | 2007-06-28 | 2011-01-26 | 罗晓民 | 脉搏波的面积重心位移分析方法 |
CN102247128A (zh) * | 2010-05-19 | 2011-11-23 | 中国科学院计算技术研究所 | 人体脉搏信息采集装置及人体健康状况监护装置 |
CN102458237A (zh) * | 2009-06-17 | 2012-05-16 | 心动力医疗公司 | 获得和处理心冲击描计图数据的方法和设备 |
CN101683260B (zh) * | 2008-09-28 | 2012-05-30 | 深圳迈瑞生物医疗电子股份有限公司 | 一种节律性生理信号的优化处理方法及装置 |
CN101866423B (zh) * | 2010-05-12 | 2012-07-25 | 重庆大学 | 一种动脉压信号逐拍分割方法 |
CN101642378B (zh) * | 2008-08-06 | 2012-10-10 | 深圳迈瑞生物医疗电子股份有限公司 | 控制谐波成像的发射脉冲频谱的方法和装置 |
CN102824166A (zh) * | 2012-07-05 | 2012-12-19 | 华东师范大学 | 一种用于脉搏波分类处理的装置 |
CN103153175A (zh) * | 2010-05-14 | 2013-06-12 | 先进计算发展中心 | 脉搏信号波形数据的诊断分类 |
CN103479334A (zh) * | 2013-09-22 | 2014-01-01 | 江苏美伦影像***有限公司 | 一种人体脉搏无线检测*** |
CN104224140A (zh) * | 2014-09-09 | 2014-12-24 | 桂林电子科技大学 | 一种利用提升小波变换滤除基线漂移的方法 |
CN106073741A (zh) * | 2016-05-30 | 2016-11-09 | 珠海脉动时代健康科技有限公司 | 一种自适应滤波及计算脉搏的***和方法 |
CN107970027A (zh) * | 2017-11-23 | 2018-05-01 | 安徽大学 | 一种桡动脉检测与人体体质识别***和方法 |
CN108042117A (zh) * | 2018-01-03 | 2018-05-18 | 刘素杰 | 一种脉搏检测***及其检测方法 |
CN108135513A (zh) * | 2016-08-09 | 2018-06-08 | 深圳市汇顶科技股份有限公司 | 一种基于心率信号的按键操作触发方法、装置及*** |
CN109833034A (zh) * | 2017-11-24 | 2019-06-04 | 深圳市岩尚科技有限公司 | 一种脉搏波信号中提取血压数据的方法及装置 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3521654B2 (ja) * | 1996-11-07 | 2004-04-19 | セイコーエプソン株式会社 | 脈拍計 |
CN1073403C (zh) * | 1997-06-28 | 2001-10-24 | 北京工业大学 | 确定体外反搏过程中的最优加压时间的方法及采用该方法的体外反搏装置 |
JP4401014B2 (ja) * | 2000-08-31 | 2010-01-20 | 株式会社デンソー | 脈波解析装置 |
SE0301201D0 (sv) * | 2003-04-24 | 2003-04-24 | St Jude Medical | Apparatus for analysing cardiac events |
CN100427029C (zh) * | 2005-01-20 | 2008-10-22 | 北京工业大学 | 基于脉搏熵的多点脉搏信息融合方法 |
CN2808060Y (zh) * | 2005-05-27 | 2006-08-23 | 北京工业大学 | 基于pci总线方式的脉搏信号采集处理器 |
-
2006
- 2006-09-22 CN CNB2006101132925A patent/CN100457023C/zh not_active Expired - Fee Related
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101332084B (zh) * | 2007-06-28 | 2011-01-26 | 罗晓民 | 脉搏波的面积重心位移分析方法 |
CN101642378B (zh) * | 2008-08-06 | 2012-10-10 | 深圳迈瑞生物医疗电子股份有限公司 | 控制谐波成像的发射脉冲频谱的方法和装置 |
CN101683260B (zh) * | 2008-09-28 | 2012-05-30 | 深圳迈瑞生物医疗电子股份有限公司 | 一种节律性生理信号的优化处理方法及装置 |
CN102458237A (zh) * | 2009-06-17 | 2012-05-16 | 心动力医疗公司 | 获得和处理心冲击描计图数据的方法和设备 |
CN101866423B (zh) * | 2010-05-12 | 2012-07-25 | 重庆大学 | 一种动脉压信号逐拍分割方法 |
CN103153175B (zh) * | 2010-05-14 | 2017-02-15 | 先进计算发展中心 | 用于做出脉搏信号波形数据的诊断分类的设备和*** |
CN103153175A (zh) * | 2010-05-14 | 2013-06-12 | 先进计算发展中心 | 脉搏信号波形数据的诊断分类 |
US9220437B2 (en) | 2010-05-14 | 2015-12-29 | Centre For Development Of Advanced Computing | Diagnostic classifications of pulse signal waveform data |
CN102247128A (zh) * | 2010-05-19 | 2011-11-23 | 中国科学院计算技术研究所 | 人体脉搏信息采集装置及人体健康状况监护装置 |
CN102247128B (zh) * | 2010-05-19 | 2014-05-07 | 中国科学院计算技术研究所 | 人体脉搏信息采集装置及人体健康状况监护装置 |
CN102824166A (zh) * | 2012-07-05 | 2012-12-19 | 华东师范大学 | 一种用于脉搏波分类处理的装置 |
CN103479334A (zh) * | 2013-09-22 | 2014-01-01 | 江苏美伦影像***有限公司 | 一种人体脉搏无线检测*** |
CN104224140A (zh) * | 2014-09-09 | 2014-12-24 | 桂林电子科技大学 | 一种利用提升小波变换滤除基线漂移的方法 |
CN106073741A (zh) * | 2016-05-30 | 2016-11-09 | 珠海脉动时代健康科技有限公司 | 一种自适应滤波及计算脉搏的***和方法 |
CN106073741B (zh) * | 2016-05-30 | 2019-03-29 | 珠海脉动时代健康科技有限公司 | 一种自适应滤波及计算脉搏的***和方法 |
CN108135513A (zh) * | 2016-08-09 | 2018-06-08 | 深圳市汇顶科技股份有限公司 | 一种基于心率信号的按键操作触发方法、装置及*** |
CN108135513B (zh) * | 2016-08-09 | 2021-02-12 | 深圳市汇顶科技股份有限公司 | 一种基于心率信号的按键操作触发方法、装置及*** |
CN107970027A (zh) * | 2017-11-23 | 2018-05-01 | 安徽大学 | 一种桡动脉检测与人体体质识别***和方法 |
CN107970027B (zh) * | 2017-11-23 | 2020-12-29 | 安徽大学 | 一种桡动脉检测与人体体质识别***和方法 |
CN109833034A (zh) * | 2017-11-24 | 2019-06-04 | 深圳市岩尚科技有限公司 | 一种脉搏波信号中提取血压数据的方法及装置 |
CN108042117A (zh) * | 2018-01-03 | 2018-05-18 | 刘素杰 | 一种脉搏检测***及其检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN100457023C (zh) | 2009-02-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100457023C (zh) | 一种基于提升小波的脉搏信号处理方法 | |
Singh et al. | Classification of unsegmented heart sound recording using KNN classifier | |
CN104706349B (zh) | 一种基于脉搏波信号的心电信号构建方法 | |
CN106691406A (zh) | 一种基于单探头光电容积脉搏波的血管弹性和血压的检测方法 | |
CN112914527B (zh) | 一种基于脉搏波光电容积描记术的动脉血压信号采集方法 | |
CN111759345B (zh) | 基于卷积神经网络的心脏瓣膜异常分析方法、***和装置 | |
CN101919704B (zh) | 一种心音信号定位、分段方法 | |
CN101259016A (zh) | 实时自动检测癫痫特征波的方法 | |
CN108567418A (zh) | 一种基于PCANet的脉搏信号亚健康检测方法及检测*** | |
CN104473660B (zh) | 一种基于子带能量包络自相关特征的异常心音识别方法 | |
Aravind Kumar et al. | Efficient FPGA‐based VLSI architecture for detecting R‐peaks in electrocardiogram signal by combining Shannon energy with Hilbert transform | |
CN104305992A (zh) | 一种交互式胎儿心电快速自动提取方法 | |
CN110141215A (zh) | 降噪自编码器的训练方法、心电信号的降噪方法及相关装置、设备 | |
Ahlstrom et al. | Heart sound cancellation from lung sound recordings using recurrence time statistics and nonlinear prediction | |
CN102240208A (zh) | 一种可在单片机中实现的心电去噪小波算法 | |
Yen et al. | Blood Pressure and Heart Rate Measurements Using Photoplethysmography with Modified LRCN. | |
CN113729653A (zh) | 一种人体脉搏波信号采集方法 | |
CN115281676B (zh) | 基于gru神经网络和ecg信号的疲劳检测方法 | |
CN116942172A (zh) | 一种基于编解码结构的小波双通道单导联心电去噪方法 | |
CN117017297A (zh) | 驾驶员疲劳的预测和识别模型建立方法及其应用 | |
CN112998676A (zh) | 基于光电阵列增强信号多特征提取的连续血压测量方法 | |
CN114145725B (zh) | 一种基于无创连续血压测量的ppg采样率估算方法 | |
Nibhanupudi | Signal denoising using wavelets | |
CN112826459B (zh) | 一种基于卷积自编码器的脉搏波波形重建的方法及*** | |
CN114668375A (zh) | 一种基于深度神经网络模型的无袖带血压预测*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090204 Termination date: 20100922 |