申请/专利权人：国立大学法人信州大学

申请日：2016-02-19

公开（公告）日：2022-12-02

公开（公告）号：CN107405114B

主分类号：A61B5/1455

分类号：A61B5/1455;A61B5/02

优先权：["20150313 JP 2015-051252"]

专利状态码：有效-授权

法律状态：2022.12.02#授权;2017.12.22#实质审查的生效;2017.11.28#公开

摘要：本发明所涉及的非侵入式血糖值测定装置1具备：脉搏波测定部2，其具备对被实验者的加速度脉搏波进行测定的FBG传感器4；以及数据处理部3，其根据进行测定所得的加速度脉搏波的波形信息，且基于预先规定的相关关系而计算出对该加速度脉搏波进行测定时的所述被实验者的血糖值。相关关系是将利用侵入式测定法进行测定所得的血糖值作为目标变量且将同时测定所得的加速度脉搏波作为说明变量来进行PLS回归分析而构建的定标线。能够实现如下非侵入式血糖值测定装置，该非侵入式血糖值测定装置能够以与侵入式血糖值测定法同等程度的测定精度而对血糖值进行测定。

主权项：1.一种非侵入式血糖值测定方法，其特征在于，所述非侵入式血糖值测定方法包括血糖值计算步骤，在该血糖值计算步骤中，根据对被实验者进行测定所得的加速度脉搏波的波形信息，且基于预先规定的加速度脉搏波与血糖值的相关关系而求出对所述被实验者的所述加速度脉搏波进行测定时的血糖值，所述相关关系是针对所述被实验者或者不同的被实验者通过侵入式测定法进行测定所得的作为血糖值的第一血糖值、与在测定该第一血糖值的同时进行测定所得的作为加速度脉搏波的第一加速度脉搏波之间的相关关系，所述相关关系是将所述第一血糖值作为目标变量且将所述第一加速度脉搏波作为说明变量来进行回归分析而构建的定标线。

全文数据：非侵入式血糖值测定方法以及非侵入式血糖值测定装置技术领域[0001]本发明涉及能够利用对被实验者进行测定所得的加速度脉搏波而测定该被实验者的血糖值的非侵入式血糖值测定方法以及装置。背景技术[0002]在糖尿病的预防、治疗这方面，对血糖值的测定是必不可少的。目前为止，一般通过穿刺而采集血液、并利用酶电极法等的血糖值传感器进行测定，由此进行血糖值的自我测定SMBG:SelfMonitoringofBloodGlucose。这种侵入式血糖值测定方法会对被采集血液的被实验者造成心理负担。另外，无法忽视穿刺、对电极等物资的一次性使用所造成的经济方面的负担，这也成为医疗费用增加的主要原因之一。[0003]目前为止，正在开展针对用于解决侵入式血糖值测定法的各种问题的非侵入式血糖值测定法的研发。作为其研发成果之一，研究了基于光谱分析的非侵入式测定法。本发明的发明人等在专利文献1中提出了通过光谱分析而对血糖值进行测定的生物体信息测定装置。另外，除了光学方法以外，还提出了考虑血液的流体力学的特性而对血糖值进行测定的系统，但该系统并未进入实用阶段。[0004]专利文献[0005]专利文献1:日本特开2012-191969号公报发明内容[0006]以往提出的光学方法等非侵入式血糖值测定方法有时在测定精度方面无法稳定地满足实用上的要求。另外，多数情况下在装置成本方面无法满足实用上的要求。[0007]本发明的课题在于，提供不利用光谱分析便能够对血糖值进行测定的非侵入式血糖值测定方法以及装置。[0008]另外，本发明的课题在于，提供能够以与侵入式血糖值测定法同等程度的测定精度而对血糖值进行测定的、廉价结构的非侵入式血糖值测定方法以及装置。[0009]本发明的发明人着眼于能够利用光纤布拉格光栅传感器（以下，称为FBG传感器。）等传感器进行观测的人的加速度脉搏波中含有依赖于动脉血的密度等而变化的成分这一点，获得了如下见解:根据测定所得的加速度脉搏波的波形图案的特征而能够提取出该测定时的血糖值信息。[0010]本发明的非侵入式血糖值测定方法以及装置是基于上述见解而完成的，其特征在于，对被实验者的加速度脉搏波进行测定，根据测定所得的加速度脉搏波的波形信息，并基于通过侵入式测定法进行测定所得的血糖值与同时测定的加速度脉搏波之间的相关关系而提取出被实验者的血糖值信息。[0011]S卩，本发明的非侵入式血糖值测定方法的特征在于，包括血糖值计算步骤，在该血糖值计算步骤中，根据对被实验者进行测定所得的加速度脉搏波的波形信息，且基于预先规定的加速度脉搏波与血糖值的相关关系而求出对被实验者的加速度脉搏波进行测定时的血糖值，所述相关关系是针对所述被实验者或者不同的被实验者通过侵入式测定法进行测定所得的作为血糖值的第一血糖值、与在测定该第一血糖值的同时进行测定所得的作为加速度脉搏波的第一加速度脉搏波之间的相关关系。[0012]本发明的发明人等将基于侵入式测定法观察血液式的测定法得到的血糖值基于遵照国际标准IS015197的侵入式血糖自我测定器得到的测定值作为目标变量、且将同时测定的加速度脉搏波作为说明变量来进行了PLS回归分析，结果能够确认：由此能够制作规定的定标精度的定标线。另外，能够确认:即使在该定标线的验证中也能够获得规定的血糖值预测精度。根据本发明的测定方法，能够实现如下廉价结构的测定装置，该测定装置包括:加速度脉搏波的测定部；以及数据处理部，其利用预先设定的定标线而进行血糖值信息的提取。[0013]在本发明中，通过使用FBG传感器，能够直接对加速度脉搏波进行测定。另外，血糖值计算步骤可以基于定标线来进行，该定标线是将利用侵入式测定法进行测定所得的血糖值作为目标变量、且将同时测定所得的加速度脉搏波作为说明变量来预先进行回归分析而构建的，特别是预先进行PLS回归分析而构建的。[0014]在该情况下，为了提高血糖值测定精度，作为用于构建定标线的说明变量，优选采用所述第一加速度脉搏波的标准化脉搏波。作为标准化脉搏波，可以采用针对第一加速度脉搏波进行其波形位移的标准化、以及其波形长度的标准化而获得的单脉冲的波形数据。同样，对于为了测定血糖值而对被实验者进行测定所得的加速度脉搏波，也同样优选采用对其实施标准化后所得的数据。[0015]本发明的发明人等将基于侵入式测定法的血糖值作为目标变量、且将同时测定所得的加速度脉搏波作为说明变量而进行了PLS回归分析，结果能够确认：由此能够在80mgd1〜178mgd1的血糖值范围内制作定标精度为±15mgd1的定标线。能够确认：即使在该定标线的验证中，也能够获得同等的血糖值预测精度。[0016]另外，还能够确认:利用FBG传感器而能够以20kHz的高时间分辨率、亚微米以下的高灵敏度而对加速度脉搏波进行检测，从而能够提高利用基于光学方法的以往的非侵入式测定法而得到的定标线精度。进而发现:检测出的定标线的因子factors在加速度脉搏波的传播速度发生变化的位置处呈现出较大的作用，从而是适当的。[0017]根据本发明，能够实现能通过使用加速度脉搏波的新方法而对血糖值进行测定的非侵入式血糖值测定方法以及装置。另外，根据本发明，能够利用加速度脉搏波实现能以与侵入式血糖值测定法同等程度的测定精度对血糖值进行测定的非侵入式血糖值测定方法以及装置。并且，根据本发明，能够利用加速度脉搏波实现与以往的侵入式血糖值测定方法、光学方法等非侵入式血糖值测定方法相比而更廉价的结构的非侵入式血糖值测定方法以及装置。附图说明[0018]图1是本发明的实施方式所涉及的非侵入式血糖值测定装置的概要结构图。[0019]图2.1是示出实验例1的基于标准化方法1的标准化脉搏波的曲线图。[0020]图2.2是示出实验例1的标准化方法1的情况下的血糖值预测结果的曲线图、以及示出定标线的构建结果以及验证结果的图表。[0021]图2.3是示出实验例1的标准化方法1的情况下的负荷loading结果的曲线图。[0022]图2.4是示出实验例1的基于标准化方法2的标准化脉搏波的曲线图。[0023]图2•5是示出实验例1的标准化方法2的情况下的血糖值的预测结果的曲线图、以及定标线的构建结果和验证结果的图表。[0024]图2.6是示出实验例1的标准化方法2的情况下的负荷结果的曲线图。[0025]图3.1是示出针对实验例2通过标准化方法1实施标准化并截取的波形、以及以最小采样点数而使得长度一致的波形的曲线图。[0026]图3•2是分别示出实验例2的各波形的情况下的血糖值的预测结果的曲线图。[0027]图3•3是分别示出实验例2的各波形的情况下的定标线的构建结果以及验证结果的图表。[0028]图3.4是分别示出实验例2的各波形的情况下的负荷结果的曲线图。[0029]图3•5是示出针对实验例2通过标准化方法2实施标准化、且以5000点以及10000点的采样点数而截取的波形的曲线图。[0030]图3.6是分别示出实验例2的各波形的情况下的血糖值的预测结果的曲线图。[0031]图3.7是分别示出实验例2的各波形的情况下的定标线的构建结果以及验证结果的图表。[0032]图3.8是分别示出实验例2的各波形的情况下的负荷结果的曲线图。[0033]图4.1是示出实验例3的被实验者的标准化脉搏波的曲线图。[0034]图4.2是示出实验例3的情况下的血糖值的预测结果的曲线图。[0035]图4.3是示出实验例3的定标线的构建以及验证结果的图表。[0036]图4.4是示出实验例3的负荷结果的曲线图。[0037]图5.1是分别示出实验例4的以不同的采样点数而截取的3种波形的曲线图。t〇〇38]图5.2是分别示出实验例4的以不同的采样点数而截取的3种波形的各情况下的血糖值的预测结果的曲线图。[0039]图5.3是分别示出实验例4的以不同的采样点数而截取的3种波形的各情况下的定标线的构建结果以及验证结果的图表。[0040]图5.4是分别示出实验例4的以不同的采样点数而截取的3种波形的各情况下的负荷结果的曲线图。[°041]图6.1是示出实验例5的通过标准化方法1、2而实施标准化后的肘部的波形的曲线图。[0042]图6.2是分别示出实验例5的通过标准化方法1、2而实施标准化后的肘部的波形的各情况下的血糖值的预测结果的曲线图。[0043]图6.3是分别示出实验例5的通过标准化方法1、2而实施标准化后的肘部的波形的各情况下的定标线的构建结果以及验证结果的图表。[0044]图6•4是分别示出实验例5的通过标准化方法1、2而实施标准化后的肘部的波形的各情况下的负荷结果的曲线图。[0045]图7•1示出了实验例6的手腕以及肘部的脉搏波的原始波形。[0046]图7.2是分别示出使用实验例6的原始波形的各情况下的血糖值的预测结果的曲线图。[0047]图7.3是分别示出使用实验例6的原始波形的各情况下的定标线的构建结果以及验证结果的图表。[0048]图7.4是分别示出使用实验例6的原始波形的各情况下的负荷结果的曲线图。具体实施方式[0049]以下，参照附图，对应用本发明的、利用加速度脉搏波的非侵入式血糖值测定装置的实施方式进行说明。[0050][测定装置的概况][0051]图1是本实施方式所涉及的非侵入式血糖值测定装置的概要结构图。非侵入式血糖值测定装置1以下，称为“血糖值测定装置1”。）具有:脉搏波测定部2，其对被实验者的加速度脉搏波进行测定;数据处理部3，其根据测定所得的加速度脉搏波的波形图案的特征而提取出被实验者的血糖值信息；以及操作显示部20。此外，在以下说明中，只要未特别声明，则脉搏波意味着加速度脉搏波。[0052]脉搏波测定部2具备:FBG传感器4;光源5，其射出向FBG传感器4入射的基准光；以及光检测器6，其对来自FBG传感器4的反射光进行检测，数据处理部3基于光检测器6的检测结果而对被实验者的血糖值进行计算。FBG传感器4安装于对被实验者的脉搏波进行测定的部位、例如手腕的内侧位置、肘部的内侧位置等而使用。[0053]在本例中，FBG传感器4具备FBG传感器1以及FBG传感器2,经由环行器7而将来自上述传感器的反射光向马赫曾德千涉仪8引导。利用光检测器6对来自马赫曾德千涉仪8的输出光进行检测。马赫曾德干涉仪8利用入射侧的分束器9而使得反射光向具有光路差的2条光路分离，并利用出射侧的分束器1〇而使得分离的2条光路再次重叠为一条光路，由此形成干涉光。由于相干光与光路差相应地产生干涉条纹，因此，通过对干涉条纹的图案进行测定，能够计算出FBG传感器4所产生的布拉格波长的变化，从而能够检测出形变的变化、即脉搏波。[0054]以下，示出脉搏波测定部2的各部分的各要素。[0055]光源5:ASEAmplifiedSpontaneousEmission:放大自发射光[0056]FGB传感器4:[0057]布拉格波长[0058]FBG传感器1:1550±0.5nm[0059]FBG传感器2:1560±0.5nm[0060]FBG传感器1、2的长度:5mm[0061]光纤直径：145wn[0062]光纤的芯径：10.5wn[0063]光纤材料:桂玻璃[0064]光检测器6:InGaAsPINPD[0065]波长分辨率：±0•1pm[0066]数据处理部3具备数据解析部11。将通过侵入式测定法观察血液式的测定法进行测定所得的血糖值作为目标变量、且将同时测定的加速度脉搏波作为说明变量而进行PLS回归分析，将由此构建的定标线预先存储保持于数据解析部11。数据解析部11基于定标线、并根据利用脉搏波测定部2进行测定所得的加速度脉搏波而对被实验者的血糖值进行预测推断）。数据处理部3构成为以微机为中心，并通过执行存储保持的解析用程序而作为数据解析部11来发挥功能。_[0067]血糖值测定装置1预先将定标线存储保持于数据解析部11。如上所述，定标线是通过将利用侵入式测定法测定所得的血糖值作为目标变量且将同时测定的加速度脉搏波作为说明变量、并进行PLS回归分析而构建的。作为定标线，使用事先根据血糖值测定对象的被实验者而获得的定标线。也可以取而代之地使用事先根据其他被实验者而获得的定标线。[0068]当测定血糖值时，例如，使血糖值测定对象的被实验者在安静状态下以使得作为测定部位的手腕处于与心脏相同的高度的方式形成仰卧的姿势。在该状态下，利用医疗用胶带等将FBG传感器4固定于被实验者的一个手腕的挠骨动脉上，以规定的采样周期在规定的期间内对加速度脉搏波进行测定。[0069]将测定所得的加速度脉搏波数据读取至数据处理部3而实施规定的数据处理。数据处理与构建定标线时针对加速度脉搏波的处理相同。[0070]例如，在数据处理部3中，首先，为了除去噪声，使加速度脉搏波数据从规定的通频带、例如0.5Hz〜fflz的通频带的带通滤波器未图示通过。接下来，为了将除去噪声后的加速度脉搏波数据分割为单脉冲的脉搏波部分，以加速度脉搏波数据中所含有的各峰值为基准点而截取1个脉冲。对截取的多个单脉冲的加速度脉搏波数据进行平均而计算出单脉冲的平均脉搏波数据。针对该平均脉搏波数据，实施其振幅波形位移的标准化以及长度采样点数的标准化。[0071]在数据处理部3的数据解析部11中，根据实施了标准化之后的单脉冲的平均脉搏波数据并利用存储保持的定标线而对测定被实验者的加速度脉搏波时的血糖值进行计算。将计算出的血糖值数据例如发送至操作显示部20、并在其显示画面上对该血糖值数据进行显不〇[0072]为了确认基于血糖值测定装置1的血糖值测定方法的有效性，本发明的发明人等进行了各种实验。以下列举其中的一部分实验例。[0073][实验例1:手腕的加速度脉搏波的解析之1][0074]实验方法以及解析方法）[0075]利用医疗用胶带将FBG传感器4粘贴于被实验者的右手腕的挠骨动脉上而对被实验者的脉搏波进行测定。在测定脉搏波的同时，利用血糖仪（产品名称为“FreestylePrecisionPro”，ABB0TTJAPAN株式会社制对被实验者的血糖值进行测定，并将测定值设为参照血糖值。被实验者是1名20多岁的男性，脉搏波的测定条件如下。[0076]采样频率:20kHz[0077]测定时间:从自动血压计开始测定起直至测定结束为止[0078]测定次数:80次[0079]测定时的被实验者的状态、姿势:作为安静状态下的仰卧的姿势，将测定部位保持于与心脏相同的高度位置[00S0]在测定时，为了除去噪声而使用通频带为0_5Hz〜5Hz的带通滤波器。另外，为了将由身体的移动而引起的噪声除去，以下述方式而使用平均脉搏波。[0081]即，针对在各测定中所获得的除去噪声之后的脉搏波，以其中出现的各峰值为基准点而针对每1个脉冲截取出多个单脉冲脉搏波，对这些脉搏波实施平均而生成单脉冲的平均脉搏波。在本实验中，测定次数为8〇次，因此生成80个单脉冲的平均脉搏波。[0082]并且，在各测定中所获得的平均脉搏波的峰值高度和长度的差异成为误差的主要原因，因此，为了将其除去，以下述方式对平均脉搏波实施标准化处理并使用实施了标准化后的平均脉搏波。[0083]接下来，为了对测定所得的脉搏波的波形标准化后的平均脉搏波的波形）、与同时测定的血糖值的相关关系进行解析，将各参照血糖值作为目标变量、并将与各参照血糖值同时测定的脉搏波作为说明变量，通过PLS回归分析而构建了回归模型（定标线）。另外，利用在回归模型的构建中未使用的数据集而对回归模型进行了验证。在回归模型的构建中使用的数据的个数为6〇，将剩余的20个数据用于对回归模型的验证。[0084]加速度脉搏波的标准化方法）[0085]作为脉搏波平均脉搏波）的标准化方法，采用了仅对波形位移实施标准化标准化方法1、以及对波形位移和采样点数波长这二者实施标准化标准化方法2的2种方法。[0086]标准化方法1:将脉搏波的峰值设为1、且将最小值设为〇。此外，以通过80次的测定而获得的8〇个单脉冲的平均脉搏波中的最小采样点数而使得这些脉搏波的长度一致。[0087]标准化方法2:除了上述的波形位移的标准化以外，还以20000点的采样点数而使得脉搏波的长度一致。[0088]实验结果卜丄）[0089]图2.1是示出通过标准化方法1对脉搏波实施标准化而获得的标准化脉搏波标准化后的平均脉搏波）的标准化波形位移的曲线图。图2.2a是示出血糖值的预测结果的曲线图，图2.2⑹是示出定标线的构建结果以及验证结果的说明图。另外，图2.3是示出负荷结果的曲线图。[0090]此处，在图2.2a所示的曲线图中，在基于误差栅格分析法EGA法）的误差栅格上标绘出由圆圈表示的定标线的构建数据和由方块表示的验证数据。误差栅格的横轴表示参照血糖值，纵轴表示预测血糖值。[0091]对角线表示参照血糖值和预想血糖值的一致性，当标绘点处于对角线的上方时，表示预测血糖值被评价为过大，当标绘点处于对角线的下方时，表示预测血糖值被评价为过小。在误差栅格中，A区是预测血糖值仅偏差20%的区域。另外，参照血糖值低于70mgdl的A区的部分是表示低血糖值（70mgdl的区域。B区表示虽然预测血糖值相对于参照血糖值向上方和下方偏差的幅度大于20%但却进行了良性的治疗的区域。C区表示对优选的血糖值水平进行了过度校正的区域。D区表示检测出错误的、会引发“危险性的失败”的区域，E区表示进行了“错误的治疗”的区域。[0092]实验结果1-2[0093]图2.4是示出通过标准化方法2对脉搏波实施标准化而获得的标准化脉搏波的标准化波形位移的曲线图，图2.5a是示出血糖值的预测结果的曲线图，图2.5b是示出定标线的构建结果以及验证结果的说明图。图2.6是示出负荷结果的曲线图。[0094]根据实验结果1-1以及1-2,标准化方法2与标准化方法1相比，能够以更高的精度预测血糖值。根据图2.3、图2•6，在标准化波形位移的峰值的下降、第二个波的上升时，负荷的值变大。可以认为：由于血糖值的变化是对血液粘度造成影响的主要原因之一，因此，血糖值的变化会对容积脉搏波的速度造成影响。可以认为：由于利用FBG传感器4测定的脉搏波是加速度脉搏波，因此，负荷值在该加速度脉搏波的倾斜部分变大。[0095][实验例2:手腕的加速度脉搏波的解析之2][0096]实验方法以及解析方法）[0097]利用医疗用胶带将FBG传感器4粘贴于与实验例1的情况相同的被实验者的右手腕的挠骨动脉上而对脉搏波进行了测定。另外，在测定脉搏波的同时，利用血糖仪产品名称为“FreestylePrecisionPro”，ABB0TTJAPAN株式会社制对血糖值进行测定，并将测定值作为参照血糖值。脉搏波的测定条件如下。[0098]采样频率:10kHz[0099]测定时间：从自动血压计开始测定起直至测定结束为止[0100]测定次数:80次[0101]测定时的被实验者的状态、姿势:取安静状态下的仰卧的姿势，将测定部位保持于与心脏相同的高度位置[0102]在测定时，为了除去噪声而使用通频带为0.5Hz〜5Hz的带通滤波器。另外，从在各测定中所获得的脉搏波中截取出单脉冲的脉搏波并获得多个单脉冲的脉搏波，对这些单脉冲的脉搏波实施平均而生成平均脉搏波。[0103]为了对测定所得的脉搏波平均脉搏波的波形与血糖值的相关关系进行解析，将各参照血糖值作为目标变量、且将与各参照血糖值同时测定所得的脉搏波作为说明变量，并通过PLS回归分析而构建了回归模型（定标线）。利用在回归模型的构建中未使用的数据集而对回归模型进行了验证。用于回归模型的构建的数据数为60,将剩余的20个数据用于对回归模型的验证。[0104]加速度脉搏波的标准化方法以及截取方法）[0105]脉搏波的标准化方法采用标准化方法1波形位移的标准化）、以及标准化方法2波形位移和采样点数的标准化这2种方法。[0106]标准化方法1:将脉搏波的峰值设为1，并将最小值设为0。此外，以最小采样点数而使得脉搏波的长度一致。[0107]标准化方法2:除了波形位移的标准化以外，以10000点的采样点数而使得脉搏波的长度一致。[0108]在脉搏波的截取方法中，在以上述方式对该波形实施了标准化之后，以5000点的采样点数而截取了标准化波形。[0109]在解析中，使用如下共计4种脉搏波：以5000点的采样点数而对通过标准化方法1实施标准化所得的标准化波形进行截取后的脉搏波、和以最小采样点数而使得长度一致的脉搏波；以及以5000点的采样点数而对通过标准化方法2实施标准化所得的标准化波形进行截取后的脉搏波、和以10000点而使得长度一致的脉搏波。[0110]实验结果2-1[0111]图3.1a、图3.1b是示出通过标准化方法1而实施标准化并以5000点截取的标准化脉搏波的标准化波形位移、以及以最小采样点数使得长度一致后的标准化波形位移的曲线图。图3•2a、图3.2b是分别示出各情况下的血糖值的预测结果的曲线图，图J3a、图3.3b是分别示出各情况下的定标线的构建结果以及验证结果的说明图。另外，图3.4a、图3.4⑹是分别示出各情况下的负荷结果的曲线图。[0112]实验结果2-2[0113]图3.5a、图3.5b是示出通过标准化方法2而实施标准化并以5000点而截取的标准化脉搏波、以及以10000的点数而使得长度一致的情况下的标准化脉搏波的曲线图。图3.6a、图3.6⑹是分别示出各情况下的血糖值的预测结果的曲线图，图3.7a、图3_7⑹是分别示出各情况下的定标线的构建结果以及验证结果的说明图。另外，图3.8a、图3.8⑹是分别示出各情况下的负荷结果的曲线图。__[0114]获得了与实验例1的情况相同的实验结果。即，根据实验结果2-1以及2-2,标准化方法2与标准化方法1相比，能够以更高的精度预测血糖值。[0115][实验例3:手腕的加速度脉搏波的解析之3][0116]实验方法、解析方法）一[0117]针对与实验例1、2的情况下的被实验者不同的20多岁的男性被实验者而进行了实验。[0118]利用医疗用胶带将FBG传感器固定于被实验者的右手腕的挠骨动脉上而对脉搏波进行了测定。另外，在测定脉搏波的同时，利用血糖仪产品名称为“FreestylePrecisionPro”，ABB0TTJAPAN株式会社制）对血糖值进行测定，并将测定值作为参照血糖值。脉搏波的测定条件如下。[0119]采样频率:10kHz[0120]测定时间：自利用血糖仪开始测定起的15秒钟[0121]测定次数:60次[0122]测定时的被实验者的状态、姿势：在安静状态以使得测定部处于与心脏相同的高度的方式形成为仰卧的姿势[0123]为了将噪声除去，使测定所得的脉搏波从通频带为〇.5Hz〜5Hz的带通滤波器通过。另外，以除去噪声后的脉搏波的峰值为基准点而截取出1个脉冲。为了降低由身体移动而引起的噪声，对截取的单脉冲的脉搏波实施了平均。由于一般的成人的脉搏数的平均值为60bpm，因此，在1次的测定中，截取出约15个单脉冲的脉搏波。对上述脉搏波实施平均而计算出60bpm时的平均脉搏波。[0124]另外，由于脉搏波的峰值高度和长度的差异成为误差的主要原因，因此，与实验例1、2的情况下的标准化方法2相同，针对平均脉搏波而实施了将平均脉搏波的最初的点（峰值设为1且将最小值设为〇的波形位移的标准化、以及将采样数设为10000点的长度的标准化。[0125]将通过标准化而获得的标准化脉搏波作为说明变量、且将参照血糖值作为目标变量，并通过PLS回归分析而构建定标线，由此计算出血糖值。[0126]实验结果）[0127]图4.1中示出了被实验者的标准化脉搏波，图4.2中示出了血糖值计算结果，图4.3的表中示出了定标线的构建以及验证结果，图4.4中示出了PLS回归分析的负荷结果。根据图4.2以及图4.3,相关度较高、且SEP、SEC均较小，验证结果全部收敛于A区中。[0128][实验例4:加速度脉搏波的截取位置的研宄][0129]根据实验例1、2、3的负荷结果（图2.3、图2.6、图3.4、图3.8、图4.4，对于目前为止的任意被实验者，负荷值在采样点数至6000点为止的标准化脉搏波的前半部分都变大。据此可以认为:标准化脉搏波的前半部分对于血糖仪的测量较为有效。[0130]为了根据负荷而确认对于血糖值的预测有效的部分，改变脉搏波的截取位置而进行了解析。将使用的脉搏波设为手腕的通过标准化方法2而实施了标准化后的标准化脉搏波的波形，并将采样频率设为20kHz。在采样点数为10000点、lf5〇〇〇点以及19000点的3个不同的位置进行了截取。除此以外的测定条件与实验例1的情况相同。[0131]图5.1a、⑹、（c中分别示出了截取后的标准化脉搏波，图5.2a、⑹、（c中分别示出了各情况下的血糖值的预测结果，图5.3a、⑹、（c中分别示出了各情况下的定标线的构建结果以及验证结果，图5.4a、〇3、（c中分别示出了各情况下的负荷结果。[0132]无论在任何情况下，都能在标准化脉搏波与血糖值之间确认到较高的相关关系，从而能够确认:能够高精度地预测血糖值。可以认为:通过以包含作为标准化脉搏波的前半部分的采样点数达到6000点左右的部分的方式进行标准化脉搏波的截取，能够高精度地预测血糖值。[0133][实验例5:肘部的加速度脉搏波的解析][0134]在实验例1、2、3中，利用手腕的脉搏波而进行了解析。为了确认利用在手腕以外的位置处进行测定所得的脉搏波对血糖值的测定的有效性，在本实验中，针对与实验1、2的情况相同的被实验者而对肘部的脉搏波进行了测定，并利用测定所得的脉搏波进行了同样的解析。[0135]与手腕的情况相同，标准化方法采用波形位移的标准化标准化方法1、以及波形位移和采样点数的标准化标准化方法2这2种方法。另外，将测定的采样频率设为10kHz。除此以外的测定条件与实验例1的情况相同。[0136]图6.1a、图6.1⑹中分别示出了通过标准化方法1、2实施标准化之后的各波形，图6.2a、图6.2⑹中分别示出了各情况下的血糖值的预测结果，图6.3a、图6.3⑹中分别示出了各情况下的定标线的构建结果以及验证结果，图6.4a、图6.4b中分别示出了各情况下的负荷结果。[0137]针对在肘部测定所得的脉搏波的标准化脉搏波，也能确认：其与血糖值之间存在规定的相关关系。另外，标准化方法2与标准化方法1相比，能够以更高的精度预测血糖值。[0138][实验例6:手腕和肘部的原始波形的解析][0139]为了确认测定所得的脉搏波的标准化对于提高血糖值测定的精度较为有效，在本实验中，利用测定所得的加速度脉搏波的原始波形对血糖值进行了测定。[0140]g卩，在本实验中，与实验例1、2同样地针对手腕和肘部这二者而测定脉搏波。利用针对测定所得的手腕和肘部的脉搏波并未进行标准化的原始波形而进行了解析。将采样频率设为10kHz。除此以外的测定条件与实验例1的情况相同。[0141]图7.1a、图7.1⑹中示出了手腕以及肘部的脉搏波的原始波形。图7.2a、图7.2b是分别示出各情况下的血糖值的预测结果的曲线图。图7.3a、图7.3b是分别示出各情况下的定标线的构建结果以及验证结果的图表。图7•4a、图7•4b是分别示出各情况下的负荷结果的曲线图。能够确认:与使用根据原始波形而获得的定标线的情况相比，使用测定所得的脉搏波的标准化脉搏波较为有效。[0142][利用加速度脉搏波的血糖值测定方法的通用性][0143]上述实验例是针对同一被实验者而构建定标线并对其进行验证。也可以针对某被实验者而构建血糖值计算用的定标线，并利用构建的定标线且根据其他被实验者的脉搏波而对血糖值进行测定。即使在该情况下，也能够以规定的精度而对血糖值进行测定。[0144][其他实施方式][0145]脉搏波传感器）[0146]上述血糖值测定装置1利用FBG传感器4对脉搏波进行测定。也可以利用FBG传感器以外的脉搏波传感器而直接对加速度脉搏波进行测定。另外，还可以对容积脉搏波进行测定并对测定脉搏波进行二次微分而求出加速度脉搏波，进而基于该加速度脉搏波对血糖值进行测定。[0147]作为脉搏波传感器，以往已知光电式结构、机械式结构、阻抗式结构、应变仪式结构这4种类型的脉搏波传感器。虽然能够使用上述脉搏波传感器，但优选使用能够更高精度地检测出脉搏波的传感器。例如，可以使用高灵敏度感压传感器L-系列传感器，芬兰EMFit公司制）、触觉传感器T40006000系列，美国PressureProfileSystems公司制）等。[0148]加速度脉搏波和血糖值的相关关系的计算方法）[0149]在上述实施方式中，作为回归分析法而采用了PLS回归分析法。也可以利用除此以外的回归分析法而求出通过侵入式测定法进行测定所得的血糖值、与同时测定的加速度脉搏波之间的相关关系。[0150]脉搏波测定部位）[0151]在上述实验例中，在手腕以及肘部的部位对脉搏波进行了测定。作为脉搏波的测定部位，也可以是除此以外的部位。

权利要求：1.一种非侵入式血糖值测定方法，其特征在于，所述非侵入式血糖值测定方法包括血糖值计算步骤，在该血糖值计算步骤中，根据对被实验者进行测定所得的加速度脉搏波的波形信息，且基于预先规定的加速度脉搏波与血糖值的相关关系而求出对所述被实验者的所述加速度脉搏波进行测定时的血糖值，所述相关关系是针对所述被实验者或者不同的被实验者通过侵入式测定法进行测定所得的作为血糖值的第一血糖值、与在测定该第一血糖值的同时进行测定所得的作为加速度脉搏波的第一加速度脉搏波之间的相关关系。2.根据权利要求1所述的非侵入式血糖值测定方法，其特征在于，所述相关关系是将所述第一血糖值作为目标变量且将所述第一加速度脉搏波作为说明变量来进行PLS回归分析而构建的定标线。3.根据权利要求2所述的非侵入式血糖值测定方法，其特征在于，_作为用于构建所述定标线的所述说明变量，采用对所述第一加速度脉搏波实施标准化所得的第一标准化脉搏波，所述第一标准化脉搏波是对所述第一加速度脉搏波进行其波形位移的标准化、以及其波形长度的标准化而获得的单脉冲的波形数据。4.根据权利要求1所述的非侵入式血糖值测定方法，其特征在于，所述非侵入式血糖值测定方法包括测定步骤，在该测定步骤中，利用光纤布拉格光栅传感器对所述被实验者测定所述加速度脉搏波，在所述血糖值计算步骤中，利用将针对所述被实验者进行测定而获得的所述加速度脉搏波实施标准化后的标准化脉搏波来求出所述血糖值，所述标准化脉搏波是对测定所得的所述加速度脉搏波进行其波形位移的标准化、以及其波形长度的标准化而获得的单脉冲的波形数据。5.—种非侵入式血糖值测定装置，其特征在于，所述非侵入式血糖值测定装置具有：脉搏波测定部，其对被实验者的加速度脉搏波进行测定；存储部，其对加速度脉搏波与血糖值之间的规定的相关关系进行存储保持;以及数据处理部，其根据进行测定所得的所述加速度脉搏波的波形信息，且利用所述相关关系来求出所述被实验者的血糖值，所述相关关系是针对所述被实验者或者不同的被实验者通过侵入式测定法进行测定所得的作为血糖值的第一血糖值、与在测定该第一血糖值的同时进行测定所得的作为加速度脉搏波的第一加速度脉搏波之间的相关关系。6.根据权利要求5所述的非侵入式血糖值测定装置，其特征在于，所述脉搏波测定部具备光纤布拉格光栅传感器。7.根据权利要求5所述的非侵入式血糖值测定装置，其特征在于，所述相关关系是将所述第一血糖值作为目标变量且将所述第一加速度脉搏波作为说明变量来进行PLS回归分析而构建的定标线，所述数据处理部具备数据解析部，该数据解析部利用所述定标线而对所述血糖值进行计算。8.根据权利要求7所述的非侵入式血糖值测定装置，其特征在于，所述数据解析部利用由所述脉搏波测定部进行测定所得的加速度脉搏波的标准化掀搏波，且根据所述定标线而计算出所述血糖值，所述标准化脉搏波是对所述加速度脉搏波进行其波形位移的标准化、以及、波形长度的标准化而获得的单脉冲的波形数据。

百度查询： 国立大学法人信州大学 非侵入式血糖值测定方法以及非侵入式血糖值测定装置

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD