CN1923133A - 利用信号变换计算血压的方法和装置 - Google Patents

Abstract

一种利用信号变换计算血压的方法和装置，所述装置包括压力传感器、压力、脉搏波放大电路、A/D转换、微处理器、控制电路、与袖套气路连接的充气泵、快放阀和慢放阀以及信号处理模块和比率法分析模块；对采样的脉搏波时间序列进行信号变换，由变换后的信号求得最大幅度及对应的频率值，最大幅度作为当前压力下的幅度值，连续对不同压力台阶下的脉搏波时间序列进行信号变换，从而获得脉搏波幅度和压力的关系曲线，从脉搏波幅度和压力的关系曲线得到平均压和对应的脉率，最后利用比率法计算出收缩压和舒张压。本发明的方法和装置能在受测者出现抖动、肌颤等情况下减少测量误差，提高血压测量的准确性。

Description

利用信号变换计算血压的方法和装置

技术领域

本发明涉及测量心和血管压力的方法和装置，尤其涉及利用对采样信号进行信号变换以提高测量精度的方法和装置。

背景技术

传统的利用振荡法测量血压的方法，是通过搜索不同压力下的脉搏波幅度，得到脉搏压力和脉搏波幅度的对应关系曲线。然后查找最大脉搏波幅度对应的压力，认为是平均压；利用平均压对应的脉搏波幅度，乘以收缩压幅度系数，得到收缩压对应的脉搏波幅度，然后根据对应关系曲线计算出脉搏压力，认为是收缩压；利用平均压对应的脉搏波幅度，乘以舒张压幅度系数，得到舒张压对应的脉搏波幅度，然后根据对应关系曲线计算出脉搏压力，认为是舒张压。传统血压测量装置的硬件构成如图1所示。从电路功能上可分为：电源电路、马达，阀门，和袖套的充气通路、压力信号采样放大处理部分、A/D转换电路、逻辑控制部分、微处理器数据处理部分、串口通讯部分，硬件保护电路。各部分的具体功能如下：

电源电路输入12V直流电源，输出5组电源——数字+5V、模拟+5V，模拟+2.5v，模拟+8v，该电路为整个血压板提供电源。

马达、阀门和袖套的充气通路，受逻辑控制部分的控制。首先关闭阀门，给马达上电，充气到指定的压力，然后停止充气，打开阀门，逐渐放气。

压力传感器将当前的压力信号转换为电信号，一路直接输出到A/D，计算后获得当前的压力；另外一路通过一个高通滤波，滤掉直流的压力信号，输出到A/D，获得交流的脉搏波信号。

A/D转换部分将经放大处理后的模拟信号转换为数字信号，然后送微处理器处理。微处理器数据处理部分将A/D转换得到的数据处理、计算，得到血压和脉率并通过串口向外发送。

硬件保护电路包括一个压力传感器，放大电路和比较器，在压力信号超过一定幅度后，引起比较器的输入超过阈值翻转，输出低电平引起系统复位。

进行无创血压测量时，先给袖带充气，至一定压力后缓慢放气。如果充气所达到的压力大于人体的收缩压，则可以认为血管已完全阻断，此时，在袖带处应该不存在脉搏的搏动。然后开始放气，当袖带压力小于收缩压时，血管已经有部分导通，此时脉搏搏动会随着袖带压力的下降而逐渐增强，由于脉搏搏动的影响，将使袖带的压力出现小范围的波动，也即是在压力信号上叠加了一个振荡信号，同样这种振荡信号随着袖带压力的下降而逐渐增强。但当袖带压力减小到一定程度后，振荡信号的幅度反而会开始下降。这主要是因为袖带压力减小，使得人体皮下组织对脉搏搏动的衰减作用逐渐增强。当袖带压力持续下降，这种衰减作用越来越明显，振荡信号的幅度下降也越明显。

测量过程袖套绑缚在受测者胳膊上并充气到一定压力，受测者不可避免的感觉到压迫，在无意识的情况下会出现抖动；或者受测者的病情导致出现肌颤等情况。这些情况下将使脉搏波的基线和幅度受到影响。因此，现有技术寻找脉搏波的波峰和波谷的方法很容易受到上述抖动、肌颤等的干扰，导致脉搏波幅度及相关的测量计算出现较大的误差。

发明内容

本发明要解决的技术问题是提出一种能明显减少由于受测者抖动、肌颤等引起的测量误差的血压测量方法和装置。

本发明采用如下技术方案：设计一种利用信号变换计算血压的方法，其特征在于该方法包括如下步骤：

a.把带有压力传感器的袖套绑缚在被测者的胳膊上并充气到一定压力后打开阀门放气，放气到指定压力；

b.维持该指定压力在DT秒内不变，在DT秒内压力传感器不间断地把当前的压力和脉搏波信号转换为电信号，该电信号经放大后被送至A/D转换变成压力和脉搏波数字信号的时间序列；

c.对所述脉搏波时间序列进行信号变换，利用信号变换计算血压计算最大幅度及对应的频率值，通过最大幅度值求取当前压力下的幅度值；最大幅度对应的频率值作为当前压力下的脉率值，继续放气至下一个指定压力，返回到步骤b；

d.步骤b和步骤c被循环执行NUM次，获得NUM个在不同压力台阶上采样的脉搏波幅度序列和压力序列，用这两个数据序列构成脉搏波幅度和压力的关系曲线；

e.利用所述脉搏波幅度和压力的关系曲线计算平均压Mean及对应的脉搏波幅度P0，然后利用比率法进一步计算出所要测量的收缩压SYS和舒张压DIA对应的脉搏波幅度PS和PD，其中PS＝α*P0，PD＝β*P0，然后根据PS和PD并利用脉搏波幅度和压力的关系曲线求出收缩压和舒张压。

所述步骤e中比率法的两个比率α和β是通过大量的实验获得的与装置和算法有关的经验系数，α的取值范围为0.2至0.9，β的取值范围为0.2至0.8。所述NUM的取值范围为6至10，所述DT的取值范围为3至10秒。

所述对脉搏波时间序列进行信号变换指对脉搏波时间序列进行CHIRPZ变换，通过变换后的信号序列计算最大幅度及对应的频率值，最大幅度作为当前压力下的幅度值。

对所述步骤d获得的NUM个在不同压力台阶上采样的的脉搏波幅度序列和压力序列进行归一化、插值、加权平均处理，使用处理后的脉搏波幅度和压力序列构成脉搏波幅度和压力的关系曲线，以保证测量结果更加准确。

本发明解决技术问题所采用的技术方案还包括：设计一种利用信号变换计算血压的装置，包括：

数据采样单元，包括位于袖套内的压力传感器、压力放大电路和脉搏波放大电路以及A/D转换；所述压力传感器把当前的压力和脉搏波信号转换为电信号，该电信号经压力放大电路和脉搏波放大电路放大后被送至A/D转换变成压力和脉搏波的数字信号；

主控单元，包括含有计算机软件的微处理器和控制电路，用于根据压力信号控制袖套放气的时间和速度以及计算被测人体的血压和脉率；

充放气单元，包括与所述袖套气路连接的充气泵、快放阀和慢放阀，用于接受主控单元的控制信号对所述袖套实施自动的充、放气操作；其特征在于还包括：

信号处理模块，用于对A/D转换输出的多个压力台阶下的脉搏波的数字信号进行信号变换，得到被分析数据的幅频特性从而进一步得到与被测人体血压紧密相关的波幅最大值；

比率法分析模块，用于调取所述信号处理模块输出的多个压力台阶下的波幅最大值以及A/D转换输出的多个压力台阶的压力值，以形成脉搏波幅度和压力的关系曲线，根据该关系曲线找出测量过程的波幅最大值，再根据系统确定的比率计算和输出被测人体的血压和脉率。

所述比率法分析模块还包括：归一化模块、插值模块和加权平均模块，用于在计算和输出被测人体的血压和脉率之前对脉搏波幅度和压力的关系曲线进行归一化、插值和加权平均处理。

与现有技术相比较，本发明利用信号变换计算血压的方法和装置具有如下优点：通过采用信号变换的方法计算脉搏波幅度，提高了抗干扰的能力，在受测者出现抖动、肌颤等情况下，可以有效提高血压测量的准确性。

附图说明

图1是现有血压测量装置硬件构成的原理框图；

图2是本发明利用信号变换计算血压的装置的原理框图；

图3是图2中信号处理模块和比率法分析模块的算法流程图；

图4是对不同脉率下的脉搏波信号进行CHIRPZ变换的示意图，其中(a)代表40脉率的情况，(b)代表80脉率的情况，(c)代表240脉率的情况；

图5是系统实现流程图。

具体实施方式

以下结合附图及附图所示之实施例对本发明装置和方法作进一步详述。

由于本发明装置区别现有技术之处在于对采样的脉搏波进行处理的算法模块，硬件部分与现有技术与现有技术基本相同，因此我们在描述本发明装置时对与现有装置相同的硬件部分尽量简化。如图2所示，本发明利用信号变换计算血压的装置包括：

数据采样单元，包括位于袖套内的压力传感器、压力放大电路和脉搏波放大电路以及A/D转换；所述压力传感器把被测者的血压和脉搏波信号转换为电信号，该电信号经压力放大电路和脉搏波放大电路放大后被送至A/D转换变成压力和脉搏波的数字信号；

主控单元，包括含有计算机软件的微处理器和控制电路，用于根据测量需要设定或根据当前的压力信号控制袖套放气的时间和速度以及根据算法软件计算被测人体的血压和脉率；

充放气单元，包括与所述袖套气路连接的充气泵、快放阀和慢放阀，用于接受主控单元的控制信号对所述袖套实施自动的充、放气操作；

辅助功能模块，用于系统复位、漏气检测、压力校准等；本发明装置有别于现有技术之处在于还包括：

信号处理模块，用于对A/D转换输出的多个压力台阶下的脉搏波的数字信号进行信号处理，得到被分析数据的幅频特性从而进一步得到与被测人体血压紧密相关的波幅最大值，本发明的实施例采用CHIRPZ变换，当然也可以采用其他形式的变换；

比率法分析模块，用于调取所述信号处理模块输出的多个压力台阶下的波幅最大值以及A/D转换输出的多个压力台阶的压力值，以形成脉搏波幅度和压力的关系曲线；

所述比率法分析模块还包括：归一化模块、插值模块和加权平均模块，用于在计算和输出被测人体的血压和脉率之前对脉搏波幅度和压力的关系曲线进行归一化、插值和加权平均处理；然后根据处理后的脉搏波幅度和压力关系曲线找出测量过程的波幅最大值，最后根据系统确定的比率计算和输出被测人体的血压和脉率。

下面较详细地说明以本发明装置为硬件基础的利用信号变换计算血压的方法的一个具体实施例，所述方法包括如下步骤：

a.把带有压力传感器的袖套绑缚在被测者的胳膊上并充气到一定压力后打开阀门放气，放气到指定压力；

b.维持该指定压力在DT秒内不变(即定时中断)，在DT秒(本实施例为3秒)内压力传感器不间断地把当前的压力(压力单位为mmHg)和脉搏波信号转换为电信号，该电信号经放大后被送至A/D转换变成压力和脉搏波数字信号的时间序列，即对采样数据进行预处理；

c.对所述脉搏波时间序列进行信号变换，刚开始采集的100点数据由于放气刚刚结束，波形还不稳定，舍弃不用。选取之后的256个采样点进行CHIRPZ变换，利用FFT算法，计算出有限序列的DFT，也即是Z变换在单位园上的全部等间隔采样值。然而，在许多场合，并不一定需要计算全部频谱值，而仅需要对某一频带内的信号频谱坐较密集的分析。另外，采样也不一定局限于单位圆上，而需要计算出某一螺旋线上的等角度间隔的采样值。CHIRPZ变换就是利用FFT算法快速计算螺旋线采样的算法，具体算法如下：

输入N点采样值：x[n]                  n＝0，1，...，N-1

螺旋线的采样点：Z_k＝AW^-k    k＝0，1，…，M-1

其中 $A=A_0{e}^{{\mathrm{j\θ}}_0}$

A决定起始角度，W决定螺旋线的走向，₀表示螺旋线上的间隔角度

(1)选择FFT的点数L，L应满足

                       L≥N+M-1

且L＝2^v    v为正整数

实际运算中选择N＝M＝256点

(2)截取一段h(n)

$\begin{array}{cc}h\left(n\right)=W^{n^2/n}& -(L-M)\≤n\≤M-1\end{array}$

(3)由h(n)构成主值序列 h(n)

$\stackrel{\‾}{h}\left(n\right)=\left\{\begin{array}{cc}h\left(n\right)=W^{-n^2/2}& 0\≤n\≤M-1\\ h(n-L)=W^{-{(n-L)}^2/2}& M\≤n\≤L-1\end{array}\right.$

(4)计算H(k)＝FFT[ h(n)]，L点

(5)对x(n)加权并补零

$g\left(n\right)=\left\{\begin{array}{cc}x\left(n\right)\·A^{-n}{W}^{-n^2/2}& 0\≤n\≤N-1\\ 0& N\≤n\≤L-1\end{array}\right.$

(6)计算G(k)＝FFT[g(n)]

(7)计算Y(k)＝H(k)G(k)

(8)计算y(n)＝IFFT[Y(k)]

(9)计算 $\begin{array}{cc}X\left(Z_k\right)=W^{n^2/2}\·y\left(n\right)& 0\≤n\≤M-1\end{array}$

经过上述对N点采样点进行CHIRPZ变换，得到M点复数X(Z_k)，这里k＝0，1，...，M-1，M点复数的模代表某一频谱范围内的的信号幅度序列。因此，对256点脉搏波信号采样值进行CHIRPZ变换，得到对应某一频率范围内的脉搏波信号幅度序列，然后根据信号幅度序列计算当前压力下的最大脉搏波幅度。如图4所示。图4右边是经CHIRPZ变换后波幅随频率的变化曲线，图4(a)、(b)、(c)都有一个明显的主峰，该主峰就是最大波幅，该最大波幅对应的频率值与对应左边的原始数据是一致的，右边的幅频曲线若干个明显比最大波幅小的次峰值，这些次峰值是由干扰信号所造成的，我们只分析有物理意义的最大波幅及所对应的频率，而把若干次峰值置之不理，这样做就等于排除了干扰信号对测量结果的影响。

通过信号变换获得最大幅度及对应的频率值，计算最大幅度和幅度对应的频率，最大幅度作为当前压力下的幅度，最大幅度对应的频率值作为当前压力下的脉率值。继续放气至下一个指定压力(下降至下一个压力台阶)，返回到步骤b；

d.步骤b和步骤c被循环执行NUM次后，本实施例NUM取8，即获得8个在不同压力台阶上采样的脉搏波幅度序列和压力序列，对这两个数据序列进行归一化、插值、加权平均处理，然后用这两个经处理的数据序列构成脉搏波幅度和压力的关系曲线；以保证测量结果更加准确，确保此基准脉搏波序列符合比率法所需要满足的一些特性。

e.利用所述脉搏波幅度和压力的关系曲线计算平均压Mean及对应的脉搏波幅度P0，然后利用比率法进一步计算出所要测量的收缩压SYS和舒张压DIA对应的脉搏波幅度PS和PD，其中PS＝α*P0，PD＝β*P0，然后根据PS和PD并利用脉搏波幅度和压力的关系曲线求出收缩压和舒张压。

上述步骤a至步骤e的测量和计算流程如图3所示。

比率法基于如下理论：当袖带压等于舒张压或收缩压时，所对应的振荡波波幅，与振荡波中波幅最大值的比值，是一个相对恒定的比率。只要先找出波幅最大值，再根据确定的比率，就可找出相应的舒张压或收缩压对应的振荡波的波幅，进而找到与之对应的舒张压或收缩压的压力值。在本发明方法中，比率值的确定是很重要的。一般情况下，这个比率与仪器设备的特性有关，是通过大量实验获得的经验系数。在本实施例中，α的取值范围为0.2至0.9，β的取值范围为0.2至0.8。

图5是本发明装置系统实现的流程图，简述如下：

(1)上电后，完成硬件初始化，系统自检，变量初始化，进入主程序；

(2)设定定时中断，在定时中断中采集压力信号和脉搏波信号；

(3)接受上位机指令，开启一次测量，充气至指定压力，然后放气，放气到指定压力，维持指定压力一段时间(3秒)，启动信号处理模块和比率法分析模块，计算搜索该压力下的脉搏波信号幅度，然后再次放气下降一个压力台阶，再次启动信号处理模块和比率法分析模块，计算搜索该压力下的脉搏波信号幅度；

(4)经过8个压力台阶的采样计算，获得了脉搏波信号幅度和压力的关系曲线，利用该脉搏波信号幅度和压力的关系曲线，分别计算平均压，收缩压，舒张压。

(5)其他辅助功能，如漏气检测功能，压力校准功能等。

经过试验室测试证明，在正常情况下模拟器脉率设定在40BPM，80BPM和240BPM，血压为120/90/80下，信号变换算法计算结果和常规算法一致；在使用BIO-TECK模拟器，模拟10级肌颤情况下，信号变换的算法效果更好。两种算法测量结果的比较如下表：

模拟器设定	信号变换算法计算结果	常规算法计算结果	比较结果
				Cufflink模拟器120/90/80 40BPM 100％增益	121/92/8240BPM	122/92/8240BPM	两种算法一致
Cufflink模拟器120/90/80 80BPM 100％增益	120/90/8380BPM	118/90/8280BPM	两种算法一致
				Cufflink模拟器120/90/80 240BPM100％	122/92/82238BPM	122/92/83238BPM	两种算法一致
BIO-TEK模拟器肌颤10级120/93/8080BPM	120/94/8382BPM	130/99/8582BPM	信号变换算法更接近模拟器结果

Claims (9)

1、一种利用信号变换计算血压的方法，其特征在于该方法包括如下步骤：

a.把带有压力传感器的袖套绑缚在被测者的胳膊上并充气到一定压力后打开阀门放气，放气到指定压力；

b.维持该指定压力在DT秒内不变，在DT秒内压力传感器不间断地把当前的压力和脉搏波信号转换为电信号，该电信号经放大后被送至A/D转换变成压力和脉搏波数字信号的时间序列；

c.对所述脉搏波时间序列进行信号变换，由变换后的信号计算最大幅度及对应的频率值，通过最大幅度值求取当前压力下的幅度值；最大幅度对应的频率值作为当前压力下的脉率值，继续放气至下一个指定压力，返回到步骤b；

d.步骤b和步骤c被循环执行NUM次，获得NUM个在不同压力台阶上采样的脉搏波幅度序列和压力序列，用这两个数据序列构成脉搏波幅度和压力的关系曲线；

e.利用所述脉搏波幅度和压力的关系曲线计算平均压Mean及对应的脉搏波幅度P0，然后利用比率法进一步计算出所要测量的收缩压SYS和舒张压DIA对应的脉搏波幅度PS和PD，其中PS＝α*P0，PD＝β*P0，然后根据PS和PD并利用脉搏波幅度和压力的关系曲线求出收缩压和舒张压。

2、根据权利要求1所述的利用信号变换计算血压的方法，其特征在于：所述步骤e中比率法的两个比率α和β是通过实验获得的经验系数，α的取值范围为0.2至0.9，β的取值范围为0.2至0.8。

3、根据权利要求2所述的利用信号变换计算血压的方法，其特征在于：所述对脉搏波时间序列进行信号变换是指对脉搏波时间序列进行CHIRPZ变换，通过变换后的信号序列计算最大幅度及对应的频率值，最大幅度作为当前压力下的幅度值。

4、根据权利要求1、2或3所述的利用信号变换计算血压的方法，其特征在于：所述NUM的取值范围为6至10，所述DT的取值范围为3至10秒。

5、根据权利要求4所述的利用信号变换计算血压的方法，其特征在于：对所述步骤d获得的NUM个在不同压力台阶上采样的的脉搏波幅度序列和压力序列进行归一化、插值、加权平均处理，使用处理后的脉搏波幅度和压力序列构成脉搏波幅度和压力的关系曲线。

6、一种利用信号变换计算血压的装置，包括：

数据采样单元，包括位于袖套内的压力传感器、压力放大电路和脉搏波放大电路以及A/D转换；所述压力传感器把当前的压力和脉搏波信号转换为电信号，该电信号经压力放大电路和脉搏波放大电路放大后被送至A/D转换变成压力和脉搏波的数字信号；

主控单元，包括含有计算机软件的微处理器和控制电路，用于根据压力信号控制袖套放气的时间和速度以及根据算法软件计算被测人体的血压和脉率；

充放气单元，包括与所述袖套气路连接的充气泵、快放阀和慢放阀，用于接受主控单元的控制信号对所述袖套实施自动的充、放气操作；其特征在于还包括：

信号处理模块，用于对A/D转换输出的多个压力台阶下的脉搏波的数字信号进行信号变换，得到被分析数据的幅频特性从而进一步得到与被测人体血压紧密相关的波幅最大值；

比率法分析模块，用于调取所述信号处理模块输出的多个压力台阶下的波幅最大值以及A/D转换输出的多个压力台阶的压力值，以形成脉搏波幅度和压力的关系曲线，根据该关系曲线找出测量过程的波幅最大值，再根据系统确定的比率计算和输出被测人体的血压和脉率。

7、根据权利要求6所述的利用信号变换计算血压的装置，其特征在于：所述信号变换是指对脉搏波时间序列进行CHIRPZ变换。

8、根据权利要求6所述的利用信号变换计算血压的装置，其特征在于：所述比率法分析模块还包括：归一化模块、插值模块和加权平均模块，用于在计算和输出被测人体的血压和脉率之前对脉搏波幅度和压力的关系曲线进行归一化、插值和加权平均处理。

9、根据权利要求6、7或8所述的利用信号变换计算血压的装置，其特征在于：所述装置还包括与微处理器双向连接的辅助功能模块，用于系统复位、漏气检测、压力校准。

Images