CN106073702B - Based on the more time-frequency scale diencephalon myoelectricity coupling analytical methods of small echo-transfer entropy - Google Patents
Based on the more time-frequency scale diencephalon myoelectricity coupling analytical methods of small echo-transfer entropy Download PDFInfo
- Publication number
- CN106073702B CN106073702B CN201610362111.6A CN201610362111A CN106073702B CN 106073702 B CN106073702 B CN 106073702B CN 201610362111 A CN201610362111 A CN 201610362111A CN 106073702 B CN106073702 B CN 106073702B
- Authority
- CN
- China
- Prior art keywords
- signal
- brain
- wavelet
- myoelectricity
- entropy
- 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
Classifications
-
- 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/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/389—Electromyography [EMG]
-
- 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
Abstract
The invention discloses a kind of based on the more time-frequency scale diencephalon myoelectricity coupling analytical methods of small echo-transfer entropy, the method includes brain electromyography signal synchronous acquisition part and signal processing, brain electromyography signal synchronous acquisition part includes eeg signal acquisition and electromyographic signal collection;Signal processing includes small echo-transfer entropy analysis method of Signal Pretreatment and brain myoelectricity.The present invention has applicability, admissibility, has important application value in medical science of recovery therapy field.
Description
Technical field
The present invention relates to neural rehabilitation project and locomotory mechanism research fields, and in particular to one kind is based on small echo-transfer entropy
More time-frequency scale diencephalon myoelectricity coupling analytical methods.
Background technique
Brain electricity (electroencephalogram, EEG) and myoelectricity (electromyographic, EMG) signal wraps respectively
The functional response information that information and muscle are intended to brain control is controlled containing somatic movement, the multi-scale coupling between brain electromyography signal
The multi-level cortex of message reflection-muscle function coupling (Functional corticomuscular coupling, FCMC)
Link information.Currently, brain myoelectricity synchronous characteristic research is based primarily upon coherence analysis, obtains the driving of brain motor mindedness and transported with muscle
Functional cohesion feature between dynamic response, but traditional coherent analysis can not embody coupling direction character.To be best understood from
Function between cerebral cortex and respective muscle interacts and information transmission characteristic, and it is same that Granger Causality analysis is applied to brain myoelectricity
Step research finds that there are two-way (downlink EEG → EMG, uplink EMG → EEG) couplings to contact between brain myoelectricity.But due to brain myoelectricity it
Between Function Coupling between coupling model is unknown and brain electromyography signal there is nonlinear causal relationship, the Glan based on set model
Outstanding causality analysis method cannot effectively describe brain myoelectricity Non-linear coupling feature.Transfer entropy has independent of set model and reality
The characteristics of existing quantitative analysis of nonlinear, it can effectively estimate the Function Coupling intensity and information transfer side between cortex-muscle
To.Therefore, transmitting entropy model is for the Function Coupling intensity and information direction of transfer feature, announcement between estimation cortex-muscle
Motion control and response mechanism in motion process between cortex and muscle have feasibility.2015, author of the present invention " had been based on
The brain electromyography signal coupling analysis of multiple dimensioned transfer entropy " in it is proposed that multiple dimensioned transfer entropy method, and not based on this method research
With brain electromyography signal coupling feature in time scale.But with the increase of coarse scale, sequence length reduces, may make
Entropy estimate inaccuracy.The it is proposed of mobile equalization overcomes this drawback, so that the length of time series of each scale keeps phase
Together.But the above research still has several drawbacks: coarse is only when carrying out brain electricity and electromyography signal with mobile equalization method
Between it is sized, can not depict brain electricity and myoelectricity time-frequency domain characteristic and different time-frequency scale diencephalon electromyography signal it is non-thread
Property coupling and information transmitting.
Summary of the invention
The purpose of the present invention is to provide one kind it can be found that non-linear dependencies between cortex muscle, further investigation brain skin
Layer muscle between coupling and information transfer characteristic based on the more time-frequency scale diencephalon myoelectricity coupling analysis of small echo-transfer entropy
Method.
The step of the method for the invention, is as follows:
Step 1, Neuroscan equipment synchronous acquisition EEG signals and electromyography signal are led in use 64;
Step 2, collected EEG signals and electromyography signal are gone respectively using Neuroscan device data processing software
Except baseline drift, spilling, eye movement and Hz noise;
Step 3, spectral decomposition, analysis are carried out to brain electricity and electromyography signal using Daubechies class db4 wavelet basis function
Synchronizing characteristics between EEG signals and electromyography signal difference time-frequency scale, Non-linear coupling and information transmitting between quantitative description brain flesh
Feature;
Step 4, motor function analysis is carried out between Non-linear coupling brain flesh and information transfer characteristic.
Further, in step 1, electrode for encephalograms uses world 10-20 system standard, using the mastoid process of ears as ginseng
It examines, the EEG signals of corresponding movement is recorded from 32 top guide skin brain wave acquisition equipment;Electromyography signal is acquired using Synamp2 equipment,
Electrode is placed at belly of muscle position along muscle fibre direction.
Further,
The specific method is as follows for the step 3:
Based on measured data pretreated in step 2 building EEG signals x (t) and two groups of time sequences of electromyography signal y (t)
Column;Wavelet transformation is carried out to EEG signals x (t), wavelet transformation is applied among partition of the scale;
Wavelet structure function, formula are as follows first:
Ψj,k(t)=2-j/2Ψ(2-jt-k) (1)
In formula, Ψ (t) is morther wavelet;Translational movement of the k for Ψ (t) ordinate direction, the number of plies of j expression signal, j, k ∈ Z,
Z is set of integers;Scale parameter is 2j, translation parameters 2jk;T is time index;
Then 7 layers of spectral decomposition are carried out to EEG signals x (t), obtains wavelet conversion coefficient
Wavelet coefficient Cj,kAccording to the sequence arrangement of frequency range from high to low, extract the 3rd, 4,5,6,7 layer coefficients reconstruct
Gamma (32~64Hz), beta (16~32Hz), alpha (8~16Hz), theta (4~8Hz) and delta (1~4Hz) frequency
The signal of section:
In formula,For the quality factor of filter;
Reconstruction signal x1(t)、x2(t)、x3(t)、x4(t) and x5(t) respectively correspond brain Electricity Functional frequency band delta, theta,
The signal of alpha, beta and gamma;
Above-mentioned identical wavelet transform procedure is carried out for electromyography signal y (t), obtains ym(t) (m=1,2,3,4,5), point
The signal of myoelectricity delta, theta, alpha, beta and gamma frequency range is not corresponded to;
Based on transfer entropy calculation method, small echo-transfer entropy WTE that x (t) arrives y (t) is constructedx→y, formula is as follows:
In formula, u is predicted time;Joint probability of the p () between variable;Respectively indicate brain electricity and myoelectricity
The delay vector of delta, theta, alpha, beta, gamma component;ForForecasting sequence;
WTEx→yThen indicate the x of EEG signals EEGi(t) y of the component to electromyography signal EMGm(t) transfer entropy between component
Value;Similarly signal y (t) arrives small echo-transfer entropy WTE of x (t)y→xExpression formula are as follows:
In formula,ForForecasting sequence;WTEy→xIndicate the y of EMGm(t) to component to the x of EEGi(t) between component
Transmit entropy;It is bigger to transmit entropy, illustrates that the coupling of cortex muscle is stronger between this frequency range;Vice versa.
The present invention has the beneficial effect that compared with prior art: the method for the present invention analyzes brain myoelectricity using small echo-transfer entropy
Signal message transmission characteristic, it is nonlinear between quantitative description EEG signals and electromyography signal frequency range to couple and information transfer characteristic,
Help to explore the functional cohesion between cerebral cortex and muscle, study movement controls feedback mechanism and dyskinesia pathology machine
System establishes the rehabilitation state evaluation index based on brain electromyography signal, constructs healing robot motion state and patient physiological condition
Evaluation mechanism can obtain considerable Social benefit and economic benefit.
Detailed description of the invention
Fig. 1 is the structure diagram of Neuroscan equipment in the present invention.
Fig. 2 is the work flow diagram of the method for the present invention.
Fig. 3 is time-frequency result figure after the wavelet decomposition of the channel subject C4 EEG signals.
Fig. 4 is time-frequency result figure after the wavelet decomposition of electromyography signal at subject's musculus flexor digitorum sublimis.
Fig. 5 is that brain electromyography signal small echo-transfer entropy of subject analyzes result figure.
Drawing reference numeral: 1- electrode for encephalograms, 2- electrode cap, 3- brain myoelectricity Acquisition Instrument, 4- myoelectricity conducting wire, 5- electromyographic electrode.
Specific embodiment
EEG signals and electromyography signal are very faint, have the characteristics that non-linear, non-stationary and frequency domain characteristic is prominent.?
In motion process, the interactive controlling mechanism between nervous system and muscle can pass through the synchronous coupling analysis body of brain electromyography signal
It is existing.Wavelet decomposition can extract the specific time-frequency data segment of brain electromyography signal, and transfer entropy can portray it is non-thread between signal
Property coupling and information transfer characteristic, the present invention by research brain myoelectricity between small echo-transfer entropy analysis, obtain different motion shape
Information transfering relation between cerebral cortex and muscle under state, and then the physiological mechanism that study movement dysfunction generates.
Embodiment 1:
As shown in Fig. 2, method and step is as follows:
Step 1, Neuroscan equipment synchronous acquisition EEG signals and electromyography signal are led in use 64.
The structure of Neuroscan equipment is as shown in Figure 1, by electrode for encephalograms, electrode cap, brain myoelectricity Acquisition Instrument, myoelectricity lead
Line, electromyographic electrode connection composition.
Eeg signal acquisition: electrode for encephalograms adopt international standards 10-20 electrode place standard, pass through electrode cap 2 realize brain
Electrode 1 is contacted with scalp.The experiment of brain electromyography signal synchronous acquisition is carried out under the output movement of hand static grip.M1, M2 are led
Mastoid process is hit exactly overhead as reference electrode, ground electrode arrangement after connection is connected respectively to left and right ear, is adopted from 32 top guide skin brain electricity
Collect the EEG signals of record corresponding movement in the selection area C3, C4 and CPZ in equipment.
Electromyographic signal collection: musculus flexor digitorum sublimis (flexor digitorum is acquired using Synamp2 equipment
Superficialis, FDS) at electromyography signal, be tested the skin surface at position with alcohol wipe first, remove skin surface
Then electromyographic electrode 5 is pasted at belly of muscle position by grease and scurf along muscle fibre direction, and myoelectricity conducting wire 4 is appropriate
It is fixed to reduce the interference that conducting wire shakes in motion process to the greatest extent.
Step 2, collected EEG signals and electromyography signal are gone respectively using Neuroscan device data processing software
Except baseline drift, spilling, eye movement and Hz noise;
Step 3, spectral decomposition is carried out to brain electricity and electromyography signal using Daubechies class db4 wavelet basis function, be based on
Pretreated measured data building EEG signals x (t) and two groups of time serieses of electromyography signal y (t) in step 2.With brain telecommunications
For number x (t), wavelet transformation is applied among partition of the scale.Wavelet structure function, formula are as follows first:
Ψj,k(t)=2-j/2Ψ(2-jt-k) (1)
In formula, Ψ (t) is morther wavelet;Translational movement of the k for Ψ (t) ordinate direction, the number of plies of j expression signal, j, k ∈ Z,
Z is set of integers;Scale parameter is 2j, translation parameters 2jk;T is time index.
Then 7 layers of spectral decomposition are carried out to EEG signals x (t), obtains wavelet conversion coefficient
The wavelet coefficient C obtained based on above formulaj,kAccording to frequency range from high to low sequence arrangement, extract the 3rd, 4,5,
6,7 layer coefficients reconstruct gamma (32~64Hz), beta (16~32Hz), alpha (8~16Hz), theta (4~8Hz) and
The signal of delta (1~4Hz) frequency range:
In formula,For the quality factor of filter;
Reconstruction signal x1(t)、x2(t)、x3(t)、x4(t) and x5(t) respectively correspond brain Electricity Functional frequency band delta, theta,
The signal of alpha, beta and gamma.
Above-mentioned identical wavelet transform procedure is carried out for electromyography signal y (t), obtains ym(t) (m=1,2,3,4,5), point
The signal of myoelectricity delta, theta, alpha, beta and gamma frequency range is not corresponded to.
Based on transfer entropy calculation method, small echo-transfer entropy WTE that x (t) arrives y (t) is constructedx→y, formula is as follows:
In formula, u is predicted time;Joint probability of the p () between variable;Respectively indicate brain electricity and myoelectricity
The delay vector of delta, theta, alpha, beta, gamma component;ForForecasting sequence;
WTEx→yThen indicate the x of EEG signals EEGi(t) y of the component to electromyography signal EMGm(t) transfer entropy between component
Value.Similarly signal y (t) arrives small echo-transfer entropy WTE of x (t)y→xExpression formula are as follows:
In formula,ForForecasting sequence;WTEy→xIndicate the y of EMGm(t) to component to the x of EEGi(t) between component
Transmit entropy;It is bigger to transmit entropy, illustrates that the coupling of cortex muscle is stronger between this frequency range;Vice versa.
It based on These parameters, calculates under the output movement of hand static grip, on difference coupling direction, between different time-frequency scales
WTE value, can non-linear synchronous coupling between quantitative description EEG signals EEG and the more time-frequency scales of electromyography signal EMG it is special
Sign.
For the feasibility and validity for verifying brain electromyography signal small echo-transfer entropy analysis method of the present invention, 8 are raised
The subject of name health carries out the output experiment of hand static grip, and subject's relevant information is as shown in table 1.According to of the present invention
Brain myoelectricity acquisition and analytic process, synchronous acquisition subject's constant force output movement under brain electromyography signal, analyzed and ground
Study carefully the coupling in subject motion's motion process between cortex muscle and information transmission mechanism.
Myoelectricity at this experiment acquisition left hand musculus flexor digitorum sublimis (flexor digitorum superficialis, FDS)
Signal and the channel opposite side C4 EEG signals, and calculate WTE value.
Fig. 3 and Fig. 4 is respectively subject's brain electricity, time-frequency domain knot of the electromyography signal after wavelet decomposition on each time-frequency scale
Fruit (left side is time-domain diagram, and right side is frequency domain figure), it can be seen that brain electricity and electromyography signal can obtain after wavelet decomposition
To the signal of delta, theta, alpha, beta and gamma frequency range.
Average value after the small echo that Fig. 5 is coupled between 8 subject's brain fleshes-transfer entropy analysis.There it can be seen that quiet
The stiffness of coupling of beta frequency range is the most significant between cortex muscle during the output of state grip, and each time-frequency on different coupling directions
There is also differences for stiffness of coupling between scale, provide theoretical study method to probe into neuromuscular function coupling mechanism.
1 subject's relevant information of table.
Claims (1)
1. a kind of based on the more time-frequency scale diencephalon myoelectricity coupling analytical methods of small echo-transfer entropy, which is characterized in that the method
The step of it is as follows:
Step 1, Neuroscan equipment synchronous acquisition EEG signals and electromyography signal are led in use 64;
Step 2, base is removed respectively to collected EEG signals and electromyography signal using Neuroscan device data processing software
Line drift, spilling, eye movement and Hz noise, for pretreated measured data building EEG signals x (t) and electromyography signal y
(t) two groups of time serieses;
Step 3, spectral decomposition is carried out to brain electricity and electromyography signal using Daubechies class db4 wavelet basis function, analyzes brain
Synchronizing characteristics between electric signal and electromyography signal difference time-frequency scale, Non-linear coupling and information transmitting are special between quantitative description brain flesh
Sign;Specific steps include:
1. EEG signals x (t) carries out wavelet transformation, wavelet transformation is applied among partition of the scale;Wavelet structure function first,
Formula is as follows:
In formula, Ψ (t) is morther wavelet;K is the translational movement of Ψ (t) ordinate direction, and j indicates the number of plies of signal, and j, k ∈ Z, Z are
Set of integers;Scale parameter is 2j, translation parameters 2jk;T is time index;
2. then carrying out 7 layers of spectral decomposition to EEG signals x (t), wavelet conversion coefficient is obtained
Wavelet coefficient Cj,kAccording to the sequence arrangement of frequency range from high to low, extract the 3rd, 4,5,6,7 layer coefficients reconstruct
The letter of gamma:32~64Hz, beta:16~32Hz, alpha:8~16Hz, theta:4~8Hz and delta:1~4Hz frequency range
Number:
In formula,For the quality factor of filter;
3. reconstruction signal x1(t)、x2(t)、x3(t)、x4(t) and x5(t) respectively correspond brain Electricity Functional frequency band delta, theta,
The signal of alpha, beta and gamma;Above-mentioned identical wavelet transform procedure is carried out for electromyography signal y (t), obtains ym(t),
M=1,2,3,4,5, respectively correspond the signal of myoelectricity delta, theta, alpha, beta and gamma frequency range;
4. being based on transfer entropy calculation method, construction x (t) arrives small echo-transfer entropy WTE of y (t)x→y, formula is as follows:
In formula, u is predicted time;Joint probability of the p () between variable;Respectively indicate brain electricity and myoelectricity delta,
The delay vector of theta, alpha, beta, gamma component;ForForecasting sequence;
WTEx→yThen indicate the x of EEG signals EEGi(t) y of the component to electromyography signal EMGm(t) the transmitting entropy between component;Similarly
Signal y (t) arrives small echo-transfer entropy WTE of x (t)y→xExpression formula are as follows:
In formula,ForForecasting sequence;WTEy→xIndicate the y of EMGm(t) to component to the x of EEGi(t) transfer entropy between component
Value;It is bigger to transmit entropy, illustrates that the coupling of cortex muscle is stronger between this frequency range;It is smaller to transmit entropy, illustrates in this frequency range mesothelium
Layer muscle coupling is weaker;
Step 4, Non-linear coupling the brain flesh under motion state and information transfer characteristic are analyzed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610362111.6A CN106073702B (en) | 2016-05-27 | 2016-05-27 | Based on the more time-frequency scale diencephalon myoelectricity coupling analytical methods of small echo-transfer entropy |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610362111.6A CN106073702B (en) | 2016-05-27 | 2016-05-27 | Based on the more time-frequency scale diencephalon myoelectricity coupling analytical methods of small echo-transfer entropy |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106073702A CN106073702A (en) | 2016-11-09 |
CN106073702B true CN106073702B (en) | 2019-05-28 |
Family
ID=57229890
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610362111.6A Active CN106073702B (en) | 2016-05-27 | 2016-05-27 | Based on the more time-frequency scale diencephalon myoelectricity coupling analytical methods of small echo-transfer entropy |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106073702B (en) |
Families Citing this family (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106901728B (en) * | 2017-02-10 | 2019-07-02 | 杭州电子科技大学 | Multichannel brain myoelectricity coupling analytical method based on mutative scale symbol transfer entropy |
CN106874589A (en) * | 2017-02-10 | 2017-06-20 | 泉州装备制造研究所 | A kind of alarm root finding method based on data-driven |
CN108733921B (en) * | 2018-05-18 | 2021-07-16 | 山东大学 | Transformer winding hot spot temperature fluctuation range prediction method based on fuzzy information granulation |
CN108742613A (en) * | 2018-05-30 | 2018-11-06 | 杭州电子科技大学 | Orient coupling analytical method between the flesh of coherence partially based on transfer entropy and broad sense |
CN109088770B (en) * | 2018-08-21 | 2020-03-31 | 西安交通大学 | Electromechanical system interactive network modeling method based on self-adaptive symbol transfer entropy |
CN109674445B (en) * | 2018-11-06 | 2021-10-08 | 杭州电子科技大学 | Inter-muscle coupling analysis method combining non-negative matrix factorization and complex network |
CN109497999A (en) * | 2018-12-20 | 2019-03-22 | 杭州电子科技大学 | Brain electromyography signal time-frequency coupling analytical method based on Copula-GC |
CN109657651A (en) * | 2019-01-16 | 2019-04-19 | 杭州电子科技大学 | A kind of continuous method for estimating of lower limb knee joint based on electromyography signal |
CN110251088A (en) * | 2019-07-03 | 2019-09-20 | 博睿康科技(常州)股份有限公司 | A kind of functional electric stimulation rehabilitation equipment for combining myoelectricity brain electricity |
CN110367974B (en) * | 2019-07-10 | 2022-10-28 | 南京邮电大学 | Brain and muscle electric coupling research method based on variational modal decomposition-transfer entropy |
CN110680315A (en) * | 2019-10-21 | 2020-01-14 | 西安交通大学 | Electroencephalogram and electromyogram signal monitoring method based on asymmetric multi-fractal detrending correlation analysis |
CN111067514B (en) * | 2020-01-08 | 2021-05-28 | 燕山大学 | Multi-channel electroencephalogram coupling analysis method based on multi-scale multivariable transfer entropy |
CN111227830B (en) * | 2020-02-14 | 2021-06-29 | 燕山大学 | Electroencephalogram and electromyographic coupling analysis method based on complex improved multi-scale transfer entropy |
CN111814390B (en) * | 2020-06-18 | 2023-07-28 | 三峡大学 | Voltage transformer error prediction method based on transfer entropy and wavelet neural network |
CN111904428A (en) * | 2020-06-29 | 2020-11-10 | 西安交通大学 | Electroencephalogram and electromyogram correlation analysis method for fine gait phase |
CN112541415B (en) * | 2020-12-02 | 2024-02-02 | 杭州电子科技大学 | Brain muscle function network motion fatigue detection method based on symbol transfer entropy and graph theory |
CN113197585B (en) * | 2021-04-01 | 2022-02-18 | 燕山大学 | Neuromuscular information interaction model construction and parameter identification optimization method |
CN113229831B (en) * | 2021-05-10 | 2022-02-01 | 燕山大学 | Movement function monitoring and management method based on myoelectricity and myooxygen signals |
CN113576403A (en) * | 2021-07-07 | 2021-11-02 | 南方科技大学 | Quantitative evaluation method for human body bidirectional coupling information conduction path and sensing system |
CN113408712A (en) * | 2021-07-16 | 2021-09-17 | 杭州电子科技大学 | Brain muscle coupling method based on time delay scale long-short term memory network and transfer entropy |
CN114041807A (en) * | 2021-12-20 | 2022-02-15 | 杭州电子科技大学 | Inter-muscle coupling analysis method based on wavelet packet-Copula transfer entropy |
CN114052750B (en) * | 2021-12-22 | 2024-04-30 | 杭州电子科技大学 | Brain muscle information transfer rule extraction method based on standard template myoelectricity decomposition |
CN114052751A (en) * | 2021-12-22 | 2022-02-18 | 杭州电子科技大学 | Movement function cortical muscle coupling method based on brain and muscle electricity |
CN115474945B (en) * | 2022-09-15 | 2024-04-12 | 燕山大学 | Multi-channel brain myoelectricity coupling analysis-oriented multi-element global synchronization index method |
CN116035597A (en) * | 2023-02-03 | 2023-05-02 | 首都医科大学宣武医院 | Electroencephalogram signal coupling analysis method, device and system |
CN116269392B (en) * | 2023-05-22 | 2023-07-18 | 华南理工大学 | Multi-parameter coupling stress level assessment method and system |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012094621A2 (en) * | 2011-01-06 | 2012-07-12 | The Johns Hopkins University | Seizure detection device and systems |
EP2535000A1 (en) * | 2011-06-17 | 2012-12-19 | Technische Universität München | Method and system for quantifying anaesthesia or a state of vigilance |
-
2016
- 2016-05-27 CN CN201610362111.6A patent/CN106073702B/en active Active
Non-Patent Citations (2)
Title |
---|
"Information transfer at multiple scales;Max Lungarella and Alex Pitti;《Physical Review》;20071230;第056117-1页第2栏第2段 |
基于多尺度传递熵的脑肌电信号耦合分析;谢平、杨芳梅等;《物理学报》;20151230(第24期);摘要,第248702-1页-248702-8页 |
Also Published As
Publication number | Publication date |
---|---|
CN106073702A (en) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106073702B (en) | Based on the more time-frequency scale diencephalon myoelectricity coupling analytical methods of small echo-transfer entropy | |
Chen et al. | A novel EEMD-CCA approach to removing muscle artifacts for pervasive EEG | |
CN110495880B (en) | Movement disorder cortical plasticity management method based on transcranial electrical stimulation brain muscle coupling | |
CN110238863B (en) | Lower limb rehabilitation robot control method and system based on electroencephalogram-electromyogram signals | |
CN104173124B (en) | A kind of upper limb healing system based on bio signal | |
CN103845137B (en) | Based on the robot control method of stable state vision inducting brain-machine interface | |
CN111227830B (en) | Electroencephalogram and electromyographic coupling analysis method based on complex improved multi-scale transfer entropy | |
Molla et al. | Artifact suppression from EEG signals using data adaptive time domain filtering | |
CN110367974B (en) | Brain and muscle electric coupling research method based on variational modal decomposition-transfer entropy | |
CN109497999A (en) | Brain electromyography signal time-frequency coupling analytical method based on Copula-GC | |
CN107822629B (en) | Method for detecting myoelectric axes on surfaces of limbs | |
CN104914994A (en) | Aircraft control system and fight control method based on steady-state visual evoked potential | |
CN102178524B (en) | Synchronization likehood-based electroencephalograph and electromyography synergistic analyzing method | |
CN103584855B (en) | Electroencephalogram and electromyogram synchronous acquisition and information transfer characteristic analysis method | |
Veer et al. | Wavelet and short-time Fourier transform comparison-based analysis of myoelectric signals | |
Huang et al. | Application and contrast in brain-computer interface between Hilbert-Huang transform and wavelet transform | |
Wang et al. | Feature extraction of brain-computer interface based on improved multivariate adaptive autoregressive models | |
CN108742613A (en) | Orient coupling analytical method between the flesh of coherence partially based on transfer entropy and broad sense | |
CN112137616B (en) | Consciousness detection device for multi-sense brain-body combined stimulation | |
CN112541415B (en) | Brain muscle function network motion fatigue detection method based on symbol transfer entropy and graph theory | |
CN105286860A (en) | Motor imagery brain electrical signal recognition method based on dual-tree complex wavelet energy difference | |
Pavlov et al. | Recognition of electroencephalographic patterns related to human movements or mental intentions with multiresolution analysis | |
Guerrero-Mendez et al. | Coherence-based connectivity analysis of EEG and EMG signals during reach-to-grasp movement involving two weights | |
CN103300850A (en) | Method for collecting and processing EEG (Electroencephalogram) signals of stroke patient | |
CN114601476A (en) | EEG signal emotion recognition method based on video stimulation |
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 |