Brevet FR3017529A1 — publié le 20150821

DOMAINE TECHNIQUE DE L'INVENTION La présente invention concerne en général un procédé et un système de surveillance d'un système nerveux autonome (SNA).L'invention propose par ailleurs, un procédé et un système d'aide à la régulation de l'activité du SNA. Elle trouve pour avantage particulier les systèmes de surveillance ou de régulation qui nécessitent un suivi du SNA en temps réel ou tout au moins avec une très bonne réactivité. ÉTAT DE LA TECHNIQUE Le système nerveux autonome (SNA) est la partie du système nerveux responsable des fonctions automatiques, non soumises au contrôle volontaire d'un être vivant tel qu'une personne. Le SNA, est composé de deux branches, respectivement le système sympathique et le système parasympathique. Les rôles des systèmes sympathique et parasympathique sont différents et ils exercent généralement des effets antagonistes sur les mêmes organes cibles. D'une manière générale, le système sympathique aide le corps à répondre à des situations d'urgence par accélérer les réactions des neurones afin d'augmenter par exemple la fréquence des battements cardiaques et de ralentir le processus de digestion. En revanche, le système parasympathique aide le corps à préserver et à restaurer son énergie. Lorsqu'une personne se détend, par exemple en se reposant dans une chaise, le système parasympathique abaisse sa pression sanguine, ralentit la fréquence des battements cardiaques et accélère le processus de digestion. Un stress prolongé génère une activité émotionnelle négative qui stimule le système sympathique en créant ainsi un déséquilibre constant qui provoque de nombreuses symptômes. Afin de fournir à des chercheurs et médecins des données précises pour effectuer des diagnostics et décider des traitements adaptés, il existe alors un besoin consistant à surveiller de manière fiable et précise les systèmes sympathique et parasympathique et à en donner une représentation aisément interprétable. Une interprétation du SNA peut permettre d'identifier une pathologie particulière comme par exemple un état dans lequel le système sympathique d'un patient est anormalement actif par rapport au système parasympathique. Ainsi, le médecin peut effectuer des diagnostics et décider des remèdes appropriés. Il existe de nombreuses méthodes pour surveiller et analyser le SNA. La variabilité de la fréquence cardiaque (VFC ; acronyme de l'anglais « Heart Rate Variability », abréviation HRV) est en général utilisée pour évaluer les fonctions physiologiques telles que les activités des systèmes sympathique et parasympathique. La VFC peut être générée par exemple en utilisant des données d'une séquence d'intervalles d'inter-battements dans un signal physiologique tel qu'un signal d'électrocardiogramme (ECG). Le plus souvent, un intervalle/décalage entre deux battements cardiaques successifs, dénommé « intervalle RR », est défini comme un intervalle de temps entre deux sommets préférentiellement successifs d'une onde R d'un signal ECG. Il existe plusieurs méthodes d'analyse et de traitement de signal utilisant la VFC, dont la majorité est basée sur une analyse dans le domaine fréquentiel. Néanmoins, les analyses fréquentielles présentent plusieurs inconvénients notamment en termes de précision et de temps d'exécution. En effet, ces analyses fréquentielles exigent de mettre en oeuvre des fonctions mathématiques telles que la transformée de Fourier rapide habituellement désignée par l'acronyme FFT (de l'anglais « Fast Fourier Transformation »,) ou/et la méthode de démodulation complexe (CDM, acronyme de l'anglais «Complex Demodulation Method »), décrites dans les documents US7079888 et US20130009779A1. Or ces fonctions mathématiques nécessitent un temps d'enregistrement long et ne peuvent donc pas être réalisées en temps réel. De plus, il est difficile de séparer correctement dans le domaine fréquentiel les données appartenant respectivement au système sympathique et au système parasympathique. Autrement dit, les deux systèmes sympathique et parasympathique ne peuvent pas facilement être analysés de manière indépendante. Les analyses fréquentielles sont en outre très sensibles aux bruits d'un signal ECG tels que les faux battements générés lors de la capture du signal par exemple dans le cas où la personne se déplace lors d'une activité sportive ou si le capteur ECG ne fonctionne pas correctement. Les analyses fréquentielles sont sensibles également aux battements atypiques tels que les extrasystoles, arythmies (bigéminées, trigéminées), etc. Ainsi, les données transformées dans le domaine fréquentiel sont très fortement faussées du fait de bruits ou de faux battements, ce qui réduit considérablement la précision de la surveillance du SNA. Les analyses fréquentielles basées sur les méthodes connues sont donc très peu fiables dès lors que le patient présente des symptômes entrainant des battements atypiques, ou/et dans le cas où de faux battements sont générés par exemple par un détecteur peu performant ou fonctionnant dans des conditions d'acquisition défavorables. D'autres solutions proposent des analyses effectuées dans le domaine temporel. Ces solutions connues prévoient, comme c'est le cas dans le document US 20130079652A1, des calculs tels qu'une mesure sélectionnée parmi un ensemble de mesures comprenant un intervalle RR moyen, un écart-type des intervalles RR, un ratio d'écarts-types suivant de différents axes d'un diagramme de dispersion des intervalles RR, une différence entre ratios des écarts-types, un écart-type d'une moyenne d'intervalles RR dans différents segments de temps etc. Les calculs d'analyses temporelles comme mentionnés ci-dessus peuvent être réalisés en temps réel. En revanche, ces calculs ne permettent pas plus d'analyser les deux systèmes sympathique et parasympathique de manière indépendante. Par exemple, si le système sympathique d'un patient reste actif plus longtemps qu'une période prédéterminée et empêche le système parasympathique de jouer son rôle de récupération, cela ne peut pas être reflété par le résultat de ces calculs d'analyses temporelles. Ces calculs ne permettent pas d'analyser l'état d'équilibre entre les deux systèmes sympathique et parasympathique du SNA et ainsi ne permettent pas de déterminer avec précision le niveau de stress d'un patient. Par conséquent, il existe un besoin consistant à proposer une solution permettant de fournir une analyse fiable, simple et rapide pour surveiller les systèmes sympathique et parasympathique de manière indépendante, avec une bonne réactivité et de préférence en temps réel et tout en limitant la sensibilité à des facteurs perturbants tels que des battements atypiques et des bruits comme de faux battements générés lors de capture de signaux ECG. RÉSUMÉ DE L'INVENTION Pour atteindre cet objectif, un aspect de la présente invention concerne un procédé de surveillance d'un système nerveux autonome (SNA) d'un sujet, comprenant : - l'acquisition d'au moins un signal d'électrocardiogramme (ECG), ce signal ECG comprenant W battements cardiaques, chaque battement cardiaque étant détecté à un instant t, avec t compris entre tx_vv+i et tx, l'instant tx étant le dernier battement des W battements cardiaques ; - la génération d'une donnée qui est fonction de la variabilité de la fréquence cardiaque (VFC) sur l'ensemble des W battements cardiaques, la VFC comprenant W-1 intervalles RR séparant deux battements cardiaques consécutifs détectés respectivement aux instants tk_i et tk, chaque intervalle RR présentant une durée d'une valeur ak = tk tk-1 avec k = (X-W+2)... X ; - le calcul à l'aide d'au moins un microprocesseur d'au moins un paramètre pris parmi : - le paramètre TAS(tx) représentant le niveau d'activité du système sympathique et qui est fonction du rapport entre d'une part une première somme obtenue en sommant, parmi les valeurs ak, uniquement celles qui sont supérieures à la valeur ak_i de l'intervalle immédiatement précédent donc uniquement les valeurs ak qui respectent la condition suivante ak > ak-1 étant la durée de l'intervalle immédiatement précédent celui de la valeur ak et 30 une deuxième somme obtenue en sommant toutes les W-1 valeurs ak; et/ou - le paramètre TAP(tx) représentant le niveau d'activité du système parasympathique et qui est fonction du rapport entre d'une part une première somme obtenue en sommant, parmi les valeurs ak, uniquement celles qui sont inférieures à la valeur ak_i de l'intervalle immédiatement précédent et une deuxième somme obtenue en sommant toutes les W-1 valeurs ak ; - le paramètre TPSP(tx) représentant le taux de pollution du système 5 sympathique vers le système parasympathique et qui est fonction du rapport entre d'une part une première somme obtenue en sommant, parmi les valeurs ak, uniquement celles qui respectent la condition suivante : ak ak-2 et ak > ak et d'autre part une deuxième somme obtenue en sommant toutes les W-1 valeurs ak ; et/ou 10 - le paramètre TPPS(tx) représentant le taux de pollution du système parasympathique vers le système sympathique et qui est fonction du rapport entre d'une part une première somme obtenue en sommant, parmi les valeurs ak, uniquement celles qui respectent la condition suivante : (ak > ak_i et ak_i Optionnellement, l'invention comprend au moins l'une quelconque des caractéristiques et étapes suivantes prises séparément ou en combinaison : Dans les équations suivantes ak_i = valeur de l'intervalle de temps immédiatement précédent l'intervalle de temps de durée ak, soit ak-1 = tk-1 tk-2, et k = (X-W+2), (X-W+3)... X. - le TAS(tx) représentant le niveau d'activité du système sympathique et dans lequel le TAS(tx) est calculé en appliquant l'équation suivante : X (tk -tk-i) si (akG ak-i) TAS (tx) =100 -W +2 (tk k-X -W +2 - le TAP(tx) représentant le niveau d'activité du système parasympathique et dans lequel le TAP(tX) est calculé en appliquant l'équation suivante : X (tk - tk-i) si (aki ak-i) TAP (tx) =100 k-X -W +2 (tk tk-1) k =X -W +2 - le TPSP(tx) étant calculé en appliquant l'équation suivante : (tk - tk_i) si (ak ak_2 et ak' > ak) TPSP(tx) =100 k =X -W +2 (tk tk-1) k-X -W +2 Ainsi, le TPSP(tx) est calculé de manière à obtenir une troisième somme partielle par accumulation des valeurs ak dont chacune perturbe la tendance monotone croissante formée par des intervalles RR adjacents à l'intervalle RR X X X correspondant à la valeur ak et à diviser ladite somme troisième somme partielle par la somme de toutes les W-1 valeurs ak ; - le TPPS(tx) étant calculé en appliquant l'équation suivante : (t, - t,_,) si (ak > a,_, et a,_, an_2 et an, > an) TPSP = 100 n=1 N -1 - le paramètre TPPS représentant le taux approximé de pollution du système parasympathique vers le système sympathique et calculé en appliquant 20 l'équation suivante : N-1 1 si (an > an_, et an_, L'invention concerne une étape de fourniture d'une représentation visuelle et/ou auditive dudit au moins un paramètre. De préférence, ledit au moins un paramètre est calculé en temps réel ou à intervalles réguliers et comprend une étape de fourniture d'une représentation visuelle de l'évolution au cours du temps dudit au moins un paramètre. De manière avantageuse, l'étape de fourniture d'une représentation visuelle comprend l'affichage d'un graphique ou d'une barre évoluant au cours du temps. Dans une autre application de l'invention, l'étape de fourniture d'une représentation visuelle de l'évolution au cours du temps dudit au moins un paramètre comprend l'affichage d'au moins un avatar et/ou d'au moins un objet animé dont l'animation ou l'évolution visuelle est fonction de l'évolution au cours du temps dudit au moins un paramètre. De préférence, on calcule au moins deux paramètres, au moins l'un des paramètres étant pris parmi le TAS(tx) et le TPSP(tx) et au moins un autre paramètres étant pris parmi le TAP(tx) et le TPPS(tx) et comprenant une étape de fourniture d'une représentation visuelle desdits au moins deux paramètres, lesdits au moins deux paramètres étant représentés indépendamment l'un de l'autre. De manière encore plus avantageuse, lesdits au moins deux paramètres sont calculés en temps réel ou à intervalles réguliers et comprennent une étape de fourniture d'une représentation visuelle de l'évolution au cours du temps 5 desdits au moins deux paramètres. L'invention concerne également un produit programme d'ordinateur comprenant des instructions, qui lorsqu'elles sont effectuées par au moins un processeur, exécute le procédé de surveillance comme ci-dessus. 10 Selon un mode de réalisation, les paramètres TAS(tx), TAP(tx), TPSP(tx), TPPS(tx), NS(tx) et NSR(tx) peuvent être utilisés dans des jeux vidéo et des simulateurs afin de montrer la condition actuelle du corps d'un joueur telle que le niveau de stress, le niveau d'énergie, ou de simuler une condition virtuelle 15 d'un personnage dans un jeu vidéo ou par un simulateur. De plus, l'invention concerne un système de surveillance du système nerveux autonome (SNA) d'un sujet, comprenant un dispositif configuré pour recevoir des données d'un signal d'électrocardiogramme (ECG), ce signal ECG 20 comprenant W battements cardiaques, chaque battement cardiaque étant détecté à un instant t, avec t compris entre tx_vv+i et tx, l'instant tx étant le dernier battement des W battements cardiaques; et comprenant au moins un traitement des données équipé d'au moins un processeur configuré pour exécuter les étapes suivantes : 25 - génération d'une donnée qui est fonction de la variabilité de la fréquence cardiaque (VFC) sur l'ensemble des W battements cardiaques, la VFC comprenant W-1 intervalles RR séparant deux battements cardiaques consécutifs détectés respectivement aux instants tk_i et tk, chaque intervalle RR présentant une durée ak = tk tk-1 d'une valeur avec k = (X-W+2)... X ; 30 - calcul d'au moins un paramètre pris parmi : - le paramètre TAS(tx) représentant le niveau d'activité du système sympathique et qui est fonction du rapport entre d'une part une première somme obtenue en sommant, parmi les valeurs ak, uniquement celles qui sont supérieures à la valeur ak_i de l'intervalle immédiatement précédent donc uniquement les valeurs ak qui respectent la condition suivante ak > ak-1 étant la durée de l'intervalle immédiatement précédent celui de la valeur ak et une deuxième somme obtenue en sommant toutes les W-1 valeurs ak; et/ou - le paramètre TAP(tx) représentant le niveau d'activité du système parasympathique et qui est fonction du rapport entre d'une part une première somme obtenue en sommant, parmi les valeurs ak, uniquement celles qui sont inférieures à la valeur ak_i de l'intervalle immédiatement précédent donc uniquement les valeurs ak qui respectent la condition suivante ak ak-2 et ak > ak et d'autre part une deuxième somme obtenue en sommant toutes les W-1 valeurs ak; et/ou - le paramètre TPPS(tx) représentant le taux de pollution du système parasympathique vers le système sympathique et qui est fonction du rapport 20 entre d'une part une première somme obtenue en sommant, parmi les valeurs ak, uniquement celles qui respectent la condition suivante : (ak > ak_i et ak_i L'invention concerne un appareil d'aide à la régulation du SNA d'un sujet comprenant le système de surveillance comme ci-dessus. L'invention concerne en outre un simulateur comprenant le système de surveillance comme ci-dessus. BRÈVE DESCRIPTION DES FIGURES Les buts, objets, ainsi que les caractéristiques et avantages de l'invention ressortiront mieux de la description détaillée d'un mode de réalisation de cette dernière qui est illustré par les dessins d'accompagnement suivants dans lesquels : La FIGURE 1 illustre un système de surveillance d'un système nerveux autonome (SNA) d'un être vivant tel qu'une personne selon un exemple de mode de réalisation de l'invention. La FIGURES 2(a) montre un signal ECG comprenant N sommets signifiant respectivement un battement cardiaque de la personne 10 mesuré à un instant to, t1 à par le capteur ECG 110 selon un mode de réalisation de l'invention. La FIGURE 2(b) montre la variabilité de la fréquence cardiaque VFC illustrée par N-1 intervalles RR a1, a2 à awi correspondant respectivement aux instants t2 à tN-1- La FIGURE 3 montre des résultats des calculs de deux paramètres TAS et TAP sur une durée d'échantillonnage TE de 100 battements cardiaques (N = 100). La FIGURE 4 montre les paramètres TAS(tx) et TAP(tx) obtenus en temps réel par calcul basé sur 30 échantillons de la fenêtre glissante à chaque instant tx (W = 30). Les FIGURES 5(a) et 5(b) montrent chacune un spectre de la VFC selon deux types de respirations. La FIGURE 6(a) montre une partie d'un premier signal ECG et sa VFC, le 30 premier signal ECG comprenant 300 battements sans ajout de faux battements cardiaques. La FIGURE 6(b) montre une partie d'un deuxième signal ECG et sa VFC, le deuxième signal ECG est obtenu par ajout de deux faux battements au premier signal ECG. Les FIGURES 6(c) et 6(d) illustrent respectivement les spectres de VFC obtenues par réalisation d'une analyse fréquentielle typique de la VFC du premier signal ECG et ceux pour le deuxième signal ECG. La FIGURE 7(a) illustre les résultats d'un exemple d'un calcul du paramètre NS. La FIGURE 7(b) illustre les résultats d'un exemple d'un calcul des paramètres NS(tx), NSR(tx) et RC(tx) pendant une phase de repos, une phase d'activité sportive et une phase de récupération. Les dessins joints sont donnés à titre d'exemples et ne sont pas limitatifs de l'invention. Ces dessins sont des représentations schématiques et ne sont pas nécessairement à l'échelle de l'application pratique. DESCRIPTION DÉTAILLÉE DE L'INVENTION La figure 1 illustre un système de surveillance 100 du système nerveux autonome (SNA) d'un sujet physique 10 selon un mode de réalisation de l'invention. Le système de surveillance 100, comprend un module de traitement des données 222 et un module d'affichage des données 130. Il est configuré pour analyser un signal physiologique reçu d'un capteur de signaux physiologiques 110, typiquement un électrocardiogramme (ECG). Le capteur de signaux physiologiques 110 est configuré pour détecter des battements cardiaques de la personne physique 10 pendant une certaine durée et envoyer ainsi un signal physiologique au module de traitement des 30 données 222. Le capteur de signaux physiologiques 110 est de préférence un capteur électrocardiogramme (ECG) configuré pour générer un signal physiologique ECG tel qu'illustré en figure 2(a) et envoyer ce signal ECG au module de traitement des données 222. Ce capteur ECG 110 est apte à être porté par la personne 10 par exemple sur son thorax ou/et son poignet. L'invention n'est pas limitée à des méthodes d'acquisition d'une fréquence / d'un rythme cardiaque. Par exemple, un signal indiquant des battements cardiaques mesurés par des dispositifs optiques peuvent également être appliqués comme source de données de calculs de l'invention. Le module de traitement des données 222 est couplé, par une communication filaire ou sans fil tel que par ondes radio, au capteur ECG 110. Dans le présent mode de réalisation, un module de communication filaire 113 est utilisé pour réaliser la communication entre le capteur ECG 110 et le module de traitement des données 222. Le module de traitement des données 222 est configuré pour exécuter un procédé d'analyse de préférence dans le domaine temporel à partir d'un signal physiologique tel qu'un signal ECG reçu du capteur ECG 110. L'invention n'est pas limitée à des moyens de communication utilisés pour la réalisation de la communication filaire ou sans fil entre le module de traitement des données 222 et le capteur ECG 110. Ledit procédé d'analyse effectué par le module de traitement des données 222 comprend un calcul d'intervalles RR et des calculs des paramètres TAS, TAP, TPSP, TPPS qui seront décrits ultérieurement. La figure 2(a) montre un signal ECG comprenant N sommets signifiant 30 respectivement un battement cardiaque de la personne 10 mesuré à un instant to, t1 à twi par le capteur 110. La figure 2(b) montre une variabilité de la fréquence cardiaque VFC illustrée par N-1 intervalles RR al, a2 à aNi_i correspondant respectivement aux instants t1, t2 à tN-1. Chaque valeur a, d'un intervalle RR est la durée (typiquement en secondes) qui s'écoule entre deux sommets adjacents du signal ECG, ces deux sommets correspondant aux deux battements cardiaques successifs mesurés aux instants t,_1 et tn. La valeur a, est représentée dans l'équation (1) ci-dessous : a, =t, -t11_1 Eq. (1) Par exemple, comme indiqué en figure 2(a), la valeur al correspond à l'intervalle de temps RR entre les deux battements cardiaques détectés aux deux instants to et t1. La valeur du premier intervalle est notée al et le nombre n des intervalles RR est égal à N-1. Le module de traitement des données 222 effectue, à partir des N-1 intervalles RR al, a2 à aNi_i de la VFC, des calculs des paramètres TAS et TAP présentant respectivement l'activité du système sympathique et celle du système parasympathique de la personne 10 surveillée durant une durée d'échantillonnage TE comprenant N battements cardiaques du signal ECG. Les calculs des paramètres TAS et TAP utilisent un phénomène physiologique : une durée d'un intervalle RR qui diminue entre deux battements cardiaques successifs indique un rythme cardiaque accéléré par le système sympathique alors qu'une durée d'un intervalle RR qui augmente entre deux battements successifs indique un rythme cardiaque ralenti par le système parasympathique. Le paramètre TAS est configuré pour représenter le taux d'activité sympathique : N -1 N -1 (tn -tn-1) si (anG an-1) an si (anG an-1) TAS =100. n=1 -100 n=1 Eq. (2) N -1 N -1 (tn - t n n=1 n=1 Le paramètre TAS est calculé, à partir des N-1 intervalles RR dérivés du 5 signal ECG comprenant N battements cardiaques. Pour obtenir ce paramètre, on accumule (c'est-à-dire que l'on somme) les valeurs a, si un segment entre deux valeurs a, et an_i est descendant comme illustré en figure 2(b) ; c'est-à-dire on accumule des valeurs a, dont chacune est inférieure à la valeur immédiatement précédente an_i. La somme des valeurs a, correspondant aux 10 segments descendants est ensuite divisée par la somme de toutes les valeurs al, a2 à aN-i. Le paramètre TAP est configuré pour représenter le taux d'activité parasympathique : 15 N -1 N -1 (tn -tn-1) si (ani an-1) an si (ani an-1) TAP =100. n=1 -100 n=1 Eq. (3) N -1 N -1 (tn-tn_1) n n=1 n=1 Le paramètre TAP est calculé, à partir des N-1 intervalles RR de la VFC dérivée du signal ECG comprenant N battements cardiaques, de manière à 20 accumuler des valeurs a, si un segment entre deux valeurs a, et an_i est montant comme illustré en figure 2(b) ; c'est-à-dire on accumule des valeurs a, dont chacune est supérieure à la valeur immédiatement précédente an_i. La somme des valeurs a, correspondant aux segments montants est ensuite divisée par la somme de toutes les valeurs al, a2 à aN-i. 25 La figure 3 montre des résultats des calculs des paramètres TAS et TAP sur une durée d'échantillonnage TE de 100 battements cardiaques (N=100). La représentation visuelle de la figure 3, présente pour avantage de pouvoir comparer très facilement les activités respectives des systèmes sympathique et parasympathique. Dans cet exemple, le calcul des paramètres TAS et TAP permet d'identifier très aisément que le système SNA de la personne 10 est dans un équilibre relatif et imparfait car le système sympathique est légèrement plus actif que celui parasympathique. Par ailleurs, cette représentation varie en temps réel afin de suivre en temps réel l'évolution des activités respectives des systèmes sympathique et parasympathique. Ce suivi en temps réel permet par exemple au patient de contrôler son SNA en effectuant des actions appropriées, notamment pour augmenter l'activité de son système parasympathique. Au niveau de la fiabilité du calcul d'un point de vue mathématique, c'est- à-dire dans une condition idéale où le signal ECG est bien capturé par le capteur ECG 110, les calculs des paramètres TAS et TAP sont, par rapport à ceux des méthodes connues d'analyse fréquentielle, beaucoup plus précis. Le nombre N de battements peut être réduit à un nombre compris entre 20 et 40. Des comparaisons des deux types de calculs (méthode de l'invention et méthodes connues d'analyse fréquentielle) seront fournies dans les figures 5(a), 5(b) et 6(a) à 6(d) illustrées ultérieurement. Néanmoins, par exemple lorsque le signal ECG de la personne 10 est capturé lors d'une activité sportive, le nombre N de battements est compris de préférence entre 200 et 350 (c'est-à-dire la durée d'enregistrement des battements d'un signal ECG est environ de trois minutes), car en pratique, 200 à 350 battements successifs permettent de refléter un état physiologique pertinent. En effet, les calculs proposés par l'invention font appel à des implémentations mathématiques simples, ce qui permet d'obtenir en temps réel un résultat de calcul, notamment dans un mode de réalisation où le nombre N de battements cardiaques est inférieur à, par exemple, 1000. Ainsi, deux paramètres TAS(tx) et TAP(tx) montrant respectivement l'activité du système sympathique et celle du système parasympathique, peuvent être calculés, en temps réel, en utilisant par exemple un certain nombre de derniers battements cardiaques mesurés par le capteur ECG 110. De manière préférentielle mais non limitative, une fenêtre glissante d'échantillonnage peut être utilisée pour échantillonner W battements cardiaques qui seront utilisés pour calculer les paramètres TAS(tx) et TAP(tx) d'un signal ECG. L'invention n'est pourtant pas limitée à des méthodes de sélection ou d'échantillonnage des battements cardiaques utilisés pour les calculs des paramètres TAS(tx) et TAP(tx). Avec une fenêtre glissante dont le nombre d'échantillons est W, le module de traitement des données 222 peut effectuer, à chaque instant tel un nouveau battement cardiaque mesuré à l'instant tx du signal ECG, des calculs des paramètres TAS(tx) et TAP(tx) pour obtenir en temps réel des résultats de calculs qui montrent l'activité du système sympathique et celle du système parasympathique entre deux battements cardiaques (i.e. le (X-W+1 )ème et le xème battements cardiaques) respectivement aux instants tx_vv+i et tx du signal ECG. Le nombre des intervalles RR entre les deux battements cardiaques aux instants tx_vv+i et tx est donc W-1. De manière préférentielle mais non limitative, les calculs des paramètres TAS(tx) et TAP(tx) sont effectués après le Wème battement cardiaque du signal ECG. Par exemple, pour obtenir un résultat reflétant les activités des systèmes sympathique et parasympathique, le nombre d'échantillons W est compris par exemple entre 10 et 30. Le calcul du paramètre TAS(tx) est représenté dans l'équation (4) ci-dessous : X (tk -tk-il si (a,G ak_i) TAS (tx) = 100 -W +2 Eq. (4), dans X (tk - tk_i) k-X -W +2 laquelle ak = tk - tk-1 De manière proche de celle du calcul du paramètre TAS, le paramètre TAS(tx) est calculé, à partir des W-1 intervalles RR de la VFC dérivée d'une durée d'échantillonnage TE du signal ECG comprenant W battements cardiaques entre les deux battements cardiaques respectivement mesurés aux instants tx_vv+i et tx du signal ECG, de manière à sommer les valeurs ak si un segment entre deux valeurs ak et ak_i est descendant ; c'est-à-dire on accumule les valeurs ak dont chacune est inférieure à la valeur immédiatement précédente ak_i. La somme des valeurs a, correspondant aux segments descendants est ensuite divisée par la somme de toutes les W-1 valeurs ax_vv+2, aX-W+3 à aX. Le calcul du paramètre TAP(tx) suit l'équation (5) ci-dessous : X (tk - tk-1) si (a,i ak_i) TAP (tx) =100 k-X -W +2 X Eq. (5) (tk tk-1) k-X -W +2 De manière proche de celle du calcul du paramètre TAP, le paramètre TAP(tx) est calculé, à partir des W-1 intervalles RR de la VFC dérivée d'une durée d'échantillonnage TE comprenant W battements cardiaques entre les deux battements cardiaques respectivement mesurés aux instants tx_vv+i et tx du signal ECG, de manière à accumuler des valeurs ak si un segment entre deux valeurs ak et ak_i est montant ; c'est-à-dire on accumule des valeurs ak dont chacune est supérieure à la valeur immédiatement précédente ak_i. La somme des valeurs a, correspondant aux segments montants est ensuite divisée par la somme de toutes les W-1 valeurs aX-W+2, aX-W+3 aX. La figure 4 illustre les paramètres TAS(tx) et TAP(tx) obtenus en temps réel par un calcul basé sur 30 échantillons de la fenêtre glissante à chaque instant tx (W = 30). Comme illustré en figure 4, les états des systèmes sympathique et parasympathique de la personne 10 sont calculés en temps réel et montrés au fur et à mesure de la durée d'échantillonnage TE entre l'instant de départ à la 32ème seconde et l'instant de fin à la 316ème seconde. Selon un mode de réalisation avantageux mais optionnel, le module de traitement des données 222 effectue en outre le calcul des paramètres 10 désignés TPSP et TPPS. Le paramètre TPSP est configuré pour représenter le taux de pollution du système sympathique vers le système parasympathique qui indique probablement une mauvaise qualité de respiration de la personne 10 à cause 15 du stress perçu ou d'un effort physique Le calcul du paramètre TPSP est calculé selon l'équation (6) ci-dessous : N-1 (tn - tn_1) si (a. an_, et an+1 > a.) TPSP =100 /1-1 N-1 Eq. (6) (tn - tn_i) n=1 20 Le paramètre TPSP est calculé, à partir des N-1 intervalles RR de la VFC dérivée du signal ECG comprenant N battements cardiaques, de manière à accumuler des valeurs a, si un segment entre deux valeurs a, et an_i est descendant et entre deux segments adjacents montants dont l'un précède ledit segment et l'autre suit ledit segment. 25 Autrement dit, ce taux de pollution TPSP prend en compte un segment descendant qui perturbe la tendance monotone croissante, également désignée tendance croissante, formée par deux segments montants dont l'un précède ledit segment descendant et l'autre suit ledit segment descendant. 30 La somme des valeurs a, satisfaisant les conditions ci-dessus est ensuite divisée par la somme de toutes les valeurs al, a2 à awi. Le paramètre TPPS est configuré pour représenter le taux de pollution du système parasympathique vers le système sympathique qui indique probablement une mauvaise qualité de respiration de la personne 10 lors d'un repos ou de la récupération après un effort physique. Le calcul du paramètre TPPS est représenté dans l'équation (7) ci-dessous : N-1 n-tn_1) si (a. >ant et ant an_2 et an, >an) TPSP =100 n=1 N -1 Eq. (8) N-1 Eq. (9) 1 si (an >an_, et an_1 ak_, et a', > a, ) TPSP(tx) = 100 k-X -W +2 X Eq. (10) (tk tk-1) k-X -W +2 De manière proche de celle du calcul du paramètre TPSP, le paramètre 15 TPSP(Tx) est calculé, à partir des W-1 intervalles RR échantillonnés, de manière à accumuler des valeurs ak si un segment donné entre deux valeurs ak et ak_i est descendant et perturbe donc la tendance monotone croissante formée par les deux segments montants et adjacents audit segment donné dont l'un précède ledit segment descendant et l'autre suit immédiatement ledit 20 segment descendant. La somme des valeurs ak satisfaisant les conditions ci-dessus est ensuite divisée par la somme de toutes les W-1 valeurs ax-VV+2, aX-VV+3 à aX. 25 Le calcul du paramètre TPPS(tx) est représenté dans l'équation (11) ci- dessous : X (tk - tk_,) si (ak > ak_, et ak_, Le module d'affichage des données 130 est couplé, par une communication filaire ou sans fil tel que par ondes radio, au module de traitement des données 222. Ce module d'affichage des données 130 est configuré pour recevoir au moins une partie des résultats du procédé d'analyse effectué par le module de traitement des données 222 tels que les paramètres TAS, TAP, TPSP, TPPS, TAS(tx), TAP(tx), TPSP(tx), TPPS(tx) et en donner une représentation visuelle sur un écran d'affichage. Par exemple, les résultats des calculs des paramètres TAS et TAP comme illustrés en figures 3 et 4 sont visuellement affichés de manière graphique par le module d'affichage des données 130. L'invention n'est pas limitée à des moyens de communication utilisés pour la réalisation de la communication filaire ou sans fil entre le module d'affichage des données 130 et le module de traitement des données 222 ni à l'implémentation du module d'affichage des données 130. Des résultats des calculs des paramètres TAS, TAP, TPSP, TPPS, TAS(tx), TAP(tx), TPSP(tx), TPPS(tx) peuvent être affichés de manière graphique ou auditive. D'autres exemples de représentation visuelle affichables sur un écran sont donnés ci-dessous : - Objet (voiture, avion), personnage humain ou animal dont le comportement ou l'animation varie en fonction du ou des paramètres calculés. - Figures géométriques dont les dimensions et/ou la couleur varient en fonction du ou des paramètres calculés (par exemple cercles, ellipses, carrés, triangles). Les calculs des paramètres TAS, TAP, TPSP, TPPS, TAS(tx), TAP(tx), TPSP(tx), TPPS(tx) du procédé d'analyse selon l'invention, effectués par le module de traitement des données 222, sont des analyses temporelles à partir de la VFC dérivée d'un signal ECG. De plus, les paramètres tels que ci-dessus permettent de mesurer et d'indiquer le niveau énergétique du corps de la personne 10 en état de repos ou/et en activité par exemple pendant une activité sportive. Le procédé d'analyse temporelle de l'invention présente ainsi plusieurs avantages ci-dessous : 1) Surveillance des systèmes sympathique et parasympathique de manière indépendante et simultanée : Comme mentionné ci-dessus, le procédé de l'invention divise les valeurs al, a2, à awi des intervalles RR en deux groupes. Plus précisément, le paramètre TAS, reflétant l'activité du système sympathique, est calculé de manière à sommer seulement les valeurs a, correspondant aux segments descendants comme explicité ci-dessus alors que le paramètre TAP, reflétant l'activité du système parasympathique, est calculé de manière à sommer seulement les valeurs a, correspondant aux segments montants. Les deux branches parasympathique et sympathique peuvent être ainsi surveillées indépendamment en temps réel, ce qui permet à une personne ou son médecin de vérifier en temps réel le niveau de stress perçu (par exemple lorsque le paramètre TAS est supérieur au paramètre TAP) et de déterminer des remèdes adaptés et réalisés de préférence sur place pour diminuer la valeur du paramètre TAS. Le résultat des remèdes appliqués peuvent être également vérifiés sur place et en temps réel. Le médecin ou l'utilisateur peuvent alors rétroagir sur les remèdes ou actions à effectuer pour ramener l'activité du SNA à la valeur souhaitée ou ramener l'équilibre entre systèmes sympathique et parasympathique à l'équilibre souhaité. Les paramètres TPSP et TPPS reprennent ces caractéristiques et permettent donc également de montrer individuellement, durant la même période de surveillance, le niveau d'interférence entre les deux branches sympathique et parasympathique du SNA, ce qui permet ainsi à une personne ou son médecin de vérifier en temps réel la qualité de sa respiration. Plus la valeur des paramètres TPSP et TPPS est grande, moins la qualité de la respiration est bonne. Dans le cas d'une respiration normale et inconsciente, les valeurs de paramètres TPSP et TPPS sont d'environ de 10%. Dans le cas d'une respiration cohérente et symétrique, par exemple 5 secondes d'inspiration et 5 secondes d'expiration, les valeurs des paramètres TPSP et TPPS sont sensiblement égales à zéro. En revanche, les méthodes d'analyse fréquentielle appliquant des calculs tels que la FFT sur la VFC ne peuvent pas séparer les données concernant la branche sympathique de celles concernant la branche parasympathique. Autrement dit, les données concernant les systèmes sympathique et parasympathique sont mélangées et de manière aléatoire. La séparation des données de ces deux systèmes n'est alors pas possible. Une personne ou son médecin ne peut alors pas déterminer, à partir des calculs obtenus par les méthodes d'analyse fréquentielle, l'activité de la branche sympathique indépendamment de celle de la branche parasympathique et inversement. Si une personne inspire pendant 5 secondes et expire pendant les 5 secondes suivantes, on peut considérer qu'il respire de façon cohérente symétrique. La figure 5(a) montre un spectre de la VFC qui comporte un seul sommet à 0.1Hz et correspond à une période de respiration de 10 secondes (5 secondes d'inspiration et 5 secondes d'expiration). Ainsi, il n'est pas possible de séparer les données des deux branches sympathique et parasympathique à partir d'un tel spectre. En effet, la seule information valide obtenue à partir du spectre est la fréquence respiratoire dont le sommet est à 0.1Hz. Dans un deuxième exemple, une personne respire de manière inconsciente (aléatoire) dans un état de repos. La figure 5(b) montre un spectre de la VFC selon le deuxième exemple. On remarque qu'une distribution spectrale aléatoire change à chaque fois en raison du changement de l'activité réalisée par la personne et/ou du changement du rythme respiratoire de la personne. La distribution spectrale de la VFC représente donc la fréquence respiratoire comportant les données des deux branches du SNA déjà cumulées et mélangées de manière inséparable dans le domaine temporel. Il n'est alors plus possible de discriminer chaque composante à partir de la distribution spectrale de la VFC. Autrement dit, il est alors impossible d'extraire les informations souhaitées à partir de la distribution spectrale de la VFC comportant les données inséparablement mélangées des deux branches du SNA. 2) Suivi en temps réel : Les calculs de l'analyse temporelle tels que ci-dessus sont des implémentations mathématiques simples. Par rapport aux méthodes d'analyse fréquentielle impliquant des calculs exigeant une complexité considérablement plus élevée comme la FFT, la méthode CDM, etc, les calculs des paramètres TAS, TAP, TPSP, TPPS, TAS(tx), TAP(tx), TPSP(tx), TPPS(tx) consomment moins de temps d'exécution et permettent donc de surveiller en temps réel les activités des systèmes sympathique et parasympathique. Il est alors possible de rétroagir en temps réel pour contrôler l'activité du SNA. Dans un exemple de réalisation des méthodes d'analyse fréquentielle en temps réel où W (le nombre d'échantillons de la fenêtre glissante FFT) est 20 et le rythme cardiaque est de 60 battements/minutes, la précision spectrale est de 0.05Hz (Hertz) =-1 20 Hz, et il en résulte qu'il faut analyser et séparer un spectre qui s'étale entre 0 et 0.5Hz avec seulement 10 points de FFT. Ainsi, pour obtenir une précision raisonnable par exemple de 0.005Hz, W devrait être 200, avec les méthodes basées sur une analyse fréquentielle, il est nécessaire d'avoir les données de 200 battements cardiaques. La technique d'insertion des zéros, de l'anglais « zero padding » peut éventuellement augmenter la résolution spectrale mais n'augmente pas la précision de l'analyse fréquentielle. Il faut donc un temps d'exécution de plus de 3 minutes avec un rythme cardiaque moyen de par exemple 60 battements/minute pour avoir un premier résultat ; autrement dit, le résultat de calculs ne peut être facilement obtenu en temps réel. 3) Précision améliorée grâce au meilleur rapport signal sur bruit et une sensibilité réduite aux faux battements cardiaques. La sensibilité aux faux battements et/ou aux diverses arythmies du procédé d'analyse temporelle de l'invention est faible en raison d'une dépendance linéaire entre l'erreur et le résultat. Dans un exemple utilisant 300 battements cardiaques dont deux sont des faux battements, au niveau des calculs des paramètres TAS et TAP de l'invention, le taux d'erreur est d'environ -2 = 300 0.0066% et est donc négligeable. A l'inverse, en mettant en oeuvre l'analyse fréquentielle et en prenant le même exemple, deux faux battements conduisent à un effet bien supérieur comme montré par la suite en référence aux figures 6(a) à 6(d). La figure 6(a) montre une partie d'un premier signal ECG et sa VFC, le premier signal ECG comprenant 300 battements sans ajout de faux battements cardiaques. La figure 6(b) montre une partie d'un deuxième signal ECG et sa VFC, le deuxième signal ECG étant obtenu par ajout de deux faux battements respectivement aux instants t = 62 sec et t = 80.5 sec au premier signal ECG. Les figures 6(c) et 6(d) illustrent respectivement les spectres de VFC obtenues par réalisation d'une analyse fréquentielle conventionnelle de la VFC du premier signal ECG (sans faux battement) et ceux pour le deuxième signal ECG (avec deux faux battements). Le spectre montré en figure 6(d) est considérablement affecté par les bruits résultant des deux faux battements. Les deux faux battements pourraient éventuellement être enlevés manuellement après la fin de l'acquisition d'un signal ECG et avant l'analyse spectrale de la VFC, mais dans ce cas toute notion de temps réel disparaît. De plus, en un cas d'arythmie, les spectres de VFC seraient encore plus sérieusement affectés. Ainsi, il est clair que le procédé d'analyse temporelle de l'invention est, par rapport aux méthodes d'analyse fréquentielle, moins sensible aux bruits. 4) Procédé non-invasif et apte à être utilisé en milieu sportif Les calculs sont basés sur des intervalles RR et n'ont pas besoin d'autres types de données. Les intervalles RR peuvent être mesurés par un capteur ECG 110 de manière non invasive et dans un milieu imposant des contraintes tel qu'un milieu sportif. Les paramètres TAS, TAP, TPSP et TPPS de l'invention montrent ainsi efficacement l'activité du système sympathique et celle du système parasympathique et permettent alors aux chercheurs et médecins de disposer d'une représentation précise et fiable du SNA, ce qui leur permet ainsi d'accroitre la pertinence de leurs diagnostics et de réduire le temps et donc le coût d'acquisition et d'obtention de la représentation. Un utilisateur peut par ailleurs vérifier son état de stress en visualisant les résultats desdits paramètres, et peut donc les corréler avec l'état émotionnel de la personne et/ou les activités physiques que la personne a effectuées. Ces paramètres peuvent être appliqués à des calculs divers tels qu'un calcul du stress perçu et accumulé par une personne, un calcul montrant l'état émotionnel lié à une situation stressante, ou/et celui de l'énergie dépensée par un joueur lors d'une activité sportive, etc. Une des applications est l'estimation du stress accumulé pendant une longue période. Par exemple, si une personne respire de manière cohérente et symétrique (i.e. 5 secondes d'inspiration suivies immédiatement de 5 secondes d'expiration, constituant une cohérence cardiaque), le résultat des calculs des paramètres TAS et TAP montre que le système SNA de la personne 10 est à l'équilibre. C'est-à-dire le TAS et le TAP sont respectivement 50% et le ratio TAS/TAP est donc égal à 1. Néanmoins, si la personne subit un stress prolongé, par exemple plusieurs semaines, mois ou années, le résultat des calculs des paramètres TAS et TAP montre un déséquilibre (valeur du TAS supérieure à celle du TAP), et reflète cet état de stress. Un paramètre NS configuré pour représenter le niveau de stress de la personne 10 est représenté dans l'équation (12) : NS=100+ TAS- TAP Eq. (12) Cette estimation du niveau de stress peut être réalisée en temps réel, par exemple, à l'instant tx du signal ECG, en utilisant les paramètres TAS(TX) et TAP(TX) tels que ci-dessus. L'estimation du niveau de stress NS(TX) est représentée dans l'équation (13): NS(tx) = 100 + TAS(tx)- TAP(tx) Eq. (13) La figure 7(a) illustre les résultats d'un exemple d'un calcul du paramètre NS. Dans cet exemple, la personne 10 a subi un stress sur le long terme. Le 20 niveau de stress NS est égal à 126.3%, ce qui montre un dysfonctionnement du système SNA avec une augmentation de l'activité sympathique et une baisse de l'activité parasympathique. Dans cet exemple, le ratio TAP/TAS est de 1.68. En revanche, une valeur du paramètre NS inférieure à 100 (soit le ratio 25 TAS/TAP étant inférieur à 1) montre que la personne ne subit pas de stress spécifique ou bien l'impact du stress est très limité par exemple en raison d'une excellente capacité de gestion du stress ou d'un environnement peu stressant. Une valeur du paramètre NS égale à 100 montre que la personne n'a 30 pas de problèmes particuliers de stress. Une autre application de ces paramètres concerne l'estimation du stress résiduel d'une personne pendant une activité sportive. Un paramètre NSR 10 15 configuré pour représenter le niveau de stress résiduel de la personne 10 est représenté dans l'équation (14) : NSR = NS RCRepos Eq. (14) RC Cette estimation du niveau de stress résiduel peut être réalisée en temps réel, par exemple, à l'instant tx du signal ECG. L'estimation du niveau de stress résiduel NSR(TX) est représentée dans l'équation (15) : NSR(tx)= NS(tx) RCRePos Eq. (15) RC(tx) dans laquelle RCRepos est la fréquence cardiaque au repos, c'est-à-dire la personne 10 est inactive pendant au moins 20 secondes et de préférence 40 secondes, RC(tx) est la fréquence cardiaque à l'instant tx. La figure 7(b) illustre les résultats d'un exemple d'un calcul des paramètres NS(tx), NSR(tx) et RC(tx) pendant une phase de repos, une phase d'activité sportive et une phase de récupération. Les paramètres TAS, TAP, TPSP, TPPS, les applications telles que NS et NSR et leurs déclinaisons en temps réel, peuvent être utilisés dans des jeux vidéo et des simulateurs afin de montrer l'état du SNA et donc l'état émotionnel d'un joueur/patient. On peut ainsi détecter le niveau de stress, le niveau d'énergie de la personne. Il faut noter qu'au niveau de fiabilité de calcul d'un point de vue mathématique, c'est-à-dire dans une condition idéal où le signal ECG est bien capturé par le capteur ECG 110, les calculs des paramètres et des applications proposés par l'invention sont, par rapport à ceux des méthodes connues d'analyse fréquentielle, beaucoup plus précis. Le nombre (N ou W) de battements peut être réduit à un nombre compris entre 20 et 40. Une application possible consiste par exemple à prendre en compte l'activité de l'un ou plusieurs de ces paramètres relatifs au SNA pour modifier le comportement d'un avatar. Tel peut être le cas dans un jeu vidéo ou dans un programme médical ou de bien-être où l'avatar, personnage, animal ou objet animé verrait sa représentation et/ou son comportement et/ou ses facultés modifiés en fonction des paramètres relatifs au SNA. Dans le cadre d'un jeu vidéo, un objectif serait d'accroître le réalisme du jeu et d'augmenter la sensation du joueur. Dans le cadre d'un programme médical ou de bien-être, un objectif serait de contrôler son SNA afin de modifier la représentation et/ou le comportement et/ou les facultés de l'avatar de manière à atteindre un objectif. Par exemple, pour les enfants, un contrôle de la respiration impacterait en temps réel la représentation d'un avatar tel qu'un avion en vol, de sorte à stabiliser l'avion sur un écran, incitant de ce fait l'enfant à mieux maîtriser sa respiration. L'invention propose ainsi un système d'aide à la régulation du SNA. Ainsi, la représentation visuelle en temps réel de l'activité du SNA d'un sujet permet de rétroagir pour modifier son SNA en temps réel. En fonction des applications, cette rétroaction peut être : commandée par la personne elle-même, par exemple pour réduire son stress ; ou commandée par une autre personne ou un système externe afin de fournir au sujet des stimuli pour modifier son état émotionnel ou son niveau de stress. En alternative à une représentation visuelle, l'invention peut également fournir une représentation sonore ou tactile des paramètres calculés et relatifs au SNA. L'invention peut également cumuler une représentation visuelle ou sonore ou tactile des paramètres calculés et relatifs au SNA. L'invention n'est pas limitée aux seuls modes et exemples de réalisation décrits ci-dessus, mais s'étend à tous les modes de réalisation entrant dans la portée des revendications. En particulier, l'invention ne s'étend pas qu'aux procédés et systèmes permettant le calcul de chacun des paramètres TAS, TAP, TPSP et TPPS. Elle s'étend également aux procédés et systèmes permettant le calcul de l'un ou certains seulement des paramètres TAS, TAP, TPSP et TPPS.5 REVENDICATIONS1. Procédé de surveillance du système nerveux autonome (SNA) d'un sujet, le procédé comprenant : - l'acquisition d'au moins un signal d'électrocardiogramme (ECG), ce signal ECG comprenant W battements cardiaques, chaque battement cardiaque étant détecté à un instant t, avec t compris entre tx_vv+i et tx, l'instant tx étant le dernier battement des W battements cardiaques ; - la génération d'une donnée qui est fonction de la variabilité de la fréquence cardiaque (VFC) sur l'ensemble des W battements cardiaques, la VFC comprenant W-1 intervalles RR séparant deux battements cardiaques consécutifs détectés respectivement aux instants tk_i et tk, chaque intervalle RR présentant une durée d'une valeur ak = tk tk-1 avec k = (X-W+2)... X ; caractérisé en ce qu'il comprend : - le calcul à l'aide d'au moins un microprocesseur d'au moins un paramètre pris parmi : - le paramètre TAS(tx) représentant le niveau d'activité du système sympathique et qui est fonction du rapport entre d'une part une première somme obtenue en sommant, parmi les valeurs ak, uniquement celles qui sont supérieures à la valeur ak_i de l'intervalle immédiatement précédent et une deuxième somme obtenue en sommant toutes les W-1 valeurs ak; et/ou - le paramètre TAP(tx) représentant le niveau d'activité du système parasympathique et qui est fonction du rapport entre d'une part une première somme obtenue en sommant, parmi les valeurs ak, uniquement celles qui sont inférieures à la valeur ak_i de l'intervalle immédiatement précédent et une deuxième somme obtenue en sommant toutes les W-1 valeurs ak; - le paramètre TPSP(tx) représentant le taux de pollution du système sympathique vers le système parasympathique et qui est fonction du rapport 30 entre d'une part une première somme obtenue en sommant, parmi les valeurs ak, uniquement celles qui respectent la condition suivante :ak a et a > a k k-1 k-1 k-2 k ak et d'autre part une deuxième somme obtenue en sommant toutes les W-1 valeurs ak ; et/ou - le paramètre TPPS(tx) représentant le taux de pollution du système parasympathique vers le système sympathique et qui est fonction du rapport 5 entre d'une part une première somme obtenue en sommant, parmi les valeurs ak, uniquement celles qui respectent la condition suivante : (ak > ak_,, et ak_,, 2. Procédé selon la revendication précédente dans lequel ledit au moins un paramètre est le TAS(tx) représentant le niveau d'activité du système sympathique et dans lequel le TAS(tx) est calculé en appliquant l'équation 15 suivante : X (tk -tk-l) si (akG ak-l) TAS (tx) = 100 -W +2 (tk - tk_i) k-X -W +2 dans laquelle : = valeur de l'intervalle de temps immédiatement précédent l'intervalle de temps de durée ak, soit ak-1 = tk-1 tk-2. 20 k = (X-W+2), (X-W+3)... X 3. Procédé de surveillance selon l'une quelconque des revendications précédentes dans lequel ledit au moins un paramètre est le TAP(tx) représentant le niveau d'activité du système parasympathique et dans lequel 25 le TAP(tx) est calculé en appliquant l'équation suivante : X (tk - tk-l) si (aki ak-l) TAP (tx) = 100 k-X -W +2 (tk tk-1) k =X -W +2 X Xdans laquelle : - = valeur de l'intervalle de temps immédiatement précédent l'intervalle de temps de durée ak, soit ak-1 = tk-1 tk-2. - k = (X-W+2), (X-W+3)... X 4. Procédé de surveillance selon l'une quelconque des revendications précédentes dans lequel ledit au moins un paramètre est le TPSP(tx) et dans lequel le TPSP(tx) est calculé en appliquant l'équation suivante : (tk - tk_,) si (ak ak_, et a, > ak) TPSP(tX) ' P(tx) =100 k =X -W +2 (tk tk-1) k =X -W +2 dans laquelle : = valeur de l'intervalle de temps immédiatement précédent l'intervalle de temps de durée ak, soit ak-1 = tk-1 tk-2. k = (X-W+2), (X-W+3)... X 5. Procédé de surveillance selon l'une quelconque des revendications précédentes dans lequel ledit au moins un paramètre est le TPPS(tx) et dans lequel le TPPS(tx) est calculé en appliquant l'équation suivante : (tk - tk_,) si (ak > ak_, et ak_, 6. Procédé de surveillance selon l'une quelconque des revendications précédentes dans lequel le calcul d'au moins un paramètre comprend au moins le calcul des paramètres TAS(tx) et TAP(tx) et comprenant le calcul d'un paramètre NS(tx) représentatif d'un niveau de stress du sujet à l'instant tx et calculé en appliquant l'équation suivante : X XNS(tx) = 100 +TAS(tx)- TAP(tx) 7. Procédé de surveillance selon l'une quelconque des revendications précédentes comprenant le calcul d'un paramètre NSR(tx) relatif un niveau de stress résiduel du sujet à l'instant tx et calculé en appliquant l'équation suivante : NSR(tx)= NS(tx) RCRepos RC(tx ) dans laquelle RCRepos est la fréquence cardiaque au repos, c'est-à-dire le sujet est inactif pendant au moins 20 secondes et de préférence 40 secondes, RC(tx) est la fréquence cardiaque à l'instant tx 8. Procédé selon l'une quelconque des revendications précédentes dans lequel on calcule le paramètre TPSP représentant le taux approximé de pollution du système sympathique vers le système parasympathique et calculé en appliquant l'équation suivante : N-1 1 si (an an_2 et an, >an) TPSP =100 n=1 N -1 9. Procédé selon l'une quelconque des revendications précédentes dans lequel on calcule le paramètre TPPS représentant le taux approximé de pollution du système parasympathique vers le système sympathique et calculé en appliquant l'équation suivante : N-1 1 si (an >an_, et an_1 10. Procédé selon l'une quelconque des revendications précédentes comprenant une étape de fourniture d'une représentation visuelle et/ou auditive dudit au moins un paramètre. 11. Procédé selon la revendication précédente dans lequel ledit au moins paramètre est calculé en temps réel ou à intervalles réguliers et comprenant une étape de fourniture d'une représentation visuelle de l'évolution au cours du temps dudit au moins un paramètre. 12. Procédé de surveillance selon l'une quelconque des revendications précédentes dans lequel on calcule au moins deux paramètres, au moins l'un des paramètres étant pris parmi le TAS(tx) et le TPSP(tx) et au moins un autre paramètres étant pris parmi le TAP(tx) et le TPPS(tx) et comprenant une étape de fourniture d'une représentation visuelle desdits au moins deux paramètres, lesdits au moins deux paramètres étant représentés indépendamment l'un de l'autre. 13. Procédé selon la revendication précédente dans lequel lesdits au moins deux paramètres sont calculés en temps réel ou à intervalles réguliers et comprenant une étape de fourniture d'une représentation visuelle de l'évolution au cours du temps desdits au moins deux paramètres. 14. Procédé selon l'une quelconque des revendications précédentes dans lequel ledit au moins paramètre est calculé en temps réel et dans lequel la représentation visuelle de l'évolution au cours du temps dudit au moins un paramètre comprend l'affichage d'un graphique ou d'une barre évoluant au cours du temps. 15. Procédé selon l'une quelconque des revendications précédentes dans lequel ledit au moins un paramètre est calculé en temps réel et dans lequel la représentation visuelle de l'évolution au cours du temps dudit au moins un paramètre comprend l'affichage d'au moins un avatar et/ou d'au moins un objet animé dont l'animation ou l'évolution visuelle est fonction de l'évolution au cours du temps dudit au moins un paramètre. 16. Produit programme d'ordinateur comprenant des instructions, qui lorsqu'elles sont effectuées par au moins un processeur, exécute le procédé selon l'une quelconque des revendications précédentes. 17. Système de surveillance (100) du système nerveux autonome (SNA) d'un sujet, comprenant un dispositif configuré pour recevoir des données d'un signal d'électrocardiogramme (ECG), ce signal ECG comprenant W battements cardiaques, chaque battement cardiaque étant détecté à un instant t, avec t compris entre tx_W+i et tx, l'instant tx étant le dernier battement des W battements cardiaques; et comprenant au moins un module de traitement des données (222) équipé d'au moins un processeur configuré pour exécuter les étapes suivantes : - génération d'une donnée qui est fonction de la variabilité de la fréquence cardiaque (VFC) sur l'ensemble des W battements cardiaques, la VFC comprenant W-1 intervalles RR séparant deux battements cardiaques consécutifs détectés respectivement aux instants tk_i et tk, chaque intervalle RR présentant une durée d'une valeur ak = tk tk-1 avec k = (X-W+2)... X ; caractérisé en ce qu'il comprend : - calcul d'au moins un paramètre pris parmi : - le paramètre TAS(tx) représentant le niveau d'activité du système sympathique et qui est fonction du rapport entre d'une part une première somme obtenue en sommant, parmi les valeurs ak, uniquement celles qui sont supérieures à la valeur ak_i de l'intervalle immédiatement précédent et une deuxième somme obtenue en sommant toutes les W-1 valeurs ak; et/ou - le paramètre TAP(tx) représentant le niveau d'activité du système parasympathique et qui est fonction du rapport entre d'une part une première somme obtenue en sommant, parmi les valeurs ak, uniquement celles qui sont inférieures à la valeur ak_i de l'intervalle immédiatement précédent et une deuxième somme obtenue en sommant toutes les W-1 valeurs ak; - le paramètre TPSP(tx) représentant le taux de pollution du système sympathique vers le système parasympathique et qui est fonction du rapport entre d'une part une première somme obtenue en sommant, parmi les valeursak, uniquement celles qui respectent la condition suivante : ak ak-2 et ak > ak et d'autre part une deuxième somme obtenue en sommant toutes les W-1 valeurs ak ; et/ou - le paramètre TPPS(tx) représentant le taux de pollution du système 5 parasympathique vers le système sympathique et qui est fonction du rapport entre d'une part une première somme obtenue en sommant, parmi les valeurs ak, uniquement celles qui respectent la condition suivante : (ak > ak_i et ak_i 18. Système (100) selon la revendication précédente comprenant un détecteur (110) configuré pour acquérir le signal ECG et un dispositif d'affichage (130) configuré pour afficher une représentation visuelle dudit au moins un 15 paramètre. 19. Système (100) selon l'une quelconque des deux revendications précédentes configuré de manière à fournir en temps réel et audit sujet une représentation visuelle dudit au moins un paramètre. 20 20.Appareil d'aide à la régulation du SNA d'un sujet comprenant un système selon l'une quelconque des trois revendications précédentes. 21. Simulateur comprenant un système selon l'une quelconque des 25 revendications 17 à 20.