CN116510183B - Transcranial magnetic stimulation device, electronic device, and storage medium - Google Patents

Transcranial magnetic stimulation device, electronic device, and storage medium Download PDF

Info

Publication number
CN116510183B
CN116510183B CN202310775449.4A CN202310775449A CN116510183B CN 116510183 B CN116510183 B CN 116510183B CN 202310775449 A CN202310775449 A CN 202310775449A CN 116510183 B CN116510183 B CN 116510183B
Authority
CN
China
Prior art keywords
pulse
transcranial magnetic
test
target
channels
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
CN202310775449.4A
Other languages
Chinese (zh)
Other versions
CN116510183A (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.)
First Affiliated Hospital of Nanchang University
Original Assignee
First Affiliated Hospital of Nanchang 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 First Affiliated Hospital of Nanchang University filed Critical First Affiliated Hospital of Nanchang University
Priority to CN202310775449.4A priority Critical patent/CN116510183B/en
Publication of CN116510183A publication Critical patent/CN116510183A/en
Application granted granted Critical
Publication of CN116510183B publication Critical patent/CN116510183B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N2/00—Magnetotherapy
    • A61N2/02—Magnetotherapy using magnetic fields produced by coils, including single turn loops or electromagnets
    • 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]
    • A61B5/372—Analysis of electroencephalograms
    • 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]
    • A61B5/377—Electroencephalography [EEG] using evoked responses
    • A61B5/383—Somatosensory stimuli, e.g. electric stimulation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/48—Other medical applications
    • A61B5/4836—Diagnosis combined with treatment in closed-loop systems or methods
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N2/00—Magnetotherapy
    • A61N2/004—Magnetotherapy specially adapted for a specific therapy
    • A61N2/006—Magnetotherapy specially adapted for a specific therapy for magnetic stimulation of nerve tissue
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Psychology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Psychiatry (AREA)
  • Neurology (AREA)
  • Magnetic Treatment Devices (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

The application provides transcranial magnetic stimulation equipment, electronic equipment and a storage medium. Determining the stimulus intensity Si of the ith transcranial magnetic test; issuing a plurality of first pulses to the subject; collecting a plurality of first electroencephalogram signals; performing effective response component analysis on the plurality of first electroencephalogram signals to determine whether effective response components exist in the ith transcranial magnetic test; if the excitation intensity Si+1 exists, determining the excitation intensity Si+1 of the i+1th transcranial magnetic test until the excitation intensity Si+1 does not exist; if not, determining the stimulation intensity Dj of the jth transcranial magnetic test; issuing a plurality of second pulses to the subject; collecting a plurality of second electroencephalogram signals; performing effective response component analysis on the plurality of second electroencephalogram signals to determine whether effective response components exist in the jth transcranial magnetic test; if not, determining the stimulation intensity Dj+1 of the j+1th transcranial magnetic test until effective response components exist; if present, the stimulus intensity Dj is taken as the target stimulus intensity.

Description

Transcranial magnetic stimulation device, electronic device, and storage medium
Technical Field
The application relates to the technical field of medical treatment, in particular to transcranial magnetic stimulation equipment, electronic equipment and a storage medium.
Background
Transcranial magnetic stimulation techniques produce magnetic fields of equal strength by nuclear magnetic resonance through coils placed on the scalp and act on brain tissue unattenuated penetrating the skull, thereby depolarizing neurons, inducing the conduction of nerve impulses. Functional reorganization of the cerebral cortex is achieved by effects on neuronal excitability, networking and plasticity. The transcranial magnetic stimulation technology has the advantages of painless, safety, low cost, strong adaptability, high plasticity and the like, and has an important application value in the nerve function rehabilitation treatment of mental and neurological diseases.
In the application of transcranial magnetic stimulation to the treatment of brain diseases, the stimulation intensity needs to be set first. Currently, the stimulus intensity is set by measuring the motor threshold of the individual. Specifically, transcranial magnetic stimulation is applied to the primary motor region of the brain, induced myoelectric signals are synchronously acquired at the motor nerve ends of the contralateral upper limbs, the minimum stimulation intensity capable of stably inducing myoelectric responses and having a peak value greater than 50 μV is referred to as the motor threshold of the person, and the stimulation intensity is set to 90% -120% of the motor threshold.
Due to the difference of the cerebral cortex structure and the scalp-to-cortex thickness distribution, when the stimulation intensity set by utilizing the motion threshold value is applied to transcranial magnetic stimulation treatment of other cortex layers, the irritancy intensity has unreasonable problems such as larger or smaller, namely, the set stimulation intensity has lower precision, and the transcranial magnetic stimulation treatment effect is unstable.
Disclosure of Invention
The application provides transcranial magnetic stimulation equipment, electronic equipment and a storage medium, which are used for obtaining the stimulation intensity which is matched with a subject with high precision through multiple transcranial magnetic tests before transcranial magnetic stimulation treatment, so that the transcranial magnetic stimulation treatment effect is improved.
In a first aspect, the present application provides a transcranial magnetic stimulation device, including an electroencephalogram acquisition module, an evoked potential online analysis module, a transcranial magnetic stimulation module, and a stimulation intensity control module;
the stimulation intensity control module is used for determining the stimulation intensity Si of the ith transcranial magnetic test based on the stimulation intensity Si-1 of the ith-1 th transcranial magnetic test, wherein the stimulation intensity Si is smaller than the stimulation intensity Si-1;
the transcranial magnetic stimulation module is used for sending a plurality of first pulses corresponding to the stimulation intensity Si to a subject;
the electroencephalogram acquisition module is used for acquiring a plurality of first electroencephalograms corresponding to each first pulse on a plurality of channels, wherein the plurality of first electroencephalograms correspond to the plurality of channels one by one;
the evoked potential online analysis module is used for carrying out effective response component analysis on a plurality of first electroencephalogram signals corresponding to each first pulse and determining whether effective response components exist in the ith transcranial magnetic test;
The stimulation intensity control module is further used for determining the stimulation intensity Si+1 of the ith+1th transcranial magnetic test based on the stimulation intensity Si if the effective response component exists, so as to perform the ith+1th transcranial magnetic test until the effective response component does not exist;
the stimulation intensity control module is further used for determining the stimulation intensity Dj of the jth transcranial magnetic test based on the stimulation intensity Dj-1 of the jth transcranial magnetic test if no effective response component exists, wherein the stimulation intensity Dj is larger than the stimulation intensity Dj-1, and when j=1, the stimulation intensity D1 is the stimulation intensity Si;
the transcranial magnetic stimulation module is further configured to deliver a plurality of second pulses to the subject corresponding to the stimulation intensity Dj;
the electroencephalogram acquisition module is further used for acquiring a plurality of second electroencephalograms corresponding to each second pulse on the plurality of channels, wherein the plurality of second electroencephalograms correspond to the plurality of channels one by one;
the evoked potential online analysis module is further used for carrying out effective response component analysis on a plurality of second brain electrical signals corresponding to each second pulse and determining whether effective response components exist in the jth transcranial magnetic test;
The stimulation intensity control module is further used for determining the stimulation intensity Dj+1 of the j+1th transcranial magnetic test based on the stimulation intensity Dj if no effective response component exists, so as to perform the j+1th transcranial magnetic test until the effective response component exists; if an effective response component exists, the stimulus intensity Dj is taken as a target stimulus intensity.
In a second aspect, the present application provides a method for setting stimulation intensity, the method is applied to a transcranial magnetic stimulation device, the transcranial magnetic stimulation device includes an electroencephalogram acquisition module, an evoked potential on-line analysis module, a transcranial magnetic stimulation module, and a stimulation intensity control module; the method comprises the following steps:
the stimulation intensity control module determines the stimulation intensity Si of the ith transcranial magnetic test based on the stimulation intensity Si-1 of the ith-1 th transcranial magnetic test, wherein the stimulation intensity Si is smaller than the stimulation intensity Si-1;
the transcranial magnetic stimulation module delivers a plurality of first pulses corresponding to the stimulation intensity Si to a subject;
the electroencephalogram acquisition module acquires a plurality of first electroencephalograms corresponding to each first pulse on a plurality of channels, wherein the plurality of first electroencephalograms correspond to the plurality of channels one by one;
The evoked potential online analysis module analyzes the effective response components of a plurality of first electroencephalograms corresponding to each first pulse, and determines whether effective response components exist in the ith transcranial magnetic test;
if the effective response component exists, the stimulation intensity control module determines the stimulation intensity Si+1 of the ith+1st transcranial magnetic test based on the stimulation intensity Si so as to perform the ith+1st transcranial magnetic test until the effective response component does not exist;
if no effective response component exists, the stimulation intensity control module determines the stimulation intensity Dj of the jth transcranial magnetic test based on the stimulation intensity Dj-1 of the jth transcranial magnetic test, wherein the stimulation intensity Dj is greater than the stimulation intensity Dj-1, and when j=1, the stimulation intensity D1 is the stimulation intensity Si;
the transcranial magnetic stimulation module delivers a plurality of second pulses to the subject corresponding to the stimulation intensity Dj;
the electroencephalogram acquisition module acquires a plurality of second electroencephalograms corresponding to each second pulse on the plurality of channels, wherein the plurality of second electroencephalograms correspond to the plurality of channels one by one;
the evoked potential online analysis module analyzes the effective response components of a plurality of second brain electrical signals corresponding to each second pulse, and determines whether effective response components exist in the jth transcranial magnetic test;
If no effective response component exists, the stimulation intensity control module determines the stimulation intensity Dj+1 of the j+1th transcranial magnetic test based on the stimulation intensity Dj so as to perform the j+1th transcranial magnetic test until the effective response component exists; if the effective response component exists, the stimulus intensity control module takes the stimulus intensity Dj as the target stimulus intensity.
In a third aspect, the present application provides an electronic device comprising: a processor and a memory, the processor being connected to the memory, the memory being for storing a computer program, the processor being for executing the computer program stored in the memory to cause the electronic device to perform the method as described in the second aspect.
In a fourth aspect, the present application provides a computer-readable storage medium storing a computer program that causes a computer to perform the method according to the second aspect.
In a fifth aspect, the present application provides a computer program product comprising a non-transitory computer readable storage medium storing a computer program, the computer being operable to cause a computer to perform the method of the second aspect.
The embodiment of the application has the following beneficial effects:
it can be seen that in the embodiment of the present application, before the transcranial magnetic stimulation treatment is performed on the subject, the subject (i.e., the patient) is subjected to multiple transcranial magnetic tests in the area to be treated of the subject, and when the transcranial magnetic tests are performed multiple times, the stimulus intensity is gradually reduced (i.e., the stimulus intensity Si is smaller than the stimulus intensity Si-1), and the analysis of the effective response component is performed based on the first electroencephalogram signal of each transcranial magnetic test until the stimulus intensity without the effective response component is found, so that the target stimulus intensity suitable for the subject should be greater than the stimulus intensity because the patient needs to respond to achieve the treatment during the transcranial magnetic stimulation treatment. Thus, in the present application, the stimulus intensity is gradually increased from the stimulus intensity (i.e., the stimulus intensity Dj is greater than the stimulus intensity Dj-1), and the analysis of the effective response component is performed based on the second electroencephalogram signal of each transcranial magnetic test until the effective response component is again generated, that is, the stimulus intensity just enables the patient to respond, and the stimulus intensity is taken as the target stimulus intensity. Therefore, through the technical scheme of the application, the set stimulation intensity can ensure that the patient can be stimulated to respond, and the stimulation intensity can not be excessively high, so that the stimulation intensity matched with the area to be treated of the patient can be set with high precision, and the transcranial magnetic stimulation treatment can be carried out on the patient by using the target stimulation intensity, and the treatment effect can be improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for the description of the embodiments will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present application, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is a schematic view of a scenario for setting stimulus intensity according to an embodiment of the present application;
FIG. 2 is a schematic diagram of a transcranial magnetic stimulation device according to an embodiment of the present application;
fig. 3 is a schematic diagram of multi-pulse and multi-channel electroencephalogram signal acquisition according to an embodiment of the present application;
FIG. 4 is a schematic diagram of determining a first covariance matrix and a second covariance matrix according to an embodiment of the application;
FIG. 5 is a schematic diagram showing a gradual decrease of stimulus intensity until no effective response component exists according to an embodiment of the present application;
FIG. 6 is a flow chart of a method for setting stimulus intensity according to an embodiment of the present application;
fig. 7 is a schematic diagram of an electronic device according to an embodiment of the present application.
Detailed Description
The following description of the embodiments of the present application will be made clearly and fully with reference to the accompanying drawings, in which it is evident that the embodiments described are some, but not all embodiments of the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
The terms "first," "second," "third," and "fourth" and the like in the description and in the claims and drawings are used for distinguishing between different objects and not necessarily for describing a particular sequential or chronological order. Furthermore, the terms "comprise" and "have," as well as any variations thereof, are intended to cover a non-exclusive inclusion. For example, a process, method, system, article, or apparatus that comprises a list of steps or elements is not limited to only those listed steps or elements but may include other steps or elements not listed or inherent to such process, method, article, or apparatus.
Reference herein to "an embodiment" means that a particular feature, result, or characteristic described in connection with the embodiment may be included in at least one embodiment of the application. The appearances of such phrases in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments. Those of skill in the art will explicitly and implicitly appreciate that the embodiments described herein may be combined with other embodiments.
In order to facilitate understanding of the technical scheme of the present application, explanation and explanation of the present application related to the related art will be first provided.
Currently, in the application of transcranial magnetic stimulation to the treatment of brain diseases, the stimulation intensity needs to be set first. Because the response of the cerebral cortex of a patient to transcranial magnetic stimulation has extremely strong variability among individuals, transcranial magnetic stimulation treatment cannot be performed directly by manually setting a fixed stimulation intensity. In this regard, the conventional approach is to set the stimulus intensity of the individual by measuring the movement threshold of the individual. Specifically, transcranial magnetic stimulation is applied to a primary motor region of the brain of an individual, induced myoelectric signals are synchronously acquired at the motor nerve ends of the contralateral upper limbs, and the minimum stimulation intensity capable of stably inducing myoelectric responses and having a peak value of more than 50 μV is called the motor threshold of the individual. Then 90% -120% of the motion threshold value is used as the stimulation intensity of the human during transcranial magnetic stimulation treatment when the human is subjected to transcranial magnetic stimulation treatment.
However, the scheme of setting the stimulus intensity by measuring the movement threshold has the following drawbacks:
first, the stimulus intensity set by measuring the motion threshold is theoretically only applicable to transcranial magnetic stimulation treatment associated with the motion zone. The differences in brain cortical structure and scalp-to-cortex thickness distribution result in transcranial magnetic stimulation treatments that cannot reasonably generalize to other cortex with the stimulation intensity set by the motion zone detection. Therefore, when the stimulation intensity set by measuring the motion threshold is applied to transcranial magnetic stimulation treatment of other cortex, the irritancy intensity has unreasonable problems such as larger or smaller stimulation intensity, namely the accuracy of the set stimulation intensity is lower, so that the transcranial magnetic stimulation treatment effect is unstable.
Setting the stimulation intensity by measuring the motor threshold is severely dependent on the motor nerve pathway integrity of the patient, so that the motor threshold of the patient cannot be accurately measured once a lesion or damage occurs on the motor nerve pathway, and would make it impossible to formulate an individualized stimulation intensity for such a patient. In addition, especially for patients with significant brain damage, such patients may fail to measure the motor threshold after the brain motor area is damaged, and it is not possible to tailor the individual stimulation intensity for such patients.
In addition, it is not possible to induce a neurological response in the patient's severely diseased cortex due to transcranial magnetic stimulation. At present, by measuring the stimulus intensity set by the exercise threshold, whether the stimulus intensity can sufficiently regulate and control the change of the nerve activity cannot be known on the lesion cortex, and whether the set stimulus intensity can sufficiently regulate and control the change of the nerve activity cannot be ensured, so that the setting accuracy of the stimulus intensity is lower.
In conclusion, in transcranial magnetic stimulation treatment, the existing method for setting the stimulation intensity of transcranial magnetic stimulation has lower precision, so that the transcranial magnetic stimulation treatment effect is poor.
Firstly, the application mainly comprises two transcranial magnetic testing iterative processes for setting the stimulation intensity, wherein the first transcranial magnetic testing iterative process mainly roughly determines the stimulation intensity, and the second transcranial magnetic testing iterative process only needs to finely adjust the stimulation intensity obtained in the first transcranial magnetic testing iterative process so as to finely determine the target stimulation intensity, and the target stimulation intensity is used for transcranial magnetic stimulation treatment of the subject. For convenience of distinction, any one transcranial magnetic test in the iterative process of the first transcranial magnetic test is called an ith transcranial magnetic test, a pulse emitted in the ith transcranial magnetic test process is called a first pulse, and an acquired electroencephalogram signal is called a first electroencephalogram signal; any one transcranial magnetic test in the second transcranial magnetic test iteration process is called a jth transcranial magnetic test, pulses emitted in the jth transcranial magnetic test process are called second pulses, and acquired brain electrical signals are called second brain electrical signals.
The first transcranial magnetic testing iteration process is a process of starting with preset stimulation intensity, gradually reducing the stimulation intensity until no effective response component is detected, ending the first transcranial magnetic testing iteration process, then starting with the stimulation intensity of the first transcranial magnetic testing iteration process, gradually increasing the stimulation intensity, performing the second transcranial magnetic testing iteration process until the effective response component is detected again, ending the second transcranial magnetic testing iteration process, and taking the stimulation intensity obtained when the second transcranial magnetic testing iteration process is ended as target stimulation intensity.
Furthermore, the stimulation intensities referred to herein may also be referred to as transcranial magnetic stimulation intensities, both being consistent in nature.
Referring to fig. 1, fig. 1 is a schematic view of a scenario for setting stimulus intensity according to an embodiment of the present application.
As shown in fig. 1, a coil is placed on the brain of a subject before the subject is subjected to transcranial magnetic stimulation treatment, and then transcranial magnetic stimulation equipment performs a plurality of transcranial magnetic tests on the subject to obtain a target stimulation intensity, and then the target stimulation intensity is used as the stimulation intensity of the subject to be subjected to transcranial magnetic stimulation treatment. Wherein, the ith transcranial stimulation comprises the following steps:
The transcranial magnetic stimulation equipment firstly determines the stimulation intensity Si of the ith transcranial magnetic test based on the stimulation intensity Si-1 of the ith-1 th transcranial magnetic test, wherein the stimulation intensity Si is smaller than the stimulation intensity Si-1; then, a plurality of first pulses corresponding to the stimulus intensity Si are given to the subject, that is, the plurality of first pulses are given to the coil, and a magnetic field is generated by the coil to act on brain tissue of the subject, so that the brain nerve impulse conduction of the subject is induced; then, a plurality of first electroencephalogram signals corresponding to each first pulse are acquired on a plurality of channels, wherein the plurality of first electroencephalogram signals are in one-to-one correspondence with the plurality of channels. Finally, carrying out effective response component analysis on a plurality of first electroencephalogram signals corresponding to each first pulse, and determining whether effective response components exist in the ith transcranial magnetic test; if the effective response component exists, determining the stimulation intensity Si+1 of the ith+1st transcranial magnetic test based on the stimulation intensity Si, so as to perform the ith+1st transcranial magnetic test until the effective response component does not exist; if no effective response component exists, determining the stimulation intensity Dj of the jth transcranial magnetic test based on the stimulation intensity Dj-1 of the jth transcranial magnetic test, wherein the stimulation intensity Dj is larger than the stimulation intensity Dj-1, and when j=1, the stimulation intensity D1 is the stimulation intensity Si; issuing a plurality of second pulses to the subject corresponding to the stimulus intensity Dj; collecting a plurality of second electroencephalogram signals corresponding to each second pulse on the channels, wherein the second electroencephalogram signals are in one-to-one correspondence with the channels; performing effective response component analysis on a plurality of second electroencephalogram signals corresponding to each second pulse, and determining whether effective response components exist in the jth transcranial magnetic test; if the effective response component does not exist, determining the stimulation intensity Dj+1 of the j+1th transcranial magnetic test based on the stimulation intensity Dj so as to perform the j+1th transcranial magnetic test until the effective response component exists; if an effective response component is present, the stimulus intensity Dj is taken as the target stimulus intensity.
Referring to fig. 2, fig. 2 is a schematic diagram of a transcranial magnetic stimulation device according to an embodiment of the present application.
As shown in fig. 2, the transcranial magnetic stimulation device comprises an electroencephalogram acquisition module, an evoked potential on-line analysis module, a transcranial magnetic stimulation module and a stimulation intensity control module.
Illustratively, in performing the ith transcranial magnetic test, the stimulation intensity control module determines the stimulation intensity Si of the ith transcranial magnetic test based on the stimulation intensity Si-1 of the ith-1 th transcranial magnetic test. Specifically, the stimulation intensity control module obtains the proportionality coefficient K of the ith-1 th transcranial magnetic test i-1 Wherein K is i-1 Less than 1; according to a preset proportionality coefficient delta K and a proportionality coefficient K i-1 Determining the scaling factor K of the ith transcranial magnetic test i Wherein K is i =K i-1 δk, δk being an integer greater than 0. Finally, the scaling factor K based on the ith transcranial magnetic testing i And determining the stimulation intensity Si of the ith transcranial magnetic test according to the maximum stimulation intensity Smax which can be emitted by the transcranial magnetic stimulation device, wherein si=k i ×Smax。
It should be noted that when i=1, the scaling factor K of the first transcranial magnetic test 1 Is a preset proportionality coefficient.
After the stimulation intensity control module sets the stimulation intensity Si at the ith transcranial magnetic test, the transcranial magnetic stimulation module may obtain the stimulation intensity Si at the ith transcranial magnetic test, e.g., the stimulation intensity control module sets the stimulation intensity Si at the ith transcranial magnetic test, and may send the stimulation intensity Si to the transcranial magnetic stimulation module. Accordingly, the transcranial magnetic stimulation module delivers a plurality of first pulses to the subject corresponding to the stimulation intensity Si based on the stimulation intensity Si. Illustratively, as shown in fig. 3, the transcranial magnetic stimulation module determines a level corresponding to an ith transcranial magnetic test based on the stimulation intensity Si, and issues a plurality of first pulses to the subject at a plurality of times through the level, as shown in fig. 3 as first pulse 1, first pulses 2, … …, first pulse n. Specifically, the coil is first placed in the region to be treated of the cerebral cortex of the subject in order to precisely determine the stimulation intensity of the region to be treated at the time of transcranial magnetic stimulation treatment. The transcranial magnetic stimulation module may then deliver a plurality of first pulses to the region to be treated through the coil at the ith transcranial magnetic test.
For example, as shown in fig. 3, the electroencephalogram acquisition module may perform electroencephalogram acquisition on multiple channels. Therefore, after the transcranial magnetic stimulation module distributes a plurality of first pulses to the subject, the electroencephalogram acquisition module acquires a plurality of first electroencephalograms corresponding to each first pulse on a plurality of channels, wherein the channels are in one-to-one correspondence with the plurality of first electroencephalograms. It should be understood that the electroencephalogram acquisition module can always acquire electroencephalogram signals in the whole process of setting the stimulus intensity, however, the electroencephalogram signal change condition of a subject only needs to be studied before and after the pulse is issued. Therefore, the first electroencephalogram signal is obtained by only taking the time of issuing each first pulse as a reference and intercepting the electroencephalogram signal acquired by the electroencephalogram acquisition module, as shown in fig. 3, the electroencephalogram signal in the time period before Deltat 1 of the time of issuing each first pulse and the electroencephalogram signal in the time period after Deltat 2 of the time of issuing each first pulse are intercepted, and a plurality of first electroencephalograms of each first pulse on a plurality of channels can be obtained.
Further, after the evoked potential online analysis module obtains a plurality of first electroencephalograms corresponding to each first pulse, effective response component analysis can be performed on the plurality of first electroencephalograms corresponding to each first pulse to determine whether effective response components exist in the ith transcranial magnetic test.
Specifically, the evoked potential online analysis module respectively interpolates a plurality of first electroencephalogram signals of each first pulse on a plurality of channels based on the moment of issuing each first pulse, so as to obtain a plurality of first test times of each first pulse on the plurality of channels. And respectively interpolating the electroencephalogram signals of each first pulse in the front delta t1 time period of the first electroencephalogram signals on each channel and interpolating the electroencephalogram signals of the rear delta t2 time period by utilizing an interpolation algorithm to obtain a plurality of first test times of each first pulse on a plurality of channels, wherein the channels correspond to the first test times one by one. And then, carrying out noise fitting on all first test times corresponding to the first pulses on the channels to obtain a noise curve. For example, the number of the first pulses is m, the number of the channels is n, and each first pulse corresponds to n first test times, and then m×n first test times are used for carrying out noise fitting on the m×n first test times in total, so as to obtain a noise curve, wherein the noise curve is used for representing a change rule of noise when electroencephalogram signal acquisition is carried out after the first pulses are issued, namely, a change rule of noise when electroencephalogram signal acquisition is carried out after each first pulse is issued. Alternatively, the noise curve may be obtained by performing noise fitting by a least square method and an attenuation fitting function.
And then, carrying out component analysis on all first test times corresponding to the first pulses on the channels by the evoked potential online analysis module to obtain a plurality of component curves corresponding to the components, wherein each component curve is used for representing the change rule of the component corresponding to the component curve after the first pulses are issued.
Specifically, n components are preset (i.e., n independent components are preset), and the first test time of each first pulse on each channel can be regarded as being formed by mixing the n independent components, and then all the first test times of the first pulses on the channels can be represented by formula (1):
wherein, the liquid crystal display device comprises a liquid crystal display device,for the first trial of the c-th first pulse on the d-th channel,/for the first trial of the c-th first pulse on the d-th channel>Is n components, wherein, c is takenThe value of d is an integer from 1 to n.
Illustratively, independent component analysis is performed based on the above formula (1) to obtain n componentsAnd a first weight matrix. Wherein, each component can be represented by a component curve corresponding to the component, and the weight matrix is composed of the weights of each component under each first test time, so the first weight matrix is a matrix of m×n rows and n columns, specifically:
Further, noise is removed from the plurality of first test times corresponding to each first pulse based on the noise curve and the plurality of component curves corresponding to the plurality of components, and a plurality of second test times corresponding to each first pulse is obtained. Illustratively, each component curve is fitted to the noise curve to obtain a root mean square error corresponding to each component curve. And determining the component curves with the root mean square error smaller than the threshold value according to the root mean square error corresponding to each component curve. And setting the weight of the component corresponding to the component curve under all first test times to be 0 for the component curve with the root mean square error smaller than the threshold value, namely setting the weight of the component under each first test time to be 0, and obtaining a second weight matrix. The weights of the components under each first test corresponding to each first pulse are obtained. For example, when the root mean square error corresponding to the first component S1 is smaller than the threshold value, the second weight matrix is:
and finally, reconstructing the first test time based on the weight of the components under the first test time of each first pulse on each channel, namely weighting the components according to the weights of the components under the first test time to obtain a second test time corresponding to the first test time, namely obtaining a third test time of each first pulse on each channel, namely obtaining a plurality of second test times of each first pulse on a plurality of channels. Then, an effective response component analysis is performed on a plurality of second trials corresponding to each first pulse, and whether an effective response component exists in the ith transcranial magnetic test is determined.
Specifically, after obtaining a plurality of second trials of each first pulse on a plurality of channels, for the second trials of each first pulse on each channel, a mean value of the magnitudes of the second trials over a preset period of time is obtained. For example, the preset time period is a previous Δt3 time period of the time when the first pulse is issued. Then, subtracting the average value from the amplitude of the second test at each time to obtain a third test corresponding to the second test, namely, performing baseline correction on the second test to obtain the third test corresponding to the second test, and further obtaining a plurality of third tests of each first pulse on a plurality of channels.
Further, as shown in fig. 4, the magnitudes of the third test times of each channel under the first pulses are averaged to obtain a fourth test time corresponding to each channel, that is, from the channel dimension, the magnitudes of the third test times of each channel under the first pulses are averaged to obtain the fourth test time corresponding to each channel. Then, according to a plurality of fourth test times corresponding to a plurality of channels, determining a first covariance matrix, namely, constructing a first amplitude matrix based on the fourth test times corresponding to each channel, wherein each row of elements in the first amplitude matrix can be constructed by the amplitude of the fourth test times corresponding to each channel at each time, and the column number is the number of the channels. A first covariance matrix is then determined based on the first magnitude matrix. Further, as shown in fig. 4, the second covariance matrix corresponding to each first pulse is determined according to the third test times of each first pulse under the plurality of channels, that is, the covariance matrix corresponding to each first pulse is determined from the pulse dimension. Specifically, a second amplitude matrix corresponding to each first pulse is constructed according to a plurality of third test times of each first pulse under a plurality of channels, wherein each row element in the second amplitude matrix is the amplitude structure of the first pulse under each channel under each time, and the column number is the number of channels. Then, a second covariance matrix corresponding to each first pulse is determined according to the second amplitude matrix corresponding to each first pulse.
And finally, filtering the first pulses according to the second covariance matrix corresponding to each first pulse and the first covariance matrix to obtain a plurality of target pulses. For example, a distance between the second covariance matrix corresponding to each first pulse and the first covariance matrix is obtained, and for example, the distance may be a Kullback Leibler distance, i.e. a KL distance. Then, performing Gaussian model fitting on a plurality of distances corresponding to the first pulses to obtain a Gaussian distribution function; and then, determining a distance threshold according to the Gaussian distribution function, determining a probability density function corresponding to each distance based on the Gaussian distribution function, and eliminating the first pulse corresponding to the distance of which the probability density function is smaller than the distance threshold to obtain a plurality of target pulses.
It can be seen that in this embodiment, the covariance matrix of the whole electroencephalogram signal is obtained by calculating the covariance matrix of the test times under all the pulses, and then the covariance matrix under each pulse is calculated, and the first pulse with a larger outlier is removed based on the KL distance between the covariance matrices, so that the interference of the abnormal first pulse is avoided, and the accuracy of analyzing the effective response components is improved.
Further, after the first pulse with the larger outlier is eliminated, effective response component analysis is performed on a plurality of third test times of each target pulse on a plurality of channels, and whether effective response components exist in the ith transcranial magnetic test is determined.
Specifically, for each target pulse, the time when each target pulse is issued is taken as a dividing line, and the third test time of each target pulse on each channel is subjected to signal division to obtain a baseline signal and an induced signal of each target pulse on each channel. Illustratively, with the time at which each target pulse is issued as a dividing line, an electroencephalogram signal in a period of Δt4 preceding the time in a third trial of the target pulse on each channel is taken as a baseline signal, and an electroencephalogram signal in a period of Δt5 following the time is taken as an evoked signal.
Then, a Gaussian fit is performed based on the magnitude of the baseline signal on each channel for each target pulse, resulting in a magnitude threshold. Illustratively, based on the amplitude of the baseline signal of each target pulse on each channel at each time instant, determining the mean and standard deviation of the amplitude of each target pulse on each channel, i.e. taking the mean and standard deviation of the amplitude at each time instant as the mean and standard deviation of the target pulse on the channel; the target mean and target standard deviation are determined based on the mean and standard deviation of the amplitude of each target pulse on each channel. For example, the average value and standard deviation of the amplitude of each target pulse on each channel may be averaged to obtain a candidate average value and a candidate standard deviation corresponding to each target pulse; and finally, respectively averaging a plurality of candidate average values and a plurality of candidate standard deviations of a plurality of target pulses to obtain the target average value and the target standard deviation. Finally, constructing Gaussian distribution corresponding to the target mean value and the target standard deviation; the magnitude threshold is determined based on the gaussian distribution and a preset confidence.
And finally, screening a plurality of target channels from the plurality of channels, for example, acquiring the distance between each channel and the coil, and taking the channel with the distance smaller than a preset threshold value as the target channel. For example, four channels around the coil may be taken as target channels. Then, the amplitude of the induction signal of each target pulse on the plurality of target channels at each time is averaged to obtain an average amplitude corresponding to each target pulse; averaging a plurality of average amplitudes corresponding to the target pulses to obtain target amplitudes; if the target amplitude is greater than or equal to the amplitude threshold, determining that an effective response component exists in the ith transcranial magnetic test; and if the target amplitude is smaller than the amplitude threshold, determining that no effective response component exists in the ith transcranial magnetic test.
Finally, if no effective response component exists, determining target stimulus intensity by a stimulus intensity control module based on the stimulus intensity Si; if the effective response component exists in the ith transcranial magnetic testing, the stimulation intensity control module determines the stimulation intensity Si+1 of the ith+1th transcranial magnetic testing based on the stimulation intensity Si, then the ith+1th transcranial magnetic testing is carried out through mutual coordination among the transcranial magnetic stimulation module, the electroencephalogram acquisition module and the evoked potential online analysis module, namely the ith+1th transcranial magnetic testing is continued until the effective response component does not exist, and then the target stimulation intensity can be determined based on the stimulation intensity when the effective response component does not exist. The i+1th transcranial magnetic test is similar to the i th transcranial magnetic test, and will not be described.
For example, fig. 5 illustrates three channels and three first pulses. As shown in fig. 5, when the stimulus intensity is k1×smax, it can be seen from the acquired brain electrical signals that the fluctuation of the brain electrical signals before and after the pulse is emitted is relatively large, the presence of the effective response component at the stimulus intensity is k1×smax can be detected, and then the stimulus intensity is gradually reduced, and as can be seen from fig. 5, when the stimulus intensity is gradually reduced from k1×smax to k2×smax, k3×smax, and k4×smax, the fluctuation of the brain electrical signals before and after the pulse is emitted is gradually reduced, and when the stimulus intensity is k4×smax, the effective response component is not detected, and the target stimulus intensity can be determined from the stimulus intensity k4×smax.
First, during the first transcranial magnetic testing, if no effective response component is present, it is indicated that the subject's cerebral cortex is unable to perform effective transcranial stimulation, and the test is terminated. In the first transcranial magnetic testing procedure, the stimulation intensity is gradually reduced if an active response component is present until no active response component is present, and the target stimulation intensity may be determined based on the stimulation intensity when no active response component is present. The present application is mainly described by taking the example that no response component exists in the ith transcranial magnetic test.
In practical use, the stimulus intensity is decreased gradually in a large proportion so that the stimulus intensity when no response component appears can be obtained quickly. Therefore, when the i-th transcranial magnetic test does not have an effective response component, the stimulus intensity (i.e., the critical value) when the effective response component is actually not just present is not necessarily the stimulus intensity Si, that is, the stimulus intensity Si cannot be directly used as the target stimulus intensity, and fine adjustment is required based on the stimulus intensity to determine the target stimulus intensity with high accuracy. The process of fine tuning is described below with respect to the absence of a significant response component from the ith transcranial magnetic test.
Illustratively, the stimulation intensity control module further determines a stimulation intensity Dj for the jth transcranial magnetic test based on the stimulation intensity Dj-1 for the jth transcranial magnetic test, wherein the stimulation intensity Dj is greater than the stimulation intensity Dj-1, and the stimulation intensity D1 is the stimulation intensity Si when j=1. In particular, the proportionality coefficient R of the j-i th transcranial magnetic test is obtained j-1 Wherein R is j-1 Less than 1; according to a preset proportionality coefficient delta R and a proportionality coefficient R j-1 Determining the scaling factor R of the jth transcranial magnetic test j Wherein R is j =R j-1 +δr. Finally, the proportionality coefficient R based on jth transcranial magnetic testing j And determining the stimulation intensity Dj of the jth transcranial magnetic test according to the maximum stimulation intensity Smax which can be emitted by the transcranial magnetic stimulation device, wherein dj=R ji ×Smax。
The δr is an integer greater than 0 and is smaller than δk, so that in the fine tuning process, from the stimulus intensity Si, the stimulus intensity is gradually increased, after the stimulus intensity is increased each time, whether effective response components exist or not is determined, the proportion of the stimulus intensity is smaller than the proportion of the stimulus intensity is reduced each time, and therefore the stimulus intensity without effective response components can be found out, and the target stimulus intensity can be determined with high precision.
Further, the transcranial magnetic stimulation module also delivers a plurality of second pulses to the subject corresponding to the stimulation intensity Dj; the electroencephalogram acquisition module is used for acquiring a plurality of second electroencephalogram signals corresponding to each second pulse on the plurality of channels, wherein the plurality of second electroencephalogram signals are in one-to-one correspondence with the plurality of channels; the evoked potential online analysis module is used for analyzing effective response components of a plurality of second brain electrical signals corresponding to each second pulse and determining whether effective response components exist in the jth transcranial magnetic test.
The manner of delivering the second pulse by the transcranial magnetic stimulation module is similar to that of delivering the first pulse, and will not be described again; the manner in which the electroencephalogram acquisition module acquires the second electroencephalogram is similar to that in which the first electroencephalogram is acquired, and is not described again; similarly, the evoked potential in-line analysis module determines whether an active response component is present in the jth transcranial magnetic test in a manner similar to that described above for determining whether an active response component is present in the ith transcranial magnetic test, and is not described again.
Further, if no effective response component exists, the stimulation intensity control module determines the stimulation intensity dj+1 of the j+1th transcranial magnetic test based on the stimulation intensity Dj to perform the j+1th transcranial magnetic test until the effective response component exists, i.e. increases the stimulation intensity again, and determines whether the effective response component exists or not until the effective response component exists. When an effective response component is present, the stimulus intensity control module takes the stimulus intensity Dj as the target stimulus intensity.
It can be seen that in the embodiment of the present application, before the transcranial magnetic stimulation treatment is performed on the subject, the subject (i.e., the patient) is subjected to multiple transcranial magnetic tests in the area to be treated of the subject, and when the transcranial magnetic tests are performed multiple times, the stimulus intensity is gradually reduced (i.e., the stimulus intensity Si is smaller than the stimulus intensity Si-1), and the analysis of the effective response component is performed based on the first electroencephalogram signal of each transcranial magnetic test until the stimulus intensity without the effective response component is found, so that the target stimulus intensity suitable for the subject should be greater than the stimulus intensity because the patient needs to respond to achieve the treatment during the transcranial magnetic stimulation treatment. Thus, in the present application, the stimulus intensity is gradually increased from the stimulus intensity (i.e., the stimulus intensity Dj is greater than the stimulus intensity Dj-1), and the analysis of the effective response component is performed based on the second electroencephalogram signal of each transcranial magnetic test until the effective response component is again generated, that is, the stimulus intensity just enables the patient to respond, and the stimulus intensity is taken as the target stimulus intensity. Therefore, through the technical scheme of the application, the set stimulation intensity can ensure that the patient can be stimulated to respond, and the stimulation intensity can not be excessively high, so that the stimulation intensity matched with the area to be treated of the patient can be set with high precision, and the transcranial magnetic stimulation treatment can be carried out on the patient by using the target stimulation intensity, and the treatment effect can be improved. In addition, the application tests the region to be treated practically, and eliminates the influence of cortex injury or brain region difference on the setting of the stimulation intensity in the traditional method. The application realizes complete automatic detection, automatic analysis and automatic control, eliminates artificial influence, does not need training for operators, does not need intensity setting for experienced personnel, and reduces labor cost.
Optionally, after determining the target stimulation intensity, the target stimulation intensity may be used to pulse the subject for transcranial magnetic stimulation treatment of the subject while the subject is being subjected to transcranial magnetic stimulation treatment.
Optionally, the transcranial magnetic stimulation device further comprises a display module, and after the stimulation intensity control module determines the target stimulation intensity, the target stimulation intensity can be displayed through the display module, so that a doctor or a professional can look over the target stimulation intensity.
It should be noted that, the electroencephalogram acquisition module, the evoked potential online analysis module, the transcranial magnetic stimulation module and the stimulation intensity control module may be actual functional modules or logically divided functional modules, which is not limited in the present application.
Referring to fig. 6, fig. 6 is a flow chart of a method for setting stimulus intensity according to an embodiment of the application. The method is applied to the transcranial stimulation device. The method includes, but is not limited to, the following steps:
s601: the stimulation intensity control module determines the stimulation intensity Si of the ith transcranial magnetic test based on the stimulation intensity Si-1 of the ith-1 th transcranial magnetic test.
S602: the transcranial magnetic stimulation module delivers a plurality of first pulses to the subject corresponding to the stimulation intensity Si.
S603: the electroencephalogram acquisition module acquires a plurality of first electroencephalograms corresponding to each first pulse on a plurality of channels, wherein the plurality of first electroencephalograms correspond to the plurality of channels one by one.
S604: and the evoked potential online analysis module analyzes the effective response components of a plurality of first electroencephalogram signals corresponding to each first pulse and determines whether the effective response components exist in the ith transcranial magnetic test.
S605: if the effective response component exists, the stimulation intensity control module determines the stimulation intensity Si+1 of the ith+1st transcranial magnetic test based on the stimulation intensity Si, so as to perform the ith+1st transcranial magnetic test until the effective response component does not exist.
S606: if no effective response component exists, the stimulation intensity control module determines the stimulation intensity Dj of the jth transcranial magnetic test based on the stimulation intensity Dj-1 of the jth transcranial magnetic test, wherein the stimulation intensity Dj is larger than the stimulation intensity Dj-1, and when j=1, the stimulation intensity D1 is the stimulation intensity Si.
S607: the transcranial magnetic stimulation module delivers a plurality of second pulses to the subject corresponding to the stimulation intensity Dj.
S608: the electroencephalogram acquisition module acquires a plurality of second electroencephalograms corresponding to each second pulse on the channels, wherein the second electroencephalograms correspond to the channels one by one.
S609: and the evoked potential online analysis module analyzes the effective response components of a plurality of second brain electrical signals corresponding to each second pulse and determines whether the effective response components exist in the jth transcranial magnetic test.
S610: if no effective response component exists, the stimulation intensity control module determines the stimulation intensity Dj+1 of the j+1th transcranial magnetic test based on the stimulation intensity Dj so as to perform the j+1th transcranial magnetic test until the effective response component exists; if a valid response component exists, the stimulus intensity control module takes the stimulus intensity Dj as the target stimulus intensity.
The specific implementation process of step S601 to step S610 may refer to specific functions of the electroencephalogram acquisition module, the evoked potential online analysis module, the transcranial magnetic stimulation module and the stimulation intensity control module, and will not be described again.
Referring to fig. 7, fig. 7 is a schematic diagram of an electronic device according to an embodiment of the application. The electronic device 700 shown in fig. 7 may be a transcranial magnetic stimulation device as described above. The electronic device 700 shown in fig. 7 includes a memory 701, a processor 702, a communication interface 703, and a bus 704. The memory 701, the processor 702, and the communication interface 703 are connected to each other by a bus 704.
The processor 702 may integrate the functions of the evoked potential online analysis module and the stimulus intensity control module, which are used for calculating the stimulus intensity and analyzing the effective response components. The communication interface 703 may integrate the functions of the stimulation intensity control module and the electroencephalogram acquisition module described above for issuing pulses to a subject and acquiring electroencephalogram signals from the subject.
The Memory 701 may be a Read Only Memory (ROM), a static storage device, a dynamic storage device, or a random access Memory (Random Access Memory, RAM). The memory 701 may store a program, and when the program stored in the memory 701 is executed by the processor 702, the processor 702 and the communication interface 703 are used to perform the respective steps of the stimulation intensity setting method of the embodiment of the present application.
The processor 702 may employ a general-purpose central processing unit (Central Processing Unit, CPU), microprocessor, application specific integrated circuit (Application Specific Integrated Circuit, ASIC), graphics processor (graphics processing unit, GPU) or one or more integrated circuits for executing associated programs to perform the functions required to be performed by the units in the transcranial magnetic stimulation device or to perform the stimulation intensity setting method of an embodiment of the method of the present application.
The processor 702 may also be an integrated circuit chip with signal processing capabilities. In implementation, various steps in the stimulus intensity setting method of the present application may be accomplished by instructions in the form of integrated logic circuits of hardware or software in the processor 702. The processor 702 may also be a general purpose processor, a digital signal processor (Digital Signal Processing, DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf programmable gate array (Field Programmable Gate Array, FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware components. The disclosed methods, steps, and logic blocks in the embodiments of the present application may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present application may be embodied directly in the execution of a hardware decoding processor, or in the execution of a combination of hardware and software modules in a decoding processor. The software modules may be located in a random access memory, flash memory, read only memory, programmable read only memory, or electrically erasable programmable memory, registers, etc. as well known in the art. The storage medium is located in the memory 701, and the processor 702 reads the information in the memory 701, and in combination with its hardware, performs the functions that the units included in the transcranial magnetic stimulation device according to the embodiment of the present application need to perform, or performs the steps in the stimulation intensity setting method according to the embodiment of the present application.
The communication interface 703 enables communication between the electronic device 700 and other devices or communication networks using transceiving means such as, but not limited to, transceivers, input-output devices, and the like. For example, the first pulse may be acquired through the communication interface 703. Specifically, when the communication interface 703 is an input device, a first pulse may be issued to the subject for the transcranial magnetic stimulation module described above; when the communication interface 703 is an output device, it may be the above-mentioned electroencephalogram acquisition module, and a first electroencephalogram signal is acquired from the subject.
A bus 704 may include a path that communicates information between various components of the electronic device 700 (e.g., memory 701, processor 702, communication interface 703).
It should be noted that while the electronic device 700 shown in fig. 7 illustrates only a memory, a processor, and a communication interface, those skilled in the art will appreciate that in a particular implementation, the electronic device 700 also includes other components necessary to achieve proper operation. Also, those skilled in the art will appreciate that the electronic device 700 may also include hardware components that perform other additional functions, as desired. Furthermore, it will be appreciated by those skilled in the art that the electronic device 700 may also include only the components necessary to implement embodiments of the present application, and not necessarily all of the components shown in FIG. 7.
It will be clear to those skilled in the art that, for convenience and brevity of description, specific working procedures of the above-described systems, apparatuses and units may refer to corresponding procedures in the foregoing method embodiments, and are not repeated herein.
In the several embodiments provided by the present application, it should be understood that the disclosed systems, devices, and methods may be implemented in other manners. For example, the apparatus embodiments described above are merely illustrative, e.g., the division of elements is merely a logical functional division, and there may be additional divisions of actual implementation, e.g., multiple elements or components may be combined or integrated into another system, or some features may be omitted, or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or units, which may be in electrical, mechanical or other form.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed over a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present application may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, comprising several instructions for causing a computer device (which may be a personal computer, a server, a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
The foregoing is merely illustrative of the present application, and the present application is not limited thereto, and any person skilled in the art will readily recognize that variations or substitutions are within the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.
The embodiment of the present application also provides a computer-readable storage medium storing a computer program that is executed by a processor to implement part or all of the steps of any one of the stimulus intensity setting methods described in the above method embodiments.
Embodiments of the present application also provide a computer program product comprising a non-transitory computer-readable storage medium storing a computer program operable to cause a computer to perform part or all of the steps of any one of the stimulation intensity setting methods described in the method embodiments above.
It should be noted that, for simplicity of description, the foregoing method embodiments are all described as a series of acts, but it should be understood by those skilled in the art that the present application is not limited by the order of acts described, as some steps may be performed in other orders or concurrently in accordance with the present application. Further, those skilled in the art will also appreciate that the embodiments described in the specification are alternative embodiments, and that the acts and modules referred to are not necessarily required for the present application.
In the foregoing embodiments, the descriptions of the embodiments are emphasized, and for parts of one embodiment that are not described in detail, reference may be made to related descriptions of other embodiments.
In the several embodiments provided by the present application, it should be understood that the disclosed apparatus may be implemented in other manners. For example, the apparatus embodiments described above are merely illustrative, such as the division of the units, merely a logical function division, and there may be additional manners of dividing the actual implementation, such as multiple units or components may be combined or integrated into another system, or some features may be omitted, or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or units, or may be in electrical or other forms.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present application may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units described above may be implemented either in hardware or in software program modules.
The integrated units, if implemented in the form of software program modules, may be stored in a computer-readable memory for sale or use as a stand-alone product. Based on this understanding, the technical solution of the present application may be embodied essentially or partly in the form of a software product, or all or part of the technical solution, which is stored in a memory, and includes several instructions for causing a computer device (which may be a personal computer, a server, a network device, or the like) to perform all or part of the steps of the method according to the embodiments of the present application. And the aforementioned memory includes: a U-disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a removable hard disk, a magnetic disk, or an optical disk, or other various media capable of storing program codes.
Those of ordinary skill in the art will appreciate that all or a portion of the steps in the various methods of the above embodiments may be implemented by a program that instructs associated hardware, and the program may be stored in a computer readable memory, which may include: flash disk, read-Only Memory (ROM), random access Memory (Random Access Memory, RAM), magnetic disk or optical disk.
The foregoing has outlined rather broadly the more detailed description of embodiments of the application, wherein the principles and embodiments of the application are explained in detail using specific examples, the above examples being provided solely to facilitate the understanding of the method and core concepts of the application; meanwhile, as those skilled in the art will have variations in the specific embodiments and application scope in accordance with the ideas of the present application, the present description should not be construed as limiting the present application in view of the above.

Claims (5)

1. The transcranial magnetic stimulation equipment is characterized by comprising an electroencephalogram acquisition module, an evoked potential on-line analysis module, a transcranial magnetic stimulation module and a stimulation intensity control module;
The stimulation intensity control module is used for determining the stimulation intensity Si of the ith transcranial magnetic test based on the stimulation intensity Si-1 of the ith-1 th transcranial magnetic test, wherein the stimulation intensity Si is smaller than the stimulation intensity Si-1; si= (K) i-1 - δk) ×smax, wherein δk is a preset scaling factor, K i-1 For the i-1 th transcranial magnetic test, smax is the coefficient of proportionality that the transcranial magnetic stimulation device can deliverMaximum stimulus intensity;
the transcranial magnetic stimulation module is used for sending a plurality of first pulses corresponding to the stimulation intensity Si to a subject;
the electroencephalogram acquisition module is used for acquiring a plurality of first electroencephalograms corresponding to each first pulse on a plurality of channels, wherein the plurality of first electroencephalograms correspond to the plurality of channels one by one;
the evoked potential online analysis module is used for carrying out effective response component analysis on a plurality of first electroencephalogram signals corresponding to each first pulse and determining whether effective response components exist in the ith transcranial magnetic test; the method is particularly used for: based on the moment of issuing each first pulse, respectively interpolating a plurality of first electroencephalogram signals of each first pulse on the channels to obtain a plurality of first test times of each first pulse on the channels; performing noise fitting on all first test times corresponding to the first pulses on the channels to obtain a noise curve, wherein the noise curve is used for representing a change rule of noise when electroencephalogram signal acquisition is performed after the first pulses are issued; performing component analysis on all first test times corresponding to the plurality of first pulses on the plurality of channels to obtain a plurality of component curves corresponding to a plurality of components, wherein the plurality of component curves are in one-to-one correspondence with the plurality of components, and each component curve is used for representing a change rule of the component corresponding to the component curve after the plurality of first pulses are issued; based on the noise curve and the component curves, removing noise from the first test times corresponding to each first pulse to obtain second test times corresponding to each first pulse;
For each second test corresponding to each first pulse, acquiring the average value of the amplitude values of the second test in a preset time period; subtracting the average value from the amplitude of the second test at each time to obtain a plurality of third test times of each first pulse on the channels; averaging the amplitudes of the third test times of each channel under the first pulses at all moments to obtain a fourth test time corresponding to each channel; determining a first covariance matrix according to a plurality of fourth test times corresponding to the plurality of channels; determining a second covariance matrix corresponding to each first pulse according to a plurality of third trials of each first pulse under the plurality of channels; filtering the first pulses according to a second covariance matrix corresponding to each first pulse and the first covariance matrix to obtain a plurality of target pulses;
taking the time of issuing each target pulse as a dividing line, and carrying out signal division on the third test time of each target pulse on each channel to obtain a baseline signal and an induced signal of each target pulse on each channel; performing Gaussian fitting on the basis of the amplitude of the baseline signal of each target pulse on each channel to obtain an amplitude threshold; screening a plurality of target channels from the plurality of channels; averaging the amplitudes of the induced signals of each target pulse on the plurality of target channels at each time to obtain an average amplitude corresponding to each target pulse; averaging a plurality of average amplitudes corresponding to the target pulses to obtain target amplitudes; if the target amplitude is greater than or equal to the amplitude threshold, determining that an effective response component exists in the ith transcranial magnetic test; if the target amplitude is smaller than the amplitude threshold, determining that no effective response component exists in the ith transcranial magnetic test;
The stimulation intensity control module is further used for determining the stimulation intensity Si+1 of the ith+1th transcranial magnetic test based on the stimulation intensity Si if the effective response component exists, so as to perform the ith+1th transcranial magnetic test until the effective response component does not exist;
the stimulation intensity control module is further used for determining the stimulation intensity Dj of the jth transcranial magnetic test based on the stimulation intensity Dj-1 of the jth transcranial magnetic test if no effective response component exists, wherein the stimulation intensity Dj is larger than the stimulation intensity Dj-1, and when j=1, the stimulation intensity D1 is the stimulation intensity Si; dj= (R j-1 +δr). Times.Smax, wherein R j-1 For the j-i th transcranial magnetic testing, δr is a preset scaling factor, and δr is smaller than δk;
the transcranial magnetic stimulation module is further configured to deliver a plurality of second pulses to the subject corresponding to the stimulation intensity Dj;
the electroencephalogram acquisition module is further used for acquiring a plurality of second electroencephalograms corresponding to each second pulse on the plurality of channels, wherein the plurality of second electroencephalograms correspond to the plurality of channels one by one;
the evoked potential online analysis module is further used for carrying out effective response component analysis on a plurality of second brain electrical signals corresponding to each second pulse and determining whether effective response components exist in the jth transcranial magnetic test;
The stimulation intensity control module is further used for determining the stimulation intensity Dj+1 of the j+1th transcranial magnetic test based on the stimulation intensity Dj if no effective response component exists, so as to perform the j+1th transcranial magnetic test until the effective response component exists; if an effective response component exists, the stimulus intensity Dj is taken as a target stimulus intensity.
2. The apparatus of claim 1, wherein the device comprises a plurality of sensors,
in the aspect of removing noise from the plurality of first test times corresponding to each first pulse based on the noise curve and the plurality of component curves to obtain a plurality of second test times corresponding to each first pulse, the evoked potential online analysis module is specifically configured to:
fitting each component curve with the noise curve to obtain a root mean square error corresponding to each component curve;
setting weights of components corresponding to component curves with root mean square errors smaller than a threshold value to 0 in all first test times to obtain weights of the components in each first test time corresponding to each first pulse;
and performing test reconstruction based on the weight of the components under each first test corresponding to each first pulse to obtain a plurality of second tests on the channels corresponding to each first pulse.
3. The apparatus according to claim 1 or 2, wherein,
in terms of performing gaussian fitting based on the amplitude of the baseline signal of each target pulse on each channel to obtain an amplitude threshold, the evoked potential online analysis module is specifically configured to:
determining a mean value and a standard deviation of the amplitude of each target pulse on each channel based on the amplitude of the baseline signal of each target pulse on each channel at each time instant;
determining a target mean value and a target standard deviation based on the mean and standard deviation of the amplitude of each target pulse on each channel;
constructing Gaussian distribution corresponding to the target mean value and the target standard deviation;
and determining the amplitude threshold based on the Gaussian distribution and a preset confidence.
4. An electronic device, comprising: the electronic device comprises a processor and a memory, wherein the processor is connected with the memory, the memory is used for storing a computer program, and the processor is used for executing the computer program stored in the memory so as to enable the electronic device to execute the following steps:
determining a stimulus intensity Si of the ith transcranial magnetic test based on the stimulus intensity Si-1 of the ith transcranial magnetic test, wherein the stimulus intensity Si is less than the stimulus intensity Si-1; si= (K) i-1 - δk) ×smax, wherein δk is a preset scaling factor, K i-1 For the i-1 th transcranial magnetic test, smax is the maximum stimulation intensity which can be given out by transcranial magnetic stimulation equipment;
issuing a plurality of first pulses corresponding to the stimulus intensity Si to a subject;
collecting a plurality of first electroencephalograms corresponding to each first pulse on a plurality of channels, wherein the plurality of first electroencephalograms are in one-to-one correspondence with the plurality of channels;
performing effective response component analysis on a plurality of first electroencephalogram signals corresponding to each first pulse, and determining whether effective response components exist in the ith transcranial magnetic test; comprising the following steps: based on the moment of issuing each first pulse, respectively interpolating a plurality of first electroencephalogram signals of each first pulse on the channels to obtain a plurality of first test times of each first pulse on the channels; performing noise fitting on all first test times corresponding to the first pulses on the channels to obtain a noise curve, wherein the noise curve is used for representing a change rule of noise when electroencephalogram signal acquisition is performed after the first pulses are issued; performing component analysis on all first test times corresponding to the plurality of first pulses on the plurality of channels to obtain a plurality of component curves corresponding to a plurality of components, wherein the plurality of component curves are in one-to-one correspondence with the plurality of components, and each component curve is used for representing a change rule of the component corresponding to the component curve after the plurality of first pulses are issued; based on the noise curve and the component curves, removing noise from the first test times corresponding to each first pulse to obtain second test times corresponding to each first pulse;
For each second test corresponding to each first pulse, acquiring the average value of the amplitude values of the second test in a preset time period; subtracting the average value from the amplitude of the second test at each time to obtain a plurality of third test times of each first pulse on the channels; averaging the amplitudes of the third test times of each channel under the first pulses at all moments to obtain a fourth test time corresponding to each channel; determining a first covariance matrix according to a plurality of fourth test times corresponding to the plurality of channels; determining a second covariance matrix corresponding to each first pulse according to a plurality of third trials of each first pulse under the plurality of channels; filtering the first pulses according to a second covariance matrix corresponding to each first pulse and the first covariance matrix to obtain a plurality of target pulses;
taking the time of issuing each target pulse as a dividing line, and carrying out signal division on the third test time of each target pulse on each channel to obtain a baseline signal and an induced signal of each target pulse on each channel; performing Gaussian fitting on the basis of the amplitude of the baseline signal of each target pulse on each channel to obtain an amplitude threshold; screening a plurality of target channels from the plurality of channels; averaging the amplitudes of the induced signals of each target pulse on the plurality of target channels at each time to obtain an average amplitude corresponding to each target pulse; averaging a plurality of average amplitudes corresponding to the target pulses to obtain target amplitudes; if the target amplitude is greater than or equal to the amplitude threshold, determining that an effective response component exists in the ith transcranial magnetic test; if the target amplitude is smaller than the amplitude threshold, determining that no effective response component exists in the ith transcranial magnetic test;
If the effective response component exists, determining the stimulation intensity Si+1 of the ith+1st transcranial magnetic test based on the stimulation intensity Si, so as to perform the ith+1st transcranial magnetic test until the effective response component does not exist;
if no effective response component exists, determining the stimulation intensity Dj of the jth transcranial magnetic test based on the stimulation intensity Dj-1 of the jth transcranial magnetic test, wherein the stimulation intensity Dj is larger than the stimulation intensity Dj-1, and when j=1, the stimulation intensity D1 is the stimulation intensity Si; dj= (R j-1 +δr). Times.Smax, wherein R j-1 For the j-i th transcranial magnetic testing, δr is a preset scaling factor, and δr is smaller than δk;
issuing a plurality of second pulses to the subject corresponding to the stimulus intensity Dj;
collecting a plurality of second electroencephalogram signals corresponding to each second pulse on the channels, wherein the second electroencephalogram signals are in one-to-one correspondence with the channels;
performing effective response component analysis on a plurality of second electroencephalogram signals corresponding to each second pulse, and determining whether effective response components exist in the jth transcranial magnetic test;
if the effective response component does not exist, determining the stimulation intensity Dj+1 of the j+1th transcranial magnetic test based on the stimulation intensity Dj so as to perform the j+1th transcranial magnetic test until the effective response component exists; if an effective response component exists, the stimulus intensity Dj is taken as a target stimulus intensity.
5. A computer readable storage medium storing a computer program, the computer program being executable by a processor to perform the steps of:
determining a stimulus intensity Si of the ith transcranial magnetic test based on the stimulus intensity Si-1 of the ith transcranial magnetic test, wherein the stimulus intensity Si is less than the stimulus intensity Si-1; si= (K) i-1 - δk) ×smax, wherein δk is a preset scaling factor, K i-1 For the i-1 th transcranial magnetic test, smax is the maximum stimulation intensity which can be given out by transcranial magnetic stimulation equipment;
issuing a plurality of first pulses corresponding to the stimulus intensity Si to a subject;
collecting a plurality of first electroencephalograms corresponding to each first pulse on a plurality of channels, wherein the plurality of first electroencephalograms are in one-to-one correspondence with the plurality of channels;
performing effective response component analysis on a plurality of first electroencephalogram signals corresponding to each first pulse, and determining whether effective response components exist in the ith transcranial magnetic test; comprising the following steps: based on the moment of issuing each first pulse, respectively interpolating a plurality of first electroencephalogram signals of each first pulse on the channels to obtain a plurality of first test times of each first pulse on the channels; performing noise fitting on all first test times corresponding to the first pulses on the channels to obtain a noise curve, wherein the noise curve is used for representing a change rule of noise when electroencephalogram signal acquisition is performed after the first pulses are issued; performing component analysis on all first test times corresponding to the plurality of first pulses on the plurality of channels to obtain a plurality of component curves corresponding to a plurality of components, wherein the plurality of component curves are in one-to-one correspondence with the plurality of components, and each component curve is used for representing a change rule of the component corresponding to the component curve after the plurality of first pulses are issued; based on the noise curve and the component curves, removing noise from the first test times corresponding to each first pulse to obtain second test times corresponding to each first pulse;
For each second test corresponding to each first pulse, acquiring the average value of the amplitude values of the second test in a preset time period; subtracting the average value from the amplitude of the second test at each time to obtain a plurality of third test times of each first pulse on the channels; averaging the amplitudes of the third test times of each channel under the first pulses at all moments to obtain a fourth test time corresponding to each channel; determining a first covariance matrix according to a plurality of fourth test times corresponding to the plurality of channels; determining a second covariance matrix corresponding to each first pulse according to a plurality of third trials of each first pulse under the plurality of channels; filtering the first pulses according to a second covariance matrix corresponding to each first pulse and the first covariance matrix to obtain a plurality of target pulses;
taking the time of issuing each target pulse as a dividing line, and carrying out signal division on the third test time of each target pulse on each channel to obtain a baseline signal and an induced signal of each target pulse on each channel; performing Gaussian fitting on the basis of the amplitude of the baseline signal of each target pulse on each channel to obtain an amplitude threshold; screening a plurality of target channels from the plurality of channels; averaging the amplitudes of the induced signals of each target pulse on the plurality of target channels at each time to obtain an average amplitude corresponding to each target pulse; averaging a plurality of average amplitudes corresponding to the target pulses to obtain target amplitudes; if the target amplitude is greater than or equal to the amplitude threshold, determining that an effective response component exists in the ith transcranial magnetic test; if the target amplitude is smaller than the amplitude threshold, determining that no effective response component exists in the ith transcranial magnetic test;
If the effective response component exists, determining the stimulation intensity Si+1 of the ith+1st transcranial magnetic test based on the stimulation intensity Si, so as to perform the ith+1st transcranial magnetic test until the effective response component does not exist;
if no effective response component exists, the method is based on the j-1 th transcranial magnetic measurementDetermining the stimulation intensity Dj of the jth transcranial magnetic test according to the tested stimulation intensity Dj-1, wherein the stimulation intensity Dj is larger than the stimulation intensity Dj-1, and when j=1, the stimulation intensity D1 is the stimulation intensity Si; dj= (R j-1 +δr). Times.Smax, wherein R j-1 For the j-i th transcranial magnetic testing, δr is a preset scaling factor, and δr is smaller than δk;
issuing a plurality of second pulses to the subject corresponding to the stimulus intensity Dj;
collecting a plurality of second electroencephalogram signals corresponding to each second pulse on the channels, wherein the second electroencephalogram signals are in one-to-one correspondence with the channels;
performing effective response component analysis on a plurality of second electroencephalogram signals corresponding to each second pulse, and determining whether effective response components exist in the jth transcranial magnetic test;
if the effective response component does not exist, determining the stimulation intensity Dj+1 of the j+1th transcranial magnetic test based on the stimulation intensity Dj so as to perform the j+1th transcranial magnetic test until the effective response component exists; if an effective response component exists, the stimulus intensity Dj is taken as a target stimulus intensity.
CN202310775449.4A 2023-06-28 2023-06-28 Transcranial magnetic stimulation device, electronic device, and storage medium Active CN116510183B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310775449.4A CN116510183B (en) 2023-06-28 2023-06-28 Transcranial magnetic stimulation device, electronic device, and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310775449.4A CN116510183B (en) 2023-06-28 2023-06-28 Transcranial magnetic stimulation device, electronic device, and storage medium

Publications (2)

Publication Number Publication Date
CN116510183A CN116510183A (en) 2023-08-01
CN116510183B true CN116510183B (en) 2023-09-19

Family

ID=87397939

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310775449.4A Active CN116510183B (en) 2023-06-28 2023-06-28 Transcranial magnetic stimulation device, electronic device, and storage medium

Country Status (1)

Country Link
CN (1) CN116510183B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116712673A (en) * 2023-08-11 2023-09-08 江西省精一医疗科技有限公司 Median nerve stimulation system, electronic device, and storage medium

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101234233A (en) * 2008-01-29 2008-08-06 武汉依瑞德医疗设备新技术有限公司 Transcranial magnetic field stimulator stimulus strength modulating method and device therefor
CN208511024U (en) * 2017-08-29 2019-02-19 南京伟思医疗科技股份有限公司 A kind of magnetic induction triggering myoelectricity collector suitable for transcranial magnetic stimulation instrument
CN110114115A (en) * 2016-10-25 2019-08-09 脑路有限公司 Device and method for predicted treatment result
CN113545792A (en) * 2021-07-20 2021-10-26 深圳市人民医院 TMS-EEG-based neural biomarker detection method
CN113827865A (en) * 2021-10-26 2021-12-24 中国医学科学院生物医学工程研究所 Transcranial magnetic stimulation position determination method and device, electronic equipment and storage medium
CN114041806A (en) * 2021-11-09 2022-02-15 电子科技大学 Self-adaptive filtering TMS evoked potential extraction device
CN114917479A (en) * 2022-06-30 2022-08-19 深圳英智科技有限公司 Transcranial magnetic pulse interference filtering method and system and electronic equipment
CN115212456A (en) * 2022-07-07 2022-10-21 四川脑科学与类脑智能研究院 Device and method for improving sleep through transcranial real-time asynchronous electrical stimulation
CN115364377A (en) * 2022-08-11 2022-11-22 江西杰联医疗设备有限公司 Frontal lobe positioning method of transcranial magnetic stimulation coil, closed-loop transcranial magnetic stimulation system and device
WO2023066020A1 (en) * 2021-10-19 2023-04-27 南京伟思医疗科技股份有限公司 Therapeutic target automatic positioning method and system for transcranial magnetism, and helmet
CN116328192A (en) * 2023-03-22 2023-06-27 中国医学科学院生物医学工程研究所 Transcranial magnetic stimulation system and method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10112056B2 (en) * 2012-05-25 2018-10-30 Monash University Optimising current direction and intensity of transcranial magnetic stimulation

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101234233A (en) * 2008-01-29 2008-08-06 武汉依瑞德医疗设备新技术有限公司 Transcranial magnetic field stimulator stimulus strength modulating method and device therefor
CN110114115A (en) * 2016-10-25 2019-08-09 脑路有限公司 Device and method for predicted treatment result
CN208511024U (en) * 2017-08-29 2019-02-19 南京伟思医疗科技股份有限公司 A kind of magnetic induction triggering myoelectricity collector suitable for transcranial magnetic stimulation instrument
CN113545792A (en) * 2021-07-20 2021-10-26 深圳市人民医院 TMS-EEG-based neural biomarker detection method
WO2023066020A1 (en) * 2021-10-19 2023-04-27 南京伟思医疗科技股份有限公司 Therapeutic target automatic positioning method and system for transcranial magnetism, and helmet
CN113827865A (en) * 2021-10-26 2021-12-24 中国医学科学院生物医学工程研究所 Transcranial magnetic stimulation position determination method and device, electronic equipment and storage medium
CN114041806A (en) * 2021-11-09 2022-02-15 电子科技大学 Self-adaptive filtering TMS evoked potential extraction device
CN114917479A (en) * 2022-06-30 2022-08-19 深圳英智科技有限公司 Transcranial magnetic pulse interference filtering method and system and electronic equipment
CN115212456A (en) * 2022-07-07 2022-10-21 四川脑科学与类脑智能研究院 Device and method for improving sleep through transcranial real-time asynchronous electrical stimulation
CN115364377A (en) * 2022-08-11 2022-11-22 江西杰联医疗设备有限公司 Frontal lobe positioning method of transcranial magnetic stimulation coil, closed-loop transcranial magnetic stimulation system and device
CN116328192A (en) * 2023-03-22 2023-06-27 中国医学科学院生物医学工程研究所 Transcranial magnetic stimulation system and method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Comparison of synchrosqueezing transform to alternative methods for time-frequency analysis of TMS-evoked EEG oscillations;Yong Wang et al.;《Biomedical Signal Processing and Control 》;1-9 *
不同强度低频重复经颅磁刺激对帕金森大鼠行为学的影响;何晓阔 等;《中国临床康复》;第9卷(第41期);11-13 *
基于脑电的意识障碍重复经颅磁刺激调控评估;白洋 等;《中国生物医学工程学报》;第38卷(第6期);687-694 *
意识障碍的脑电分析与神经调控评估;白洋;《中国博士学位论文全文数据库》;1-148 *

Also Published As

Publication number Publication date
CN116510183A (en) 2023-08-01

Similar Documents

Publication Publication Date Title
Horvath et al. Evidence that transcranial direct current stimulation (tDCS) generates little-to-no reliable neurophysiologic effect beyond MEP amplitude modulation in healthy human subjects: a systematic review
CN101677775B (en) System and method for pain detection and computation of a pain quantification index
Casula et al. Evidence for interhemispheric imbalance in stroke patients as revealed by combining transcranial magnetic stimulation and electroencephalography
US6195576B1 (en) Quantitative magnetoencephalogram system and method
US9826916B2 (en) Device and method for examining a phase distribution used to determine a pathological interaction between different areas of the brain
US8014870B2 (en) Method and apparatus for the treatment of tinnitus
Baur et al. Induction of LTD-like corticospinal plasticity by low-frequency rTMS depends on pre-stimulus phase of sensorimotor μ-rhythm
WO2017189757A1 (en) Brain stimulation treatment in depression
WO2004100765A2 (en) System and method of assessment of the efficacy of treatment of neurological disorders using the electroencephalogram
Peters et al. On the feasibility of concurrent human TMS-EEG-fMRI measurements
US20200281517A1 (en) Stimulative Electrotherapy Using Autonomic Nervous System Control
CN116510183B (en) Transcranial magnetic stimulation device, electronic device, and storage medium
Rich et al. Determining electrode placement for transcranial direct current stimulation: a comparison of EEG-versus TMS-guided methods
US20180256912A1 (en) Method for monitoring treatment of neuropsychiatric disorders
Wu et al. Correlated disruption of resting-state fMRI, LFP, and spike connectivity between area 3b and S2 following spinal cord injury in monkeys
Bashir et al. The number of pulses needed to measure corticospinal excitability by navigated transcranial magnetic stimulation: eyes open vs. close condition
Piitulainen et al. Spatial variability in cortex-muscle coherence investigated with magnetoencephalography and high-density surface electromyography
Nguyen et al. Principal component regression on motor evoked potential in single-pulse transcranial magnetic stimulation
Zazio et al. Perceptual and physiological consequences of dark adaptation: a TMS-EEG study
Livingston et al. Motor evoked potential response latencies demonstrate moderate correlations with height and limb length in healthy young adults
Vink et al. EEG functional connectivity is a weak predictor of causal brain interactions
Zhu et al. The relationship between magnetic and electrophysiological responses to complex tactile stimuli
Kasten et al. Integrating electric field modelling and neuroimaging to explain inter-individual variability of tACS effects
Wei et al. [Retracted] Simulation Study on the Spatiotemporal Difference of Complex Neurodynamics between P3a and P3b
Schmid et al. Simultaneous EEG and fMRI in the macaque monkey at 4.7 Tesla

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant