CN114366126A - A method and system for closed-loop EEG regulation - Google Patents
A method and system for closed-loop EEG regulation Download PDFInfo
- Publication number
- CN114366126A CN114366126A CN202210105120.2A CN202210105120A CN114366126A CN 114366126 A CN114366126 A CN 114366126A CN 202210105120 A CN202210105120 A CN 202210105120A CN 114366126 A CN114366126 A CN 114366126A
- Authority
- CN
- China
- Prior art keywords
- eeg
- signal
- theta rhythm
- ultrasonic
- eeg signal
- 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
- 230000033228 biological regulation Effects 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000033764 rhythmic process Effects 0.000 claims abstract description 113
- 230000000638 stimulation Effects 0.000 claims abstract description 85
- 238000012545 processing Methods 0.000 claims abstract description 22
- 230000008569 process Effects 0.000 claims abstract description 20
- 238000007781 pre-processing Methods 0.000 claims abstract description 12
- 210000004556 brain Anatomy 0.000 claims abstract description 10
- 238000011217 control strategy Methods 0.000 claims description 13
- 238000001914 filtration Methods 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 7
- 230000003321 amplification Effects 0.000 claims description 3
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 8
- 239000003550 marker Substances 0.000 abstract description 3
- 230000004007 neuromodulation Effects 0.000 description 6
- 238000002604 ultrasonography Methods 0.000 description 6
- 230000001537 neural effect Effects 0.000 description 4
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 230000013016 learning Effects 0.000 description 3
- 210000005036 nerve Anatomy 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000000971 hippocampal effect Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 241000283984 Rodentia Species 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000003925 brain function Effects 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000003920 cognitive function Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 210000001153 interneuron Anatomy 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 230000006386 memory function Effects 0.000 description 1
- 210000000537 nasal bone Anatomy 0.000 description 1
- 230000007383 nerve stimulation Effects 0.000 description 1
- 230000007230 neural mechanism Effects 0.000 description 1
- 230000010004 neural pathway Effects 0.000 description 1
- 210000000118 neural pathway Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 210000001176 projection neuron Anatomy 0.000 description 1
- 230000010831 regulation of synaptic plasticity Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 210000003625 skull Anatomy 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/369—Electroencephalography [EEG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/369—Electroencephalography [EEG]
- A61B5/372—Analysis of electroencephalograms
- A61B5/374—Detecting the frequency distribution of signals, e.g. detecting delta, theta, alpha, beta or gamma waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/369—Electroencephalography [EEG]
- A61B5/377—Electroencephalography [EEG] using evoked responses
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Animal Behavior & Ethology (AREA)
- Psychiatry (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Signal Processing (AREA)
- Psychology (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本申请提供一种闭环脑电调控方法及系统,所述方法包括如下步骤:接收实验对象的脑电信号;预处理所述脑电信号,得到预处理脑电信号;计算处理所述预处理脑电信号,得到θ节律脑电信号和脑电特征值;对所述θ节律脑电信号和所述脑电特征值进行处理,得到超声刺激参数和θ节律设定刺激相位;根据所述超声刺激参数和所述θ节律设定刺激相位,对所述实验对象发出超声刺激以进行脑电调控。本申请通过实时采集分析实验对象的脑电信号,使用闭环超声刺激技术,实现以脑电θ节律作为标志物对实验对象进行脑电调控。
The present application provides a closed-loop EEG regulation method and system. The method includes the following steps: receiving an EEG signal of an experimental subject; preprocessing the EEG signal to obtain a preprocessed EEG signal; calculating and processing the preprocessed brain signal electrical signal, obtain the theta rhythm EEG signal and EEG characteristic value; process the theta rhythm EEG signal and the EEG characteristic value to obtain ultrasonic stimulation parameters and theta rhythm to set the stimulation phase; according to the ultrasonic stimulation The parameters and the theta rhythm set the stimulation phase, and ultrasonic stimulation was sent to the subject for EEG modulation. In the present application, the EEG signal of the experimental object is collected and analyzed in real time, and the closed-loop ultrasonic stimulation technology is used to realize the EEG regulation of the experimental object by using the EEG theta rhythm as a marker.
Description
技术领域technical field
本申请涉及超声刺激神经调控技术领域,尤其涉及一种闭环脑电调控方法及系统。The present application relates to the technical field of ultrasonic stimulation nerve regulation, and in particular, to a closed-loop EEG regulation method and system.
背景技术Background technique
脑神经调控技术是利用光、磁、电、超声等物理或化学等外部技术手段改变脑部内源神经信号传递,从而引起脑功能变化的重要方法。其中超声神经调控是一种新型的脑刺激技术,相对于其他传统的电、磁、光神经刺激技术而言,具有方向性强、穿透深度大、靶点控制准确和无创等优点,因而近年来超声无创脑神经调控技术倍受神经科学领域学者的关注。Brain nerve regulation technology is an important method to change the transmission of endogenous nerve signals in the brain by using external technical means such as light, magnetism, electricity, ultrasound and other physical or chemical means, thereby causing changes in brain function. Among them, ultrasonic neuromodulation is a new type of brain stimulation technology. Compared with other traditional electrical, magnetic and optical nerve stimulation technologies, it has the advantages of strong directionality, large penetration depth, accurate target control and non-invasiveness. Ultrasound non-invasive brain neuromodulation technology has attracted the attention of scholars in the field of neuroscience.
超声作为神经调控技术,可以广泛应用于细胞、啮齿类动物、非人灵长类动物及人类的神经调控研究中,是一种理想的神经调控工具。然而,和电刺激神经调控技术不同,在之前的经颅超声刺激研究中,还很少出现闭环超声系统的使用,在调控时不能在刺激的过程中根据对象的特征情况实时改变刺激参数。As a neuromodulation technology, ultrasound can be widely used in neuromodulation research in cells, rodents, non-human primates and humans, and is an ideal neuromodulation tool. However, unlike electrical stimulation neuromodulation techniques, in previous studies of transcranial ultrasound stimulation, the use of closed-loop ultrasound systems was rarely used, and the stimulation parameters could not be changed in real time according to the characteristics of the object during the stimulation process.
此外,在哺乳动物的脑电中一个很重要的节律就是θ节律,θ节律是由海马中间神经元和锥体神经元在局部的相互作用产生的,θ节律对多个类型的学习和记忆过程起到关键的作用,同时对突触可塑性的调节作用也至关重要。然而,现有的闭环脑电调控方法均未考虑θ节律。为此,本申请提出一种闭环脑电调控方法及系统。In addition, a very important rhythm in mammalian EEG is theta rhythm. Theta rhythm is generated by the local interaction of hippocampal interneurons and pyramidal neurons. Theta rhythm is important for many types of learning and memory processes. It also plays a key role in the regulation of synaptic plasticity. However, none of the existing closed-loop EEG regulation methods consider theta rhythm. To this end, the present application proposes a closed-loop EEG regulation method and system.
发明内容SUMMARY OF THE INVENTION
本申请的目的是针对以上问题,提供一种闭环脑电调控方法及系统。The purpose of this application is to provide a closed-loop EEG regulation method and system for the above problems.
第一方面,本申请提供一种闭环脑电调控方法,所述方法包括如下步骤:In a first aspect, the present application provides a closed-loop EEG regulation method, the method includes the following steps:
接收实验对象的脑电信号;Receive the EEG signals of the experimental subjects;
预处理所述脑电信号,得到预处理脑电信号;Preprocessing the EEG signal to obtain a preprocessed EEG signal;
计算处理所述预处理脑电信号,得到θ节律脑电信号和脑电特征值;Calculating and processing the preprocessed EEG signal to obtain the theta rhythm EEG signal and EEG characteristic value;
对所述θ节律脑电信号和所述脑电特征值进行处理,得到超声刺激参数和θ节律设定刺激相位;processing the theta rhythm EEG signal and the EEG characteristic value to obtain ultrasonic stimulation parameters and theta rhythm setting stimulation phase;
根据所述超声刺激参数和所述θ节律设定刺激相位,对所述实验对象发出超声刺激以进行脑电调控。The stimulation phase is set according to the ultrasonic stimulation parameters and the theta rhythm, and ultrasonic stimulation is sent to the experimental subject for EEG regulation.
根据本申请某些实施例提供的技术方案,预处理所述脑电信号,得到预处理脑电信号,具体包括如下步骤:According to the technical solutions provided by some embodiments of the present application, preprocessing the EEG signal to obtain the preprocessed EEG signal specifically includes the following steps:
对所述脑电信号进行降采样、50Hz陷波去噪和4-200Hz滤波处理,得到预处理脑电信号。The EEG signal is subjected to down-sampling, 50Hz notch denoising and 4-200Hz filtering to obtain a preprocessed EEG signal.
根据本申请某些实施例提供的技术方案,计算处理所述预处理脑电信号,得到θ节律脑电信号和脑电特征值,具体包括如下步骤:According to the technical solutions provided by some embodiments of the present application, calculating and processing the preprocessed EEG signal to obtain the theta rhythm EEG signal and EEG characteristic value, which specifically includes the following steps:
对所述预处理脑电信号进行4-8Hz滤波,得到θ节律脑电信号;4-8Hz filtering is performed on the preprocessed EEG signal to obtain the theta rhythm EEG signal;
计算处理所述预处理脑电信号和所述θ节律脑电信号,得到所述脑电特征值。Calculate and process the preprocessed EEG signal and the theta rhythm EEG signal to obtain the EEG characteristic value.
根据本申请某些实施例提供的技术方案,对所述θ节律脑电信号和所述脑电特征值进行处理,得到超声刺激参数和θ节律设定刺激相位,具体包括如下步骤:According to the technical solutions provided by some embodiments of the present application, the theta rhythm EEG signal and the EEG characteristic value are processed to obtain ultrasonic stimulation parameters and theta rhythm setting stimulation phase, which specifically includes the following steps:
计算所述脑电特征值和预设期望值的误差,并按照控制策略得到超声调控参数;Calculate the error between the EEG characteristic value and the preset expected value, and obtain ultrasonic regulation parameters according to the control strategy;
采用希尔波尔变换计算所述θ节律脑电信号的θ节律相位;Calculate the theta rhythm phase of the theta rhythm EEG signal by using Hillbore transform;
根据θ节律相位和预测模型预测所述θ节律脑电信号的θ节律设定刺激相位。The theta rhythm setting stimulation phase of the theta rhythm EEG signal is predicted according to the theta rhythm phase and the prediction model.
根据本申请某些实施例提供的技术方案,所述控制策略为广义最小方差控制策略。According to the technical solutions provided by some embodiments of the present application, the control strategy is a generalized minimum variance control strategy.
根据本申请某些实施例提供的技术方案,根据所述超声刺激参数和所述θ节律设定刺激相位,对所述实验对象发出超声刺激以进行脑电调控,具体包括如下步骤:According to the technical solutions provided by some embodiments of the present application, the stimulation phase is set according to the ultrasonic stimulation parameters and the theta rhythm, and ultrasonic stimulation is sent to the experimental subject to perform EEG regulation, which specifically includes the following steps:
根据所述超声刺激参数和所述θ节律设定刺激相位生成第一控制波形信号;generating a first control waveform signal according to the ultrasonic stimulation parameters and the theta rhythm setting stimulation phase;
放大处理所述第一控制波形信号,得到第二控制波形信号;amplifying and processing the first control waveform signal to obtain a second control waveform signal;
将所述第二控制波形信号使能为超声波信号输出;enabling the second control waveform signal to be output as an ultrasonic signal;
将所述超声波信号限定在设定空间范围内并向所述实验对象发射。The ultrasonic signal is confined within a set spatial range and emitted to the experimental subject.
第二方面、本申请提供一种上述闭环脑电调控方法的闭环脑电调控系统,包括:In the second aspect, the present application provides a closed-loop EEG regulation system of the above-mentioned closed-loop EEG regulation method, including:
处理模块,配置用于接收并预处理所述实验对象的脑电信号,以得到预处理脑电信号,还配置用于计算处理所述预处理脑电信号,以得到θ节律脑电信号和脑电特征值;The processing module is configured to receive and preprocess the EEG signal of the experimental subject to obtain the preprocessed EEG signal, and is also configured to calculate and process the preprocessed EEG signal to obtain the theta rhythm EEG signal and the brain electrical characteristic value;
控制模块,配置用于对所述θ节律脑电信号和所述脑电特征值进行处理,以得到超声刺激参数和θ节律设定刺激相位;a control module, configured to process the theta rhythm EEG signal and the EEG characteristic value to obtain ultrasonic stimulation parameters and theta rhythm setting stimulation phase;
超声刺激模块,配置用于根据所述超声刺激参数和所述θ节律设定刺激相位,对所述实验对象发出超声刺激以进行脑电调控。The ultrasonic stimulation module is configured to set a stimulation phase according to the ultrasonic stimulation parameters and the theta rhythm, and send ultrasonic stimulation to the experimental subject to perform EEG regulation.
根据本申请某些实施例提供的技术方案,所述处理模块包括:According to the technical solutions provided by some embodiments of the present application, the processing module includes:
预处理模块,配置用于接收所述脑电信号,并对所述脑电信号进行降采样、陷波去噪和滤波处理,以得到预处理脑电信号;a preprocessing module, configured to receive the EEG signal, and perform down-sampling, notch denoising and filtering processing on the EEG signal to obtain a preprocessed EEG signal;
θ节律脑电计算模块,配置用于对所述预处理脑电信号进行4-8Hz滤波,得到θ节律脑电信号;The theta rhythm EEG calculation module, configured to filter the preprocessed EEG signal at 4-8 Hz to obtain the theta rhythm EEG signal;
脑电特征值计算模块,配置用于计算处理所述预处理脑电信号和所述θ节律脑电信号,得到所述脑电特征值。The EEG feature value calculation module is configured to calculate and process the preprocessed EEG signal and the theta rhythm EEG signal to obtain the EEG feature value.
根据本申请某些实施例提供的技术方案,所述控制模块包括:According to the technical solutions provided by some embodiments of the present application, the control module includes:
θ节律相位预测模块,配置用于采用希尔波尔变换计算所述θ节律脑电信号的θ节律相位,还配置用于根据θ节律相位和预测模型预测所述θ节律脑电信号的θ节律设定刺激相位;The theta rhythm phase prediction module, configured to calculate the theta rhythm phase of the theta rhythm EEG signal by adopting Hillbore transform, and also configured to predict the theta rhythm of the theta rhythm EEG signal according to the theta rhythm phase and the prediction model Set the stimulus phase;
超声调控参数计算模块,配置用于计算所述脑电特征值和预设期望值的误差,并按照控制策略得到超声调控参数。The ultrasonic regulation parameter calculation module is configured to calculate the error between the EEG characteristic value and the preset expected value, and obtain the ultrasonic regulation parameter according to the control strategy.
根据本申请某些实施例提供的技术方案,所述超声刺激模块包括:According to the technical solutions provided by some embodiments of the present application, the ultrasonic stimulation module includes:
信号发生模块,配置用于根据所述超声刺激参数和所述θ节律设定刺激相位生成第一控制波形信号;a signal generating module, configured to generate a first control waveform signal according to the ultrasonic stimulation parameters and the theta rhythm setting stimulation phase;
功率放大模块,配置用于放大处理所述第一控制波形信号,得到第二控制波形信号;a power amplification module, configured to amplify and process the first control waveform signal to obtain a second control waveform signal;
超声换能模块,配置用于将所述第二控制波形信号使能为超声波信号输出;an ultrasonic transducer module, configured to enable the second control waveform signal to be output as an ultrasonic signal;
准直模块,配置用于将所述超声波信号限定在设定空间范围内并向所述实验对象发射。The collimation module is configured to confine the ultrasonic signal within a set spatial range and transmit the ultrasonic signal to the experimental object.
与现有技术相比,本申请的有益效果:本申请提供的闭环脑电调控方法包括:接收实验对象的脑电信号;预处理所述脑电信号,得到预处理脑电信号;计算处理所述预处理脑电信号,得到θ节律脑电信号和脑电特征值;对所述θ节律脑电信号和所述脑电特征值进行处理,得到超声刺激参数和θ节律设定刺激相位;根据所述超声刺激参数和所述θ节律设定刺激相位,对所述实验对象发出超声刺激以进行脑电调控;本申请通过实时采集分析实验对象的脑电信号,使用闭环超声刺激技术,实现以脑电θ节律作为标志物对实验对象进行脑电调控。Compared with the prior art, the beneficial effects of the present application: the closed-loop EEG regulation method provided by the present application includes: receiving the EEG signal of the experimental object; preprocessing the EEG signal to obtain the preprocessed EEG signal; The preprocessed EEG signal is described to obtain the theta rhythm EEG signal and the EEG characteristic value; the theta rhythm EEG signal and the EEG characteristic value are processed to obtain the ultrasonic stimulation parameters and the theta rhythm to set the stimulation phase; The ultrasonic stimulation parameters and the theta rhythm are set to stimulate the phase, and ultrasonic stimulation is sent to the experimental object for EEG regulation; the present application collects and analyzes the EEG signal of the experimental object in real time, and uses closed-loop ultrasonic stimulation technology to achieve The EEG theta rhythm was used as a marker to regulate the EEG of the experimental subjects.
θ节律的频段在4-8Hz,周期在125-250ms,方便在存在ms级设备运算延迟情况下,保证脑电闭环调控中的刺激准确度;脑电θ节律作用于学习和记忆过程,以脑电θ节律作为标志物进行闭环脑电调控,能够进一步调节学习记忆相关的神经作用机制和可塑性神经通路,有利于研究改善实验对象的学习记忆和认知功能。The frequency band of theta rhythm is 4-8Hz, and the period is 125-250ms, which is convenient to ensure the stimulation accuracy in the closed-loop regulation of EEG in the presence of ms-level equipment operation delay. The closed-loop EEG regulation of electrical theta rhythm as a marker can further regulate the neural mechanisms and plastic neural pathways related to learning and memory, which is beneficial to study and improve the learning and memory and cognitive functions of experimental subjects.
附图说明Description of drawings
图1为本申请实施例1提供的闭环脑电调控系统的结构示意图;1 is a schematic structural diagram of a closed-loop EEG regulation system provided in Embodiment 1 of the present application;
图2为本申请实施例2提供的闭环脑电调控方法的流程图;2 is a flowchart of a closed-loop EEG regulation method provided in Embodiment 2 of the present application;
图3为图2中步骤S3的具体流程图;Fig. 3 is the concrete flow chart of step S3 in Fig. 2;
图4为图2中步骤S4的具体流程图;Fig. 4 is the concrete flow chart of step S4 in Fig. 2;
图5为图2中步骤S5的具体流程图。FIG. 5 is a specific flowchart of step S5 in FIG. 2 .
具体实施方式Detailed ways
为了使本领域技术人员更好地理解本申请的技术方案,下面结合附图对本申请进行详细描述,本部分的描述仅是示范性和解释性,不应对本申请的保护范围有任何的限制作用。In order to make those skilled in the art better understand the technical solutions of the present application, the present application will be described in detail below with reference to the accompanying drawings. The description in this part is only exemplary and explanatory, and should not have any limiting effect on the protection scope of the present application. .
实施例1Example 1
请参考图1,本实施例提供一种闭环脑电调控系统,所述系统包括采集模块、处理模块、控制模块和超声刺激模块。Referring to FIG. 1 , this embodiment provides a closed-loop EEG regulation system, which includes an acquisition module, a processing module, a control module, and an ultrasonic stimulation module.
所述采集模块,连接在实验对象上,用于实时采集实验对象的脑电信号;其中,实验对象设置为实验室进行试验时常用的小黑鼠;所述采集模块包括脑电电极和神经信号处理器;所述脑电电极包括三路,其中一路钨微电极插在实验对象颅骨下的海马CA1区记录该脑区的脑电,另外两路电极插接在实验对象鼻骨上作为参比电极和接地电极;所述神经信号处理器与脑电电极相连接,配置用于将脑电信号转化为相应的数字信号,并将转换后的数字信号发送至脑电信号处理模块。The acquisition module is connected to the experimental object and is used for real-time acquisition of the EEG signal of the experimental object; wherein, the experimental object is set as a black mouse commonly used in experiments in the laboratory; the acquisition module includes EEG electrodes and neural signals processor; the EEG electrodes include three paths, one of which is inserted into the hippocampal CA1 area under the skull of the experimental object to record the EEG of the brain area, and the other two electrodes are inserted into the nasal bone of the experimental object as reference electrodes and the ground electrode; the neural signal processor is connected with the EEG electrode, and is configured to convert the EEG signal into a corresponding digital signal, and send the converted digital signal to the EEG signal processing module.
在本实施例中,所述神经信号处理器的型号可选的采用Bio-SignalTechnologies公司的ApolloⅡ神经信号处理器,采样频率为30kHz,通过USB接口与PC机相连接。In this embodiment, the model of the neural signal processor can be selected by the Apollo II neural signal processor of Bio-Signal Technologies, the sampling frequency is 30 kHz, and it is connected to a PC through a USB interface.
所述处理模块,与所述采集模块相连接,用于接收采集模块所采集的实验对象的脑电信号,并对所述脑电信号进行预处理以得到预处理脑电信号,还配置用于计算处理所述预处理脑电信号,以得到θ节律脑电信号和脑电特征值。The processing module, connected with the acquisition module, is used for receiving the EEG signal of the experimental object collected by the acquisition module, and preprocessing the EEG signal to obtain the preprocessed EEG signal, and is also configured to The pre-processed EEG signals are processed computationally to obtain theta rhythm EEG signals and EEG characteristic values.
具体地,所述处理模块包括:预处理模块、θ节律脑电计算模块和脑电特征值计算模块。Specifically, the processing module includes: a preprocessing module, a theta rhythm EEG calculation module, and an EEG feature value calculation module.
所述预处理模块与所述采集模块相连接,配置用于接收所述脑电信号,并对所述脑电信号进行降采样、50Hz陷波去噪和4-200Hz滤波处理,以得到预处理脑电信号,还配置用于将所述预处理脑电信号发送给所述θ节律脑电计算模块和所述脑电特征值计算模块。The preprocessing module is connected to the acquisition module, and is configured to receive the EEG signal, and perform down-sampling, 50Hz notch denoising and 4-200Hz filtering processing on the EEG signal to obtain preprocessing The EEG signal is further configured to send the preprocessed EEG signal to the theta rhythm EEG calculation module and the EEG feature value calculation module.
所述θ节律脑电计算模块,配置用于对所述预处理脑电信号进行4-8Hz滤波,以得到θ节律脑电信号,还配置用于将所述θ节律脑电信号发送至所述脑电特征值计算模块和所述控制模块。The theta rhythm EEG calculation module is configured to perform 4-8 Hz filtering on the preprocessed EEG signal to obtain the theta rhythm EEG signal, and is also configured to send the theta rhythm EEG signal to the EEG eigenvalue calculation module and the control module.
所述脑电特征值计算模块,配置用于计算处理所述预处理脑电信号和所述θ节律脑电信号,以得到所述脑电特征值,还配置用于将所述脑电特征值发送至所述控制模块。其中,脑电特征值包括预处理脑电信号的绝对功率、θ节律脑电信号相对于预处理脑电信号的相对功率、样本熵等。The EEG feature value calculation module is configured to calculate and process the preprocessed EEG signal and the theta rhythm EEG signal to obtain the EEG feature value, and is also configured to calculate the EEG feature value sent to the control module. The EEG characteristic value includes the absolute power of the preprocessed EEG signal, the relative power of the theta rhythm EEG signal relative to the preprocessed EEG signal, and the sample entropy.
所述控制模块,配置用于对所述θ节律脑电信号和所述脑电特征值进行处理,以得到超声刺激参数和θ节律设定刺激相位。The control module is configured to process the theta rhythm EEG signal and the EEG characteristic value to obtain ultrasonic stimulation parameters and theta rhythm setting stimulation phase.
具体地,所述控制模块包括:θ节律相位预测模块和超声调控参数计算模块。Specifically, the control module includes: a theta rhythm phase prediction module and an ultrasonic regulation parameter calculation module.
所述超声调控参数计算模块,与所述脑电特征计算模块相连接,配置用于计算所述脑电特征值和预设期望值的误差,并按照控制策略得到超声调控参数,还配置用于将所述超声调控参数发送至所述超声刺激模块。其中,所述超声调控参数包括超声大小、时长等参数。The ultrasonic regulation parameter calculation module is connected to the EEG characteristic calculation module, and is configured to calculate the error between the EEG characteristic value and the preset expected value, obtain ultrasonic regulation parameters according to the control strategy, and is also configured to The ultrasonic regulation parameters are sent to the ultrasonic stimulation module. Wherein, the ultrasonic regulation parameters include parameters such as ultrasonic size and duration.
所述θ节律相位预测模块,与所述θ节律脑电计算模块相连接,配置用于接收所述θ节律脑电信号,并采用希尔波尔变换计算所述θ节律脑电信号的θ节律相位,还配置用于通过预测模型预测θ节律设定刺激相位的时间得到超声刺激时间,并根据超声刺激时间将超声输出指令发送至所述超声刺激模块;所述θ节律预测模型可以根据历史脑电θ节律对未来一段时间的θ节律相位进行预测。The theta rhythm phase prediction module is connected to the theta rhythm EEG calculation module, and is configured to receive the theta rhythm EEG signal, and calculate the theta rhythm of the theta rhythm EEG signal using Hillbore transformation phase, and is also configured to predict the time of the theta rhythm to set the stimulation phase through the prediction model to obtain the ultrasonic stimulation time, and send the ultrasonic output instruction to the ultrasonic stimulation module according to the ultrasonic stimulation time; the theta rhythm prediction model can be based on the historical brain. The electrical theta rhythm predicts the theta rhythm phase in the future.
所述控制策略能够根据当前脑电特征值与预设期望值的误差对超声参数进行调控,达到减小甚至消除误差的效果;在本实施例中,所述控制策略使用的是广义最小方差控制策略,θ节律相位预测模型是AR模型,所述θ节律设定刺激相位为波峰。The control strategy can adjust the ultrasonic parameters according to the error between the current EEG characteristic value and the preset expected value, so as to achieve the effect of reducing or even eliminating the error; in this embodiment, the control strategy uses a generalized minimum variance control strategy , the theta rhythm phase prediction model is an AR model, and the theta rhythm sets the stimulus phase as a peak.
在本实施例中,所述处理模块和所述控制模块分别由PC机中C++程序与MATLAB程序实现,所述C++程序与MATLAB程序之间的数据传输由UDP数据包实现。In this embodiment, the processing module and the control module are respectively implemented by a C++ program and a MATLAB program in the PC, and the data transmission between the C++ program and the MATLAB program is implemented by a UDP data packet.
所述超声刺激模块,配置用于根据所述超声刺激参数和所述θ节律设定刺激相位,对所述实验对象发出超声刺激以进行脑电调控。The ultrasonic stimulation module is configured to set a stimulation phase according to the ultrasonic stimulation parameters and the theta rhythm, and send ultrasonic stimulation to the experimental subject to perform EEG regulation.
具体地,所述超声刺激模块包括:信号发生模块、功率放大模块、超声换能模块和准直模块。Specifically, the ultrasonic stimulation module includes: a signal generation module, a power amplification module, an ultrasonic transducer module and a collimation module.
所述信号发生模块,配置用于接收所述θ节律相位预测模块发送的输出指令、和所述调控参数计算模块发送的超声刺激参数,并生成第一控制波形信号。The signal generation module is configured to receive the output instruction sent by the theta rhythm phase prediction module and the ultrasonic stimulation parameters sent by the regulation parameter calculation module, and generate a first control waveform signal.
所述信号发生模块是可以实时远程接收超声频率、脉冲重复频率、脉冲周期数目、脉冲串长度、脉冲串间隔、串循环次数、开始刺激时间等参数的函数/任意波形发生器;在本实施例中,所述信号发生模块可选的采用RIGOL品牌的DG2080型号函数/任意波形发生器,通过USB接口与PC机连接,通过BNC线将信号输出给超声换能器,所接收的控制指令由MATLAB程序发送。The signal generating module is a function/arbitrary waveform generator that can remotely receive parameters such as ultrasonic frequency, pulse repetition frequency, pulse cycle number, pulse train length, pulse train interval, train cycle times, and start stimulation time in real time. , the signal generation module can optionally use the DG2080 model function/arbitrary waveform generator of the RIGOL brand, which is connected to the PC through the USB interface, and outputs the signal to the ultrasonic transducer through the BNC line, and the received control commands are generated by MATLAB program to send.
所述功率放大模块,配置用于将所述第一控制波形信号进行放大处理,得到第二控制波形信号,并将所述第二控制波形信号发送至所述超能换能模块。The power amplifying module is configured to amplify the first control waveform signal to obtain a second control waveform signal, and send the second control waveform signal to the super energy transducer module.
在本实施例中,所述功率放大模块由射频放大器实现,本优选实施例中,功率放大器的型号为240L的ENI线性功率放大器。In this embodiment, the power amplifying module is implemented by a radio frequency amplifier. In this preferred embodiment, the power amplifier is a 240L ENI linear power amplifier.
所述超声换能模块,配置用于将所述第二控制波形信号使能为超声波信号输出,并与所述准直模块相连接;所述超声换能模块为超声换能器,超声换能器是一种有源器件,可以实现将电功率转换为机械功率输出,即将电能转换为声能的一种装置,而且这个过程只有很小的功率损耗(超声波为一种机械波)。The ultrasonic transducer module is configured to enable the second control waveform signal to be output as an ultrasonic signal, and is connected to the collimation module; the ultrasonic transducer module is an ultrasonic transducer, and the ultrasonic transducer A device is an active device that can realize the conversion of electrical power into mechanical power output, that is, a device that converts electrical energy into sound energy, and this process has only a small power loss (ultrasonic is a mechanical wave).
在本实施例中,所述超声换能器产生超声,其中心频率设为500kHz,选用直径31mm的奥利巴斯超声换能器,型号为V301-SU。In this embodiment, the ultrasonic transducer generates ultrasonic waves, and its center frequency is set to 500 kHz, and an Olibas ultrasonic transducer with a diameter of 31 mm is selected, and the model is V301-SU.
所述准直模块固定连接在所述超声换能器的输出端,配置用于将所述超声波信号限定在设定空间范围内并向所述实验对象发射,同时通过在其中填注耦合液,减少超声在传播过程中的衰减;在本实施例中,所述准直模块为准直器。The collimation module is fixedly connected to the output end of the ultrasonic transducer, and is configured to limit the ultrasonic signal within a set space range and transmit it to the experimental object, and at the same time, by filling the coupling fluid therein, Reduce the attenuation of ultrasound during propagation; in this embodiment, the collimation module is a collimator.
实施例2Example 2
本实施例提供一种应用实施例1所述闭环脑电调控系统的闭环脑电调控方法,图2为该方法的流程图,所述方法包括如下步骤:This embodiment provides a closed-loop EEG regulation method using the closed-loop EEG regulation system described in Embodiment 1. FIG. 2 is a flowchart of the method, and the method includes the following steps:
S1、接收实验对象的脑电信号。S1. Receive the EEG signal of the experimental object.
在本实施例中,实验对象为实验室进行试验时常用的小黑鼠。In this embodiment, the experimental object is a black mouse commonly used in laboratory experiments.
S2、预处理所述脑电信号,得到预处理脑电信号。S2. Preprocess the EEG signal to obtain a preprocessed EEG signal.
本步骤具体包括:对所述脑电信号进行降采样、50Hz陷波去噪和4-200Hz滤波处理,得到预处理脑电信号。This step specifically includes: performing down-sampling, 50Hz notch denoising and 4-200Hz filtering processing on the EEG signal to obtain a pre-processed EEG signal.
S3、计算处理所述预处理脑电信号,得到θ节律脑电信号和脑电特征值。S3. Calculate and process the preprocessed EEG signal to obtain the theta rhythm EEG signal and EEG characteristic value.
如图3所示,本步骤具体包括:As shown in Figure 3, this step specifically includes:
S31、对所述预处理脑电信号进行4-8Hz滤波,得到θ节律脑电信号;S31, performing 4-8 Hz filtering on the preprocessed EEG signal to obtain the theta rhythm EEG signal;
S32、计算处理所述预处理脑电信号和所述θ节律脑电信号,得到所述脑电特征值。S32. Calculate and process the preprocessed EEG signal and the theta rhythm EEG signal to obtain the EEG feature value.
S4、对所述θ节律脑电信号和所述脑电特征值进行处理,得到超声刺激参数和θ节律设定刺激相位。S4. Process the theta rhythm EEG signal and the EEG characteristic value to obtain ultrasonic stimulation parameters and theta rhythm setting stimulation phase.
如图4所示,本步骤具体包括:As shown in Figure 4, this step specifically includes:
S41、计算所述脑电特征值和预设期望值的误差,并按照控制策略得到超声调控参数;S41, calculating the error between the EEG characteristic value and the preset expected value, and obtaining ultrasonic regulation parameters according to the control strategy;
S42、采用希尔波尔变换计算所述θ节律脑电信号的θ节律相位;S42, using Hillbore transform to calculate the theta rhythm phase of the theta rhythm EEG signal;
S43、根据θ节律相位和预测模型预测所述θ节律脑电信号的θ节律设定刺激相位。S43. Predict the theta rhythm setting stimulation phase of the theta rhythm EEG signal according to the theta rhythm phase and the prediction model.
S5、根据所述超声刺激参数和所述θ节律设定刺激相位,对所述实验对象发出超声刺激以进行脑电调控。S5. Set a stimulation phase according to the ultrasonic stimulation parameters and the theta rhythm, and send ultrasonic stimulation to the experimental subject to perform EEG regulation.
如图5所示,本步骤具体包括:As shown in Figure 5, this step specifically includes:
S51、根据所述超声刺激参数和所述θ节律设定刺激相位生成第一控制波形信号;S51, setting a stimulation phase according to the ultrasonic stimulation parameters and the theta rhythm to generate a first control waveform signal;
S52、放大处理所述第一控制波形信号,得到第二控制波形信号;S52. Amplify and process the first control waveform signal to obtain a second control waveform signal;
S53、将所述第二控制波形信号使能为超声波信号输出;S53, enabling the second control waveform signal to be output as an ultrasonic signal;
S54、将所述超声波信号限定在设定空间范围内并向所述实验对象发射。S54. Confine the ultrasonic signal within a set space range and transmit it to the experimental object.
本文中应用了具体个例对本申请的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本申请的方法及其核心思想。以上所述仅是本申请的优选实施方式,应当指出,由于文字表达的有限性,而客观上存在无限的具体结构,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进、润饰或变化,也可以将上述技术特征以适当的方式进行组合;这些改进润饰、变化或组合,或未经改进将发明的构思和技术方案直接应用于其他场合的,均应视为本申请的保护范围。Specific examples are used herein to illustrate the principles and implementations of the present application, and the descriptions of the above embodiments are only used to help understand the methods and core ideas of the present application. The above are only the preferred embodiments of the present application. It should be pointed out that due to the limited expression of words, there are objectively unlimited specific structures. For those of ordinary skill in the art, without departing from the principles of the present invention However, some improvements, modifications or changes can also be made, and the above-mentioned technical features can also be combined in an appropriate manner; these improvements, modifications, or combinations, or the ideas and technical solutions of the invention are directly applied to other occasions without improvement. , shall be regarded as the protection scope of this application.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210105120.2A CN114366126B (en) | 2022-01-28 | 2022-01-28 | A closed-loop EEG control method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210105120.2A CN114366126B (en) | 2022-01-28 | 2022-01-28 | A closed-loop EEG control method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114366126A true CN114366126A (en) | 2022-04-19 |
| CN114366126B CN114366126B (en) | 2025-03-18 |
Family
ID=81146661
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210105120.2A Active CN114366126B (en) | 2022-01-28 | 2022-01-28 | A closed-loop EEG control method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114366126B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115344122A (en) * | 2022-08-15 | 2022-11-15 | 中国科学院深圳先进技术研究院 | Acoustic non-invasive brain-computer interface and control method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060149337A1 (en) * | 2005-01-21 | 2006-07-06 | John Michael S | Systems and methods for tissue stimulation in medical treatment |
| US20130197401A1 (en) * | 2011-12-30 | 2013-08-01 | Tomo Sato | Optimization of ultrasound waveform characteristics for transcranial ultrasound neuromodulation |
| CN103394161A (en) * | 2013-07-16 | 2013-11-20 | 山西大学 | Cerebral alpha wave feedback regulation acupoint magnetic stimulation system |
| CN113181569A (en) * | 2021-04-27 | 2021-07-30 | 燕山大学 | Closed-loop transcranial brain stimulation system and method |
| CN113208610A (en) * | 2020-01-17 | 2021-08-06 | 广西医科大学 | Phase locus estimation method and device and closed-loop neural stimulation system |
| KR20210136627A (en) * | 2020-05-08 | 2021-11-17 | 광운대학교 산학협력단 | A method for predicting vagus nerve stimulation effect in pediatric epilepsy patients by brain disease verification eeg biomarker analysis algorithm |
-
2022
- 2022-01-28 CN CN202210105120.2A patent/CN114366126B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060149337A1 (en) * | 2005-01-21 | 2006-07-06 | John Michael S | Systems and methods for tissue stimulation in medical treatment |
| US20130197401A1 (en) * | 2011-12-30 | 2013-08-01 | Tomo Sato | Optimization of ultrasound waveform characteristics for transcranial ultrasound neuromodulation |
| CN103394161A (en) * | 2013-07-16 | 2013-11-20 | 山西大学 | Cerebral alpha wave feedback regulation acupoint magnetic stimulation system |
| CN113208610A (en) * | 2020-01-17 | 2021-08-06 | 广西医科大学 | Phase locus estimation method and device and closed-loop neural stimulation system |
| KR20210136627A (en) * | 2020-05-08 | 2021-11-17 | 광운대학교 산학협력단 | A method for predicting vagus nerve stimulation effect in pediatric epilepsy patients by brain disease verification eeg biomarker analysis algorithm |
| CN113181569A (en) * | 2021-04-27 | 2021-07-30 | 燕山大学 | Closed-loop transcranial brain stimulation system and method |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115344122A (en) * | 2022-08-15 | 2022-11-15 | 中国科学院深圳先进技术研究院 | Acoustic non-invasive brain-computer interface and control method |
| WO2024036794A1 (en) * | 2022-08-15 | 2024-02-22 | 中国科学院深圳先进技术研究院 | Acoustic-wave non-invasive brain-computer interface and control method |
| CN115344122B (en) * | 2022-08-15 | 2024-06-18 | 中国科学院深圳先进技术研究院 | Acoustic non-invasive brain-computer interface and control method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114366126B (en) | 2025-03-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113181569B (en) | Closed-loop transcranial brain stimulation system and method | |
| CN108865881B (en) | Cell function regulation and control system and method based on acoustomagnetic coupling electrical stimulation principle | |
| CN104826243B (en) | A kind of device of ultrasound stimulation nerve fiber | |
| CN205041400U (en) | Novel it is portable through cranium ultrasound stimulation and EEG signal collection device | |
| CN1089011C (en) | Functional neuromyo-electric signal identification and stimulation apparatus | |
| CN109924976A (en) | The stimulation of mouse TCD,transcranial Doppler and brain electromyography signal synchronous | |
| CN105251141A (en) | Closed-loop transcranial magnetoacoustic stimulation device | |
| CN105251125A (en) | Integrated device of transcranial magnetoacoustic stimulation and electroencephalogram detection | |
| CN108144199A (en) | Device, method and system are motivated by ultrasound | |
| CN104857640A (en) | Closed-loop type transcranial ultrasonic brain stimulation apparatus | |
| Qin et al. | Optimizing frequency and pulse shape for ultrasound current source density imaging | |
| CN106039600B (en) | A kind of wearable ultrasonic nerve stimulation device of wireless control | |
| CN115813344A (en) | Deep brain stimulation noninvasive monitoring system based on sound computer imaging | |
| CN114366126A (en) | A method and system for closed-loop EEG regulation | |
| CN119158176A (en) | A closed-loop transcranial alternating current stimulation system based on phase-amplitude coupled stimulation | |
| CN110975153B (en) | A deep brain stimulation configuration method and system, electronic device, and storage medium | |
| CN221998650U (en) | Transcranial magnetic stimulation system with electroencephalogram signal feedback regulation | |
| CN109513123B (en) | High-resolution three-dimensional passive cavitation imaging method based on hemispherical array | |
| CN212327192U (en) | A mental fatigue intervention device | |
| CN108355246A (en) | Optimal closed loop DBS treatment systems based on phase response curve | |
| CN111973880A (en) | Mental fatigue intervention equipment and intervention method | |
| CN116637020A (en) | A real-time sleep rhythm closed-loop control system based on ultrasonic stimulation of Shenmen acupoint | |
| CN115153586A (en) | A method and system for regulating sleep rhythm | |
| CN104548392A (en) | Transcranial ultrasonic stimulation device and stimulation method | |
| CN119215327B (en) | Directional deep brain stimulation closed-loop regulation and control system based on contact current multidimensional feedback |
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 |