CN102579040A - Acupuncture point selecting method for acquiring brain waves - Google Patents

Acupuncture point selecting method for acquiring brain waves Download PDF

Info

Publication number
CN102579040A
CN102579040A CN2012100710656A CN201210071065A CN102579040A CN 102579040 A CN102579040 A CN 102579040A CN 2012100710656 A CN2012100710656 A CN 2012100710656A CN 201210071065 A CN201210071065 A CN 201210071065A CN 102579040 A CN102579040 A CN 102579040A
Authority
CN
China
Prior art keywords
passage
acupuncture point
max
brain
information
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.)
Pending
Application number
CN2012100710656A
Other languages
Chinese (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN2012100710656A priority Critical patent/CN102579040A/en
Publication of CN102579040A publication Critical patent/CN102579040A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

The invention discloses an acupuncture point selecting method for acquiring brain waves. A correct acquisition channel is obtained finally. A plurality of brain acupuncture point signals are obtained via electroencephalogram equipment and are analyzed by the aid of a phase space reconstruction method, then embedding dimensions of the various acupuncture point signals are compared, a channel containing more information and a relatively independent channel are used as follow-up researched objects, information of other channels can be reconfigured from a channel with a large embedding dimension, and complexity of acquisition of the brain wave acupuncture point signals is simplified. Real-time analysis of the brain wave signals is strengthened. A mathematical model is simple, a physical significance is clear and definite, and the acupuncture point selecting method is easy to implement.

Description

The acupuncture point system of selection that a kind of brain wave is gathered
Technical field
The present invention relates to the electron medicine field, relate in particular to a kind of eeg signal method of sampling.
Background technology
The brain electricity be cranial nerve cell cerebral cortex and with the bioelectrical activity performance of scalp surface, it can reflect the functional status of brain physiology, pathological condition and brain effectively.Different thought state and emotion changes can both reflect different EEG signals in different cerebral cortex positions, so EEG signals contain abundant useful information.In biomedicine, the brain electricity is used as the effective means of medical diagnosis and disease treatment; In Cognitive Study, the brain electricity is as the important tool of human thinking's origin.That at first the discovery mankind have the brain electrical acti characteristic is the sick scholar Hans of Germanism Berger, and he has delivered the paper of a first piece of writing about human brain electricity in nineteen twenty-nine.
Traditional method is mainly from time domain, frequency spectrum and statistical angle research brain.Time-domain analysis mainly is waveform, the amplitude of research brain electricity, characteristics such as phase place, and frequency-domain analysis method is to analyze the distribution of brain wave on different frequency range.They can describe the Partial Feature of brain electricity, but poor sensitivity, and can't analyze the substitutive characteristics of cerebral activity.Methods such as wavelet transformation, neutral net and nonlinear kinetics also begin to be applied to the analysis of EEG signals in recent years, and they have represented the new direction of electroencephalogramsignal signal analyzing.
Become a very important field with chaology research brain electricity, brain is considered to a nonlinear dynamic system.Divide from the research purpose, mainly contain following two fields: the working mechanism of brain is disclosed in (1), and researching human body is in the Nonlinear Dynamical Characteristics under different physiological statuss, the different cerebral functional status.(2) researching human body is in the variation of nonlinear kinetics under the pathological state, is the clinical foundation of analyzing and diagnosing that provides.The phase space reconfiguration of Whitney delay reconstruct theorem theoretical and Takens provides from the theoretical basis of experimental period sequence reconfiguration system phase space.Based on this theory, Nonlinear Time Series is analyzed and studied, can very effective description goal systems, this also makes the research to the complication system of the unknown become possibility.
When non-linear sequence is carried out phase space reconfiguration; Two major parameters that influence the reconstruct effect are time delay and embedding dimension; Total principle is the nonlinear characteristic that guarantees that the sampling time sequence can be reappeared original system to greatest extent, keeps the information of system as much as possible.Nonlinear Dynamics is the new approaches of research EEG signals, all will play enormous function to deepening the mankind to understanding and disease detection, diagnosis and the treatment of brain.
" Applied Mathematics and Computation " (2009; 207 (1): 63-74) nonlinear dynamic characteristic of EEG signals has been discussed; " IEEE Transactions on Biomedical Engineering " (2011,58 (4): 1084-1093) utilization phase space reconfiguration theoretical research the brain signal situation of change of people when being in different condition.They handle the EEG signals of having gathered; But the collection to EEG signals is but also not mentioned; And we think; The collection of EEG signals is of paramount importance links of EEG Processing, if the signal that collects can comprise maximum quantity of information, will facilitate for subsequent treatment so.This mainly is because brain is not a single node, and the difference of the signal that collects at the brain different parts is very big, how therefrom to select a road optimum signal, makes its can maximum ground characterize the state of brain, comprises maximum quantity of information.Traditional method is generally chosen the signal of head acupoint as candidate signal, but can't make optimized choice to the acupuncture point signal of being gathered.
Summary of the invention
In order to overcome prior art when selecting brain to observe the acupuncture point; The imperfect complexity of bringing with multi-channel measurement of brain wave information; The present invention proposes a kind of new brain wave and observe the method at acupuncture point; Can when guaranteeing that the brain wave information gathering is complete, simplify the complexity that brain wave is gathered, for the plenty of time has been practiced thrift in follow-up brain wave real-time analysis.
The technical solution adopted for the present invention to solve the technical problems may further comprise the steps:
At first, obtain eeg signal, promptly gather the eeg signal of human brain; 32 acupuncture points with brain are acquisition target, make up the acupuncture point passage respectively, through setting-up time length, obtain the eeg signal of 32 acupuncture point passages; The data of 32 acupuncture point passages are stored respectively, obtain the array of 32 equal lengths;
Then, utilize the phase space reconfiguration theory respectively each array to be studied as an independent systematic component, confirm the embedding dimension m and the delay time T of system reconfiguration; At first adopt the FNN method (" Phys.Rev.A ", 1992,45:3403-3411) confirm to embed dimension m, adopt then the AD method (" Physica D ", first extreme point 1994:73:82-98) is as the numerical value of delay time T;
Secondly, the analysis one by one through to each array can obtain m 1... m i..m 32, m wherein iThe embedding dimension of representing the data of i degree passage; Through to embedding the screening of relatively studying passage of dimension: if m Max>32 (m MaxRepresent that the maximum in all passages embeds dimension), increase acquisition target, promptly increase and gather the acupuncture point number, make it equal m MaxIf, m Max≤32, then choose and have m MaxPassage be object of study, keep m simultaneously i=1 passage; Comprise m MaxPassage be to comprise the maximum passage of information, m i=1 passage is the relatively independent passage of E.E.G information;
At last, obtain brain wave reconstruction model information: with m MaxFor embedding dimension, be time delay with its corresponding τ, the reconfiguration system model:
X(0)=[x(0),x(τ),x(2τ),L,x((m-1)τ)]
X(1)=[x(1),x(τ+1),x(2τ+1),L,x((m-1)τ+1)]
...................................
X(n)=[x(n),x(n+τ),x(n+2τ),L,x(n+(m-1)τ)]
...............................................................
X(m-1)=[x(m-1),x(m-1+τ),x(m-1+2τ),L,x((m-1)(τ+1))]
This moment X (0) ... all data of .X (m-1) all are from having the maximum dimension m that embeds MaxPassage, but but can reflect the information of other passage, the information of X (0) expression 1 passage for example, X (1) represent 2 channel informations, X (n) representes the n+1 channel information.Therefore comprise m MaxPassage comprise except that m iThe relevant information of other all passages beyond=1 passage, subsequent analysis only need be gathered m MaxWith the m=1 channel data, need not data acquisition is being carried out at all acupuncture points.
The invention has the beneficial effects as follows: the present invention is from the collection of eeg signal; Utilized the theory of phase space reconfiguration; Embedding dimension through to each channel signal reconstruct compares; Selected to comprise relatively independently passage of maximum passage of information and information,, can from comprise the maximum passage of information, restore the information of other passage according to the phase space reconfiguration theory.Simultaneously, the present invention proposes the unified algorithm of phase space reconfiguration, when guaranteeing computational accuracy, improved the speed of reconstruct greatly, for the real-time analysis of eeg signal provides assurance; Mathematical model of the present invention is simple, explicit physical meaning, and it is accurate to be easy to realization and to calculate effect.
Below in conjunction with accompanying drawing and embodiment the present invention is further specified.
Description of drawings
The brain wave observation acupuncture point figure of Fig. 1 initial selected;
Fig. 2 is a phase space reconfiguration brain wave analysis operation flow chart of the present invention;
The specific embodiment
The present invention will find out quantitative foundation, gathers acupoint information at first in a large number, has proposed a kind of standard then and has weighed the size of each acupoint information amount, and finally selects one-channel signal and carry out processed steps.
Present embodiment is according to the operating process shown in the accompanying drawing 2, and a simple embodiment of gathering based on the human brain electric wave signal shown in the accompanying drawing 1.
At first can the data that obtain be transferred in the PC through the usb data line, in PC, the transmission type of data be changed, according to acquisition channel storage { x (t through the signal at initially selected 32 acupuncture points of eeg signal acquisition equipment collection i).
Compare the embedding dimension m value of each passage then, obtain maximum embedding dimension value m MaxIf, m Max>32, increase acquisition target, promptly increase and gather the acupuncture point number, make it equal m MaxIf, m Max≤32, m then MaxThe contained quantity of information of corresponding passage is maximum, and to choose this road be final preferred passage, keeps the passage of m=1 simultaneously; m MaxBe to comprise the maximum passage of information, m=1 is the relatively independent passage of E.E.G information;
At last, then choose and have the maximum dimension m that embeds MaxPassage be object of study, keep to embed the passage of dimension m=1 simultaneously; Because m MaxBe to comprise the maximum passage of information, m=1 is the relatively independent passage of E.E.G information
Proper collection acupuncture point during the final human brain signals collecting that obtains of the present invention; Simplify the complexity of signals collecting through the phase space reconfiguration theory, strengthened the effect of eeg signal real-time analysis, can reconstruct whole brain signal output accurately; The mathematical model of system reconfiguration is simple; Explicit physical meaning is easy to realize that effect is accurate.

Claims (1)

1. the acupuncture point system of selection that brain wave is gathered is characterized in that comprising the steps:
At first, be acquisition target with 32 acupuncture points of brain, make up the acupuncture point passage respectively, through setting-up time length, obtain the eeg signal of 32 acupuncture point passages; The data of 32 acupuncture point passages are stored respectively, obtain the array of 32 equal lengths;
Then, utilize the phase space reconfiguration theory respectively each array to be studied as an independent systematic component, adopt the FNN method to confirm to embed dimension m, adopt the numerical value of first extreme point of AD method as delay time T;
Secondly, each array is analyzed one by one, obtained m 1... m i..m 32, m wherein iThe embedding dimension of representing the data of i degree passage; Through to embedding the screening of relatively studying passage of dimension: if m Max>32, m MaxRepresent that the maximum in all passages embeds dimension, increase and gather the acupuncture point number, make it equal m MaxIf, m Max≤32, then choose and have m MaxPassage be object of study, keep m simultaneously i=1 passage; Comprise m MaxPassage be to comprise the maximum passage of information, m i=1 passage is the relatively independent passage of E.E.G information;
At last, obtain brain wave reconstruction model information: with m MaxFor embedding dimension, be time delay with its corresponding τ, the reconfiguration system model:
X(0)=[x(0),x(τ),x(2τ),L,x((m-1)τ)]
X(1)=[x(1),x(τ+1),x(2τ+1),L,x((m-1)τ+1)]
................................
X(n)=[x(n),x(n+τ),x(n+2τ),L,x(n+(m-1)τ)]
...............................................................
X(m-1)=[x(m-1),x(m-1+τ),x(m-1+2τ),L,x((m-1)(τ+1))]
Subsequent analysis is only gathered m MaxWith the m=1 channel data.
CN2012100710656A 2012-03-18 2012-03-18 Acupuncture point selecting method for acquiring brain waves Pending CN102579040A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012100710656A CN102579040A (en) 2012-03-18 2012-03-18 Acupuncture point selecting method for acquiring brain waves

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012100710656A CN102579040A (en) 2012-03-18 2012-03-18 Acupuncture point selecting method for acquiring brain waves

Publications (1)

Publication Number Publication Date
CN102579040A true CN102579040A (en) 2012-07-18

Family

ID=46468685

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012100710656A Pending CN102579040A (en) 2012-03-18 2012-03-18 Acupuncture point selecting method for acquiring brain waves

Country Status (1)

Country Link
CN (1) CN102579040A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108594989A (en) * 2018-03-26 2018-09-28 广东欧珀移动通信有限公司 Acquiring brain waves method and relevant device
CN111387976A (en) * 2020-03-30 2020-07-10 西北工业大学 Cognitive load assessment method based on eye movement and electroencephalogram data

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
MATTHEW B. KENNEL ET AL.: "Determining embedding dimension for phase-space reconstruction using a geometrical construction", 《PHYSICAL REVIEW A》 *
张菁 等: "相空间重构中延迟时间选取的新算法", 《计算物理》 *
杨大鹏: "特定心理作业多通道脑电信号识别系统研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108594989A (en) * 2018-03-26 2018-09-28 广东欧珀移动通信有限公司 Acquiring brain waves method and relevant device
CN111387976A (en) * 2020-03-30 2020-07-10 西北工业大学 Cognitive load assessment method based on eye movement and electroencephalogram data
CN111387976B (en) * 2020-03-30 2022-11-29 西北工业大学 Cognitive load assessment method based on eye movement and electroencephalogram data

Similar Documents

Publication Publication Date Title
Hramov et al. Wavelets in neuroscience
CN101596101B (en) Method for determining fatigue state according to electroencephalogram
CN104720797B (en) One kind is based on myoelectricity noise cancellation method in single pass EEG signals
CN102488516A (en) Nonlinear electroencephalogram signal analysis method and device
Zawawi et al. Electromyography signal analysis using spectrogram
CN103610447A (en) Mental workload online detection method based on forehead electroencephalogram signals
CN105326499B (en) A kind of portable brain electric acquisition method
CN112137616B (en) Consciousness detection device for multi-sense brain-body combined stimulation
CN100998503A (en) Method for automatically recogniting and eliminating ophthalmogyric interference in electroencephalo-signals
CN105125186B (en) A kind of method and system of definite therapeutic intervention mode
CN113274037A (en) Method, system and equipment for generating dynamic brain function network
Zhang et al. An improved method to calculate phase locking value based on Hilbert–Huang transform and its application
CN101433460A (en) Spatial filtering method of lower limb imaginary action potential
CN102579040A (en) Acupuncture point selecting method for acquiring brain waves
CN111543983B (en) Electroencephalogram signal channel selection method based on neural network
Zhang An improved QRS wave group detection algorithm and matlab implementation
CN117860271A (en) Classifying method for motor imagery electroencephalogram signals
CN106679799A (en) Thunder signal generation system and thunder signal simulation method
Li et al. Subject-based dipole selection for decoding motor imagery tasks
Principe et al. Representing and decomposing neural potential signals
Zhou et al. A smart universal single-channel blind source separation method and applications
CN115299960A (en) Electric signal decomposition method and electroencephalogram signal decomposition device based on short-time varying separate modal decomposition
Banta et al. A novel convolutional neural network for reconstructing surface electrocardiograms from intracardiac electrograms and vice versa
Singh et al. EEG Signals: Current trends and Future Aspects
CN118303887B (en) Electrophysiological signal processing method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20120718