ES2232223B1

ES2232223B1 - METHOD FOR THE DETERMINATION OF THE FREQUENCY OF WAVE FORM OF AN ECG SIGNAL.

Info

Publication number: ES2232223B1
Application number: ES200200710A
Authority: ES
Inventors: Jesus Ruiz Oieda; Elisabete Aramendi; Andoni Lazkano; Luis A. Leturiondo; Silvia Almaraz
Original assignee: OSATU S COOP
Current assignee: OSATU S COOP
Priority date: 2002-03-26
Filing date: 2002-03-26
Publication date: 2006-02-16
Anticipated expiration: 2022-03-26
Also published as: ES2232223A1

Abstract

El método (1) está destinado a la detección de la frecuencia de forma de onda (FFO) de señales ECG generadas por ritmos cardiacos tal como la fibrilación ventricular cuya FFO es de alto valor, y está basado en el análisis (22) de la función de autocorrelación estimada de una ventana de señal de larga duración. El método comprende una segunda etapa (25,26,27) de estimación de la autocorrelación, ahora mediante dos ventanas (Tw1,Tw2) de señal parciales, que es ejecutada para resolver un suceso de incertidumbre (23a) cuando el primer pico detectado no es el de mayor amplitud, y asegurar la transmisión (24) de un valor cierto de FFO como parámetro de actuación a un aparato de desfibrilación automático (DEA).The method (1) is intended for the detection of the waveform frequency (FFO) of ECG signals generated by heart rhythms such as ventricular fibrillation whose FFO is of high value, and is based on the analysis (22) of the Estimated autocorrelation function of a long-lasting signal window. The method comprises a second stage (25,26,27) of autocorrelation estimation, now by means of two partial signal windows (Tw1, Tw2), which is executed to solve an uncertainty event (23a) when the first detected peak does not it is the one with the greatest amplitude, and ensure the transmission (24) of a certain FFO value as an actuation parameter to an automatic defibrillation apparatus (DEA).

Description

Método para la determinación de la frecuencia de forma de onda de una señal de ECG.Method for determining the frequency of Waveform of an ECG signal.

La presente invención se relaciona con la detección de la frecuencia de la forma de onda (FFO) de la señal de un electrocardiograma (ECG), determinando un valor cierto de FFO mediante un cálculo estadístico de autocorrelación de la señal ECG para su aplicación a ritmos incoherentes de fibrilación ventricular.The present invention relates to the waveform frequency (FFO) signal detection an electrocardiogram (ECG), determining a true FFO value by means of a statistical calculation of autocorrelation of the ECG signal for application to inconsistent fibrillation rhythms ventricular

Prior state of the art

Las arritmias ventriculares como la fibrilación ventricular FV, y la taquicardia ventricular rápida TV, producen señales de ECG cuya frecuencia de forma de onda (FFO) es de valor alto, comparado con la forma de onda de otras arritmias.Ventricular arrhythmias such as fibrillation ventricular PV, and rapid ventricular tachycardia TV, produce ECG signals whose waveform frequency (FFO) is of value high, compared to the waveform of other arrhythmias.

La fibrilación ventricular se caracteriza por producir un ritmo incoherente de latidos cuyo ECG es una señal irregular sin la presencia de las ondas Q-R-S cíclicas de una señal normal, y con unas oscilaciones de amplitud y frecuencia aleatorias, en las cuales aunque no se presente una frecuencia cardiaca, se puede detectar una forma de onda a la que se le puede atribuir una frecuencia en latidos por minuto, cuyo valor es determinable para ser utilizado como un parámetro entre otros para la discriminación de una arritmia ventricular rápida.Ventricular fibrillation is characterized by produce an incoherent rhythm of beats whose ECG is a signal irregular without the presence of waves Cyclic Q-R-S of a normal signal, and with random oscillations of amplitude and frequency, in the which although a heart rate is not present, you can detect a waveform that can be attributed a frequency in beats per minute, whose value is determinable for be used as a parameter among others for discrimination of a rapid ventricular arrhythmia.

El valor de la FFO de la señal ECG suele medirse mediante cálculos de la correlación de la señal ECG registrada, una vez filtrada y digitalizada, y delimitada por una ventana de tiempo, que al menos comprende un latido. Mientras que la FFO de un ECG lento como el correspondiente a una bradicardia sinusal es de 30 latidos/min, un ritmo desfibrilable puede tener una FFO de más de 500 latidos/min, por lo que han de ser analizadas ventanas de larga duración.The FFO value of the ECG signal is usually measured by calculations of the correlation of the recorded ECG signal, a once filtered and digitized, and delimited by a time window, which at least comprises a beat. While the FFO of an ECG slow as that corresponding to sinus bradycardia is 30 beats / min, a defibrilable rhythm can have an FFO of more than 500 beats / min, so long windows have to be analyzed duration.

US-A-4403184 se divulga un método para la determinación de la frecuencia o el período de la oscilación básica de una señal aproximadamente periódica con componentes espectrales distribuidos estadísticamente, tal como una señal biológica. El valor de la FFO es determinado mediante un análisis de la similitud de la estructura de la señal de un ciclo a otro, usando el cálculo de la función de autocorrelación. El intervalo de tiempo medido desde el pico de correlación inicial hasta el siguiente pico sucesivo, es representativo del periodo o frecuencia de la señal investigada. En U-4463425 el método descrito para la determinación del período o la FFO, comprende una sola etapa de estimación de la función de autocorrelación, y siempre es juzgado el pico de mayor amplitud como pico verdadero. Previamente es determinado un nivel umbral del pico potencialmente verdadero, igual al 50% de un pico verdadero medido recientemente. Aunque en este método conocido, sean analizados solamente los picos positivos de correlación más grandes, esto es que exceden de un valor umbral alto, ocurren sucesos en la detección de señales de fibrilación ventricular, en los que la determinación de la FFO más representativa de la señal por el pico más alto detectado puede ser incierta.US-A-4403184 SE discloses a method for determining the frequency or the period of the basic oscillation of a signal approximately periodically with statistically distributed spectral components, such as a biological signal. The value of the FFO is determined through an analysis of the similarity of the signal structure from one cycle to another, using the calculation of the function of autocorrelation. The time interval measured from the peak of initial correlation to the next successive peak, is representative of the period or frequency of the investigated signal. In U-4463425 the method described for the determination of the period or the FFO, includes a single stage of estimation of the autocorrelation function, and the highest peak is always judged amplitude as true peak. Previously a level is determined potentially true peak threshold, equal to 50% of a peak true recently measured. Although in this known method, they are analyzed only the largest positive correlation peaks, this is that they exceed a high threshold value, events occur in the detection of ventricular fibrillation signals, in which the determination of the most representative FFO of the signal by the peak Higher detected may be uncertain.

Existen ritmos no ventriculares en los que la separación entre latidos no es uniforme o la forma de un latido es discordante con el resto de latidos, dando lugar a una curva de autocorrelación con un primer pico deformado y de poca amplitud, por lo que la posición relativa del pico más alto detectado, puede inducir a la determinación errónea del valor de la FFO representativo. También ocurre en algunas señales de ECG generadas por ritmos de fibrilación ventricular, que la curva de autocorrelación calculada desde una ventana de larga duración como 2,4 segundos, presenta al menos dos picos poco diferenciados uno del otro en amplitud. La determinación del pico representativo entre ambos picos por su posición relativa respecto al origen en tiempo, puede llevar a un valor de la FFO erróneo.There are non-ventricular rhythms in which the beat separation is not uniform or the shape of a beat is jarring with the rest of beats, resulting in a curve of autocorrelation with a first deformed peak and low amplitude, by what the relative position of the highest peak detected, can induce the wrong determination of the FFO value representative. It also occurs in some generated ECG signals by ventricular fibrillation rhythms, which the curve of autocorrelation calculated from a long-term window such as 2.4 seconds, has at least two slightly differentiated peaks one of the other in breadth. The determination of the representative peak between both peaks for their relative position with respect to the origin in time, It can lead to a wrong FFO value.

El valor de la FFO ha de ser determinado con certeza, puesto que será utilizado para la actuación de un aparato desfibrilador externo automático (DEA), en donde el shock eléctrico emitido por el actuador del aparato sería perjudicial para un paciente cuya arritmia no fuera del tipo desfibrilable. En la estimación de valores FFO de señales ECG se presenta una gama de valores tan amplia que obliga a delimitar ventanas de detección de larga duración, tal como 2,4 - 4 segundos. Consecuentemente las curvas de correlación obtenidas desde una ventana pueden resultar con un grado de incertidumbre que es necesario resolver antes de la transmisión de un valor de FFO definitivo al aparato DEA.The value of the FFO must be determined with certainty, since it will be used for the performance of a device automatic external defibrillator (DEA), where the electric shock issued by the device actuator would be detrimental to a patient whose arrhythmia was not the defibrilable type. In the FFO value estimation of ECG signals presents a range of values so wide that it forces to delimit detection windows of Long duration, such as 2.4-4 seconds. Consequently the correlation curves obtained from a window can result with a degree of uncertainty that needs to be resolved before the transmission of a definitive FFO value to the DEA device.

Exhibition of the invention

El objeto de la invención es un método para la determinación de un valor de la frecuencia de la forma de onda (FFO) de una señal de electrocardiograma (ECG), mediante el cálculo de la función de autocorrelación de una ventana de señal de larga duración, y la determinación de un valor cierto de la FFO como parámetro representativo de la señal ECG, a partir de la posición relativa en tiempo desde el origen a uno de los picos de gran amplitud detectados, según se define en la reivindicación 1.The object of the invention is a method for determination of a waveform frequency value (FFO) of an electrocardiogram (ECG) signal, by calculating the autocorrelation function of a long signal window duration, and the determination of a true FFO value as representative parameter of the ECG signal, starting from the position relative in time from the origin to one of the great peaks amplitude detected, as defined in claim 1.

Mediante el método proporcionado se asegura que la medida de la FFO transmitida a un aparato desfibrilador externo automático (DEA), es representativa de la señal para la actuación del aparato DEA, en todos los sucesos ocurribles incluidos los de señales ECG generadas por arritmias ventriculares. El valor FFO transmitido es obtenido en correspondencia con la posición desde el origen en tiempo a un primer pico de autocorrelación de máxima amplitud en dicha primera estimación de autocorrelación, o cuando un suceso es discriminado como incierto, el valor FFO se determina en una segunda etapa de estimación de la autocorrelación por el primer pico detectado de amplitud mayor que el 50% del pico máximo posterior.The method provided ensures that the measurement of the FFO transmitted to an external defibrillator device automatic (DEA), is representative of the signal for the performance of the DEA apparatus, in all the possible events including those of ECG signals generated by ventricular arrhythmias. FFO value transmitted is obtained in correspondence with the position from the origin in time at a maximum peak autocorrelation peak amplitude in said first autocorrelation estimate, or when an event is discriminated as uncertain, the FFO value is determined in a second stage of autocorrelation estimation by the first peak detected of amplitude greater than 50% of the maximum peak later.

El método comprende un análisis comparativo de los sucesivos picos de la curva de autocorrelación de una señal de larga duración. Puesto que el resultado del análisis de los picos de correlación puede presentarse como incierto, cuando son detectados dos o más picos grandes próximos, el método comprende además un algoritmo de discriminación de la certeza de uno de los picos de gran amplitud detectados para la determinación de la frecuencia de la señal investigada. La discriminación del pico de correlación cierto entre dos picos sucesivos de gran amplitud se resuelve mediante una segunda etapa de análisis, ahora de curvas de autocorrelación estimadas desde ventanas parciales. Así se asegura que la medida de la FFO de algunos ritmos particulares de fibrilación ventricular y taquicardia ventricular es ciertamente la más representativa para su transmisión al aparato DEA.The method comprises a comparative analysis of the successive peaks of the autocorrelation curve of a signal from Long duration. Since the result of the analysis of the peaks of correlation can be presented as uncertain, when they are detected two or more nearby large peaks, the method further comprises a discrimination algorithm of the certainty of one of the peaks of large amplitude detected for the determination of the frequency of The signal investigated. The discrimination of the correlation peak true between two successive peaks of great amplitude is resolved through a second stage of analysis, now of curves of Estimated autocorrelation from partial windows. This ensures that the FFO measure of some particular rhythms of ventricular fibrillation and ventricular tachycardia is certainly the more representative for transmission to the DEA device.

Description of the drawings

Las figuras 1A y 2A muestran dos ejemplos de ventanas de señal de electrocardiograma (ECG) correspondientes a arritmias ventriculares.Figures 1A and 2A show two examples of electrocardiogram (ECG) signal windows corresponding to ventricular arrhythmias

Las figuras 1B y 2B muestran gráficamente las funciones de autocorrelación calculadas desde las ventanas de señal de las figuras 1A y 2A respectivamente.Figures 1B and 2B graphically show the autocorrelation functions calculated from the signal windows of figures 1A and 2A respectively.

La figura 3 es el diagrama de etapas de un método seguido para la determinación del valor de la frecuencia de la forma de onda de las señales ECG de las figuras 1A y 2A.Figure 3 is the step diagram of a method followed for the determination of the value of the frequency of the form waveform of the ECG signals of Figures 1A and 2A.

Las figuras 4 y 5 son diagramas parciales de algoritmos incluidos en el diagrama de la figura 3.Figures 4 and 5 are partial diagrams of algorithms included in the diagram in figure 3.

Detailed description of the preferred embodiment

En referencia a las figuras 1-5, un método 1 para la determinación de un valor FFO de la frecuencia de forma de onda de dos ejemplos de diferentes señales 2,3 de ECG (fig. 1A y 2A), es representado en la figura 3. El valor FFO determinado como cierto es transmitido finalmente a un aparato desfibrilador externo automático DEA, el cual actúa generando shocks eléctricos en función del valor FFO transmitido. La señal ECG 2,3 es filtrada antialiasing y muestreada en un convertidor A/D, etapa 20, procesada preciamente mediante un filtro paso banda por ejemplo 0,7-35 Hz, en ventanas 5,6 de tiempo "t" de duración 2,4 s, siendo suprimidas las interferencias que acompañan a la señal ECG debidas a movimientos del paciente o al de los electrodos y a la red de 50 Hz. De una ventana 5,6 se realiza, etapa 21, una estimación sesgada de la función de autocorrelación de la señal 2,3 adquirida, representada por una de las curvas de autocorrelación 7,8 mostradas en las figuras 1B y 2B, siendo detectados y analizados los picos positivos de autocorrelación, etapa 22, que superan un valor umbral prefijado o "umbral positivo" y del resultado obtenido se resuelve la condición de certeza, etapa 24.Referring to figures 1-5, a method 1 for the determination of a frequency FFO value waveform of two examples of different 2.3 ECG signals (fig. 1A and 2A), is represented in figure 3. The FFO value determined as true is finally transmitted to a device DEA automatic external defibrillator, which acts by generating shocks electrical depending on the transmitted FFO value. The ECG 2.3 signal is filtered antialiasing and sampled in an A / D converter, stage 20, processed precisely by a bandpass filter for example 0.7-35 Hz, in windows 5.6 of time "t" of duration 2.4 s, the accompanying interference being suppressed to the ECG signal due to patient movements or to those of the electrodes and the 50 Hz network. From a window 5.6 is performed, stage 21, a biased estimate of the autocorrelation function of the acquired signal 2.3, represented by one of the curves of autocorrelation 7.8 shown in Figures 1B and 2B, being positive autocorrelation peaks detected and analyzed, step 22, which exceed a preset threshold value or "threshold positive "and from the result obtained the condition of certainty, stage 24.

El pico 9 del origen "t0" de tiempo es el de máxima amplitud o altura, pues corresponde siempre al índice de correlación máxima. En la etapa 22 de análisis de la estimación de la autocorrelación, se detectan los picos 10,11,12, de la curva de autocorrelación 7, 8 estimada, y se registran con su amplitud y su posición desde el origen "t0" de la ventana 5,6. Son posibles cuatro sucesos diferentes (fig. 3), "CASO 1", "CASO 2", "CASO 3", e "INCERTIDUMBRE"(23a):Peak 9 of the origin "t0" of time is that of maximum width or height, since it always corresponds to the index of maximum correlation In step 22 of analysis of the estimate of autocorrelation, peaks 10,11,12, of the curve of autocorrelation 7, 8 estimated, and recorded with its amplitude and its position from the origin "t0" of the window 5,6. They're possible four different events (fig. 3), "CASE 1", "CASE 2", "CASE 3", and "UNCERTAINTY" (23rd):

CASO 1: Ha sido detectado un pico único en el origen t0. Se determina que FFO=0 pues la ventana 5,6 de la muestra de señal ECG presenta un latido o ninguno.CASE 1: A single peak has been detected in the origin t0. It is determined that FFO = 0 as window 5.6 of the sample ECG signal presents a beat or none.

CASO 2: Asociado a una curva 7 como la de fig. 1B. El pico de mayor altura excluido el pico 9 del origen es el pico 10 primero, consecutivo al pico 9 del origen. El tiempo de la posición "t1" relativa al origen del pico 10 más alto determina la FFO de la señal 2 sometida a medición.CASE 2: Associated with a curve 7 like the one in fig. 1 B. The highest peak excluding peak 9 of the origin is the peak 10 first, consecutive to peak 9 of the origin. The time of the position "t1" relative to the origin of the highest peak 10 determines the FFO of the signal 2 subjected to measurement.

CASO 3: (no representado en los dibujos). El pico más alto es el segundo o uno posterior, y no hay otro pico intermedio de amplitud mayor que el 50% de aquél.CASE 3: (not shown in the drawings). Peak higher is the second or later one, and there is no other peak intermediate amplitude greater than 50% of that.

INCERTIDUMBRE: asociada a una curva 8 como la de fig. 2B: El pico más alto es el pico 12 segundo o uno posterior, y se ha detectado algún pico intermedio 11 de amplitud mayor que el 50% de aquél.UNCERTAINTY: associated with a curve 8 like that of fig. 2B: The highest peak is the 12 second or later peak, and some intermediate peak 11 of amplitude greater than the 50% of that.

En los sucesos CASO 2 y CASO 3 se determina la FFO de la señal 2,3 por la posición t1, relativa al origen del tiempo "t" del pico 10 más alto, siendo transmitido, etapa 23, al aparato DEA como valor cierto de la FFO. Si el resultado del análisis de la autocorrelación es de "INCERTIDUMBRE", no se tiene la certeza de que la posición del pico 12 más alto corresponda al valor cierto de la FFO. Para su discriminación, etapa 23, ahora se continúa con una segunda etapa (fig. 5) de estimación de la autocorrelación siendo realizados dos estimaciones parciales, etapas 25,26, de cada uno de dos segmentos Tw1,Tw2 o mitades de ventana de duración 1,2 s, y en consecuencia un análisis detallado, etapa 27, de los picos detectados en ambas estimaciones de autocorrelación (no representadas en los dibujos), para buscar en su caso, etapa 28, otro valor de FFO más cierto que el valor previamente determinado en el suceso de INCERTIDUMBRE.In the CASE 2 and CASE 3 events, the FFO of signal 2.3 by position t1, relative to the origin of the time "t" of the highest peak 10, being transmitted, stage 23, to the DEA device as the true value of the FFO. If the result of Autocorrelation analysis is "UNCERTAINTY", it is not you are certain that the position of the highest peak 12 corresponds to the true value of the FFO. For your discrimination, stage 23, now a second stage (fig. 5) of estimation of the autocorrelation being made two partial estimates, stages 25.26, of each of two segments Tw1, Tw2 or window halves of duration 1.2 s, and consequently a detailed analysis, step 27, of the peaks detected in both autocorrelation estimates (not shown in the drawings), to look for the case, step 28, another FFO value more true than the previously determined value in the UNCERTAINTY event.

Durante el segundo análisis 27 (fig. 5) se detectan los picos positivos de correlación cuya amplitud excede de un primer valor umbral prefijado diferente que el "umbral positivo" de la primera estimación (etapa 21 de fig. 3) de la autocorrelación, y son clasificados como picos de gran amplitud o picos grandes excluido el pico de origen 9, que exceden de un segundo umbral prefijado. Se registra la amplitud de cada uno de los picos grandes y su posición desde el origen "t0" del tiempo, en cada segmento Tw1,Tw2 de la ventana 6 de señal.During the second analysis 27 (fig. 5), detect positive correlation peaks whose amplitude exceeds a first preset threshold value different than the "threshold positive "of the first estimate (step 21 of fig. 3) of the autocorrelation, and are classified as wide-ranging peaks or large peaks excluding the peak of origin 9, which exceed a Second preset threshold. The amplitude of each of the large peaks and their position from the origin "t0" of time, in each segment Tw1, Tw2 of the signal window 6.

Los valores umbrales de los picos de gran amplitud detectados en cada una de las dos estimaciones, son prefijados experimentalmente a partir del registro y el procesado de gran cantidad de ritmos cardíacos, incluyendo un número elevado de fibrilaciones ventriculares. La posición del pico más alto proporcionará una medida FF01.The threshold values of large peaks amplitude detected in each of the two estimates, are experimentally preset from the registration and processing of large number of heart rhythms, including a high number of ventricular fibrillations The highest peak position will provide a measure FF01.

En la curva de autocorrelación de una primera mitad Tw1 de ventana de señal (no representada en los dibujos), si el pico más alto es el primero tras el origen t0, su posición proporciona una medida FF01. Si el pico más alto es el segundo o posterior, se compara con los picos vecinos, y si ninguno de estos últimos es mayor que el 50%, la posición del pico más alto proporciona una medida FF01, mientras que si hay un pico vecino mayor que el 50% del pico mas alto, es la posición del pico vecino la que proporcionará la medida FF01. Si no hay picos de gran amplitud en la estimación de la autocorrelación de la primera mitad Tw1 de la ventana, se analiza si el número de picos excluido el origen es igual o superior a cuatro. Si esta condición se cumple, se mide un valor FF02 correspondiente a la posición del pico más cercano al origen t0.In the autocorrelation curve of a first Tw1 half of signal window (not shown in the drawings), if the highest peak is the first after the origin t0, its position provides a measure FF01. If the highest peak is the second or later, it is compared with the neighboring peaks, and if none of these last is greater than 50%, the highest peak position provides a measure FF01, while if there is a neighboring peak greater than 50% of the highest peak, is the position of the neighboring peak which will provide measure FF01. If there are no great peaks amplitude in the first half autocorrelation estimate Tw1 of the window, it is analyzed if the number of peaks excluded the origin is equal to or greater than four. If this condition is met, it will be measures a value FF02 corresponding to the position of the peak plus close to the origin t0.

Este mismo proceso lógico es seguido en el análisis de la curva de autocorrelación de la segunda mitad Tw2 de la ventana 6, proporcionando respectivamente una medida FF03 o FF04, dependiendo de la detección de uno o dos picos de gran amplitud.This same logical process is followed in the analysis of the autocorrelation curve of the second half Tw2 of window 6, respectively providing a measure FF03 or FF04, depending on the detection of one or two peaks of great amplitude.

Después de calculados FF01, FF02, FF03 y FF04 pueden presentarse dos sucesos adicionales "CASO 4" y "CASO 5" para la determinación de un valor cierto de FFO:After calculated FF01, FF02, FF03 and FF04 two additional events "CASE 4" and "CASE can be presented 5 "for the determination of a true FFO value:

CASO 4: Si ha sido medido al menos uno de los dos valores FF01, FF03, el valor cierto de FFO para ser transmitido al aparato DEA es el mayor de ellos dos. Si no hay picos detectados se determina como valor cierto de FFO el más alto de FF02 y FF04.CASE 4: If at least one of the two has been measured values FF01, FF03, the true value of FFO to be transmitted to DEA apparatus is the largest of them two. If there are no peaks detected, determine the highest value of FFO as the highest of FF02 and FF04.

CASO 5: Si ninguna de las medidas FF01-FF04 ha sido posible, el valor cierto de FFO se determina desde la curva 8 de autocorrelación calculada de la ventana 6 previa de 2.4 s de la primera estimación (21 de la fig. 4). La posición del pico 11 detectado más próximo al origen t0, cuya amplitud era mayor que el 50% del pico 12 más alto posterior (fig. 2B) corresponde al valor cierto de la FFO.CASE 5: If none of the measures FF01-FF04 has been possible, the true value of FFO is determined from the calculated autocorrelation curve 8 of the previous window 6 of 2.4 s of the first estimate (21 of fig. 4). The position of the detected peak 11 closest to the origin t0, whose amplitude was greater than 50% of the highest posterior peak 12 (fig. 2B) corresponds to the true value of the FFO.

Claims

1. Method for determining the waveform frequency (FFO) of an electrocardiogram signal (2,3) (ECG), as a measure of an actuation parameter transmitted (24) to an automatic external defibrillator device (AED) , which includes the acquisition (20) of the digitized (2,3) ECG digitized and filtered bandpass, the estimation (21) of a signal autocorrelation function (2,3) in at least one window (5,6 ) of long duration (t), the detection and recording (22) of the sequence of the positive (9-12) correlation curves of the estimated curve (7.8), and the determination (24) of a representative FFO value of the signal (2,3) from the position (t1) relative to the time source (t0) of the window of a large peak (9-12) previously discriminated as representative of the FFO of other nearby peaks also detected, characterized in that it includes a discrimination process (23a) with certainty of said peak (10-12) of correlation corresponding to the representative FFO, by means of a second stage (25,26) of estimation of the autocorrelation of the signal (2,3) carried out in partial windows (Tw1, Tw2) of the window (6) of long duration.

2. Method for determining the frequency according to claim 1, characterized in that in said second stage (25,26,27) of estimation and autocorrelation analysis the positive correlation peaks detected are recorded and also separately the large amplitude peaks that exceed a preset threshold value, where the true value of the FFO is determined between four possible FFO values (FF01-FF04) calculated from both partial windows (Tw1, Tw2), by the position of a large amplitude correlation peak if any, or the position of the first positive peak after the origin (t0) of time between at least three positive peaks detected successively.

3. Method for determining the frequency according to claim 1, characterized in that in the second stage (25,26) of autocorrelation estimation the positive correlation peaks are recorded and separately the large amplitude peaks exceeding a predetermined threshold value , and when a correlation peak is not detected for the calculation of a representative FFO value of the signal (3), between two of said partial windows (Tw1, Tw2), a certain FFO value is determined from the first estimate ( 21) of long-term autocorrelation, by a first peak (11) with respect to the origin (t0) whose amplitude is close to that of the subsequent peak (12) of maximum amplitude.