专利名称：一种胃排空参数和胃蠕动参数同步检测方法和装置的制作方法
技术领域：
本发明涉及生物医学信息提取和处理领域，特别涉及胃排空参数与 胃蠕动参数的同步检测方法和装置。
背景技术：
胃排空是一个复杂的过程，对固体和液体食物排空方式也不尽相同， 容易受到多种因素，如神经调节、激素、进食量、酸碱度、食物的化学 构成、食物的热值、情绪等的影响，其不同区域功能的控制和协调通过 固有的和外源-的神经及激素调节实现。胃排空的过程伴随着 一 系列的胃 蠕动，胃蠕动对食物的消化和吸收至关重要，胃蠕动的情况直接影响胃
排空的过程。胃排空和胃蠕动是胃动力过程的两个环节，两者相互依赖； 胃排空和胃蠕动信息的同步提取有助于增强对胃动力功能的认识。
目前采用放射性核素法来同步测量胃蠕动和胃排空过程，但是因为 核素法存在有害性、操作困难等问题，而不太被人们接受。发明专利CN 96236620和CN 00202383等公开了 一种胃动力检测装置，包括恒流源电 路、信号调理电路、模/数转换器、计算机等实现胃运动信号的检测， 这类发明利用胃排空时胃的形态、容积及其内容物组成情况的改变较大， 电特性变化非常明显，表征胃内容物和容积状态的信号可以通过其相对 应的胃阻抗变化获得的原理，检测出阻抗信号，但是其没有将排空信号 与蠕动信号相结合进行测量分析，使得胃动力信息未得到充分利用，未 能为胃动力学研究提供可靠的参考依据。
发明内容
为解决以上问题，本发明提出 一种胃排空参数和胃蠕动参数同步检 测方法和装置。一种胃排空参数和胃蠕动参数同步检测方法，其特征在于将拾取的 生物阻抗信号做高精度模数转换， 一 方面通过高精度模数转换数据生成 拟合排空曲线，计算胃半排空时间，另一方面对高精度模数转换数据做 小波分解和小波重构，计算胃蠕动信号特定参数。
所述拾取生物阻抗信号，由恒流激励源产生电流注入人体，通过电 极拾取相应电压信号，通过放大和滤波等电路调理信号，完成生物阻抗 信号的拾取。
所述对生物阻抗信号做高精度模数转换，是指对拾取的生物阻抗信
号转换为数字信号，模数转换的采样频率为5Hz，高精度指采用24比特 以上量化位数。
所述生成拟合胃排空曲线包括以下步骤
(1) 确定基础阻抗R0,可以是空腹检测时间内的阻抗平均值或者采 样起点到饮水前时刻阻抗平均值；
(2) 确定排空时间起点，可以饮水后达到的阻抗峰值Rm点为排空 时间起点；
(3) 确定排空终点T,以饮水后达到的阻抗值降为基础阻抗到基础 阻抗1.5倍的任一时刻为排空终点，该阻值为排空终点阻抗值 Ri;
(4) 计算拟合排空曲线参数k,根据经验，拟合排空曲线表示为 jK0 = Mxe-to，根据排空终点阻抗值R1、排空终点时间T和阻 抗峰值Rm,计算出拟合排空曲线参数k;
其中，t为排空过间，oasr, W)是阻抗值，M是幅度，取值为饮 水后达到的阻抗峰值
所述计算胃半排空时间，包括以下步骤
(1)寻找快速排出相切线A，以阻抗峰值Rm与基线阻抗RO的差为单位，确定排空信号下降1/4单位和1/2单位的时间点，由
峰值点结合以上两点确定一段圆弧，快速排出相切线A就是在 峰值点的圓弧切线； (2)寻找排空终端切线，由基础阻抗基线向后延展，与阻抗信号相 交的交点为基点，连接阻抗峰值点和该点，形成的切线称为排 空终值切线B;
(3 )确定半排空时间，由快速排出相切线A和排空终值切线B得到 其角平分线C，在角平分线C上，阻抗峰值与饮水前的基础阻 抗值之差的 一 半时所对应的时间确定为半排空时间。
所述小波分解和小波重构的小波基可根据胃蠕动信号特点构造，也
可以采用通用DB4-DB7小波基；
对高精度模数转换数据做小波分解与重构后获得胃蠕动信号，提取 胃蠕动节律和幅度信号并计算胃蠕动信号特定参数
所述的节律，每分钟胃蠕动次数，通过对蠕动信号周期计数获取。 所述的幅度，胃蠕动信号的强度，单位为电阻欧姆，从提取的蠕动
信号的纵坐标数值直接获得。
所述胃蠕动信号特定参数，包括功率谱密度、变异系数IC等； 所述的功率谱密度为对胃蠕动信号做快速傅立叶变换后计算平均功
率，以一分钟为周期计算功率谱，获得时间上延续的动态谱；
<formula>formula see original document page 9</formula>其中，x(n)表示胃蠕动信号，X(k)表示快速傅立叶变换的离散频语，
N表示一分钟内信号点数，P(K)表示功率谱密度，conj表示复共轭运算；
所述变异系数/C计算式为<formula>formula see original document page 9</formula>其中，D/^，画。/为所有主频在2-4CPM范围的动态谱^跌的支配频率
均值，SDpre。由下式获得
其中，Z)尸(^是第k个动态i普段的支配频率，即主频，DF,是所有 动态谱分段的支配频率均值，其中，S是用于计算的动态谱段数目。
本发明另 一 目的在于提供一种胃排空参数与胃蠕动参数的同步检测 装置，包括生物阻抗信号拾取单元，其特征在于，还包括高精度模数转 换单元、胃排空曲线拟合单元、胃半排空时间计算单元和胃蠕动信号特 定参数计算单元；
所述生物阻抗信号拾取单元，包括信号调理电路、放大电路、滤波 电路、经过信号调理电路、放大电路、滤波电路，信号输入到信号采集
单元；
所述高精度模数转换电路单元，用于信号经高精度模数转换电路后
获取生物阻抗信号；模数转换的采样频率为5Hz，高精度指采用24比特 以上量化位数。
所述胃排空曲线生成单元，用于从生物阻抗信号中提取胃排空信息， 生成胃排空曲线解析式；
所述生成胃排空曲线解析式，包括以下步骤
(1) 确定基础阻抗RO，可以是空腹检测时间内的阻抗平均值或者采 样起点到饮水前时刻阻抗值；本实施例以空腹检测10分钟内 的阻抗平均值为基础阻抗值；
(2) 确定排空时间起点，饮用纯净水后胃内阻抗快速增加，本实施 例以饮水后达到的阻抗峰值点Rm为排空时间起点；
(3) 确定排空终点T,以饮水后达到的阻抗值降为基础阻抗到基础阻抗1.5倍的任一时刻为排空终点，该阻值为排空终点阻抗值
Rl ，本实施例以饮水后达到的阻抗值降为基础阻抗时刻为排空 终点；
作为另一优选实施方式，以饮水后达到的阻抗值降为基础阻抗1.5 倍的时刻为排空终点；
(4)计算拟合排空曲线参数k,根据经验，拟合排空曲线表示为 WhMxe-to,根据排空终点阻抗值R1、排空终点时间T和阻 抗峰值，计算出拟合排空曲线参数k; 其中，t为排空过间，oasr, XO是阻抗值，M是幅度，取值为饮 水后达到的阻抗峰值。
所述胃半排空时间测量单元，用于计算胃半排空时间； 所述计算胃半排空时间，包括以下步骤
(1) 寻找快速排出相切线A，以阻抗峰值Rm与基线阻抗RO的差 为单位，确定排空信号下降1/4单位和1/2单位的时间点，由 峰值点结合以上两点确定一段圓弧，快速排出相切线A就是在 峰值点的圆弧切线；
(2) 寻找排空终端切线，由基础阻抗基线向后延展，与阻抗信号相 交的交点为基点，连接阻抗峰值点和该点，形成的切线称为排 空终值切线B;
(3 )确定半排空时间，由快速排出相切线A和排空终值切线B得到 其角平分线C,在角平分线C上，阻抗峰值与饮水前的基础阻 抗值之差的 一半时所对应的时间确定为半排空时间。 所述胃蠕动信号特定参数计算单元，包括小波分解单元、小波重构 单元、节律提取单元、幅度提取单元、功率谱计算单元和变异系数计算 单元等，分别提取蠕动信号中包含的节律、幅度信息和计算胃蠕动信号特定参数。
所述胃蠕动信号特定参数，包括功率谱密度、变异系数IC等； 所述的功率谱密度为对胃蠕动信号做快速傅立叶变换后计算平均功
率，以一分钟为周期计算功率谱，获得时间上延续的动态谱；
X(K) = FFT{x(n)}， K = 0,1".."N-1 P(K) = 2*( X(K).*conj(X(K)))/N 其中，x(n)表示胃蠕动信号，X(k)表示快速傅立叶变换的离散频谱， N表示一分钟内信号点数，P(K)表示功率谱密度，conj表示复共轭运算；
所述变异系数/C计算式为/C、，
其中，D7^wm。/为所有主频在2-4CPM范围的动态谱段的支配频率 均值，SD^q由下式获得
其中，DF(W是第k个动态谱段的支配频率，即主频，Z)F。g是所有 动态谱分段的支配频率均值，其中，S是用于计算的动态i普段数目。
本发明有效实现了胃排空参数与胃蠕动参数的同步检测，为胃动力 学研究提供了一种更便捷的检测方式，并为胃动力学研究提供可靠的参 考依据。
图1为本发明一种胃排空参数和胃蠕动参数同步检测方法流程图 图2为本发明计算半排空时间示意图
图3为本发明一种胃排空参数和胃蠕动参数同步检测装置框图
具体实施例方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明的 一种胃排空参数和胃蠕动参数同步检测方法及 装置作进一步说明。应当理解，此处所描述的具体实施例^f又用以解释本 发明，并不用于限定本发明。
测试者空腹进入系统，进行测量信号，启动计时，空腹检测10分钟
后，饮用纯净水400ml，继续测量信号，饮水后测量时间30分钟。
一种胃排空参数和胃蠕动参数同步检测方法，其特征在于将拾取的
生物阻抗信号做高精度模数转换， 一 方面通过高精度模数转换数据生成
拟合排空曲线，并计算胃半排空时间，另一方面对高精度模数转换数据
做小波分解和小波重构，计算胃蠕动信号特定参数
所述拾取生物阻抗信号，由恒流激励源产生电流注入人体，通过电
极拾取相应电压信号，通过放大和滤波等电路调理信号，完成生物阻抗
信号的拾取。
所述恒流激励源、信号调理电路、信号放大电路、信号解调电路为 本领域公知技术，本实施不再详述。但使用时可以选用数字合成恒流源， 数字解调等新的芯片方式。
所述对生物阻抗信号做高精度模数转换，为对拾取的生物阻抗信号 转换为数字信号，本实施例模数转换的采样频率为5Hz,采用24比特量 化位数。
所述生成拟合胃排空曲线包括以下步骤
(1 )确定基础阻抗RO，可以是空腹检测时间内的阻抗平均值或者采 样起点到饮水前时刻阻抗值；本实施例以空腹^f全测10分钟内 的阻抗平均值为基础阻抗值；
(2) 确定排空时间起点，饮用纯净水后胃内阻抗快速增加，本实施 例以饮水后达到的阻抗峰值点Rm为排空时间起点；
(3) 确定排空终点T，以饮水后达到的阻抗值降为基础阻抗到基础阻抗1.5倍的任一时刻为排空终点，该阻值为排空终点阻抗值
Rl,本实施例以饮水后达到的阻抗值降为基础阻抗时刻为排空
终点；
作为另一优选实施方式，以饮水后达到的阻抗值降为基础阻抗1.5
倍的时刻为排空终点；
(4)计算拟合排空曲线参数k,根据经验，拟合排空曲线表示为 少(0-Mx^'，根据排空终点阻抗值R1、排空终点时间T和阻 抗峰值，计算出拟合排空曲线参数k;
其中，t为排空过间，ogsr, y(/)是阻抗值，M是幅度，取值为饮 水后达到的阻抗峰值。
由于完整的胃排空时间(Gastric Emptying Time,简称GET)较长， 而且排空时间终点往往难以准确确定，临床通常以半排空时间(GET / 2) 进行胃排空功能的量度。
所述计算胃半排空时间，包括以下步骤
(1) 寻找快速排出相切线A，以阻抗峰值Rm与基线阻抗R0的差 为单位，确定排空信号下降1/4单位和1/2单位的时间点，由 峰值点结合以上两点确定一段圓弧，快速排出相切线A就是在 峰值点的圓弧切线；
(2) 寻找排空终端切线，由基础阻抗基线向后延展，与阻抗信号相 交的交点为基点，连接阻抗峰值点和该点，形成的切线称为排 空终值切线B;
(3 )确定半排空时间，由快速排出相切线A和排空终值切线B得到 其角平分线C,在角平分线C上，阻抗峰值与饮水前的基础阻 抗值之差的 一 半时所对应的时间确定为半排空时间。
所述小波分解和小波重构的小波基可根据胃蠕动信号特点构造，也可以采用通用DB4-DB7小波基。
本实施例，小波分解的小波基系数如下 h
:=0.2303778133088964; h[l]:=0.7148465705529154; h[2]:=0.6308807679398597; h[3]:=-0.0279837694168599; h[4〗:=-0.1870348 U 7190931; h[5]:=0.0308413818355607; h[6]:=0.0328830116668852; h[7] :=-0.0105974017850890;
小波重构小波基系数由以上分解系数变化得到，为分解系数倒序相
反数，即
h[O]:=0.0105974017850890;
h[l]:=-0.0328830116668852
h[2]=-0.0308413818355607
h[3]:=0.1870348117190931;
h[4]=0.0279837694168599.
h[5]=-0.6308807679398597
h[6]=-0.7148465705529154
h[7]=-0.2303778133088964
对高精度模数转换数据做小波分解与重构后获得胃蠕动信号，提取 胃蠕动节律和幅度信号并计算胃蠕动信号特定参数
所述的节律，每分钟胃蠕动次数，通过对蠕动信号周期计数获取。 所述的幅度，胃蠕动信号的强度，单位为电阻欧姆，从提取的蠕动
信号的纵坐标数值直接获得；
所述胃蠕动信号特定参数，包括功率谱密度、变异系数IC等； 所述的功率谱密度为对胃蠕动信号做快速傅立叶变换后计算平均功率，以一分钟为周期计算功率谱，获得时间上延续的动态谱；
X(K) = FFT{x(n)}， K = 0,1，....，N-1 P(K) = 2*( X(K).*conj(X(K)))/N 其中，x(n)表示胃蠕动信号，X(k)表示快速傅立叶变换的离散频谱， N表示一分钟内信号点数，P(K)表示功率谱密度，conj表示复共轭运算；
所述变异系数/C计算式为/C =，"
其中，D/^g,^为所有主频在2-4CPM范围的动态谱段的支配频率
均值，SD^q由下式获得
幼一
其中，Z)^T"是第k个动态谱段的支配频率，即主频，D尸^是所有
动态谱分段的支配频率均值，其中，s是用于计算的动态谱段数目。
一种胃排空曲线和胃运动信号同步检测的装置，包括生物阻抗信号 拾取单元，其特征在于，还包括高精度模数转换单元、胃排空曲线拟合 单元、胃半排空时间计算单元和胃蠕动信号特定参数计算单元。
所述生物阻抗信号拾取单元，包括信号调理电路、放大电路、滤波 电路、经过信号调理电路、放大电路、滤波电路，信号输入到信号采集
单元；
所述高精度模数转换电路单元，用于对生物阻抗信号做高精度模数
转换，本实施例模数转换的采样频率为5Hz，采用24比特量化位数。 所述胃排空曲线生成单元，用于从生物阻抗信号中提取生成胃排空
信息，生成胃排空曲线解析式；
所述生成胃排空曲线解析式，包括以下步骤(1 )确定基础阻抗RO，可以是空腹检测时间内的阻抗平均值或者采
样起点到饮水前时刻阻抗值；本实施例以空腹检测10分钟内 的阻抗平均值为基础阻抗值；
(2) 确定排空时间起点，饮用纯净水后胃内阻抗快速增加，本实施 例以^:水后达到的阻抗峰值点Rm为排空时间起点；
(3) 确定排空终点T，以饮水后达到的阻抗值降为基础阻抗到基础 阻抗1.5倍的任一时刻为排空终点，该阻值为排空终点阻抗值 Rl,本实施例以饮水后达到的阻抗值降为基础阻抗时刻为排空 终点；
作为另一优选实施方式，以饮水后达到的阻抗值降为基础阻抗1.5 倍的时刻为排空终点；
(4) 计算拟合排空曲线参数k，根据经验，拟合排空曲线表示为 _Kr) = Mxe-to,根据排空终点阻抗值R1、排空终点时间T和阻 抗峰值，计算出拟合排空曲线参数k;
其中，t为排空过间，oasr,少(f)是阻抗值，M是幅度，取值为饮 水后达到的阻抗峰值。
所述胃半排空时间测量单元，用于计算胃半排空时间； 所述计算胃半排空时间，具体包括以下步骤
(1 )寻找快速排出相切线A,以阻抗峰值Rm与基线阻抗R0的差 为单位，确定排空信号下降1/4单位和1/2单位的时间点，由 峰值点结合以上两点确定一段圓弧，快速排出相切线A就是在 峰值点的圆弧切线；
(2)寻找排空终端切线，由基础阻抗基线向后延展，与阻抗信号相 交的交点为基点，连接阻抗峰值点和该点，形成的切线称为排 空终值切线B;(3 )确定半排空时间，由快速排出相切线A和排空终值切线B得到 其角平分线C,在角平分线C上，阻抗峰值与饮水前的基础阻 抗值之差的 一半时所对应的时间确定为半排空时间。 所述胃蠕动信号相关参数计算单元，包括小波分解单元、小波重构 单元、节律提取单元、幅度提取单元、功率谱计算单元和变异系数计算 单元等，分别用于提取蠕动信号中包含的节律、幅度信息和计算胃蠕动 信号特定参数。
所述小波分解和小波重构的小波基可根据胃蠕动信号特点构造，也 可以采用通用DB4-DB7小波基。
本实施例，小波分解的小波基系数如下
h[O] :=0.2303778133088964;
h[l]:=0.7148465705529154;
h[2] :=0.6308807679398597;
h[3] :=-0.0279837694168599;
h[4]:=-0.1870348117190931;
h[5]:=0.0308413818355607;
h[6] :=0.0328830116668852;
h[7]:=-0.0105974017850890;
小波重构小波基系数由以上分解系数变化得到，为分解系数倒序相
反数，即
h[O]:=0.0105974017850890;
h[l]'=-0.0328830116668852;
h[2]=-0.0308413818355607;
h[3]'=0.1870348117190931;
h[4]=0.0279837694168599;
h[5]=-0.6308807679398597;
h[6]=-0.7148465705529154;h[7]:=-0.2303778133088964;
对高精度模数转换数据做小波分解与重构后获得胃蠕动信号，提取 胃蠕动节律和幅度信号并计算胃蠕动信号特定参数
所述的节律，每分钟胃蠕动次数，通过对蠕动信号周期计数获取。 所述的幅度，胃蠕动信号的强度，单位为电阻欧姆，从提取的蠕动
信号的纵坐标数值直接获得；
所述胃蠕动信号特定参数，包括功率谱密度、变异系数IC等； 所述的功率谱密度为对胃蠕动信号做快速傅立叶变换后计算平均功
率，以一分钟为周期计算功率语，获得时间上延续的动态谱；
X(K) = FFT{x(n)}， K-0，1，…"N-1 P(K) = 2*( X(K).*conj(X(K)))/N 其中，x(n)表示胃蠕动信号，X(k)表示快速傅立叶变换的离散频谱，
N表示一分钟内信号点数，P(K)表示功率谱密度，conj表示复共轭运算；
所述变异系数/C计算式为/C=幼
其中，Z)&，。^。/为所有主频在2-4CPM范围的动态谱段的支配频率 均值，SD^由下式获得
其中，D尸f"是第k个动态i普段的支配频率，即主频，D/^g是所有 动态谱分段的支配频率均值，其中，S是用于计算的动态谱段数目。
对于本发明各个实施例中所阐述的方法和装置，凡在本发明的精神 和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明 的保护范围之内。
权利要求
1、一种胃排空参数和胃蠕动参数同步检测方法，其特征在于将拾取的生物阻抗信号做高精度模数转换，一方面通过高精度模数转换数据生成拟合排空曲线，计算胃半排空时间，另一方面对高精度模数转换数据做小波分解和小波重构，计算胃蠕动信号特定参数。
2、 如权利要求1所述的一种胃排空参数和胃蠕动参数同步检测方法，其特征在于，所述对生物阻抗信号做高精度模数转换，是指将拾取的生物阻抗信号转换为数字信号，模数转换的采样频率为5Hz, 24 比特以上量化位数。
3、 如权利要求1所述的一种胃排空参数和胃蠕动参数同步检测方 法，其特征在于，所述生成拟合胃排空曲线包括以下步骤(1) 确定基础阻抗R0，可以是空腹^r测时间内的阻抗平均值或 者采样起点到饮水前时刻阻抗平均值；(2) 确定排空时间起点，可以饮水后达到的阻抗峰值Rm点为排 空时间起点；(3) 确定排空终点T,以饮水后达到的阻抗值降为基础阻抗到基 础阻抗1.5倍的任一时刻为排空终点，该阻值为排空终点阻 抗值Rl;(4) 计算拟合排空曲线参数k，根据经验，拟合排空曲线表示为 _K0 = Mxe-to,根据排空终点阻抗值R1、排空终点时间T和阻 抗峰值Rm,计算出拟合排空曲线参数k;其中，t为排空过间，oasr, XO是阻抗值，M是幅度，取值为 饮水后达到的阻抗峰值。
4、 如权利要求1所述的一种胃排空参数和胃蠕动参数同步检测方 法，其特征在于，所述计算胃半排空时间，包括以下步骤(1 )寻找快速排出相切线A,以阻抗峰值Rm与基线阻抗RO的差为单位，确定排空信号下降1/4单位和1/2单位的时间点， 由峰值点结合以上两点确定一段圆弧，快速排出相切线A就 是在峰值点的圓弧切线；(2) 寻找排空终端切线，由基础阻抗基线向后延展，与阻抗信号 相交的交点为基点，连接阻抗峰值点和该点，形成的切线称 为排空终值切线B;(3) 确定半排空时间，由快速排出相切线A和排空终值切线B 得到其角平分线C,在角平分线C上，阻抗峰值与饮水前的 基础阻抗值之差的 一半时所对应的时间确定为半排空时间。
5、 如权利要求1所述的一种胃排空参数和胃蠕动参数同步检测方 法，其特征在于，所述小波分解和小波重构的小波基可根据胃蠕动信 号特点构造，也可以采用通用DB4-DB7小波基。
6、 如权利要求1所迷的一种胃排空参数和胃蠕动参数同步检测方 法，其特征在于，所述胃蠕动信号特定参数包括功率谱密度、变异系 数IC;所述的功率i普密度为对胃蠕动信号做快速傅立叶变换后计算平 均功率，以一分钟为周期计算功率谱，获得时间上延续的动态谱；X(K) = FFT{x(n)}， K = 0,1,,...， N-l P(K) = 2*( X(K).*conj(X(K)))/N 其中，x(n)表示胃蠕动信号，X(k)表示快速傅立叶变换的离散频 谱，N表示一分钟内胃蠕动信号点数，P(K)表示功率谱密度，conj表 示复共轭；所述变异系数/C计算式为/C=其中，/)尸勿 ,为所有主频在2-4CPM范围的动态谱段的支配频率均值，SD^由下式获得幼—二其中，DiYig是第k个动态i普段的支配频率，即主频，DF^是所 有动态谱分段的支配频率均值，其中，S是用于计算的动态谱段数目。
7、 一种胃排空曲线和胃蠕动信号同步检测装置，包括生物阻抗信号 拾取单元，其特征在于，还包括高精度模数转换单元、胃排空曲线拟 合单元、胃半排空时间计算单元和胃蠕动信号特定参数计算单元；所述生物阻抗信号拾取单元，包括信号调理电路、放大电路、滤 波电路、经过信号调理电路、放大电路、滤波电路，信号输入到信号 采集单元；所述高精度模数转换电路单元，用于信号经高精度模数转换电路 后获取生物阻抗信号；模数转换的采样频率为5Hz， 24比特以上量化 位数；所述胃排空曲线生成单元，用于从生物阻抗信号中提取胃排空信 息，生成胃排空曲线解析式；所述胃半排空时间测量单元，用于计算胃半排空时间； 所述胃蠕动信号相关参数计算单元，包括小波分解单元、小波重 构单元、功率谱计算单元和用于变异系数计算单元，分别用于提取蠕 动信号中包含的节律、幅度信息和计算胃蠕动信号特定参数。
8、 如权利要求7所述的一种胃排空曲线和胃蠕动信号同步检测装 置，其特征在于，所述生成胃排空曲线解析式，包括以下步骤(1)确定基础阻抗R0,可以是空腹检测时间内的阻抗平均值或 者采样起点到饮水前时刻阻抗平均值；(2) 确定排空时间起点，可以饮水后达到的阻抗峰值Rm点为排 空时间起点；(3) 确定排空终点T，以饮水后达到的阻抗值降为基础阻抗到基 础阻抗1.5倍的任一时刻为排空终点，该阻值为排空终点阻 抗值Rl;(4) 计算拟合排空曲线参数k，根据经验，拟合排空曲线表示为 少( )-Mxe-to,根据排空终点阻抗值R1、排空终点时间T和阻 抗峰值Rm，计算出拟合排空曲线参数k;其中，t为排空过间，owsr， y(r)是阻抗值，M是幅度，取值为 饮水后达到的阻抗峰值。
9、 如权利要求7所述的一种胃排空曲线和胃蠕动信号同步检测装 置，其特征在于，所述计算胃半排空时间，包括以下步骤(1) 寻找快速排出相切线A,以阻抗峰值Rm与基线阻抗RO的 差为单位，确定排空信号下降1/4单位和1/2单位的时间点， 由峰值点结合以上两点确定一段圆弧，快速排出相切线A就 是在峰值点的圆弧切线；(2) 寻找排空终端切线，由基础阻抗基线向后延展，与阻抗信号 相交的交点为基点，连接阻抗峰值点和该点，形成的切线称 为排空终值切线B;(3) 确定半排空时间，由快速排出相切线A和排空终值切线B 得到其角平分线C，在角平分线C上，阻抗峰值与饮水前的 基础阻抗值之差的 一半时所对应的时间确定为半排空时间。
10、 如权利要求7所述的一种胃排空曲线和胃蠕动信号同步检测装 置，其特征在于，所述胃蠕动信号特定参数，包括功率谱密度、变异 系数IC;所述的功率谱信息为对胃蠕动信号做快速傅立叶变换后的数据，以一分钟为周期计算功率谱，获得时间上延续的动态谱；X(K) = FFT{x(n)} ， K-0，1，.…，N-1 P(K) = 2*( X(K).*conj(X(K)))/N 其中，x(n)表示胃蠕动信号，X(k)表示快速傅立叶变换的离散频 谱，N表示一分钟内信号点数，P(K)表示功率谱密度，conj表示复共 轭运算；所述变异系数/C计算式为/C=叫啤其中，D/^g輕则/为所有主频在2-4CPM范围的动态谱段的支配频 率均值，SD—由下式获得幼细<formula>formula see original document page 6</formula>其中，Z^Y"是第k个动态谱段的支配频率，即主频，DF。^是所 有动态谱分段的支配频率均值，其中，S是用于计算的动态谱段数目。
全文摘要
本发明涉及生物医学信息提取和处理领域，特别涉及一种胃排空参数与胃蠕动参数的同步检测方法，其特征在于将拾取的生物阻抗信号做高精度模数转换，一方面通过高精度模数转换数据生成拟合排空曲线，计算胃半排空时间，另一方面对高精度模数转换数据做小波分解和小波重构，计算胃蠕动信号特定参数；还涉及一种胃排空参数与胃蠕动参数的同步检测装置，包括生物阻抗信号拾取单元、高精度模数转换单元、胃蠕动信号特定参数计算单元、胃排空曲线拟合单元、胃半排空时间计算单元；本发明实现了胃排空参数与胃蠕动参数的同步检测，为胃动力学研究提供便捷的检测方式及可靠的参考依据。
文档编号A61B5/00GK101579231SQ20091010411
公开日2009年11月18日 申请日期2009年6月17日 优先权日2009年6月17日
发明者任超世, 追 徐, 李章勇, 伟 王, 娜 胡, 赵德春, 魏进民 申请人:重庆邮电大学