FR2495330A1 - Systeme de mesure de periode - Google Patents

Systeme de mesure de periode Download PDF

Info

Publication number
FR2495330A1
FR2495330A1 FR8113883A FR8113883A FR2495330A1 FR 2495330 A1 FR2495330 A1 FR 2495330A1 FR 8113883 A FR8113883 A FR 8113883A FR 8113883 A FR8113883 A FR 8113883A FR 2495330 A1 FR2495330 A1 FR 2495330A1
Authority
FR
France
Prior art keywords
period
autocorrelation function
value
peak
biosignal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
FR8113883A
Other languages
English (en)
Other versions
FR2495330B1 (fr
Inventor
Toshinori Hirano
Masakazu Murase
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Terumo Corp
Original Assignee
Terumo Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from JP9782180A external-priority patent/JPS5722737A/ja
Priority claimed from JP9782080A external-priority patent/JPS5722736A/ja
Application filed by Terumo Corp filed Critical Terumo Corp
Publication of FR2495330A1 publication Critical patent/FR2495330A1/fr
Application granted granted Critical
Publication of FR2495330B1 publication Critical patent/FR2495330B1/fr
Granted legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/43Detecting, measuring or recording for evaluating the reproductive systems
    • A61B5/4306Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
    • A61B5/4343Pregnancy and labour monitoring, e.g. for labour onset detection
    • A61B5/4362Assessing foetal parameters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02411Detecting, measuring or recording pulse rate or heart rate of foetuses
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/02Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • GPHYSICS
    • G04HOROLOGY
    • G04FTIME-INTERVAL MEASURING
    • G04F10/00Apparatus for measuring unknown time intervals by electric means
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S128/00Surgery
    • Y10S128/92Computer assisted medical diagnostics

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Cardiology (AREA)
  • Pregnancy & Childbirth (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Gynecology & Obstetrics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Pediatric Medicine (AREA)
  • Reproductive Health (AREA)
  • Physiology (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

Système de mesure de période. Ce système est prévu pour échantillonner un biosignal à une période prédéterminée d'échantillonnage, trouver une fonction d'autocorrélation pour une variable tau à partir du biosignal échantillonné, et trouver ensuite une fonction d'autocorrélation qui correspond à la valeur de la variable de déphasage obtenue en modifiant la variable tau le long de l'axe des temps. La fonction d'autocorrélation obtenue de cette manière est stockée dans une mémoire, puis comparée avec la fonction d'autocorrélation suivante. La comparaison est répétée pour plusieurs fonctions successives d'autocorrélation ; de ce fait, on peut trouver une pointe des fonctions d'autocorrélation, ce qui permet de mesurer la période du biosignal. Application : notamment pour mesurer la période d'un signal représentant la pulsation du cîoeur d'un fîoetus. (CF DESSIN DANS BOPI)

Description

Système de mesure de période.
La présente invention se rapporte à un systèmede mesure de la pé-
riode d'un biosignal, et particulièrement d'un signal représentant la pulsa-
tion du coeur d'un foetus. Un système classique pour mesurer la période d'un biosignal est
basé sur un sytème de corrélation adapté pour dériver une fonction d'auto-
corrélation du biosignal et pour mesurer la période du biosignal en se ba-
sant sur la fonction d'autocorrélation.
Le système de mesure de période qui est basé sur le système de cor-
relation fonctionne en échantillonnant un biosignal sur une période appro-
priée d'échantillonnage, en calculant la fonction d'autocorrélation du bio-
signal à partir des données échantillonnées, et en détectant les pointes du
biosignal à partir de l'autocorrélation calculée afin d'en obtenir la pé-
riode.
La fonction d'autocorrélation indique la similarité entre deux par-
ties de la forme d'onde d'un biosignal à deux moments différents séparés par un certain intervalle de temps. En d'autres termes, elle représente le degré de similarité de la forme d'onde répétitrice du biosignal. On comprendra mieux cette notion en se reportant à la figure 1 sur laquelle on peut voir que si une partie M1, qui se répète à une certaine période T, est décalée sur l'axe des temps d'un intervalle de temps qui est égal à la période T, cette partie M1 se superposera avec le maximum de précision- à la partie M2
qui lui fait suite immédiatement.
Afin d'obtenir la fonction d'autocorrélation à partir du biosignal,
on peut écrire cette fonction d'autocorrélation A (e) en fonction du bio-
signal f(t) qui est lui-même fonction du temps t. Ainsi, A ('r) peut s'écri-
re: i IT
A( ')=im -
A ('): - 2 T f(t) f(t +Y) dt..... (1)
expression dans laquelle T représente la période du biosignal et 'un inter-
valle donné de temps entre deux points séparés dans le temps, le premier
point dans le temps étant un temps de référence en liaison avec le biosignal.
En d'autres termes, t-est une variable qui applique un déphasage au biosi-
gnal f(t) sur l'axe des temps.
On se reportera maintenant à la figure 2 pour décrire le système classique de mesure de période qui est basé sur la fonction de corrélation
pour mesurer la période d'un biosignal, en particulier d'un signal repré-
sentant la pulsation du coeur d'un foetus, lequel signal sera appelé ci-
après "signal de pulsation du coeur".
Comme le représente la figure 2, on place une sonde 2 en contact avec l'abdomen d'un sujet féminin pour détecter le signal de pulsation du coeur du foetus afin de le mesurer. La forme d'onde du signal de pulsation du coeur ainsi détectée est traitée de manière appropriée dansun circuit 3 de pré-traitement, et le signal est ensuite échantillonné à une période prédéterminée d'échantillonnage dans un circuit 4 d'échantillonnage. Les données obtenues par l'échantillonnage du signal de pulsation du coeur sont stockées dans une mémoire 6 de données composée de plusieurs registres à décalage. Au fur et à mesure que chaque élément nouveau d'information est introduit dans la mémoire 6 de données, les éléments d'information déjà
stockés jusqu' alors sont décalés vers le registre immédiatement con-
tigu, de sorte que les données sont décalées séquentiellement d'un regis-
tre à l'autre, les éléments les plus anciens d'information stockés dans le
dernier registre étant perdus à chaque nouvelle entrée de données. Un multi-
plicateur 8 et un additionneur 10 constituent un circuit de calcul de fonc-
tions d'autocorrélation, qui est prévu pour calculer une fonction d'auto-
corrélation à l'aide des données stockées dans la mémoire 6 de données. Une mémoire 12 de corrélation stocke les résultats du calcul, c'est-àdire la fonction d'autocorrélation calculée. La fonction d'autocorrélation est donc
calculée par le multiplicateur 8 et l'additionneur 10 sur la base des don-
nées stockées dans la mémoire 6 de données. Le calcul est effectué sur la
base de divisions pour un seul cycle d'échantillonnage et, pour chaque élé-
ment d'information Xi, X2, X3 --., il se poursuit de la manière suivante (Xi, XSA+ A1 --> AI, (X1I XS2)+ A2 --- A2.*. (XI.XS+m)+ Am -- Am, le résultat de chaque calcul étant stocké séquentiellement dans la mémoire 12 de corrélation. En répétant ces calculs et ces opérations de stockage pour
n cycles, les données définissant la fonction d'autocorrélation sont stoc-
kées dans la mémoire 12 de corrélation. Les pointes représentant la pério-
dicité de la fonction d'autocorrélation stockée dans la mémoire 12 de corré-
lation sont détectées par un détecteur 14 de pointes afin de pouvoir obte-
nir la période du biosignal.
Dans le système classique de mesure du type décrit, cependant, l'arrangement est tel que la variable r de déphasage varie dans chaque
cycle d'échantillonnage. Il faut donc stocker dans la mémoire 12 de corré-
lation les résultats de chaque calcul de fonction d'autocorrélation englo-
bant l'ensemble des données qui couvrentla plage dans laquelle la variable G' varie dans chaque cycle d'échantillonnage. Ceci signifie que la mémoire de corrélation doit posséder-une très grande capacité de stockage. En outre, même lorsqu'on mesure un signal présentant une période courte, les calculs décrits ci-dessus couvrent un intervalle de temps-correspondant à deux ou trois fois la longueur de la période, de sorte qu'une grande partie de ce
calcul est sans signification sensible. Ce fait nécessite également une mé-
moire de corrélation avec une grande capacité de stockage et il est également défavorable au point de vue du traitement en temps réel du fait qu'un
grand nombre de calculs sensiblement sans signification sont effectués.
En conséquence, la présente invention a pour objet de fournir un système pour mesurer la période d'un biosignal, système qui ne présente pas
les défauts mentionnés ci-dessus de sorte qu'il permet de mesurer la pério-
de avec une mémoire de corrélation dont la capacité de stockage est plus
petite et avec un temps de calcul qui est réduit au minimum possible.
Un autre objet de la présente invention est de fournir un système
de mesure de période qui permet de mesurer correctement la période en détec-
tant les pointes réelles qui correspondent à la période d'un biosignal, à
partir de plusieurs pointes obtenues depuis une fonction d'autocorréla-
tion. Pour réaliser ces objets, la présente invention fournit un système de mesure de période qui-comprend un moyen pour extraire un biosignal, un moyen pour calculer la fonction d'autocorrélation du biosignal, un moyen pour détecter une pointe à partir de la fonction d'autocorrêlation, et un moyen pour-calculer la période du biosignal à partir du point sur un axe de
corrélation auquel est détectée une pointe par le moyen de détection de poin-
tes, le calcul de la-fonction d'autocorrêlation étant continu dans un inter-
valle qui correspond sensiblement à la valeur minimale de la période de me-
sure, cet intervalle commençant avec la détection d'une pointe, étant bien
établi qu'aucune pointe supérieure à la pointe détectée n'existe dans l'in-
tervalle qui correspond à-la valeur minimale et qui commence avec la détec-
tion de la pointe, afin que la pointe détectée soit une pointe réelle. La présente invention étant considérée sous un autre aspect, la fonction de corrélation, qui est donnée par l'équation A <S) = nL f(k + 't) pour
une certaine valeur d'une variable e qui applique un déphasage au biosi-
gnal sur l'axe des temps, est calculée par le moyen de calcul des fonctions d'autocorrélation pour une valeur spécifique de la variable Cde déphasage, la valeur spécifique de la variable Dt de déphasage est avancée sur l'axe des temps pour se conformer à la progression des cycles d'échantillonnage, de ce fait le moyen de calcul des fonctions d'autocorrélation calcule les
fonctions d'autocorrélation les unes après les autres conformément aux nou-
velles valeurs spécifiques de la variable de déphasage, et la valeur cal-
culée de la fonction d'autocorrélation est stockée dans la mémoire et com-
parée à la valeur calculée la plus récente de la fonction d'autocorrélation,
afin de détecter une pointe..
La présente invention sera bien comprise à la lecture de la des-
cription suivante faite en relation avec les dessins ci-joints, dans lesquels: - la figure 1 est un diagramme représentant la forme d'onde d'un biosignal, ce diagramme étant utilisé pour décrire la mesure d'une période au moyen d'un système d'autocorrélation; - la figure 2 est un schéma de fonctionnement représentant, sous une forme simplifiée, la construction d'un appareil de mesure de période auquel s'applique le système classique de mesure de période; - la figure 3 est une vue illustrative utilisée pour décrire la manière dont est calculée une fonction d'autocorrélation dans un système de mesure de période suivant la présente invention;
- la figure 4 est un diagramme montrant la forme d'onde d'un si-
gnal représentant la pulsation du coeur d'un foetus, ce diagramme étant utilisé pour décrire un cas o le signal de mesure de période suivant la
présente invention s'applique à la mesure de la période entre les pulsa-
tions du coeur d'un foetus; - la figure 5 est un diagramme représentant une forme d'onde, ce diagramme étant utilisé pour décrire un système prévu pour continuer le calcul des fonctions d'autocorrélation pendant une période fixe de temps qui suit la détection d'une pointe, afin d'établir si la pointe détectée est une pointe réelle ou non; - la figure 6 est un schéma de fonctionnement représentant, sous une forme simplifiée, la construction d'un appareil de mesure de période auquel s'applique le système de mesure de période qui fait l'objet de la présente invention; la figure 7 est un schéma de fonctionnement qui est utilisé pour décrire le stockage de données d'échantillonnage dans une mémoire de données, ainsi que la lecture et, ensuite, le traitement de ces données; - la figure 8 est un schéma de fonctionnement représentant la construction détaillée du détecteur de pointes, du circuit de contrôle du niveau de pointe et du circuit de confirmation de pointe incorporés à l'appareil de mesure de période représenté à la figure 6; et - la figure 9 est un schéma de fonctionnement qui est utilisé
pour décrire les détails d'un générateur de niveau de référence.
On se-reportera maintenant à la figure 3 qui est utilisée pour décrire le système de mesure de période suivant la présente invention et
qui illustre le système utilisé pour calculer la fonction d'autocorréla-
tion d'un biosignal.
Soit f (k) (k étant égal à 1, 2, 3... n) une fonction qui re-
présente les données obtenues par des opérations respectives d'échantillon- nage appliquées à un biosignal à une période fixe T5 d'échantillonnage, la
fonction d'autocorrélation A(t) du biosignal s'exprimera alors par l'é-
quation suivante: n A(t) =n i f(k). f(k +).(2) dans laquelle "I " représente une variable qui applique un déphasage au biosignal sur l'axe des temps, "n" le nombre total de multiplications ou d'additions dans un cycle d'échantillonnage et "k" un nombre ordinal d'échantillonnage. En développant l'équation (2), on obtient: A(%^) = 1t f(1). f(1 +l) + f(2). f(2 +T<) + f(3).f(3 + ' +... + f(n). f(n +Yt)}.. (3)
Dans l'équation (3), f(1) représente les données les plus ré-
centes. L'équation (3) signifie que la fonction d'autocorrélation d'un
biosignal s'obtient en additionnant "n" fois au total le produit f(k).
f(k +1e) et en changeant la valeur de k, le produit f(k). f(k +1t) étant le produit des données échantillonnées f(k) et f(k +tj) en deux
points du temps séparés sur l'axe des temps par la variable t:de déphasage.
En se reportant plus précisément à la figure 3, on supposera que
plusieurs éléments d'information sont obtenus par des opérations d'échan-
tillonnage effectuées à des intervalles égaux à la période d'échantillon-
nage TS représentée sur l'axe des temps, et que la variable -c de déphasage est donnée par "m". Pour calculer la fonction d'autocorrélation A (m),deux
éléments d'information d'échantillonnage séparés l'un de l'autre de la va-
leur "m", tels que f(1) et f(m + 1), f(2) et f(m + 2), f(n) et f(m + n) , sont multipliés entre eux pour donner les produits f(1). f(m + 1),
f(2). f(m + 2)..... f(n). f(m + n). Ces produits sont ensuite addition-
nés ensemble pour les "n" opérations d'échantillonnage effectuées pendant le
cycle d'échantillonnage, afin d'obtenir la fonction d'autocorrélation A(m) .
Le système adopté dans la présente invention calcule la fonction d'autocor-
rélation pour une certaine valeur de la variable 'tqui applique le déphasa-
ge au biosignal sur l'axe des temps pendant un cycle d'échantillonnage du biosignal, fait varier la valeur de la variable -û de déphasage sur l'axe des temps conformément à la progression des cycles d'échantillonnage, et
calcule ensuite la fonction d'autocorrélation qui correspond à chaque cy-
cle d'échantillonnage. Les résultats du dernier calcul de la fonction
d'autocorrélation sont stockés dans la mémoire, ce qui permet de déter-
miner les pointes et la période du signal.
On décrira maintenant plus en détail le système de mesure de pé-
riode en prenant comme exemple un cas o la présente invention est appliquée à la mesure de la période du signal représentant la pulsation du coeur
d'un foetus.
La période de pulsation du coeur d'un foetus est de l'ordre de
300 à 1 500 millisecondes environ. Par conséquent, pour calculer une fonc-
tion d'autocorrélation couvrant la période complète du signal de pulsa-
tion du coeur, il faut trouver la fonction d'autocorrélation en faisant varier la période de mesure depuis la valeur minimale de 300 millisecondes jusqu'à la valeur de 1 500 millisecondes. En d'autres termes, il faut faire varier la variable Itde déphasage dans l'intervalle de 3OO/Ts à 1500/Ts dans l'équation (2). Du fait que la fonction d'autocorrélation présentera une pointe maximum dans cet intervalle lorsque la variable jt de déphasage sera fixée à la période T du signal de pulsation du coeur ou à une période de temps qui sera un multiple de la période T, on pourra déterminer la période
réelle du signal de pulsation du coeur si la pointe correspondant à la pé-
riode qr est détectée.
Conformément au système de mesure de période qui fait l'objet de
la présente invention, le calcul de la fonction d'autocorrélation est effec-
tué en utilisant chaque cycle d'échantillonnage comme division unique. Géné-
ralement, la période la plus courte du signal de pulsation du coeur d'un
foetus est d'environ 300 millisecondes. Comme le montrera clairement l'ex-
plication donnée ci-dessous, le calcul de la fonction d'autocorrélation
commence par la plus petite valeur possible de la période de mesure, c'est-
à-dire 300 millisecondes, afin d'extraire les résultats de mesure qui cou-
vrent un intervalle de temps équivalent à la période. Ce qui revient à dire
que, pendant le premier cycle d'échantillonnage, la fonction d'autocorréla-
tion est obtenue d'abord pour ce qui concerne l'intervalle de 300 millise-
condes correspondant à la valeur minimale de la période de pulsation du coeur d'un foetus. Dans ce cas, la variable Et de déphasage est obtenue à partir de C: = 3OO/Ts, de sorte que cette variable Ir sera égale à 60 si On fixe la période TS d'échantillonnage à cinq millisecondes. Ensuite, avec une période TS d'échantillonnage de cinq millisecondes, le temps autorisé pour un calcul concernant les données échantillonnées sera inférieur à cinq millisecondes énviron. Par conséquent, "n" opérations d'échantillonnage seront effectuées aux conditions I2 = 60 et période d'échantillonnage TS millisecondes, et la fonction d'autocorrelation A(60) sera obtenue
pour D = 60. La fonction d'autocorrelation A(60) est obtenue par le pro-
cédé utilisé pour obtenir la fonction d'autocorrélation A('/) à la figu-
re 3.
On reprendra maintenant la description précédente en se reportant
à la figure 4 qui représente un signal de pulsation du coeur. L'échantillon-
nage s'effectue jusqu'à un total de "n" fois à intervalles de cinq milli-
secondes, ce qui équivaut à la période Ts d'échantillonnage (c'est-a-dire
à des intervalles définis par TS = 5 millisecondes). Les éléments d'infor-
mation f(1), f(2), f(3), f(4).... f(n) obtenus par chaque opération
d'échantillonnage sont stockés dans la mémoire. Ensuite, deux éléments d'in-
formation f(k) et f(k + 60) obtenus à deux moments différents d'échantil-
lonnage espacés l'un de l'autre de la variable de déphasage t = 60, sont
multipliés l'un par l'autre et une série de ces produits, tels que f(1).
f(1 + 60), f(2). f(2 + 60)........ sont additionnés ensemble pour donner
la somme totale des produits. On peut donc obtenir la fonction d'autocorré-
lation A(60) pour le cas oO la variable 'lde déphasage est réglee à 60.
La valeur de A(60) indique le degré de périodicité en rapport avec:= 60 (c'est-à-dire pour une période de 300 millisecondes). La valeur de A(60)
est stockée dans la mémoire aux fins de comparaison, jusqu'à ce que la fonc-
tion d'autocorrelation soit obtenue dans le cycle suivant d'échantillonnage.
Le calcul est effectué ensuite pour le deuxième cycle d'échantillon-
nage dans lequel la valeur de la variable de déphasage est avancée de "un"
c'est-a-dire à A(61). En d'autres termes, dans le deuxième cycle d'échan-
tillonnage, la fonction d'autocorrélation est calculée pour une période de
305 millisecondes. Le calcul de la fonction d'autocorrelation A(61) s'effec-
tue sensiblement de la même manière que le calcul de la fonction d'autocorré-
lation A(60) et on ne le décrira pas de nouveau. La fonction d'autocorrela-
tion A(61) obtenue par le calcul portant sur une période de 305 millisecon-
des, est comparée à la fonction d'autocorrelation A(60) portant sur la pé-
riode de 300 millisecondes et précédemment calculée et stockée dans la mée-
moire. Ainsi, le système adopté calcule la fonction d'autocorrelation pour
une certaine valeur de la variable SSde déphasage dans un cycle d'échan-
tillonnage, stocke dans la mémoire uniquement le résultat de ce calcul et
compare ensuite ce résultat avec le résultat du calcul de la fonction d'au-
tocorrélation pour une variable de déphasage dont la valeur est avancée d'un
nombre dans le cycle suivant d'échantillonnage. En conséquence, seul le ré-
sultat du calcul de la fonction d'autocorrélation dans le cycle le plus ré-
cent d'échantillonnage doit être stocke dans la mémoire. Le système qui fait
l'objet de la présente invention permet donc de réduire la capacité né-
cessaire de la mémoire de corrélation, comparativement au système classi-
que qui oblige la mémoire de corrélation à stocker les résultats de cha-
que calcul de fonction d'autocorrélation englobant l'ensemble des données qui couvrent l'intervalle dans laquelle la variable 'Y de déphasage varie dans
chaque cycle d'échantillonnage.
Afin de détecter les pointes de signal conformément à la présente invention, la valeur qui a été calculée et stockée précédemment et qui se rapporte au cycle précédent d'échantillonnage, est comparée à la valeur calculée pour le cycle suivant d'échantillonnage. On détecte ensuite les pointes de signal en répétant ce processus de comparaison et en examinant
les changements d'état. Lorsqu'on observe un changement d'état d'une va-
leur élevée à une valeur faible entre deux cycles continus d'échantillon-
nage, ceci indique la détection d'une pointe dans le premier des deux cy-
cles. En effectuant l'opération de détection de pointe, on fait la compa-
raison uniquement avec la valeur calculée qui précède immédiatement, con-
formément à la description donnée ci-dessus. Cependant, on peut évidemment
stocker des valeurs calculées qui concernent plusieurs cycles, et faire
une comparaison de ces valeurs entre elles si on le désire.
Dans l'exemple de réalisation décrit ci-dessus, on peut utiliser un microprocesseur du fait de la réduction de la capacité nécessaire de stockage et de la réduction du nombre de calculs. On peut donc effectuer
des calculs précis de fonctions d'autocorrélation et la commande du sys-
tème. Il faut bien noter cependant que le fonctionnement ci-dessus détecte malheureusement non seulement une pointe intrinsèque correspondant à la période du signal, mais aussi d'autres pointes qui tendent généralement à se trouver au voisinage de la pointe intrinsèque. Par conséquent, afin de mesurer la période avec un ordre élevé de précision, on doit prévoir un moyen pour détecter la pointe intrinsèque ou réelle, qui correspond à la
période du signal, parmi les diverses pointes qui peuvent exister.
Afin de déterminer si une pointe détectée a la possibilité d'être
une pointe réelle, deux opérations sont nécessaires. Tout d'abord, une opé-
ration de contrôle du niveauest.effectuée sur la base d'un niveau minimal déterminé pour servir de valeur de seuil et, ensuite, lorsqu'une pointe a
été détectée, le calcul de la fonction d'autocorrélation est poursuivi pen-
dant un intervalle de temps qui correspond à la plus courte période de me-
sure, pour établir qu'aucune pointe supérieure à la pointe détectée n'exis-
te dans l'intervalle pendant lequel le calcul a été poursuivi. Ces deux
opérations permettent de détecter une pointe réelle.
L'opération de contrôle du niveau comprend les stades consis-
tant à déterminer la valeur de seuil d'un niveau utilisé en estimant si une pointe a la possibilité d'être une pointe réelle, et ensuite à apprécier si le niveau d'une pointe est supérieur à cette valeur de seuil, à la suite de quoi il est décidé si la pointe détectée, qui a la possibilité d'être une pointe réelle, doit être considérée en fait comme en étant une. Dans le présent exemple de réalisation, la valeur de seuil est fixée à la moitié de la valeur d'une pointe utilisée dans la mesure qui précède immédiatement, c'est-à-dire à la moitié de la valeur de la pointe réelle la plus récente, et seule la pointe dont le niveau est supérieur à
la valeur fixée de seuil, est considérée comme une pointe qui a la possi-
bilité d'être une pointe réelle.
La valeur de seuil ne doit pas nécessairement être fixée à la moitié de la valeur de la pointe réelle la plus récente, mais elle doit
être fixée à la valeur optimale choisie conformément à la qualité du si-
gnal reçu à ce moment-là. En général, bien que la valeur de la pointe réelle qui indique la période du signal soit influencée par l'intensité
et la forme d'onde du signal, le bruit constitue une difficulté parti-
culière. En particulier, plus le bruit sera faible, plus les pointes réel-
les se présenteront de manière distincte et apparaîtront grandes, tandis
que plus le bruit sera fort, plus les pointes réelles apparaîtront petites.
En fait, la valeur d'une pointe réelle en présence d'un bruit important peut même être inférieure à celle d'une fausse pointe au voisinage d'une
pointe réelle lorsqu'il y a un bruit faible.
C'est la raison pour laquelle la valeur de seuil doit être fixée
conformément à la qualité du signal reçu pendant la détection d'une pointe.
Dans le présent exemple de réalisation, en plus du contrôle du niveau dé-
crit ci-dessus, le calcul de la fonction d'autocorrélation est poursuivi pendant un intervalle fixe de temps faisant suite à la détection d'une pointe, et un contrôle est effectué pour déterminer si une pointe plus
grande que celle qui a été détectée, se trouve dans cet intervalle fixe.
On a indiqué ci-dessus que les pointes obtenues à partir d'une fonction d'autocorrélation comprennent, en plus d'une pointe réelle qui correspond à la période du signal, plusieurs pointes situées au voisinage de la pointe réelle. Cette dernière doit être détectée parmi les diverses pointes afin de permettre une mesure correcte de la période. Puisque les pointes au voisinage de la pointe réelle sont généralement situées très prés de celle-ci, on peut empêcher que ces pointes soient détectées comme
étant des pointesyéellesenpoursuivant le calcul de la fonction d'autocor-
rélation pendant un intervalle fixe de temps faisant suite à la détection d'une pointe, et en vérifiant ensuite si une pointe supérieure à la pointe détectée existe dans cet intervalle fixe. On doit bien noter qu'il suffit que l'intervalle fixe soit fixé à une valeur qui correspond à la période minimale de mesure. En conséquence, dans le présent exemple de réalisation, dès qu'une pointe a été détectée, le calcul de la fonction d'autocorrélation est poursuivi pendant un intervalle de temps qui correspond sensiblement à
la valeur minimale de la période de mesure, c'est-à-dire 300 microsecondes.
On se reportera maintenant à la figure 5. Si on suppose que la
pointe P1 est détectée au temps t1l (temps présent), le calcul de la fonc-
tion d'autocorrélation sera poursuivi pendant 300 millisecondes après le temps t1l, c'est-à-dire jusqu'au temps tl2. Comme le représente la figure , une pointe P. supérieure à la pointe PI est détectée au temps t21 dans l'intervalle de 300 millisecondes entre le temps t1l et le temps 12 Dans
ce cas, la pointe P1 est écartée et le calcul de la fonction d'autocorréla-
tion est poursuivi pendant un autre intervalle de 300 millisecondes com-
mençant à partir de la nouvelle pointe P2, c'est-à-dire jusqu'au temps t22.
La pointe P2 détectée est considérée comme étant la pointe réelle si l'exis-
tence d'aucune pointe supérieure à P2 n'est décelée dans ce dernier inter-
valle de 300 millisecondes. On doit bine noter que, à la figure 5, une est éetectee pointe P3 d'amplitude inférieure à la pointe P2 /à un certain temps t31 compris dans l'intervalle de 300 millisecondes entre le temps t12 auquel la pointe P2 est détectée, et le temps t22. Cependant, la pointe P dont l'amplitude-est inférieure à celle de la pointe P2, n'est pas considérée comme une pointe ayant la possibilité d'être une pointe réelle. Ainsi, la pointe P2 détectée au temps t21, est considérée comme étant une pointe
réelle indiquant la période lorsque 300 millisecondes se sont écoulées de-
puis le temps t21, c'est-à-dire lorsque le temps t22 a été atteint. A ce point du temps, le calcul de la fonction d'autocorrélation prend fin et lapériode est calculée. La valeur de la variable lkfde déphasage de la pointe réelle trouvée de cette manière correspond à la période. En admettant que TS soit la période d'échantillonnage des données, la période T est obtenue
par la formule de calcul T = t x Ts. La nouvelle mesure de période com-
mence de nouveau à partir de 'k = 60 (correspondant à la période de 300
millisecondes) et elle se poursuit de la même manière.
Ainsi,la période correcte du biosignal est mesurée de la manière
décrite ci-dessus.
Dans ce qui précède, le fait que la fonction d'autocorrélation commence à partir de 300 millisecondes sur l'axe >' d'autocorrélation et se
termine en un point équivalent à la période du biosignal T + 300 millisecon-
des, est extrêmement important en ce qui concerne la détection des pointes
réelles et le temps de sortie des résultats de la mesure.
Tout d'abord, en ce qui concerne la détection des pointes réelles, une pointe réelle ne peut exister en dessous de la période la plus courte possible du biosignal faisant l'objet de la mesure, et une pointe réelle ne
peut exister non plus dans un intervalle situé dans la période la plus cour-
te. Par conséquent, on peut dire que les pointes qui sont confirmées de cet-
te manière sont des pointes qui n'ont absolument aucune possibilité de re-
présenter des pointes dont la période est égale à deux fois la période réel-
le. En ce qui concerne le temps de sortie des résultats de mesure, la disposition mentionnée ci-dessus a pour effet de permettre la sortie des
* résultats de mesure en synchronisation avec la période réelle du biosignal.
Plus précisément, la mesure de la période commence à partir de 300 milli-
secondes qui constituent la période la plus courte possible. D'autre part,
le temps de 300 millisecondes, équivalent à la période la plus courte pos-
sible, est fixé comme étant l'intervalle de confirmation d'une pointe réel-
le, de sorte que les résultats de la mesure peuvent être obtenus dans un
intervalle de temps qui est équivalent à la période réelle du biosignal.
Par exemple, si la période réelle est de 500 millisecondes, la sortie des résultats de la mesure aura lieu toutes les 500 millisecondes. Lorsque la période varie, les intervalles de sortie varient pareillement. Ceci est dû au fait que le calcul de la fonction d'autocorrélation se poursuit en temps
réel sur l'axe de corrélation si l'intervalle du calcul de la fonction d'au-
tocorrélation coîncide avec la période d'échantillonnage des données, c'est-
à-dire que le calcul de corrélation, pendant une durée de temps allant de la
période la plus courte du biosignal jusqu'à un temps représenté par la som-
me de la période la plus courte et de la période réelle, est effectué dans
un temps équivalent à la période réelle du biosignal.
La figure 6 représente, sous une forme simplifiée, la construction d'un appareil de mesure de période pour la mise en pratique du système de
mesure de période décrit ci-dessus en rapport avec les figures 3 à 5.
Comme le représente la figure 6, un transducteur est mis en contact avec l'abdomen W d'un sujet féminin pour détecter le signal de pulsation du
coeur du foetus. Un circuit 24 d'échantillonnage est connecté au transduc-
teur 22 par l'intermédiaire d'un circuit 23 de pré-traitement. Le signal de
pulsation du coeur est détecté par le transducteur 22 et, après que sa for-
me d'ondea été conformée de manière appropriée parle circuit 23 de pré-
traitement, ce signal est échantillonné par le circuit 24 d'échantillon-
nage suivant une période prédéterminée d'échantillonnage et il est sou-
mis à une conversion d'analogique en numérique par le circuit d'échantil-
lonnage. Le signal de pulsation du coeur sort donc du circuit 24 d'échantillonnage sous forme de signal numérique. Une mémoire 26 de données est
reliée au circuit 24 d'échantillonnage et elle stocke les données échan-
tillonnées provenant de ce circuit. La mémoire 26 de données se compose
de plusieurs registres à décalage et elle fonctionne de la manière suivan-
te. Au fur et à mesure que chaque élément nouveau d'information est intro-
duit dans la mémoire de données, les éléments d'information déjà stockés
jusqu' alors sont décalés par octet, les élémentsles plus anciens d'in-
formation étant perdus à chaque nouvelle entrée de données. Un multiplica-
teur 28 est relié à la mémoire 26 de données, et un additionneur est relié
au multiplicateur 28. Plus précisément, la mémoire 26 de données ou le re-
gistre à décalage comprend un registre parallèle d'un octet (8 bits) qui est prévu pour introduire par décalage les données échantillonnées sous
forme numérique. Il est construit de telle sorte que les données arbitrai-
res de position précisées par la ligne analogique-numérique du signal peu-
vent en être extraites. La mémoire 26 de données comprend une mémoire à
accès sélectif avec capacité de lecture et d'écriture, et un organe de com-
mande de la mémoire à accès sélectif.
Le multiplicateur 28 et un additionneur 30 constituent un circuit de calcul de la fonction d'autocorrélation. Ce circuit calcule la fonction d'autocorrélation d'un biosignal, c'est-à-dire le signal représentant la
pulsation du coeur d'un foetus, en effectuant le calcul indiqué essentielle-
ment par l'équation (3) à l'aide des données stockées dans la mémoire 26 de données. En d'autres termes, le calcul de la fonction d'autocorrélation est
effectué en liaison avec une variable Or de déphasage ayant une certaine va-
leur dans chaque cycle d'échantillonnage. Pour être plus précis, deux élé-
ments d'information qui représentent deux positions sur l'axe des temps sé-
parées l'une de l'autre par la variable %-de déphasage, sont produits par
un circuit 42 de commande de la manière décrite plus loin et les deux élé-
ments d'information sont stockés dans deux adresses de la partie mémoire
de la mémoire 26 de données (les adresses donnant les emplacements de mé-
moire qui sont indiqués par des hachures dans le bloc 26 de la figure 7).
Pour calculer la fonction d'autocorrélation, les deux éléments stockés d'in-
formation sont multipliés et le produit est implanté dans un accumulateur situé dans l'additionneur 30. Le nombre d'opérations de multiplication pour une variable % de déphasage est représenté par "n" dans l'équation (3),
comme on le comprendra facilement d'après la description qui précède, de
sorte que le nombre d'additions est également égal à "n". Les "n' addi-
tions terminées, on obtient en fait la variable n, de déphasage dont la valeur représente "n" fois la fonction de corrélation. Cependant, puisque
"n" est constant, les données calculées sont proportionnelles à la fonc-
tion d'autocorrélation dans l'équation (3) de sorte que, au fond, la fonc-
tion d'autocorrélation est calculée.
Un détecteur 32 de pointes,connecté à l'additionneur 30, permet de stocker une faible quantité de données et d'effectuer une opération de
comparaison. L'entrée du détecteur 32 de pointes est la valeur de la fonc-
tion d'autocorrélation calculée par le circuit de calcul composé du multi-
plicateur 28 et de l'additionneur 30. Comme on le décrira avec plus de dé-
tails ci-dessous, le détecteur 32 de pointes stocke la valeur calculée pré-
cédemment de la fonction d'autocorrélation pour un cycle d'échantillonnage, et il compare cette valeur avec la valeur calculée nouvellement introduite
de la fonction d'autocorrélation pour le cycle suivant d'échantillonnage.
Le détecteur de pointesstocke ensuite la valeur calculée nouvellement in-
troduite, si elle est supérieure à la valeur calculée stockée précédemment.
Du fait que le détecteur 32 de pointes ne doit stocker que la valeur calcu-
lée de la fonction d'autocorrélation pour le cycle d'échantillonnage le plus
récent et la valeur de la variable 1û de déphasage à ce moment-là, une fai-
ble capacité de mémoire suffira. Ainsi, la valeur calculée stockée pour un
cycle d'échantillonnage est comparée au moyen d'un comparateur avec la va-
leur calculée de la fonction d'autocorrélation pour le cycle suivant d'échan-
tillonnage, ce qui permet de rechercher la variation des valeurs pour les
deux cycles d'échantillonnage. Lorsque le résultat de l'opération de com-
paraison fait apparaître un passage d'une valeur élevée à une valeur basse, cela indique l'existence d'une pointe dans le premier des deux cycles d'échantillonnage. Le détecteur 32 depointes établit une comparaison entre un signal de détection de pointe et un niveau de référence. Pour fixer le
niveau de référence, on peut utiliser un niveau qui représente, par exem-
ple, la moitié de la valeur de pointe réelle précédemment mesurée, comme on l'a décrit plus haut. Si la pointe détectée est supérieure au niveau de référence et s'il se confirme qu'aucune pointe supérieure à la pointe détectée n'est présente dans un intervalle fixe de temps mesuré à partir de l'instant o la pointe détectée est supérieure au niveau de référence (cet intervalle fixe de temps est de 300 millisecondes dans le présent exemple de réalisation), le détecteur 32 de pointes considère alors que la pointe détectée est une pointe réelle et il émet un signal de détection
de pointe réelle.
Au détecteur 32 de pointes est connecté un circuit 38 de calcul
de la période qui, lorsqu'il reçoit le signal de détection de pointe réel-
le venant du détecteur 32 de pointes, calcule la période en se basant sur la valeur de la variable de déphasage dans la fonction d'autocorrélation au moment o la pointe est obtenue, cette valeur étant conservée dans un
registre situé dans le détecteur de pointes.
Au circuit 38 de calcul de la période est connecté un circuit 40 de calcul des pulsations du coeur qui calcule le nombre de pulsations du coeur en se basant sur la période calculée par le circuit 38 de calcul de
la période.
Le circuit 40 de calcul des pulsations du coeur est connecté à un circuit 42 de commande auquel est relié un dispositif 44 d'affichage, tel qu'un agencement de diodes électroluminescentes. Ce dispositif 44 d'affichage affiche le nombre de pulsations du coeur contenues dans le signal de pulsation du coeur en se basant sur le signal obtenu depuis le
circuit 40 de calcul des pulsations du coeur par l'intermédiaire du cir-
cuit 42 de commande. Il peut arriver, dans certains cas, que le signal ve-
nant du circuit 40 de calcul des pulsations du coeur comprenne une com-
posante de bruit, ou que la sonde de détection des pulsations du coeur
glisse. Le circuit 42 de commande est prévu, par conséquent, pour conmman-
der le signal venant du circuit 40 de calcul des pulsations du coeur,de ma-
nière à l'empêcher d'entrer dans le dispositif 44 d'affichage dans les cas mentionnés ci-dessus, afin d'assurer qu'un nombre erroné de pulsations du
coeur ne sera pas affiché.
Le circuit 42 de commande est prévu en outre pour fournir des im-
pulsions d'horloge au circuit 24 d'échantillonnage, afin de commander la synchronisation de l'opération d'échantillonnage effectuée par le circuit d'échantillonnage. En outre, lors de chaque opération d'échantillonnage, le circuit de commande envoie au multiplicateur 28 un signal qui représente
la valeur de la variable de déphasage. La valeur de la variable de déphasa-
ge avance au fur et à mesure que progressent les cycles d'échantillonnage, en partant d'un temps qui correspond sensiblement à la valeur minimale de la période du signal représentant la pulsation du coeur. Le multiplicateur 28 est prévu pour lire, à partir de la mémoire 26 de données, deux éléments d'information séparés par la valeur de la variable de déphasage désignée
par le signal venant du circuit 42 de commande, et pour calculer le pro-
duit de ces deux éléments d'information. Le circuit 42 de commande envoie un signal de synchronisation à l'additionneur 30 qui, en se basant sur ce signal de synchronisation, additionne ensemble les résultats des calculs effectués par le multiplicateur 28. En d'autres termes, le multiplicateur 28 et l'additionneur 30, régis par le circuit 42 de commande, lisent les données provenant de la mémoire de données et calculent la fonction d'auto-
corrélation sensiblement comme l'indique l'équation (3).
Un détecteur 46 de niveau de référence est connecté au circuit 42
de commande. Conformément à un signal de synchronisation fourni par le cir-
cuit de commande à un intervalle approprié de temps, le détecteur 46 de ni-
veau de référence est prévu pour détecter le niveau optimum de référence (niveau zéro) dans le but d'appliquer un signal positif (+) ou un signal négatif (-) aux données échantillonnées, et pour envoyer au circuit 24
d'échantillonnage un signal représentant le niveau optimum de référence.
Dans l'application des signes aux données, plus la polarité des données sera équilibrée, plus la périodicité de la fonction d'autocorrélation sera faible. Le détecteur 46 de niveau de référence est prévu pour découvrir la
valeur optimale afin de réaliser cet objectif. Plus précisément, le détec-
teur 46 découvre la valeur optimale du niveau de référence en détectant la valeur maximale et la valeur minimale, ou la valeur moyenne, des données
pendant l'échantillonnage.
Le détecteur 32 de pointes peut être construit comme le représente la figure 8. Dans ce cas, une mémoire 52 comprend deux unités de mémoire, une pour stocker la valeur de la fonction de corrélation et l'autre pour stocker la valeur de la variable de déphasage. Plus précisément, la mémoire 52, commandée par un signal d'écriture venant d'un comparateur 54, stocke la valeur de la fonction d'autocorrélation calculée par l'additionneur 30, et la valeur de la variable de déphasage obtenue à partir du circuit 42 de commande. Le comparateur 54 est prévu pour comparer la valeur nouvellement calculée de la fonction-d'autocorrélation obtenue à partir de l'additionneur 30 avec la valeur calculée la plus élévée et la plus récente de la fonction d'autocorrélation stockée précédemment dans la mémoire 52, et pour envoyer le signal d'écriture à la mémoire 52 si la valeur nouvellement calculée de
la fonction d'autocorrélation est la plus élevée des deux valeurs; en con-
séquence, les contenus de la mémoire 52 sont remplacés respectivement par la valeur nouvellement calculée de la fonction d'autocorrélation et par la
valeur de la variable de déphasage obtenue à partir du circuit 42 de com-
mande. Lorsque la valeur de la fonction d'autocorrélation cesse d'être croissante pour devenir décroissante lors de la répétition de l'opération de comparaison mentionnée ci-dessus, le comparateur 54 considère qu'une pointe a été détectée et il émet donc un signal. La valeur calculée de la fonction d'autocorrélation implantée dans la mémoire 52 est envoyée
à un comparateur 56 en vue du contrôle du niveau de pointe. Le com-
parateur 56 compare cette valeur avec un niveau de référence reçu d'un générateur 58 de niveau de référence. Ce dernier est réglé par la synchronisation de sortie d'un compteur 62 avrec le temps auquel la pointe réelle précédente est détectée et, de ce fait, il stocke un niveau égal, par exemple, à la moitié de la valeur de la pointe réelle détectée par la
mesure précédente. C'est ce niveau que le générateur de niveau de ré-
férence fournit comme niveau de référence. La moitié de la valeur d'une pointe réelle est obtenue par la technique représentée à la figure 9. Plus précisément, on obtient ce résultat en décalant les données de sortie de la mémoire 52 d'un bit vers le côté du bit de poids faible et en connectant les données au comparateur 56 qui est un comparateur de grandeur. Si le résultat de la comparaison est tel que la valeur calculée de la fonction d'autocorrélation stockée dans la mémoire 52 est d'un niveau supérieur au niveau de référence, le comparateur 56 émet un signal. Une porte ET reçoit le produit logique des sorties des comparateurs 54 et 56. Une transition dans le sens positif à la sortie de la porte ET 60 remet à zéro le compteur 62 et fixe la valeur de la variable t de déphasage, qui a été stockée dans la mémoire 52, dans un registre 64. Lorsque les impulsions
d'horloge comptées par le compteur 62 atteignent un nombre qui corres-
pond à une période fixe de temps, telle que 300 millisecondes, le comp-
teur émet un signal. Le signal de sortie du compteur 62 indique qu'une pointe réelle a été détectée, de sorte que la valeur de ? qui a été fixée dans le registre 64, est envoyée au circuit 38 de calcul de la période. Ce dernier circuit calcule la période en prenant le produit de la variable -V et de la période d'échantillonnage venant du circuit 42 de commande par
la ligne SP' des signaux. A titre d'exemple, si la période d'échantillon-
nage est de cinq millisecondes et si te est égal à 60 millisecondes, la
période est calculée comme étant de 300 millisecondes. La période obte-
nue est envoyée au circuit 40 de calcul des pulsations du coeur o le nom-
bre de pulsations du coeur pendant une période d'une minute est obtenu en divisant 60 x 103 (ms) par la période (ms). Le nombre de pulsations du coeur trouvé de cette manière est envoyé ensuite au circuit 42 de commande et affiché sur le dispositif 44 d'affichage régi par le circuit
de commande.
Ainsi, les pointes sont détectées et contrôlées grâce au dispo-
sitif et aux opérations décrites ci-dessus afin d'assurer l'extraction de
pointes qui sont réelles.
Suivant la présente invention décrite ci-dessus, on mesure la période d'un biosignal en effectuant les opérations suivantes: calcul de la fonction d'autocorrélation pour une certaine valeur de la variable à', de déphasage dans un cycle d'échantillonnage du biosignal; changement
de la valeur de la variable -t de déphasage sur l'axe des temps confor-
mément à la progression des cycles d'échantillonnage; calcul de la fonc-
tion d'autocorrélation dans chaque cycle d'échantillonnage; stockage uni-
quement du résultat du calcul de la fonction d'autocorrélation pour le cycle initial de deux cycles consécutifs d'échantillonnage; comparaison de ce résultat avec le résultat du calcul de la fonction d'autocorrélation pour le cycle suivant; et détection d'une pointe à partir de la croissance et de la décroissance du résultat de la comparaison, ce qui permet de mesurer la période du biosignal. Cette disposition permet de réduire
fortement la capacité de stockage des résultats des calculs de la fonc-
tion d'autocorrélation, et d'éliminer les calculs dénués de sens de la fonction d'autocorrélation pendant de longs intervalles de temps dont la
longueur peut représenter deux ou trois fois la période réelle du biosi-
gnal, ce qui permet de traiter les données approximativement en temps
réel.
En outre, selon une autre caractéristique de la présente inven-
tion, on peut mesurer la période correcte d'un biosignal en effectuant
les opérations suivantes commencement du calcul de la fonction d'auto-
corrélation sensiblement à partir de la valeur minimale de la mesure
de la période du biosignal; continuation du calcul de la fonction d'auto-
corrélation dans un intervalle correspondant à cette valeur minimale faisant suite à la détection d'une pointe; et confirmation qu'il n'existe
pas de pointe supérieure à la pointe initiale dans cet intervalle corres-
pondant à la valeur minimale mesuré à partir du point de détection de la pointe initiale, afin de déterminer que la pointe initiale est une pointe
réelle. Il est donc possible de détecter fiablement uniquement une poin-
te réelle qui indique la période intrinsèque du biosignal, ce qui permet donc de mesurer la période correcte. En outre, du fait que la portée du calcul de la fonction d'autocorrélation est limitée à une zone comprise sensiblement entre la valeur minimale mentionnée ci-dessus et une gamme de valeurs représentées par la somme de la période réelle du biosignal et de l'intervalle de confirmation (tel que la valeur minimale),
la présente invention a pour effet d'éliminer les calculs sans significa-
tion et de permettre le traitement en temps réel. De plus, les résultats
des mesures peuvent être fournis à un intervalle de temps qui est équi-
valent à la période du signal qui fait l'objet de la mesure.
La présente invention n'est pas limitée aux exemples de réali-
sation qui viennent d'être décrits, elle est au contraire susceptible de
variantes et de modifications qui apparaîtront à l'homme de l'art.
REV ENDICATIONS
1. Système de mesure de période, caractérisé en ce qu'il com-
prend un moyen d'échantillonnage pour échantillonner un biosignal à une
période prédéterminée d'échantillonnage; un moyen de calcul de la fonc-
tion d'autocorrélation pour calculer la fonction d'autocorrélation du biosignal en se basant sur le biosignal fourni par le moyen d'échantillonnage; un moyen de détection de pointe pour détecter une pointe à partir de la fonction d'autocorrélation calculée; et un moyen de calcul de période pour calculer la période du biosignal à partir du point sur l'axe des
temps de corrélation auquel une pointe est détectée par le moyen de dé-
tection de pointe; la fonction d'autocorrélation, donnée par l'équation n A (1) =- 'Z f (k). f (k +t) pour une certaine valeur d'une variable n k=l
qui applique un déphasage au biosignal sur l'axe des temps, étant calcu-
lée par le moyen de calcul de la fonction d'autocorrélation pour une va-
leur spécifique de la variable F de déphasage; cette valeur spécifique de la variable t de déphasage étant avancée sur l'axe des temps pour se conformer à la progression des cycles d'échantillonnage, ce qui permet
au moyen de calcul de la fonction d'autocorrélation de calculer des fonc-
tions d'autocorrélation les unes après les autres correspondant aux nou-
velles valeurs spécifiques de la variable de déphasage; la valeur calcu-
lée de la fonction d'autocorrélation étant stockée en mémoire et compa-
rée ensuite à la valeur calculée la plus récente de la fonction d'autocor-
rélation, ce qui permet de détecter une pointe.
2. Système de mesure de période suivant la revendication 1,
caractérisé en ce que la valeur de la variable t qui applique le déphasa-
ge est choisie de manière à commencer le calcul de la fonction d'auto-
corrélation sensiblement à partir de la valeur minimale de la période
du biosignal qui fait l'objet de la mesure.
3. Système de mesure de période suivant la revendication 1,
caractérisé en ce que le calcul de la fonction d'autocorrélation se pour-
suit pendant une période fixe de temps qui fait suite à la détection d'une certaine pointe; et en ce qu'il est confirmé qu'aucune pointe supérieure à cette pointe n'existe dans cette période fixe de temps, ce qui permet de
déterminer que cette pointe est une pointe réelle.
4. Système de mesure de période suivant la revendication 3, caractérisé en ce que cette période fixe de temps est sensiblement égale
à la valeur minimale de la.période du biosignal qui fait l'objet de la me-
sure. 5. Système de mesure de période suivant la revendication 4,
caractérisé en ce que la valeur de la variable It qui applique les dépha-
sages commence sensiblement à la valeur minimale ' de la période de
mesure de la période du biosignal; et en ce que,lorsqu'une pointe est dé-
tectée, le calcul de la fonction d'autocorrélation se poursuit pendant un
intervalle de temps t qui correspond à la période minimale de la me-
sure afin de confirmer que la pointe est une pointe réelle, ce qui permet
de sortir le résultat de la mesure à un intervalle de temps qui est équi-
valent à la période du biosignal.
6. Système de mesure de période, caractérisé en ce qu'il com-
prend un moyen pour extraire un biosignal; un moyen de calcul de la fonction d'autocorrélation pour calculer la fonction d'autocorrélation du biosignal; un moyen de détection pour détecter une pointe à partir de la
fonction d'autocorrélation; et un moyen de calcul de période pour calcu-
ler la période du biosignal à partir du point sur l'axe de corrélation au-
quel une pointe est détectée par le moyen de détection de pointe; le cal-
cul de la fonction d'autocorrélation se poursuivant dans un intervalle
correspondant sensiblement à la valeur minimale de la période de me-
sure, cet intervalle commençant avec la détection d'une pointe; le fait étant confirmé qu'aucune pointe supérieure à la pointe détectée n'existe
dans cet intervalle qui correspond à la valeur minimale et qui commen-
ce avec la détection de la pointe, ce qui permet de déterminer que la
pointe détectée est une pointe réelle.
7. Système de mesure de période suivant la revendication 6, caractérisé en ce que la valeur de seuil du niveau de détection de pointes du moyen de détection de pointes est fixée à la valeur d'une pointe réelle détectée servant de référence, le moyen de détection de pointes détectant comme étant des pointes les seules pointes qui dépassent la valeur de seuil.
FR8113883A 1980-07-17 1981-07-16 Systeme de mesure de periode Granted FR2495330A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP9782180A JPS5722737A (en) 1980-07-17 1980-07-17 Cycle measuring system
JP9782080A JPS5722736A (en) 1980-07-17 1980-07-17 Cycle measuring system

Publications (2)

Publication Number Publication Date
FR2495330A1 true FR2495330A1 (fr) 1982-06-04
FR2495330B1 FR2495330B1 (fr) 1985-01-25

Family

ID=26438963

Family Applications (1)

Application Number Title Priority Date Filing Date
FR8113883A Granted FR2495330A1 (fr) 1980-07-17 1981-07-16 Systeme de mesure de periode

Country Status (4)

Country Link
US (1) US4463425A (fr)
CA (1) CA1174733A (fr)
DE (1) DE3128171A1 (fr)
FR (1) FR2495330A1 (fr)

Families Citing this family (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4562846A (en) * 1983-09-15 1986-01-07 Duke University System and process for monitoring myocardial integrity
US4732157A (en) * 1986-08-18 1988-03-22 Massachusetts Institute Of Technology Method and apparatus for quantifying beat-to-beat variability in physiologic waveforms
EP0613652B1 (fr) * 1990-02-15 1997-04-16 Hewlett-Packard GmbH Appareil et méthode pour la mesure non-invasive du taux d'oxygène
US5241473A (en) * 1990-10-12 1993-08-31 Ken Ishihara Ultrasonic diagnostic apparatus for displaying motion of moving portion by superposing a plurality of differential images
US5558096A (en) * 1995-07-21 1996-09-24 Biochem International, Inc. Blood pulse detection method using autocorrelation
WO1997005821A1 (fr) * 1995-08-10 1997-02-20 Pentavox Kft. Procede et appareil pour mesurer le rythme cardiaque d'un foetus et sonde electroacoustique pour capter les sons provenant du coeur d'un foetus
US6115624A (en) * 1997-07-30 2000-09-05 Genesis Technologies, Inc. Multiparameter fetal monitoring device
US6937696B1 (en) * 1998-10-23 2005-08-30 Varian Medical Systems Technologies, Inc. Method and system for predictive physiological gating
US6279579B1 (en) 1998-10-23 2001-08-28 Varian Medical Systems, Inc. Method and system for positioning patients for medical treatment procedures
US6973202B2 (en) * 1998-10-23 2005-12-06 Varian Medical Systems Technologies, Inc. Single-camera tracking of an object
US7158610B2 (en) * 2003-09-05 2007-01-02 Varian Medical Systems Technologies, Inc. Systems and methods for processing x-ray images
US6980679B2 (en) * 1998-10-23 2005-12-27 Varian Medical System Technologies, Inc. Method and system for monitoring breathing activity of a subject
US6621889B1 (en) * 1998-10-23 2003-09-16 Varian Medical Systems, Inc. Method and system for predictive physiological gating of radiation therapy
DE69916871T2 (de) 1998-10-23 2005-03-31 Varian Medical Systems Technologies, Inc., Palo Alto Verfahren und system zur physiologischen steuerung von radiotherapie
US8788020B2 (en) 1998-10-23 2014-07-22 Varian Medical Systems, Inc. Method and system for radiation application
US6477476B1 (en) * 1999-12-06 2002-11-05 Koninklijke Philips Electronics N.V. Periodic-signal analysis via correlation
US7769430B2 (en) * 2001-06-26 2010-08-03 Varian Medical Systems, Inc. Patient visual instruction techniques for synchronizing breathing with a medical procedure
ES2232223B1 (es) * 2002-03-26 2006-02-16 Osatu, S. Coop. Metodo para la determinacion de la frecuencia de forma de onda de una señal ecg.
US7254502B2 (en) * 2002-08-26 2007-08-07 Infineon Technologies Ag Method and device for detecting period length fluctuations of periodic signals
US20050053267A1 (en) * 2003-09-05 2005-03-10 Varian Medical Systems Technologies, Inc. Systems and methods for tracking moving targets and monitoring object positions
US8571639B2 (en) * 2003-09-05 2013-10-29 Varian Medical Systems, Inc. Systems and methods for gating medical procedures
US20060074305A1 (en) * 2004-09-30 2006-04-06 Varian Medical Systems Technologies, Inc. Patient multimedia display
US9119541B2 (en) * 2005-08-30 2015-09-01 Varian Medical Systems, Inc. Eyewear for patient prompting
US20070167847A1 (en) * 2006-01-19 2007-07-19 Guglielmino Michael F Method and device for using a physiological parameter to express evolution
US7953247B2 (en) * 2007-05-21 2011-05-31 Snap-On Incorporated Method and apparatus for wheel alignment
JP5453861B2 (ja) * 2008-03-31 2014-03-26 Jfeスチール株式会社 周期性欠陥検出装置及びその方法
US10667727B2 (en) * 2008-09-05 2020-06-02 Varian Medical Systems, Inc. Systems and methods for determining a state of a patient
US20100061596A1 (en) * 2008-09-05 2010-03-11 Varian Medical Systems Technologies, Inc. Video-Based Breathing Monitoring Without Fiducial Tracking
JP5516428B2 (ja) 2010-10-14 2014-06-11 株式会社村田製作所 拍動周期算出装置およびこれを備えた生体センサ
JP6129082B2 (ja) * 2011-03-11 2017-05-17 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 生理学的信号を監視する監視装置
KR102655669B1 (ko) * 2016-07-20 2024-04-05 삼성전자주식회사 생체신호의 특징 추출 장치 및 방법과, 생체정보 검출 장치
CN112869724B (zh) * 2021-01-19 2022-04-22 西安交通大学 一种基于多通道被动式采集信号的胎儿健康监测仪

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS567592B2 (fr) * 1974-10-31 1981-02-18
US4054862A (en) * 1975-10-28 1977-10-18 Raytheon Company Ranging system with resolution of correlator ambiguities
US4037151A (en) * 1976-02-27 1977-07-19 Hewlett-Packard Company Apparatus for measuring the period and frequency of a signal
DE2818768C2 (de) * 1978-04-28 1986-07-24 Hewlett-Packard GmbH, 7030 Böblingen Verfahren und Vorrichtung zum Messen der Frequenz bzw. der Periodendauer der Grundschwingung eines annähernd periodischen Eingangssignales
US4239048A (en) * 1979-02-06 1980-12-16 Multitronics Corporation Cardiotachometer using autocorrelation techniques

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
IEE PROCEEDINGS A, vol. 128, no. 8, novembre 1971, HITCHIN HERTS (GB) *
IEEE TRANSACTIONS ON BIO-MEDICAL ENGINEERING, vol. BME-13, no. 1, janvier 1966, NEW YORK (US) *
IEEE TRANSACTIONS ON BIO-MEDICAL ENGINEERING, vol.BME-15, no. 1, janvier 1968, NEW YORK (US) *

Also Published As

Publication number Publication date
US4463425A (en) 1984-07-31
CA1174733A (fr) 1984-09-18
FR2495330B1 (fr) 1985-01-25
DE3128171A1 (de) 1982-04-22

Similar Documents

Publication Publication Date Title
FR2495330A1 (fr) Systeme de mesure de periode
FR2487075A1 (fr) Systeme de mesure de periode
US5008902A (en) Automatic baud rate detection
WO2018010107A1 (fr) Procédé et dispositif pour déterminer si une onde pulsatile est valide
FR2626120A1 (fr) Generateur d&#39;impulsions en correlation
FR2562675A1 (fr) Procede de levee d&#39;ambiguite de la mesure par effet doppler de la vitesse d&#39;un mobile
WO2005107584A1 (fr) Procede d’obtention de parametres hemodynamiques derives d’un signal de pression arterielle
EP0052041B1 (fr) Procédé de détection de la fréquence de mélodie dans un signal de parole, et dispositif destiné à la mise en oeuvre de ce procédé
FR2559908A1 (fr) Circuit numerique de calcul d&#39;une frequence moyenne
WO2015110398A1 (fr) Procédé et système d&#39;extraction dynamique d&#39;impulsions dans un signal temporel bruité
EP0224315A1 (fr) Discriminateur numérique de fréquence, et application à un signal composite VOR
JPS624971B2 (fr)
EP0200596B1 (fr) Système de suppression du bruit et de ses variations pour la détection d&#39;un signal pur dans signal discret bruité mesuré
US20080284409A1 (en) Signal Recognition Method With a Low-Cost Microcontroller
EP0309320A1 (fr) Système de suppression du bruit et de ses variations pour la détection d&#39;un signal pur dans un signal discret bruité mesuré
FR2550083A1 (fr) Appareil et procede pour mesurer la tension arterielle suivant une methode indirecte
EP1580893B1 (fr) Dispositif et procédé de détection de phase d&#39;un signal
FR2498071A1 (fr) Procede de mesure de la pression arterielle et dispositif pour mettre en oeuvre ce procede
CN108761441A (zh) 一种雷达信号相位编码规律识别方法及装置
EP0727893B1 (fr) Méthode et dispositif pour régénérer, par détection des maxima, un signal d&#39;horloge rythmant la transmission de signaux numérisés reçus
FR2475323A1 (fr) Procede et dispositif pour la determination d&#39;un signal utile
FR2983055A1 (fr) Detection et estimation du complexe qrs pour le suivi d&#39;une activite cardiaque et pulmonaire
EP0566928A2 (fr) Procédé de détection de tons pour appareil téléphonique
JP2556369B2 (ja) アナログ波形のオーバーシュート位置検出方式
JPH021566A (ja) 周期測定装置