CA2797886C - Procede et systeme d'analyse de l'activite respiratoire d'un patient et applications correspondantes - Google Patents
Procede et systeme d'analyse de l'activite respiratoire d'un patient et applications correspondantes Download PDFInfo
- Publication number
- CA2797886C CA2797886C CA2797886A CA2797886A CA2797886C CA 2797886 C CA2797886 C CA 2797886C CA 2797886 A CA2797886 A CA 2797886A CA 2797886 A CA2797886 A CA 2797886A CA 2797886 C CA2797886 C CA 2797886C
- Authority
- CA
- Canada
- Prior art keywords
- signal
- phi
- elementary
- respiratory
- phase
- 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.)
- Active
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/087—Measuring breath flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7271—Specific aspects of physiological measurement analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7271—Specific aspects of physiological measurement analysis
- A61B5/7275—Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7271—Specific aspects of physiological measurement analysis
- A61B5/7278—Artificial waveform generation or derivation, e.g. synthesising signals from measured signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
- A61M16/00—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
- A61M16/021—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes operated by electrical means
- A61M16/022—Control means therefor
- A61M16/024—Control means therefor including calculation means, e.g. using a processor
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
- A61M16/00—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
- A61M16/06—Respiratory or anaesthetic masks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Hematology (AREA)
- Emergency Medicine (AREA)
- Anesthesiology (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Discrete Mathematics (AREA)
- Power Engineering (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
Ce procédé d'analyse de l'activité respiratoire d'un patient comprend les étapes d'acquisition (20) d'au moins un signal d'activité respiratoire comprenant au moins un signal élémentaire correspondant à un cycle respiratoire, dont la forme générale peut être exprimée par x(t) = x 0 + x 1 cos(F(t)), où F(t) est la phase dudit signal élémentaire, et d'analyse (24, 26, 28) dudit signal d'activité respiratoire. Cette analyse (24, 26, 28) comprend les étapes d'extraction (24), à partir dudit signal d'activité respiratoire, dudit signal élémentaire, de détermination (26) d'une expression d'une équation de phase (I) dudit signal élémentaire et de détermination (26) d'une expression de la phase F(t) dudit signal élémentaire en fonction de paramètres mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcosn et psinn définies par : (II)
Description
Procédé et système d'analyse de l'activité respiratoire d'un patient et applications correspondantes La présente invention concerne un procédé et un système d'analyse de l'activité
respiratoire d'un patient, le procédé comprenant les étapes d'acquisition d'au moins un signal d'activité respiratoire comprenant au moins un signal élémentaire correspondant à
un cycle respiratoire, dont la forme générale peut être exprimée par x(t)= x0 + x1 cos(c13(t)), où ci(t) est la phase dudit signal élémentaire, et d'analyse dudit signal d'activité respiratoire.
Elle concerne également des applications de ceux-ci à la commande d'un dispositif d'assistance respiratoire et à un dispositif de surveillance respiratoire.
Elle s'applique en particulier à la détection de troubles respiratoires ou à
la commande de dispositifs de ventilation artificielle.
Un signal d'activité respiratoire est un signal de mesure de la variation d'une grandeur liée à l'activité respiratoire du patient, telle que le débit et la pression d'air ou la concentration en oxygène et en dioxyde de carbone à l'entrée de ses voies respiratoires, ou la concentration en oxygène dans le sang. Ces grandeurs peuvent être mesurées par des appareils de mesure non invasifs, par exemple un capteur de débit ou de pression intégré à un masque placé devant la bouche du patient ou un oxymètre, ou par des capteurs internes, par exemple des capteurs de pression placés dans le circuit respiratoire du patient. De tels signaux peuvent également être déduits de signaux d'électrocardiogramme.
L'activité respiratoire est constituée d'une succession de cycles respiratoires, comprenant une phase inspiratoire et une phase expiratoire, à une fréquence appelée fréquence respiratoire. Par conséquent, les signaux d'activité respiratoire sont des signaux quasi-périodiques, comprenant une succession de signaux élémentaires, chacun de ces signaux élémentaires étant caractéristique d'un cycle respiratoire.
L'analyse de ces signaux permet de détecter des troubles ou anomalies respiratoires, tels que l'apnée du sommeil ou l'asthme. Cependant, cette analyse est généralement limitée à la détermination de la fréquence respiratoire et de sa variabilité et de l'amplitude de ces signaux, et aucune analyse de la forme d'onde de ces signaux n'est réalisée.
Or, la forme d'onde de signaux d'activité respiratoire est caractéristique de cette activité respiratoire, et leur analyse peut permettre de détecter efficacement d'éventuelles anomalies respiratoires.
On connaît de nombreuses méthodes d'analyse et de caractérisation d'un signal périodique. En particulier, l'analyse fréquentielle d'un signal permet de décrire ce signal
respiratoire d'un patient, le procédé comprenant les étapes d'acquisition d'au moins un signal d'activité respiratoire comprenant au moins un signal élémentaire correspondant à
un cycle respiratoire, dont la forme générale peut être exprimée par x(t)= x0 + x1 cos(c13(t)), où ci(t) est la phase dudit signal élémentaire, et d'analyse dudit signal d'activité respiratoire.
Elle concerne également des applications de ceux-ci à la commande d'un dispositif d'assistance respiratoire et à un dispositif de surveillance respiratoire.
Elle s'applique en particulier à la détection de troubles respiratoires ou à
la commande de dispositifs de ventilation artificielle.
Un signal d'activité respiratoire est un signal de mesure de la variation d'une grandeur liée à l'activité respiratoire du patient, telle que le débit et la pression d'air ou la concentration en oxygène et en dioxyde de carbone à l'entrée de ses voies respiratoires, ou la concentration en oxygène dans le sang. Ces grandeurs peuvent être mesurées par des appareils de mesure non invasifs, par exemple un capteur de débit ou de pression intégré à un masque placé devant la bouche du patient ou un oxymètre, ou par des capteurs internes, par exemple des capteurs de pression placés dans le circuit respiratoire du patient. De tels signaux peuvent également être déduits de signaux d'électrocardiogramme.
L'activité respiratoire est constituée d'une succession de cycles respiratoires, comprenant une phase inspiratoire et une phase expiratoire, à une fréquence appelée fréquence respiratoire. Par conséquent, les signaux d'activité respiratoire sont des signaux quasi-périodiques, comprenant une succession de signaux élémentaires, chacun de ces signaux élémentaires étant caractéristique d'un cycle respiratoire.
L'analyse de ces signaux permet de détecter des troubles ou anomalies respiratoires, tels que l'apnée du sommeil ou l'asthme. Cependant, cette analyse est généralement limitée à la détermination de la fréquence respiratoire et de sa variabilité et de l'amplitude de ces signaux, et aucune analyse de la forme d'onde de ces signaux n'est réalisée.
Or, la forme d'onde de signaux d'activité respiratoire est caractéristique de cette activité respiratoire, et leur analyse peut permettre de détecter efficacement d'éventuelles anomalies respiratoires.
On connaît de nombreuses méthodes d'analyse et de caractérisation d'un signal périodique. En particulier, l'analyse fréquentielle d'un signal permet de décrire ce signal
2 dans l'espace de Fourier. La décomposition de Fourier consiste en effet à
décomposer un signal périodique de fréquence f en une somme infinie de fonctions sinusoïdales de fréquences multiples de f, pondérées par les coefficients de Fourier. Ces coefficients de Fourier, qui constituent un codage du signal analysé, sont des paramètres caractéristiques de ce signal. En pratique, le nombre de coefficients de Fourier conservés est limité, et seuls les premiers termes de la décomposition de Fourier sont gardés. Ces termes doivent cependant être en nombre suffisant pour caractériser efficacement le signal.
Or, les signaux d'activité respiratoire sont des signaux anharmoniques, c'est-à-dire non-linéaires, et la décomposition de Fourier de tels signaux nécessite de conserver un grand nombre de coefficients, coefficients auxquels il est difficile de donner un sens physique. La décomposition de Fourier est donc inadaptée pour l'analyse de ces signaux.
L'invention a donc pour but de permettre l'analyse des formes d'onde de signaux d'activité respiratoire au moyen d'un petit nombre de paramètres porteurs d'un sens physique et constituant une signature simple et explicite de la forme de ces signaux.
A cet effet, l'invention a pour objet un procédé d'analyse du type précité, caractérisé en ce que l'analyse du signal d'activité respiratoire comprend les étapes suivantes :
- extraction, à partir dudit signal d'activité respiratoire, dudit signal élémentaire;
ck13 - détermination d'une expression d'une équation de phase F(c13) = ¨dt dudit signal élémentaire ; et - détermination d'une expression de la phase c13(t) dudit signal élémentaire en fonction de paramètres (r, rk, 430, pk) mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcos, et psinr, définies paroc :
pcos, (t, r)= Ecos(e)¨rk k et psin, (t, r)= sin(kt)¨r .
k=1 kn k=1 k"
Le procédé selon l'invention comporte également les caractéristiques suivantes, prises séparément ou en combinaison :
- l'équation de phase est exprimée sous la forme :
dc13 1+ r2 + 2r cos() dt 1¨r2 dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité
dudit signal élémentaire ;
- le signal élémentaire est exprimé au moyen de deux paramètres r et (130, sous la forme :
décomposer un signal périodique de fréquence f en une somme infinie de fonctions sinusoïdales de fréquences multiples de f, pondérées par les coefficients de Fourier. Ces coefficients de Fourier, qui constituent un codage du signal analysé, sont des paramètres caractéristiques de ce signal. En pratique, le nombre de coefficients de Fourier conservés est limité, et seuls les premiers termes de la décomposition de Fourier sont gardés. Ces termes doivent cependant être en nombre suffisant pour caractériser efficacement le signal.
Or, les signaux d'activité respiratoire sont des signaux anharmoniques, c'est-à-dire non-linéaires, et la décomposition de Fourier de tels signaux nécessite de conserver un grand nombre de coefficients, coefficients auxquels il est difficile de donner un sens physique. La décomposition de Fourier est donc inadaptée pour l'analyse de ces signaux.
L'invention a donc pour but de permettre l'analyse des formes d'onde de signaux d'activité respiratoire au moyen d'un petit nombre de paramètres porteurs d'un sens physique et constituant une signature simple et explicite de la forme de ces signaux.
A cet effet, l'invention a pour objet un procédé d'analyse du type précité, caractérisé en ce que l'analyse du signal d'activité respiratoire comprend les étapes suivantes :
- extraction, à partir dudit signal d'activité respiratoire, dudit signal élémentaire;
ck13 - détermination d'une expression d'une équation de phase F(c13) = ¨dt dudit signal élémentaire ; et - détermination d'une expression de la phase c13(t) dudit signal élémentaire en fonction de paramètres (r, rk, 430, pk) mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcos, et psinr, définies paroc :
pcos, (t, r)= Ecos(e)¨rk k et psin, (t, r)= sin(kt)¨r .
k=1 kn k=1 k"
Le procédé selon l'invention comporte également les caractéristiques suivantes, prises séparément ou en combinaison :
- l'équation de phase est exprimée sous la forme :
dc13 1+ r2 + 2r cos() dt 1¨r2 dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité
dudit signal élémentaire ;
- le signal élémentaire est exprimé au moyen de deux paramètres r et (130, sous la forme :
3 x(t)= x, a,h sin(t,r)+ cos(t,r) où a, = x, cos(c1)0) et b, = -x sin(d) 0) , les fonctions hsin et hcos étant définies par:
\ (1+ r2)cos(t)- 2r (1- r2)sin(t) h cos :(t,r --> et h sin : (t, ->
1 + r2 - 2r cosU 1+r2 - 2r cos(t) - l'équation de phase est exprimée sous la forme :
P4) Fc1))=
Q((1)) dans laquelle P(0) et Q((I)) sont des polynômes trigonométriques ;
- l'expression de la phase 1100 est déterminée sous la forme :
t (C) = '51) akp (cb -Pk,rk )- bk P cosi (cp -Pk ,rk) k dans laquelle les fonctions psini et pcosi sont définies par:
r p cosi(tr) cos((t)¨ rk et p sin, (t, r) = sinkt)¨
k =1 k=1 Ainsi réalisé, le procédé selon l'invention permet d'analyser les signaux d'activité
respiratoire et de caractériser ces signaux au moyen d'un petit nombre de paramètres, par rapport aux procédés d'analyse de signaux périodiques selon l'état de la technique.
De plus, ces paramètres ont un sens physique, et sont caractéristiques des formes d'onde de ces signaux.
Selon un autre aspect, l'invention a également pour objet un système d'analyse de l'activité respiratoire d'un patient comprenant des moyens pour acquérir un signal d'activité respiratoire comprenant au moins un signal élémentaire correspondant à un cycle respiratoire dont la forme générale peut être exprimée par x(t)- x, + xl cos(a)(0), où
(DO est la phase dudit signal élémentaire, et des moyens pour analyser ledit signal d'activité respiratoire, caractérisé en ce que les moyens pour analyser ledit signal d'activité respiratoire comprennent :
- des moyens pour extraire, à partir dudit signal d'activité respiratoire, ledit signal élémentaire ;
dict.
- des moyens pour déterminer une expression d'une équation de phase F()=
di dudit signal élémentaire ; et - des moyens pour déterminer une expression de la phase D(t) dudit signal élémentaire en fonction de paramètres (r, rk, (1)0, pk) mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcosn et psin, définies par:
\ (1+ r2)cos(t)- 2r (1- r2)sin(t) h cos :(t,r --> et h sin : (t, ->
1 + r2 - 2r cosU 1+r2 - 2r cos(t) - l'équation de phase est exprimée sous la forme :
P4) Fc1))=
Q((1)) dans laquelle P(0) et Q((I)) sont des polynômes trigonométriques ;
- l'expression de la phase 1100 est déterminée sous la forme :
t (C) = '51) akp (cb -Pk,rk )- bk P cosi (cp -Pk ,rk) k dans laquelle les fonctions psini et pcosi sont définies par:
r p cosi(tr) cos((t)¨ rk et p sin, (t, r) = sinkt)¨
k =1 k=1 Ainsi réalisé, le procédé selon l'invention permet d'analyser les signaux d'activité
respiratoire et de caractériser ces signaux au moyen d'un petit nombre de paramètres, par rapport aux procédés d'analyse de signaux périodiques selon l'état de la technique.
De plus, ces paramètres ont un sens physique, et sont caractéristiques des formes d'onde de ces signaux.
Selon un autre aspect, l'invention a également pour objet un système d'analyse de l'activité respiratoire d'un patient comprenant des moyens pour acquérir un signal d'activité respiratoire comprenant au moins un signal élémentaire correspondant à un cycle respiratoire dont la forme générale peut être exprimée par x(t)- x, + xl cos(a)(0), où
(DO est la phase dudit signal élémentaire, et des moyens pour analyser ledit signal d'activité respiratoire, caractérisé en ce que les moyens pour analyser ledit signal d'activité respiratoire comprennent :
- des moyens pour extraire, à partir dudit signal d'activité respiratoire, ledit signal élémentaire ;
dict.
- des moyens pour déterminer une expression d'une équation de phase F()=
di dudit signal élémentaire ; et - des moyens pour déterminer une expression de la phase D(t) dudit signal élémentaire en fonction de paramètres (r, rk, (1)0, pk) mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcosn et psin, définies par:
4 rk rk p cos, (t, = )¨ et p sin, (t, r) E sin ---k k kn Le système selon l'invention comporte également les caractéristiques suivantes, prises séparément ou en combinaison :
- le système comporte des moyens pour exprimer l'équation de phase sous la forme :
el) 1 -L r 2 + 2r cos (cl)) dt 1-r2 dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité
dudit signal élémentaire ;
- le système comporte des moyens pour exprimer ledit signal élémentaire au moyen de deux paramètres r et (Do, sous la forme :
x(t)-= x0 + ai h sin(t, r)+ cos(t, r) où al =x1 cos(c1)0) et b1 = -xl sin(11)0), les fonctions hsin et hcos étant définies par:
\ + r )cos(t) - 2r , - T2 )sin(t) h cos : t, r õ et h sin : r) 1 + r2 -2r cos(t ) 1 + r 2 - 2r cos(t) - le système comporte des moyens) pour exprimer l'équation de phase sous la forme Q() dans laquelle P(1:1)) et Q(cI)) sont des polynômes trigonométriques ;
- le système comporte des moyens pour exprimer la phase (1)(t) sous la forme :
01= (D+ ak p sin, (.1) - pk ,rk)- b k p cosi(a) - pk , rk) k=1 dans laquelle les fonctions psini et pcosi sont définies par:
rk p cos (t, r) cos(kt)¨ et p sin, (t,T)-sin(kt)¨rk k =1 k=1 Selon d'autres aspects, l'invention a également pour objet un dispositif d'assistance respiratoire et un système d'analyse de l'activité respiratoire comprenant un système d'analyse de l'activité respiratoire selon l'invention.
L'invention sera mieux comprise au regard d'un exemple de réalisation de l'invention qui va maintenant être décrit en faisant référence aux figures annexées parmi lesquelles :
- la figure 1 représente de manière schématique un système d'analyse selon un mode de réalisation de l'invention ;
- la figure 2 représente un signal d'activité respiratoire ; et - la figure 3 est un schéma synoptique illustrant le procédé selon un mode de
- le système comporte des moyens pour exprimer l'équation de phase sous la forme :
el) 1 -L r 2 + 2r cos (cl)) dt 1-r2 dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité
dudit signal élémentaire ;
- le système comporte des moyens pour exprimer ledit signal élémentaire au moyen de deux paramètres r et (Do, sous la forme :
x(t)-= x0 + ai h sin(t, r)+ cos(t, r) où al =x1 cos(c1)0) et b1 = -xl sin(11)0), les fonctions hsin et hcos étant définies par:
\ + r )cos(t) - 2r , - T2 )sin(t) h cos : t, r õ et h sin : r) 1 + r2 -2r cos(t ) 1 + r 2 - 2r cos(t) - le système comporte des moyens) pour exprimer l'équation de phase sous la forme Q() dans laquelle P(1:1)) et Q(cI)) sont des polynômes trigonométriques ;
- le système comporte des moyens pour exprimer la phase (1)(t) sous la forme :
01= (D+ ak p sin, (.1) - pk ,rk)- b k p cosi(a) - pk , rk) k=1 dans laquelle les fonctions psini et pcosi sont définies par:
rk p cos (t, r) cos(kt)¨ et p sin, (t,T)-sin(kt)¨rk k =1 k=1 Selon d'autres aspects, l'invention a également pour objet un dispositif d'assistance respiratoire et un système d'analyse de l'activité respiratoire comprenant un système d'analyse de l'activité respiratoire selon l'invention.
L'invention sera mieux comprise au regard d'un exemple de réalisation de l'invention qui va maintenant être décrit en faisant référence aux figures annexées parmi lesquelles :
- la figure 1 représente de manière schématique un système d'analyse selon un mode de réalisation de l'invention ;
- la figure 2 représente un signal d'activité respiratoire ; et - la figure 3 est un schéma synoptique illustrant le procédé selon un mode de
5 réalisation de l'invention.
On a représenté sur la figure 1, un système d'acquisition et d'analyse de signaux d'activité respiratoire.
Ce système comprend des moyens 1 d'acquisition d'un signal d'activité
respiratoire, un boîtier 3 d'acquisition, relié aux moyens 1 d'acquisition, et des moyens 5 d'analyse d'un signal d'activité respiratoire, par exemple un processeur, reliés au boîtier 3 d'acquisition.
Les moyens 1 d'acquisition sont aptes à recueillir un signal d'activité
respiratoire. A
titre d'exemple, ces moyens 1 comprennent un pneumotachographe 7, placé à la sortie d'un masque respiratoire 9 couvrant le nez et la bouche d'un patient 11. Le pneumotachographe 7 est apte à mesurer de manière continue le débit d'air inspiré et expiré par le patient 11 et à transmettre au boîtier 3 d'acquisition un signal électrique analogique, caractéristique de ce débit.
Le boîtier 3 d'acquisition comprend un convertisseur analogique/numérique, apte à
convertir un signal analogique en un signal numérique, appelé par la suite signal d'activité
respiratoire, par échantillonnage et quantification du signal analogique.
Le processeur 5 est apte à analyser un signal d'activité respiratoire numérique de manière à en extraire des paramètres caractéristiques de la forme de ce signal.
On a représenté sur la figure 2 un tracé illustrant la forme d'un signal 13 d'activité
respiratoire tel que mesuré par le pneumotachographe 7. Sur ce tracé, le temps est représenté en abscisse, et le débit d'air, depuis l'extérieur vers l'appareil respiratoire du patient 11, en ordonnée. On reconnaît sur ce tracé quatre signaux élémentaires 15, correspondant chacun à un cycle respiratoire, comprenant une phase inspiratoire 16, durant laquelle le débit est positif, et une phase expiratoire 17, durant laquelle le débit est négatif.
La figure 3 est un schéma synoptique illustrant l'acquisition et l'analyse d'un signal d'activité respiratoire au moyen du système décrit en référence à la figure 1, selon un mode de réalisation de l'invention.
Dans une première étape 20 d'acquisition, le débit d'air inspiré et expiré par le patient 11 est mesuré par le pneumotachographe 7, qui transmet de manière continue au boîtier 3 d'acquisition un signal électrique analogique, dont l'amplitude instantanée est proportionnelle au débit mesuré.
On a représenté sur la figure 1, un système d'acquisition et d'analyse de signaux d'activité respiratoire.
Ce système comprend des moyens 1 d'acquisition d'un signal d'activité
respiratoire, un boîtier 3 d'acquisition, relié aux moyens 1 d'acquisition, et des moyens 5 d'analyse d'un signal d'activité respiratoire, par exemple un processeur, reliés au boîtier 3 d'acquisition.
Les moyens 1 d'acquisition sont aptes à recueillir un signal d'activité
respiratoire. A
titre d'exemple, ces moyens 1 comprennent un pneumotachographe 7, placé à la sortie d'un masque respiratoire 9 couvrant le nez et la bouche d'un patient 11. Le pneumotachographe 7 est apte à mesurer de manière continue le débit d'air inspiré et expiré par le patient 11 et à transmettre au boîtier 3 d'acquisition un signal électrique analogique, caractéristique de ce débit.
Le boîtier 3 d'acquisition comprend un convertisseur analogique/numérique, apte à
convertir un signal analogique en un signal numérique, appelé par la suite signal d'activité
respiratoire, par échantillonnage et quantification du signal analogique.
Le processeur 5 est apte à analyser un signal d'activité respiratoire numérique de manière à en extraire des paramètres caractéristiques de la forme de ce signal.
On a représenté sur la figure 2 un tracé illustrant la forme d'un signal 13 d'activité
respiratoire tel que mesuré par le pneumotachographe 7. Sur ce tracé, le temps est représenté en abscisse, et le débit d'air, depuis l'extérieur vers l'appareil respiratoire du patient 11, en ordonnée. On reconnaît sur ce tracé quatre signaux élémentaires 15, correspondant chacun à un cycle respiratoire, comprenant une phase inspiratoire 16, durant laquelle le débit est positif, et une phase expiratoire 17, durant laquelle le débit est négatif.
La figure 3 est un schéma synoptique illustrant l'acquisition et l'analyse d'un signal d'activité respiratoire au moyen du système décrit en référence à la figure 1, selon un mode de réalisation de l'invention.
Dans une première étape 20 d'acquisition, le débit d'air inspiré et expiré par le patient 11 est mesuré par le pneumotachographe 7, qui transmet de manière continue au boîtier 3 d'acquisition un signal électrique analogique, dont l'amplitude instantanée est proportionnelle au débit mesuré.
6 Ce signal analogique est numérisé dans une étape 22 par le convertisseur analogique/numérique du boîtier 3 d'acquisition, par échantillonnage et quantification, et le signal ainsi numérisé, ou signal d'activité respiratoire, est transmis au processeur 5.
Ce signal 13 d'activité respiratoire est composé d'une succession de signaux élémentaires 15, correspondant chacun à un cycle respiratoire. Il n'est cependant pas rigoureusement périodique, en raison notamment de la variabilité de la fréquence respiratoire et du débit d'air inspiré et expiré pendant chacun des cycles. On définit donc la période respiratoire instantanée comme l'intervalle de temps entre le début d'un cycle respiratoire, c'est-à-dire le début de sa phase d'inspiration, et le début du cycle respiratoire précédent.
Dans une étape 24, le processeur décompose le signal 13 d'activité
respiratoire en signaux élémentaires 15. Cette décomposition peut par exemple être effectuée en détectant les instants auxquels le signal passe d'une valeur négative à une valeur positive, ce qui correspond à une transition entre la phase expiratoire d'un cycle et la phase inspiratoire du cycle suivant. Ainsi, lors de l'étape 24, le processeur détermine la fréquence respiratoire instantanée du patient et sa variabilité.
Puis, le processeur 5 analyse dans une étape 26 chacun des signaux élémentaires 15 issus de cette décomposition.
Chaque signal élémentaire x(t) est un signal anharmonique, qui peut être décrit sous la forme suivante :
x(t) = xo + X, cos(cP(t)) (1) dans laquelle toute la dépendance temporelle est contenue dans la fonction de phase (1).
Ce signal élémentaire x(t) est considéré comme un signal périodique de période T, T étant égal à l'inverse de la fréquence respiratoire instantanée.
Or, dans un signal anharmonique, la principale contribution à l'anharmonicité
provient de la brisure de symétrie de la dynamique de phase. Ainsi, toute l'information dynamique pertinente est exprimée par la dynamique de phase. Lors de l'analyse du signal x(t), il convient donc d'étudier cette phase 43(t), et en particulier la dynamique de phase exprimée par la fonction F, dérivée de la fonction (13 par rapport au temps t:
F(43) = cl(1dt 3 (2) Ainsi, la morphologie du signal x(t) est complètement déterminée par la connaissance de F.
L'étape d'analyse 26 du procédé selon l'invention consiste donc à décrire cette fonction F au moyen d'un très petit nombre de paramètres. On entendra par petit nombre
Ce signal 13 d'activité respiratoire est composé d'une succession de signaux élémentaires 15, correspondant chacun à un cycle respiratoire. Il n'est cependant pas rigoureusement périodique, en raison notamment de la variabilité de la fréquence respiratoire et du débit d'air inspiré et expiré pendant chacun des cycles. On définit donc la période respiratoire instantanée comme l'intervalle de temps entre le début d'un cycle respiratoire, c'est-à-dire le début de sa phase d'inspiration, et le début du cycle respiratoire précédent.
Dans une étape 24, le processeur décompose le signal 13 d'activité
respiratoire en signaux élémentaires 15. Cette décomposition peut par exemple être effectuée en détectant les instants auxquels le signal passe d'une valeur négative à une valeur positive, ce qui correspond à une transition entre la phase expiratoire d'un cycle et la phase inspiratoire du cycle suivant. Ainsi, lors de l'étape 24, le processeur détermine la fréquence respiratoire instantanée du patient et sa variabilité.
Puis, le processeur 5 analyse dans une étape 26 chacun des signaux élémentaires 15 issus de cette décomposition.
Chaque signal élémentaire x(t) est un signal anharmonique, qui peut être décrit sous la forme suivante :
x(t) = xo + X, cos(cP(t)) (1) dans laquelle toute la dépendance temporelle est contenue dans la fonction de phase (1).
Ce signal élémentaire x(t) est considéré comme un signal périodique de période T, T étant égal à l'inverse de la fréquence respiratoire instantanée.
Or, dans un signal anharmonique, la principale contribution à l'anharmonicité
provient de la brisure de symétrie de la dynamique de phase. Ainsi, toute l'information dynamique pertinente est exprimée par la dynamique de phase. Lors de l'analyse du signal x(t), il convient donc d'étudier cette phase 43(t), et en particulier la dynamique de phase exprimée par la fonction F, dérivée de la fonction (13 par rapport au temps t:
F(43) = cl(1dt 3 (2) Ainsi, la morphologie du signal x(t) est complètement déterminée par la connaissance de F.
L'étape d'analyse 26 du procédé selon l'invention consiste donc à décrire cette fonction F au moyen d'un très petit nombre de paramètres. On entendra par petit nombre
7 de paramètres un nombre de paramètres réduit par rapport au nombre de paramètres nécessaire à la décomposition de la même fonction, au moyen des séries de Fourier, avec un niveau de précision équivalent.
Cette étape d'analyse 26 comprend ainsi une première étape consistant à
exprimer la phase cl), et en particulier la fonction F, dérivée de (I) par rapport au temps.
Dans le cas le plus simple, et pour un signal de période 27c, la dynamique de phase peut être écrite sous la forme :
F(1)) = doI) = 1+r2 + 2r cos(d)) (3) dt 1- r2 appelée équation de phase.
La fonction F présente dans ce cas une symétrie de réflexion par rapport à
l'axe GD=0. Cette expression de la dynamique de phase ne contient qu'un seul paramètre, r, qui varie dans l'intervalle [0,1]. La limite r=0 correspond à un signal harmonique, la limite r=1 à un signal infiniment anharmonique.
Le signal x(t), qui peut être écrit sous la forme :
x(t)-x0 + x1 cos(d)(t, r) - (Do ) (4) où (1)0 est une origine de phase, est décomposé et réécrit sous une forme faisant intervenir les paramètres r et (Do :
x(t) = x0 + aih sin(t, r)+ b cos(t, r) (5) avec ai = x1 cos(d0) et b, = 1 sin(d),:,), et dans laquelle on a défini les fonctions hcos et hsin suivantes :
(1+ r2)cos(t)- 2r h cos : r) --->(6) 1+r - 2r cos(t) (1- r2 )sin(t) hsin: r)-> (7) 1+ r 2 - 2r cos(t) Ainsi, la décomposition du signal x(t) fait intervenir seulement deux paramètres, r et 1.
r, appelé paramètre d'anharmonicité, mesure le degré d'anharmonicité du signal, la limite r=0 correspond à un signal harmonique, la limite r=1 à un signal infiniment anharmonique. Par ailleurs, le paramètre (1)2, qui définit la composition du signal dans les deux fonctions hcos et hsin, est un paramètre de morphologie, qui correspond à
l'angle de symétrie de réflexion de la dynamique de phase.
Dans le cas général, c'est-à-dire pour un signal périodique quelconque, l'équation de phase peut s'écrire sous la forme :
Cette étape d'analyse 26 comprend ainsi une première étape consistant à
exprimer la phase cl), et en particulier la fonction F, dérivée de (I) par rapport au temps.
Dans le cas le plus simple, et pour un signal de période 27c, la dynamique de phase peut être écrite sous la forme :
F(1)) = doI) = 1+r2 + 2r cos(d)) (3) dt 1- r2 appelée équation de phase.
La fonction F présente dans ce cas une symétrie de réflexion par rapport à
l'axe GD=0. Cette expression de la dynamique de phase ne contient qu'un seul paramètre, r, qui varie dans l'intervalle [0,1]. La limite r=0 correspond à un signal harmonique, la limite r=1 à un signal infiniment anharmonique.
Le signal x(t), qui peut être écrit sous la forme :
x(t)-x0 + x1 cos(d)(t, r) - (Do ) (4) où (1)0 est une origine de phase, est décomposé et réécrit sous une forme faisant intervenir les paramètres r et (Do :
x(t) = x0 + aih sin(t, r)+ b cos(t, r) (5) avec ai = x1 cos(d0) et b, = 1 sin(d),:,), et dans laquelle on a défini les fonctions hcos et hsin suivantes :
(1+ r2)cos(t)- 2r h cos : r) --->(6) 1+r - 2r cos(t) (1- r2 )sin(t) hsin: r)-> (7) 1+ r 2 - 2r cos(t) Ainsi, la décomposition du signal x(t) fait intervenir seulement deux paramètres, r et 1.
r, appelé paramètre d'anharmonicité, mesure le degré d'anharmonicité du signal, la limite r=0 correspond à un signal harmonique, la limite r=1 à un signal infiniment anharmonique. Par ailleurs, le paramètre (1)2, qui définit la composition du signal dans les deux fonctions hcos et hsin, est un paramètre de morphologie, qui correspond à
l'angle de symétrie de réflexion de la dynamique de phase.
Dans le cas général, c'est-à-dire pour un signal périodique quelconque, l'équation de phase peut s'écrire sous la forme :
8 F()= Pr, (4.) Q ni (43) (8) dans laquelle Pr, et Qm sont des polynômes trigonométriques de degrés respectifs n et m.
La forme générale d'un polynôme trigonométrique de degré n est:
13, (cl)) = a 0 En ak cos,k(13, + bk sin(kcP) (9) k=1 L'analyse du signal x(t) comprend alors la détermination d'une expression de c13 faisant intervenir un petit nombre de paramètres, ce qui permet de déterminer une expression du signal x(t) en fonction de ces paramètres.
Avantageusement, l'équation de phase (2) peut être réécrite sous la forme :
1 dt Qm (cD) F(cI)) P n (cI)) (10) La factorisation du polynôme F',.,() permet de transformer F(1,) en une somme de termes simples, ce qui permet de réécrire l'équation de phase sous la forme :
d t ak cos(c13 ¨ pk)+ bk sin(c13 + pk) (11) cel) k=1 (1+ rk2 ¨ 2rk cos(c1) + pk )) dans laquelle les paramètres rk, compris entre 0 et 1, mesurent l'anharmonicité du signal x(t), et les paramètres pk caractérisent sa morphologie.
La période T du signal peut être déterminée en intégrant cette équation par rapport à (13, entre 0 et 2t:
<I>=27c ce ( r ak T= F(E13.) = 2A- a + Ek 1 ¨k __________ (12) ,D=o A partir de ce résultat, et des contraintes selon lesquelles la période est égale à 27c et le signal est harmonique lorsque les coefficients rk sont tous nuls, l'équation de phase peut être exprimée ainsi :
d t =1+EDk(cP¨pk) (13) k=1 Où la fonction Dk est définie par:
rk(ak cos(*) + bk sin()¨ ak) Dk : C13 (14) (1+ ¨ 2rk cos(<13)) Et vérifie :
<I>=2z f D k(C0C1c13 = 0 (15) d¨o
La forme générale d'un polynôme trigonométrique de degré n est:
13, (cl)) = a 0 En ak cos,k(13, + bk sin(kcP) (9) k=1 L'analyse du signal x(t) comprend alors la détermination d'une expression de c13 faisant intervenir un petit nombre de paramètres, ce qui permet de déterminer une expression du signal x(t) en fonction de ces paramètres.
Avantageusement, l'équation de phase (2) peut être réécrite sous la forme :
1 dt Qm (cD) F(cI)) P n (cI)) (10) La factorisation du polynôme F',.,() permet de transformer F(1,) en une somme de termes simples, ce qui permet de réécrire l'équation de phase sous la forme :
d t ak cos(c13 ¨ pk)+ bk sin(c13 + pk) (11) cel) k=1 (1+ rk2 ¨ 2rk cos(c1) + pk )) dans laquelle les paramètres rk, compris entre 0 et 1, mesurent l'anharmonicité du signal x(t), et les paramètres pk caractérisent sa morphologie.
La période T du signal peut être déterminée en intégrant cette équation par rapport à (13, entre 0 et 2t:
<I>=27c ce ( r ak T= F(E13.) = 2A- a + Ek 1 ¨k __________ (12) ,D=o A partir de ce résultat, et des contraintes selon lesquelles la période est égale à 27c et le signal est harmonique lorsque les coefficients rk sont tous nuls, l'équation de phase peut être exprimée ainsi :
d t =1+EDk(cP¨pk) (13) k=1 Où la fonction Dk est définie par:
rk(ak cos(*) + bk sin()¨ ak) Dk : C13 (14) (1+ ¨ 2rk cos(<13)) Et vérifie :
<I>=2z f D k(C0C1c13 = 0 (15) d¨o
9 La définition des fonctions des fonctions polycos et polysin, notées pcosn et psinn, qui s'expriment par:
rk pcosn(t,r)=cos(kt)¨ (16) k=1 kn psinn(t,r), sin(kt)4- (17) k=1 kn et possèdent entre autre les propriétés suivantes :
r(cos(t)¨r) pcoso(t, r) = (18) 1+r 2 ¨ 2r cos(t) r sin(t) psino(t,r)= (19) 1+r 2 ¨ 2r cos(t) pcosl (t, = --1In(1+ r2 ¨ 2r cos(t)) (20) psini(t,r), tan'r rsin(t) (21) r cos(t)) permet de réécrire l'équation de phase sous la forme :
¨dt=1 + En akpcoso (c13¨ pk , rk )+ bk psin 0 (.13 ¨ pk , rk ) (22) k=1 La résolution de cette équation permet d'accéder à une expression analytique de t(<13), qui s'exprime par:
t(c1)), + En akpsinl(d)¨ pk,rk )¨ bkpcosi(c1D¨ pk, rk) (23) k=1 Le temps t est donc exprimé en fonction de la phase <13, et de manière duale la phase <D est exprimée en fonction du temps t, à l'aide de paramètres indépendants clairement définis, qui mesurent l'anharmonicité (paramètres r ou rk), et la morphologie (paramètres c1)0 ou pk).
Ainsi, lors de l'étape 26 d'analyse, le processeur 5 code chaque signal élémentaire x(t) au moyen d'un petit nombre de paramètres. Selon un mode de réalisation, chaque signal élémentaire x(t) est décrite de manière quasi-exacte par une amplitude, une harmonicité r et une morphologie .1)0. Selon un autre mode de réalisation, chaque signal élémentaire x(t) est décrit de manière encore plus précise par deux couples de paramètres (ri, Pi) et (r2, P2), complétés de leurs poids respectifs.
Chacun des signaux élémentaires, donc chaque cycle respiratoire, est donc caractérisé par un nombre restreint de paramètres, porteurs d'un sens physique car représentatifs de la non-linéarité et de la morphologie de ce signal.
Puis, dans une étape 28, le processeur compare les valeurs des paramètres déterminés pour chaque signal élémentaire, c'est-à-dire la fréquence respiratoire instantanée, l'amplitude et les paramètres de morphologie et d'harmonicité du signal, à
des valeurs tabulées de paramètres préalablement enregistrées, de manière à
détecter 5 d'éventuelles anomalies respiratoires.
Le procédé selon l'invention permet ainsi d'analyser l'activité respiratoire, et d'extraire d'un signal d'activité respiratoire un nombre restreint de paramètres, permettant une représentation compacte et pertinente de la forme d'onde de ce signal, et la détection précise d'éventuelles anomalies respiratoires.
rk pcosn(t,r)=cos(kt)¨ (16) k=1 kn psinn(t,r), sin(kt)4- (17) k=1 kn et possèdent entre autre les propriétés suivantes :
r(cos(t)¨r) pcoso(t, r) = (18) 1+r 2 ¨ 2r cos(t) r sin(t) psino(t,r)= (19) 1+r 2 ¨ 2r cos(t) pcosl (t, = --1In(1+ r2 ¨ 2r cos(t)) (20) psini(t,r), tan'r rsin(t) (21) r cos(t)) permet de réécrire l'équation de phase sous la forme :
¨dt=1 + En akpcoso (c13¨ pk , rk )+ bk psin 0 (.13 ¨ pk , rk ) (22) k=1 La résolution de cette équation permet d'accéder à une expression analytique de t(<13), qui s'exprime par:
t(c1)), + En akpsinl(d)¨ pk,rk )¨ bkpcosi(c1D¨ pk, rk) (23) k=1 Le temps t est donc exprimé en fonction de la phase <13, et de manière duale la phase <D est exprimée en fonction du temps t, à l'aide de paramètres indépendants clairement définis, qui mesurent l'anharmonicité (paramètres r ou rk), et la morphologie (paramètres c1)0 ou pk).
Ainsi, lors de l'étape 26 d'analyse, le processeur 5 code chaque signal élémentaire x(t) au moyen d'un petit nombre de paramètres. Selon un mode de réalisation, chaque signal élémentaire x(t) est décrite de manière quasi-exacte par une amplitude, une harmonicité r et une morphologie .1)0. Selon un autre mode de réalisation, chaque signal élémentaire x(t) est décrit de manière encore plus précise par deux couples de paramètres (ri, Pi) et (r2, P2), complétés de leurs poids respectifs.
Chacun des signaux élémentaires, donc chaque cycle respiratoire, est donc caractérisé par un nombre restreint de paramètres, porteurs d'un sens physique car représentatifs de la non-linéarité et de la morphologie de ce signal.
Puis, dans une étape 28, le processeur compare les valeurs des paramètres déterminés pour chaque signal élémentaire, c'est-à-dire la fréquence respiratoire instantanée, l'amplitude et les paramètres de morphologie et d'harmonicité du signal, à
des valeurs tabulées de paramètres préalablement enregistrées, de manière à
détecter 5 d'éventuelles anomalies respiratoires.
Le procédé selon l'invention permet ainsi d'analyser l'activité respiratoire, et d'extraire d'un signal d'activité respiratoire un nombre restreint de paramètres, permettant une représentation compacte et pertinente de la forme d'onde de ce signal, et la détection précise d'éventuelles anomalies respiratoires.
10 Les étapes 20, 22, 24, 26 et 28 peuvent être réalisées au fur et à
mesure de l'acquisition du signal d'activité respiratoire, de manière à surveiller en continu l'activité
respiratoire du patient.
Selon un autre mode de réalisation, le système et le procédé selon l'invention peuvent être mis en oeuvre dans un respirateur artificiel.
Un respirateur artificiel est un dispositif visant à pallier l'insuffisance respiratoire d'un patient, en insufflant de l'air vers les poumons du patient par l'intermédiaire d'un masque, au moyen d'un ventilateur. Pour que cette ventilation soit efficace, il est nécessaire que les efforts inspiratoires du patient et le déclenchement du ventilateur soient synchronisés, c'est-à-dire que l'air soit insufflé dans les poumons par le ventilateur seulement pendant la phase inspiratoire du cycle respiratoire du patient.
La mesure du débit d'air inspiré et expiré par le patient permet de détecter le début et la fin de cette phase inspiratoire, et ainsi de synchroniser le fonctionnement du ventilateur sur le cycle respiratoire du patient.
Cette détection est généralement réalisée en fixant deux seuils de débit ou de pression, le ventilateur étant actionné dès que le débit ou la pression d'air inspiré par le patient dépasse un premier seuil, et arrêté dès que ce débit ou cette pression devient inférieur à un deuxième seuil.
Cette méthode de détection peut s'avérer inefficace, en raison notamment de la variabilité des cycles respiratoires, conduisant à une inadéquation entre les seuils de débit fixés et ces cycles respiratoires. Cette inadéquation peut entraîner une mauvaise synchronisation entre le ventilateur et le patient, ou même un non déclenchement du ventilateur.
Ainsi, selon un mode de réalisation de l'invention, le procédé d'analyse de l'activité
respiratoire selon l'invention est mis en oeuvre pour déterminer, préalablement à
l'utilisation du respirateur, les valeurs de seuils des paramètres respiratoires les mieux adaptés pour un fonctionnement optimal du respirateur. Cette analyse est également
mesure de l'acquisition du signal d'activité respiratoire, de manière à surveiller en continu l'activité
respiratoire du patient.
Selon un autre mode de réalisation, le système et le procédé selon l'invention peuvent être mis en oeuvre dans un respirateur artificiel.
Un respirateur artificiel est un dispositif visant à pallier l'insuffisance respiratoire d'un patient, en insufflant de l'air vers les poumons du patient par l'intermédiaire d'un masque, au moyen d'un ventilateur. Pour que cette ventilation soit efficace, il est nécessaire que les efforts inspiratoires du patient et le déclenchement du ventilateur soient synchronisés, c'est-à-dire que l'air soit insufflé dans les poumons par le ventilateur seulement pendant la phase inspiratoire du cycle respiratoire du patient.
La mesure du débit d'air inspiré et expiré par le patient permet de détecter le début et la fin de cette phase inspiratoire, et ainsi de synchroniser le fonctionnement du ventilateur sur le cycle respiratoire du patient.
Cette détection est généralement réalisée en fixant deux seuils de débit ou de pression, le ventilateur étant actionné dès que le débit ou la pression d'air inspiré par le patient dépasse un premier seuil, et arrêté dès que ce débit ou cette pression devient inférieur à un deuxième seuil.
Cette méthode de détection peut s'avérer inefficace, en raison notamment de la variabilité des cycles respiratoires, conduisant à une inadéquation entre les seuils de débit fixés et ces cycles respiratoires. Cette inadéquation peut entraîner une mauvaise synchronisation entre le ventilateur et le patient, ou même un non déclenchement du ventilateur.
Ainsi, selon un mode de réalisation de l'invention, le procédé d'analyse de l'activité
respiratoire selon l'invention est mis en oeuvre pour déterminer, préalablement à
l'utilisation du respirateur, les valeurs de seuils des paramètres respiratoires les mieux adaptés pour un fonctionnement optimal du respirateur. Cette analyse est également
11 réalisée de manière continue pendant le fonctionnement du respirateur, de manière à
adapter ces seuils à la respiration du patient.
Selon un autre mode de réalisation, le système et le procédé selon l'invention peuvent être mis en oeuvre en imagerie tomodensitométrique (TMD), l'analyse des signaux d'activités respiratoires permettant de corriger, notamment sur les images des poumons, les artéfacts dus aux mouvements respiratoires.
Selon un autre mode de réalisation, le procédé de mesure de l'activité
respiratoire d'un patient comprend les étapes suivantes :
- placer une interface sur les voies respiratoires du patient, ladite interface comprenant un pneumotachographe apte à mesurer de manière continue le débit d'air inspiré et expiré par le patient ;
- acquérir au moins un signal d'activité respiratoire du patient comprenant au moins un signal élémentaire correspondant à un cycle respiratoire, dont la forme générale peut être exprimée par x(t) = xo +x1 cos(cI)(t)), où t est le temps, IO est la phase dudit signal élémentaire, xo et x, des coefficients, - analyser ledit signal d'activité respiratoire, où l'analyse dudit signal d'activité
respiratoire comprend les étapes suivantes :
- extraction, à partir dudit signal d'activité respiratoire, dudit signal élémentaire ;
\ de) - détermination d'une expression d'une équation de phaseF(c13) = ¨t dudit d signal élémentaire ; et - détermination d'une expression de la phase MO dudit signal élémentaire en fonction de paramètres mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcos, et psin,, définies par:
p cos n(t , r) .--Icos(kt) ¨r et p sin, (t, r) sin (k) rk , k" k"
k=1 k =1 où r est un paramètre d'anharmonicité mesurant le degré d'anharmonicité du signal, qui varie dans l'intervalle [0,1], la limite r=0 correspondant à un signal harmonique, la limite r=1 à un signal infiniment anharmonique, et où k et n sont des entiers.
Selon un autre mode de réalisation, le système de mesure de l'activité
respiratoire d'un patient comprend une interface pour recouvrir les voies respiratoires du patient,comprenant un pneumotachographe apte à mesurer de manière continue le débit 1 1 a d'air inspiré et expiré par le patient ; des moyens pour acquérir un signal d'activité
respiratoire du patient comprenant au moins un signal élémentaire correspondant à un cycle respiratoire dont la forme générale peut être exprimée par x(t)=. xo +
x1 cos(1)(0), où
t est le temps, (1)(t) est la phase dudit signal élémentaire, xo et xi des coefficients ; des moyens pour analyser ledit signal d'activité respiratoire, où les moyens pour analyser ledit signal d'activité respiratoire comprennent : des moyens pour extraire, à
partir dudit signal d'activité respiratoire, ledit signal élémentaire ; des moyens pour déterminer une \ ch) expression d'une équation de phase F((1)r ________________________________ dudit signal élémentaire ; et des moyens dt pour déterminer une expression de la phase (DO dudit signal élémentaire en fonction de paramètres (r, rk, *0, pk) mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcos, et psin, définies par:
p cos, (t, r) cos(kt) ' et p sin, (t,r)= sin(kt)¨rk , k=1 kn k =1 kn où r est un paramètre d'anharmonicité mesurant le degré d'anharmonicité du signal, qui varie dans l'intervalle [0,1], la limite r=0 correspondant à un signal harmonique, la limite r=1 à un signal infiniment anharmonique, et où k et n sont des entiers.
Il devra toutefois être compris que les exemples de réalisation présentés ci-dessus ne sont pas limitatifs.
Notamment, le signal d'activité respiratoire n'est pas nécessairement relatif au débit d'air inspiré ou expiré.
adapter ces seuils à la respiration du patient.
Selon un autre mode de réalisation, le système et le procédé selon l'invention peuvent être mis en oeuvre en imagerie tomodensitométrique (TMD), l'analyse des signaux d'activités respiratoires permettant de corriger, notamment sur les images des poumons, les artéfacts dus aux mouvements respiratoires.
Selon un autre mode de réalisation, le procédé de mesure de l'activité
respiratoire d'un patient comprend les étapes suivantes :
- placer une interface sur les voies respiratoires du patient, ladite interface comprenant un pneumotachographe apte à mesurer de manière continue le débit d'air inspiré et expiré par le patient ;
- acquérir au moins un signal d'activité respiratoire du patient comprenant au moins un signal élémentaire correspondant à un cycle respiratoire, dont la forme générale peut être exprimée par x(t) = xo +x1 cos(cI)(t)), où t est le temps, IO est la phase dudit signal élémentaire, xo et x, des coefficients, - analyser ledit signal d'activité respiratoire, où l'analyse dudit signal d'activité
respiratoire comprend les étapes suivantes :
- extraction, à partir dudit signal d'activité respiratoire, dudit signal élémentaire ;
\ de) - détermination d'une expression d'une équation de phaseF(c13) = ¨t dudit d signal élémentaire ; et - détermination d'une expression de la phase MO dudit signal élémentaire en fonction de paramètres mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcos, et psin,, définies par:
p cos n(t , r) .--Icos(kt) ¨r et p sin, (t, r) sin (k) rk , k" k"
k=1 k =1 où r est un paramètre d'anharmonicité mesurant le degré d'anharmonicité du signal, qui varie dans l'intervalle [0,1], la limite r=0 correspondant à un signal harmonique, la limite r=1 à un signal infiniment anharmonique, et où k et n sont des entiers.
Selon un autre mode de réalisation, le système de mesure de l'activité
respiratoire d'un patient comprend une interface pour recouvrir les voies respiratoires du patient,comprenant un pneumotachographe apte à mesurer de manière continue le débit 1 1 a d'air inspiré et expiré par le patient ; des moyens pour acquérir un signal d'activité
respiratoire du patient comprenant au moins un signal élémentaire correspondant à un cycle respiratoire dont la forme générale peut être exprimée par x(t)=. xo +
x1 cos(1)(0), où
t est le temps, (1)(t) est la phase dudit signal élémentaire, xo et xi des coefficients ; des moyens pour analyser ledit signal d'activité respiratoire, où les moyens pour analyser ledit signal d'activité respiratoire comprennent : des moyens pour extraire, à
partir dudit signal d'activité respiratoire, ledit signal élémentaire ; des moyens pour déterminer une \ ch) expression d'une équation de phase F((1)r ________________________________ dudit signal élémentaire ; et des moyens dt pour déterminer une expression de la phase (DO dudit signal élémentaire en fonction de paramètres (r, rk, *0, pk) mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcos, et psin, définies par:
p cos, (t, r) cos(kt) ' et p sin, (t,r)= sin(kt)¨rk , k=1 kn k =1 kn où r est un paramètre d'anharmonicité mesurant le degré d'anharmonicité du signal, qui varie dans l'intervalle [0,1], la limite r=0 correspondant à un signal harmonique, la limite r=1 à un signal infiniment anharmonique, et où k et n sont des entiers.
Il devra toutefois être compris que les exemples de réalisation présentés ci-dessus ne sont pas limitatifs.
Notamment, le signal d'activité respiratoire n'est pas nécessairement relatif au débit d'air inspiré ou expiré.
Claims (14)
1.- Procédé de mesure de l'activité respiratoire d'un patient, comprenant :
- placer une interface sur les voies respiratoires du patient, ladite interface comprenant un pneumotachographe apte à mesurer de manière continue le débit d'air inspiré et expiré par le patient ;
- acquérir au moins un signal d'activité respiratoire du patient comprenant au moins un signal élémentaire correspondant à un cycle respiratoire, dont la forme générale peut être exprimée par x(t)=x0+ x1, cos(.PHI.(t)), où t est le temps, (.PHI.(t) est la phase dudit signal élémentaire, x0 et x1 des coefficients, - analyser ledit signal d'activité respiratoire, où l'analyse dudit signal d'activité
respiratoire comprend les étapes suivantes :
- extraction, à partir dudit signal d'activité respiratoire, dudit signal élémentaire ;
- détermination d'une expression d'une équation de phase F(.PHI.)=~ dudit signal élémentaire ; et - détermination d'une expression de la phase .PHI.(t) dudit signal élémentaire en fonction de paramètres mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcosn et psinn définies par :
où r est un paramètre d'anharmonicité mesurant le degré d'anharmonicité du signal, qui varie dans l'intervalle [0,1], la limite r=0 correspondant à un signal harmonique, la limite r=1 à un signal infiniment anharmonique, et où k et n sont des entiers.
- placer une interface sur les voies respiratoires du patient, ladite interface comprenant un pneumotachographe apte à mesurer de manière continue le débit d'air inspiré et expiré par le patient ;
- acquérir au moins un signal d'activité respiratoire du patient comprenant au moins un signal élémentaire correspondant à un cycle respiratoire, dont la forme générale peut être exprimée par x(t)=x0+ x1, cos(.PHI.(t)), où t est le temps, (.PHI.(t) est la phase dudit signal élémentaire, x0 et x1 des coefficients, - analyser ledit signal d'activité respiratoire, où l'analyse dudit signal d'activité
respiratoire comprend les étapes suivantes :
- extraction, à partir dudit signal d'activité respiratoire, dudit signal élémentaire ;
- détermination d'une expression d'une équation de phase F(.PHI.)=~ dudit signal élémentaire ; et - détermination d'une expression de la phase .PHI.(t) dudit signal élémentaire en fonction de paramètres mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcosn et psinn définies par :
où r est un paramètre d'anharmonicité mesurant le degré d'anharmonicité du signal, qui varie dans l'intervalle [0,1], la limite r=0 correspondant à un signal harmonique, la limite r=1 à un signal infiniment anharmonique, et où k et n sont des entiers.
2.- Procédé selon la revendication 1, dans lequel l'équation de phase est exprimée sous la forme : dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité
dudit signal élémentaire.
dudit signal élémentaire.
3.- Procédé selon la revendication 2 dans lequel le signal élémentaire est exprimé
au moyen de deux paramètres r et .PHI.0, sous la forme :
x(t)= x0+ a1h sin(t,r) + b1h cos(t,r) où a1 = x1 cos(.PHI.0 ) et b1 = ¨x1 sin(.PHI.0), les fonctions hsin et hcos étant définies par :
au moyen de deux paramètres r et .PHI.0, sous la forme :
x(t)= x0+ a1h sin(t,r) + b1h cos(t,r) où a1 = x1 cos(.PHI.0 ) et b1 = ¨x1 sin(.PHI.0), les fonctions hsin et hcos étant définies par :
4.- Procédé selon la revendication 1, dans lequel l'équation de phase est exprimée sous la forme :
dans laquelle P(.PHI.) et Q(.PHI.) sont des polynômes trigonométriques.
dans laquelle P(.PHI.) et Q(.PHI.) sont des polynômes trigonométriques.
5.- Procédé selon la revendication 4, dans lequel l'expression de la phase .PHI.(t) est déterminée sous la forme :
dans laquelle les fonctions psin1 et pcos1 sont définies par :
dans laquelle les fonctions psin1 et pcos1 sont définies par :
6.- Procédé selon l'une quelconque des revendications 1 à 5, dans lequel ladite interface comprend un masque respiratoire et qui est placé sur le nez et la bouche du patient.
7.- Système de mesure de l'activité respiratoire d'un patient comprenant :
- une interface pour recouvrir les voies respiratoires du patient, comprenant un pneumotachographe apte à mesurer de manière continue le débit d'air inspiré et expiré
par le patient ;
- des moyens pour acquérir un signal d'activité respiratoire du patient comprenant au moins un signal élémentaire correspondant à un cycle respiratoire dont la forme générale peut être exprimée par x(t) = x0 + x1 cos(.PHI.(t)), où t est le temps, .PHI.(t) est la phase dudit signal élémentaire, x0 et x1 des coefficients ;
- des moyens pour analyser ledit signal d'activité respiratoire, où les moyens pour analyser ledit signal d'activité respiratoire comprennent :
- des moyens pour extraire, à partir dudit signal d'activité respiratoire, ledit signal élémentaire ;
- des moyens pour déterminer une expression d'une équation de phase F(.PHI.) = ~ dudit signal élémentaire ; et - des moyens pour déterminer une expression de la phase (.PHI.)(t) dudit signal élémentaire en fonction de paramètres (r, rk, .PHI.0, pk) mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcosn et psinn définies par :
où r est un paramètre d'anharmonicité mesurant le degré d'anharmonicité du signal, qui varie dans l'intervalle [0,1], la limite r=0 correspondant à un signal harmonique, la limite r=1 à un signal infiniment anharmonique, et où k et n sont des entiers.
- une interface pour recouvrir les voies respiratoires du patient, comprenant un pneumotachographe apte à mesurer de manière continue le débit d'air inspiré et expiré
par le patient ;
- des moyens pour acquérir un signal d'activité respiratoire du patient comprenant au moins un signal élémentaire correspondant à un cycle respiratoire dont la forme générale peut être exprimée par x(t) = x0 + x1 cos(.PHI.(t)), où t est le temps, .PHI.(t) est la phase dudit signal élémentaire, x0 et x1 des coefficients ;
- des moyens pour analyser ledit signal d'activité respiratoire, où les moyens pour analyser ledit signal d'activité respiratoire comprennent :
- des moyens pour extraire, à partir dudit signal d'activité respiratoire, ledit signal élémentaire ;
- des moyens pour déterminer une expression d'une équation de phase F(.PHI.) = ~ dudit signal élémentaire ; et - des moyens pour déterminer une expression de la phase (.PHI.)(t) dudit signal élémentaire en fonction de paramètres (r, rk, .PHI.0, pk) mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcosn et psinn définies par :
où r est un paramètre d'anharmonicité mesurant le degré d'anharmonicité du signal, qui varie dans l'intervalle [0,1], la limite r=0 correspondant à un signal harmonique, la limite r=1 à un signal infiniment anharmonique, et où k et n sont des entiers.
8.- Système selon la revendication 7, qui comporte des moyens pour exprimer l'équation de phase sous la forme :
dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité
dudit signal élémentaire.
dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité
dudit signal élémentaire.
9.- Système selon la revendication 8, qui comporte des moyens pour exprimer ledit signal élémentaire au moyen de deux paramètres r et .PHI.0, sous la forme :
x(t) = x0 + a1hsin(t,r)+ b1h cos(t,r) où a1 = x1 cos(.PHI.0) et b1 = ¨x1 sin(.PHI.0), les fonctions hsin et hcos étant définies par :
x(t) = x0 + a1hsin(t,r)+ b1h cos(t,r) où a1 = x1 cos(.PHI.0) et b1 = ¨x1 sin(.PHI.0), les fonctions hsin et hcos étant définies par :
10.- Système selon la revendication 9, qui comporte des moyens pour exprimer l'équation de phase sous la forme:
dans laquelle P(.PHI.) et Q(.PHI.) sont des polynômes trigonométriques.
dans laquelle P(.PHI.) et Q(.PHI.) sont des polynômes trigonométriques.
11.- Système selon la revendication 10, qui comporte des moyens pour exprimer la phase .PHI.(t) sous la forme :
dans laquelle les fonctions psin1 et pcos1 sont définies par :
dans laquelle les fonctions psin1 et pcos1 sont définies par :
12. ¨ Système selon l'une quelconque des revendications 7 à 11, dans lequel ladite interface comprend un masque respiratoire pour être placé sur le nez et la bouche du patient.
13.- Dispositif d'assistance respiratoire comprenant un système de mesure de l'activité respiratoire selon l'une quelconque des revendications 7 à 12.
14.- Dispositif de surveillance respiratoire comprenant un système de mesure de l'activité respiratoire selon l'une quelconque des revendications 7 à 12.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1053357A FR2959407A1 (fr) | 2010-04-30 | 2010-04-30 | Procede et systeme d'analyse de l'activite respiratoire d'un patient et applications correspondantes |
FR1053357 | 2010-04-30 | ||
PCT/FR2011/050957 WO2011135257A2 (fr) | 2010-04-30 | 2011-04-28 | Procédé et système d'analyse de l'activité respiratoire d'un patient et applications correspondantes |
Publications (2)
Publication Number | Publication Date |
---|---|
CA2797886A1 CA2797886A1 (fr) | 2011-11-03 |
CA2797886C true CA2797886C (fr) | 2018-07-03 |
Family
ID=43661798
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA2797886A Active CA2797886C (fr) | 2010-04-30 | 2011-04-28 | Procede et systeme d'analyse de l'activite respiratoire d'un patient et applications correspondantes |
Country Status (6)
Country | Link |
---|---|
US (1) | US9649051B2 (fr) |
EP (1) | EP2564326A2 (fr) |
JP (1) | JP5932769B2 (fr) |
CA (1) | CA2797886C (fr) |
FR (1) | FR2959407A1 (fr) |
WO (1) | WO2011135257A2 (fr) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NZ736250A (en) * | 2012-04-13 | 2019-04-26 | ResMed Pty Ltd | Apparatus and methods for ventilatory treatment |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5020540A (en) | 1987-10-09 | 1991-06-04 | Biometrak Corporation | Cardiac biopotential analysis system and method |
US5490502A (en) * | 1992-05-07 | 1996-02-13 | New York University | Method and apparatus for optimizing the continuous positive airway pressure for treating obstructive sleep apnea |
US6306088B1 (en) * | 1998-10-03 | 2001-10-23 | Individual Monitoring Systems, Inc. | Ambulatory distributed recorders system for diagnosing medical disorders |
US6644312B2 (en) * | 2000-03-07 | 2003-11-11 | Resmed Limited | Determining suitable ventilator settings for patients with alveolar hypoventilation during sleep |
JP5451013B2 (ja) * | 2007-09-06 | 2014-03-26 | 英次 麻野井 | 睡眠の質を評価するために用いる表示又は印刷を行う装置、方法、コンピュータプログラム、呼吸補助装置、慢性心疾患患者を対象とした呼吸補助装置、睡眠導入装置、マッサージ装置、検査装置 |
JP5207172B2 (ja) * | 2008-03-21 | 2013-06-12 | 国立大学法人 大分大学 | 波形解析装置及び波形解析プログラム |
FR2955187A1 (fr) * | 2010-01-08 | 2011-07-15 | Centre Nat Rech Scient | Procede de decomposition d'un signal periodique anharmonique et programme d'ordinateur correspondant |
FR2958529A1 (fr) * | 2010-04-13 | 2011-10-14 | Centre Nat Rech Scient | Procede et systeme d'analyse de l'activite cardiaque d'un patient et applications correspondantes |
-
2010
- 2010-04-30 FR FR1053357A patent/FR2959407A1/fr active Pending
-
2011
- 2011-04-28 WO PCT/FR2011/050957 patent/WO2011135257A2/fr active Application Filing
- 2011-04-28 EP EP11723517A patent/EP2564326A2/fr not_active Ceased
- 2011-04-28 CA CA2797886A patent/CA2797886C/fr active Active
- 2011-04-28 US US13/695,266 patent/US9649051B2/en active Active
- 2011-04-28 JP JP2013506717A patent/JP5932769B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
WO2011135257A2 (fr) | 2011-11-03 |
EP2564326A2 (fr) | 2013-03-06 |
FR2959407A1 (fr) | 2011-11-04 |
CA2797886A1 (fr) | 2011-11-03 |
US20130096452A1 (en) | 2013-04-18 |
JP5932769B2 (ja) | 2016-06-08 |
WO2011135257A3 (fr) | 2012-01-12 |
US9649051B2 (en) | 2017-05-16 |
JP2013528415A (ja) | 2013-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2719169C (fr) | Dispositif de detection de l'observance par un patient d'un traitement d'oxygenotherapie | |
US9706946B2 (en) | Spirometer system and methods of data analysis | |
US8425428B2 (en) | Nitric oxide measurements in patients using flowfeedback | |
JP6199312B2 (ja) | チェーンストークス呼吸パターンの判別 | |
JP2008525060A5 (fr) | ||
US20130190642A1 (en) | Apparatus and method for diagnosing obstructive sleep apnea | |
EP3110321A1 (fr) | Procede, dispositif, systeme et programme informatique de filtrage d'une serie rr obtenue a partir d'un signal cardiaque avec controle automatique de la qualite de la serie rr | |
EP2590701B1 (fr) | Système de détection électroencéphalographique d'une inadéquation entre l'état d'un patient placé sous assistance ventilatoire et le réglage de la machine utilisée pour cette assistance, et utilisation de cette détection pour l'adaptation du réglage | |
EP0892618A1 (fr) | Dispositif de determination de la profondeur d'anesthesie | |
EP2364108A1 (fr) | Systeme de detection de l'activite musculaire respiratoire d'un patient sous assistance respiratoire | |
CA2797886C (fr) | Procede et systeme d'analyse de l'activite respiratoire d'un patient et applications correspondantes | |
EP2364174B1 (fr) | Système de quantification du desaccord entre un patient sous assistance respiratoire et un appareil d'assistance correspondant | |
JP2007190082A (ja) | 肺音診断装置及び肺音診断方法 | |
EP3600010A1 (fr) | Systeme de determination d'un ensemble d'au moins un descripteur cardio-respiratoire d'un individu pendant son sommeil | |
Sa et al. | Validation of the Forced Oscillation Technique in the diagnostic of respiratory changes in patients with silicosis | |
WO2015128568A1 (fr) | Procede et dispositif de contrôle automatique de la qualite d'une serie rr obtenue a partir d'un signal cardiaque | |
WO2007097634A1 (fr) | Procédé et appareil de surveillance d'un patient sous sédation | |
CA2796172C (fr) | Procede et systeme d'analyse de l'activite cardiaque d'un patient et applications correspondantes | |
Galli et al. | Accurate ECG monitoring by Gaussian feature streaming | |
Terrill et al. | Investigating parameters participating in the infant respiratory control system attractor | |
Uchida et al. | An ICA approach to reject noise from pressure changes of pillow | |
Ko et al. | An EEG Spike Prescreening Algorithm Using RBF ANN with Multi-Channel Correlation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
EEER | Examination request |
Effective date: 20160321 |