System and method for analyzing selected signal components in electrocardiographic signals, particularly late potentials in electrocardiograms

• US 5,211,179 A
• Filed: 07/13/1990
• Issued: 05/18/1993
• Est. Priority Date: 07/14/1989
• Status: Expired due to Fees

• Current Assignee: Hans F. Schels, Ralph Haberl
• Original Assignee: Hans F. Schels, Ralph Haberl
• Inventors: Schels, Hans F., Haberl, Ralph
• Primary Examiner(s): Howell, Kyle L.
• Assistant Examiner(s): SCHAETZLE, KENNEDY
• Application Number:
  US07/553,500
• Time in Patent Office: 1,040 Days
• Field of Search: 128/696, 128/702, 128/703, 128/704, 364/413.06, 364/413.02, 364/413.03, 364/413.05
• US Class Current: 600/515
• CPC Class Codes: A61B 5/349 Detecting specific parameters of the electrocardiograph cycle
• CPC Class Codes: A61B 5/413 Monitoring transplanted tissue or organ, e.g. for possible rejection reactions after a transplant
• CPC Class Codes: A61B 5/7257 using Fourier transforms

Abstract

Disclosed is a system and a method for computer-supported analysis of arrhythmic potentials in ECG signals, particularly those of late potentials. Interference discrimination advantages of frequency domain analysis are combined with temporal localization advantages of time domain analysis to determine the accurate location of arrhythmic potentials. Several small signal segments are selected in an ECG waveform. A determination is made of parameters corresponding to extended signals which closely match fluctuations of each respective small signal, allowing more information is discerned about the small signals than is possible with more conventional techniques. A comparison is made with respect to extended signals rather than small signals. Two autoregressive models are used, the maximum entropy method and adaptive filter determination. Area integrals of the frequency characteristics of small signal segments are recorded successively with respect to the frequency range of the arrhythmic potentials analyzed. The area integrals for first segments are calculated, then the remaining area integrals are calculated to determine the spectral manifestation of the selected signal region. A normality factor is determined as quotient of the summed area integrals of the first and second number of time segments analyzed.

Claims

1. A system for analyzing selected signal components in physiological measurement signals, comprising:

means for selecting a portion of a physiological measurement signal to be analyzed;

means for selecting a plurality of substantially overlapping mutually offset time segments of said portion;

means for deriving a plurality of low order frequency functions representing the low-band spectral components in each respective time segment;

means for deriving a plurality of high order frequency functions representing the broadband spectral components in each respective time segment;

means for deriving a plurality of difference frequency functions by subtracting said low order frequency functions from respective high order frequency functions; and

means for comparing said difference frequency functions, the comparison indicating detailed analysis of selected signal components in physiological measurement signals.

2. A system in accordance with claim 1 wherein said means for deriving low order and high order frequency functions each comprise means for transforming time domain data values into frequency domain data values.

3. A system in accordance with claim 2 wherein said means for transforming time domain data values into frequency domain data values comprises means for autoregressive computation.

4. A system in accordance with claim 1 wherein said means for comparing comprises:

means for calculating a plurality of area integrals over a particular frequency range, one for each respective difference frequency function; and

means for comparing said area integrals, so as to indicate parametric changes associated with said time segments.

5. A system in accordance with claim 1 wherein said portion of a physiological measurement signal to be analyzed is a portion of an electrocardiographic signal having a QRS complex obtained from a patient.

6. A system in accordance with claim 5 wherein said means for selecting a portion of an electrocardiographic signal comprises:

means for determining a velocity vector function for the electrocardiographic signal;

means for determining the rate of change of said velocity vector function, said rate of change indicating generally the beginning and end of the QRS complex; and

means for determining an autocorrelation peak of said velocity vector function in a particular time interval, said autocorrelation peak representing the end of the QRS complex.

7. A system in accordance with claim 6 wherein said means for selecting a plurality of mutually offset time segments comprises means for determining where said plurality of mutually offset time segments are to be positioned relative to the location of said autocorrelation peak.

8. A method for analyzing selected signal components in physiological measurement signals, comprising the steps of:

selecting a portion of a physiological measurement signal to be analyzed;

selecting a plurality of substantially overlapping mutually offset time segments of said portion;

deriving a plurality of low order frequency functions representing the low-band spectral components in each respective time segment;

deriving a plurality of high order frequency functions representing the broadband spectral components in each respective time segment;

deriving a plurality of difference frequency functions by subtracting said low order frequency functions from respective high order frequency functions; and

comparing said difference frequency functions, the comparison indicating detailed analysis of selected signal components in physiological measurement signals.

9. A method in accordance with claim 8 wherein said portion of a physiological measurement signal to be analyzed is a portion of an electrocardiographic signal having a P wave, a QRS complex, and an ST portion obtained from a patient.

10. A method in accordance with claim 9 wherein said time segments are at the trailing edge of the QRS complex and the ST portion for detecting late potentials.

11. A method in accordance with claim 9 wherein said time segments are within the range of the P wave up to the beginning of the QRS complex for detecting potentials of the atrium and of the bundle of His.

12. A method in accordance with claim 9 wherein said time segments are within the range of the P wave, the QRS complex itself and the ST portion for detecting rejection reactions after a heart transplant.

13. A system for analyzing selected signal components in physiological measurement signals, comprising:

means for selecting a portion of a physiological measurement signal to be analyzed;

means for selecting a plurality of substantially overlapping mutually offset time segments of said portion;

means for deriving a plurality of low order frequency functions representing the low-band spectral components in each respective time segment;

means for deriving a plurality of high order frequency functions representing the broadband spectral components in each respective time segment;

means for deriving a plurality of difference frequency functions by subtracting said low order frequency functions from respective high order frequency functions;

means for comparing said difference frequency functions, the comparison indicating detailed analysis of selected signal components in physiological measurement signals; and

means for outputting a three-dimensional plot of said frequency functions indicating detailed analysis of selected signal components in physiological measurement signals.

14. A system in accordance with claim 13 wherein said portion of a physiological measurement signal to be analyzed is a portion of an electrocardiographic signal having a P wave, a QRS complex, and an ST portion obtained from a patient.

15. A system in accordance with claim 14 wherein said time segments are at the trailing edge of the QRS complex and the ST portion for detecting late potentials.

16. A system in accordance with claim 14 wherein said time segments are within the range of the P wave up to the beginning of the QRS complex for detecting potentials of the atrium and of the bundle of His.

17. A system in accordance with claim 14 wherein said time segments are within the range of the P wave, the QRS complex itself and the ST portion for detecting rejection reactions after a hear transplant.

18. A system in accordance with claim 13 wherein said means for deriving low order and high order frequency functions each comprise means for autoregressive computation.

19. A system in accordance with claim 13 wherein said means for deriving low order and high order frequency functions each comprise means for maximum entropy method (MEM) computation as a method of autoregressie computation for carrying out spectral analysis of the individual time segments.

20. A system in accordance with claim 19 wherein said means for MEM computation comprises:

means for deriving a plurality of MEM coefficients; and

means for tapering said plurality of MEM coefficients, thereby reducing frequency shifting and splitting-up of the frequency peaks.

21. A system in accordance with claim 13 wherein said means for deriving low order and high order frequency functions each comprise means for adaptive filter determination (AFD) computation as a method of autoregressive computation for carrying out spectral analysis of the individual time segments.

22. A system in accordance with claim 21 wherein said means for AFD computation comprises:

means for deriving a plurality of AFD coefficients; and

means for predetermining the optimal number of said AFD coefficients.