基于阻抗谱法的无创血糖检测系统设计-李洋.pdf

《基于阻抗谱法的无创血糖检测系统设计-李洋.pdf》由会员分享，可在线阅读，更多相关《基于阻抗谱法的无创血糖检测系统设计-李洋.pdf（4页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、2016年第35卷第3期 传感器与微系统(Transducer and Microsystem Technologies) 77DOI：1013873J1000-9787(2016)03-0077-03基于阻抗谱法的无创血糖检测系统设计木李 洋，唐飞，李曙哲，王晓浩(清华大学精密仪器系精密测试技术及仪器国家重点实验室。北京100084)摘要：以生物阻抗检测技术为基础，设计了一种频谱范围在10-60 MHz的高频阻抗检测系统用于人体阻抗测试，验证了系统的测试性能。根据阻抗谱法无创血糖测试的原理，对健康人体进行测试后得到112组数据，采用主成分分析回归模型预测血糖浓度，对比同一测试者的预测结果与有

2、创检测数据，相关系数为08058，且较好预测了动态血糖变化趋势，为进一步研究无创连续监测血糖浓度提供了参考。关键词：无创血糖检测；阻抗谱；主成分分析中图分类号：TM 930 文献标识码：A 文章编号：1000_9787(2016)03-0077-03一 n - Design 0I a non-lnVaSlVe glucose detecting systembased 0n impedance pectroscopy methodspectroscopymethod半0nLI Yang，TANG Fei，LI Shuzhe，WANG Xiaohao(State Key Laboratory o

3、f Precision Measurement Technology and Instruments，Department of Precision Instrument，Tsinghua University，Beijing 100084，China)Abstract：A high frequency impedance monitoring system based on impedance spectroscopy has been developed，whose frequency ranges of 1060 MHz，applied for human body testAccord

4、ing to principle of non-invasiveglucose measuring based on impedance spectroscopy，after testing on healthy body，1 12 groups of data arecollected；use principal component analysis regression model to predict glucose concentration，compare predictingresults and invasive data from the s,锄tne tester，corre

5、lation coefficient is 0805 8，dynamic glucose change trend canbe well predicted，provide reference for further study of noninvasive continuous monitoring of glucose levelsKey words：noninvasive glucose detecting；impedance spectroscopy；principal component analysisO引 言糖尿病是以高血糖为特征、可引起多种并发症的慢性代谢疾病，其治疗需要不断对

6、人体的血糖水平进行检测。目前常用的血糖检测方法多为有创或微创方法，给患者造成疼痛的同时，增加了感染风险，限制了血糖浓度检测的频率和效果。因此，以光谱法、反离子渗透法、代谢热整合法等为代表的无创血糖检测，因其便捷、安全、实时等优势而广受关注。o。近年来，国外多个研究机构将生物阻抗检测技术与无创血糖测试相结合，开辟了无创血糖检测的新方向。瑞士的Caduff研究小组研究了1200 MHz范围的人体阻抗变化，证明在较高频率段，血糖浓度与阻抗值之间存在较为明显的关系H1；Harry Rechardson Elden等人设计的个人监护仪采用单频率测量，得到血糖浓度=031阻抗模值+O24相角的关系【5j。

7、韩国的Kiseok Song研究组综合阻收稿日期：2015-05-26女基金项目：国家“863”计划资助项目(2011AA040406)抗谱法与红外光谱法进行无创血糖测试，研究的频率范围在IO76kHz6|。基于生物阻抗测试原理，本文设计了基于STM32F405芯片的人体阻抗高频测试系统，该系统功耗低、性能稳定，能输出10-60 MHz的激励信号对人体进行阻抗检测，得到人体阻抗频谱参数，在此基础上，通过主成分分析回归模型表征血糖浓度的变化。1 阻抗谱法无创检测原理血液与组织中的血糖浓度改变会引起人体内相应的电解质浓度变化，从而打破血液与组织液之间的电平衡。当穿过细胞膜的离子浓度发生变化时，细胞

8、膜的电特性也发生相应改变，宏观上体现为人体阻抗变化。在100 kHz100MHz频段(13频散)，细胞的电性能变化主要由膜电容的容性短路和细胞内外大分子的转运引起71，研究表明，此时病变组织与健康组织介电常数对比较为显著。糖尿病万方数据78 传感器与微系统 第35卷患者的血样阻抗幅频曲线向低频平移，明显降低；相频曲线中负相移明显增加。人体组织既包含有阻性成分又包含有容性成分，常用Cole-Cole模型来等效，其经验公式表示为式(1)z巩+高篑，式中z为复电阻抗，R。为频率无穷大时的电阻，凡为频率为0 Hz时电阻，为频率Z为特征频率，a为散射系数。阻抗频谱表示在复平面上为第一象限上的一段圆弧，通

9、过拟合圆弧的半径、圆心坐标，可得到人体电特性参数。在本高频阻抗测试系统中，为了得到较为明显的电特性参数，系统通过一个外加的电感器与人体构成谐振电路，等效阻抗测试电路如图1。 ， ；宁上图1等效测试电路Fig 1 Equivalent test circuit当激励信号为正弦信号时，输入输出分别表示为：Vi。=Asin tot，vo。=Asin(tot+。两信号之间的幅值比k和相移6可通过幅相检测器得到。通过式(2)求得阻抗模值IzI和相位妒=万菰kR；io丽,妒-8-tani丽ksin 6(2)2测试系统设计系统以STM32F405芯片作为控制核心，控制DDS(direct di画tM syn

10、thesizer)在特定频率范围内按照预定间隔扫频；收集幅相检测器得到的电极测试信号并进行AD转换；通过串口与上位机通信。总体结构如图2所示。、r板l乜惝 歧耿图2测试系统硬件结构示意图Fig 2 Hardware structure of measuring system21硬件设计211主板设计阻抗测试模块的主电路包括以STM32为核心的主控模块、以AD9913为核心的DDS信号源、以幅相测试芯片AD8302为核心的检测模块以及电源、放大滤波模块和数据存储模块。系统采用5V供电，主芯片STM32F405拥有32位内核，运行时最高频率可达168 MHz，内部自带3个12位带有DMA(dire

11、ct memoryaccess)功能的AD转换，丰富的端口方便地为微处理器与其他外设通信提供了基础。扫描信号的产生依赖于AD9913，频率分辨率达0058 Hz，相位精度达0022。幅相测试芯片ADS302可测量从低频到27GHz的频率范围下，增益比为-30+30dB、相位差为-180000或0。-1800范围内的两信号之间的幅值比和相位差。DDS发出的初始波形通过后级的滤波电路滤除噪声后再进入高速放大器OPA691进行放大，参考信号与电极返回的测试信号之间通过电压跟随器0PA2652隔离，分别输入到幅相检测芯片AD8302进行检测，检测结果以电压信号的形式通过STM32自带的AD转换控制器转

12、换后输出结果。212测试电极系统的测试电极采用两端子的指状电极，基底尺寸为30mm x20咖，采用柔性材料以便与人体皮肤全面接触，指状电极采用表面镀金工艺，利于微弱电流的传导。电极片表面覆膜，为谐振电路提供了较为稳定的电容，提高测试系统的测试灵敏度和准确性H 0|，并能对人体起到一定保护作用。22软件设计221阻抗检测流程设计阻抗检测控制程序在EWARM(IAR E删ded Workbench for ARM)环境下用C语言编写，包括上位机通信程序与系统检测主控程序。通信程序通过串口接收上位机命令，执行相应的测试和数据传输等操作，检测主控程序的流程如图3。图3 系统测试程序流程图Fig 3 F

13、low chart of system test program222数据转换与处理系统每次扫频循环连续测量3组数据，对3组数据取平均以减小误差，所得阻抗频谱滤波后，对谐振频率附近数万方数据第3期 李洋，等：基于阻抗谱法的无创m糖检测系统设计 79据进行多项式拟合，避免了离散点之间数值跳动大、取值不稳定的问题，得到更为准确的谐振点处的阻抗信息(z赫，五，吼)。23系统性能测试采用高精度的无感电阻、电容模拟人体阻抗与选定的村圈电感器组成谐振电路分别进行测试。在阻抗最低点，实部X=z。e口，虚部Y=Z。i。sin p，等效的电容C和电阻尺通过式(3)计算得到c=了：iL二-7y陌oI，尺=王竺等手

14、兰芷+x (3)其中，=2砺，为阻抗最低点对应的角频率，L为谐振电感。系统测试表明，按照仿真得到的理论阻抗曲线与测试所得频谱曲线基本重合(图4)。计算所得等效电阻值与实际接入电路的无感电阻相对误差在05以内，电容值相对误差在15以内，满足人体阻抗测试的要求。25 30 3 5 40 4 5 5 0 5 5频率101tz圈4谐振电路阻抗频谱曲线与理论曲线对比圈Fig 4 Contrast最郾砖of妇叩酬la矗传spectrum almaf rcsonmt circuit and theoretical cum3人体血糖罚试实验31 实验方法以健康人群为实验对象，用设计的系统对20一30岁之问的l

15、O名志愿者分别进行多次了测试。被测者早晨8点前空腹到达室内实验地点，静坐l h后口服定量的葡萄糖溶液(509葡萄糖溶于250mL水)，并对其阻抗与血糖浓度进行连续监测。参考血糖浓度采用AecuChek Performa罗氏血糖仪指尖采血得到。实验时被测试者保持静坐姿势，被测手臂无大幅运动，通过弹性带将柔性电极片固定在小臂上，并保持指状电极轴向平行于小臂轴向。32 实验结果实验共得到112组有效阻抗数据及其相应的血糖浓度。采用统计分析软件SPSS(statistical pdtJct and servicesolutions)将包含有复阻抗特征参数(电导c0=llRo，G。=1R。，Ol Z)、

16、阻抗最低点Z。及其对应的频率点，等效的电阻R、电容C等多个参数在内的信息初步筛选后进行主成分分析。相应的主成分因子的载荷矩阵如表l所示。表1主成分园子载荷矩阵Tab 1 Load matrix of mBin compo越nt factor主成分l中足和电导G。，Go所占比重较大，该主成分侧重于表现阻抗信息中的阻性成分；主成分2中五所占比重较大，决定大小的是人体阻抗的容性成分。对主成分变量做最小二乘回归，还原到原始变量表示为式(4)glucose=14810气-026R+242912Go一529724G。-14109 (4)用该模型对同一测试对象25组数据进行分析，阻抗谱法的血糖估计值与参考血

17、糖浓度之间的相关系数为08058，二者对比如图5所示。表2为测试数据对比。无创检测能够较好地预测血糖变化趋势，为动态血糖监测的进一步研究提供了基础。图5Fig 5f川20ram计算得到的血糖浓度与有创血糖测试结果对比Contrast of calculated glucose concentration andinvasive data表2测试数据对比结果1Ib 2 Comparison result of test data序标称值 测试值绝对误差相对误差 序标称值 测试值绝对误差相对误差号(耐dL)(吖dL)(叫dL)喁 号(mgdL)(rrdL)(mdL) 1 7200 9762 256

18、2 35 59 lO 72 00 7230 030 0422 7740 9485 1745 22 54 Il 77 40 7292 4 48 5793 l皓00 125昭 1783 165l 12 10260 105 74 3 14 3064 14760 l蛤59 240l 1627 13 11160 10757 4 03 3 615 12960 12122 8鹄 646 14 9720 10494 774 7966 11340 9304 2036 1795 15 95如 10789 1249 13097 12960 113 76 15 84 12228 l 眈70 810 803 24 66

19、60 7894 1234 18539 9540 9563 023 024 25 66舳 81强 14n 221l4结束语本文设计了一种高频条件下人体阻抗的测试系统，验证了系统阻抗测试的性能。通过对健康人体的测试，利用主成分分析回归建立了阻抗谱法血糖预测模型，与有创测(下转第82页)如如如如m加如钙弱如笛加万方数据82 传感器与微系统 第35卷图4 SI值计算流程图4 Calculation flow chart of SI value噬l芝。图5 SI值计算误差比较Fig 5 calculation error value comparison of SI value震阀门控制器并通过无线传输

20、到监控中心。整套系统安装于唐山燃气集团天元社区实验站呻。在进行试运行期间，在测得sI值超过40 cms时，传感器发出报警信号，作为判断地震动强度的参考数据。4结论地震动sI传感器采用z变换方法得到的差分方程计算地震相对速度反应，简化的SI计算方法可以保证在一定精度的前提下实时计算SI值。传感器实时输出的最大加速度值和SI值可以作为判断地震动强弱的参考数据。参考文献：IGBT17742中国地震烈度表S20082 郝敏，谢礼力集集地震等震线和PGA、PGV等值关系的研究J地震工程与工程振动，2006，26(1)：19-213胡聿贤地震工程学M2版北京：地震出版社，2006：135-1374 Kar

21、im Kazi R，Yamazaki FumioCorrelation of JMA instrumentalseismic intensity with strong motion parametersJEarthquakeEngineering and Structural Dynamics，2002，31：1 191-12125谢礼立，于双久强震观测与分析原理M北京：地震出版社。1982：219-2216 马强，金星，李山有单自由度系统地震动力反应的实时计算方法J地震工程与工程振动，2003，船(5)：61-637许亮华，郭永刚，杜修力水工建筑物强震动加速度计录的分析处理J地震工程与工程

22、振动，2012，32(5)：28-308杨学山，马树林，匙庆磊，等城市燃气管网地震安全监控和紧急处理系统及其应用研究R哈尔滨：中国地震局工程力学研究所，2015：128作者简介：王雷(1972一)，男，黑龙江哈尔滨人，博士，副研究员，主要研究方向为数字强震仪、地震计标定与测试。q妒pu)pqpt声t声q妒t户pu、u，q妒户t)q)PPqp、妒声q声t声u妒q妒t声q)声pu)ppq)pq)pq声qq妒qppk)pu妒)u妒)pk妒qpup(上接第79页)量结果作对比，得到相关系数为08058，验证了血糖浓度与人体阻抗参数之间的相关性，反映了血糖浓度的动态变化趋势，为进一步研究无创连续监测血糖

23、浓度提供了参考。参考文献：1Malin S F，Ruchti T L，Blank T B，et a1Noninvasive prediction ofglucose by nearinfrared diffuse reflectance spectmscopyJClinical Chemistry，1999，45(9)：1651-16582Tiemey M J，Tamada J A，Ports R O，et a1Clinical evaluation ofthe GlueoWatchbiographer：A continual，noninvasive glueosemonitor for pa

24、tients with diabetesJBiosensors and Bioeleetmnics，2001，16(9)：621-6293 唐飞，王晓浩。王东生，等代谢热整合法无创血糖检测技术研究J仪器仪表学报，2007(10)：1858-18594Caduff A，Hilt E，Feldman Y，et a1First human experiments witha novel noninvasive，nonoptical continuous glucose monitoringsystemJBiosensors and Bioeleetmnics，2003，19(3)：2092175El

25、den H R，Wiekett R R，Ollmar SMethod and apparatus fornon-invasive determination of glucoin body fluids：US。6，517，482P12003-02-116Song K，Ha U，Park s，et a1An impedance and multi-wavelengthneaY-infrared spectroscopy IC for noninvasive blood glucose esti-mationCff2014 Symposium on VLSI Circuits，IEEE，2014：

26、l一27Pethig RDielectric properties of body tissuesJClinical Physics and Physiological Measurement，1987(8)：58 胡茂清，黄 华，袁支润，等血液组成交流电特性的研究J生物医学工程学杂志，2006(23)：36-409刘锐岗，史学涛，付峰，等多频电阻抗扫描的ColeCole模型分析J航天医学与医学工程，2006(6)：1110Caduff A，Dewarrat F，Talary M，et a1Non-invasive glueose monitoting in patients with diabetes：A novel system based on impedance spectroscopyJBiosensors and Bioelectronies，2006(22)：598-604作者简介：李洋(1989)，女，甘肃平凉人，硕士研究生，研究方向为科学仪器、无创检测系统。唐飞。通讯作者，Email：tangfmailtsinghuaeduca。瑚m舳。万方数据