Power frequency filtering method for reducing influence of TMS electromagnetic pulse on rat movement evoked potential
Technical Field
The invention relates to a TMS motion evoked potential power frequency interference filtering method. In particular to a power frequency filtering method for reducing the influence of TMS electromagnetic pulse on the rat movement evoked potential.
Background
Transcranial Magnetic Stimulation (TMS) technology is a magnetic stimulation technology that uses a time-varying pulsed magnetic field to act on the central nervous system (mainly the brain) to change the membrane potential of cortical nerve cells, so that induced current is generated to affect the metabolism in the brain and the neuroelectrical activity, thereby triggering a series of physiological and biochemical reactions. Motor Evoked Potentials (MEPs) are motor compound potentials recorded at the target muscle when a single TMS pulse is applied to the motor cortex. The existing rat myoelectricity acquisition mode comprises a needle electrode and a surface electrode, wherein the needle electrode is high in acquisition precision and good in muscle positioning, but needs to be implemented and recorded in an anesthesia state, and the surface electrode is very suitable for recording the movement evoked potential in a waking state. However, it is still difficult to accurately extract the motor-induced potential, latency and amplitude recorded by the surface electrode in the awake state. Compared with the human body, the rat limb surface has fine villi which are difficult to completely eliminate, the impedance is large, and the resistance reducing and resisting capacity of the scrub cream and the conductive cream is limited, so that the power frequency interference of the motion-induced potential recorded by the surface electrode is great.
At present, the existing power frequency interference filtering method is mainly to perform filtering by a way of constructing a digital wave trap, and a surface electromyographic signal recorded at a target muscle also contains a TMS electromagnetic pulse with a very large amplitude besides a motion evoked potential, and the interval of the TMS electromagnetic pulse is only a few milliseconds before the motion evoked potential.
Disclosure of Invention
The invention aims to solve the technical problem of providing a power frequency filtering method which can accurately extract the movement evoked potential waveform, amplitude and latency of a rat under the action of TMS and reduce the influence of TMS electromagnetic pulses on the movement evoked potential of the rat.
The technical scheme adopted by the invention is as follows: a power frequency filtering method for reducing the influence of TMS electromagnetic pulses on rat movement evoked potentials comprises the following steps:
1) acting a single TMS pulse on a motor cortex of a rat, and recording a hindlimb surface electromyographic signal of the rat in an awake state by adopting a sampling frequency of more than 1000Hz, wherein the signal is represented by x (N), N is 1,2, …, and N is the data point number of the electromyographic signal;
2) carrying out Fourier transform on the electromyographic signal x (N), setting Fourier coefficients corresponding to power frequency interference and direct current components to zero, keeping the other Fourier coefficients unchanged, carrying out inverse Fourier transform, realizing filtering of the power frequency interference, and obtaining a filtered electromyographic signal which is represented as y (N), wherein N is 1,2, …, and N is the data point number of the electromyographic signal;
3) extracting motion-evoked potential from the filtered electromyographic signal y (n) to obtain the motion-evoked potential P with power frequency interference being filteredmM is 1,2, …, M is the number of movement-induced potentials in the recorded electromyographic signals.
The single TMS pulse is acted on the motor cortex of the rat in the step 1), and Magstim rapid is adopted2The 8-shaped coil stimulates a rat movement region, a TMS stimulation mode is single-pulse stimulation, the stimulation intensity is 60% of the maximum output, 30 pulses are continuously output, and the pulse interval is larger than 1 s.
The step 2) comprises the following steps:
(1) fourier transform is carried out on the electromyographic signals x (n) to obtain Fourier coefficients X (k) corresponding to W frequency points,
wherein, W is the frequency point number, and W is N;
(2) calculating the real frequency values f (k) corresponding to the W frequency points,
f(k)=(k-1)·Fs/W,
wherein Fs is sampling frequency which is more than 1000 Hz;
(3) setting Fourier coefficients corresponding to 50Hz power frequency interference and 50Hz frequency multiplication and direct current components to zero, keeping the rest Fourier coefficients unchanged to obtain new Fourier coefficients Y (k),
wherein f is0Is a DC component, 50Hz and 50Hz frequency multiplication, f0={0,50,100,150,200,250,300,350,400,450,500};
(4) Carrying out inverse Fourier transform on the new Fourier coefficient Y (k) to obtain an electromyographic signal which is shown as y (n) after power frequency interference is removed,
the step 3) comprises the following steps:
(1) identifying TMS electromagnetic pulse in the electromyographic signals according to the electromyographic signals y (n) with the power frequency interference removed, wherein the amplitude of the electromyographic signals y (n) is larger than 20000 mu V and is a maximum value, and the time point corresponding to the TMS electromagnetic pulse is represented as TlL is 1,2, …, L is the number of TMS electromagnetic pulses,
Tl=n/Fs if y(n)>20000&y(n)>y(n-1)&y(n)>y(n+1),n=2,3,…,N-1,
wherein Fs is sampling frequency which is more than 1000 Hz;
(2) at a time point TlFor marking, the movement evoked potentials of the filtered myoelectric signals y (n) are identified within 150ms after the TMS electromagnetic pulse, the movement evoked potential peak-to-peak value is larger than 100 μ V, M movement evoked potentials are identified in total, and the time point corresponding to the TMS electromagnetic pulse is taken as a mark and is marked as T'm,T'm∈Tl,m=1,2,…,M,M≤L;
(3) Time point T 'is intercepted from electromyographic signal y (n) after power frequency interference is filtered out'mFirst 100ms to T'mAnd the electromyographic signal of the last 200ms is the movement evoked potential containing the TMS electromagnetic pulse, and the waveform, the amplitude and the latency of the movement evoked potential are read out.
The invention relates to a power frequency filtering method for reducing the influence of TMS electromagnetic pulse on rat movement evoked potential, wherein single TMS pulse is acted on rat motor cortex, and a higher sampling frequency is adopted to record the surface electromyographic signal of the target muscle of the rat in a waking state; carrying out Fourier transform on the electromyographic signals, setting Fourier coefficients corresponding to power frequency interference to zero, keeping the rest Fourier coefficients unchanged, and carrying out inverse Fourier transform to obtain the electromyographic signals after power frequency interference is filtered; and identifying the movement evoked potential from the electromyographic signals after the power frequency interference is filtered out to obtain the movement evoked potential after the power frequency interference is filtered out. The method for filtering power frequency interference by inverse Fourier transform does not need to establish a digital wave trap, is insensitive to sharp pulse, can greatly reduce the influence of TMS electromagnetic pulse on the movement evoked potential in the filtering process, can better reserve the self-owned characteristics of the movement evoked potential, and has important significance for accurately extracting the movement evoked potential waveform, amplitude and latency of a rat under the action of TMS.
Drawings
FIG. 1 is a flow chart of a power frequency filtering method for reducing the influence of TMS electromagnetic pulses on rat movement evoked potentials according to the present invention;
FIG. 2a is a diagram of a raw electromyography signal;
FIG. 2b is a diagram of motion-evoked potentials filtered using the method of the present invention;
fig. 2c is a diagram of motion-induced potentials filtered using a conventional digital trap.
Detailed Description
The following describes the power frequency filtering method for reducing the influence of the TMS electromagnetic pulse on the rat motor evoked potential according to the present invention in detail with reference to the following embodiments and the accompanying drawings.
The invention discloses a power frequency filtering method for reducing the influence of TMS electromagnetic pulses on rat movement evoked potentials, which particularly filters power frequency interference by adopting inverse Fourier transform. Firstly, transiently anaesthetizing a rat by using a gas anaesthesia machine, binding and fixing the rat, shaving the periphery of a target muscle of the rat, and sticking a myoelectric measuring electrode and a grounding electrode; secondly, acting a single TMS pulse on the motor cortex of the rat, and recording the surface electromyographic signals of the target muscles of the rat in the waking state by adopting higher sampling frequency; then carrying out Fourier transform on the electromyographic signals, setting Fourier coefficients corresponding to power frequency interference to zero, keeping the rest Fourier coefficients unchanged, and carrying out inverse Fourier transform to obtain the electromyographic signals after power frequency interference is filtered; and finally, identifying the movement evoked potential from the electromyographic signals after the power frequency interference is filtered out, and obtaining the movement evoked potential after the power frequency interference is filtered out.
The invention discloses a power frequency filtering method for reducing influence of TMS electromagnetic pulse on rat movement evoked potential. Then the following steps are carried out:
1) acting a single TMS pulse on a motor cortex of a rat, and recording a hindlimb surface electromyographic signal of the rat in an awake state by adopting a sampling frequency of more than 1000Hz, wherein the signal is represented by x (N), N is 1,2, …, and N is the data point number of the electromyographic signal;
the single TMS pulse acts on the motor cortex of the rat by adopting Magstimrapid2The 8-shaped coil stimulates a rat movement region, a TMS stimulation mode is single-pulse stimulation, the stimulation intensity is 60% of the maximum output, 30 pulses are continuously output, and the pulse interval is larger than 1 s.
Although the rats were subjected to transient gas anesthesia at the initial stage, the rats were awake after preparations such as rat restraint fixation, shaving, and electrode adhesion. The sampling frequency of more than 1000Hz ensures the frequency domain resolution of the electromyographic signals after Fourier transform, and is beneficial to distinguishing the Fourier coefficients of the power frequency band and the Fourier coefficients of other frequency bands.
2) Carrying out Fourier transform on the electromyographic signal x (N), setting Fourier coefficients corresponding to power frequency interference and direct current components to zero, keeping the other Fourier coefficients unchanged, carrying out inverse Fourier transform, realizing filtering of the power frequency interference, and obtaining a filtered electromyographic signal which is represented as y (N), wherein N is 1,2, …, and N is the data point number of the electromyographic signal; the method comprises the following steps:
(1) fourier transform is carried out on the electromyographic signals x (n) to obtain Fourier coefficients X (k) corresponding to W frequency points,
wherein, W is the frequency point number, and W is N;
(2) calculating the real frequency values f (k) corresponding to the W frequency points,
f(k)=(k-1)·Fs/N,
wherein Fs is sampling frequency which is more than 1000 Hz;
(3) setting Fourier coefficients corresponding to 50Hz power frequency interference and 50Hz frequency multiplication and direct current components to zero, and removing power frequency interference
Removing low-frequency drift while disturbing, keeping the rest Fourier coefficients unchanged to obtain new Fourier coefficients Y (k),
wherein f is0Is a DC component, 50Hz and 50Hz frequency multiplication, f0={0,50,100,150,200,250,300,350,400,450,500};
(4) Carrying out inverse Fourier transform on the new Fourier coefficient Y (k) to obtain an electromyographic signal which is shown as y (n) after power frequency interference is removed,
3) extracting motion-evoked potential from the filtered electromyographic signal y (n) to obtain the motion-evoked potential P with power frequency interference being filteredmM is 1,2, …, M is the number of movement-induced potentials in the recorded electromyographic signals. The method comprises the following steps:
(1) identifying TMS electromagnetic pulse in the electromyographic signals according to the electromyographic signals y (n) with the power frequency interference removed, wherein the amplitude of the electromyographic signals y (n) is larger than 20000 mu V and is a maximum value, and the time point corresponding to the TMS electromagnetic pulse is represented as TlL is 1,2, …, L is the number of TMS electromagnetic pulses,
Tl=n/Fs if y(n)>20000&y(n)>y(n-1)&y(n)>y(n+1),n=2,3,…,N-1,
fs is sampling frequency which is more than 1000 Hz;
(2) at a time point TlFor marking, the movement evoked potentials of the filtered myoelectric signals y (n) are identified within 150ms after the TMS electromagnetic pulse, the movement evoked potential peak-to-peak value is larger than 100 μ V, M movement evoked potentials are identified in total, and the time point corresponding to the TMS electromagnetic pulse is taken as a mark and is marked as T'm,T'm∈Tl,m=1,2,…,M,M≤L;
(3) Time point T 'is intercepted from electromyographic signal y (n) after power frequency interference is filtered out'mFirst 100ms to T'mAnd the electromyographic signal of the last 200ms is the movement evoked potential containing the TMS electromagnetic pulse, and the waveform, the amplitude and the latency of the movement evoked potential are read out.
The rat hind limb muscle group was used as the target muscle, monopulse TMS was applied to the motor cortex, and the recorded raw myoelectric signal is shown in fig. 2 a. The inverse Fourier transform provided by the invention is used for filtering the power frequency interference in the movement evoked potential under the rat TMS action, and the obtained movement evoked potential is shown in figure 2 b. The movement evoked potential power frequency interference under rat TMS action is filtered by using a traditional digital wave trap, and the obtained movement evoked potential is shown in figure 2 c. Compared with the power frequency interference filtering result, the method provided by the invention can better reserve the amplitude and the latency of the motion evoked potential.