CN102389313B - 一种方波调制光电容积脉搏波测量方法 - Google Patents

Abstract

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

Description

一种方波调制光电容积脉搏波测量方法

技术领域

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

背景技术

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

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

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

发明内容

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

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

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

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

本发明提供了一种方波调制光电容积脉搏波测量方法，所述方法包括以下步骤： 

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

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

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

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

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

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

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

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

附图说明

图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：运算放大器； 

PY口：I/O口。 

具体实施方式

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

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

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

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

$A=-1g\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，则动脉充盈时的吸光度和动脉收缩时的吸光度差值ΔA为： 

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

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

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

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

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

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

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

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

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

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

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

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

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

为简便说明起见，以4种波长发光二极管2为例进行说明，假定λ1、λ2、λ3和λ4波长的发光二极管的驱动方波频率分别为8倍f、4倍f、2倍f、1倍f。假定模数转换器5的采样频率为f_S，且f_S＝2f，并保证在λ1驱动信号高、低电平中间采样。 

数字信号序列 

可以表示为： 

$D_i^t=D_i^{\mathrm{\λ}1}+D_i^{\mathrm{\λ}2}+D_i^{\mathrm{\λ}3}+D_i^{\mathrm{\λ}4}+D_i^B---\left(7\right)$

其中， 和 

分别为波长λ1、λ2、λ3和λ4的光电容积脉搏波， 

为背景光和光敏器件3的暗电流、电流/电压转换放大器4的失调电压的总和 信号(简称背景信号)。 

假定采样频率f_S远高于调制方波信号和背景光的变化频率，在最低驱动信号频率的一个周期可以近似认为各路方波信号的幅值和背景光信号的幅值不变。以最前16个采样数据为例： 

$\begin{array}{c}\left\{\begin{array}{c}{D}_1^{\mathrm{\λ}1}=D_3^{\mathrm{\λ}1}=D_5^{\mathrm{\λ}1}=D_7^{\mathrm{\λ}1}=D_9^{\mathrm{\λ}1}=D_{11}^{\mathrm{\λ}1}=D_{13}^{\mathrm{\λ}1}=D_{15}^{\mathrm{\λ}1}=D_A^{\mathrm{\λ}1}\\ D_2^{\mathrm{\λ}1}=D_4^{\mathrm{\λ}1}{=D}_6^{\mathrm{\λ}1}=D_8^{\mathrm{\λ}1}=D_{10}^{\mathrm{\λ}1}=D_{12}^{\mathrm{\λ}1}=D_{14}^{\mathrm{\λ}1}=D_{16}^{\mathrm{\λ}1}=0\\ D_1^{\mathrm{\λ}2}=D_2^{\mathrm{\λ}2}=D_5^{\mathrm{\λ}2}=D_6^{\mathrm{\λ}2}=D_9^{\mathrm{\λ}2}=D_{10}^{\mathrm{\λ}2}=D_{13}^{\mathrm{\λ}2}=D_{14}^{\mathrm{\λ}2}=D_A^{\mathrm{\λ}2}\\ D_3^{\mathrm{\λ}2}=D_4^{\mathrm{\λ}2}=D_7^{\mathrm{\λ}2}=D_8^{\mathrm{\λ}2}=D_{11}^{\mathrm{\λ}2}=D_{12}^{\mathrm{\λ}2}=D_{15}^{\mathrm{\λ}2}=D_{16}^{\mathrm{\λ}2}=0\\ D_1^{\mathrm{\λ}3}=D_2^{\mathrm{\λ}3}=D_3^{\mathrm{\λ}3}=D_4^{\mathrm{\λ}3}=D_9^{\mathrm{\λ}3}=D_{10}^{\mathrm{\λ}3}=D_{11}^{\mathrm{\λ}3}=D_{12}^{\mathrm{\λ}3}=D_A^{\mathrm{\λ}3}\\ D_5^{\mathrm{\λ}3}=D_6^{\mathrm{\λ}3}=D_7^{\mathrm{\λ}3}=D_8^{\mathrm{\λ}3}=D_{13}^{\mathrm{\λ}3}=D_{14}^{\mathrm{\λ}3}=D_{15}^{\mathrm{\λ}3}=D_{16}^{\mathrm{\λ}3}=0\\ D_1^{\mathrm{\λ}4}=D_2^{\mathrm{\λ}4}=D_3^{\mathrm{\λ}4}=D_4^{\mathrm{\λ}4}=D_5^{\mathrm{\λ}4}=D_6^{\mathrm{\λ}4}=D_7^{\mathrm{\λ}4}=D_8^{\mathrm{\λ}4}=D_A^{\mathrm{\λ}4}\\ D_9^{\mathrm{\λ}4}=D_{10}^{\mathrm{\λ}4}=D_{11}^{\mathrm{\λ}4}=D_{12}^{\mathrm{\λ}4}=D_{13}^{\mathrm{\λ}4}=D_{14}^{\mathrm{\λ}4}=D_{15}^{\mathrm{\λ}4}=D_{16}^{\mathrm{\λ}4}=0\\ D_1^B=D_2^B=D_3^B=D_4^B=D_5^B=D_6^B=D_7^B=D_8^B=D_9^B=D_{10}^B=D_{11}^B=D_{12}^B=D_{13}^B=D_{14}^B\end{array}\right.---\left(8\right)\\ {=D}_{15}^B=D_{16}^B=D_A^B\end{array}$

其中， 

和 

分别为波长λ1、λ2、λ3和λ4的光信号和背景信号的幅值。 

换言之，以顺序每16个数字信号为一组进行运算： 

$D_{16n+1}-D_{16n+2}+D_{16n+3}-D_{16n+4}+D_{16n+5}-D_{16n+6}+D_{16n+7}-D_{16n+8}+D_{16n+9}-D_{16n+10}+$

$D_{16n+11}-D_{16n+12}+D_{16n+13}-D_{16n+14}+D_{16n+15}-D_{16n+16}=8D_{\mathrm{An}}^{\mathrm{\λ}1},n=\mathrm{0,1,2}......---\left(9\right)$

即得到8倍的波长λ1的光电容积脉搏波 

而且完全消除了背景信号 

的影响。 

$D_{16n+1}+D_{16n+2}-D_{16n+3}-D_{16n+4}+D_{16n+5}+D_{16n+6}-D_{16n+7}-D_{16n+8}+D_{16n+9}+D_{16n+10}-$

$D_{16n+11}-D_{16n+12}+D_{16n+13}+D_{16n+14}-D_{16n+15}-D_{16n+16}=8D_{\mathrm{An}}^{\mathrm{\λ}2},n=\mathrm{0,1,2}......---\left(10\right)$

即得到8倍的波长λ2的光电容积脉搏波 而且完全消除了背景信号 

的影响。 

$D_{16n+1}+D_{16n+2}+D_{16n+3}+D_{16n+4}-D_{16n+5}-D_{16n+6}-D_{16n+7}-D_{16n+8}+D_{16n+9}+D_{16n+10}+$

$D_{16n+11}+D_{16n+12}-D_{16n+13}-D_{16n+14}-D_{16n+15}-D_{16n+16}=8D_{\mathrm{An}}^{\mathrm{\λ}3},n=\mathrm{0,1,2}......---\left(11\right)$

即得到8倍的波长λ3的光电容积脉搏波 

而且完全消除了背景信号 

的影响。 

$D_{16n+1}+D_{16n+2}+D_{16n+3}+D_{16n+4}+D_{16n+5}+D_{16n+6}+D_{16n+7}+D_{16n+8}-D_{16n+9}-D_{16n+10}-$

$D_{16n+11}-D_{16n+12}-D_{16n+13}-D_{16n+14}-D_{16n+15}-D_{16n+16}=8D_{\mathrm{An}}^{\mathrm{\λ}4},n=\mathrm{0,1,2}......---\left(12\right)$

即得到8倍的波长λ4的光电容积脉搏波 

而且完全消除了背景信号 

的影响。 

分别计算波长λ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； 

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

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

如图5所示，一种方波调制光电容积脉搏波测量装置采用了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频率的二倍速度将预设幅值电压信号转换成数字信号通过PY口送入到微处理器1。数字信号在微处理器1先分离出不同波长光电容积脉搏波：每顺序获取的16个数字信号为一组，按照 

$D_{16n+1}-D_{16n+2}+D_{16n+3}-D_{16n+4}+D_{16n+5}-D_{16n+6}+D_{16n+7}-D_{16n+8}+D_{16n+9}-D_{16n+10}+$

$D_{16n+11}-D_{16n+12}+D_{16n+13}-D_{16n+14}+D_{16n+15}-D_{16n+16}=8D_{\mathrm{An}}^{\mathrm{\λ}1},n=\mathrm{0,1,2}......$

$D_{16n+1}+D_{16n+2}-D_{16n+3}-D_{16n+4}+D_{16n+5}+D_{16n+6}-D_{16n+7}-D_{16n+8}+D_{16n+9}+D_{16n+10}-$

$D_{16n+11}-D_{16n+12}+D_{16n+13}+D_{16n+14}-D_{16n+15}-D_{16n+16}=8D_{\mathrm{An}}^{\mathrm{\λ}2},n=\mathrm{0,1,2}......$

$D_{16n+1}+D_{16n+2}+D_{16n+3}+D_{16n+4}-D_{16n+5}-D_{16n+6}-D_{16n+7}-D_{16n+8}+D_{16n+9}+D_{16n+10}+$

$D_{16n+11}+D_{16n+12}-D_{16n+13}-D_{16n+14}-D_{16n+15}-D_{16n+16}=8D_{\mathrm{An}}^{\mathrm{\λ}3},n=\mathrm{0,1,2}......$

$D_{16n+1}+D_{16n+2}+D_{16n+3}+D_{16n+4}+D_{16n+5}+D_{16n+6}+D_{16n+7}+D_{16n+8}-D_{16n+9}-D_{16n+10}-$

$D_{16n+11}-D_{16n+12}-D_{16n+13}-D_{16n+14}-D_{16n+15}-D_{16n+16}=8D_{\mathrm{An}}^{\mathrm{\λ}4},n=\mathrm{0,1,2}......$

即分别得到8倍的波长λ1、λ2、λ3和λ4的光电容积脉搏波 

和 

而且完全消除了背景信号 

的影响。 

得到各个波长的光电容积脉搏波，据此计算出λ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倍比率关系的方波驱动至少2种发光二极管；

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

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

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

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

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

其中，所述根据所述光电容积脉搏波获取谷值和峰值的操作具体为：

λ1、λ2、λ3和λ4波长的发光二极管的驱动方波频率分别为8倍f、4倍f、2倍f、1倍f；模数转换器的采样频率为f_S，且f_S＝2f，并保证在λ1驱动信号高、低电平中间采样；

数字信号序列

表示为：

$D_i^t=D_i^{\mathrm{\λ}1}+D_i^{\mathrm{\λ}2}+D_i^{\mathrm{\λ}3}+D_i^{\mathrm{\λ}4}+D_i^B$

其中，和

分别为波长λ1、λ2、λ3和λ4的光电容积脉搏波，

为背景光和光敏器件的暗电流、电流/电压转换放大器的失调电压的总和信号；

采样频率f_S远高于调制方波信号和背景光的变化频率，在最低驱动信号频率的一个周期近似认为各路方波信号的幅值和背景光信号的幅值不变，以最前16个采样数据为例：

$\left\{\begin{array}{c}{D}_1^{\mathrm{\λ}1}=D_3^{\mathrm{\λ}1}=D_5^{\mathrm{\λ}1}=D_7^{\mathrm{\λ}1}=D_9^{\mathrm{\λ}1}=D_{11}^{\mathrm{\λ}1}=D_{13}^{\mathrm{\λ}1}=D_{15}^{\mathrm{\λ}1}=D_A^{\mathrm{\λ}1}\\ D_2^{\mathrm{\λ}1}=D_4^{\mathrm{\λ}1}=D_6^{\mathrm{\λ}1}=D_8^{\mathrm{\λ}1}=D_{10}^{\mathrm{\λ}1}=D_{12}^{\mathrm{\λ}1}=D_{14}^{\mathrm{\λ}1}=D_{16}^{\mathrm{\λ}1}=0\\ D_1^{\mathrm{\λ}2}=D_2^{\mathrm{\λ}2}=D_5^{\mathrm{\λ}2}=D_6^{\mathrm{\λ}2}=D_9^{\mathrm{\λ}2}=D_{10}^{\mathrm{\λ}2}=D_{13}^{\mathrm{\λ}2}=D_{14}^{\mathrm{\λ}2}=D_A^{\mathrm{\λ}2}\\ D_3^{\mathrm{\λ}2}=D_4^{\mathrm{\λ}2}=D_7^{\mathrm{\λ}2}=D_8^{\mathrm{\λ}2}=D_{11}^{\mathrm{\λ}2}=D_{12}^{\mathrm{\λ}2}=D_{15}^{\mathrm{\λ}2}=D_{16}^{\mathrm{\λ}2}=0\\ D_1^{\mathrm{\λ}3}=D_2^{\mathrm{\λ}3}=D_3^{\mathrm{\λ}3}=D_4^{\mathrm{\λ}3}=D_9^{\mathrm{\λ}3}=D_{10}^{\mathrm{\λ}3}=D_{11}^{\mathrm{\λ}3}=D_{12}^{\mathrm{\λ}3}=D_A^{\mathrm{\λ}3}\\ \\ \\ \\ \\ \end{array}\right.---\left(8\right)$

$\left[\begin{array}{c}{D}_5^{\mathrm{\λ}3}=D_6^{\mathrm{\λ}3}=D_7^{\mathrm{\λ}3}=D_8^{\mathrm{\λ}3}=D_{13}^{\mathrm{\λ}3}=D_{14}^{\mathrm{\λ}3}=D_{15}^{\mathrm{\λ}3}=D_{16}^{\mathrm{\λ}3}=0\\ D_1^{\mathrm{\λ}4}=D_2^{\mathrm{\λ}4}=D_3^{\mathrm{\λ}4}=D_4^{\mathrm{\λ}4}=D_5^{\mathrm{\λ}4}=D_6^{\mathrm{\λ}4}=D_7^{\mathrm{\λ}4}=D_8^{\mathrm{\λ}4}=D_A^{\mathrm{\λ}4}\\ D_9^{\mathrm{\λ}4}=D_{10}^{\mathrm{\λ}4}=D_{11}^{\mathrm{\λ}4}=D_{12}^{\mathrm{\λ}4}=D_{13}^{\mathrm{\λ}4}=D_{14}^{\mathrm{\λ}4}=D_{15}^{\mathrm{\λ}4}=D_{16}^{\mathrm{\λ}4}=0\\ D_1^B=B_2^B=D_3^B=D_4^B=D_5^B=D_6^B=D_7^B=D_8^B=D_9^B=D_{10}^B=D_{11}^B=D_{12}^B=D_{13}^B{=D}_{14}^B\\ =D_{15}^B=D_{16}^B=D_A^B\end{array}\right]$

其中，

和分别为波长λ1、λ2、λ3和λ4的光信号和背景信号的幅值；

以顺序每16个数字信号为一组进行运算：

$\begin{array}{c}{D}_{16n+1}-D_{16n+2}+D_{16n+3}-D_{16n+4}+D_{16n+5}-D_{16n+6}+D_{16n+7}-D_{16n+8}+D_{16n+9}-D_{16n+10}+\\ D_{16n+11}-D_{16n+12}+D_{16n+13}-D_{16n+14}+D_{16n+15}-D_{16n+16}=8D_{\mathrm{An}}^{\mathrm{\λ}1}\begin{array}{}\end{array}n=\mathrm{0,1,2}......\end{array}---\left(9\right)$

即得到8倍的波长λ1的光电容积脉搏波

消除了背景信号

的影响；

$\begin{array}{c}{D}_{16n+1}+D_{16n+2}-D_{16n+3}-D_{16n+4}+D_{16n+5}+D_{16n+6}-D_{16n+7}-D_{16n+8}+D_{16n+9}+D_{16n+10}-\\ D_{16n+11}-D_{16n+12}+D_{16n+13}+D_{16n+14}-D_{16n+15}-D_{16n+16}=8D_{\mathrm{An}}^{\mathrm{\λ}2}\begin{array}{}\end{array}n=\mathrm{0,1,2}......\end{array}---\left(10\right)$

即得到8倍的波长λ2的光电容积脉搏波消除了背景信号

的影响；

$\begin{array}{c}{D}_{16n+1}+D_{16n+2}+D_{16n+3}+D_{16n+4}-D_{16n+5}-D_{16n+6}-D_{16n+7}-D_{16n+8}+D_{16n+9}+D_{16n+10}+\\ D_{16n+11}+D_{16n+12}-D_{16n+13}-D_{16n+14}-D_{16n+15}-D_{16n+16}=8D_{\mathrm{An}}^{\mathrm{\λ}3}\begin{array}{}\end{array}n=\mathrm{0,1,2}......\end{array}---\left(11\right)$

即得到8倍的波长λ3的光电容积脉搏波

消除了背景信号

的影响；

$\begin{array}{c}{D}_{16n+1}+D_{16n+2}+D_{16n+3}+D_{16n+4}+D_{16n+5}+D_{16n+6}+D_{16n+7}+D_{16n+8}-D_{16n+9}-D_{16n+10}-\\ D_{16n+11}-D_{16n+12}-D_{16n+13}-D_{16n+14}-D_{16n+15}-D_{16n+16}=8D_{\mathrm{An}}^{\mathrm{\λ}4}\begin{array}{}\end{array}n=\mathrm{0,1,2}......\end{array}---\left(12\right)$

即得到8倍的波长λ4的光电容积脉搏波

消除了背景信号

的影响；

分别计算波长λ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。

Images