CN103494609A - Method for indicating electromyographic signal interference and method for removing interference - Google Patents
Method for indicating electromyographic signal interference and method for removing interference Download PDFInfo
- Publication number
- CN103494609A CN103494609A CN201310486036.0A CN201310486036A CN103494609A CN 103494609 A CN103494609 A CN 103494609A CN 201310486036 A CN201310486036 A CN 201310486036A CN 103494609 A CN103494609 A CN 103494609A
- Authority
- CN
- China
- Prior art keywords
- signal
- electromyographic
- electromyographic signal
- annoyance level
- disturbed
- 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.)
- Granted
Links
Images
Abstract
The invention discloses a method for indicating electromyographic signal interference. The method comprises the following steps that (1) in a time frame t, the sampling value sequence x (t) of electromyographic signals is obtained, and t=1-10s; (2) the power spectrum P (f) and the total power Tp of the electromyographic signals are calculated; (3) in the bandwidth range of the electromyographic signals, the power values P50, P (n*50) of 50Hz and the frequency doubling (n*50) Hz are determined, and n is a natural number larger than 1; (4) the percent K of the sum of the power values of 50Hz and the frequency doubling (n*50) Hz accounted for in the total power Tp is calculated, and Q=100-K; (5) the value of K or the value of Q is indicated for a user. The invention further discloses a method for removing the electromyographic signal interference according to the indicating method. The method for indicating the electromyographic signal interference helps a doctor to judge the reliability of the signals and provides a basis for interference removal, and therefore the reliable signals can be obtained under the condition of removing 50Hz and the harmonic energy as less as possible.
Description
Technical field
The present invention relates to the indicating means that a kind of electromyographic signal is disturbed, for collection, analysis and the biofeedback training equipment of the human body electromyographic signal of medical apparatus and instruments.
Background technology
Due to the human biological signal, very faint, such as, between EEG signals amplitude 5uV-100uV, electrocardiosignal is between 10uV-10mV, between electromyographic signal 2uV-10mV, all these bioelectrical signals all easily are subject to the interference of external electromagnetic ripple, the particularly interference of 50Hz power frequency and harmonic wave thereof.The try every possible means capacity of resisting disturbance of raising equipment of equipment vendors, such as the common mode rejection ratio that improves equipment, increase 50Hz notch filter circuit.But the interference signal intensity under varying environment is different, still can't break away from the impact of interference.Electromyographic signal particularly, the broader bandwidth of signal, the surface electromyogram signal bandwidth is at 20Hz-500Hz, pin electromyographic signal bandwidth is at 10Hz-10KHz, 50Hz and harmonic wave thereof are all in the bandwidth range of signal, more easily be interfered, because also there is useful signal at 50Hz and harmonic wave place thereof, so should not pass through the method for trap 50Hz and the whole filterings of harmonic energy thereof yet, but in some cases, in order to improve the capacity of resisting disturbance of equipment, can carry out the trap of 50Hz, and seldom the harmonic wave of 50Hz be processed.
From primary signal and RMS value signal, the doctor is difficult to judgement, the signal collected be by patient contract muscles cause or because interference causes.Due in the equipment such as muscle electricity biofeedback instrument, seldom there is the signal spectral analysis function, be all to carry out feedback training with signal amplitude, the doctor has no idea to distinguish at all, and signal is muscle contraction signal or interfering signal.Even in thering is the equipment of spectrum analysis, due to the non-industry science background of doctor, sometimes also be difficult to distinguish interference by spectrogram.
Summary of the invention
Goal of the invention: the object of the invention is to for the deficiencies in the prior art, a kind of method that can indicate intuitively the electromyographic signal annoyance level is provided, to help the credibility of doctor's discernible signal.
Technical scheme: the indicating means that electromyographic signal of the present invention is disturbed comprises the steps:
(1), in time period t, obtain the sampled value sequence x (t) of electromyographic signal, t=1-10s;
(2) calculate this section electromyographic signal power spectrum P (f) and general power Tp;
(3) in the electromyographic signal bandwidth range, determine the performance number of 50Hz and frequency multiplication thereof (n*50) Hz: P50, P (n*50), n gets the natural number that is greater than 1;
(4) the performance number sum of calculating 50Hz and frequency multiplication (n*50) Hz accounts for the percentage ratio K of general power Tp, the size of K value representation annoyance level, and K=0-100, wherein 0 means without any interfering signal, 100 mean it is interfering signal entirely;
Mean the electromyographic signal quality with Q=100-K, Q=0-100, wherein, 0 means without any useful signal, Q=100 means it is useful signal entirely;
(5) K or Q-value are indicated to the user.
Preferably, in step (1), t=2s.
In formula, the fast Fourier transform that FFT (X) is x (t), the number that N is data;
Utilize indicating means of the present invention to carry out electromyographic signal and disturb the method for removing, comprise the steps:
(1) disturb indicated value judgement annoyance level according to signal: when annoyance level is low or when noiseless, without removing interfering signal, step finishes; When annoyance level is high, enter step (2);
(2) utilize digital filter to remove the signal of 50Hz frequency;
(3) again according to signal, disturb indicated value judgement annoyance level: when annoyance level is low or when noiseless, without removing interfering signal, step finishes; When annoyance level is high, enter step (4);
(4) repetitive cycling step (2), (3), utilize digital filter to remove the signal of (n*50) Hz frequency, and n gets the natural number that is greater than 1 successively, until, without removing interfering signal, step finishes.
Normally, when K=0-20 or Q=80-100, without removing interfering signal; Be that annoyance level is high when K>20 or Q<80, need to be removed interfering signal.Also can further accurately judge whether to need to remove to disturb according to special requirement.
The present invention compared with prior art, its beneficial effect is: the indicating means that 1, electromyographic signal disclosed by the invention is disturbed, can help the doctor judge the signal collected be by patient contract muscles cause or because interference causes, can indicate the order of severity of disturbing in electromyographic signal, the credibility of judgement signal, and then provide accurate foundation for doctor's further diagnosis; 2, the present invention disturbs indicating means, for removing to disturb, provides foundation, thereby can be in the situation that the least possible removal 50Hz and harmonic energy thereof obtain reliable signal; 3, due to 50Hz and harmonic wave place thereof, a small amount of useful signal is also arranged, the present invention is by disturbing indicating means to remove interfering signal, removal useful signal that can be as far as possible few, and the data that obtain approach real data more.
The accompanying drawing explanation
Fig. 1 is the surface myoelectric analytical system in the situation that disturb the very little parent mass peak peak signal collected, peak value 40uV left and right.
Fig. 2 is the surface myoelectric analytical system in the situation that disturb the very little rms signal collected, the 13uV left and right.
Fig. 3 is the surface myoelectric analytical system in the situation that disturb the very little signal spectrum figure collected.
Fig. 4 is the surface myoelectric analytical system in the situation that the primary signal figure disturbed without contracting muscles and gathering is arranged.
Fig. 5 is the surface myoelectric analytical system in the situation that the rms signal figure disturbed without contracting muscles and gathering is arranged.
Fig. 6 is the surface myoelectric analytical system in the situation that the signal spectrum figure disturbed without contracting muscles and gathering is arranged.
The specific embodiment
Below technical solution of the present invention is elaborated, but protection scope of the present invention is not limited to described embodiment.
Embodiment 1: followingly take the surface myoelectric analytical system as example, the existence of disturbing is described.This system has been used the 50Hz trap function, and the 50Hz power frequency is disturbed by filtering.Fig. 1 is this system in the situation that disturb the very little parent mass peak peak signal collected, peak value 40uV left and right, and Fig. 2 is rms signal, the 13uV left and right, Fig. 3 is spectrogram.In the situation that interference is very little, the surface myoelectric RMS value collected is 10uV left and right, signal of telecommunication when this 10uV is the real muscle contraction of patient.From spectrogram, can find out, the signal energy value in the 20-50Hz scope, has wherein been concentrated most energy in the 30-250Hz scope respectively, and this is typical surface electromyogram signal frequency spectrum.In noisy situation, patient does not contract muscles, but produces signal amplitude owing to disturbing, and as shown in Figure 4, Fig. 5 is rms signal 10uV left and right to the primary signal figure of collection, and Fig. 6 is signal spectrum.From spectrogram, can find out, signal capabilities mainly concentrates on the harmonic wave place of the 50Hz such as 100Hz, 200Hz, 300Hz, 400Hz, 150Hz, 250Hz, 350Hz, 450Hz, and the energy of 50Hz is by the filtering of 50Hz wave trap.
Therefore, the indicating means that electromyographic signal of the present invention is disturbed, comprise the steps:
(1) in time period t, obtain the sampled value sequence x (t) of electromyographic signal, t=1-10s, comparative optimization be 2s;
(2) calculate this section electromyographic signal power spectrum P (f) and general power Tp;
(3) in the electromyographic signal bandwidth range, determine the performance number of 50Hz and frequency multiplication thereof (n*50) Hz: P50, P (n*50), n gets the natural number that is greater than 1;
In formula, the fast Fourier transform that FFT (X) is x (t), the number that N is data;
(4) the performance number sum of calculating 50Hz and frequency multiplication (n*50) Hz accounts for the percentage ratio K of general power Tp, the size of K value representation annoyance level, and K=0-100, wherein 0 means without any interfering signal, 100 mean it is interfering signal entirely;
Mean the electromyographic signal quality with Q=100-K, Q=0-100, wherein, 0 means without any useful signal, Q=100 means it is useful signal entirely;
(5) K or Q-value are indicated to the user.
Utilize above-mentioned indicating means to carry out electromyographic signal and disturb the method for removing, comprise the steps:
(1) disturb indicated value judgement annoyance level according to signal: when annoyance level is low or when noiseless, without removing interfering signal, step finishes; When annoyance level is high, enter step (2);
(2) utilize digital filter to remove the signal of 50Hz frequency;
(3) again according to signal, disturb indicated value judgement annoyance level: when annoyance level is low or when noiseless, without removing interfering signal, step finishes; When annoyance level is high, enter step (4);
(4) repetitive cycling step (2), (3), utilize digital filter to remove the signal of (n*50) Hz frequency, and n gets the natural number that is greater than 1 successively, until, without removing interfering signal, step finishes.
Normally, when K=0-20 or Q=80-100, think noiseless or disturb low, without removing interfering signal; Be that annoyance level is high when K>20 or Q<80, need to be removed interfering signal.Also can further accurately judge whether to need to remove to disturb according to special requirement.
As mentioned above, although meaned and explained the present invention with reference to specific preferred embodiment, it shall not be construed as the restriction to the present invention self.Under the spirit and scope of the present invention prerequisite that does not break away from the claims definition, can make in the form and details various variations to it.
Claims (5)
1. the indicating means that electromyographic signal is disturbed, is characterized in that comprising the steps:
(1), in time period t, obtain the sampled value sequence x (t) of electromyographic signal, t=1-10s;
(2) calculate this section electromyographic signal power spectrum P (f) and general power Tp;
(3), in the electromyographic signal bandwidth range, determine the performance number of 50Hz and frequency multiplication thereof (n*50) Hz: P50, P (n*50), n gets the natural number that is greater than 1;
(4) the performance number sum of calculating 50Hz and frequency multiplication (n*50) Hz thereof accounts for the percentage ratio K of general power Tp, the size of K value representation annoyance level, and K=0-100, wherein 0 means without any interfering signal, 100 mean it is interfering signal entirely;
Mean the electromyographic signal quality with Q=100-K, Q=0-100, wherein, 0 means without any useful signal, Q=100 means it is useful signal entirely;
(5) K or Q-value are indicated to the user.
2. the indicating means that electromyographic signal according to claim 1 is disturbed is characterized in that:
In step (1), t=2s.
3. the indicating means that electromyographic signal according to claim 1 is disturbed, it is characterized in that: the computational methods of electromyographic signal power spectrum P (f) are:
In formula, the fast Fourier transform that FFT (X) is x (t), the number that N is data;
4. carry out electromyographic signal according to the described indicating means of claim 1-3 and disturb the method for removing, it is characterized in that comprising the steps:
(1) disturb indicated value judgement annoyance level according to signal: when annoyance level is low or when noiseless, without removing interfering signal, step finishes; When annoyance level is high, enter step (2);
(2) utilize digital filter to remove the signal of 50Hz frequency;
(3) again according to signal, disturb indicated value judgement annoyance level: when annoyance level is low or when noiseless, without removing interfering signal, step finishes; When annoyance level is high, enter step (4);
(4) repetitive cycling step (2), (3), utilize digital filter to remove the signal of (n*50) Hz frequency, and n gets the natural number that is greater than 1 successively, until, without removing interfering signal, step finishes.
5. electromyographic signal according to claim 4 is disturbed the method for removing, and it is characterized in that: when K=0-20 or Q=80-100, without removing interfering signal; Be that annoyance level is high when K>20 or Q<80, need to be removed interfering signal.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310486036.0A CN103494609B (en) | 2013-10-16 | 2013-10-16 | Method for indicating electromyographic signal interference and method for removing interference |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310486036.0A CN103494609B (en) | 2013-10-16 | 2013-10-16 | Method for indicating electromyographic signal interference and method for removing interference |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103494609A true CN103494609A (en) | 2014-01-08 |
CN103494609B CN103494609B (en) | 2015-05-06 |
Family
ID=49859942
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310486036.0A Active CN103494609B (en) | 2013-10-16 | 2013-10-16 | Method for indicating electromyographic signal interference and method for removing interference |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103494609B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016057823A1 (en) * | 2014-10-08 | 2016-04-14 | MAD Apparel, Inc. | Method and system for measuring beat parameters |
CN109901711A (en) * | 2019-01-29 | 2019-06-18 | 西安交通大学 | By the asynchronous real-time brain prosecutor method of the micro- expression EEG signals driving of weak Muscle artifacts |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002504405A (en) * | 1998-02-27 | 2002-02-12 | ユニベルシテ ドゥ モントリオール | Method and apparatus for activating ventilatory assistance in response to myoelectric activity |
US20050125473A1 (en) * | 2003-12-04 | 2005-06-09 | Lars Lindstrom | Signal filtering using orthogonal polynomials and removal of edge effects |
US20070038382A1 (en) * | 2005-08-09 | 2007-02-15 | Barry Keenan | Method and system for limiting interference in electroencephalographic signals |
CN101548885A (en) * | 2009-04-17 | 2009-10-07 | 南京大学 | Method for eliminating power frequency interfering signals in electrophysiological signals |
CN101773393A (en) * | 2010-01-14 | 2010-07-14 | 中国科学院合肥物质科学研究院 | Low-cost multi-channel surface myoelectric amplification module |
CN201578241U (en) * | 2010-01-18 | 2010-09-15 | 段峰 | Surface electromyogram signal filter |
-
2013
- 2013-10-16 CN CN201310486036.0A patent/CN103494609B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002504405A (en) * | 1998-02-27 | 2002-02-12 | ユニベルシテ ドゥ モントリオール | Method and apparatus for activating ventilatory assistance in response to myoelectric activity |
US20050125473A1 (en) * | 2003-12-04 | 2005-06-09 | Lars Lindstrom | Signal filtering using orthogonal polynomials and removal of edge effects |
US20070038382A1 (en) * | 2005-08-09 | 2007-02-15 | Barry Keenan | Method and system for limiting interference in electroencephalographic signals |
CN101548885A (en) * | 2009-04-17 | 2009-10-07 | 南京大学 | Method for eliminating power frequency interfering signals in electrophysiological signals |
CN101773393A (en) * | 2010-01-14 | 2010-07-14 | 中国科学院合肥物质科学研究院 | Low-cost multi-channel surface myoelectric amplification module |
CN201578241U (en) * | 2010-01-18 | 2010-09-15 | 段峰 | Surface electromyogram signal filter |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016057823A1 (en) * | 2014-10-08 | 2016-04-14 | MAD Apparel, Inc. | Method and system for measuring beat parameters |
CN109901711A (en) * | 2019-01-29 | 2019-06-18 | 西安交通大学 | By the asynchronous real-time brain prosecutor method of the micro- expression EEG signals driving of weak Muscle artifacts |
Also Published As
Publication number | Publication date |
---|---|
CN103494609B (en) | 2015-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102368951B (en) | Nerve monitoring during electrosurgery | |
Parák et al. | ECG signal processing and heart rate frequency detection methods | |
Li et al. | Denoising and R-peak detection of electrocardiogram signal based on EMD and improved approximate envelope | |
Nayak et al. | Filtering techniques for ECG signal processing | |
CN107041743A (en) | A kind of real-time R wave detecting methods of electrocardiosignal | |
CN108888259A (en) | A kind of real-time QRS wave detection method of electrocardiosignal | |
CN105054926A (en) | Electrocardiosignal feature information extraction method and device | |
CN105796096A (en) | Heart rate variability analysis method, heart rate variability analysis system and terminal | |
CN101449973A (en) | Judgment index generation method and device for cardiac interference signal identification | |
CN103417206B (en) | ECG removes the method and system of Hz noise | |
CN112754444A (en) | Radar-based non-contact pig respiration detection method | |
CN102429655B (en) | A kind of interference of the self adaptation myoelectricity based on electrocardiosignal minimizing technology | |
CN106974617A (en) | The Signal Pre-Processing Method and signal wave crest detection method of a kind of efficiently and accurately | |
CN105852850A (en) | Method and related device for extracting respiratory signals from ECG (electrocardiogram) signals | |
CN109009087B (en) | Rapid detection method for electrocardiosignal R wave | |
CN105266800A (en) | Fetal electrocardiogram blind separation method based on low signal-to-noise ratio | |
CN107997754A (en) | A kind of intelligence mattress system and human body physiological characteristics data extraction method | |
CN103494609B (en) | Method for indicating electromyographic signal interference and method for removing interference | |
CN112115856A (en) | Electroencephalogram quality evaluation method, storage medium and system | |
CN103169469B (en) | A kind of real-time brain electrical impedance detection method based on high frequency pumping | |
Xu et al. | Digital filter design for peak detection of surface EMG | |
CN112861328B (en) | Generator damping evaluation device and method based on random response signals | |
KR101048763B1 (en) | Apparauts and method for detecting signal | |
Dai et al. | Removal of baseline wander from dynamic electrocardiogram signals | |
CN111803054B (en) | Method and instrument for evaluating ventricular fibrillation signal quality and defibrillation success rate in real time |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20170418 Address after: 48 No. 210000 Jiangsu province Nanjing city Yuhuatai District software Avenue (SoHo International Plaza A building two B zone) Patentee after: Nanjing Weisi medical Polytron Technologies Inc Address before: 211100 Jiangning, Jiangsu Province, the District of Nanjing Road No. 18 (Science Park) Patentee before: Nanjing Mailande Medical Technology Co., Ltd. |
|
TR01 | Transfer of patent right |