Background
Electroencephalography (EEG) is a non-invasive method of measuring voltage fluctuations due to ionic currents in neurons in the brain, which are directly reflected in the electrophysiological activity of brain neurons in the cerebral cortex. Electroencephalography is a non-linear time-domain representation of electroencephalogram, which is recorded by a plurality of electrodes placed on the scalp and contains a large amount of physiological and pathological information. The electroencephalogram signal is the comprehensive expression of a large number of neuron cluster postsynaptic currents in brain tissues on the surface of a cerebral cortex or a scalp, and can be regarded as the result of the superposition of different oscillation frequency components on different time scales. As an objective evaluation index of brain function, the electroencephalogram signal has very high time precision, can dynamically observe the state change of the brain, and provides a basis for real-time diagnosis and treatment of brain diseases. And the human brain health care product has rich contents including feelings, thinking, spirit and psychological activities, so that the human brain health care product becomes an important research method for high-level functions of human brain learning, memory, understanding and the like.
The frequency cross coupling is an important method for researching the electric potential of the electroencephalogram related event, and can be used for researching the mutual relation among signals with different frequencies. Frequency cross-coupling of synchronous phases describes the phase independence of low frequency signals from high frequency signals, which has some relationship to recognized cognitive processing mechanisms such as language, memory, etc. Brain diseases such as alzheimer's, epilepsy, and parkinson's cause impairment of cognitive function in the human brain. Therefore, the study of the nonlinear characteristics of the phase cross frequency of the electroencephalogram signals has important significance for the study and detection of brain diseases.
However, the existing nonlinear synchronization algorithm only surrounds the phase synchronization of the electroencephalogram signal under a certain specific frequency, neglects the phase synchronization coupling characteristic among frequencies, and if the algorithm is used for constructing a frequency cross synchronization network of the electroencephalogram signal, the calculated phase coupling has errors, and an accurate phase synchronization relation cannot be obtained.
Disclosure of Invention
The invention aims to provide an electroencephalogram signal feature extraction method which can obtain an accurate phase synchronism relation by considering the phase synchronism coupling characteristic among frequencies. The technical scheme is as follows:
an electroencephalogram signal feature extraction method based on a complex network comprises the following steps:
(1) wavelet packet decomposition
Collecting X-channel electroencephalogram signals through an electroencephalogram machine, carrying out eight-level wavelet packet decomposition on each channel electroencephalogram signal, decomposing the signals into a low-frequency part and a high-frequency part at the first level, decomposing the decomposed low-frequency part into a low-frequency part and a high-frequency part at the second level, decomposing the high-frequency part into a low-frequency part and a high-frequency part, and so on, averagely dividing the frequency of the original signals into 2 parts through the eighth-level wavelet packet decomposition8And (4) decomposing each lead brain electrical signal into Q sub-band brain electrical signals according to frequency band components.
(2) Calculating the frequency cross-synchronization of EEG signals
Two coefficients n and m are defined,
wherein f is
x,f
yThe center frequencies of the sub-bands where the electroencephalogram signals x and y are located are respectively, n and m respectively take the minimum positive integer which meets the proportional relation, and therefore the phase difference delta phi between the electroencephalogram signals of any two sub-bands is calculated to be n phi
x-mφ
yAnd calculating two EEG signal synchronism indexes by using the obtained phase difference delta phi, wherein the synchronism index range is between 0 and 1, the larger the value is, the higher the phase synchronism is, and a weighting matrix M containing the cross phase synchronism degree of the frequencies of the same sub-band and different sub-bands is obtained.
(3) Constructing a frequency crossover network
Performing proportional thresholding on the weighting matrix M obtained in thestep 2, wherein the matrix is subjected to threshold value T, T is more than 0 and less than 10, T represents an element which retains the maximum value of T × 10/100 in the matrix, so as to obtain a thresholded matrix, each element in the matrix represents the synchronism between two leads, and the leads represent nodes in the network; and constructing a frequency cross network according to the obtained thresholded matrix.
(4) Network feature parameter extraction
Extracting characteristic parameters of the frequency crossing network obtained in thestep 3, setting all node sets as G for the network with N nodes, and defining the length of the shortest node path of the network as:
lijdefined as the number of edges on the shortest path connecting two nodes i and j.
The local efficiency of the network is defined as:
global efficiency is the average of the local efficiency sums of all nodes:
calculating the frequency cross network obtained in thestep 3 according to the formula (2), the formula (3) and the formula (4) to obtain a phase
The corresponding network parameter characteristics are used for revealing the EEG frequency cross-coupling characteristics.
The method has the effects of quantifying the frequency band cross synchronization relation of the electroencephalogram signals, constructing a frequency cross network of a comprehensive frequency domain and a space domain, and revealing the frequency cross coupling characteristics of the electroencephalogram signals through extraction and analysis of network parameters.
Detailed Description
As shown in FIG. 1, the method for extracting the electroencephalogram signal features comprises the following steps:
1. wavelet packet decomposition
The original brain wave data is collected, and in the embodiment, 16 brain wave signals are collected for 20 Alzheimer patients and 20 age-matched normal persons. All experimenters lie on the special bed for electroencephalogram acquisition quietly, and eyes are closed in the whole process.
Wavelet packet decomposition reconstruction is carried out on the electroencephalogram signals by utilizing a one-dimensional wavelet packet decomposer (wavelet packet decomposition 1-D), and 4 sub-band physiological rhythms are extracted. The wpdec function in MATLAB 2015b used for data analysis is used for carrying out 8-level wavelet packet decomposition on an original signal, and the sub-frequency band is realized: delta (0.5-4Hz), theta (4-8Hz), alpha (8-16Hz) and beta (16-30Hz) rhythms. The wavelet reconstruction principle is shown in fig. 2.
1. Calculating the frequency cross-synchronization of EEG signals
For the extracted sub-band multi-lead electrical time sequence signal (4 sub-bands with 16 lead electrical data in each band). Firstly, the phase phi of each pilot sequence signal is calculated through Hilbert
changei1,2,3, 64. Calculating phase difference between brain electrical signals (n phi)
x-
mφy1,2,3, 64 and x ≠ y, wherein
f
x,f
yIs the center frequency of the sub-band where the electroencephalogram signals x and y are located. Then, the phase difference between the two electroencephalogram signals is converted into a phase synchronism index by formula (1), thereby obtaining a frequency cross adjacency matrix M, as shown in fig. 3. The matrix comprises the following parts: { M
ijI, j ═ 1,2,3.., 16} represents the delta in-band adjacency matrix; { M
ijI, j ═ 17,18,19.., 32} represents the θ in-band synchronization matrix; { M
ijI, j 33,34,35., 48} represents an α -in-band synchronization matrix; { M
ijI, j ═ 49,50,51.., 64} represents the β in-band synchronization matrix; {
Mij1,2,3, 16, j 17,18,19, 32 represents a frequency cross synchronization matrix of δ - θ; {
Mij1,2,3, 16, j 33,34,35, 48 represents a frequency cross synchronization matrix of δ - α; {
Mij1,2,3, 16, j 49,50,51, 64 represents a frequency cross synchronization matrix of δ - β; { M
ij17,18,19, 32, j 33,34,35, 48 represents a frequency cross-synchronization matrix of θ - α; { M
ij17,18,19, 32, j 49,50,51, 64 represents a frequency cross-synchronization matrix of θ - β; { M
ij33,34,35., 48, j 49,50,51., 64 represents a frequency cross-synchronization matrix of α - β;
3. constructing a frequency crossover network
And (3) taking a 0.2 proportional threshold value for the frequency cross adjacency matrix M obtained in the step (2) to obtain a two-dimensional threshold value matrix. The matrix is used to construct a frequency crossover network, as shown in fig. 4.
4. Network feature parameter extraction
And (4) extracting characteristic parameters of the frequency cross network obtained in the step (3), and setting all node sets as G.
The shortest path length of a node of the network is defined as:
lijdefined as the number of edges on the shortest path connecting two nodes i and j.
The local efficiency of the network is defined as:
global efficiency is the average of the local efficiency sums of all nodes:
and (3) calculating the local efficiency and the global efficiency of the frequency cross network according to the formula (2), the formula (3) and the formula (4), and finishing the feature extraction of the electroencephalogram signals. The extraction results are shown in fig. 5 and 6. FIG. 4 vertical axis is local efficiency E of the networkiThe abscissa represents two sub-bands corresponding to the frequency cross network, and the network with the delta band crossing other frequency bands has a significant increase in local efficiency compared with the normal control group under the influence of the alzheimer. FIG. 5 ordinate is the global efficiency E of the networkgThe abscissa is two frequency sub-bands corresponding to the frequency cross network, and in the alzheimer group, the global efficiency of all the frequency cross networks is increased relative to the normal control group. Under the influence of the Alzheimer, the extracted local efficiency and the global efficiency are obviously improved, so that the characteristic value can be used for distinguishing the Alzheimer patient from a healthy person.