CN104921736B - 一种包含参数估计功能滤波模块的连续血糖监测设备 - Google Patents
一种包含参数估计功能滤波模块的连续血糖监测设备 Download PDFInfo
- Publication number
- CN104921736B CN104921736B CN201510310303.8A CN201510310303A CN104921736B CN 104921736 B CN104921736 B CN 104921736B CN 201510310303 A CN201510310303 A CN 201510310303A CN 104921736 B CN104921736 B CN 104921736B
- Authority
- CN
- China
- Prior art keywords
- blood glucose
- value
- time
- blood sugar
- filtering
- 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.)
- Active
Links
- 239000008280 blood Substances 0.000 title claims abstract description 115
- 210000004369 blood Anatomy 0.000 title claims abstract description 115
- 238000001914 filtration Methods 0.000 title claims abstract description 52
- 238000012806 monitoring device Methods 0.000 title claims abstract description 10
- 239000008103 glucose Substances 0.000 claims abstract description 85
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 claims abstract description 82
- 238000000034 method Methods 0.000 claims abstract description 19
- 238000012544 monitoring process Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 15
- 125000002791 glucosyl group Chemical group C1([C@H](O)[C@@H](O)[C@H](O)[C@H](O1)CO)* 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 32
- 230000006870 function Effects 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 101000802640 Homo sapiens Lactosylceramide 4-alpha-galactosyltransferase Proteins 0.000 claims description 3
- 102100035838 Lactosylceramide 4-alpha-galactosyltransferase Human genes 0.000 claims description 3
- GZCGUPFRVQAUEE-SLPGGIOYSA-N aldehydo-D-glucose Chemical compound OC[C@@H](O)[C@@H](O)[C@H](O)[C@@H](O)C=O GZCGUPFRVQAUEE-SLPGGIOYSA-N 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000013256 coordination polymer Substances 0.000 claims description 3
- 238000013500 data storage Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 238000013459 approach Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 8
- 208000013016 Hypoglycemia Diseases 0.000 abstract description 5
- 230000002218 hypoglycaemic effect Effects 0.000 abstract 2
- 239000007787 solid Substances 0.000 abstract 1
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000000366 juvenile effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/14532—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring glucose, e.g. by tissue impedance measurement
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/746—Alarms related to a physiological condition, e.g. details of setting alarm thresholds or avoiding false alarms
-
- G06F19/30—
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Emergency Medicine (AREA)
- Optics & Photonics (AREA)
- Investigating Or Analysing Biological Materials (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明公开了一种包含参数估计功能滤波模块的连续血糖监测设备,该设备针对连续血糖监测信号进行滤波处理,以便进行准确的高低血糖报警。由于不同病人和不同传感器的情况不同,相应的血糖监测信号的噪声水平也不相同,从而需要设置不同的滤波参数进行血糖滤波处理。本发明设备中的滤波参数能够随着不同病人、不同传感器的改变而调整。本发明中所采用的滤波参数的估计方法能及时准确的估计出卡尔曼滤波所需要的滤波参数,从而更好地对血糖信号进行滤波处理,为提高高低血糖实时报警的准确性打下了坚实的基础,具有重要的作用。
Description
技术领域
本发明属于血糖数据处理和分析的研究领域,特别是涉及一种包含参数估计功能的滤波模块的连续血糖监测设备。
背景技术
为了管理监测血糖水平,必须要对血糖水平进行测量,目前采用的是连续血糖监测设备。随着连续血糖监测(CGM)设备的快速发展,使得更高精度的控制变为了可能,实时的CGM***对高/低血糖的提前检测起到了重要的作用。通过比较当前的测量值与高/低血糖阀值之间的关系就可以产生报警,而及时的高低血糖报警尤其是夜间的报警对于糖尿病患者尤为重要。但是通常来说,CGM数据都含有噪声,也正因为此,会对高低血糖的报警产生影响。影响CGM数据准确性的原因有很多,首先是对于自我血糖监测(SMBG)和CGM之间的校准偏差,其次是传感器中的物理、化学和电子过程中的随机误差,最后是在测量过程中产生的高频测量噪声。在一些公开的专利中一般采用滑动平均滤波的算法,但是滑动平均滤波的滤波效果不如卡尔曼滤波。
然而,在卡尔曼滤波的实际应用中,并不知道随机噪声的多少,并且信噪比会随着不同对象和不同传感器的变化而变化(个体间差异)。这时候,如果卡尔曼滤波的参数不能跟随信噪比的变化而变化,那么滤波效果就会是次优的。
发明内容
本发明的目的在于针对现有滤波方法的不足,提供一种包含参数估计功能滤波模块的连续血糖监测设备。
本发明的目的是通过以下技术方案来实现的:一种包含参数估计功能滤波模块的连续血糖监测设备,该设备包括:用于采集人体血糖信息,输出血糖可用信号的传感器;用于对传感器的输出信号进行放大处理的信号放大器;用于对信号放大器输出的模拟信号进行数字转换的单片机;用于对单片机输出的数字信号进行数据处理的滤波器,滤波器可以集成在单片机中,也可以单独使用;用于将滤波器输出的滤波后的血糖值进行数值和波形显示的显示器;用于数据存储的存储器;所述滤波器进行数据处理的过程包括以下步骤:
(1)血糖数据预处理:将以一定采样周期Δt获得的单片机输出的连续血糖监测信号组合成一维时序数据y1×l,其中,y表示检测到的血糖信号,l为采样个数,去除其中的尖峰信号;
(2)对于连续血糖监测信号建立二阶模型:用u(k)表示k时刻实际的血糖值,那么u(k)满足以下公式:
u(k)=2u(k-1)-u(k-2)+w1(k) (1)
其中w1(k)是均值为0、方差为λ2的高斯白噪声;令x1(k)=u(k),x2(k)=u(k-1),x(k)=[x1(k) x2(k)]T,由此得:
令
***观测方程如下:
y(k)=Cx(k)+v(k) (3)
其中y(k)是k时刻实际监测的血糖值,v(k)是均值为0、方差为σ2的高斯白噪声,C=[1 0];
(3)通过基于EM的方法进行模型中***噪声w(k)的协方差矩阵Q和测量噪声v(k)的协方差矩阵R的估计;
(4)根据步骤(3)的参数估计,采用卡尔曼滤波对连续血糖监测信号进行数据处理,具体包括以下子步骤:
(4.1)确定对状态提前一步的预测:
P′(k)=AP(k-1)AT+Q(k-1) (5)
其中x′(k)是对于当前时刻血糖实际值的估测,是血糖前一时刻滤波后的值,而P′(k)是估测误差的协方差矩阵,P(k-1)是血糖前一时刻实际值和滤波值之间误差的协方差矩阵;
(4.2)确定滤波后的***状态:
K=P′(k)CT[CP′(k)CT+R]-1 (6)
P(k)=(I-KC)P′(k) (8)
其中,K是卡尔曼滤波的修正矩阵,I为二阶单位矩阵,P(k)为血糖当前时刻实际值和滤波值之间误差的协方差矩阵,为血糖当前时刻滤波后的值,即最终估计的血糖值。
进一步地,所述步骤(3)具体包括以下子步骤:
(3.1)对于某个特定病人选取n个采样点的血糖信号用于EM算法的参数估计,为了下面公式表示方便,用yk表示k时刻的血糖采样值,xk表示k时刻的血糖真实值,这里认为血糖的初始值x0服从均值为μ,协方差矩阵为Σ的分布;
(3.2)根据步骤(2)建立的二阶模型采用EM算法估计参数,具体步骤如下:
(3.2.1)对于μ,Q和R赋初值μ(0),Q(0)和R(0);
(3.2.2)根据式(9)-(17)计算k=1,2,...,n
其中是k时刻的血糖滤波值,为血糖真值和滤波值之间误差的协方差矩阵,是k-1时刻对于k时刻血糖的估计值,为k时刻血糖的估计值与真实值的协方差矩阵,Kk是卡尔曼滤波增益矩阵;其中,为了计算和采用向后递归,k=n,n-1,...,1,
其中 为了下一步的计算,需要获得的值,采用向后递归的方法,k=n,n-1,...,1,
其中,
(3.2.3)根据式(18)-(20)计算Q(1)和R(1):
Q(r+1)=(W-VU-1VT)/n (19)
其中:
(3.2.4)重复步骤(3.2.2)和(3.2.3),直到估计值Q、R和对数似然函数logL稳定,其中,对数似然函数为:
与现有技术相比,本发明的有益效果是:本发明所提出的包含参数估计功能滤波模块的连续血糖监测设备能够根据不同的病人和不同的传感器估计出适当的***噪声与过程噪声,以便用于卡尔曼滤波中,其滤波效果明显好于滑动平均滤波,并能提高高低血糖的报警精度。本发明易于实施,为血糖处理和分析的研究指明了新的方向。
附图说明
图1是本发明连续血糖监测设备的结构框图;
图2是本发明连续血糖监测设备中滤波器的实现流程图;
图3是采用滑动平均滤波后的血糖数据和真实血糖值的比较图(分析对象为儿童组第五人的血糖信号);
图4是采用基于EM算法的参数估计方法的卡尔曼滤波后的血糖值和真实血糖值的比较图(分析对象为儿童组第五人的血糖信号);
图5是将两种滤波方法所得的血糖值与真实血糖值的局部放大比较图(分析对象为儿童组第五人的血糖信号)。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细说明。
如图1所示,本发明一种包含参数估计功能滤波模块的连续血糖监测设备,包括:用于采集人体血糖信息,输出血糖可用信号的传感器;用于对传感器的输出信号进行放大处理的信号放大器;用于对信号放大器输出的模拟信号进行数字转换的单片机;用于对单片机输出的数字信号进行数据处理的滤波器,滤波器可以集成在单片机中,也可以单独使用;用于将滤波器输出的滤波后的血糖值进行数值和波形显示的显示器;用于数据存储的存储器;所述滤波器进行数据处理的过程包括以下步骤:
步骤1:血糖数据预处理
将以采样周期为5分钟所获得的连续血糖监测信号组合成一维时序数据y1×l,其中,y表示检测到的血糖信号,l为样个数,去除其中的尖峰信号。本实例中,共有三组对象的采样信号,采样周期为5分钟,第1组为青少年组,第2组为成人组,第3组为儿童组,每组10人,三组共30人,每个对象的采样信号包括五天的数据。
步骤2:对于连续血糖监测信号建立二阶模型
用u(k)表示k时刻实际的血糖值,那么u(k)满足以下公式:
u(k)=2u(k-1)-u(k-2)+w1(k) (1)
其中w1(k)是均值为0、方差为λ2的高斯白噪声;令x1(k)=u(k),x2(k)=u(k-1),x(k)=[x1(k) x2(k)]T,由此得:
令
***观测方程如下:
y(k)=Cx(k)+v(k) (3)
其中y(k)是k时刻实际监测的血糖值,v(k)是均值为0、方差为σ2的高斯白噪声,C=[1 0];
步骤3:通过基于EM的方法进行模型中***噪声w(k)的协方差矩阵Q和测量噪声v(k)的协方差矩阵R的估计,具体包括以下子步骤:
(3.1)对于某个特定病人选取n个采样点的血糖信号用于EM算法的参数估计,为了下面公式表示方便,用yk表示k时刻的血糖采样值,xk表示k时刻的血糖真实值,这里认为血糖的初始值x0服从均值为μ,协方差矩阵为Σ的分布;这里选取病人第一天的连续血糖监测数据,即n=288;
(3.2)根据步骤2建立的二阶模型采用EM算法估计参数,具体步骤如下:
(3.2.1)对于μ,Q和R赋初值μ(0),Q(0)和R(0);
(3.2.2)根据式(4)-(12)计算k=1,2,...,n;
其中是k时刻的血糖滤波值,为血糖真值和滤波值之间误差的协方差矩阵,是k-1时刻对于k时刻血糖的估计值,为k时刻血糖的估计值与真实值的协方差矩阵,Kk是卡尔曼滤波增益矩阵,其中,为了计算和采用向后递归,k=n,n-1,...,1;
其中 为了下一步的计算,需要获得的值,采用向后递归的方法,k=n,n-1,...,1,
其中,
(3.2.3)根据式(13)-(15)计算Q(1)和R(1):
Q(r+1)=(W-VU-1VT)/n (14)
其中:
(3.2.4)重复步骤(3.2.2)和(3.2.3),直到估计值Q、R和对数似然函数logL稳定,其中,对数似然函数为:
对于不同的对象而言,每次进行滤波之前都要先选取该对象一天的数据用于参数估计。
步骤4:由于现有的血糖仪多采用滑动平均滤波的方法,所以对于血糖信号分别进行滑动平均滤波和卡尔曼滤波处理,以便进行对比。
(4.1)滑动平均滤波(MA)处理:
其中为第k次采样值经过滤波后的输出,y(k-i)为未经过滤波的第k-i次采样值;N为滑动平均的项数;ci为常数。其中,随着N的增大,对过去数据的“记忆”也越大,即所使用的过去数据也就越多,滤波效果也就越明显,但同时,也会造成信号有明显的延迟,使得滤波信号不能快速的跟随真实信号。
对于固定的阶数N,对于权重ci有许多不同的选法,对于血糖信号的处理,一般选取指数权重,即ci=μi,其中μ称作遗忘因子,且0<μ<1。对于滑动平均滤波来说,最大的缺点在于一旦选定了阶数和权重,对于所有的时间序列的处理都是相同的,不会因为不同的传感器和不同人造成的信噪比不同而改变,这也会造成在处理不同的CGM信号时造成次优滤波。这里设置N=5,μ=0.65。
(4.2)卡尔曼滤波(KF)处理,具体包括以下子步骤:
(4.2.1)确定对状态提前一步的预测:
P′(k)=AP(k-1)AT+Q(k-1) (22)
其中x′(k)是对于当前时刻血糖实际值的估测,是血糖前一时刻滤波后的值,而P′(k)是估测误差的协方差矩阵,P(k-1)是血糖前一时刻实际值和滤波值之间误差的协方差矩阵;
(4.2.2)确定滤波后的***状态:
K=P′(k)CT[CP′(k)CT+R]-1 (23)
P(k)=(I-KC)P′(k) (25)
其中,K是卡尔曼滤波的修正矩阵,I为二阶单位矩阵,P(k)为血糖当前时刻实际值和滤波值之间误差的协方差矩阵,为血糖当前时刻滤波后的值,即最终估计的血糖值。
(4.3)采用以下两个指标作为滤波性能的评价指标:
(4.3.1)均方根误差
均方误差RMSE的计算公式为:
其中,是k时刻血糖滤波后的值,x(k)是k时刻血糖的实际值,N是样本的总量,均方根误差RMSE越小代表滤波后的值与实际值的偏差越小,滤波效果越好;
(4.3.2)延迟时间
血糖滤波的目的是为了最终用于高/低血糖报警,及时的报警有利于病人及时采取相应的措施,从而减弱或者避免异常血糖对于患者造成的影响。然而实际情况中,经过滤波后的值与实际值或者未经过滤波的值相比总会存在一定程度的延迟,因此引入延迟时间TL作为评价延迟的指标。对于实际应用而言,时延要求在30分钟以内,当采样周期Δt为5分钟时,6个采样时刻才有意义,否则认为数据是无效的。
延迟时间TL采用以下公式计算:
其中y(k)代表k时刻实际监测的血糖值,是k时刻血糖滤波后的值,t为采样时间间隔,使得最小的t值即为所求的延迟时间。
从表1中可以看出,基于EM算法参数估计的KF要明显好于滑动平均滤波算法,对于均方误差RMSE而言,基于EM算法参数估计的KF比MA平均减少36.2%;而在延迟时间TL方面,基于EM算法参数估计的KF的性能更加优越,比MA平均减少94.9%。为了更清晰的展示两种滤波效果的好坏,图2和图3展示了对象5采用两种滤波方法所得的血糖值与真实血糖值的比较图(噪声标准差为2)。
表1针对3组(青少年组、成人组和儿童组)共30个对象的血糖采样数据分别采用滑动平均滤波和基于EM算法的卡尔曼滤波的滤波结果对比(结果用均值±标准差表示)
为了更清晰的显示两种方法滤波性能的好坏,选择图4进行展示。
Claims (2)
1.一种包含参数估计功能滤波模块的连续血糖监测设备,其特征在于,该设备包括:用于采集人体血糖信息,输出血糖可用信号的传感器;用于对传感器的输出信号进行放大处理的信号放大器;用于对信号放大器输出的模拟信号进行数字转换的单片机;用于对单片机输出的数字信号进行数据处理的滤波器,滤波器可以集成在单片机中,也可以单独使用;用于将滤波器输出的滤波后的血糖值进行数值和波形显示的显示器;用于数据存储的存储器;所述滤波器进行数据处理的过程包括以下步骤:
(1)血糖数据预处理:将以一定采样周期Δt获得的单片机输出的连续血糖监测信号组合成一维时序数据y1×l,其中,y表示检测到的血糖信号,l为采样个数,去除其中的尖峰信号;
(2)对于连续血糖监测信号建立二阶模型:用u(k)表示k时刻实际的血糖值,那么u(k)满足以下公式:
u(k)=2u(k-1)-u(k-2)+w1(k) (1)
其中w1(k)是均值为0、方差为λ2的高斯白噪声;令x1(k)=u(k),x2(k)=u(k-1),x(k)=[x1(k) x2(k)]T,由此得:
令
***观测方程如下:
y(k)=Cx(k)+v(k) (3)
其中y(k)是k时刻实际监测的血糖值,v(k)是均值为0、方差为σ2的高斯白噪声,C=[10];
(3)通过基于EM的方法进行模型中***噪声w(k)的协方差矩阵Q和测量噪声v(k)的协方差矩阵R的估计;
(4)根据步骤(3)的参数估计,采用卡尔曼滤波对连续血糖监测信号进行数据处理,具体包括以下子步骤:
(4.1)确定对状态提前一步的预测:
P′(k)=AP(k-1)AT+Q(k-1) (5)
其中x′(k)是对于当前时刻血糖实际值的估测,是血糖前一时刻滤波后的值,而P′(k)是估测误差的协方差矩阵,P(k-1)是血糖前一时刻实际值和滤波值之间误差的协方差矩阵;
(4.2)确定滤波后的***状态:
K=P′(k)CT[CP′(k)CT+R]-1 (6)
P(k)=(I-KC)P′(k) (8)
其中,K是卡尔曼滤波的修正矩阵,I为二阶单位矩阵,P(k)为血糖当前时刻实际值和滤波值之间误差的协方差矩阵,为血糖当前时刻滤波后的值,即最终估计的血糖值。
2.根据权利要求1所述一种包含参数估计功能滤波模块的连续血糖监测设备,其特征在于,所述步骤(3)具体包括以下子步骤:
(3.1)对于某个特定病人选取n个采样点的血糖信号用于EM算法的参数估计,为了下面公式表示方便,用yk表示k时刻的血糖采样值,xk表示k时刻的血糖真实值,这里认为血糖的初始值x0服从均值为μ,协方差矩阵为Σ的分布;
(3.2)根据步骤(2)建立的二阶模型采用EM算法估计参数,具体步骤如下:(3.2.1)对于μ,Q和R赋初值μ(0),Q(0)和R(0);
(3.2.2)根据式(9)-(17)计算
其中是k时刻的血糖滤波值,为血糖真值和滤波值之间误差的协方差矩阵,是k-1时刻对于k时刻血糖的估计值,为k时刻血糖的估计值与真实值的协方差矩阵,Kk是卡尔曼滤波增益矩阵;其中,为了计算和采用向后递归,k=n,n-1,…,1,
其中为了下一步的计算,需要获得的值,采用向后递归的方法,k=n,n-1,…,1,
其中,
(3.2.3)根据式(18)-(20)计算Q(1)和R(1):
Q(r)=(W-VU-1VT)/n (19)
其中:
r表示迭代次数,从1开始;
(3.2.4)重复步骤(3.2.2)和(3.2.3),直到估计值Q、R和对数似然函数logL稳定,其中,对数似然函数为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510310303.8A CN104921736B (zh) | 2015-06-08 | 2015-06-08 | 一种包含参数估计功能滤波模块的连续血糖监测设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510310303.8A CN104921736B (zh) | 2015-06-08 | 2015-06-08 | 一种包含参数估计功能滤波模块的连续血糖监测设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104921736A CN104921736A (zh) | 2015-09-23 |
CN104921736B true CN104921736B (zh) | 2017-08-04 |
Family
ID=54109421
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510310303.8A Active CN104921736B (zh) | 2015-06-08 | 2015-06-08 | 一种包含参数估计功能滤波模块的连续血糖监测设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104921736B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105403705B (zh) * | 2015-12-05 | 2017-04-12 | 浙江大学 | 一种包含血糖分类功能故障检测模块的连续血糖监测设备 |
CN105608317B (zh) * | 2015-12-18 | 2018-06-26 | 上海集成电路研发中心有限公司 | 一种基于线性***的数字滤波装置及方法 |
CN106805943B (zh) * | 2016-11-23 | 2018-06-05 | 舒糖讯息科技(深圳)有限公司 | 基于分数阶微分方程的血糖数据处理方法及装置 |
CN106859666B (zh) * | 2017-02-15 | 2018-12-04 | 舒糖讯息科技(深圳)有限公司 | 一种血糖检测装置及其检测方法 |
CN107137093B (zh) * | 2017-04-20 | 2019-06-07 | 浙江大学 | 一种包含异常血糖概率报警器的连续血糖监测设备 |
CN108805011B (zh) * | 2018-04-24 | 2022-01-11 | 长江大学 | 一种数字滤波方法及*** |
CN111603151B (zh) * | 2020-06-17 | 2023-05-16 | 深圳智领人工智能健康科技有限公司 | 一种基于时频联合分析的无创血液成分检测方法及*** |
CN111991003A (zh) * | 2020-08-12 | 2020-11-27 | 上海萌草科技有限公司 | 基于Savitzky-Golay滤波的连续血糖平滑方法、装置、设备及存储介质 |
CN112042371B (zh) * | 2020-10-13 | 2021-07-20 | 中国农业大学 | 一种玉米籽粒清选损失监测***及方法 |
CN113796858B (zh) * | 2021-11-18 | 2022-02-08 | 湖州美奇医疗器械有限公司 | 用于血糖数据监测的参比偏移校准算法*** |
CN114052688B (zh) * | 2021-12-07 | 2022-12-20 | 山东大学 | 基于单路脉搏波的血压监测装置、存储介质及电子设备 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7022072B2 (en) * | 2001-12-27 | 2006-04-04 | Medtronic Minimed, Inc. | System for monitoring physiological characteristics |
CN101138500A (zh) * | 2007-08-29 | 2008-03-12 | 中南大学 | 基于能量守恒法的无创血糖检测仪 |
JP4571220B2 (ja) * | 2009-05-13 | 2010-10-27 | ネルコー ピューリタン ベネット エルエルシー | データ信号適応平均化方法及び装置 |
CN102293654B (zh) * | 2011-06-17 | 2014-07-16 | 清华大学 | 基于代谢热-光学方法的无创血糖检测仪 |
CN103310113B (zh) * | 2013-06-24 | 2016-03-30 | 浙江大学 | 一种基于频带分离和数据建模的通用血糖预测方法 |
CN104665842A (zh) * | 2013-11-28 | 2015-06-03 | 程昊 | 一种家用血糖监护*** |
-
2015
- 2015-06-08 CN CN201510310303.8A patent/CN104921736B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104921736A (zh) | 2015-09-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104921736B (zh) | 一种包含参数估计功能滤波模块的连续血糖监测设备 | |
US11896371B2 (en) | Method and apparatus for improving lag correction during in vivo measurement of analyte concentration with analyte concentration variability and range data | |
EP1886113B1 (en) | Use of multiple data points and filtering in an analyte sensor | |
US7650244B2 (en) | Method and device for monitoring analyte concentration by determining its progression in the living body of a human or animal | |
US9936888B2 (en) | Pulse period calculation device and biosensor equipped with the same | |
Sparacino et al. | Continuous glucose monitoring time series and hypo/hyperglycemia prevention: requirements, methods, open problems | |
JP6795980B2 (ja) | 加速度の測定値を処理して、ユーザによるステップを識別する方法及び対応するコンピュータプログラム、並びに加速度の測定値を処理して、ユーザによるステップを識別する装置及びそのような装置を備えるデバイス | |
KR102095959B1 (ko) | 인공지능 딥러닝 학습을 이용한 생체측정물 농도 측정방법 | |
CN107137093B (zh) | 一种包含异常血糖概率报警器的连续血糖监测设备 | |
JP2002532694A (ja) | 測定値を分析するための方法及び装置 | |
CN111540436B (zh) | 一种自适应葡萄糖和胰岛素的浓度预测***及方法 | |
MX2021008013A (es) | Procesamiento de se?ales de sensores con calibracion basada en kalman. | |
CN105806898A (zh) | 一种用于气体传感器的气体浓度标定方法 | |
JP2004089267A (ja) | 睡眠深度推定装置及びこれを具えた布団 | |
KR101276973B1 (ko) | 맥박수 측정 방법 및 장치 | |
CN113288133A (zh) | 用于预测血糖值的方法、装置、存储介质及处理器 | |
CN106618570B (zh) | 一种基于生物介电谱的皮肤生化指标检测方法及*** | |
CN115639335A (zh) | 水质监测数据校准方法、***及计算机可读存储介质 | |
Pernice et al. | A validity and reliability study of Conditional Entropy Measures of Pulse Rate Variability | |
US20230218204A1 (en) | Personalized calibration for glucose sensing | |
CN113100735A (zh) | 一种基于静脉血流图的血栓检测装置及方法 | |
US10024709B2 (en) | Method and system to quickly determine a weight | |
CN112932424A (zh) | 一种数据获取方法及*** | |
Poliński et al. | Estimation of blood pressure parameters using ex-Gaussian model | |
US20170363565A1 (en) | Analyte Measurement |
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 | ||
GR01 | Patent grant |