专利名称：一种方波调制光电容积脉搏波测量装置和测量方法
技术领域：
本发明涉及一种方波调制光电脉搏容积波测量装置和测量方法。
背景技术：
光电容积脉搏波(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 1/0 口 ；
PX2 1/0 口 ；PX. η 1/0 口 ；
PX3 1/0 口 ；PX. 4 1/0 口 ；
Rl第一电阻；VCC 电源；
R2第二电阻；R3 第三电阻；
R4第四电阻；R5 第五电阻；
R6第六电阻；Cl 第一电容；
C2第二电容；Dl 第一发光二极管；
D2第二发光二极管；D3 第三发光二极管；
D4第四发光二极管；Al 运算放大器；
PY 口 1/0 口。
具体实施例方式为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。由于动脉的脉动现象，使血管中血流量呈周期性变化，而血液是高度不透明液体， 因此脉搏搏动的变化必然引起吸光度的变化，如图1所示。考虑动脉血管充盈度最低状态，来自光源的入射光没有被脉动动脉血液吸收，此时的出射光强Imax最强，可视为脉动动脉血液的入射光I ；而动脉血管充盈度最高状态对应光电脉搏波的谷点，即脉动动脉血液作用最大的时刻，此时的出射光强Imin最弱，为脉动动脉血液的最小出射光强I。所以，通过记录动脉充盈至最大与动脉收缩至最小时的吸光度值，就可以消除皮肤组织、皮下组织等一切具有恒定吸收特点的人体成分对于吸光度的影响。根据修正的朗伯-比尔定律，设I” I分别为入射光强和出射光强，α为分子消光系数，c为各成分浓度，1为光在组织中的平均光路长，G是由散射引起的光损失，则吸光度 A可表示为
A = -lg— = -2.303acl + G(1)
Iq设生物组织的吸收系数为μ a，则μ a = α c，代入式(1)可得A = -2. 303 μ al+G(2)在光透射检测中，吸光度主要由被透射组织的吸收与散射构成，其中血液散射相对较小，可忽略不计。这样，G仅仅由除了脉动动脉血外的组织贡献，在测量过程中保持不变。设除脉动动脉血外的被透射组织共η层，第i层的吸收系数为μ ti，动脉血的吸收系数为μ ab，一个光电脉搏波周期上动脉充盈时最大光路长为Imax，动脉收缩时的最小光路长为 Imin，则动脉充盈时吸光度A1和动脉收缩时吸光度A2可分别表示为
A = -2.303^ ^Zmax - 2.303&/max + G(3)
A = -2.303^^/mm -2303Uablmm +G(4)
i=\设1为Imax与Imin之差。由于除了脉动动脉血液以外的其他组织基本稳定，不进行周期变化，因此该部分在动脉充盈和收缩时对吸光度没有影响，即式C3)和式(4)中的第一个分量相等。则动脉充盈时的吸光度和动脉收缩时的吸光度之差为Δ A = A1-A2 = -2. 303 μ ab (Imax-Imin) = _2· 303 μ abl(5)在上面的推导过程中，非脉动血液和各层组织的吸收和散射的吸光度分量都被消掉了，动脉充盈时和动脉收缩时的吸光度差值ΔΑ仅由动脉血的脉动吸收部分贡献，主要反映脉动的动脉血的吸收变化。在本质上相当于在被透射组织中、皮肤、肌肉以及静脉血液等除脉动动脉血液外的其他组织的影响都被去除了，只留下纯粹的脉动动脉血部分来进行吸光度差值ΔΑ的测量。这样一来，皮肤、骨骼和肌肉等个体差异的影响都被去除了。设入射光强为Itl，动脉充盈时检测光强和动脉收缩时检测光强分别为Imin和Imax， 则动脉充盈时的吸光度和动脉收缩时的吸光度差值ΔΑ为
M = A1-A2=- lg(-^) = Igi^l(6)
minmax\ min J测量各个光电容积脉搏波的谷值Imin和峰值Imax即可得到光电容积脉搏波所对应的吸光度差值ΔΑ，可以得到由ΔΑλ1、ΔΑλ2……成的光谱值。一种方波调制光电容积脉搏波测量装置，参见图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、l倍f。假定模数转换器5的采样频率为fs，且fs = 2f，并保证在λ 1驱动信号高、低电平中间采样。数字信号序列Ai可以表示为
Dti = Df1 + Df2 + Df3 + Df4 + Df(7)其中，忧1、Df2、D尸和^4分别为波长λ 1、λ 2、λ 3禾Π λ 4的光电容积脉搏波， As为背景光和光敏器件3的暗电流、电流/电压转换放大器4的失调电压的总和信号(简称背景信号)。假定采样频率fs远高于调制方波信号和背景光的变化频率，在最低驱动信号频率的一个周期可以近似认为各路方波信号的幅值和背景光信号的幅值不变。以最前16个采样数据为例r DM = DM = DM = DM = DM = Dxi = Dxi = Dxi = Dxi
(MMl M 3 M 5 1^A
Dm = Dm = Dm = Dm = Dm = Dm = Dm =Dm =0
ljI ljA U6 J^s -L^10 I-J11-Ly16 yJ
Dx2 = Dx2 = Dx2 = Dx2 = Dx2 = Dx2 = Dx2 = Dx2 = Dx2
ljI ljI ljS U6 I-J9 -tvIO "13 "14 ljA
Df = Df = Df = Dl1 = D^ = Df22 = D1]2 = Df62 = 0 < Df3 = Df = Df = Df = Df = D^ = D^ = Df23 = Df(8)
Df = Df = Df = Df = Df33 = Df43 = D^ = D^ = 0 D1"4 = Df = D3"4 = Df = D》4 = D164 = Df = D^4 = Df Dm - ηΛ4 - ηΛ4 - ηΛ4 - ηΛ4 - ηΛ4 - ηΛ4 - ηΛ4 - η
"9 ~ ljV) - ljW - ljM _"13 _"14 _ "15 _ "16
\ DB _ dB _ dB _ dB _ dB _ dB _ dB _ dB _ dB _ dB _ dB _ dB _ dB _ dB 、“1 - ljI _ "3 - ljA _ "5 _ "6 - ljI _ "8 _ "9 _ Mo _ Ml ~ M2 _ "13 _ "14
—Db _ dB _ dB
--- LJa其中，Dj1、Df、Df、Dj4和叱分别为波长λ 1、λ 2、λ 3禾口 λ 4的光信号和
背景信号的幅值。换言之，以顺序每16个数字信号为一组进行运算Dl6n+l — Dl6n+2 + ~ ^ 6η+4 ^16n+5 _ ^16n+6 ^16n+7 _ ^16n+8 ^16n+9 _ ^16n+10 +
D16n+U ~Dl6n+l2 +D16n+13 -D16n+14 +D16n+15 -D^6=Wfn =0,1,2...... (9)即得到8倍的波长λ 1的光电容积脉搏波D=，而且完全消除了背景信号Af的影响。-Dl6n+l + Dl6n+2 ~ Dl6n+ ) — D16n+4 + Dl6n+i + Dl6n+6 — Dl6n+1 — Dl6n+i + D16n+9 + Dl6n+W —
Dl6n+ll ~Dl6n+u +D16n+13 +D16n+14 -D16n+15 -D16n+16 = ^Dfn =0,1,2...... (10)即得到8倍的波长λ 2的光电容积脉搏波Df，而且完全消除了背景信号Af的影响。Dl6n+l + Dl6n+2 +^ 6η+4 ~ ^16n+5 _ ^16n+6 _ ^16n+7 _ ^16n+8 ^16n+9 Μδη+ΙΟ +
D16n+U +D16n+U ~D16n+13 -D16n+14 -D16n+15 -UDfn =0,1,2...... (11)即得到8倍的波长λ 3的光电容积脉搏波Di3，而且完全消除了背景信号Af的影响。Dl6n+l + Dl6n+2 + ^ 6η+4 ^16n+5 ^16n+6 ^16n+7 ^16n+8 _ ^16n+9 _ ^16n+10 ~~
Αδη+Π _ Αδη+12 _ Αδη+13 _ Αδη+14 _ Αδη+15 _ Αδη+16 = ^^An W—0,1,2...... (12)即得到8倍的波长λ 4的光电容积脉搏波Di4，而且完全消除了背景信号Af的影响。分别计算波长λ 1、λ 2、λ 3和λ 4的光电容积脉搏波的谷值和峰值=Imin；u、Imax；u、
IfliinA 2Λ ImaxA 2 Λ IfliinA 3Λ ^maxA 3Λ IfliinA 4 和 ^maxA 4 ‘106:对谷值和峰值进行计算得到吸光度差值，通过吸光度差值获取光谱值。
采用公式(6)计算各个波长的吸光度差值ΔΑλ1、ΔΑλ2.......ΔΑλη，并由吸光度
差值构成光谱值。如图5所示，一种方波调制光电容积脉搏波测量装置采用了 4种发光二极管2，微处理器1的四个I/O 口 PX. UPX. 2、PX. 3和PX. 4分别通过第一电阻Rl、第二电阻R2、第三电阻R3和第四电阻R4驱动第一发光二极管Dl、第二发光二极管D2、第三发光二极管D3和第四发光二极管D4，第一发光二极管D1、第二发光二极管D2、第三发光二极管D3和第四发光二极管D4发出的光透光被测手指被光敏传感器3所接收，光敏传感器3所接收的信号经过由运算放大器Al、第一电容Cl、第二电容C2、第五电阻R5和第六电阻R6所组成的电流/ 电压转换放大器4转换成预设幅值电压信号，然后模数转换器5以最高驱动发光二极管2 频率的二倍速度将预设幅值电压信号转换成数字信号通过PY 口送入到微处理器1。数字信号在微处理器1先分离出不同波长光电容积脉搏波每顺序获取的16个数字信号为一组， 按照- Dl6n+2 + ^16^+3 _ ΜδΜ+4 η+5 ~ ^16η+6 ^\βη+ ~^\βη+9 ~+
-^\6η+\\ ~ ΜδΜ+12 +Dl6M+13 _ ^\6n+\4 ^\6n+\5 ~ ^\6n+\6 ~ ^^An/7—0,1,2......^16^+1 ΜδΜ+2 _ ^16η+3 ~^\βη+5 ^16η+6 ~ ^\βη+ ~^\βη+9 Μδ^+ΙΟ -
-^\6η+\\ ~ ΜδΜ+12 +Dl6M+13 ^\6n+\4 ~ ^\6n+\5 ~ ^\6n+\6 ~ ^^An/7—0,1,2......^16^+1 ^\6n+2 +Dl6M+3~ ^\βη+5 ~ ^16η+6 ~ ^\βη+ ~^\βη+9 Μδ^+ΙΟ +
-^\6n+\\ ΜδΜ+12 _ ^16m+13 _ ^\6n+\4 ~ ^\6n+\5 ~ ^\6n+\6 ~ ^^An/7—0,1,2......^16^+1 ^\6n+2 +Dl6M+3^\βη+5 ^16η+6 ^\βη+ ~ ^\βη+9 ~ ^I6n+10 -
-^\6n+\\ ~ -^\6n+\2 ~ ^16m+13 _ ^\6n+\4 ~ ^\6n+\5 ~ ^\6n+\6 ~ ^^An/7—0,1,2......即分别得到8倍的波长λ 1、λ 2、λ 3和λ 4的光电容积脉搏波、Dfn、Ζ)丄3和 daa，而且完全消除了背景信号Af的影响。得到各个波长的光电容积脉搏波，据此计算出λ 1、λ 2、λ 3和λ 4的光电容积脉搏波白勺存值和峰值Ιη ηλ1、Ι·λ1、IminA2、ImaxA2、Ιπ ηλ3、ImaxA3、Ιπ ηλ4 和 ImaxA4°再由IminM、ImaxM、IminA2、ΙΜχλ2、Ι-λ3、Ι_λ3、IminA4 和 Ι_λ4 计算各个波长所对应的吸光度差值ΔΑ，可以得到由吸光度差值ΔΑλ1、ΔΑλ2……成的光谱值。综上所述，本发明实施例提供了一种方波调制光电容积脉搏波测量装置和测量方法，本发明实施例依据朗伯-比尔定律，采用方波频分调制和数字解调技术设计一种方波调制多波长发光二极管的光电容积脉搏波及其谷值和峰值，通过谷值和峰值得到光谱值的装置和测量方法，具有测量精确、电路简单、无需调试、工艺性好以及成本低廉的特点。本领域技术人员可以理解附图只是一个优选实施例的示意图，上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
权利要求
1.一种方波调制光电容积脉搏波测量装置，其特征在于，所述光电容积脉搏波测量装置包括微处理器、至少2种发光二极管、光敏器件、电流/电压转换放大器和模数转换器，所述微处理器输出不同频率且成2倍比率关系的方波，所述方波驱动所述至少2种发光二极管，所述发光二极管发出的光经被测手指后被所述光敏器件接收，所述光敏器件转换成电压信号，所述电压信号经所述电流/电压转换放大器转换成预设幅值电压信号，所述模数转换器将所述预设幅值电压信号转换成数字信号，所述微处理器对所述数字信号进行处理，获取光电容积脉搏波及其谷值和峰值，通过所述谷值和所述峰值得到光谱值。
2.根据权利要求1所述的一种方波调制光电容积脉搏波测量装置，其特征在于，所述微处理器采用MCU、ARM、DSP或FPGA中的任意一种。
3.一种用于权利要求1所述的一种方波调制光电容积脉搏波测量装置的测量方法，其特征在于，所述方法包括以下步骤(1)微处理器采用不同频率且成2倍比率关系的方波驱动至少2种发光二极管；(2)所述发光二极管发出的光经过被测手指后由光敏器件接收转换成电压信号，所述电压信号经过电流/电压转换放大器放大成预设幅值电压信号；(3)所述预设幅值电压信号经模数转换器转换成数字信号送入所述微处理器；(4)所述微处理器对所述数字信号进行分离处理得到光电容积脉搏波并消除背景光的干扰；(5)根据所述光电容积脉搏波获取谷值和峰值；(6)对所述谷值和所述峰值进行计算得到吸光度差值，通过所述吸光度差值获取光谱值。
全文摘要
本发明公开了一种方波调制光电容积脉搏波测量装置和测量方法，微处理器输出不同频率且成2倍比率关系的方波，方波驱动至少2种发光二极管，发光二极管发出的光经被测手指后被光敏器件接收，光敏器件转换成电压信号，电压信号经电流/电压转换放大器转换成预设幅值电压信号，模数转换器将预设幅值电压信号转换成数字信号，微处理器对数字信号进行处理，获取光电容积脉搏波及其谷值和峰值，通过谷值和峰值得到光谱值；测量方法包括微处理器对数字信号进行分离处理得到光电容积脉搏波并消除背景光的干扰；根据光电容积脉搏波获取谷值和峰值；对谷值和峰值进行计算得到吸光度差值，通过吸光度差值获取光谱值。具有测量精确、电路简单以及成本低廉。
文档编号A61B5/1455GK102389313SQ20111023639
公开日2012年3月28日 申请日期2011年8月17日 优先权日2011年8月17日
发明者刘近贞, 周梅, 李刚, 林凌, 郝丽玲 申请人:天津大学