CN1036175C - Automatic diagnostic instrument for ECG information - Google Patents
Automatic diagnostic instrument for ECG information Download PDFInfo
- Publication number
- CN1036175C CN1036175C CN94103588A CN94103588A CN1036175C CN 1036175 C CN1036175 C CN 1036175C CN 94103588 A CN94103588 A CN 94103588A CN 94103588 A CN94103588 A CN 94103588A CN 1036175 C CN1036175 C CN 1036175C
- Authority
- CN
- China
- Prior art keywords
- signal
- unit
- output
- mind
- frequency domain
- 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.)
- Expired - Fee Related
Links
- 238000012545 processing Methods 0.000 claims abstract description 77
- 238000004458 analytical method Methods 0.000 claims abstract description 35
- 230000006870 function Effects 0.000 claims description 27
- 238000001228 spectrum Methods 0.000 claims description 17
- 238000003786 synthesis reaction Methods 0.000 claims description 14
- 230000001427 coherent effect Effects 0.000 claims description 10
- 238000002955 isolation Methods 0.000 claims description 8
- 238000012546 transfer Methods 0.000 claims description 8
- 238000005314 correlation function Methods 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 5
- 230000001575 pathological effect Effects 0.000 claims description 5
- 230000005622 photoelectricity Effects 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 5
- 238000005311 autocorrelation function Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 238000005316 response function Methods 0.000 claims 2
- 230000000052 comparative effect Effects 0.000 claims 1
- 238000004148 unit process Methods 0.000 claims 1
- 238000003745 diagnosis Methods 0.000 abstract description 12
- 208000024172 Cardiovascular disease Diseases 0.000 abstract description 5
- 238000000034 method Methods 0.000 description 18
- 230000008569 process Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 7
- 208000029078 coronary artery disease Diseases 0.000 description 6
- 238000001514 detection method Methods 0.000 description 6
- 230000002159 abnormal effect Effects 0.000 description 5
- 210000004556 brain Anatomy 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 5
- 230000004069 differentiation Effects 0.000 description 5
- 201000010099 disease Diseases 0.000 description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 208000010125 myocardial infarction Diseases 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000001360 synchronised effect Effects 0.000 description 4
- 230000005856 abnormality Effects 0.000 description 3
- 206010003119 arrhythmia Diseases 0.000 description 3
- 230000006793 arrhythmia Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 210000000056 organ Anatomy 0.000 description 3
- 230000001225 therapeutic effect Effects 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000003759 clinical diagnosis Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 230000006855 networking Effects 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 230000035790 physiological processes and functions Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 206010006578 Bundle-Branch Block Diseases 0.000 description 1
- 208000031229 Cardiomyopathies Diseases 0.000 description 1
- 206010052895 Coronary artery insufficiency Diseases 0.000 description 1
- 206010020880 Hypertrophy Diseases 0.000 description 1
- 208000009525 Myocarditis Diseases 0.000 description 1
- 206010042434 Sudden death Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000028161 membrane depolarization Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000005405 multipole Effects 0.000 description 1
- 210000004165 myocardium Anatomy 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000000661 pacemaking effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 208000011580 syndromic disease Diseases 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The present invention relates to an automatic diagnostic instrument for omnidirectional stereo ECG signal, which comprises a plurality of detectors capable of synchronously detecting an ECG signal and a vectorcardiogram signal at different parts of the surface of a human body. The detected signals are collected by a signal collecting device. The collected signals are comprehensively analyzed from a time domain, a spatial domain, a frequency domain and a spatial domain frequency domain by a signal processing device. A processed result is output by a signal output device composed of a monitor and a drawing printer. The diagnostic instrument also comprises external storage equipment, a keyboard and an alarm device. The diagnostic instrument of the present ivnention can synchronously detect, collect and treat cardiovascular diseases at multiple domains and can carry out waveform analysis, parameter table output and automatic comprehensive diagnosis.
Description
The present invention relates to omnibearing stereo ecg information automatic diagnostic instrument.
As everyone knows, the cardiovascular diseases is the highest a kind of disease of middle-aged and elderly people M ﹠ M, and wherein, coronary heart disease, arrhythmia are again important topics highly visible in the cardiovascular diseases field.Current, mainly be the detection that depends on electrocardiosignal to cardiovascular disease diagnosis.Electrocardiosignal is a kind of quasi-periodic bioelectrical signals, in each cardiac electrical cycle, electrocardiosignal has all experienced from the electrocardio pace-making, pacing signal conduction and cardiac muscle be the strict process of finishing depolarization and multipole in an orderly manner under the regulation and control of conducted signal, form the complete cycle that electrocardiosignal changes thus, it has characteristic elements such as amplitude, frequency, time limit and phase place.At present, the detection to electrocardiosignal mainly is to realize by electrocardiogram and vectorcardiogram.
Electrocardiogram is on the special leads system, as wilson's lead2s, measure the time dependent function curve of electrocardiosignal, only can analyze the feature of electrocardiosignal in the time domain scope on the routine electrocardiogram, it has certain specificity to the diagnosis of coronary heart disease, myocardial infarction, cardiomyopathy and arrhythmia, but sensitivity is not high.Statistics shows: about 60% in the coronary heart disease is in the not significantly reaction experimentally of routine electrocardiogram application of load; Electrocardiographic diagnosis rate before the infraction sudden death also has only about 55%.
As replenishing and development of routine electrocardiogram, vectorcardiogram also is widely used clinically, the electrocardiosignal that vectorcardiogram is measured on X, Y, Z three lead system, form vectorcardiogram, and with the variation of ring-type pictorial representation on the size and Orientation that Descartes's rectangular coordinate is fastened, show in the cardiac electrical cycle each moment electrocardio-activity phase relation.Its graph direct, the dimensional orientation definite conception, careful accurately to time and spatial relationship expression.It is generally acknowledged: vectorcardiogram gives the diagnosis that swashs syndrome etc. to myocardial infarction, coronary insufficiency, bundle branch retardance, fascicular block, chamber hypertrophy and is better than electrocardiogram, once be difficult to provide Useful Information to diseases such as arrhythmia.
Although electrocardiogram and vectorcardiogram have certain effect aspect diagnosis of cardiovascular diseases,, as mentioned above, their effect all is subjected to restriction to a certain degree, but has so just formed clinical demand and the contradiction between the supplying apparatus in fact.How to solve this contradiction, the problem of pendulum in face of numerous medical personnels and scientific worker that Here it is.On the other hand, electronic technology in recent years, fast development of computer technology particularly, microcomputer being extensive use of in every field, can the enlightenment that give people be incorporated into them in the checkout gear of bioelectrical signals to remove (as electrocardio, brain), by these traditional signals that checkout gear detected are carried out analyzing and processing, obtain more signal or analysis processing result and be used for clinical again.In this respect, present inventor and Feng Genquan etc. are through years of researches, with regard to the collection of the electric physiological process of electrocardio, brain electricity and other vital organs, analyze and contribute, the document is disclosed on the Chinese patent CN87-102381A, and title is " biological signal detection processing unit and method " (on October 21st, 1988 is open).
CN87-102381A discloses a kind of relevant multiple signals of electric physiological process that utilize a plurality of electrodes synchronously to gather the heart, brain and other vital organs on the different parts of people or other animal body surface, in signal processing apparatus, each road signal is carried out the comprehensive analysis of time domain, frequency domain, spatial domain, amplitude domain, time difference domain dynamic process then, utilize each road correlation between signals and not the complementary action of the analytical method of same area carry out multifactorial comparison, confirmation and dynamic tracking, improve the accuracy of diagnosis and the reliability of differentiate judging with this.This device can also be exported the result who detects and handle by monitor and drawing printing equipment, so that for the clinical use of medical worker.This device comprises that also the comprehensive pathological index of multiple domain with result and the aforementioned various diseases that obtains through the clinical medicine statistics compares, can carry out automatic diagnosis to various diseases, can also start corresponding alert program and device according to diagnostic result to the critical patient, send to operator and report to the police and prompting.
The disclosed apparatus and method of CN87-102381A, the signal of telecommunication that has solved the heart, brain and other vital organ that can not synchronous acquisition preferably can detect in same device, by signal processing apparatus it is analyzed synchronously then, thereby the required data of clinical diagnosis, curve and report are provided.But the disclosed apparatus and method of CN87-102381A are that detection, processing and the analysis of synchronous electrocardiosignal seems then and can not adapt to that mainly there is the defective of two aspects in it being applied to originally.The first, it does not have enough measures to solve the shielding between the signal of being gathered, and therefore, it exists high frequency noise and 50HZ noise inevitably.The second, the signal of gathering is not utilized fully, it can only provide electrocardiogram (ECG) respectively, vectorcardiogram (VCG) and ecg frequency spectrum (FCG), so, also just limited its range of application clinically.
Therefore, the objective of the invention is to overcome the shortcoming in the prior art, a kind of omnibearing stereo ecg information automatic diagnostic instrument of novelty is provided, this diagnostic apparatus energy synchronous acquisition multichannel electrocardiosignal, and filter undesirable high frequency noise and 50HZ noise effectively, and the ground of prestage and the ground of microcomputer can be kept apart, thereby reduced interfering noise widely, so that guarantee that signal is undistorted; This diagnostic apparatus can not only provide ECG, VCG and the FCG in the prior art, and the analysis-by-synthesis figure (SCG) of mind-set amount spectrogram (FVG) and electrocardiogram and vectorcardiogram can also further be provided; This diagnostic apparatus can also be sent into the analysis result of above-mentioned five aspects specialist system and carry out analysis-by-synthesis and state of an illness classification, and then diagnostic result and therapeutic scheme are provided; This diagnostic apparatus can also start corresponding alert program and device to the critical patient, sends alarm to operator.
The objective of the invention is to solve by following proposal.Different parts on the body surface of people or other animal is placed the electrode that is no less than 16 lead systems, they should comprise the orthogonal mind-set amount X that 12 ECG signal and 3 of leading system are led at least, Y, the lead signals of Z, above-mentioned signal is synchronously by the signal pickup assembly collection, the signal that signal pickup assembly is gathered each road separates them with isolation level, in case the interference between stop signal, therefrom select then various differences lead or lead the combination signal by preposition amplification, filtering, the photoelectricity isolation level is sent into sampling holder, signal through above-mentioned harvester collection is sent into signal processing apparatus, (comprise frequency domain by a series of algorithm process, time domain, spatial domain, amplitude domain, the processing of time difference domain and time-space domain) and D/A conversion, draw electrocardiogram (ECG), ecg frequency spectrum (FCG), vectorcardiogram (VCG), mind-set amount spectrogram (FVG) and electrocardiogram and vectorcardiogram analysis-by-synthesis figure (SCG); Then, above-mentioned analysis result input specialist system is carried out analysis-by-synthesis and state of an illness classification, print diagnostic result and corresponding treatment scheme, simultaneously, can also start alert program and device, realize reporting to the police the critical patient.
Be not difficult to find out that from solution of the present invention it has done important improvement to the bioelectric signal measuring apparatus of prior art, make it be fit to detection, processing and the analysis of electrocardiosignal more, maximally utilised obtained information, provide clinical diagnosis required information, thereby improved the accuracy and the reliability of diagnosis greatly.
ECG, FCG, VCG, FVG and SCG are that electrocardiosignal is comprehensive, three-dimensional, the mutual integrated dynamic analysis of multiple domain (time domain, frequency domain, spatial domain, amplitude domain, time difference domain and time-space domain), thereby reach the advantage of giving full play to separately, replenish deficiency separately again, realize optimum organization, corroborate each other.It all has special diagnostic value for early coronary disease, myocarditis, the myocardial infarction that is not true to type, the conduction block of myocardial infarction merging bundle branch etc.
Diagnostic apparatus of the present invention can also not have the electrophysiological function that wound ground detects heart, and heart degree of aging with age.On the index basis of diagnosis, grade assessment, thus the normal person is carried out classification to the patients with coronary heart disease order of severity, the development to coronary heart disease has had a quantitative notion like this.
Purpose of the present invention and other purposes and feature and advantage are by will be more readily apparent from below in conjunction with the explanation of accompanying drawing to the preferred embodiments of the present invention, in the accompanying drawings:
Fig. 1 is the schematic block diagram according to omnibearing stereo ecg information automatic diagnostic instrument of the present invention;
Fig. 2 is an embodiment according to the signal pickup assembly of omnibearing stereo ecg information automatic diagnostic instrument of the present invention shown in Figure 1;
Fig. 3 is an embodiment according to the signal processing apparatus of omnibearing stereo ecg information automatic diagnostic instrument of the present invention shown in Figure 1;
Fig. 4 is the oscillogram of electrocardio time domain processing signals;
Fig. 6 is the oscillogram of mind-set quantity space territory signal processing;
Fig. 7 A, 7B, 7C are the oscillograms that mind-set amount frequency-region signal is handled;
Fig. 8 A, 8B, 8C, 8D, 8E, 8F, 8G, 8H are the oscillograms of complex electrocardio figure;
Fig. 9 is the drawing principle of description time vectorcardiogram;
Figure 10 is the workflow diagram according to omnibearing stereo ecg information automatic diagnostic instrument of the present invention;
Figure 11 is the detailed description flow chart of specialist system analysis-by-synthesis step 180 among Figure 10;
Figure 12 is to be the flow chart of the manual intervention program that provides of example with step 134 among Figure 10;
Figure 13 is to be the sketch map of used sign in the example explanation manual intervention program with the electrocardiogram;
Figure 14 is the detailed description flow chart of analysis scan step 187 in critical region among Figure 11;
Figure 15 is the schematic appearance of a ride formula configuration embodiment of omnibearing stereo ecg information automatic diagnostic instrument of the present invention.
Referring to Fig. 1, be depicted as the schematic structure block diagram of omnibearing stereo ecg information automatic diagnostic instrument of the present invention, number in the figure 11-In represents N detecting electrode, wherein, N can select arbitrarily according to clinical needs, but can not be less than 16, promptly Electrocardiographic 12 lead system and orthogonal vectorcardiogram X, Y, Z, three lead.Certainly, also can be as the disclosed electrode vest of CN1043621A.Label 2 is electrical signal collection devices, the preferred embodiment that its detailed structure provides referring to Fig. 2.Label 3 is signal processing apparatus, and its handles the signal of collecting from electrical signal collection device 2, by calculating and waveform recognition draws ECG, FCG, VCG, FVG and five kinds of curves of SCG and relevant parameters C, F, V, VF, TV.Print by signal output apparatus then, if above-mentioned curve and parameter have special abnormality, the alarm device of being represented by label 5 sends warning with sound equipment, light or their combination.Keyboard of label 6 expressions, the operator can import the various information relevant with testee by the operation of the whole diagnostic apparatus of Keyboard Control, for example, the name of testee, sex, age, medical record number, test time started or the like.It can also be intervened by keyboard 6 waveform recognition of signal processing apparatus 3.Obviously, keyboard 6 also can be substituted by other instruction inputting devices.External memory equipment of label 7 expressions, it can adopt the external equipment of any routine, be used to store treated various signal and the information relevant with tested person person, follow the tracks of so that carry out the accumulation of case history, the extension of data is handled and later medical statistics.
Referring to Fig. 2, the figure shows out the schematic block diagram of a preferred embodiment of the electrocardiogram signal acquisition device 2 among Fig. 1.The present invention is in order to take the needed signal of telecommunication from body surface, designed the system 1 of leading that is not less than 16 lines that lead, the signal of each path electrode isolation level by separately 22 respectively arrives resistance and leads and select networking 23, exportable electrocardiosignal and the orthogonal mind-set amount signal that contains conventional leads in this networking, these signals are sent to the soft switch 24 of electronic simulation, according to clinical needs, by the soft switch 24 of signal processing apparatus 3 control electronic simulations, variously lead or lead composite signal to select, these signals are sent to preamplifier 25 separately, preamplifier 25 is selected typical three amplifier difference amplifiers for use, it is amplified to the volt level without distortion with the signal of the microvolt level of body surface, so that post processing.Deliver to wave filter 26 through preamplifier 25 amplifying signals, wave filter is made up of low pass filter and 50HZ line trap, is used to filter high frequency noise and 50HZ noise.Filtered signal is again by one-level gain-adjusted level 27, to satisfy various needs.Then, signal is sent to photoelectricity isolation level 28, uses the effect of photoelectricity isolation level to be to make the ground of prestage and the ground of signal processing apparatus 3 to keep apart, thereby has reduced earth-return widely, has also lowered interfering noise widely.Simultaneously, it has more effectively guaranteed the safety of tested human body also with prestage and isolated from power.Signal is sent into sampling holder 29 after by the photoelectricity isolation level.So just finished signals collecting, so that use for subsequent treatment.
Referring to Fig. 3, it is an embodiment of the signal processing apparatus shown in Fig. 1.Signal processing apparatus 3 comprises a time domain processing unit 31, a frequency domain processing unit 32, a spatial domain processing unit 33, a spatial domain frequency domain processing unit 34, a complex electrocardio processing unit 35,36, one indexs in a waveform recognition unit compare to determine unit 37 and a control unit 38.
During signal processing apparatus 3 work, when above-mentioned 5 processing unit 31-35 worked by control unit 3 signal processing apparatus 3, above-mentioned 5 processing unit 31-35 passed through control unit 38 29 extraction data of the sampling holder in signal pickup assembly 2 respectively, wherein:
31 pairs of Electrocardiographic sample frequencys of time domain processing unit are 250HZ, and it is through digital filtering and output device 4 drawing are delivered in the compression processing respectively and external memory storage 7 carries out record.
Frequency domain processing unit 32 and spatial domain frequency domain processing unit 34 need extraction data respectively in order to carry out fourier transform, and every segment data is 2
nIndividual point, its frequency can be decided with the need, to ecg frequency spectrum (FCG) or mind-set amount spectrogram usually between 50HZ-500HZ, the two paths of signals that collects is carried out fast fourier transform as function X (t), Y (t) and Z (t) respectively, so calculate they power spectrum, self correlation, cross-correlation, transfer function, pulsion phase should, correlation function etc.Each function curve and mathematical derivation illustrate when key diagram 5A, 5B and Fig. 7 A, 7B and 7C below.Above-mentioned more detailed description is asked for an interview " principle and the application of heart and brain electrograph computer-aided analysis " (Science Press, publish in October, 1986) that Feng Genquan shows, and this book is introduced as a reference at this.The operation result of above-mentioned function is sent into output device 4 drawing respectively and input external memory 7 carries out record.
Also face is fixed with the need for the data fetch frequency of spatial domain processing unit 33, usually be between the 250HZ-1000HZ to vectorcardiogram, this unit carries out digital filtering to corresponding X, Y, Z three road signals, by intercepting this relevant paragraph of three tunnel, as the P in the vectorcardiogram, QRS, T ripple, constitute face amount (X, Y), plane of structure (X, Z), side (Y, Z) and corresponding three-dimensional stereo model are sent into output device 4 respectively and are drawn and send into external memory 7 and carry out record.
Complex electrocardio processing unit 35 is to 12 tunnel conventional electrocardiosignaies, and three lead the unit that mind-set amount signal X, Y, Z carry out comprehensive analysis processing.In order to carry out electrocardio, the multicycle of mind-set amount signal handles synchronously, this unit electricity of coring, and the minimum common frequency 500HZ of mind-set amount is a benchmark, automatically gets one section more than or equal to 5 electrocardios, mind-set amount cycle from the multiple signals of gathering.And it is carried out the identification of digital filtering and periodic characteristic, and for electrocardio and mind-set amount signal extraction periodic characteristic mark P, QRS, the T ripple, and indicated.Send into output device 4 drawing and external memory 7 at last and carry out record.
Waveform recognition unit 36 in the signal processing apparatus 3 comprises corresponding five parts with each processing unit 31-35, they carry out waveform recognition to the output of five processing unit 31-35 respectively, location Calculation, and can carry out manual intervention (this also will be described in detail in the back) as required, so that produce five relevant parameters tables, output device 4 is delivered in the output of waveform recognition unit 36, prints by the printing equipment in the output device 4.
The index of signal processing apparatus 3 compares to determine five parameter lists that 36 outputs of waveform recognition unit are accepted in unit 37, and with compare by the comprehensive pathological index of the multiple domain that sets in advance, and on this basis testing result is judged, and final report (FR) comprised that the final result of diagnosis and therapeutic scheme, this result send into output device 4 and print.In addition, can also send instruction to alarm device according to judged result reports to the police.
In addition, control unit 38 can offer this device and discerns and handle depositing relevant patient's data in the memorizer 7 externally in.Equally also can be with the input store as a result 7 of each processing unit.
The result that the waveform of time domain processing unit 31 output prints through output device 4 among the Fig. 3 of being shown in Figure 4, the curve shown in the figure is 12 routine electrocardiograms that lead.
Fig. 5 A, 5B and Fig. 7 A, 7B, 7C are respectively the oscillograms of electrocardio frequency spectrum processing unit 32 among Fig. 3 and 34 outputs of mind-set quantity space territory frequency domain processing unit.
Shown in Fig. 5 A, the 5B be with the curve V5 among Fig. 4 as function X (t), curve II by fourier transform, can be transformed into frequency-region signal as Function Y (t) with time-domain signal, the present invention at first utilizes formula in the frequency domain processing procedure:
To curve V5 shown in Figure 4 and curve II, promptly X (t) and Y (t) carry out fast fourier transform respectively, obtain frequency curve F (ω) and F (ω), then, and according to the computing formula of power spectrum:
Just can draw the power spectrum G of the curve V5 among Fig. 4 and the power spectrum G of curve II (shown in Fig. 5 A).
The computing formula of crosspower spectrum is:
Can draw curve V5 among Fig. 4 and the crosspower spectrum Gxx (seeing shown in Fig. 5 A) of curve II thus.
Utilize the computing formula of coherent function:
Thereby draw curve V5 shown in the curve RF among Fig. 5 A and the coherent function curve chart of curve II.
Pass through formula
Conversion, by transfer function H and the Q shown in the also exportable picture 5A of frequency domain processing unit 32, the 5B.The mould of HxyHxy (ω) wherein, Qxy is the argument of Hxy (ω), that is: Hxy=|Hxy (ω) |=| Gxy (ω) |/Gxx (ω) (7)
Wherein:
Curve PIH shown in Fig. 5 B is that curve V5 in 32 couples of Fig. 4 of frequency domain processing unit and curve II handle the two the pulse respond figure that exports the back, and this impulse response is the fourier transform of transfer function, utilizes formula
H (t)=F
-1[H (ω)] (10) carry out following conversion by means of frequency domain processing unit 32,
Get final product:
Rx(τ)=F
-1[Gxx(ω)] (12)
Ry(τ)=F
-1[Gyy(ω)] (13)
Rx (τ) and Ry (τ) are autocorrelator trace Vxx and the Vyy among Fig. 5 B, curve V5 in their difference presentation graphs 4 and the auto-correlation function curve chart of curve II.
The V of Fig. 5 B is curve V5 among Fig. 4 and the cross-correlation function curve chart of curve II, and it carries out the cross-correlation function computing by frequency processing unit 32 and draws, and its operational formula is as follows:
Hence one can see that,
Rxy(τ)=F
-1[Gxy(ω)] (15)
The Rxy here (t) just is equivalent to the Vxy among Fig. 5 B.
Shown in Figure 6 is the oscillogram after spatial domain is handled of spatial domain processing unit 33 outputs among Fig. 3, and that provide among the figure is X, Y, Z vectorcardiogram and face amount, plane of structure and lateral heart vector loop figure LS, H and F.
Fig. 7 A, 7B and 7C are the mind-set amount spectrograms (FVG) after the spatial domain frequency domain processing unit 34 in Fig. 3 is handled, the algorithm of each waveform of FVG is identical with FCG, difference be the fast fourier transform (FFT) of FCG be with the curve V5 shown in Fig. 4 and II as input quantity, VFG carries out computing with the variable that three of mind-set amount X, Y, Z shown in Figure 6 lead as FFT.
Power spectrum Gxx, Gyy, Gzz that X, Y, Z lead are:
The curve of Gxx, Gyy, Gzz is shown in Fig. 7-A.
Self correlation Vxx, Vyy, Vzz are:
Wherein Ux (t), Uy (t), Uz (t) are the mind-set amount signal that X, Y, Z lead, and t is a time difference.The curve of Vxx, Vyy, Vzz is shown in Fig. 7-B.
Cross-correlation Vxy is:
Wherein Ux (t), Uy (t) are X, the Y mind-set amount signal that leads, and τ is a time difference.The curve of Vxy is shown in Fig. 7-8.
Transfer function H xy (ω)
Hxy=|Hxy (ω) |=| Gxy (ω) |/Gxx (ω) wherein:
Wherein:
Hxy in the formula is the mould of Hxy (ω), and Qxy is the argument of Hxy (ω).The figure of Hxy, Qxy is shown in Fig. 7-C.
Impulse response PIH is:
h(t)=F
-1[Hxy(ω)]
The figure of P1H is shown in Fig. 7-C.
Coherent function RF is:
R(ω)=|Gxy(ω)|
2/Gxx(ω)·Gyy(ω)
The figure of coherent function is shown in Fig. 7-C.
Fig. 8 A, 8B, 8C, 8D, 8E, 8F, 8G, 8H are the oscillogram that the complex electrocardio processing unit 35 among Fig. 3 is exported by output device 4.As previously mentioned, complex electrocardio figure (SCG) is made up of 12 lead electrocardiogram and quadrature vectorcardiogram two parts.Wherein:
Fig. 8 A and 8B are 12 conventional lead electrocardiogram, i.e. I, II, III, aVR, aVL, aVF (above at Fig. 8 A), and V1-V6 (above at Fig. 8 B).
Fig. 8 C is an orthogonal ECG, and it is to be traced by electrocardiogram methods and obtained by the lead signal of X, Y, three faces of Z of input of Frank.
T-VCG among Fig. 8 D, 8E, 8F, 8G, the 8H is a timed vectorcardiogram, the graphy method of this figure be with Frank lead X, Y, the Z mind-set amount signal of input in plane coordinate system by certain direction, the successive projection that moves with certain speed.Its mathematical model is as follows:
With reference to Fig. 9, the coordinate that Pn is ordered is (Xn, Yn), and when the coordinate system OXY at Pn place along continuous straight runs moved with even velocity V, then at a time the position of some Pn (Xn, Yn) in former coordinate system of tn correspondence was:
Xn′=Xn+V(tn-t。)
The position of Yn '=Yn in former coordinate system is:
Xn′=Xn+V(tn-t
o)
Yn′=Yn
Be Pn ' (Xn+V Δ t, Yn).
So, by
Set constitute time heart vector loop, the formation of visible time heart vector loop
Being by parameter V Δ t having occurred in the abscissa value of a Pn, is a normal value if get V Δ t, then
Pn(Xn+n,Yn)
Can draw out timed vectorcardiogram T-VCG according to above-mentioned formula.
C-TVCG among Fig. 8 D is the transposition figure of timed vectorcardiogram T-VCG.
C-VCG among Fig. 8 D is continuous mind-set figure, and it is the mind-set amount circular chart that XY, YZ, three groups of signals of XZ are traced continuously by polar coordinate system.
Shown in Figure 10 is workflow diagram according to omnibearing stereo ecg information diagnostic apparatus of the present invention.
At first, start diagnostic apparatus, by the keyboard 6 inputs information relevant, as name, sex, age, detection time etc. with the person under inspection in step 110.Utilize a plurality of electrodes (being no less than 16) that the different parts of person under inspection's body surface is detected in step 120 then.To store through the digital information of over-sampling in step 130, use so that carry out time domain, spatial domain, frequency domain, spatial domain frequency domain and analysis-by-synthesis respectively for subsequent step 131,141,151,161 and 171.
Respectively through the data of 131,141,151,161 and 171 step process, produce the time domain deal with data of external memory record respectively in step 132,142,152,162 and 172, the spatial domain deal with data, electrocardio frequency domain deal with data, spatial domain frequency domain deal with data and complex electrocardio deal with data; And time domain deal with data, electrocardio frequency domain deal with data can be made electrocardio time-domain diagram (ECG), electrocardio frequency domain figure (FCG), and by the spatial domain frequency domain data draw mind-set amount frequency domain figure (FVG) this shown in step 133,153 and 163.Carry out waveform recognition in step 134,143,154,164 and 173 respectively through the data that handle in above-mentioned five territories.According to design of the present invention, above-mentioned steps can be carried out automatically, also can add manual intervention where necessary so that difficult waveform is carried out artificial aid identification, below seeing also for the detailed step of manual intervention aid identification to the description of Figure 12 and Figure 13.Produce C, V, F, VF and five tables of data of TV (Table I-Table V) thus.Wherein:
The C parameter list mainly contains the amplitude and the time width of each waveform of time domain figure, i.e. the amplitude and the parameter values such as time width and heart rate of the P ripple in the electrocardiogram, Q ripple, R ripple, S ripple and T ripple.
The V parameter list mainly comprises the orbiting motion direction of spatial domain figure, and angle and area parameters are as vectorcardiogram all quadrants area ratio, vector angle, initial sum termination locations, parameters such as the direction of rotation of vector loop.
The F parameter list mainly comprises the form and the position of electrocardio frequency domain figure, peak value and the correspondent frequency of the preceding four peak g1-g4 of power spectrum curve G shown in Fig. 5 A, 5B in the electrocardio frequency domain figure, the peak value at the main peak of impulse response PIH and negative sense peak and position, height and the position of the γ 1 of self correlation and cross-correlation curve V, V and V, γ 2 and γ 3, among the coherent function curve RF, with the corresponding coherent point f1 of frequency location of power spectrum curve first peak g1; Parameters such as the amplitude of transfer function amplitude curve Hxy peak h and frequency location.
The VF parameter list mainly comprises the form and the position of VFG figure, peak and the correspondent frequency value of power spectrum G, G in the mind-set amount spectrogram shown in Fig. 7 A, 7B, 7C, preceding four peak g1-g4 of G, the amplitude at the main peak of impulse response PIH and negative sense peak and position, height and the position of the γ 1 of self correlation and cross-correlation VXX, V, V, V, γ 2, γ 3, among the coherent function RF with parameters such as the amplitude of the coherent value f1 of the frequency location corresponding point of power spectrum first peak g1 and transfer function H xy peak h and frequency locations.
The position of TV parameter list essential record X, Y, the P of each phase of Z orthogonal ECG, QRS, T ripple, argument is used to instruct the output of complex electrocardio figure.
Carry out time-domain analysis, spatial domain analysis, electrocardio frequency-domain analysis, spatial domain frequency-domain analysis and analysis-by-synthesis respectively in step 135,145,155,165 and 175 respectively according to C, V, F, VF, five parameter lists of TV, and provide the analysis report of time domain, spatial domain, electrocardio frequency domain, spatial domain frequency domain and analysis-by-synthesis, shown in step 136,146,156,166 and 176 among the figure.And then carry out expert system analysis in step 180, will do more detailed discussion below to this.To carry out state of an illness classification in step 210 after expert system analysis, print diagnostic result and therapeutic scheme in step 300 afterwards, whole flow process is to this termination.
In addition, can draw vectorcardiogram according to the result of the spatial domain waveform recognition of step 143, shown in step 144.Equally, can draw the relevant figure of analysis-by-synthesis in the waveform recognition result of step 173 pair analysis-by-synthesis, this is shown in step 174.
Figure 11 is the detail flowchart of the various processes of specialist system analysis-by-synthesis step 180 representative among Figure 10.As shown in the figure, after having produced C, V, F, VF, TV parameter list, finish in step 181 prompting individual event analysis, thereby change step 182 over to, the comprehensive case of multiple domain is scanned, in step 183 scanning result and the given comprehensive pathological index of multiple domain are compared differentiation, if meet These parameters, can provide analysis report in step 184, show that scanning result is in the scope of the comprehensive pathological index of multiple domain.If but distinguish not within scope, then respectively the analysis result in each territory is checked in step 185, if when the check result of step 186 is not found positive value, then change step 187 over to and carry out the critical zone analysis scan, will go through in the back for the critical zone analysis scan.By the comprehensive distinguishing of critical zone index, if meet the critical zone discriminant criterion, then give, as shown in the step 189 by normal report.If belong to abnormal conditions, then provide exception reporting or prompting explanation at step 190 place through differentiation.When finding that it is positive that two territories or the analysis result more than two territories are arranged, shown in step 192, then at step 193 output analysis report.If not so, then continue execution in step 195, determine whether and can provide definite analysis result according to the positive index of existing single domain, if can, then provide analysis result in step 196.Otherwise execution in step 198, whether in other several territories, exist corresponding relation so that determine the positive findings of this single domain analysis, if there is this corresponding relation, then provide definite analysis report, otherwise will judge whether this territory positive findings is false positive, and result of determination is indicated as true positives in step 201 in step 199, then provide the prompting explanation in step 202, otherwise will think that this positive findings is a false positive, and be vetoed, shown in step 204.With in analysis result among the figure or the explanation input state of an illness classification 210, carry out subsequent treatment at last.
Referring to Figure 12, be the detailed process of step 134 time-domain analysis among Figure 10 shown in this figure, should be noted that step 143,173 analytic process are identical with step 134, describe step 134 at this.In the identification of step 1341 beginning time domain waveform, step 1342 need to determine whether manual intervention, should need not manual intervention for the waveform of routine, can directly enter step 1346 and discern automatically, produces the C parameter list.For difficult waveform, enter waveform or markings that step 1343 demonstrates needs manual intervention on monitor, as shown in figure 13, provide markings code name (as the Pb line among Figure 13) and displacement thereof by the operator in step 1344, make it reach the desired position.Determine whether to finish manual intervention in step 1345,, then turn back to step 1344, move another markings, after finishing manual intervention, enter step 1346, to the waveform generation C parameter after the identification if do not finish.Should be understood that, because multiple signals are synchronous acquisition, therefore, only one road signal is wherein carried out artificial cognition, can correspondingly determine the waveform position of other each road signal, so a certain road signal that only needs the operator to select shows separately on monitor, and, can greatly improve accuracy in this way to difficult waveform recognition by with the local amplification of waveform.
Figure 13 shows that electrocardiographic wave passes through the initial sum final position of markings Pb, Pe, Qb, Se, the determined P ripple of Tx, Te, QRS ripple and T ripple.Dotted line Pb ' has represented the displacement process of markings Pb.
Figure 14 shows that the detailed description of the critical zone analysis scan step 187 among Figure 11.When the differentiation of the single domain analysis result of step 186 is all negative, program enters step 187, and the related parameter that respectively has in C, F, V, VF, five parameter lists of TV is scanned, and differentiates that parameter and is within the critical zone that sets in advance.Should be understood that, the differentiation result in each territory can divide positive, negative and a critical zone between the two in the device of the present invention, this critical zone is to determine according to a large amount of clinical case statistical results and expert's experience, mainly be meant one section span between positive index and negative index, in this scope, be difficult to normal condition is distinguished mutually with the mile abnormality situation.For example, normal index is below the 2.40mv, and abnormal index is more than the 2.50mv, then the critical zone between the 2.41-2.49.In step 182, if there is not the interior parameter of critical range in above-mentioned five parameters, then in the normal report of step 189 output.If parameter of critical region (step 1872) is arranged, then carry out Fuzzy Processing in step 1872, its method is that the parameter of critical zone is carried out the Euclidean distance computing, and compare differentiation in step 1873 pair The above results, if its spatial distribution point is not in predetermined abnormal area, then turn back to the normal report of step 189 output through step 1874, otherwise in report of step 190 output abnormality or relevant prompting.
Figure 15 shows that the schematic appearance of ride formula configuration of a specific embodiment of diagnostic apparatus of the present invention, wherein label 403 is a main frame, it comprises the signal of telecommunication gathering-device 2 of omnibearing stereo ecg information diagnostic apparatus of the present invention, signal processing apparatus 3, alarm device 5, keyboard 6, external memory 7.404 is monitor, and it constitutes an integral body with external monitor 7 and main frame, and 405 is high speed temperature-sensitive drawing printer, and 406 is ride.Owing to adopt the ride formula to join, be convenient to the use of in hospital, flowing, and also be suitable for the needs of promptly paying a home visit.
The above is a specific embodiments of the present invention, can make various improvement and conversion on this basis for those skilled in the art, but these are all within protection scope of the present invention.
Table I
I II III aVR aVL aVF V1 V2 V3 V4 V5 V6
P1a .05 .09 .02 -.07 .00 .06 .02 .05 .04 .04 .05 .05
P2a .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00
p1d .09 .11 .08 .10 .03 .10 .08 .08 .08 .08 .09 .08
p2d .00 .00 .00 .00 .05 .00 .00 .00 .00 .00 .00 .00
Qa .00 .00 -.20 -.56 .00 -.00 .00 .00 .00 .00 .00 -.03
Qd .00 .00 .04 .06 .00 .00 .00 .00 .00 .00 .00 .02
Ra .68 .52 .76 .03 .43 .60 .17 .68 .85 1.68 1.77 1.11
Rd .05 .11 .08 .04 .05 .10 .04 .04 .04 .05 .04 .05
Sa -.43 .00 .00 .00 -.56 .00 -.72 -2.23 -2.22 -1.95 -.78 -.07
Sd .07 .00 .00 .00 .07 .00 .08 .08 .08 .07 .06 .05
R″ .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00
STa .00 .00 -.06 .00 .00 .00 .14 .42 .33 .11 .00 -.07
STd .00 .00 .00 .00 .00 .00 1.00 1.00 1.00 1.00 .00 .00
T1a .20 .10 -.13 -.14 .15 .04 .00 .91 .71 .33 .20 .14
T2a .00 .00 .00 .00 .00 .00 -.02 .00 .00 .00 .00 .00
P-R .15 .16 .15 .16 .15 .16 .15 .15 .15 .15 .15 .15
Q-T .37 .35 .36 .36 .36 .35 .36 .36 .37 .37 .37 .36
QRS .12 .11 .12 .10 .12 .11 .12 .12 .12 .12 .11 .12
PP min: PP max: PP ave: A.HR:
RR min: RR max: RR ave: V.HR:85.0
QRS: .12 PR Int: .16 QT Int: .36 QTC:
AXIS=74.00CLOCKWISE ROTATION
ECG FINDING:
NORMAL SINUS RHYTHM
ST ELEVATION:V2,V3,
SIGNATURE: Table II
***** FCG DIAGNOSIS *****GXX Mark 1/2 HG 1-N
+ - -
3/4N TU 5/10
5 - - +GYY Mark 1/2 HG 1-N
+ - -
3/4N TU 5/10
5 - - +QXY Mark D W D+W
5 + - -PIH Mark PV M1 M2 M3
2 - + - -VXX.VYY Mark RH RL FPX
FPY TVX TVY
0 - - -VXY Mark RV RD RF NW
5 - - - +RF Mark Cp Ct Cb1
0 - - -TOTAL: 8LAST SUCGESTION: ABNORMAL FCG: +
***** FCG PARAMETER TABLE *****
Rxy Rx Ry Pmax Pmin
7.75 4.53 19.63 213.76-189.39Gxx 1 2 3 4A 2.5 5.0 6.9 3.6F 5.6 11.2 16.8 21.6Gyy 1 2 3 4A 1.9 17.1 5.7 4.1F 5.6 11.2 16.8 22.4
Table III.
Table IV
| F LS H | |
| P | rot direction CCW CCW CCW max.v-amp(mv) .00 .00 .00 |
| Q R S | rot direction CW CCW CCW 0.01″ amp (mv) .14 .28 .31 0.01″ ang (degree) 105.9 174.60 63.19 0.02″ amp (mv) .69 .38 .75 0.02″ ang (degree) 14.94 151.80 26.46 0.03″ amp (mv) 1.10 .38 1.05 0.03″ ang (degree) 18.88 68.73 -7.58 0.04″ amp (mv) .40 .83 .73 0.04″ ang (degree) 89.42 28.48 -89.69 max v-amp(mv) 1.14 1.05 1.19 max ang(degress) 16.44 17.33 -122.22 duration(sec) 1.23 area L/R .43 .45 .30 ″ S/I .00 .00 2.19 ″ P/A Total Area: 368.00 267.00 927.00 area I% 70.10 68.90 31.10 ″ II% 29.80 31.00 .20 ″ III% .00 .00 22.60 ″ IV% .00 .00 45.90 |
| T | rot direction CW CCW CCW max v amp(mv) .26 .31 .36 ang(degree) 29.16 169.1 58.32 |
| QRS-T angle(degree) 12.12 152.58 180.53 QRS/T ratio 4.43 3.35 3.33 | |
| VCG FINDING: | |
GXX Mark 1/2 HG 1-N 3/4N TU
+ - - - -
O U E P
5 + - - +
GYY Mark 1/2 HG 1-N 3/4N TU
+ - - - -
O U E P
5 - - - -
GZZ Mark 1/2 HG 1-N 3/4N TU
+ - - - -
O U E P
5 + - - -
QXY Mark D W DW T S
5 + - - - -
PIH Mark PV M1 M2 M3 MZ
5 - + + - - v/x/y/z Mark RH RLX RLY RLZ FPX FPY
- - - - - +
FPZ TVX TVY TVZ
5 - - - +
VXY Mark RV RD NW+ NW- VP
5 - + - + -
RF Mark Cp Ct V RW
2 - - - +
TOTAL:10 LAST SUGGESTION:SERIOUS ABNORMAL VFG:++ Table V
I II III aVR aVL aVF V1 V2 V3 V4 V5 V6 X Y Z PB1 .05 .05 .06 .05 .051 .052 .05 .051 .051 .051 .051 .052 .05 .05 .05 PE1 .02 .02 .02 .02 .019 .02 .018 .02 .02 .02 .02 .02 .02 .02 .02 QB1 .14 .13 .14 .13 .14 .14 .15 .14 .14 .15 .14 .15 .14 .14 .15 QE1 .08 .08 .07 .08 .08 .08 .08 .08 .08 .07 .08 .08 .08 .07 .08 SB1 .32 .31 .32 .33 .32 .33 .32 .32 .32 .33 .31 .32 .31 .31 .31 SE1 .11 .10 .10 .11 .12 .11 .11 .11 .10 .11 .10 .12 .10 .11 .12 PB2 .45 .44 .44 .45 .41 .45 .44 .45 .45 .43 .45 .45 .43 .45 .44 PE2 .02 .03 .02 .02 .04 .03 .03 .02 .03 .02 .02 .02 .02 .03 .03 QB2 .51 .50 .51 .50 .53 .50 .51 .51 .52 .52 .50 .50 .51 .51 .52 QE2 .09 .10 .09 .09 .09 .10 .10 .09 .08 .09 .10 .09 .10 .09 .09 SB2 .67 .66 .65 .66 .65 .66 .65 .66 .67 .67 .67 .66 .65 .67 .68 SE2 .12 .11 .12 .12 .13 .11 .12 .12 .12 .13 .12 .11 .12 .13 .12 PB3 .85 .86 .87 .85 .84 .85 .85 .85 .84 .83 .84 .84 .84 .85 .84 PE3 .03 .03 .04 .04 .03 .03 .02 .03 .03 .04 .03 .04 .03 .04 .03 QB3 1.0 1.0 1.1 1.0 1.0 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.0 1.1 1.0 QE3 .09 .08 .09 .09 .08 .09 .09 .09 .09 .08 1.0 .09 0.9 0.8 0.8 SB3 1.4 1.4 1.3 1.4 1.4 1.3 1.3 1.4 1.4 1.4 1.3 1.4 1.3 1.4 1.3 SE3 .11 .10 .10 .10 .10 .10 .11 .10 .10 .10 .11 .11 .11 .12 .11 pB4 1.6 1.5 1.6 1.6 1.6 1.5 1.5 1.4 1.6 1.5 1.6 1.5 1.5 1.5 1.6 PE4 .03 .03 .03 .04 .04 .04 .04 0.3 .04 .03 .04 .03 .04 .03 .03 QB4 1.7 1.7 1.6 1.7 1.6 1.5 1.6 1.6 1.5 1.7 1.6 1.6 1.5 1.6 1.6 QE4 .12 .13 .12 .12 .13 .13 .12 .13 .13 .14 .12 .12 .12 .12 .13 SB4 1.9 1.9 2.0 1.8 1.9 1.9 1.9 1.8 1.9 1.9 1.9 1.8 2.0 1.9 1.8 SE4 .13 .12 .13 .12 .12 .13 .13 .13 .12 .12 .13 .13 .12 .12 .13 PB5 2.2 2.2 2.3 2.2 2.2 2.2 2.3 2.2 2.3 2.3 2.3 2.3 2.2 2.2 2.2 PE5 .03 .03 .03 .04 .04 .04 .04 0.3 .04 .03 .04 .03 .04 .03 .03 QB5 2.5 2.6 2.6 2.4 2.6 2.5 2.6 2.6 2.5 2.4 2.6 2.6 2.5 2.6 2.6 QE5 .12 .13 .12 .12 .13 .13 .12 .13 .13 .14 .12 .12 .12 .12 .13 SB5 2.8 2.9 2.8 2.8 2.9 2.9 2.9 2.8 2.9 2.9 2.9 2.8 2.8 1.8 2.8 SE5 .13 .12 .12 .12 .12 .13 .13 .13 .12 .12 .13 .13 .13 .12 .13
Claims (12)
1. omnibearing stereo ecg information automatic diagnostic instrument, it is by a plurality of electrodes, and electrical signal collection device and signal processing apparatus and signal output apparatus are formed, and it is characterized in that:
Said a plurality of electrode is no less than 16, comprising 12 electrocardioelectrodes, and three mind-set amount electrodes;
Said electrical signal collection device is synchronously to measure multichannel electrocardio and mind-set amount signal by above-mentioned a plurality of detecting electrodes from the different parts of body surface;
Said signal processing apparatus is that the multiple signals of above-said current signal harvester output are handled, and comprising:
An electrocardio time domain processing unit, this unit carry out time domain to each road electrocardiosignal of described signal pickup assembly output to be handled;
An electrocardio frequency domain processing unit, this unit carry out the electrocardio frequency domain to each road electrocardiosignal of described signal pickup assembly output to be handled;
A mind-set quantity space territory processing unit, this unit carries out the spatial domain processing to each road mind-set amount signal of described electrical signal collection device output;
A mind-set quantity space territory frequency domain processing unit, this unit carries out the processing of spatial domain frequency domain to the two-way of appointing in the mind-set amount signal of described electrical signal collection device output;
An integerated analytic unit, this unit is handled the electrocardio and the mind-set amount signal of described electrical signal collection device output;
Said signal output apparatus is used to export the result of said signal processing device, comprising:
The time-domain curve of each road electrocardiosignal of electrocardio time domain processing unit;
Electrocardio frequency domain processing unit is handled the frequency domain of appointing the two-way electrocardiosignal, comprises the realization fast fourier transform, the power spectrum function of calculating, auto-correlation function, cross-correlation function, transfer function, impulse response function, the respective function curve of coherent function;
The X-Y face amount of the mind-set amount signal of mind-set quantity space territory processing unit, the plane of structure that X-Z forms, the side that Y-Z forms and the corresponding geometric locus of three-dimensional stereo model;
Mind-set quantity space territory frequency domain processing unit is handled appointing two-way mind-set amount signal to carry out the spatial domain frequency domain, comprises the realization fast fourier transform, calculates their power spectrum function, auto-correlation function, cross-correlation function, transfer function, impulse response function, the corresponding function curve of coherent function;
The electrocardio of integerated analytic unit and mind-set amount analysis-by-synthesis figure.
2, according to the described omnibearing stereo ecg information of claim 1 automatic diagnostic instrument, it is characterized in that described electrical signal collection device carries out isolated isolation level by the signal that each road is gathered, the resistance selection network that leads, the soft switch of electronic simulation, preamplifier state, filtering stage, gain-adjusted, photoelectricity isolation level and sampling holder are formed.
3,, it is characterized in that described preamplifier is three amplifier difference amplifiers according to the described omnibearing stereo ecg information of claim 2 automatic diagnostic instrument.
4,, it is characterized in that described wave filter is made up of low pass filter and 50HZ line trap according to the described omnibearing stereo ecg information of claim 2 automatic diagnostic instrument.
5, according to the described automatic diagnostic instrument for ECG information of claim 1, after it is characterized in that integerated analytic unit is analyzed, the curve of drawing by output device comprises twelve-lead electrocardiogram, orthogonal ECG, timed vectorcardiogram and continuous vectorcardiogram.
6,, it is characterized in that described signal processing apparatus further comprises according to the described automatic diagnostic instrument for ECG information of claim 1:
An electrocardio time domain waveform recognition unit, relevant parameters is discerned and obtained in this unit to the wave amplitude of each waveform of the curve of described time domain processing unit output and ripple are wide;
An electrocardio frequency domain recognition unit, it discerns and obtains relevant parameters to the form and the position of described electrocardio frequency domain processing unit curve of output;
A mind-set quantity space territory recognition unit, it discerns and obtains relevant parameters to area under the curve, angle and the track trend of the output of described mind-set quantity space territory processing unit;
A mind-set quantity space territory frequency domain recognition unit, it discerns and obtains relevant parameters to the form and the position of the curve of described mind-set quantity space territory frequency domain unit output;
An analysis-by-synthesis recognition unit, it discerns and obtains relevant parameters to the curve of described comprehensive analysis processing unit output;
Through above-mentioned five parameter lists that recognition unit processes produced, through the output device output print.
7,, it is characterized in that described signal processing apparatus also comprises according to the described automatic diagnostic instrument for ECG information of claim 6:
A specialist system integerated analytic unit, comparing the parameter of described five waveform recognition unit output in this unit with the multiple domain pathological index value that presets, judges according to comparative result, and differentiate the result by output device output.
8,, it is characterized in that described signal processing apparatus also comprises according to the described automatic diagnostic instrument for ECG information of claim 7:
A state of an illness stage unit, this unit carries out state of an illness classification to specialist system integerated analytic unit discrimination result;
Described automatic diagnostic instrument for ECG information also comprises:
A warning system, this system starts warning system according to the state of an illness classification results of state of an illness stage unit, so that send seriously ill alarm signal to the operator.
9, according to the described automatic diagnostic instrument for ECG information of claim 1, it is characterized in that:
Described signal processing apparatus can be made up of a computer, is used for the signal of signal pickup assembly is sampled, and carries out the processing of time domain, electrocardio frequency domain, spatial domain, spatial domain frequency domain and analysis-by-synthesis, thereby produce the output result in described territory.
10, according to the described automatic diagnostic instrument for ECG information of claim 9, it is characterized in that: this diagnostic apparatus also comprises a keyboard, can supply the operator to this diagnostic apparatus input operation instruction, data, and waveform recognition is intervened.
11, according to the described automatic diagnostic instrument for ECG information of claim 9, it is characterized in that also comprising an external memory, be used to store each experimenter's diagnostic message.
12,, it is characterized in that described signal output apparatus comprises a monitor and a drawing printer according to the described automatic diagnostic instrument for ECG information of claim 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN94103588A CN1036175C (en) | 1994-04-14 | 1994-04-14 | Automatic diagnostic instrument for ECG information |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN94103588A CN1036175C (en) | 1994-04-14 | 1994-04-14 | Automatic diagnostic instrument for ECG information |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1099254A CN1099254A (en) | 1995-03-01 |
| CN1036175C true CN1036175C (en) | 1997-10-22 |
Family
ID=5031153
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN94103588A Expired - Fee Related CN1036175C (en) | 1994-04-14 | 1994-04-14 | Automatic diagnostic instrument for ECG information |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN1036175C (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100488448C (en) * | 2007-01-10 | 2009-05-20 | 赵峰 | Three-dimensional image electrocardiograph and implementation method |
| EP2429643B1 (en) * | 2009-04-30 | 2017-08-23 | Medtronic, Inc. | Patient state detection based on support vector machine based algorithm |
| CN103876729B (en) * | 2014-03-28 | 2016-05-18 | 浙江大学 | A kind of cardioelectric monitor system based on photoelectric sensing |
| CN116196013B (en) * | 2023-04-25 | 2023-08-15 | 毕胜普生物科技有限公司 | Electrocardiogram data processing method, device, computer equipment and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN87102381A (en) * | 1987-03-30 | 1988-10-12 | 创建基立有限公司 | Bioelectric signal detection and processing device and method |
| CN1041692A (en) * | 1988-10-07 | 1990-05-02 | 赵峰 | Multifunctional automatic diagnosis and treatment instrument for complex electrocardio-activity |
-
1994
- 1994-04-14 CN CN94103588A patent/CN1036175C/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN87102381A (en) * | 1987-03-30 | 1988-10-12 | 创建基立有限公司 | Bioelectric signal detection and processing device and method |
| CN1041692A (en) * | 1988-10-07 | 1990-05-02 | 赵峰 | Multifunctional automatic diagnosis and treatment instrument for complex electrocardio-activity |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1099254A (en) | 1995-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN87102381A (en) | Bioelectric signal detection and processing device and method | |
| CN107358196B (en) | Heart beat type classification method and device and electrocardiograph | |
| US20230329948A1 (en) | System, method, and computer algorithm for characterization and classification of electrophysiological evoked potentials | |
| US6607480B1 (en) | Evaluation system for obtaining diagnostic information from the signals and data of medical sensor systems | |
| Wartzek et al. | ECG on the road: Robust and unobtrusive estimation of heart rate | |
| CA2784267C (en) | Method and device for point-of-care neuro-assessment and treatment guidance | |
| CN1301080C (en) | Reduced lead set for detecting symptoms of acute cardiac ischemia | |
| US8948860B2 (en) | Field-deployable concussion detector | |
| JP4386235B2 (en) | Method and apparatus for sequential comparison of electrocardiograms | |
| US20070135727A1 (en) | Detection of artifacts in bioelectric signals | |
| EP0489010A1 (en) | Cardiac biopotential analysis system and method | |
| CN109288515B (en) | Method and device for periodic monitoring of premature beats in wearable ECG signals | |
| CN112057087B (en) | Method and device for data processing of autonomic nervous function in high-risk population of schizophrenia | |
| CN112426162A (en) | Fatigue detection method based on electroencephalogram signal rhythm entropy | |
| JP5422513B2 (en) | Derived electrocardiogram generation system and derived electrocardiogram generation method | |
| KR102802002B1 (en) | Ecg measurement service providing method and system using potable ecg device | |
| JP2007514488A (en) | System and method for analyzing electrocardiogram curvature and drug effects in long QT syndrome | |
| CN1032288A (en) | Method and device for multi-domain processing and correlation analysis of electrocardiographic signals | |
| CN106859673A (en) | A kind of depression Risk Screening system based on sleep cerebral electricity | |
| CN1036175C (en) | Automatic diagnostic instrument for ECG information | |
| US20130035604A1 (en) | Frequency Analysis of 12-Lead Cardiac Electrical Signals to Detect and Identify Cardiac Abnormalities | |
| CN105144225B (en) | Multi-patient EEG monitoring | |
| CN117297621A (en) | An AI intelligent analysis method and system for ECG data | |
| CN1265869A (en) | Heart disease testing and positioning system | |
| EP4609784A1 (en) | Apparatus and method for electrocardiogram ("ecg") signal analysis and artifact detection |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C06 | Publication | ||
| PB01 | Publication | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C19 | Lapse of patent right due to non-payment of the annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |