一种基于正则化最小二乘回归学习的qrs波起点终点定位方法
【专利摘要】一种基于正则化最小二乘回归学习的QRS波起点终点定位方法，该方法是：基于函数在Hilbert空间中正交分解的思想结合正则化最小二乘回归学习算法，构建一系列检测方法、检测准则和阈值参数，最终检测出QRS波并定位出其起点、终点位置；首先将心电信号与高斯函数的一阶导函数做内积运算得到内积序列，选择由高斯核函数生成的再生核Hilbert空间作为正则化最小二乘回归学习算法的逼近空间，并用实对称矩阵的平方根分解的办法求解算法，得到新的内积序列。
【专利说明】
一种基于正则化最小二乘回归学习的QRS波起点终点定位 方法
技术领域
[0001] 本发明涉及动态心电图、常规心电图中QRS波起点及终点的定位测量分析，具体为 一种基于正则化最小二乘回归学习的QRS波起点终点定位方法，属于心电图（ECG)自动诊断 分析技术领域。
【背景技术】
[0002] 心电图是临床应用最广泛的检查手段之一，能够帮助诊断心律失常，帮助诊断心 肌缺血、心肌梗死、判断梗死部位，判断药物或电解质对心脏的影响，以及帮助植入人工心 脏起搏器等。另外，描记心电图操作简单，费用低廉，特别适合医疗条件不发达的地区使用。
[0003] 目前，我国专业的心电图技师数量匮乏，对于一些较为疑难的心电图，很多临床医 生都不能准确判读。因此，基于动态心电图和常规心电图自动诊断技术应运而生，其中最关 键的莫过于QRS波起点、终点位置的准确定位。用分规人工测量误差较大，即便是一个毫秒 的误差在临床上都可能造成两种截然不同的诊断结果，比如完全性右束支传导阻滞与不完 全性右束支传导阻滞只有一个毫秒的差别。在采样率足够高的情形下，利用先进的算法，借 助计算机自动分析确定QRS波起点、终点的位置就能更好的解决该问题。但是，心电图千变 万化，QRS波形态各异，这在客观上要求算法具有良好的泛化性，同时能够保证较高的处理 效率。
[0004] 心电图自动分析技术正在逐步发展，对QRS波（心拍）的定位技术已经较为成熟，但 对于QRS波起点终点位置定位的准确性亟待提高。目前，市场上主流的QRS波起点终点定位 方法大多只适用于某一类或某几类QRS波，对于一些特殊情形如T波高尖(图1.2)、基线严重 漂移（图1.3)等波形的QRS波分析定位准度较低，对于信号质量较差的波形更难做出准确的 分析判断。
【发明内容】
[0005] 本发明的目的是克服现有技术存在的不足之处，降低基线漂移、肌电干扰、ST段改 变、QRS波的形态等因素对分析定位的影响，进而提供一种具有良好泛化性和准确度的QRS 波起点、终点定位算法。本发明衍生于回归学习，具有坚实的数学基础，心拍检测尤其是QRS 波起点、终点定位较为准确，为动态心电图、常规心电图自动诊断技术提供了一种非常可靠 的定位算法。
[0006] 本发明是通过如下技术方案来完成的，借鉴Hilbert空间中函数正交分解的思想， 将原始心电信号映射到一个新的Hilbert空间中，最后在上述新空间中完成QRS波的检测及 其起点和终点的定位工作。上述映射是通过将原始心电信号与高斯函数的一阶导函数做内 积实现的。为了降低P波、T波对检测QRS波的影响，同时使得QRS波起点和终点的特征更加明 显，在检波和定位工作开始前先用正则化最小二乘回归学习算法处理内积序列。
[0007] 本发明具有良好的泛化性，对各种不同形态的QRS波都具有良好的适用性，极大提 高了心电图分析的准确率，可应用于各种动态心电图机、常规心电图机、心电监护仪等设备 的心电图自动诊断分析。
【附图说明】
[0008] 图1.1为原始心电信号（正常心电图）和未经正则化最小二乘回归学习算法处理的 内积序列。
[0009] 图1.2为原始心电信号(T波高尖的心电图）和未经正则化最小二乘回归学习算法 处理的内积序列。
[0010] 图1.3为原始心电信号(基线严重漂移的心电图）和未经正则化最小二乘回归学习 算法处理的内积序列。
[0011 ] 图2.1 (正常心电图，与图1.1对应)为未经处理的内积序列和经过正则化最小二乘 回归学习算法处理之后的内积序列。
[0012]图2.2(T波高尖，与图1.2对应)为未经处理的内积序列和经过正则化最小二乘回 归学习算法处理之后的内积序列。
[0013]图2.3(基线严重漂移，与图1.3对应)为未经处理的内积序列和经过正则化最小二 乘回归学习算法处理之后的内积序列。
[0014]图3.1(正常心电图，与图1.1对应)为心电信号QRS波起点、终点定位效果图。
[0015 ]图3.2 (Τ波高尖，与图1.2对应)为心电信号QRS波起点、终点定位效果图。
[0016]图3.3(基线严重漂移，与图1.3对应)为心电信号QRS波起点、终点定位效果图。 [0017]图4.1为基于正则化最小二乘回归学习算法的QRS波起点终点定位方法流程图。 [0018]图4.2为用遗传算法优化正则化参数lambda的流程图。
【具体实施方式】
[0019]以下结合技术方案和附图以实例详细说明本发明的【具体实施方式】。
[0020]本发明解决心电图自动化诊断系统中QRS波检测及其起点、终点定位问题。基于回 归学习与信号在Hilbert空间中的正交分解理论构建一系列检测方法、检测准则和阈值参 数，检测和识别QRS波，定位QRS波的起点、终点位置，具体实施步骤是如下所述。
[0021 ] 在以下论述中，若h表示函数，则h (X)表示函数h在点X处的函数值;若h表示数组 (整型或浮点型），则h[i]表示数组h的第i个元素，其中i为非负整数。
[0022] 1.心电信号变换
[0023]由于心电信号复杂多变，为了提高算法的泛化性，我们需要将心电信号投影到一 个维数相对较低的空间中。由于新空间维数较低，所以在新空间中研究会大大减少问题的 复杂性。
[0024]众所周知，Hilbert空间是完备化的度量空间，Hilbert空间具有标准正交基。在一 组标准正交基下，Hilbert空间中的元素具有正交分解的形式。
[0025]即若Η是Hilbert空间，feH，丨^=1是《的一组标准正交基，则f具有如下正交分解 形式：
[0027]其中，
[0029] 由连续函数组成的Hilbert空间中的元素具有如（1)所示的正交分解形式，同样信 号空间中的元素也具有如（1)所示的正交分解形式，只是(2)式需要做下相应变形。
[0030] 记所有心电信号组成的空间为EjlJ_Ee/2，12表示所有均方收敛数列组成的空间， 这是一个Hilbert空间。设^^二是12的一组标准正交基，gel2,则g具有如下正交分解形式：
[0033]当然，我们可以通过选择一组合适的标准正交基将E投影到一个较低维的空间中， 我们也可以选择一组普通的基底达到同样的效果。
[0034]本发明中采用的是由函数组成的一组基底，其中
[0036] 经投影变换后的心电图如图1.1、1.2、1.3所示，这种变换消除了T波高尖、基线漂 移等因素对定位QRS波的影响，原因在于与T波、P波相比，QRS波在这样一组基底下投影的模 更大。
[0037] 记原始心电信号为e，内积序列为S，e与S都是单精度浮点型数据，长度为HKKKLS 的计算过程如下所示。
[0039] 2.用正则化最小二乘回归学习算法处理内积序列S
[0040]设{^·,,乃为一组样本，正则化最小二乘回归学习算法形式如下
[0042]其中，Ηκ是由核函数K所生成的再生核Hilbert空间，| | · | |κ表示空间Ηκ中的范数， λ表示正则化参数。这里我们取
[0044]作为核函数，其中I I · I |2表示欧式空间中的2范数。如(5)所示的正则化最小二乘 回归学习算法具有如下形式的解析解：
[0046]其中，
[0047] a = (ai,a2,---,am)T0
[0048] a是方程
[0049] (Aml+G)a = y (6)的解，其中
[0050] y = (yi,y2,---,ym)T〇
[0051] G是一个mXm维实对称正定矩阵。矩阵G的第i行，第j列元素通过如下计算得到
[0053] 在实际操作过程中，为了提高效率，首先需要将(2)所示线性系统转化为
[0054] LLTa = y
[0055] 再求解，这里L是矩阵λπιΙ+G的平方根分解。
[0056] 在本发明中，取
[0057] m= 100,
[0059] 使用正则化最小二乘回归学习算法最关键的步骤是确定正则化参数λ的值。首先， 我们利用CV(Cross Validation)的办法确定λ，然后使用遗传算法(Genetic Algorithms) 对λ进行优化。CV是用来验证学习型算法性能的一种统计分析方法，其基本思想是在某种意 义下将原始数据进行分组，一部分作为训练集，另一部分作为验证集:先用训练集训练，再 利用验证集测试训练得到的模型，以得到的均方误差作为评价学习型算法性能的指标。
[0060] 遗传算法是1962年由美国密歇根大学Holland教授提出的模拟自然界遗传机制和 生物进化论的一种并行随机搜索最优化方法。它把自然界优胜劣汰，适者生存的生物进化 原理引入优化参数形成的编码串联群体中，按照所选择的适应度函数并通过遗传中的选 择，使适应度值好的个体被保留，适应度差的个体被淘汰，新的群体既继承了上一代的信 息，又优于上一代。这样反复循环，直至满足条件。
[0061 ]遗传算法的基本操作依次为:选择操作，交叉操作，变异操作。将对训练集进行CV 意义下的最小均方误差作为GA中的适应度函数值，则利用GA对正则化最小二乘回归学习算 法进行优化的整体算法过程如图4.2所示。定义单精度浮点型数组Si，长度为10000,将经正 则化最小二乘回归学习算法处理后的内积序列存储在Si中。
[0062]经正则化最小二乘回归学习算法处理前后的内积序列如图2.1、2.2、2.3所示。 [0063] 3.检测 QRS 波
[0064] 该过程可分解为下列几个小步骤进行：
[0065] 3.1定义单精度浮点型数组&，长度为10000。对31的每个元素都取绝对值然后按顺 序存储在3 2中。
[0066] 3.2定义单精度浮点型数组S3,长度为10000;定义整型数组N，长度为10000。对S 2做 降序排列，得到的新序列按顺序存储在&中，将&中的元素在&中所处的位置存储在N中，BP
[0067] S3[i] = S2[N[i]]，i = 0,1，···，9999。
[0068] 3 · 3构造二维整型数组Μ，维数为200 X 200。其中，
[0069] M[i][j]=|N[i]_N[j]|，i，j = 0，l，."，199。
[0070] 经过3.2的排序操作之后，所有QRS波位置处所对应的内积的模一定都排在了前 面。如果Μ的维数过大，QRS波的个数可能会多出来;反之，可能会漏掉QRS波。大量测试表明 将Μ的维数定为200 X 200是最合适的。
[0071] 这里| · |表示绝对值运算。
[0072] 3.4 重置 Μ[0][0]为 1
[0073] Μ[0][0]存储的是内积序列中模最大的元素所在的位置，这个位置一定存在一个 QRS波，但是在3.3中已经将该位置置0,为了在接下来的搜索中不至于漏掉一个QRS波，所以 要将该位置重置为0。
[0074] 3.5若紙1][」]彡200，则将第]_行的值全部置0，1 = 0，1，~，198;]_ = 1+1，...，199。
[0075] 200毫秒是心肌细胞绝对不应期的长度，也就是说在200毫秒内绝对不可能出现两 个心搏，所以我们选择200毫秒作为步长筛选QRS波。
[0076] 3.6定义整型数组A，长度为50，初始化为0。若M[ i ] [0]乒0，i = 0，1，…，199，则将Ν
[i]存储在A[i]中。定义整型变量Num，计算A[i]中所有不为0的元素的个数，并将其存储在 Num 中。
[0077] 3.7定义单精度浮点型数组S4，长度为Num，按照以下过程初始化
[0078] S4[i] = S2[A[i]]，i = 0，l，···，Num_l。
[0079] 3.8定义单精度浮点型变量Median，计算数组S4的中位数，并存储在Median中，若
[0080] S4[i ]<0 · 5 · Median，i = 0，1，…，Num_l，
[0081 ]则将A[ i ]重置为0，同时Num自减1。施行这一步操作的目的是剔除可能错检的QRS 波。
[0082] 3.9定义整型数组处8，长度为灿111，并将4[1]中所有不为0的数存储在处8中，〇^中 存储的就是所有QRS波位置。
[0083] 4.定位QRS波起点，该过程分两种情形进行。定位QRS波起点的基本思想是：
[0084]通过寻找内积序列的突变点最终定位起点的准确位置。从信号处理的角度来讲 QRS波起点前后心电信号中占主要成分的频率是不同的，也就是说QRS波起点处是信号频率 分布发生突变的点。这在内积序列上表现为QRS波起点所在位置附近对应的内积会出现一 个小的波峰。本发明就是通过定位这个小波峰的位置来定位QRS波起点位置的。
[0085]首先，定义整型数组QrsStart用于存储QRS波起点位置，长度为Num，初始化为0。
[0086]当51[(^[;[]]>〇4 = 0，1，"，1111111-1时，(^波起点定位过程如下：
[0087] 4.1定义整型变量j，并用Qrs[ i ]初始化，i = 0，1，…，Num-Ι。
[0088] 4.2从j开始向前搜索，若Si [ j ] -Si [ j-1 K0，1 < j < 9999，则j即第i个QRS波起点所 处位置，将j存储在QrsStart [ i ]中。
[0089 ]当51[(^[;[]]<〇4 = 0，1，~，1111111-1时，(^波起点定位过程如下：
[0090] 4.3定义整型变量j，初始化为Qrs [ i ]，i = 0，1，…，Num-Ι。
[0091] 4 · 4定义整型变量temp，初始化为0。
[0092] 4 · 5从j开始向前搜索，若Si [ j hSi [ j-Ι ]彡0，1 < j <9999，则将j存储在temp中。
[0093] 4 · 6从temp 开始向前搜索，若Si [temp ]_Si[ temp-1 ] <0，K temp <9999，则 temp 即 第i个QRS波起点所处位置，将temp存储在QrsStart[i]中。
[0094] 5.定位QRS波终点，定位QRS波终点的思想与起点相同。该过程分两种情形进行。
[0095]首先定义整型数组QrsFinal用于存储QRS波终点位置，长度为Num，初始化为0。
[0096] 当51[(^[;[]]>〇4 = 0，1，"，1111111-1时，(^波终点定位过程如下：
[0097] 5.1定义整型变量j，并用Qrs[i]初始化，i = 0，l，，"，Num-l。
[0098] 5.2定义整型变量temp，初始化为0。
[0099] 5.3从」开始向后搜索，若51[」]-51口+1]彡0,0彡衫9998，则将」存储在仏111?中。
[0100] 5 · 4从 temp开始向后搜索，Si[temp]_Si[temp+l ]多0，CXtemp^i9998，则 temp 即第 i 个QRS波终点所处位置，将其存储在QrsFinal [ i ]中。
[0101] 当51[(^[;[]]<〇4 = 0，1，~，11111]1-1时，(^波终点定位过程如下：
[0102] 5 · 5定义整型变量j，初始化为Qrs[ i ]，i = 0，1，…，Num-1。
[0103] 5 · 6从j开始向后搜索，若Si [ j ] -Si [ j+1 ]彡0，0彡j < 9998，则j即第i个QRS波终点所 处位置，将其存储在QrsFinal [ i ]中。
[0104] QRS波起点终点定位效果如图3.1、3.2、3.3所示。
[0105] 以上各步骤中所述的具体内容(包括公式、变换方法、技术路线等），均可通过计算 机代码来实现，本领域技术人员应该能够理解，对于任何一个特定的方法和步骤采用不同 的方式来实现所述功能，或者通过替换其中的一些参数以及推演出类似的公式方法和判定 规则，只要是不脱离本专利构思的前提下，均应当视为属于本专利的保护范围。
【主权项】
1. 一种基于正则化最小二乘回归学习的QRS波起点终点定位方法，该方法是:基于函数 在Hilbert空间中正交分解的思想结合正则化最小二乘回归学习算法，构建一系列检测方 法、检测准则和阈值参数，最终检测出QRS波并定位出其起点、终点位置;其特征在于： 首先将心电信号与高斯函数的一阶导函数做内积运算得到内积序列，选择由高斯核函 数生成的再生核Hilbert空间作为正则化最小二乘回归学习算法的逼近空间，并用实对称 矩阵的平方根分解的办法求解算法，得到新的内积序列； 然后用正则化最小二乘回归学习算法处理内积序列，具体是通过对内积序列依照绝对 值大小进行排序，用特定的步长筛选内积序列，实现对QRS波位置的定位； 最后通过分析内积序列检测出QRS波并确定其起点、终点位置；即通过在QRS波位置一 定邻域内寻找内积序列的突变点确定QRS起点、终点位置。2. 根据权利要求1所述的基于正则化最小二乘回归学习的QRS波起点终点定位方法，其 特征在于所述的QRS波起点终点定位方法包括如下实施步骤： 在以下步骤中，若h表示函数，则h (X)表示函数h在点X处的函数值;若h表示数组(整型 或浮点型），则h[i]表示数组h的第i个元素，其中i为非负整数； 步骤1:心电信号变换;将心电信号投影到一个维数相对较低的空间中；若Η是Hilbert 空间，f e H，是!1的一组标准正交基，则f具有如下正交分解形式： / = Σ?,·η ⑴ /-1 其中： ?,=〈/- *Λ〉"' =厂 /(命(·ν)办 12、 由连续函数组成的Hilbert空间中的元素具有如式（1)所示的正交分解形式，同样信号 空间中的元素也具有如式(1)所示的正交分解形式，只是式(2)需要做下相应变形： 记所有心电信号组成的空间为E，则五g; f，I2表示所有均方收敛数列组成的空间，这是 一个Hilbert空间；设丨￡〇^是12的一组标准正交基，gel2,则g具有如下正交分解形式： s=i^iei ⑶ ?~1 β； =(s^,)r=Y,s{.i)-eX.i) H) 采用由函数组成的一组基底，其中：经投影变换后的心电图消除了 T波高尖、基线漂移等因素对定位QRS波的影响，原因在 于与T波、P波相比，QRS波在这样一组基底下投影的模更大； 记原始心电信号为e，内积序列为S，e与S都是单精度浮点型数据，长度为lOOOOd的计 算过程如下所示。 ,Γ,Τ |Σ^(/>"[./]^· = 4·^··-9995 冲 4: τ Ο, / - 0,1, 2, 3,9996,9997,9998, 9999 1) 用正则化最小二乘回归学习算法处理内积序列S， 设为一组样本，正则化最小二乘回归学习算法形式如下(5) 其中，Ηκ是由核函数Κ所生成的再生核Hilbert空间，| | · | |κ表示空间Ηκ中的范数，λ表 示正则化参数，这里我们取作为核函数，其中I I · I |2表示欧式空间中的2范数;如式(5)所示的正则化最小二乘回 归学习算法具有如下形式的解析解： χ：(·ν)^Σβλ'(-ν，-ν) > i-l. 其中， α= (αι,α2, ··· ,am)T； a是方程 (Aml+G)a = y (6) 的解，其中， y=(yi，y2，."，ym)T， G是一个mXm维实对称正定矩阵，矩阵G的第i行，第j列元素通过如下计算得到：在实际操作过程中，为了提高效率，首先需要将式(2)所示线性系统转化为 LLTa = y 再求解，这里L是矩阵λπιΙ+G的平方根分解， 取 m= 100， {(W, )}: = ;(，-外·))}:〉 使用正则化最小二乘回归学习算法最关键的步骤是确定正则化参数λ的值，首先，利用 CV(Cross Validation)的办法确定λ，然后使用遗传算法(Genetic Algorithms)对λ进行优 化； 遗传算法的基本操作依次为:选择操作，交叉操作，变异操作，将对训练集进行CV意义 下的最小均方误差作为GA中的适应度函数值，则利用GA对正则化最小二乘回归学习算法进 行优化的整体算法，定义单精度浮点型数组Si，长度为10000,将经正则化最小二乘回归学 习算法处理后的内积序列存储在Si中； 2) 检测QRS波，该过程可分解为下列几个小步骤进行： 第一，定义单精度浮点型数组&，长度为10000;对31的每个元素都取绝对值然后按顺序 存储在S2中； 第二，定义单精度浮点型数组S3，长度为10000;定义整型数组N，长度为10000;对32做降 序排列，得到的新序列按顺序存储在S3中，将S3中的元素在S2中所处的位置存储在N中，BP S3[i]=S2[N[i]]，i = 0，l，...，9999。 第三，构造二维整型数组Μ，维数为200 X 200，其中， M[i][j]=|N[i]-N[j]|，i，j = 0，l，…，199。 经过第二步的排序操作之后，所有QRS波位置处所对应的内积的模一定都排在了前面， 如果Μ的维数过大，QRS波的个数可能会多出来;反之，可能会漏掉QRS波;大量测试表明将Μ 的维数定为200 X 200是最合适的。 这里I · I表示绝对值运算； 第四，重置Μ[0][0]为1， Μ[0][0]存储的是内积序列中模最大的元素所在的位置，这个位置一定存在一个QRS 波，但是在第三步骤中已经将该位置置0,为了在接下来的搜索中不至于漏掉一个QRS波，所 以要将该位置重置为1; 第五，若M[i][j]彡200,则将第j行的值全部置0，i = 0，l，…，198;j = i+l，."，199;200 毫秒是心肌细胞绝对不应期的长度，也就是说在200毫秒内绝对不可能出现两个心搏，所以 我们选择200毫秒作为步长筛选QRS波； 第六，定义整型数组A，长度为50，初始化为0;若M[ i ] [0]辛0，i = 0，1，…，199，则将N[ i ] 存储在A[ i ]中；定义整型变量Num，计算A[ i ]中所有不为0的元素的个数，并将其存储在Num 中； 第七，定义单精度浮点型数组S4,长度为Num，按照以下过程初始化： S4[i] = S2[A[i] ]，i = 0,1，…，Num_l 〇 第八，定义单精度浮点型变量Median，计算数组S4的中位数，并存储在Median中，若 S4[i]<0.5 · Median, i = 0，1，…，Num_l， 则将A[i ]重置为0，同时Num自减1。施行这一步操作的目的是剔除误检的QRS波； 第九，定义整型数组Qrs，长度为Num，并将A[ i ]中所有不为0的数存储在Qrs中，Qrs中存 储的就是所有QRS波位置； 3)定位QRS波起点，该过程分两种情形进行； 首先，定义整型数组QrsStart用于存储QRS波起点位置，长度为Num，初始化为0; 当Si[Qrs[i]] >0, i = 0,1，…，num-l时，QRS波起点定位过程如下： 第一，定义整型变量j，并用Qrs[i]初始化，i = 0，l，···，Num-l; 第二，从j开始向前搜索，若Si [ j ] -Si [ j-1 ]彡0，1彡j彡9999，则j即第i个QRS波起点所处 位置，将j存储在QrsStart [ i ]中； 当Si[Qrs[i]]<0, i = 0,1，…，num-1时，QRS波起点定位过程如下： 第三，定义整型变量j，初始化为Qrs[i ]，i = 0，1，…，Num-l; 第四，定义整型变量temp，初始化为0; 第五，从j开始向前搜索，若Si [ j ] -Si [ j-Ι ]》0，1 < j <9999，则将j存储在temp中； 第六，从 temp 开始向前搜索，若 Si [temp]-Si [temp-1 ]<0, K temp< 9999,则temp 即第i 个QRS波起点所处位置，将temp存储在QrsStart [ i ]中； 4)定位QRS波终点，定位QRS波终点的思想与起点相同；该过程分两种情形进行： 首先定义整型数组QrsFinal用于存储QRS波终点位置，长度为Num，初始化为0; 当Si[Qrs[i]] >0, i = 0,1，…，num-l时，QRS波终点定位过程如下： 第一，定义整型变量j，并用Qrs[i]初始化，i = 0，l，···，Num-l; 第二，定义整型变量temp，初始化为0; 第三，从j开始向后搜索，若Si [ j ] -Si [ j+Ι ]彡0，0彡j彡9998，则将j存储在temp中； 第四，从temp开始向后搜索，51|^11^1]-51|^11^1+1]多0,0^^11^1<9998，则七611^1即第;[个 QRS波终点所处位置，将其存储在QrsFinal [ i ]中； 当31[(^[;[]]<〇4 = 0，1，~，1111111-1时，(^波终点定位过程如下： 第五，定义整型变量j，初始化为Qrs[i]，i = 0，l，···，Num-l; 第六，从j开始向后搜索，若Si [ j ] -Si [ j+1 ]彡0，0彡j彡9998，则j即第i个QRS波终点所处 位置，将其存储在QrsFinal [ i ]中。
【文档编号】A61B5/0472GK106063704SQ201610369281
【公开日】2016年11月2日
【申请日】2016年5月27日 公开号201610369281.7, CN 106063704 A, CN 106063704A, CN 201610369281, CN-A-106063704, CN106063704 A, CN106063704A, CN201610369281, CN201610369281.7
【发明人】朱玉奎, 符灵建
【申请人】浙江好络维医疗技术有限公司