CN104905787A - 心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置及方法 - Google Patents
心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置及方法 Download PDFInfo
- Publication number
- CN104905787A CN104905787A CN201510218250.7A CN201510218250A CN104905787A CN 104905787 A CN104905787 A CN 104905787A CN 201510218250 A CN201510218250 A CN 201510218250A CN 104905787 A CN104905787 A CN 104905787A
- Authority
- CN
- China
- Prior art keywords
- impedance
- data
- electrical impedance
- module
- parameter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置及方法，包括第一电阻抗测量模块、第二电阻抗测量模块和上位机；上位机控制第一电阻抗测量模块实时动态地采集心脏附近的单频电阻抗数据，由上位机识别出心阻抗模值的极值，该极值包括极大值和极小值，并基于该极值控制第二电阻抗测量模块按照设定的参数对外周组织进行电阻抗频谱数据测量，上位机对电阻抗频谱数据进行处理得到电阻抗特征参数，上位机计算在心阻抗模值的极大值和极小值两次驱动下测得电阻抗特征参数的平均值作为外周组织电阻抗频谱特性的评估参数，并显示该评估参数。本发明能够实现稳定地评估外周组织的电阻抗频谱特性。
Description
技术领域
本发明属于电阻抗频谱特性测量技术，具体涉及一种心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置及方法。
背景技术
生物组织的电阻抗频谱特性是指生物组织的电阻抗的阻性和容性成份随着所施加激励信号的频率的变化而发生明显的变化。而电阻抗频谱特性可以通过在被测目标表面以扫频的方式施加具有一定幅值的交流电流，并测量相应检测电极上的响应电压来获得被测目标的电阻抗实虚部的频谱数据，然后可绘制出电阻抗虚部和实部的散点图，通过曲线拟合获得Cole-Cole圆弧，借助于几何学的运算可计算出四个电阻抗特征参数R0、R∞、α和fc(R0代表直流阻抗，R∞代表频率为无穷大时的阻抗，α代表散射因子，fc代表特征频率)，从而实现定量评估。现有评估人体外周组织电阻抗频谱特性的技术方法没有考虑到心脏射血周期的影响，只通过任意时间点的外周组织电阻抗频谱特性来评估(如：US 8750978B2)，其评估结果的重复性和稳定性不佳，因此需要将心阻抗信号的引入到外周组织电阻抗频谱测量的过程中，以实现稳定地评估外周组织的电阻抗频谱特性。
发明内容
本发明的目的是提供一种心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置及方法，以实现稳定地评估外周组织的电阻抗频谱特性。
本发明所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置，包括第一电阻抗测量模块、第二电阻抗测量模块和上位机；
所述第一电阻抗测量模块基于上位机的控制实时动态地采集心脏附近的单频电阻抗数据，并利用单频电阻抗数据的实部和虚部计算心阻抗模值后传输给上位机，由上位机识别出心阻抗模值的极值，该极值包括极大值和极小值，并基于该极值控制第二电阻抗测量模块按照设定的参数对外周组织进行电阻抗频谱数据测量，并将电阻抗频谱数据传给上位机进行存储，上位机对所述电阻抗频谱数据进行处理得到电阻抗特征参数，并将该电阻抗特征参数进行存储，上位机计算在心阻抗模值的极大值和极小值两次驱动下测得电阻抗特征参数的平均值作为外周组织电阻抗频谱特性的评估参数，并显示该评估参数。
所述上位机包括第一数据采集控制模块、极值识别模块、极值驱动模块、第二数据采集控制模块、数据存储模块、电阻抗特征参数提取模块、参数存储模块、两次测量参数平均值计算模块、平均值显示模块；
所述第一数据采集控制模块将单频电阻抗数据采集的参数传递给第一电阻抗测量模块，该单频电阻抗数据采集的参数包括激励电流的幅值Am、频率f1和数据采集频率f_acq；
所述第一电阻抗测量模块包括第一电阻抗数据采集子模块、第一接口和第一电阻抗电极，第一电阻抗电极与第一电阻抗数据采集子模块连接，第一电阻抗数据采集子模块通过第一接口与第一数据采集控制模块连接；所述第一接口接受第一数据采集控制模块的控制，将施加给目标的激励电流的幅值Am、频率f1和数据采集频率f_acq传递给第一电阻抗数据采集子模块，所述第一电阻抗数据采集子模块通过第一电阻抗电极将激励电流施加到被测目标的心脏附近位置，第一电阻抗电极中的电压信号被第一电阻抗数据采集子模块按照设定的数据采集频率采集，第一电阻抗数据采集子模块利用所采集的单频电阻抗数据的实部和虚部计算出心阻抗模值后通过第一接口发送给极值识别模块；
所述极值识别模块将第一电阻抗测量模块传递过来的实测心阻抗模值与心阻抗模值的极大或极小先验值进行比较，将差异最小的实测心阻抗模值作为心阻抗极值，并将极值驱动标志f_drv置为真，该极值识别模块与第一接口连接；
所述极值驱动模块在极值驱动标志f_drv为真时启动第二数据采集控制模块，该极值驱动模块与极值识别模块连接；
所述第二数据采集控制模块将电阻抗频谱数据测量的参数传递给第二电阻抗测量模块，该电阻抗频谱数据测量的参数包括激励电流的幅值Am、起始频率f0、频率间隔Δf和激励次数N；
所述第二电阻抗测量模块包括第二电阻抗数据采集子模块、第二接口和第二电阻抗电极；第二电阻抗电极与第二电阻抗数据采集子模块连接，第二电阻抗数据采集子模块通过第二接口分别与第二数据采集控制模块连接；所述第二接口接受第二数据采集控制模块的控制，将施加给目标的激励电流的幅值Am、起始频率f0、频率间隔Δf和激励次数N传递给第二电阻抗数据采集子模块，第二电阻抗数据采集子模块通过第二电阻抗电极按照扫频测量方式将激励电流逐个施加到被测目标的肢端附近位置，第二电阻抗电极中的电压信号被第二电阻抗数据采集子模块采集，并通过第二接口发送给数据存储模块；
所述数据存储模块将第二电阻抗测量模块传递过来的一次扫频测量测得的N个频点的电阻数据的实部和虚部分别存储于4行N列的矩阵Ma1中前两行或后两行，该数据存储模块与第二接口连接；
所述电阻抗特征参数提取模块利用矩阵Ma1中的前两行或后两行数据作为N个电阻抗数据的实部和虚部，将N个实部作为X变量，将N个虚部作为Y变量，利用最小二乘法拟合出两者之间的圆弧状曲线，并计算该曲线的几何参数得到电阻抗特征参数，该电阻抗特征参数包括R0、R∞、α和fc，其中，R0代表直流阻抗，R∞代表频率为无穷大时的阻抗，α代表散射系数，fc代表特征频率，该电阻抗特征参数提取模块与数据存储模块连接；
所述参数存储模块将电阻抗特征参数提取模块所获得的电阻抗特征参数存储于2行4列的矩阵Ma2中，该参数存储模块与电阻抗特征参数提取模块连接；
所述两次测量参数平均值计算模块用于求心阻抗模值取极大值时所对应的电阻抗参数与心阻抗模值取极小值时所对应的电阻抗参数的平均值，该两次测量参数平均值计算模块与参数存储模块连接；
所述平均值显示模块用于显示四个电阻抗特征参数在两次心阻抗极值驱动中外周组织电阻抗特征参数的平均值，平均值显示模块与两次测量参数平均值计算模块连接。
所述数据采集频率f_acq是第一电阻抗测量模块在一秒钟内所采集单频心阻抗数据的次数。
所述心阻抗模值的极大先验值为从所述第一电阻抗测量模块在一分钟内所采集的心阻抗数据中提取出的电阻抗模值的极大值；
所述心阻抗模值的极小先验值极为从所述第一电阻抗测量模块在一分钟内所采集的心阻抗数据中提取出的电阻抗模值的极小值。
所述扫频测量是第二电阻抗数据采集子模块按照预设的电阻抗数据采集的参数将N个激励信号逐个施加给目标，每施加一个激励信号，测得一组电阻抗实部和虚部数据。
所述电阻抗频谱数据是指一次扫频测量获得的所有电阻抗数据的实部和虚部。
所述第一接口和第二接口采用USB接口。
本发明所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估方法，采用本发明所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置，所述第一电阻抗测量模块基于上位机的控制实时动态地采集心脏附近的单频电阻抗数据，并利用单频电阻抗数据的实部和虚部计算心阻抗模值后传输给上位机，由上位机识别出心阻抗模值的极值，该极值包括极大值和极小值，并基于该极值控制第二电阻抗测量模块按照设定的参数对外周组织进行电阻抗频谱数据测量，并将电阻抗频谱数据传给上位机进行存储，上位机对所述电阻抗频谱数据进行处理得到电阻抗特征参数，并将该电阻抗特征参数进行存储，上位机计算在心阻抗模值的极大值和极小值两次驱动下测得电阻抗特征参数的平均值作为外周组织电阻抗频谱特性的评估参数，并显示该评估参数。
包括以下步骤：
步骤1、第一数据采集控制模块将单频电阻抗数据采集的参数传递给第一电阻抗测量模块，该单频电阻抗数据采集的参数包括激励电流的幅值Am、频率f1和数据采集频率f_acq；
步骤2、第一接口接受第一数据采集控制模块的控制，将施加给目标的激励电流的幅值Am、频率f1和数据采集频率f_acq传递给第一电阻抗数据采集子模块，所述第一电阻抗数据采集子模块通过第一电阻抗电极将激励电流施加到被测目标的心脏附近位置，第一电阻抗电极中的电压信号被第一电阻抗数据采集子模块按照设定的数据采集频率采集，第一电阻抗数据采集子模块利用所采集的单频电阻抗数据的实部和虚部计算心阻抗模值后通过第一接口发送给极值识别模块；
步骤3、极值识别模块将第一电阻抗测量模块传递过来的实测心阻抗模值与心阻抗模值的极大或极小先验值进行比较，将差异最小的实测心阻抗模值作为心阻抗极值，并将极值驱动标志f_drv置为真；
步骤4、极值驱动模块在极值驱动标志f_drv为真时启动第二数据采集控制模块；
步骤5、第二数据采集控制模块将电阻抗频谱数据测量的参数传递给第二电阻抗测量模块，该电阻抗频谱数据测量的参数包括激励电流的幅值Am、起始频率f0、频率间隔Δf和激励次数N；
步骤6、第二接口接受第二数据采集控制模块的控制，将施加给目标的激励电流的幅值Am、起始频率f0、频率间隔Δf和激励次数N传递给第二电阻抗数据采集子模块，第二电阻抗数据采集子模块通过第二电阻抗电极按照扫频测量方式将激励电流逐个施加到被测目标的肢端附近位置，第二电阻抗电极中的电压信号被第二电阻抗数据采集子模块采集，并通过第二接口发送给数据存储模块；
步骤7、数据存储模块将第二电阻抗测量模块传递过来的一次扫频测量测得的N个频点的电阻数据的实部和虚部分别存储于4行N列的矩阵Ma1中前两行或后两行；
步骤8、电阻抗特征参数提取模块利用矩阵Ma1中的前两行或后两行数据作为N个电阻抗数据的实部和虚部，将N个实部作为X变量，将N个虚部作为Y变量，利用最小二乘法拟合出两者之间的圆弧状曲线，并计算该曲线的几何参数得到电阻抗特征参数，该电阻抗特征参数包括R0、R∞、α和fc，其中R0代表直流阻抗，R∞代表频率为无穷大是的阻抗，α代表散射系数，fc代表特征频率；
步骤9、参数存储模块将电阻抗特征参数提取模块所获得的电阻抗特征参数存储于2行4列的矩阵Ma2中；
步骤10、两次测量参数平均值计算模块用于求心阻抗模值取极大值时所对应的电阻抗参数与心阻抗模值取极小值时所对应的电阻抗参数的平均值；
步骤11、平均值显示模块用于显示四个电阻抗特征参数在两次心阻抗极值驱动中外周组织电阻抗特征参数的平均值。
所述第一电阻抗电极、第二电阻抗电极采用两电极法或四电极法。
本发明具有以下优点：利用心阻抗极值来驱动外周组织电阻抗频谱特性的测量和评估，能够确保评估结果的重复性和稳定性较佳。利用本发明对人体外周组织(前臂)的电阻抗频谱特性进行20次重复测试和评估，四个特征参数的变异系数(即标准差/均数)均低于0.2％。
附图说明
图1是本发明的原理框图。
图2是本发明的方法的流程图。
具体实施方式
下面结合附图对本发明作进一步说明。
如图1所示的心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置，包括第一电阻抗测量模块1、第二电阻抗测量模块2和上位机3。第一电阻抗测量模块1基于上位机3的控制实时动态地采集心脏附近的单频电阻抗数据，并利用单频电阻抗数据的实部和虚部计算心阻抗模值后传输给上位机3，由上位机3识别出心阻抗模值的极值，该极值包括极大值和极小值，并基于该极值控制第二电阻抗测量模块2按照设定的参数对外周组织进行电阻抗频谱数据测量，并将电阻抗频谱数据传给上位机3进行存储，上位机3对所述电阻抗频谱数据进行处理得到电阻抗特征参数，并将该电阻抗特征参数进行存储，上位机3计算在心阻抗模值的极大值和极小值两次驱动下测得电阻抗特征参数的平均值作为外周组织电阻抗频谱特性的评估参数，并显示该评估参数。
如图1所示，上位机3包括第一数据采集控制模块31、极值识别模块32、极值驱动模块33、第二数据采集控制模块34、数据存储模块35、电阻抗特征参数提取模块36、参数存储模块37、两次测量参数平均值计算模块38、平均值显示模块39。
第一数据采集控制模块31将单频电阻抗数据采集的参数传递给第一电阻抗测量模块1，该单频电阻抗数据采集的参数包括激励电流的幅值Am、频率f1和数据采集频率f_acq，该数据采集频率f_acq是第一电阻抗测量模块1在一秒钟内所采集单频心阻抗数据的次数，本实例中幅值Am取值1.25mA，频率f1取值50kHz，数据采集频率f_acq取值3Hz。
第一电阻抗测量模块1包括第一电阻抗数据采集子模块12、第一接口11和第一电阻抗电极13，第一电阻抗电极13与第一电阻抗数据采集子模块12连接，第一电阻抗数据采集子模块12通过第一接口11与第一数据采集控制模块31连接。本实例中，第一电阻抗测量模块1利用基于AD5933的复阻抗测量仪(AnalogDevice Inc.,ADI)的快速模式实现，第一接口11采用USB接口。所述第一接口11接受第一数据采集控制模块31的控制，将施加给目标的激励电流的幅值Am、频率f1和数据采集频率f_acq传递给第一电阻抗数据采集子模块12，所述第一电阻抗数据采集子模块12通过第一电阻抗电极13将激励电流施加到被测目标的心脏附近位置，第一电阻抗电极13中的电压信号被第一电阻抗数据采集子模块12按照设定的数据采集频率采集，第一电阻抗数据采集子模块12利用所采集的单频电阻抗数据的实部和虚部计算出心阻抗模值后通过第一接口11发送给极值识别模块32。
极值识别模块32将第一电阻抗测量模块1传递过来的实测心阻抗模值与心阻抗模值的极大或极小先验值进行比较，将差异最小的实测心阻抗模值作为心阻抗极值，并将极值驱动标志f_drv置为真，该极值识别模块32与第一接口11连接。所述心阻抗模值的极大先验值为从所述第一电阻抗测量模块1在一分钟内所采集的心阻抗数据中提取出的电阻抗模值的极大值。所述心阻抗模值的极小先验值极为从所述第一电阻抗测量模块1在一分钟内所采集的心阻抗数据中提取出的电阻抗模值的极小值。
极值驱动模块33在极值驱动标志f_drv为真时启动第二数据采集控制模块34，该极值驱动模块33与极值识别模块32连接。
第二数据采集控制模块34将电阻抗频谱数据测量的参数传递给第二电阻抗测量模块2，该电阻抗频谱数据测量的参数包括激励电流的幅值Am、起始频率f0、频率间隔Δf和激励次数N。本实例中幅值Am取值1.25mA，起始频率f0取值1kHz，频率间隔Δf取值1kHz，激励次数N取值100。
第二电阻抗测量模块2包括第二电阻抗数据采集子模块22、第二接口21和第二电阻抗电极23；第二电阻抗电极23与第二电阻抗数据采集子模块22连接，第二电阻抗数据采集子模块22通过第二接口21分别与第二数据采集控制模块34连接。本实例中，第二电阻抗测量模块2利用AD5933的复阻抗测量仪(AnalogDevice Inc.,ADI)的正常模式实现，第二接口21采用USB接口。第二接口21接受第二数据采集控制模块34的控制，将施加给目标的激励电流的幅值Am、起始频率f0、频率间隔Δf和激励次数N传递给第二电阻抗数据采集子模块22，第二电阻抗数据采集子模块22通过第二电阻抗电极23按照扫频测量方式(即上述正常模式)将激励电流逐个施加到被测目标的肢端附近位置，所述扫频测量是第二电阻抗数据采集子模块22按照预设的电阻抗数据采集的参数将N个(本实例中N取值100)激励信号逐个施加给目标，每施加一个激励信号，测得一组电阻抗实部和虚部数据。第二电阻抗电极23中的电压信号被第二电阻抗数据采集子模块22采集，并通过第二接口21发送给数据存储模块35。
第一电阻抗电极13是第一电阻抗测量模块1向被测目标施加激励电流信号和采集响应电压信号的传感器。第二电极阻抗23是第二电阻抗测量模块2向被测目标施加激励电流信号和采集响应电压信号的传感器。第一电阻抗电极13和第二电极阻抗23可以基于两电极法或四电极法分配电极的功能，两电极法有两个电极既作为激励电极又作为测量电极，四电极法有四个电极，位于外侧的两个电极作为激励电极，而位于内侧两个电极作为测量电极。本实例中的第一电阻抗测量模块1和第二电阻抗测量模块2均采用两电极法。第一电阻抗电极13和第二电极阻抗23均由摁扣式心电电极线和两片一次性心电电极组成，其中电极线的摁扣端与一次性心电电极片连接，另外一端与AD5933复阻抗测量仪的一个测量端口接线柱连接，组成了电阻抗测量模块。第一电阻抗电极13的两个电极贴于胸骨外的皮肤上，它们的位置分别位于颈根到剑突的线段中点往上或往下5cm。第二电极阻抗23的两个电极贴于上臂中段的皮肤上，两个电极中心连线长4cm，线段的中点位于前臂的中点。
数据存储模块35将第二电阻抗测量模块2传递过来的一次扫频测量测得的N个频点的电阻数据的实部和虚部分别存储于4行N列(N为100)的矩阵Ma1中前两行或后两行，该数据存储模块35与第二接口21连接。
电阻抗特征参数提取模块36利用矩阵Ma1中的前两行或后两行数据作为N个(N为100)电阻抗数据的实部和虚部，将N个实部作为X变量，将N个虚部作为Y变量，利用最小二乘法拟合出两者之间的圆弧状曲线，并计算该曲线的几何参数得到电阻抗特征参数，该电阻抗特征参数包括R0、R∞、α和fc，其中，R0代表直流阻抗，R∞代表频率无穷大时的阻抗，α代表散射因子，fc代表特征频率，该电阻抗特征参数提取模块36与数据存储模块35连接。
参数存储模块37将电阻抗特征参数提取模块36所获得的电阻抗特征参数存储于2行4列的矩阵Ma2中，该参数存储模块37与电阻抗特征参数提取模块36连接。
两次测量参数平均值计算模块38用于求心阻抗模值取极大值时所对应的电阻抗参数与心阻抗模值取极小值时所对应的电阻抗参数的平均值，该两次测量参数平均值计算模块38与参数存储模块37连接。
平均值显示模块39用于显示四个电阻抗特征参数在两次心阻抗极值驱动中外周组织电阻抗特征参数的平均值，平均值显示模块39与两次测量参数平均值计算模块38连接。
本发明所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估方法，采用本发明所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置，第一电阻抗测量模块1基于上位机3的控制实时动态地采集心脏附近的单频电阻抗数据，并利用单频电阻抗数据的实部和虚部计算心阻抗模值后传输给上位机3，由上位机3识别出心阻抗模值的极值，该极值包括极大值和极小值，并基于该极值控制第二电阻抗测量模块2按照设定的参数对外周组织进行电阻抗频谱数据测量，并将电阻抗频谱数据传给上位机3进行存储，上位机3对所述电阻抗频谱数据进行处理得到电阻抗特征参数，并将该电阻抗特征参数进行存储，上位机3计算在心阻抗模值的极大值和极小值两次驱动下测得电阻抗特征参数的平均值作为外周组织电阻抗频谱特性的评估参数，并显示该评估参数。
如图2所示，本发明所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估方法，具体包括以下步骤：
步骤1、第一数据采集控制模块31将单频电阻抗数据采集的参数传递给第一电阻抗测量模块1，该单频电阻抗数据采集的参数包括激励电流的幅值Am、频率f1和数据采集频率f_acq。
步骤2、第一接口11接受第一数据采集控制模块31的控制，将施加给目标的激励电流的幅值Am、频率f1和数据采集频率f_acq传递给第一电阻抗数据采集子模块12，所述第一电阻抗数据采集子模块12通过第一电阻抗电极13将激励电流施加到被测目标的心脏附近位置，第一电阻抗电极13中的电压信号被第一电阻抗数据采集子模块12按照设定的数据采集频率采集，第一电阻抗数据采集子模块12第一电阻抗数据采集子模块12利用所采集的单频电阻抗数据的实部和虚部计算出心阻抗模值后通过第一接口11发送给极值识别模块32。
步骤3、极值识别模块32将第一电阻抗测量模块1传递过来的实测心阻抗模值与心阻抗模值的极大或极小先验值进行比较，将差异最小的实测心阻抗模值作为心阻抗极值，并将极值驱动标志f_drv置为真。
步骤4、极值驱动模块33在极值驱动标志f_drv为真时启动第二数据采集控制模块34。
步骤5、第二数据采集控制模块34将电阻抗频谱数据测量的参数传递给第二电阻抗测量模块2，该电阻抗频谱数据测量的参数包括激励电流的幅值Am、起始频率f0、频率间隔Δf和激励次数N。
步骤6、第二接口21接受第二数据采集控制模块34的控制，将施加给目标的激励电流的幅值Am、起始频率f0、频率间隔Δf和激励次数N传递给第二电阻抗数据采集子模块22，第二电阻抗数据采集子模块22通过第二电阻抗电极23按照扫频测量方式将激励电流逐个施加到被测目标的肢端附近位置，第二电阻抗电极23中的电压信号被第二电阻抗数据采集子模块22采集，并通过第二接口21发送给数据存储模块35。
步骤7、数据存储模块35将第二电阻抗测量模块2传递过来的一次扫频测量测得的N个频点的电阻数据的实部和虚部分别存储于4行N列的矩阵Ma1中前两行或后两行。
步骤8、电阻抗特征参数提取模块36利用矩阵Ma1中的前两行或后两行数据作为N个电阻抗数据的实部和虚部，将N个实部作为X变量，将N个虚部作为Y变量，利用最小二乘法拟合出两者之间的圆弧状曲线，并计算该曲线的几何参数得到电阻抗特征参数，该电阻抗特征参数包括R0、R∞、α和fc，其中，R0代表直流阻抗，R∞代表频率为无穷大时的阻抗，α代表散射因子，fc代表特征频率。
步骤9、参数存储模块37将电阻抗特征参数提取模块36所获得的电阻抗特征参数存储于2行4列的矩阵Ma2中。
步骤10、两次测量参数平均值计算模块38用于求心阻抗模值取极大值时所对应的电阻抗参数与心阻抗模值取极小值时所对应的电阻抗参数的平均值。
步骤11、平均值显示模块39用于显示四个电阻抗特征参数在两次心阻抗极值驱动中外周组织电阻抗特征参数的平均值。
所述第一电阻抗电极13、第二电阻抗电极23采用两电极法或四电极法。
以一健康志愿者为测量对象，利用本发明方法的实施例对其外周电阻抗频谱特性的进行20次重复评估，四个特征参数的变异系数(即标准差/均数)均低于0.2％。
Claims (10)
1.一种心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置，其特征在于：包括第一电阻抗测量模块(1)、第二电阻抗测量模块(2)和上位机(3)；
所述第一电阻抗测量模块(1)基于上位机(3)的控制实时动态地采集心脏附近的单频电阻抗数据，并利用单频电阻抗数据的实部和虚部计算出心阻抗模值后传输给上位机(3)，由上位机(3)识别出心阻抗模值的极值，该极值包括极大值和极小值，并基于该极值控制第二电阻抗测量模块(2)按照设定的参数对外周组织进行电阻抗频谱数据测量，并将电阻抗频谱数据传给上位机(3)进行存储，上位机(3)对所述电阻抗频谱数据进行处理得到电阻抗特征参数，并将该电阻抗特征参数进行存储，上位机(3)计算在心阻抗模值的极大值和极小值两次驱动下测得电阻抗特征参数的平均值作为外周组织电阻抗频谱特性的评估参数，并显示该评估参数。
2.根据权利要求1所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置，其特征在于：所述上位机(3)包括第一数据采集控制模块(31)、极值识别模块(32)、极值驱动模块(33)、第二数据采集控制模块(34)、数据存储模块(35)、电阻抗特征参数提取模块(36)、参数存储模块(37)、两次测量参数平均值计算模块(38)、平均值显示模块(39)；
所述第一数据采集控制模块(31)将单频电阻抗数据采集的参数传递给第一电阻抗测量模块(1)，该单频电阻抗数据采集的参数包括激励电流的幅值Am、频率f1和数据采集频率f_acq；
所述第一电阻抗测量模块(1)包括第一电阻抗数据采集子模块(12)、第一接口(11)和第一电阻抗电极(13)，第一电阻抗电极(13)与第一电阻抗数据采集子模块(12)连接，第一电阻抗数据采集子模块(12)通过第一接口(11)与第一数据采集控制模块(31)连接；所述第一接口(11)接受第一数据采集控制模块(31)的控制，将施加给目标的激励电流的幅值Am、频率f1和数据采集频率f_acq传递给第一电阻抗数据采集子模块(12)，所述第一电阻抗数据采集子模块(12)通过第一电阻抗电极(13)将激励电流施加到被测目标的心脏附近位置，第一电阻抗电极(13)中的电压信号被第一电阻抗数据采集子模块(12)按照设定的数据采集频率采集，第一电阻抗数据采集子模块(12)利用所采集的单频电阻抗数据的实部和虚部计算出心阻抗模值后通过第一接口(11)发送给极值识别模块(32)；
所述极值识别模块(32)将第一电阻抗测量模块(1)传递过来的实测心阻抗模值与心阻抗模值的极大或极小先验值进行比较，将差异最小的实测心阻抗模值作为心阻抗极值，并将极值驱动标志f_drv置为真，该极值识别模块(32)与第一接口(11)连接；
所述极值驱动模块(33)在极值驱动标志f_drv为真时启动第二数据采集控制模块(34)，该极值驱动模块(33)与极值识别模块(32)连接；
所述第二数据采集控制模块(34)将电阻抗频谱数据测量的参数传递给第二电阻抗测量模块(2)，该电阻抗频谱数据测量的参数包括激励电流的幅值Am、起始频率f0、频率间隔Δf和激励次数N；
所述第二电阻抗测量模块(2)包括第二电阻抗数据采集子模块(22)、第二接口(21)和第二电阻抗电极(23)；第二电阻抗电极(23)与第二电阻抗数据采集子模块(22)连接，第二电阻抗数据采集子模块(22)通过第二接口(21)分别与第二数据采集控制模块(34)连接；所述第二接口(21)接受第二数据采集控制模块(34)的控制，将施加给目标的激励电流的幅值Am、起始频率f0、频率间隔Δf和激励次数N传递给第二电阻抗数据采集子模块(22)，第二电阻抗数据采集子模块(22)通过第二电阻抗电极(23)按照扫频测量方式将激励电流逐个施加到被测目标的肢端附近位置，第二电阻抗电极(23)中的电压信号被第二电阻抗数据采集子模块(22)采集，并通过第二接口(21)发送给数据存储模块(35)；
所述数据存储模块(35)将第二电阻抗测量模块(2)传递过来的一次扫频测量测得的N个频点的电阻数据的实部和虚部分别存储于4行N列的矩阵Ma1中前两行或后两行，该数据存储模块(35)与第二接口(21)连接；
所述电阻抗特征参数提取模块(36)利用矩阵Ma1中的前两行或后两行数据作为N个电阻抗数据的实部和虚部，将N个实部作为X变量，将N个虚部作为Y变量，利用最小二乘法拟合出两者之间的圆弧状曲线，并计算该曲线的几何参数得到电阻抗特征参数，该电阻抗特征参数包括R0、R∞、α和fc，其中，R0代表直流阻抗，R∞代表频率为无穷大时的阻抗，α代表散射系数，fc代表特征频率，该电阻抗特征参数提取模块(36)与数据存储模块(35)连接；
所述参数存储模块(37)将电阻抗特征参数提取模块(36)所获得的电阻抗特征参数存储于2行4列的矩阵Ma2中，该参数存储模块(37)与电阻抗特征参数提取模块(36)连接；
所述两次测量参数平均值计算模块(38)用于求心阻抗模值取极大值时所对应的电阻抗参数与心阻抗模值取极小值时所对应的电阻抗参数的平均值，该两次测量参数平均值计算模块(38)与参数存储模块(37)连接；
所述平均值显示模块(39)用于显示四个电阻抗特征参数在两次心阻抗极值驱动中外周组织电阻抗特征参数的平均值，平均值显示模块(39)与两次测量参数平均值计算模块(38)连接。
3.根据权利要求2所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置，其特征在于：所述数据采集频率f_acq是第一电阻抗测量模块(1)在一秒钟内所采集单频心阻抗数据的次数。
4.根据权利要求2所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置，其特征在于：所述心阻抗模值的极大先验值为从所述第一电阻抗测量模块(1)在一分钟内所采集的心阻抗数据中提取出的电阻抗模值的极大值；
所述心阻抗模值的极小先验值极为从所述第一电阻抗测量模块(1)在一分钟内所采集的心阻抗数据中提取出的电阻抗模值的极小值。
5.根据权利要求2所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置，其特征在于：所述扫频测量是第二电阻抗数据采集子模块按照预设的电阻抗数据采集的参数将N个激励信号逐个施加给目标，每施加一个激励信号，测得一组电阻抗实部和虚部数据。
6.根据权利要求2所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置，其特征在于：所述电阻抗频谱数据是指一次扫频测量获得的所有电阻抗数据的实部和虚部。
7.根据权利要求1或2所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置，其特征在于：所述第一接口(11)和第二接口(21)采用USB接口。
8.一种心阻抗极值驱动的外周组织电阻抗频谱特性的评估方法，其特征在于：采用如权利要求1至7任一所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置，所述第一电阻抗测量模块(1)基于上位机(3)的控制实时动态地采集心脏附近的单频电阻抗数据，并利用单频电阻抗数据的实部和虚部计算心阻抗模值后传输给上位机(3)，由上位机(3)识别出心阻抗模值的极值，该极值包括极大值和极小值，并基于该极值控制第二电阻抗测量模块(2)按照设定的参数对外周组织进行电阻抗频谱数据测量，并将电阻抗频谱数据传给上位机(3)进行存储，上位机(3)对所述电阻抗频谱数据进行处理得到电阻抗特征参数，并将该电阻抗特征参数进行存储，上位机(3)计算在心阻抗模值的极大值和极小值两次驱动下测得电阻抗特征参数的平均值作为外周组织电阻抗频谱特性的评估参数，并显示该评估参数。
9.根据权利要求8所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估方法，其特征在于，包括以下步骤：
步骤1、第一数据采集控制模块(31)将单频电阻抗数据采集的参数传递给第一电阻抗测量模块(1)，该单频电阻抗数据采集的参数包括激励电流的幅值Am、频率f1和数据采集频率f_acq；
步骤2、第一接口(11)接受第一数据采集控制模块(31)的控制，将施加给目标的激励电流的幅值Am、频率f1和数据采集频率f_acq传递给第一电阻抗数据采集子模块(12)，所述第一电阻抗数据采集子模块(12)通过第一电阻抗电极(13)将激励电流施加到被测目标的心脏附近位置，第一电阻抗电极(13)中的电压信号被第一电阻抗数据采集子模块(12)按照设定的数据采集频率采集，第一电阻抗数据采集子模块(12)利用所采集的单频电阻抗数据的实部和虚部计算心阻抗模值后通过第一接口(11)发送给极值识别模块(32)；
步骤3、极值识别模块(32)将第一电阻抗测量模块(1)传递过来的实测心阻抗模值与心阻抗模值的极大或极小先验值进行比较，将差异最小的实测心阻抗模值作为心阻抗极值，并将极值驱动标志f_drv置为真；
步骤4、极值驱动模块(33)在极值驱动标志f_drv为真时启动第二数据采集控制模块(34)；
步骤5、第二数据采集控制模块(34)将电阻抗频谱数据测量的参数传递给第二电阻抗测量模块(2)，该电阻抗频谱数据测量的参数包括激励电流的幅值Am、起始频率f0、频率间隔Δf和激励次数N；
步骤6、第二接口(21)接受第二数据采集控制模块(34)的控制，将施加给目标的激励电流的幅值Am、起始频率f0、频率间隔Δf和激励次数N传递给第二电阻抗数据采集子模块(22)，第二电阻抗数据采集子模块(22)通过第二电阻抗电极(23)按照扫频测量方式将激励电流逐个施加到被测目标的肢端附近位置，第二电阻抗电极(23)中的电压信号被第二电阻抗数据采集子模块(22)采集，并通过第二接口(21)发送给数据存储模块(35)；
步骤7、数据存储模块(35)将第二电阻抗测量模块(2)传递过来的一次扫频测量测得的N个频点的电阻数据的实部和虚部分别存储于4行N列的矩阵Ma1中前两行或后两行；
步骤8、电阻抗特征参数提取模块(36)利用矩阵Ma1中的前两行或后两行数据作为N个电阻抗数据的实部和虚部，将N个实部作为X变量，将N个虚部作为Y变量，利用最小二乘法拟合出两者之间的圆弧状曲线，并计算该曲线的几何参数得到电阻抗特征参数，该电阻抗特征参数包括R0、R∞、α和fc，其中，R0代表直流阻抗，R∞代表频率为无穷大时的阻抗，α代表散射系数，fc代表特征频率；
步骤9、参数存储模块(37)将电阻抗特征参数提取模块(36)所获得的电阻抗特征参数存储于2行4列的矩阵Ma2中；
步骤10、两次测量参数平均值计算模块(38)用于求心阻抗模值取极大值时所对应的电阻抗参数与心阻抗模值取极小值时所对应的电阻抗参数的平均值；
步骤11、平均值显示模块(39)用于显示四个电阻抗特征参数在两次心阻抗极值驱动中外周组织电阻抗特征参数的平均值。
10.根据权利要求8或9所述的心阻抗极值驱动的外周组织电阻抗频谱特性的评估方法，其特征在于：所述第一电阻抗电极(13)、第二电阻抗电极(23)采用两电极法或四电极法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510218250.7A CN104905787B (zh) | 2015-04-30 | 2015-04-30 | 心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510218250.7A CN104905787B (zh) | 2015-04-30 | 2015-04-30 | 心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置及方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104905787A true CN104905787A (zh) | 2015-09-16 |
| CN104905787B CN104905787B (zh) | 2017-06-16 |
Family
ID=54075548
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510218250.7A Expired - Fee Related CN104905787B (zh) | 2015-04-30 | 2015-04-30 | 心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104905787B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105534524A (zh) * | 2016-02-05 | 2016-05-04 | 思澜科技（成都）有限公司 | 用于甲状腺术中甲状旁腺的快速识别设备及其方法 |
| CN108852352A (zh) * | 2018-03-26 | 2018-11-23 | 深圳市麦德安医疗设备有限公司 | 生物阻抗Cole方程的参数估计方法、智能设备及存储介质 |
| CN111225727A (zh) * | 2017-07-26 | 2020-06-02 | 弗洛设计声能学公司 | 声换能器自动起动及运行 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1305777A (zh) * | 2000-01-20 | 2001-08-01 | 深圳市辉大高科技发展有限公司 | 胃动力检测方法以及采用该方法构造的检测仪 |
| US20040049237A1 (en) * | 2002-07-12 | 2004-03-11 | Larson Dennis E. | Minute ventilation sensor with dynamically adjusted excitation current |
| CN1543912A (zh) * | 2003-11-18 | 2004-11-10 | 华中科技大学 | 一种生物组织多频率阻抗测量方法及其装置 |
| CN101156776A (zh) * | 2007-09-17 | 2008-04-09 | 中国人民解放军第四军医大学 | 实时多信息提取的电阻抗扫描检测系统及方法 |
| CN102008304A (zh) * | 2010-11-10 | 2011-04-13 | 中国医学科学院生物医学工程研究所 | 一种用于电阻抗胃动力检测的信号模拟方法 |
| US7988634B1 (en) * | 2007-01-10 | 2011-08-02 | Pacesetter, Inc. | Adaptive cancellation of a signal component |
| US20120172730A1 (en) * | 2011-01-04 | 2012-07-05 | College Of William And Mary | Method and system for detecting apnea |
| US20140148712A1 (en) * | 2012-11-29 | 2014-05-29 | Siemens Medical Solutions Usa, Inc. | System and Method To Calculate Cardiac Characteristics |
| CN104114087A (zh) * | 2012-02-15 | 2014-10-22 | 皇家飞利浦有限公司 | 生物阻抗摄谱系统和方法 |
-
2015
- 2015-04-30 CN CN201510218250.7A patent/CN104905787B/zh not_active Expired - Fee Related
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1305777A (zh) * | 2000-01-20 | 2001-08-01 | 深圳市辉大高科技发展有限公司 | 胃动力检测方法以及采用该方法构造的检测仪 |
| US20040049237A1 (en) * | 2002-07-12 | 2004-03-11 | Larson Dennis E. | Minute ventilation sensor with dynamically adjusted excitation current |
| CN1543912A (zh) * | 2003-11-18 | 2004-11-10 | 华中科技大学 | 一种生物组织多频率阻抗测量方法及其装置 |
| US7988634B1 (en) * | 2007-01-10 | 2011-08-02 | Pacesetter, Inc. | Adaptive cancellation of a signal component |
| CN101156776A (zh) * | 2007-09-17 | 2008-04-09 | 中国人民解放军第四军医大学 | 实时多信息提取的电阻抗扫描检测系统及方法 |
| CN102008304A (zh) * | 2010-11-10 | 2011-04-13 | 中国医学科学院生物医学工程研究所 | 一种用于电阻抗胃动力检测的信号模拟方法 |
| US20120172730A1 (en) * | 2011-01-04 | 2012-07-05 | College Of William And Mary | Method and system for detecting apnea |
| CN104114087A (zh) * | 2012-02-15 | 2014-10-22 | 皇家飞利浦有限公司 | 生物阻抗摄谱系统和方法 |
| US20140148712A1 (en) * | 2012-11-29 | 2014-05-29 | Siemens Medical Solutions Usa, Inc. | System and Method To Calculate Cardiac Characteristics |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105534524A (zh) * | 2016-02-05 | 2016-05-04 | 思澜科技（成都）有限公司 | 用于甲状腺术中甲状旁腺的快速识别设备及其方法 |
| CN111225727A (zh) * | 2017-07-26 | 2020-06-02 | 弗洛设计声能学公司 | 声换能器自动起动及运行 |
| CN108852352A (zh) * | 2018-03-26 | 2018-11-23 | 深圳市麦德安医疗设备有限公司 | 生物阻抗Cole方程的参数估计方法、智能设备及存储介质 |
| CN108852352B (zh) * | 2018-03-26 | 2021-07-30 | 深圳市麦德安医疗设备有限公司 | 生物阻抗Cole方程的参数估计方法、智能设备及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104905787B (zh) | 2017-06-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8918170B2 (en) | Electrophysiological analysis system | |
| JP5161772B2 (ja) | インピーダンスパラメータ値 | |
| AU2008286194B2 (en) | Impedance measurement process | |
| JP2017529977A (ja) | 生体組織型の決定 | |
| CN107847174A (zh) | 起搏信号处理方法、系统和心电监护仪 | |
| US5339827A (en) | Acupuncture system and method | |
| US8594782B2 (en) | Apparatus and methods for determining the location of the apex of a dental root canal | |
| WO2015123603A1 (en) | Electrical impedance myography | |
| CA2692795A1 (en) | Cardiac monitoring system | |
| US20120323134A1 (en) | Method and system for determining a location of nerve tissue in three-dimensional space | |
| US20090312666A1 (en) | Skin conduction measuring apparatus | |
| CN104905787A (zh) | 心阻抗极值驱动的外周组织电阻抗频谱特性的评估装置及方法 | |
| JP2008514327A (ja) | 歯の構造をモニタリングする方法および装置 | |
| CN103230273B (zh) | 早期亚临床无症状糖尿病的无创检测系统及其方法 | |
| US20150038873A1 (en) | Apparatus and method for stimulator on-skin short detection | |
| CN104523271A (zh) | 一种临床用肌电测痛方法及装置 | |
| CN112423837A (zh) | 电极接触监测 | |
| WO2024067956A1 (en) | System and method for neuromonitoring based on a biosignal | |
| CN203263391U (zh) | 早期亚临床无症状糖尿病的无创检测系统 | |
| JP3984332B2 (ja) | 身体組成推計装置及び身体組成推計プログラムを記録したコンピュータ読み取り可能な記録媒体 | |
| Iusan et al. | The voltage effect on electrical parameters of the skin-ECG electrodes interface | |
| JP2017520768A (ja) | インビトロ組織における物質作用のモニタリング | |
| CN204351830U (zh) | 慢性病早期电生理检测设备和系统 | |
| JP6207130B2 (ja) | 心臓電気刺激装置 | |
| CN106419916B (zh) | 一种人体电生理参数测定装置及方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170616 Termination date: 20200430 |