CN102258365B - 一种正弦波调制光电容积脉搏波测量装置和测量方法 - Google Patents

Abstract

本发明公开了一种正弦波调制光电容积脉搏波测量装置和测量方法，微处理器输出不同频率且成2倍比率关系的正弦波，正弦波驱动至少4种发光二极管，发光二极管发出的光经被测手指后被光敏器件接收，光敏器件转换成电压信号，电压信号经电流/电压转换放大器转换成预设幅值电压信号，模数转换器将预设幅值电压信号转换成数字信号，微处理器对数字信号进行处理，获取光电容积脉搏波及其谷值和峰值，通过谷值和峰值得到光谱值；方法包括：微处理器对数字信号进行锁相计算、分离处理得到光电容积脉搏波并消除背景光的干扰；根据光电容积脉搏波获取谷值和峰值；对谷值和峰值进行计算得到吸光度差值，通过吸光度差值获取光谱值。电路简单、数字信号简洁。

Description

一种正弦波调制光电容积脉搏波测量装置和测量方法

技术领域

本发明涉及一种正弦波调制光电容积脉搏波测量装置和测量方法。

背景技术

光电容积脉搏波(Photo Plethysmo Graphy，以下简称PPG)是一种重要的生理信号，广泛地应用对心血管系统和血液成分进行分析。如对血氧饱和度的测量中就是采用2种或2种以上的LED(发光二极管)测量PPG而实现的。在这些测量中通常采用时分方式采集PPG并消除背景光的干扰。

发明人在实现本发明的过程中发现，现有技术中至少存在以下缺点和不足：

现有的多波长PPG的测量方法存在电路结构复杂、器件和工艺要求高、调试困难、可靠性低、计算量大以及结果不够准确等缺点。

发明内容

本发明要解决的技术问题在于提供一种正弦波调制光电容积脉搏波测量装置和测量方法，该测量装置和测量方法可以实现高精度测量，且电路结构简单、器件和工艺要求低、调试容易、可靠性高、计算量小等优点，详见下文描述：

一种正弦波调制光电容积脉搏波测量装置，所述光电容积脉搏波测量装置包括：微处理器、至少4种发光二极管、光敏器件、电流/电压转换放大器和模数转换器，

所述微处理器输出不同频率且成2倍比率关系的正弦波，所述正弦波驱动至少所述4种发光二极管，所述发光二极管发出的光经被测手指后被所述光敏器件接收，所述光敏器件转换成电压信号，所述电压信号经所述电流/电压转换放大器转换成预设幅值电压信号，所述模数转换器将所述预设幅值电压信号转换成数字信号，所述微处理器对所述数字信号进行处理获取光电容积脉搏波及其谷值和峰值，通过所述谷值和所述峰值得到光谱值。

所述微处理器采用MCU、ARM、DSP或FPGA中的任意一种。

一种正弦波调制光电容积脉搏波测量方法，所述方法包括以下步骤：

(1)微处理器采用不同频率且成2倍比率关系的正弦波驱动至少4种发光二极管；

(2)所述发光二极管发出的光经过被测手指后由光敏器件接收转换成电压信号，所述电压信号经过电流/电压转换放大器放大成预设幅值电压信号；

(3)所述预设幅值电压信号经模数转换器转换成数字信号送入所述微处理器；

(4)所述微处理器对所述数字信号进行锁相计算、分离处理得到光电容积脉搏波；

(5)根据所述光电容积脉搏波获取谷值和峰值；

(6)对所述谷值和所述峰值进行计算得到吸光度差值，通过所述吸光度差值获取光谱值。

步骤(4)中所述微处理器对所述数字信号进行锁相计算、分离处理得到光电容积脉搏波并消除背景光的干扰，具体包括：

1)假设所述微处理器控制所述模数转换器以驱动所述发光二极管的最高频率f_max的4M倍速度对所述数字信号进行采样f_s＝4M×f_max获取采样信号x(m)，其中M为大于等于1的正整数；

$x\left(m\right)=\underset{i=0}{\overset{M-1}{\mathrm{\Σ}}}x(4\mathrm{lm}+i),$ l＝0，1，2，......

2)所述微处理器将采样频率f_s＝4_M×f_max下抽样M倍到4×f_max；

3)根据所述采样信号x(m)、所述微处理器产生的正交参考序列y_s(k)和y_c(k)计算两个正交分量R_S和R_C；

$R_s=\frac{1}{Q}\underset{m=0}{\overset{N-1}{\mathrm{\Σ}}}[x(4m+1)-x(4m+3)];$

$R_C=\frac{1}{Q}\underset{m=0}{\overset{N-1}{\mathrm{\Σ}}}[x\left(4m\right)-x(4m+2)];$

4)根据正交分量R_S和R_C通过低通滤波获取所述数字信号的幅值A；

β为下抽样M倍而取的某个常数值；

5)对所述采样信号x(m)下抽样2倍，重复第3)步、4)步继续计算频率低一半所述数字信号的幅值A，直到计算、分离完全部频率的光电容积脉搏波。

本发明提供的一种正弦波调制光电容积脉搏波测量装置和测量方法，与现有技术相比具有如下的优点：

本发明依据朗伯-比尔定律，采用正弦波频分调制和数字解调技术设计一种正弦波调制多波长发光二极管的光电容积脉搏波及其谷值和峰值，通过谷值和峰值得到光谱值的装置和测量方法，具有测量精确、电路简单、无需调试、工艺性好以及成本低廉的特点。

附图说明

图1为本发明提供的计算吸光度的原理示意图；

图2为本发明提供的一种正弦波调制光电容积脉搏波测量装置的结构示意图；

图3为本发明提供的分离不同波长光电容积脉搏波的示意图；

图4为本发明提供的一种正弦波调制光电容积脉搏波测量方法的流程图；

图5为本发明提供的一种正弦波光电容积脉搏波测量装置的另一结构示意图。

附图中各标号所代表的部件列表如下：

1：微处理器；           2：发光二极管；

3：光敏器件；           4：电流/电压转换放大器；

5：模数转换器；         PX.1：I/O口；

PX.2：I/O口；           PX.n：I/O口；

PX.3：I/O口；           PX.4：I/O口；

R1：第一电阻；          VCC：电源；

R2：第二电阻；          R3：第三电阻；

R4：第四电阻；          R5：第五电阻；

R6：第六电阻；          C1：第一电容；

C2：第二电容；          D1：第一发光二极管；

D2：第二发光二极管；    D3：第三发光二极管；

D4：第四发光二极管；    A1：运算放大器；

P Y 口：I/O 口。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。

由于动脉的脉动现象，使血管中血流量呈周期性变化，而血液是高度不透明液体，因此脉搏搏动的变化必然引起吸光度的变化，如图1所示。

考虑动脉血管充盈度最低状态，来自光源的入射光没有被脉动动脉血液吸收，此时的出射光强I_max最强，可视为脉动动脉血液的入射光I；而动脉血管充盈度最高状态对应光电脉搏波的谷点，即脉动动脉血液作用最大的时刻，此时的出射光强I_min最弱，为脉动动脉血液的最小出射光强I。所以，通过记录动脉充盈至最大与动脉收缩至最小时的吸光度值，就可以消除皮肤组织、皮下组织等一切具有恒定吸收特点的人体成分对于吸光度的影响。

根据修正的朗伯-比尔定律，设I₀、I分别为入射光强和出射光强，α为分子消光系数，c为各成分浓度，l为光在组织中的平均光路长，G是由散射引起的光损失，则吸光度A可表示为：

$A=-\lg \frac{I}{I_0}=-2.303\mathrm{\αcl}+G---\left(1\right)$

设生物组织的吸收系数为μ_a，则μ_a＝αc，代入式(1)可得：

A＝-2.303μ_al+G    (2)

在光透射检测中，吸光度主要由被透射组织的吸收与散射构成，其中血液散射相对较小，可忽略不计。这样，G仅仅由除了脉动动脉血外的组织贡献，在测量过程中保持不变。设除脉动动脉血外的被透射组织共n层，第i层的吸收系数为μ_ti，动脉血的吸收系数为μ_ab，一个光电脉搏波周期上动脉充盈时最大光路长为l_max，动脉收缩时的最小光路长为l_min，则动脉充盈时吸光度A₁和动脉收缩时吸光度A₂可分别表示为：

$A_1=-2.303\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_{\mathrm{ti}}{l}_{\max }-2.303{\mathrm{\μ}}_{\mathrm{ab}}{l}_{\max }+G---\left(3\right)$

$A_2=-2.303\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_{\mathrm{ti}}{l}_{\min }-2.303{\mathrm{\μ}}_{\mathrm{ab}}{l}_{\min }+G---\left(4\right)$

设l为l_max与l_min之差。由于除了脉动动脉血液以外的其他组织基本稳定，不进行周期变化，因此该部分在动脉充盈和收缩时对吸光度没有影响，即式(3)和式(4)中的第一个分量相等。则动脉充盈时的吸光度和动脉收缩时的吸光度之差为：

ΔA＝A₁-A₂＝-2.303μ_ab(l_max-l_min)＝-2.303μ_abl    (5)

在上面的推导过程中，非脉动血液和各层组织的吸收和散射的吸光度分量都被消掉了，动脉充盈时和动脉收缩时的吸光度差值ΔA仅由动脉血的脉动吸收部分贡献，主要反映脉动的动脉血的吸收变化。在本质上相当于在被透射组织中、皮肤、肌肉以及静脉血液等除脉动动脉血液外的其他组织的影响都被去除了，只留下纯粹的脉动动脉血部分来进行吸光度差值ΔA的测量。这样一来，皮肤、骨骼和肌肉等个体差异的影响都被去除了。

设入射光强为I₀，动脉充盈时检测光强和动脉收缩时检测光强分别为I_min和I_max，则动脉充盈时的吸光度和动脉收缩时的吸光度差值为：

$\mathrm{\ΔA}=A_1-A_2=\lg \left(\frac{I_0}{I_{\min }}\right)-\lg \left(\frac{I_0}{I_{\max }}\right)=\lg \left(\frac{I_{\max }}{I_{\min }}\right)---\left(6\right)$

测量各个光电容积脉搏波的谷值I_min和峰值I_max即可得到光电容积脉搏波所对应的吸光度差值ΔA，可以得到由ΔA_λ1、ΔA_λ2……ΔA_λn组成的光谱值。

一种正弦波调制光电容积脉搏波测量装置，参见图2，该正弦波调制光电容积脉搏波测量装置包括：微处理器1、至少4种发光二极管2、光敏器件3、电流/电压转换放大器4和模数转换器5，

微处理器1输出不同频率且成2倍比率关系的正弦波，正弦波驱动至少4种发光二极管2，发光二极管2发出的光经被测手指后被光敏器件3接收，光敏器件3转换成电压信号，电压信号经电流/电压转换放大器4转换成预设幅值电压信号，模数转换器5将预设幅值电压信号转换成数字信号，微处理器1对数字信号进行处理，获取光电容积脉搏波及其谷值和峰值，通过谷值和峰值得到光谱值。

其中，发光二极管2的数量大于等于4。具体实现时，发光二极管2的数量根据实际应用中的需要进行设定，本发明实施例对此不做限制。

其中，预设幅值根据实际应用中的需要进行设定，具体实现时，本发明实施例对此不做限制。

其中，微处理器1可以采用MCU、ARM、DSP或FPGA中的任意一种。

一种正弦波调制光电容积脉搏波测量方法，参见图3和图4，该方法包括以下步骤：

101：微处理器1采用不同频率且成2倍比率关系的正弦波驱动至少4种发光二极管2；

102：发光二极管2发出的光经过被测手指后由光敏器件3接收转换成电压信号，电压信号经过电流/电压转换放大器4放大成预设幅值电压信号；

103：预设幅值电压信号经模数转换器5转换成数字信号送入微处理器1；

104：微处理器1对数字信号进行锁相计算、分离处理得到光电容积脉搏波并消除背景光的干扰：

该步骤具体为：

1、假设微处理器1控制模数转换器5以驱动发光二极管2的最高频率f_max的4M倍速度对数字信号进行采样f_s＝4M×f_max获取采样信号x(m)，其中M为大于等于1的正整数；

$x\left(m\right)=\underset{i=0}{\overset{M-1}{\mathrm{\Σ}}}x(4\mathrm{lm}+i),$ l＝0，1，2，......

2、微处理器1将采样频率f_s＝4M×f_max下抽样M倍到4×f_max；

3、根据采样信号x(m)、微处理器1产生的正交参考序列y_s(k)和y_c(k)计算两个正交分量R_S和R_C；

$R_s=\frac{1}{Q}\underset{m=0}{\overset{N-1}{\mathrm{\Σ}}}[x(4m+1)-x(4m+3)];$

$R_C=\frac{1}{Q}\underset{m=0}{\overset{N-1}{\mathrm{\Σ}}}[x\left(4m\right)-x(4m+2)];$

4、根据正交分量R_S和R_C通过低通滤波获取数字信号的幅值A，

β为下抽样M倍而取的某个常数值。

得到该数字信号的幅值A，也即分离出最高频率f_max的PPG信号。

5、对采样信号x(m)下抽样2倍，重复第3步、4步继续计算频率低一半数字信号的幅值A，直到计算、分离完全部频率的PPG信号。

105：根据光电容积脉搏波获取谷值和峰值；

106：对谷值和峰值进行计算得到吸光度差值，通过吸光度差值获取光谱值。

采用公式(6)计算各个吸光度差值ΔA_λ1、ΔA_λ2、......ΔA_λn，并由吸光度差值构成光谱值。

如图4所示，一种正弦波调制光电容积脉搏波测量装置采用了4种发光二极管2，微处理器1的四个I/O口PX.1、PX.2、PX.3和PX.4分别通过第一电阻R1、第二电阻R2、第三电阻R3和第四电阻R4驱动第一发光二极管D1、第二发光二极管D2、第三发光二极管D3和第四发光二极管D4，第一发光二极管D1、第二发光二极管D2、第三发光二极管D3和第四发光二极管D4发出的光透光被测手指被光敏传感器3所接收，光敏传感器3所接收的信号经过由运算放大器A1、第一电容C1、第二电容C2和第五电阻R5、第六电阻R6所组成的电流/电压转换放大器4转换成预设幅值电压信号，然后模数转换器5以最高驱动发光二极管2频率的4M倍速度将预设幅值电压信号转换成数字信号通过PY口送入到微处理器1。数字信号在微处理器1先计算、分离出不同波长光电容积脉搏波：得到各个波长的光电容积脉搏波，据此计算出波长λ1、λ2、λ3和λ4的光电容积脉搏波的谷值和峰值：I_minλ1、I_maxλ1、I_minλ2、I_maxλ2、I_minλ3、I_maxλ3、I_minλ4和I_maxλ4。再由I_minλ1、I_maxλ1、I_minλ2、I_maxλ2、I_minλ3、I_maxλ3、I_minλ4和I_maxλ4计算各个波长所对应的吸光度差值ΔA，可以得到由吸光度差值ΔA_λ1、ΔA_λ2……ΔA_λn组成的光谱值。

综上所述，本发明实施例提供了一种正弦波调制光电容积脉搏波测量装置和测量方法，本发明实施例依据朗伯-比尔定律，采用正弦波频分调制和数字解调技术设计一种正弦波调制多波长发光二极管的光电容积脉搏波及其谷值和峰值，通过谷值和峰值得到光谱值的装置和测量方法，具有测量精确、电路简单、无需调试、工艺性好以及成本低廉的特点。

本领域技术人员可以理解附图只是一个优选实施例的示意图，上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (1)

1.一种正弦波调制光电容积脉搏波测量装置的测量方法，其特征在于，所述方法包括以下步骤：

(1)微处理器采用不同频率且成2倍比率关系的正弦波驱动至少4种发光二极管；

(2)所述发光二极管发出的光经过被测手指后由光敏器件接收转换成电压信号，所述电压信号经过电流/电压转换放大器放大成预设幅值电压信号；

(3)所述预设幅值电压信号经模数转换器转换成数字信号送入所述微处理器；

(4)所述微处理器对所述数字信号进行锁相计算、分离处理得到光电容积脉搏波并消除背景光的干扰；

(5)根据所述光电容积脉搏波获取谷值和峰值；

(6)对所述谷值和所述峰值进行计算得到吸光度差值，通过所述吸光度差值获取光谱值；

其中，步骤(4)中所述微处理器对所述数字信号进行锁相计算、分离处理得到光电容积脉搏波并消除背景光的干扰，具体包括：

1）假设所述微处理器控制所述模数转换器以驱动所述发光二极管的最高频率f_max的4M倍速度对所述数字信号进行采样f_s=4M×f_max，获取采样信号x(m)，其中M为大于等于1的正整数；

$x\left(m\right)=\underset{i=0}{\overset{M-1}{\mathrm{\Σ}}}(4\mathrm{lm}+i),l=\mathrm{0,1,2},\·\·\·\·\·\·$

2）所述微处理器将采样频率f_s=4M×f_max下抽样M倍到4×f_max；

3）根据所述采样信号x(m)、所述微处理器产生的正交参考序列y_s(k)和y_c(k)计算两个正交分量R_S和R_C；

$R_S=\frac{1}{Q}\underset{m=0}{\overset{N-1}{\mathrm{\Σ}}}\[x(4m+1)-x(4m+3)\];$

$R_C=\frac{1}{Q}\underset{m=0}{\overset{N-1}{\mathrm{\Σ}}}[x\left(4m\right)-x(4m+2)];$

4）根据正交分量R_S和R_C通过低通滤波获取所述数字信号的幅值A；

β为下抽样M倍而取的某个常数值；

5）对所述采样信号x(m)下抽样2倍，重复第3）步、4）步继续计算频率低一半所述数字信号的幅值A，直到计算、分离完全部频率的光电容积脉搏波。

Images