CN102499675B - Feedback system random resonance intensification method of electro-corticogram signal - Google Patents

Feedback system random resonance intensification method of electro-corticogram signal Download PDF

Info

Publication number
CN102499675B
CN102499675B CN 201110331467 CN201110331467A CN102499675B CN 102499675 B CN102499675 B CN 102499675B CN 201110331467 CN201110331467 CN 201110331467 CN 201110331467 A CN201110331467 A CN 201110331467A CN 102499675 B CN102499675 B CN 102499675B
Authority
CN
China
Prior art keywords
signal
eeg signals
cortex
cortex eeg
neuron
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN 201110331467
Other languages
Chinese (zh)
Other versions
CN102499675A (en
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.)
Yunfeng Pipe Industry Co Ltd
Original Assignee
Hangzhou Dianzi 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 Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN 201110331467 priority Critical patent/CN102499675B/en
Publication of CN102499675A publication Critical patent/CN102499675A/en
Application granted granted Critical
Publication of CN102499675B publication Critical patent/CN102499675B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention relates to a feedback system random resonance intensification method of an electro-corticogram signal. Different from a traditional intensification method of the elctro-corticogram signal, the method is advantageous to improvement of the performance of intensifying the electro-corticogram signal through a random resonance mechanism under noise interference with certain strength. In consideration of the characteristics of electrophysiological signals such as the electro-corticogram signal, a FitzHugh-Nagumo model in accordance with the neuron electrophysiological characteristics is adopted, and a feedback element is introduced to a one-way neuron network structure, so that the method is more in accordance with the complicated interconnection relationship among neurons in a neuron system. Compared with the traditional noise filtering method, the method can better restore and intensify the electro-corticogram signal.

Description

A kind of feedback system accidental resonance Enhancement Method of cortex EEG signals
Technical field
The invention belongs to biomedical engineering field, relate to the processing method that a kind of cortex EEG signals strengthens, be specifically related to a kind of faint cortex EEG signals real time enhancing method based on neuron Feedback network model random resonance mechanism.
Background technology
Because the cortex EEG signals easily is subject to the interference of various electricity physiological signals and other noises, the cortex EEG signals that therefore detects is fainter, and Noise is larger, and this studies EEG signals for us and has a great impact.Effectively eliminate the adverse effect of noise, the EEG signals that remains with simultaneously usefulness becomes particularly important.Traditional noise for harmful is mainly taked the method for filtering, the priori of interfering signal, the cortex EEG signals signal to noise ratio that collects is lower, if removal noise, EEG signals can be subject to very large infringement, even EEG signals can be used as together filtering of noise, so that the EEG signals distortion can't be restored.Because random resonance mechanism can be in harmonious proportion relation between nonlinear system, signal and the noise, noise energy is shifted to signal energy, therefore random resonance mechanism thinks that the existence of noise is significant, the cortex EEG signals that suitable signal to noise ratio is lower in some sense.The nonlinear system of present existing accidental resonance derives from the physics abstract models such as bistable system usually, and it realizes the enhancing of EEG signals by EEG signals, noise information and the triangular accidental resonance of nonlinear system.But these physics abstract models are too idealized, and whether its electro-physiological signals to this quasi-representative of brain electricity is fit to, and does not have sufficient foundation.Therefore the present invention proposes to adopt FitzHugh-Nagumo (FHN) neuron models that can truly reflect the neuron electrophysiological characteristics, and forms the network-feedback structure with being connected to each other between the simulation cerebral neuron; To be subject to the input of the low signal-to-noise ratio cortex EEG signals of various interference as aforementioned feedback neural metanetwork model, utilize random resonance mechanism to realize the enhancing of weak signal.
Summary of the invention
The present invention is based on the FHN neuron models, a kind of random resonance mechanism based on the neuron Feedback network model is provided, realize the enhancing of faint cortex EEG signals with this.Feedback element has been added on basis at the FHN neural network model, avoids unidirectional control and the unstability of neural network model with this, and the realization output signal is improved the accidental resonance performance of FHN neural network model for the regulating action of input signal.
The inventive method may further comprise the steps:
Step (1) obtains the cortex EEG signals, and is divided into several nonoverlapping specific duration windows by the cortex eeg collection system.Cortex EEG signals in each window is asked for reference value, with the intermediate value between its maximum and minima, as the reference value of this cortex EEG Processing.
Step (2) is to the cortex EEG signals in each window, carry out respectively the bipolar processes of amplitude: the cortex EEG signals value in each window is deducted the reference value that step (1) is tried to achieve, acquisition has ambipolar cortex EEG signals, makes it satisfy the neuron models input signal and has ambipolar requirement.
The value of feedback of the cortex EEG signals that step (3) obtains step (2), the noise signal of interpolation and output signal is as the input signal based on FHN neuron Feedback network model, utilize the accidental resonance effect of FHN neuron Feedback network model, obtain the cortex brain electrical output signal that strengthens.
Step (4) under the varying strength that adds noise signal, is calculated respectively the signal to noise ratio of FHN neuron Feedback network model response in the step (3).Utilize random resonance mechanism, when signal to noise ratio reached maximum, the optimum that the cortex EEG signals will obtain under the signal to noise ratio evaluation index this moment strengthened, with this cortex EEG signals as output signal.
The reference value that step (5) is asked for above-mentioned output signal and step (1) is sued for peace, and former amplitude scope is returned in inverse mapping, thereby obtains the cortex EEG signals of the enhancing after the amplitude reduction.
Beneficial effect of the present invention:
1, because the cortex EEG signals has transient characteristic, the cortex EEG signals that the present invention arrives dynamic acquisition, its intermediate value is asked in setting in short-term long window, obtains the reference value that follow-up bipolarity mapping is processed, this reference value has dynamic characteristic, and the real-time processing that utilizes the transition EEG signals is arranged.
2, the present invention is based on the random resonance mechanism of FHN neuron Feedback network model, different from traditional signal Enhancement Method based on noise filtering, its passive power conversion with noise is the positive energy of signal, thereby realizes the enhancing of faint cortex EEG signals.
3, the present invention has given up the abstract models such as bistable system commonly used in the accidental resonance of EEG signals strengthens, but adopts the neuron models that meet true neuron electrophysiological characteristics; Network models in the unidirectional connection of neuron has increased feedback element simultaneously, more meets the interconnected relationship between the neuron in the nervous system, is conducive to improve the stability in the faint cortex EEG signals enhancing process.
Description of drawings
Fig. 1 is FHN neuron feedback double-layer network model structure schematic diagram.
The specific embodiment
Step (1) is used cortex eeg signal acquisition system, gathers one section continuous cortex EEG signals, and it is divided into several nonoverlapping windows, and the window duration is designated as
Figure DEST_PATH_IMAGE002
, wherein NSampling number in the expression window, TThe expression sampling period.Therefore the cortex EEG signals in the window can be designated as
Figure DEST_PATH_IMAGE004
, (
Figure DEST_PATH_IMAGE006
).It is asked for maximum
Figure DEST_PATH_IMAGE008
And minima
Figure DEST_PATH_IMAGE010
, their average as the dynamic benchmark value of this window EEG signals, is designated as
Figure DEST_PATH_IMAGE012
Step (2) is carried out the bipolar processes of amplitude with the cortex EEG signals in each window of step (1).Each cortex EEG signals sampled value that is about in the window deducts respectively the reference value that step (1) is tried to achieve, and obtains to have ambipolar cortex EEG signals
Figure DEST_PATH_IMAGE014
, (
Figure 178967DEST_PATH_IMAGE006
), make it satisfy the input signal bipolarity characteristics of neuron models.
Cortex EEG signals, the white Gaussian noise of interpolation and the value of feedback of output signal that step (3) obtains step (2) are as the input signal of FHN neuron Feedback network model.
Be specifically described as an example of bilayer feedback FHN neural network model example, wherein the model structure schematic diagram as shown in Figure 1, among the figure
Figure DEST_PATH_IMAGE016
Be the cortex EEG signals of current time after the bipolarity mapping is processed;
Figure DEST_PATH_IMAGE018
Figure DEST_PATH_IMAGE020
It is the independent noise item of same noise intensity;
Figure DEST_PATH_IMAGE022
Figure 846840DEST_PATH_IMAGE020
Be the feedback regulation parameter, wherein network system is negative feedback;
Figure DEST_PATH_IMAGE024
Figure 564260DEST_PATH_IMAGE020
The of expression ground floor
Figure DEST_PATH_IMAGE026
The neuronic coefficient of connection of individual neuron and the second layer;
Figure DEST_PATH_IMAGE028
Figure 831906DEST_PATH_IMAGE020
Be of ground floor
Figure 583961DEST_PATH_IMAGE026
Individual FHN neuron;
Figure DEST_PATH_IMAGE030
Be second layer FHN neuron;
Figure DEST_PATH_IMAGE032
Be the neuronic output membrane voltage of the second layer, i.e. output signal.
Mathematical model is suc as formula shown in (1) and the formula (2):
Ground floor:
Figure DEST_PATH_IMAGE034
(1)
The second layer:
Figure DEST_PATH_IMAGE036
(2)
In formula (1) and the formula (2), Be the sampling period;
Figure DEST_PATH_IMAGE040
Be output signal;
Figure DEST_PATH_IMAGE042
Be time constant, determined neuronic ignition rate;
Figure DEST_PATH_IMAGE044
Be marginal value, impel neuron regularly to light a fire;
Figure DEST_PATH_IMAGE046
For the signal level average with Difference;
Figure DEST_PATH_IMAGE048
,
Figure DEST_PATH_IMAGE050
Be the equation group constant;
Figure DEST_PATH_IMAGE052
Be the input signal of bilayer feedback FHN neural network model,
Figure DEST_PATH_IMAGE054
,
Figure DEST_PATH_IMAGE056
Be the cortex EEG signals after the bipolarity mapping processing, The noise item in the input signal, usually by average be 0, auto-correlation function is
Figure DEST_PATH_IMAGE060
White Gaussian noise simulated, wherein
Figure DEST_PATH_IMAGE062
Be noise intensity,
Figure DEST_PATH_IMAGE064
The expression impulse function;
Figure 2619DEST_PATH_IMAGE022
Figure 139203DEST_PATH_IMAGE020
Be the feedback regulation parameter; Be the first floor iIndividual neuronic output membrane voltage;
Figure 676811DEST_PATH_IMAGE032
Be the neuronic output membrane voltage of the second layer, i.e. output signal;
Figure DEST_PATH_IMAGE068
Slowly become the recovery variable for the second layer is neuronic; Coefficient of connection is
Utilize the value of feedback of the noise signal of cortex EEG signals to be strengthened, interpolation and response as excitation, random resonance mechanism by FHN neuron Feedback network model, thereby realize the enhancing of faint cortex EEG signals, obtain the cortex brain electrical output signal that strengthens.
Step (4) is for the noise signal in the step (3)
Figure 563974DEST_PATH_IMAGE058
, get respectively varying strength DNoise figure.Calculate the signal to noise ratio of FHN neuron Feedback network model response under the varying strength noise.According to random resonance mechanism, in the certain noise strength range, along with the increase of noise intensity, the response signal to noise ratio is with monotone increasing; And when noise intensity was increased to certain value, if continue to increase noise intensity, the response signal to noise ratio can descend on the contrary, until noise floods signal fully.Therefore when the response signal to noise ratio reached maximum, the optimum that the cortex EEG signals will obtain on the signal to noise ratio meaning this moment strengthened.Wherein signal to noise ratio is defined as:
(3)
Wherein,
Figure DEST_PATH_IMAGE072
,
Figure DEST_PATH_IMAGE074
Represent respectively the cortex EEG signals output signal of corresponding certain window duration in power spectral density and the noise signal of interpolation, the unit of signal to noise ratio is decibel (dB).Consider the randomness of noise, therefore in power spectral density is calculated, adopt power spectral density cumulative mean method.Namely under the cortex EEG signals and noise intensity effect of identical input, repeat to ask for the response value of model, to every group of response difference rated output spectrum density, again the power spectral density of all response values is carried out cumulative mean.
The reference value that step (5) is asked for above-mentioned output signal and step (1) is sued for peace, and former amplitude scope is returned in inverse mapping, namely obtains the cortex EEG signals of the enhancing of amplitude reduction.

Claims (1)

1. the feedback system accidental resonance Enhancement Method of a cortex EEG signals is characterized in that the method comprises the steps:
Step 1, by the cortex eeg collection system, obtain the cortex EEG signals, and be divided into several nonoverlapping specific duration windows, cortex EEG signals in each window is asked for reference value, the method of asking for is with the maximum of the cortex EEG signals in each window and the intermediate value between minima, as the described reference value of this cortex EEG Processing;
Step 2, to the cortex EEG signals in each window, carry out respectively the bipolar processes of amplitude, specifically: the cortex EEG signals value in each window is deducted the reference value that step 1 is tried to achieve, acquisition has ambipolar cortex EEG signals, makes it satisfy the neuron models input signal and has ambipolar requirement;
The noise signal of step 3, the cortex EEG signals that step 2 is obtained, interpolation and the value of feedback of output signal are as the input signal based on FHN neuron Feedback network model, utilize the accidental resonance effect of FHN neuron Feedback network model, obtain the cortex brain electrical output signal that strengthens;
Step 4, in the step 3, adding under the varying strength of noise signal, calculate respectively the signal to noise ratio of FHN neuron Feedback network model response; Utilize random resonance mechanism, when signal to noise ratio reached maximum, the optimum that the cortex EEG signals will obtain under the signal to noise ratio evaluation index this moment strengthened, with this cortex EEG signals as output signal;
Step 5, the reference value that the output signal in the step 4 and step 1 are asked for are sued for peace, and former amplitude scope is returned in inverse mapping, thereby obtains the cortex EEG signals of the enhancing after the amplitude reduction.
CN 201110331467 2011-10-27 2011-10-27 Feedback system random resonance intensification method of electro-corticogram signal Active CN102499675B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110331467 CN102499675B (en) 2011-10-27 2011-10-27 Feedback system random resonance intensification method of electro-corticogram signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110331467 CN102499675B (en) 2011-10-27 2011-10-27 Feedback system random resonance intensification method of electro-corticogram signal

Publications (2)

Publication Number Publication Date
CN102499675A CN102499675A (en) 2012-06-20
CN102499675B true CN102499675B (en) 2013-09-18

Family

ID=46211857

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110331467 Active CN102499675B (en) 2011-10-27 2011-10-27 Feedback system random resonance intensification method of electro-corticogram signal

Country Status (1)

Country Link
CN (1) CN102499675B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104951082B (en) * 2015-07-09 2018-01-12 浙江大学 A kind of brain-machine interface method for strengthening EEG signals using accidental resonance
CN110840411B (en) * 2019-12-06 2022-03-11 深圳市德力凯医疗设备股份有限公司 Measuring device, storage medium and electronic equipment of anesthesia degree of depth
CN113158746B (en) * 2021-02-02 2024-04-30 杭州电子科技大学 Weak signal sensing method based on random resonance of small-world network of nerve cells
CN113180706B (en) * 2021-04-19 2023-08-15 西安交通大学 FHN stochastic resonance-based SSVEP characteristic frequency extraction method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6032072A (en) * 1998-01-30 2000-02-29 Aspect Medical Systems, Inc. Method for enhancing and separating biopotential signals
CN101076281A (en) * 2004-06-10 2007-11-21 荷兰联合利华有限公司 Apparatus and method for reducing interference

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003235823A (en) * 2002-02-15 2003-08-26 Naohiro Toda External ac noise eliminating system for biological electric signal
US20040092801A1 (en) * 2002-11-13 2004-05-13 Budimir Drakulic System for, and method of, acquiring physiological signals of a patient

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6032072A (en) * 1998-01-30 2000-02-29 Aspect Medical Systems, Inc. Method for enhancing and separating biopotential signals
CN101076281A (en) * 2004-06-10 2007-11-21 荷兰联合利华有限公司 Apparatus and method for reducing interference

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JP特开2003-235823A 2003.08.26

Also Published As

Publication number Publication date
CN102499675A (en) 2012-06-20

Similar Documents

Publication Publication Date Title
CN102499675B (en) Feedback system random resonance intensification method of electro-corticogram signal
CN108256697B (en) Prediction method for short-term load of power system
CN104809434B (en) A kind of sleep stage method removed based on single channel EEG signals eye electricity artefact
CN103576060B (en) Based on the partial discharge signal denoising method of wavelet adaptive threshold
CN107844755A (en) A kind of combination DAE and CNN EEG feature extraction and sorting technique
CN106963369B (en) Electroencephalogram relaxation degree identification method and device based on neural network model
EP2877958B1 (en) Neuron circuit and method
CN102542133B (en) Short-time wind speed forecasting method and system for wind power plant
CN110334580A (en) The equipment fault classification method of changeable weight combination based on integrated increment
CN105719029B (en) A kind of wind power combination forecasting method based on fluctuations in wind speed feature extraction
CN109038675B (en) Modeling method based on wind power fluctuation multi-scale decomposition
CN103808509A (en) Fan gear box fault diagnosis method based on artificial intelligence algorithm
CN103168306A (en) Method for determining a feature of the circadian rhythm of a subject
WO2022165737A1 (en) Method for diagnosing early weak fault signal features of marine machinery
CN109904878B (en) Multi-wind-field power generation time sequence simulation scene construction method
Wu et al. Image denoising with rectified linear units
CN105046322A (en) Method for diagnosing leading screw faults
CN107280659A (en) The processing method and system of a kind of electrocardiosignal
CN116340824A (en) Electromyographic signal action recognition method based on convolutional neural network
CN102693529B (en) Grayscale image enhancement method based on stochastic resonance mechanism of delayed self-feedback FHN (fitzhugh-nagumo) model
CN113469332A (en) Virtual synchronous generator inertia damping self-adaptive control method based on fuzzy nerves
CN111541484A (en) Optical fiber communication system optical signal-to-noise ratio monitoring method based on time delay sampling
CN108113667A (en) A kind of cardioelectric monitor method and monitoring system
CN102327116B (en) Network system random resonance restoration method of cortex electroencephalographic signal
CN102087311B (en) Method for improving measurement accuracy of power mutual inductor

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20181126

Address after: 313000 industrial zone of Donglin Town, Wuxing District, Huzhou, Zhejiang

Patentee after: Zhejiang far new material Co., Ltd.

Address before: 310018 2 street, Xiasha Higher Education Park, Hangzhou, Zhejiang

Patentee before: Hangzhou Electronic Science and Technology Univ

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20191211

Address after: No. 988, yujiahui East Road, Yaozhuang Town, Jiashan County, Jiaxing City, Zhejiang Province

Patentee after: YUNFENG PIPE INDUSTRY CO., LTD.

Address before: 313000 Zhejiang Province, Huzhou city Wuxing District East Town Industrial Zone

Patentee before: Zhejiang far new material Co., Ltd.

CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: No. 988, yujiahui East Road, Yaozhuang Town, Jiashan County, Jiaxing City, Zhejiang Province

Patentee after: Yunfeng Pipe Industry Co., Ltd

Address before: No. 988, yujiahui East Road, Yaozhuang Town, Jiashan County, Jiaxing City, Zhejiang Province

Patentee before: YUNFENG PIPE INDUSTRY Co.,Ltd.