CN110840411B - 一种麻醉深度的测量装置、存储介质及电子设备 - Google Patents

一种麻醉深度的测量装置、存储介质及电子设备 Download PDF

Info

Publication number
CN110840411B
CN110840411B CN201911245011.5A CN201911245011A CN110840411B CN 110840411 B CN110840411 B CN 110840411B CN 201911245011 A CN201911245011 A CN 201911245011A CN 110840411 B CN110840411 B CN 110840411B
Authority
CN
China
Prior art keywords
signal
sub
electroencephalogram
features
predicted
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
Application number
CN201911245011.5A
Other languages
English (en)
Other versions
CN110840411A (zh
Inventor
韩如泉
温鹏
熊飞
任冠清
李兴
周赤宜
戴仁泉
王筱毅
李明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Delikai Medical Electronics Co ltd
Original Assignee
Shenzhen Delica Medical Equipment Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen Delica Medical Equipment Co ltd filed Critical Shenzhen Delica Medical Equipment Co ltd
Priority to CN201911245011.5A priority Critical patent/CN110840411B/zh
Publication of CN110840411A publication Critical patent/CN110840411A/zh
Priority to PCT/CN2020/106066 priority patent/WO2021109601A1/zh
Application granted granted Critical
Publication of CN110840411B publication Critical patent/CN110840411B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4821Determining level or depth of anaesthesia
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/369Electroencephalography [EEG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7246Details of waveform analysis using correlation, e.g. template matching or determination of similarity
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7275Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor

Abstract

本发明公开了一种麻醉深度的测量装置、存储介质及电子设备,其包括:获取脑电信号,并提取所述脑电信号的信号特征;获取所述脑电信号对应的预测脑电信号,并提取所述预测脑电信号的预测信号特征;根据所述信号特征以及预测信号特性生成所述脑电信号对应的模拟脑电信号,并根据所述模拟脑电信号确定麻醉深度值。本发明通过将通过网络脑模型生成的预测脑电信号与实际测量得到的脑电信号相结合得到模拟脑电信号,再根据该模拟脑电信号计算麻醉深度,这样可以通过预测脑电信号对实际测量脑电信号进行验证,可以提高测量得到的麻醉深度的准确性,提高麻醉深度的可靠性,并且提高麻醉深度测量的动态实时跟踪能力和抗干扰能力。

Description

一种麻醉深度的测量装置、存储介质及电子设备
技术领域
本发明涉及超声技术领域,特别涉及一种麻醉深度的测量装置、存储介质及电子设备。
背景技术
在进行重大手术前,需要对患者进行全身麻醉(俗称全麻),而全麻具有极高风险,例如,麻醉过深会引起麻醉并发症,甚至危及患者生命;麻醉过浅易发生“术中知晓”,产生痛苦和恐惧及精神后遗症;因此如何精确地估算麻醉深度,使得医生可以根据麻醉深度确定麻醉用量,以提高麻醉安全性。
目前基于脑电信号展开的麻醉深度监测的方法主要包括双频指数、听觉诱发电位指数、脑功能状态指数,熵指数,复杂度和小波分析法等。例如,美国Aspect公司(now partof Covidien)主推出的BIS监护仪,其实采用双谱分析方便,并通过0-100的指数来反映麻醉意识深度。然而,上述方法中普遍采用的是患者的真实脑电信号,而真实脑电信号在采集过程中会存在干扰信号等,造成真实脑电信号异常。此时,医生需要根据血压、心率、呼吸频率、肌松程度等间接指标来判断病人的麻醉状态,这需要依赖医生的业务水平,并且不同医生的判断结果也会存在不同。
发明内容
鉴于现有技术的不足,本发明旨在提供一种麻醉深度的测量装置、存储介质及电子设备。
本发明所采用的技术方案如下:
一种麻醉深度的测量方法,其包括:
获取脑电信号,并提取所述脑电信号的信号特征;
获取所述脑电信号对应的预测脑电信号,并提取所述预测脑电信号的预测信号特征,其中,所述预测脑电信号为通过预设矩阵型脑网络生成;
根据所述信号特征以及预测信号特性生成所述脑电信号对应的模拟脑电信号,并根据所述模拟脑电信号确定麻醉深度值。
所述麻醉深度的测量方法,其中,所述预测脑电信号的生成过程具体包括:
以预先建立的神经元群模块为元素建立矩阵型脑网络;
根据预设联通参数驱动所述矩阵型脑网络产生脑电信号,其中,所述预设联通参数用于控制所述矩阵型脑网络中各神经元群模块状态。
所述麻醉深度的测量方法,其中,所述神经元群模块包括激励神经细胞网络、锥体细胞网络以及抑制神经细胞网络;所述激励神经细胞网络的激励信号、抑制神经细胞网络的抑制信号、以及锥体细胞网络的外部激励信号形成脑电波,并将所述脑电波输出至外部的神经元群模块,并且所述脑电波分别反馈至该神经元群模块的激励神经细胞网络以及抑制神经细胞网络。
所述麻醉深度的测量方法,其中,所述获取脑电信号,并提取所述脑电信号的信号特征具体包括:
获取脑电信号,并对所述脑电信号进行去噪处理;
提取去噪处理后的脑电信号的信号特征。
所述麻醉深度的测量方法,其中,所述根据所述信号特征以及预测信号特性生成所述脑电信号对应的模拟脑电信号具体包括:
将所述信号特征中各子信号特征分别与所述预测信号特征中对应的预测信号特征进行匹配;
根据所述匹配结果确定脑电信号对应的模拟信号特征,并根据所述模拟信号特征生成所述模拟脑电信号。
所述麻醉深度的测量方法,其中,所述根据所述匹配结果确定脑电信号对应的模拟信号特征具体包括:
对于匹配结果为相同的子信号特征,将该子信号特征作为该子信号特征对应的子模拟信号特征;
对于匹配结果为不相同的子信号特征,将该子信号特征进行处理,并将处理后的子信号特征作为该子信号特征对应的子模拟信号特征;
根据确定所有子模拟信号特征生成所述脑电信号对应的模拟信号特征。
所述麻醉深度的测量方法,其中,所述对于匹配结果为不相同的子信号特征,将该子信号特征进行处理,并将处理后的子信号特征作为该子信号特征对应的子模拟信号特征具体包括:
对于匹配结果为不相同的子信号特征,判断是否通过线性插值方式将该子信号特征匹配至其对应的子预测信号特征;
若是,则将该子信号特征进行线性插值处理以得到该子信号特征对应的子模拟信号特征。
所述麻醉深度的测量方法,其中,所述对于匹配结果为不相同的子信号特征,将该子信号特征进行处理,并将处理后的子信号特征作为该子信号特征对应的子模拟信号特征包括:
若否,则获取所述脑电信号的信噪比;
当所述信噪比满足预设条件时,采用移动窗平移算法确定该子信号特征对应的子模拟信号特征;
当所述信噪比未满足预设条件时,将该子信号特征对应的子预测信号特征作为该子信号特征对应的子模拟信号特征;
当未检测到信噪比时,将该子信号特征与该子信号特征对应的子预测信号特征进行加权处理,以得到该子信号特征对应的子模拟信号特征。
一种计算机可读存储介质,所述计算机可读存储介质存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现如上任一所述的麻醉深度的测量方法中的步骤。
一种电子设备,其包括:处理器、存储器及通信总线;所述存储器上存储有可被所述处理器执行的计算机可读程序;
所述通信总线实现处理器和存储器之间的连接通信;
所述处理器执行所述计算机可读程序时实现如上任一所述的麻醉深度的测量方法中的步骤。
有益效果:与现有技术相比,本发明提供了一种麻醉深度的测量装置、存储介质及电子设备,所述方法包括获取脑电信号,并提取所述脑电信号的信号特征;获取所述脑电信号对应的预测脑电信号,并提取所述预测脑电信号的预测信号特征;根据所述信号特征以及预测信号特性生成所述脑电信号对应的模拟脑电信号,并根据所述模拟脑电信号确定麻醉深度值。本发明通过将通过网络脑模型生成的预测脑电信号与实际测量得到的脑电信号相结合得到模拟脑电信号,再根据该模拟脑电信号计算麻醉深度,这样可以通过预测脑电信号对实际测量脑电信号进行验证,可以提高测量得到的麻醉深度的准确性,提高麻醉深度的可靠性,并且提高麻醉深度测量的动态实时跟踪能力和抗干扰能力。
附图说明
图1为本发明提供的麻醉深度的测量方法的流程图。
图2为本发明提供的麻醉深度的测量方法中神经元群模块的示意图。
图3为本发明提供的麻醉深度的测量方法中人脑模型的一个角度的示意图。
图4为本发明提供的麻醉深度的测量方法中人脑模型的另一个角度的示意图。
图5为本发明提供的电子设备的结构原理图。
具体实施方式
本发明提供一种麻醉深度的测量装置、存储介质及电子设备,为使本发明的目的、技术方案及效果更加清楚、明确,以下参照附图并举实施例对本发明进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本技术领域技术人员可以理解,除非特意声明,这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是,本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件,但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。应该理解,当我们称元件被“连接”或“耦接”到另一元件时,它可以直接连接或耦接到其他元件,或者也可以存在中间元件。此外,这里使用的“连接”或“耦接”可以包括无线连接或无线耦接。这里使用的措辞“和/或”包括一个或更多个相关联的列出项的全部或任一单元和全部组合。
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语),具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语,应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样被特定定义,否则不会用理想化或过于正式的含义来解释。
下面结合附图,通过对实施例的描述,对发明内容作进一步说明。
本实施提供了一种麻醉深度的测量方法,如图1和图2所示,所述方法包括:
S10、获取脑电信号,并提取所述脑电信号的信号特征。
具体地,所述脑电信号可以通过EEG采集设备采集得到,也可以通过电极片获取得到,还可以通过传感器采集得到,并且通过传感器采集到所述脑电信号可以是模拟脑电信号,并且在采集到所述模拟脑电信号可以将模拟脑电信号转换为数字脑电信号,并将所述数字脑电信号作为所述脑电信号,这样可以便于脑电信号的存储。当然,值得说明的,EEG采集设备可以为本领域技术人员所公知的各种采集脑电信号的设备。此外,所述脑电信号所述脑电信号可以是患者在麻醉全程中各段脑电信号,例如,每段脑电信号的信号长度为10秒,并通过滑动窗口技术,通过步长控制每秒脑电信号数据的更新。
所述信号特征可以包括时域特征、频域特征和非频域特征中的至少一种。在本实施例的一个实现方式中,所述信号特征包括时域特征或频域特征,即在获取到脑电信号后,提取该脑电信号的时域特征或频域特征。
进一步,在本实施例的一个实现方式中,所述获取脑电信号,并提取所述脑电信号的信号特征具体包括:
S11、获取脑电信号,并对所述脑电信号进行去噪处理;
S12、提取去噪处理后的脑电信号的信号特征。
具体地,所述去噪处理为对生理干扰信号和非生理干扰信号进行识别和处理。对脑电信号进行去噪处理,可以提高脑电信号的指令,由此,去除干扰信号能够更加准确的麻醉深度值。在本实施例的一个实现方式中,所述生理干扰信号可以包括眼动干扰信号和肌电干扰信号,所述非生理干扰信号可以包括信号幅度异常、信号斜率异常以及电刀干扰等。由此,为了去除生理干扰信号和非生理干扰信号,可以对所述脑电信号进行数字滤波,其中,所述数字滤波可以包括低通滤波,高通滤波以及特征波删除等方式。
S20、获取所述脑电信息对应的预测脑电信息,并提取所述预测脑电信号的预测信号特征,其中,所述预测脑电信息为通过预设矩阵型脑网络生成。
具体地,所述预设矩阵型脑网络为预先建立,并通过矩阵型脑网络可以生成预测脑电信号。所述预测脑电信号可以以根据预设麻醉方案确定脑电信号对应的麻醉数据,并根据所述麻醉数据实时改变模型参数动态地获取脑电信号对应的预测脑电信号。所述预测脑电信号可以根据预设麻醉方案中麻醉数据,通过改变模型参数提前离线地获取并存储备用麻醉过程中的预测脑电信号。在本实施例的一个可能实现方式中,所述预测脑电信号的获取过程为:在确定麻醉方案后,根据麻醉方案确定各时刻对应的模型参数,之后通过各时刻对应的模型参数按照时间顺序输入至矩阵型脑网,可以得到麻醉过程中的预测脑电信号,在得到预测脑电信号后,可以存储该麻醉过程对应的预测脑电信号。而所述脑电信号对应的预测脑电信号,可以是根据所述脑电信号对应的时间段在预设矩阵型脑网络生成的预测脑电信号中选取该时间段对应的预测脑电信号,以得到所述脑电信号对应的预测脑电信号,这样可以减少预测脑电信号的获取时间,提高麻醉深度获取的实时性。此外,所述麻醉数据可以包括麻醉时间以及进入麻醉状态的速度等。
进一步,在本实施例的一个实现方式中,所述预测脑电信息的生成过程具体包括:
M10、以预先建立的神经元群模块为元素建立矩阵型脑网络;
M20、根据预设联通参数驱动所述矩阵型脑网络产生脑电信号,其中,所述预设联通参数用于控制所述矩阵型脑网络中各神经元群模块状态。
具体地,所述神经元群模块为预选建立的,每个神经元群模块表示一脑皮层区域,所述矩阵脑网络包括多个神经元群模块,并且各神经元群模块对应的脑皮层区域构成脑皮层,即脑皮层可以划分为多个脑皮层区域,每个脑皮层对于一个神经元群模块,并且该神经元群模块包含于形成所述矩阵型脑网络的多个神经元群模块中。同时所述多个脑皮层区域中的每个脑皮层区域负责不同的功能,例如,枕叶区负责视觉功能;顶叶区负责触觉和空间;颞叶区负责听觉和综合感知;额叶运动皮层负责协调肢体运动;前额皮层区负责理解、记忆和判断等。在本实施例的一个可能实现方式中,所述神经元群模块包括76个,即将左右脑皮层共划分为76个区域,相应的,所述矩阵型脑网络具有76个元素,目前76个元素是我们能够达到的最好结果,当然更多元素会使模型更加准确。理论上越多的JR模型能模拟出来的脑电信号越接近真实脑电信号,但会带来系统效率低下。目前选择76这个值只是在保证模拟的脑电信号有效的前提下在系统运行性能与结果复杂度之间平衡的结果。
进一步,脑电信号EEG是大量神经元活动的反映,特别是在麻醉过程中,脑电信号EEG变化具有很强的非线性特征。其中,麻醉过程中脑电信号EEG是根据丘脑和脑皮层间的大量神经元会产生振荡耦合形成,并且且随着麻醉深度的变化,耦合强度会有明显不同,以使得脑电信号EEG变化具有很强的非线性特征,此外,在丘脑和脑皮层间的振荡耦合过程中,丘脑的功能很像是一个知觉“开关”或整体脑活动的显示器,并且丘脑和脑皮层神经网络的规则性活动是麻醉和知觉的关键。例如,在开始进入睡眠状态和维持睡眠的过程中,有一个正向的抑制过程,该抑制过程是由前脑底部和丘脑下部经由GaBa+离子传递到觉醒核心,同时觉醒核心同时向失觉核心发出以反馈过程,这说明当觉醒核心受到抑制时,该反馈正过程加强了失觉核心的神经活动,从而引发一个类似于单稳触发电路一样的有知觉和无知觉状态,并且人脑始终处于一个有知觉或无知觉的状态。由此,基于预设神经元群模块构建一个分部平滑的过程模型,即矩阵型脑网络,所述矩阵型脑网络是以预设建立的神经元群模块为元素,并按照各神经元群模块之间的耦合强度建立矩阵型脑网络,同时所述矩阵型脑网络可以把微观的脑细胞活动,同介观的脑皮层区活动和宏观的脑电图信号号结合起来,以使得通过所述矩阵型脑网络可以模拟脑电信号。
进一步,在本实施例的一个实现方式中,所述神经元群模块为JANSEN-RIT(J-R)模块,所述JANSEN-RIT模块包括激励神经细胞网络、锥体细胞网络以及抑制神经细胞网络;所述激励神经细胞网络的激励信号、抑制神经细胞网络的抑制信号、以及锥体细胞网络的外部激励信号形成脑电波,并将所述脑电信号输出至外部的神经元群模块。如图2所示,所述激励神经细胞网络、锥体细胞网络以及抑制神经细胞网络均可以表示为神经环路。所述激励神经细胞网络包括第一脉冲支路以及第一反馈支路,第一脉冲支路的电压信号与第一反馈支路的电压信号相加后,通过电压-脉冲编码器将相加后的电压信号转换为脉冲信号后传输至锥体细胞网络;所述抑制神经细胞网络包括第二脉冲支路以及第二反馈支路,第二脉冲支路的电压信号与第二反馈支路的电压信号相加后,通过电压-脉冲编码器将相加后的电压信号转换为脉冲信号后传输至锥体细胞网络;锥体细胞网络包括局部激励支路、局部抑制支路以及外部激励支路,局部激励支路接收激励神经细胞网络传输的脉冲信号,局部抑制支路接收抑制神经网络传输的脉冲信号,外部激励支路接收输入脉冲,并且所述局部激励支路的、局部抑制支路以及外部激励支路汇合后形成脑电波信号,所述脑电波信号通过电压-脉冲编码器后输出至其他神经元群模块,通过分别反馈至该JANSEN-RIT模块的激励神经细胞网络以及抑制神经细胞网络。此外,如图2所示,在所述Si,i=1,2,3为电压-脉冲编码器,用于将细胞内的电压信号转换成轴突丘上的脉冲信号;
Figure GDA0003390699740000091
以及
Figure GDA0003390699740000092
均为是脉冲-电压译码器,用于把收到的脉冲信号进行加权卷积后转换成突触后的电压,其中,其中,1T、2T以及3T均为外部输出信号的神经元群模块的编号;
Figure GDA0003390699740000093
为加法器,用于把突触后的树突电压进行整合转换成细胞内的电压信号;C1T,C13,C2T,C23,C31,C32以及C3T均为神经网络群之间平均突触数;m1T(t),m1(t),m2T(t),m2(t),m3(t)以及m3T(t)均为脉冲信号,v1T(t),v13(t),v1(t),v2T(t),v23(t),v2(t),v3T(t),v31(t),v32(t)以及v3(t)均为电压信号,其中,v3(t)为可被检测到的脑电信号。
进一步,在本实施例的一个实现方式中,所述以预先建立的神经元群模块为元素建议具有空间特性的矩阵型脑网络具体包括:
S11、获取用户的脑部数据,并根据所述脑部数据生成结构型脑网络;
S12、获取该结构型脑网络中神经网络节点的耦合强度系数,并根据获取到的耦合强度系数形成矩阵序列,以得到矩阵型脑网络。
具体地,所述脑部数据为根据解剖学获取得到,根据所述脑部数据可以建议一个结构脑网络,例如,所述结构脑网络包括76个区域。在建立所述结构脑网络后,对于结构脑网络的每个区域,获取该区域中各神经网络节点之间的耦合强度系数,并根据获取到所有耦合强度系数确定该区域的矩阵,这样结构型脑网络的动态演变就会产生一系列脑网络矩阵,而在不同时间点的动态脑网络矩阵形成矩阵型脑网络。此外,所述脑部数据包括脑皮层区域和各脑皮层区域的关联性,其中,各脑皮层区域的关联强度可以用0,1和-1表示,也可以是归一化后从-1到+1的任何值,而各种脑活动则可以用一组描述大规模脑神经网络的动态时域方程表示,这些方程可以根据研究的目的进行选择,例如,我们可以选择如下的方程组对常规的脑电信号进行模拟:
Figure GDA0003390699740000094
其中,
Figure GDA0003390699740000095
代表节点i在时间t时基于局部动态函数f(xi(t))的平均电场,wij是连接节点i和j的各向同质电导矩阵,g是全局耦合函数,(t-Δtij)代表时间延迟。
此外,上述动态时域方程中的参数可以由已知生理参数值,包括不限于模型中各种脑组织的导电率等来确定。
进一步,所述预设联通参数用于表示各J-R模块的连通状态,其中,所述联通参数包括未联通、弱联通、正常联通以及强联通,例如,分别采用0,1,2,和3表示。此外,所述J-R模块是由如图2所示的三个神经网络环路组成,并且能够根据预设参数而产生脑电波,其中,所述预设参数可以为脑组织导电系数等,所述脑组织导电系数可以包括头皮,头盖骨,脑组织液,白色组织,灰色组织以及血管等。
进一步,每一个神经网路环路都是用其关键状态变量来表示,例如,平均膜电位,平均激活率和它们之间的相互转换(脉冲-波电位及波电位-脉冲)等,这样可以根据震荡频率和锥体的平均PSP基准线确定各J-R模块的连通状态,即根据预先选定用于控制状态转移的联通参数,可以使得各个脑皮区之间的相互作用将会产生诸如多稳、同步以及直至协同等状态,从而产生相应的脑电波信号。
进一步,在本实施例的一个实现方式中,所述根据预设联通参数驱动所述矩阵型脑网络产生脑电信号具体包括:
M21、获取丘脑以及脑皮层的神经活动;
M22、根据所述连通参数驱动所述矩阵型脑网络,通过所述矩阵型脑网络确定所述神经活动对应的电位信息,并根据所述电位信息得到脑电信号。
具体地,所述丘脑以及脑皮层的神经活动为已知的,可以通过CT或者磁共振等途径获取到,在获取到丘脑和脑皮层的神经活动后,根据丘脑和脑皮层的神经活动确定各神经元群模块的连通参数,并通过所述连通参数驱动所述矩阵型脑网络,以使得所述通过所述矩阵型脑网络确定所述神经活动对应的电位信息,并根据所述电位信息得到脑电信号。
进一步,在本实施例的一个实现方式中,所述通过所述矩阵型脑网络确定所述神经活动对应的电位信息,并根据所述电位信息得到脑电信号;
对丘脑以及脑皮层对应的脑部图像进行识别以区分脑部组分,并将各脑部组分与其对应的电质参数相结合,以得到人脑模型;
根据所述人脑模型确定脑部头皮的各电极点单元的电位,并根据确定到的所有电位生成所述脑电信号。
具体地,所述脑部图像可以是通过核磁共振获取到,并且所述脑部图像包括多张连续核磁共振图像,以便于可以获取到脑部各部位的图像。在获取到脑部图像后,对所述脑部图像进行识别以确定脑部组分,其中,所述脑部组分包括盖骨,脑液,脑灰色物质,脑白色物质等几十种人头组织。在识别到脑部组分后,可以将脑部组分与已知电质参数,例如,头皮,头盖骨,脑组织液,白色组织,灰色组织以及血管等,相结合而建立以数字人头模型,即得到人脑模型,如图3和4所示。此外,在建立人脑模型后,将该数字化的人脑模型抽象成几百万个小的单元,每一个单元代表人脑的不同部位和组织,并被赋予相应的电质参数,这样每一个单元不仅可以代表不同的人头组织,也可以代表同样组织的不同电质参数,之后确定数字化的人脑模型中各个小单元的点位,然后从中找出头皮上对应电极点的单元的点位,以得到脑电信号。
S30、根据所述信号特征以及预测信号特性生成所述脑电信号对应的模拟脑电信号,并根据所述模拟脑电信号确定麻醉深度值。
具体地,所述根据所述信号特征以及预测信号特性生成所述脑电信号对应的模拟脑电信号指的是将所述信号特征与预测信号特征进行匹配,以将所述脑电信号进行划分,分别划分为匹配信号集、线性插值信号集以及完全不匹配信号。其中,所述匹配信号集中的脑电信号的信号特征与其对应的预测信号特征相匹配;所述线性插值信号集中的脑电信号的信号特征可以通过线性插值的方式匹配至其对应的预测信号特征;所述完全不匹配信号中的脑电信号的信号特征与其对应的预测信号特征不匹配,并且无法通过线性插值的方式匹配至其对应的预测信号特征。此外,为了快速计算麻醉深度,匹配信号集、线性插值信号集以及完全不匹配信号中的脑电信号可以配置不同的麻醉估算标识,例如,匹配信号集中脑电信号配置有自动标识,线性插值信号集中脑电信号配置有半自动标识,完全不匹配信号中的脑电信号在根据预测脑电信号进行修正后,配置有半自动标识。这样在对信号特征与预测信号特征匹配后,可以根据麻醉估算标识确定麻醉深度。
进一步,在本实施例的一个实现方式中,所述根据所述信号特征以及预测信号特性生成所述脑电信号对应的模拟脑电信号具体包括:
将所述信号特征中各子信号特征分别与所述预测信号特征中对应的预测信号特征进行匹配;
根据所述匹配结果确定脑电信号对应的模拟信号特征,并根据所述模拟信号特征生成所述模拟脑电信号。
具体地,所述匹配结果包括匹配、通过线性插值后匹配以及不匹配且通过线性插值也无法匹配。对于不同的匹配结果的脑电信号的处理方式不同,特征是对于不匹配且通过线性插值也无法匹配的脑电信号,可以通过预测脑电信号对所述脑电信号进行修正,并根据修正后的脑血信号生成模拟脑血信号,这样可以提高模拟脑电信号的准确性,从而提高麻醉深度测量的准确性,
进一步,在本实施例的一个实现方式中,所述根据所述匹配结果确定脑电信号对应的模拟信号特征具体包括:
对于匹配结果为相同的子信号特征,将该子信号特征作为该子信号特征对应的子模拟信号特征;
对于匹配结果为不相同的子信号特征,将该子信号特征进行处理,并将处理后的子信号特征作为该子信号特征对应的子模拟信号特征;
根据确定所有子模拟信号特征生成所述脑电信号对应的模拟信号特征。
具体地,所述匹配结果为相同,说明信号特征与预测信号特征一样,由此,可以确定该信号特征对应的脑血信号为准确脑电信号,从而可以将该子信号特征作为该子信号特征对应的子模拟信号特征。而对于匹配结果为不相同的子信号特征,则需要对该子信号特征进行处理,以使得所述子信号特征与其对应的预测信号特征匹配,而当无法使得子信号特征与其对应的预测信号特征匹配是,根据所述预测信号特征对所述该子信号特征进行修正。相应的,在本实施例的一个实现方式中,所述对于匹配结果为不相同的子信号特征,将该子信号特征进行处理,并将处理后的子信号特征作为该子信号特征对应的子模拟信号特征具体包括:
对于匹配结果为不相同的子信号特征,判断是否通过线性插值方式将该子信号特征匹配至其对应的子预测信号特征;
若是,则将该子信号特征进行线性插值处理以得到该子信号特征对应的子模拟信号特征;
若否,则获取所述脑电信号的信噪比;
当所述信噪比满足预设条件时,采用移动窗平移算法确定该子信号特征对应的子模拟信号特征;
当所述信噪比未满足预设条件时,将该子信号特征对应的子预测信号特征作为该子信号特征对应的子模拟信号特征;
当未检测到信噪比时,将该子信号特征与该子信号特征对应的子预测信号特征进行加权处理,以得到该子信号特征对应的子模拟信号特征。
具体地,所述对非吻合部分将用线性插值的方法进行数据整合后产生一个半自动的标识,对于不匹配的但又无法进行插值整合的部分,分以下三种情况进行处理:1)如果实测脑电信噪比在预设范围内,应用移动窗平滑算法将实测脑电信号的信号特征转化为子模拟信号特征,2)如果实测脑电信噪比在预设范围外,将预测脑电信号的预测信号特征转化为子模拟信号特征,3)如果得不到实测脑电信噪比,则用信号特征A和预测信号特征B两者的加权值(αA+(1-α)B),以期得到子模拟信号特征。其中,所述预测信噪比对应的预设范围可以随选取的麻醉特征而定,例如,在10%到50%之间。在此,在匹配是一个过程,由于麻醉是一个受控制的渐进的连续过程,中间不可以有间断,从而在不匹配时,可以进行线性插值以使得信号特征匹配至预测信号特征。
进一步,在确定模拟脑电信号后,可以根据模拟脑电信号确定麻醉深度,其中,所述麻醉深度可以经过一个类似百分比的归一化处理给出一个0到100的深度指数,例如,所述指数可以如下图表所示,反映病人的麻醉深度。
麻醉深度指数 麻醉状态 生理/物理特征
100 清醒 各项反应正常
80 轻度麻醉 血压心律升高,对强刺激有反应
60 全麻醉 血压稳定,无知觉无反应
40 深度麻醉 血压高不稳,心律不齐
20 爆发抑制 间歇性脑电活动
0 无脑电活动 无脑电活动
基于上述麻醉深度的测量方法,本实施例提供了一种计算机可读存储介质,所述计算机可读存储介质存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现如上述实施例所述的麻醉深度的测量方法中的步骤。
基于上述麻醉深度的测量方法,本发明还提供了一种电子设备,如图5所示,其包括至少一个处理器(processor)20;显示屏21;以及存储器(memory)22,还可以包括通信接口(Communications Interface)23和总线24。其中,处理器20、显示屏21、存储器22和通信接口23可以通过总线24完成相互间的通信。显示屏21设置为显示初始设置模式中预设的用户引导界面。通信接口23可以传输信息。处理器20可以调用存储器22中的逻辑指令,以执行上述实施例中的方法。
此外,上述的存储器22中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。
存储器22作为一种计算机可读存储介质,可设置为存储软件程序、计算机可执行程序,如本公开实施例中的方法对应的程序指令或模块。处理器20通过运行存储在存储器22中的软件程序、指令或模块,从而执行功能应用以及数据处理,即实现上述实施例中的方法。
存储器22可包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序;存储数据区可存储根据电子设备的使用所创建的数据等。此外,存储器22可以包括高速随机存取存储器,还可以包括非易失性存储器。例如,U盘、移动硬盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(Random Access Memory,RAM)、磁碟或者光盘等多种可以存储程序代码的介质,也可以是暂态存储介质。
此外,上述存储介质以及移动终端中的多条指令处理器加载并执行的具体过程在上述方法中已经详细说明,在这里就不再一一陈述。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (6)

1.一种麻醉深度的测量装置,其特征在于,其包括:
采集模块,用于获取脑电信号,并提取所述脑电信号的信号特征;
神经元群模块,用于获取所述脑电信号对应的预测脑电信号,并提取所述预测脑电信号的预测信号特征,其中,所述预测脑电信号为通过预设矩阵型脑网络生成;
确定模块,用于根据所述信号特征以及所述预测信号特征生成所述脑电信号对应的模拟脑电信号,并根据所述模拟脑电信号确定麻醉深度值;
所述根据所述信号特征以及所述预测信号特征生成所述脑电信号对应的模拟脑电信号具体包括:
将所述信号特征中各子信号特征分别与所述预测信号特征中对应的预测信号特征进行匹配;
根据匹配结果确定脑电信号对应的模拟信号特征,并根据所述模拟信号特征生成所述模拟脑电信号;
所述根据所述匹配结果确定脑电信号对应的模拟信号特征具体包括:
对于匹配结果为相同的子信号特征,将该子信号特征作为该子信号特征对应的子模拟信号特征;
对于匹配结果为不相同的子信号特征,将该子信号特征进行处理,并将处理后的子信号特征作为该子信号特征对应的子模拟信号特征;
根据确定所有子模拟信号特征生成所述脑电信号对应的模拟信号特征;
所述对于匹配结果为不相同的子信号特征,将该子信号特征进行处理,并将处理后的子信号特征作为该子信号特征对应的子模拟信号特征具体包括:
对于匹配结果为不相同的子信号特征,判断是否通过线性插值方式将该子信号特征匹配至其对应的子预测信号特征;
若是,则将该子信号特征进行线性插值处理以得到该子信号特征对应的子模拟信号特征;
所述对于匹配结果为不相同的子信号特征,将该子信号特征进行处理,并将处理后的子信号特征作为该子信号特征对应的子模拟信号特征包括:
若否,则获取所述脑电信号的信噪比;
当所述信噪比满足预设条件时,采用移动窗平移算法确定该子信号特征对应的子模拟信号特征;
当所述信噪比未满足预设条件时,将该子信号特征对应的子预测信号特征作为该子信号特征对应的子模拟信号特征;
当未检测到信噪比时,将该子信号特征与该子信号特征对应的子预测信号特征进行加权处理,以得到该子信号特征对应的子模拟信号特征。
2.根据权利要求1所述麻醉深度的测量装置,其特征在于,所述预测脑电信号的生成过程具体包括:
以预先建立的神经元群模块为元素建立矩阵型脑网络;
根据预设联通参数驱动所述矩阵型脑网络产生脑电信号,其中,所述预设联通参数用于控制所述矩阵型脑网络中各神经元群模块状态。
3.根据权利要求2所述麻醉深度的测量装置,其特征在于,所述神经元群模块包括激励神经细胞网络、锥体细胞网络以及抑制神经细胞网络;所述激励神经细胞网络的激励信号、抑制神经细胞网络的抑制信号、以及锥体细胞网络的外部激励信号形成脑电波,并将所述脑电波输出至外部的神经元群模块,并且所述脑电波分别反馈至该神经元群模块的激励神经细胞网络以及抑制神经细胞网络。
4.根据权利要求1所述麻醉深度的测量装置,其特征在于,所述获取脑电信号,并提取所述脑电信号的信号特征具体包括:
获取脑电信号,并对所述脑电信号进行去噪处理;
提取去噪处理后的脑电信号的信号特征。
5.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现如权利要求1~4任意一项所述的麻醉深度的测量装置中的步骤。
6.一种电子设备,其特征在于,包括:处理器、存储器及通信总线;所述存储器上存储有可被所述处理器执行的计算机可读程序;
所述通信总线实现处理器和存储器之间的连接通信;
所述处理器执行所述计算机可读程序时实现如权利要求1-4任意一项所述的麻醉深度的测量装置中的步骤。
CN201911245011.5A 2019-12-06 2019-12-06 一种麻醉深度的测量装置、存储介质及电子设备 Active CN110840411B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201911245011.5A CN110840411B (zh) 2019-12-06 2019-12-06 一种麻醉深度的测量装置、存储介质及电子设备
PCT/CN2020/106066 WO2021109601A1 (zh) 2019-12-06 2020-07-31 一种麻醉深度的测量方法、存储介质及电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911245011.5A CN110840411B (zh) 2019-12-06 2019-12-06 一种麻醉深度的测量装置、存储介质及电子设备

Publications (2)

Publication Number Publication Date
CN110840411A CN110840411A (zh) 2020-02-28
CN110840411B true CN110840411B (zh) 2022-03-11

Family

ID=69608120

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911245011.5A Active CN110840411B (zh) 2019-12-06 2019-12-06 一种麻醉深度的测量装置、存储介质及电子设备

Country Status (2)

Country Link
CN (1) CN110840411B (zh)
WO (1) WO2021109601A1 (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110840411B (zh) * 2019-12-06 2022-03-11 深圳市德力凯医疗设备股份有限公司 一种麻醉深度的测量装置、存储介质及电子设备
CN113133744A (zh) * 2021-04-30 2021-07-20 鹤壁市人民医院 一种多功能麻醉科用麻醉深度监测装置
CN114521904B (zh) * 2022-01-25 2023-09-26 中山大学 一种基于耦合神经元群的脑电活动模拟方法及系统
CN116491960B (zh) * 2023-06-28 2023-09-19 南昌大学第一附属医院 脑瞬态监测设备、电子设备及存储介质
CN117281534B (zh) * 2023-11-22 2024-03-22 广东省人民医院 一种多指标的麻醉状态监测方法及系统

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101259015A (zh) * 2007-03-06 2008-09-10 李小俚 一种脑电信号分析监测方法及其装置
CN101488162A (zh) * 2008-01-17 2009-07-22 复旦大学 一种用于脑电信号自动评估的脑电信号特征提取方法
JP4836140B2 (ja) * 2007-03-23 2011-12-14 独立行政法人産業技術総合研究所 脳活動解析方法および装置
WO2011119224A3 (en) * 2010-03-24 2012-03-08 Steven John Schiff Model based control of parkinson's disease
CN102499675A (zh) * 2011-10-27 2012-06-20 杭州电子科技大学 一种皮层脑电信号的反馈系统随机共振增强方法
CN102509282A (zh) * 2011-09-26 2012-06-20 东南大学 一种融合结构连接的各脑区间的效能连接分析方法
JP2012146116A (ja) * 2011-01-12 2012-08-02 Kyushu Institute Of Technology スピーチ内容を識別する装置及び方法
CN103006211A (zh) * 2013-01-17 2013-04-03 西安电子科技大学 一种基于脑电网络分析的地形图描绘装置
KR20180059985A (ko) * 2016-11-28 2018-06-07 참엔지니어링(주) 히든 마르코프 모델을 이용한 마취 심도 측정 방법 및 장치
KR101911506B1 (ko) * 2016-12-27 2018-10-25 울산과학기술원 뇌신호 기반 3차원 상지운동 외부 보조기기 제어를 위한 가상 뇌파 생성장치와 시뮬레이션 장치 및 방법
CN109065128A (zh) * 2018-09-28 2018-12-21 郑州大学 一种加权图正则化稀疏脑网络构建方法
CN109701160A (zh) * 2019-01-23 2019-05-03 中国人民解放军总医院 影像引导下可见光定位导航装置及方法
CN110114115A (zh) * 2016-10-25 2019-08-09 脑路有限公司 用于预测治疗结果的设备及方法
CN110393525A (zh) * 2019-06-18 2019-11-01 浙江大学 一种基于深度循环自编码器的大脑活动检测方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6731975B1 (en) * 2000-10-16 2004-05-04 Instrumentarium Corp. Method and apparatus for determining the cerebral state of a patient with fast response
KR101248055B1 (ko) * 2011-05-24 2013-03-26 한국과학기술원 뇌활성도 및 마취심도 변화에 따른 eeg신호 모델 및 시뮬레이터
CN105769146A (zh) * 2016-03-24 2016-07-20 美合实业(苏州)有限公司 一种多监测指标的麻醉深度监护仪
US11241186B2 (en) * 2017-01-03 2022-02-08 Myndlift Ltd. Systems and methods for processing EEG signals of a neurofeedback protocol
CN107595247A (zh) * 2017-08-29 2018-01-19 深圳市德力凯医疗设备股份有限公司 一种基于脑电信号的麻醉深度的监测方法及系统
US11000224B2 (en) * 2017-09-05 2021-05-11 Korea University Research And Business Foundation Method and apparatus of monitoring anaesthesia and consciousness depth through brain network analysis
WO2019127558A1 (zh) * 2017-12-29 2019-07-04 深圳迈瑞生物医疗电子股份有限公司 基于脑电的麻醉深度监测方法和装置
CN109645989B (zh) * 2018-12-10 2021-01-08 燕山大学 一种麻醉深度估计系统
CN110840411B (zh) * 2019-12-06 2022-03-11 深圳市德力凯医疗设备股份有限公司 一种麻醉深度的测量装置、存储介质及电子设备

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101259015A (zh) * 2007-03-06 2008-09-10 李小俚 一种脑电信号分析监测方法及其装置
JP4836140B2 (ja) * 2007-03-23 2011-12-14 独立行政法人産業技術総合研究所 脳活動解析方法および装置
CN101488162A (zh) * 2008-01-17 2009-07-22 复旦大学 一种用于脑电信号自动评估的脑电信号特征提取方法
WO2011119224A3 (en) * 2010-03-24 2012-03-08 Steven John Schiff Model based control of parkinson's disease
JP2012146116A (ja) * 2011-01-12 2012-08-02 Kyushu Institute Of Technology スピーチ内容を識別する装置及び方法
CN102509282A (zh) * 2011-09-26 2012-06-20 东南大学 一种融合结构连接的各脑区间的效能连接分析方法
CN102499675A (zh) * 2011-10-27 2012-06-20 杭州电子科技大学 一种皮层脑电信号的反馈系统随机共振增强方法
CN103006211A (zh) * 2013-01-17 2013-04-03 西安电子科技大学 一种基于脑电网络分析的地形图描绘装置
CN110114115A (zh) * 2016-10-25 2019-08-09 脑路有限公司 用于预测治疗结果的设备及方法
KR20180059985A (ko) * 2016-11-28 2018-06-07 참엔지니어링(주) 히든 마르코프 모델을 이용한 마취 심도 측정 방법 및 장치
KR101911506B1 (ko) * 2016-12-27 2018-10-25 울산과학기술원 뇌신호 기반 3차원 상지운동 외부 보조기기 제어를 위한 가상 뇌파 생성장치와 시뮬레이션 장치 및 방법
CN109065128A (zh) * 2018-09-28 2018-12-21 郑州大学 一种加权图正则化稀疏脑网络构建方法
CN109701160A (zh) * 2019-01-23 2019-05-03 中国人民解放军总医院 影像引导下可见光定位导航装置及方法
CN110393525A (zh) * 2019-06-18 2019-11-01 浙江大学 一种基于深度循环自编码器的大脑活动检测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A permutation Lempel-Ziv complexity measure for EEG analysis;YangBai,等;《Biomedical Signal Processing and Control》;20150531;第19卷;102-114 *
ELECTROENCEPHALOGRAM AND VISUAL-EVOKED POTENTIAL GENERATION IN A MATHEMATICAL-MODEL OF COUPLED CORTICAL COLUMNS;JANSEN, BH,等;《BIOLOGICAL CYBERNETICS》;19950930;第73卷(第4期);357-366 *
Neural mass model-based tracking of anesthetic brain states;LevinKuhlmann,等;《NeuroImage》;20160630;1743–1755 *
Time rescaling reproduces EEG behavior during transition from propofol anesthesia-induced unconsciousness to consciousness;Boussen, S等;《SCIENTIFIC REPORTS》;20180416;第8卷;1-17 *

Also Published As

Publication number Publication date
WO2021109601A1 (zh) 2021-06-10
CN110840411A (zh) 2020-02-28

Similar Documents

Publication Publication Date Title
CN110840411B (zh) 一种麻醉深度的测量装置、存储介质及电子设备
CN110859600A (zh) 一种生成脑电信号的方法、存储介质及电子设备
CN111920408B (zh) 结合虚拟现实的脑电神经反馈系统的信号分析方法及组件
EP3014290A1 (en) Systems and methods for tracking non-stationary spectral structure and dynamics in physiological data
Ortiz et al. Analysis of the EEG rhythms based on the empirical mode decomposition during motor imagery when using a lower-limb exoskeleton. A case study
AU2001213824A1 (en) Monitoring auditory evoked potentials
EP1272106A1 (en) Monitoring auditory evoked potentials
Bashashati et al. Fractal dimension-based EEG biofeedback system
Nair et al. Chaotic analysis of the electroretinographic signal for diagnosis
Yan et al. Removal of transcranial alternating current stimulation eeg artifacts using blind source separation and wavelets
CN110931123B (zh) 一种矩阵型脑网络及其构建方法
KR20100027460A (ko) 다중 윈도우를 이용한 생체신호의 특성 파라미터 추출 장치및 방법
Wang et al. A Novel Rapid Assessment of Mental Stress by Using PPG Signals Based on Deep Learning
Liu et al. Fractal dimension based neurofeedback training to improve cognitive abilities
CN116313029B (zh) 一种数字针灸动态控制优化的方法、系统和装置
Hegde et al. A review on ECG signal processing and HRV analysis
CN115363523A (zh) 一种癫痫发作区定位系统、设备及介质
CN109124605B (zh) 一种减少icu内误报警的方法、装置及设备
RU2766409C2 (ru) Система и способ для оценивания объема крови в мозге, и/или кровотока в мозге, и/или глубины анестезии пациента
Chen et al. Research on Working Memory States Based on Weighted-Order Propagation Number Algorithm: An EEG Perspective
CN117017234B (zh) 多参数集成麻醉监测与分析系统
Jefri et al. EYE BLINK IDENTIFICATION AND REMOVAL FROM SINGLE-CHANNEL EEG USING EMD WITH ENERGY THRESHOLD AND ADAPTIVE FILTER
CN116491958B (zh) 靶点确定设备、电子设备及存储介质
Piaggi et al. Adaptive filtering and random variables coefficient for analyzing functional magnetic resonance imaging data
Kilby Wavelet analysis and classification surface electromyography signals

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 518000, 18th Floor, Building B, High tech Innovation Center, Guangqiao Road, Tianliao Community, Yutang Street, Guangming District, Shenzhen, Guangdong Province

Patentee after: Shenzhen Delikai Medical Electronics Co.,Ltd.

Address before: 6 / F, building 10, Guanlong Village Second Industrial Zone, Xili street, Nanshan District, Shenzhen City, Guangdong Province

Patentee before: SHENZHEN DELICA MEDICAL EQUIPMENT Co.,Ltd.