FR 3095546 A1 20201030 FR 1904422 A 20190426 Description Titre de l'invention : Dispositif et procédé d'analyse de l'état d'un système en contexte bruité Domaine de l'invention [0001] L'invention est dans le domaine de l'analyse de l'état d'un système, et en particulier, l'invention concerne un dispositif et un procédé permettant de compenser l'influence d'un bruit additif dans l'analyse de régularité d'un signal. Technique antérieure [0002] Dans un environnement donné, il peut être nécessaire de disposer d'informations sur l'état psycho-physiologique d'une personne ou l'état d'endommagement d'un équipement (structure métallique, pièce de moteur, pont, etc.). Pour ce faire, un ensemble de capteurs recueillent des informations de la personne ou de l'équipement et délivrent des signaux qui après analyse permettent d'attester et classifier l'état de la personne ou de l'équipement. [0003] A partir des signaux issus des capteurs, des signatures sont extraites. Ces dernières caractérisent les propriétés statistiques, temporelles ou fréquentielles des séries temporelles étudiées. [0004] Cependant, les signaux issus des capteurs sont entachés d'un bruit de mesure qui peut modifier l'information d'intérêt et la valeur des signatures extraites, et de fait entrainer des conséquences néfastes sur leur exploitation. [0005] Une signature particulière porte sur l'analyse temporelle du signal, dont notamment sa régularité notée 'a '. Cette dernière retranscrit une invariance de comportement du signal par changement d'échelle temporelle. Pour connaître la régularité d'un signal, Il est commun d'estimer le coefficient ou index de Hurst, noté 'H'. [0006] Or, les méthodes connues d'estimation du coefficient de Hurst, comme la méthode dite « Detrended Fluctuation Analysis » (DFA) ou la méthode dite « Detrending Moving Average » (DMA) fournissent une estimation erronée de la signature lorsque les signaux sont perturbés par un bruit. Un bruit peut par exemple être un bruit blanc ou un bruit additif. [0007] Pour pallier ce problème, une approche peut être de rehausser le signal dans le but. de supprimer le bruit, puis d'estimer la régularité sur le signal réhaussé. Cependant, les techniques de débruitage classiquement employées, par exemple en traitement de la parole qui reposent sur des approches de type Kalman, de soustraction spectrale ou qui se fondent sur les ondelettes ne donnent pas de résultats satisfaisants. En effet, l'étape de débruitage augmente la complexité calculatoire et le temps de calcul, ce qui est dommageable dans des applications de monitoring où les aspects temps réel sont es- 2 sentiels. De plus, cette étape de réhaussement (ou débruitage) altère la régularité du signal, en ce que la valeur estimée de cc sur le signal rehaussé s'en trouve donc faussée. Les techniques de réhaussement connues ne permettent pas de conserver la régularité du signal. [0008] Aussi, il est nécessaire de mettre en place des méthodes d'estimation de signatures qui soient robustes au bruit additif. La présente invention répond à ce besoin. Résumé de l'invention [0009] Un objet de la présente invention est un procédé d'estimation de la régularité d'un signal en contexte bruité. [0010] Un autre objet de la présente invention est un dispositif qui comprend des moyens permettant de mettre en oeuvre un procédé d'estimation de la régularité d'un signal en contexte bruité. [0011] L'invention trouvera des applications avantageuses dans tous les domaines dont par exemple le milieu (hio)médical où l'analyse de l'état psycho-physiologique d'une personne est requise, mais aussi dans tous les domaines industriels nécessitant le contrôle de pièces, que ce soit en aéronautique, automobile, ferroviaire, ou toute autre industrie lourde. [0012] Pour obtenir les résultats recherchés, il est revendiqué des procédés, dispositifs et produit programme d'ordinateur selon différents modes de réalisation. [0013] En particulier, il est proposé un procédé mis en oeuvre par ordinateur pour déterminer l'état d'un système. Le procédé comprend des étapes consistant à : - collecter des données relatives à un système, lesdites données étant des données bruitées comprenant des données utiles et du bruit ; - générer un signal à analyser à partir des données collectées, ledit signal étant un signal bruité comprenant un signal d'intérêt et du bruit ; - effectuer une analyse de régularité du signal bruité en compensant l'influence du bruit dans le calcul de la puissance du résidu du signal bruité; et - déterminer l'état dudit système en fonction du résultat de l'analyse de régularité du signal bruité. [0014] Selon des modes de réalisation du procédé, alternativement ou en combinaison : - l'étape d'analyse de régularité du signal bruité comprend les étapes de : - définir un jeu de valeurs de 'N' : - pour chaque valeur de 'N' : - calculer la valeur de la puissance du résidu du signal bruité relatif à la valeur 'N' à partir d'une fonction d'autocorrélation du signal bruité ; - calculer la valeur de la puissance du résidu du signal d'intérêt à partir d'une fonction d'autocorrélation du signal d'intérêt ; et 3 - estimer la régularité du signal bruité. - le procédé comprend après l'étape de générer un signal à analyser, une étape de calculer une fonction d'autocorrélation du signal bruité. - l'étape de calculer la valeur de la puissance du résidu du signal d'intérêt comprend les étapes de : - estimer les propriétés du bruit ; - calculer une fonction d'autocorrélation du bruit ; - calculer une fonction d'autocorrélation du signal d'intérêt à partir de la fonction d'autocorrélation du bruit et de la fonction d'autocorrélation du signal bruité ; et - compenser mathématiquement par les fonctions d'autocorrélation l'influence du bruit dans le calcul de la valeur de la puissance du résidu du signal d'intérêt. - l'étape d'estimation des propriétés du bruit consiste à estimer la stationnarité et/ou la nature colorée ou blanche du brui et/ou la modélisation paramétrique. - l'étape de définir un jeu de valeurs de 'N' consiste à découper le signal bruité en segments de taille 'N' et l'étape d'estimer la régularité du signal bruité consiste à estimer la tendance dudit signal pour le jeu de valeurs de 'N'. - l'étape de définir un jeu de valeurs de 'N' consiste à appliquer un filtre passe-bas d'ordre 'N' au signal bruité intégré et l'étape d'estimer la régularité du signal d'intérêt consiste à estimer la tendance dudit signal. - l'étape de collecter des données, consiste à collecter des données depuis des capteurs couplés au système. - les capteurs sont configurés pour relever des données relatives à la réponse électrique du cerveau (EEG) « ElectroEncephaloGramme » et/ou à la réponse électrique du coeur (ECG) « ElectroCardioGranune » et/ou à la réponse électrodermale et/ou à la pupillométrie. - les capteurs comprennent au moins un capteur inertiel qui comprend des accéléromètres et des gyromètres. - le procédé comprend de plus une étape consistant à afficher l'état dudit système en fonction du résultat de l'analyse de régularité. [0015] L'invention couvre aussi un produit programme d'ordinateur, ledit programme d'ordinateur comprenant des instructions de code permettant d'effectuer les étapes du procédé revendiqué, lorsque le programme est exécuté sur un ordinateur. [0016] L'invention couvre de plus un dispositif pour déterminer l'état d'un système, le dispositif comprenant : - des moyens pour collecter des données relatives à un système, lesdites données étant des données bruitées comprenant des données utiles et du bruit ; - des moyens pour générer un signal à analyser à partir des données collectées, ledit signal étant un signal bruité comprenant un signal d'intérêt et du bruit ; - des moyens pour effectuer une analyse de régularité du signal bruité en compensant l'influence du bruit dans le calcul de la puissance du résidu du signal bruité ; et - des moyens pour déterminer l'état dudit système en fonction du résultat de l'analyse de régularité. [0017] Selon un mode de réalisation, le dispositif comprend de plus une interface utilisateur configure° pour fournir l'état du système. Description des figures [0018] Différents aspects et avantages de l'invention vont apparaitre en appui de la des- cription (l'un mode préféré d'implémentation de l'invention mais non limitatif, avec référence aux figures ci-dessous, où : [0019] [fig.1] est un schéma bloc illustrant une implémentation du dispositif de l'invention selon un mode dc réalisation ; [0020] W8.21 illustre un enchainemcnt d'étapes du procédé de l'invention selon un mode de réalisation ; [0021] [fig.3] illustre un enchainemcnt d'étapes pour estimer la régularité d'un signal d'intérêt, selon un mode dc réalisation de l'invention ; [0022] [fig.4] illustre un enchainemcnt d'étapes pour compenser l'influence du bruit selon un mode de réalisation de l'invention ; [0023] [fig.51 est une représentation de la pente de la régression dc points obtenus à partir de différentes tendances d'un même signal ; [0024] [fig.6] illustre sur un graphe log(F(N)) en fonction dc log(N), l'évolution de dif- férentes valeurs de A 2(N) pour différentes valeurs dc N par la méthode DFA et par F. la méthode de l'invention clans un cas d'un signal non bruité et de ce même signal bruité; et 100251 Ifig.71 illustre des comparaisons de résultats obtenus par une approche de ré- haussement connue et par la méthode de l'invention dans l'estimation de la régularité de signaux synthétiques de régularité connue et maitrisée, qui ont été bruités. Description détaillée de l'invention [0026] La figure 1 montre une implémentation du dispositif (100) de l'invention, selon un mode de réalisation. Le dispositif (100) comprend généralement un module d'acquisition de données (102) configure pour recevoir des données d'au moins un capteur apte à enregistrer des informations sur un système surveillé (101), un module de génération de signal (104) configure pour générer un signal à analyser à partir des données acquises, un module (l'analyse (106) configure selon les principes de l'invention pour estimer la régularité d'un signal bruité et un module de résultat (108) couplé au module (l'analyse configure pour déterminer l'état d'un système à partir des résultats de l'analyse. Le système comprend de plus une interface utilisateur (110) configure° pour fournir le résultat de l'analyse. L'interface utilisateur peut avantageusement être configurée pour paramétrer au besoin le module d'acquisition des données. [0027] La désignation d'un système surveillé dans le cadre de la présente invention doit se comprendre comme pouvant être une personne ou un équipement à surveiller. [0028] La désignation d'au moins un capteur dans le cadre de la présente invention n'est pas limitative en ternies de nombre et de nature de capteurs. Ainsi, un ou une pluralité de capteurs peuvent être mis en oeuvre pour enregistrer des informations sur le système surveillé. [0029] Dans le cas où le système surveillé est une personne, dans le but de connaître par exemple son niveau d'attention, de somnolence, les capteurs peuvent être des électrodes installées sur la personne afin d'enregistrer des informations relatives à la réponse électrique du cerveau (EEG) « ElectroEncephaloGramme » et/ou à la réponse électrique du coeur (ECG) « ElectroCardioGramme » et/ou à la réponse électrodermale et/ou à la pupillométrie. [0030] Dans le cas où le système surveillé est un équipement, dans le but de connaître par exemple son niveau d'endommagement, les capteurs peuvent permettre d'étudier les forces qui lui sont appliquées et des vibrations qu'il subit. Dans un mode de réalisation, le capteur peut être un capteur inertiel de type « Incrtial Measurement Unit » (IMU) qui comprend des accéléromètres et des gyromètres, qui permettent aussi de suivre l'évolution de la position de l'objet dans l'espace. [0031] Le module d'acquisition de données (102) est configuré pour recevoir et enregistrer des données issues des capteurs. La transmission des données depuis les capteurs peut être soit filaire, soit être sans fil. Comme expliqué dans la section de l'état de l'art, les signaux issus des capteurs sont généralement entachés d'un bruit de mesure, bruit lié au contexte de mesure (environnement, capteurs). Les données acquises sont alors des données bruitées comprenant des données utiles et du bruit. [0032] Le module de génération de signal (104) est configuré pour générer une série temporelle à partir des données acquises. Selon des modes de réalisation, les données peuvent être formatées, filtrées, combinées afin de générer un signal à analyser. Les données acquises étant entachées de bruit de mesure additif, le signal à analyser qui est généré, est un signal bruité, composé du signal d'intérêt et du bruit additionnel. [0033] Le module d'analyse (106) est configuré pour estimer selon les principes de l'invention, la régularité du signal bruité. Avantageusement, le procédé de l'invention permet de compenser l'influence du bruit sur l'estimation de la régularité du signal bruité analysé. [0034] Le module de détermination de l'état du système (108) est couplé au module d'estimation de régularité pour déterminer selon le résultat fourni, l'état du système 6 surveillé. [0035] Le dispositif de l'invention comprend une interface utilisateur (110). L'interface uti- lisateur peut être configurée pour paramétrer le module d'acquisition de données (102). L'interface utilisateur peut être graphique pour visualiser les résultats obtenus par le module de détermination de l'état du système (108). L'interface utilisateur peut comprendre des moyens additionnels pour générer par exemple les résultats de manière audio. [0036] Le dispositif de la présente invention peut s' implémenter à partir d'éléments matériels et/ou logiciel. Le procédé opéré peut être disponible en tant que produit programme d'ordinateur sur un support lisible par ordinateur. Le support peut être électronique, magnétique, optique ou électromagnétique. Dans le cas où l'invention est implantée sur une machine de calcul reprogrammable (par exemple un circuit FPGA), le programme correspondant (c'est-à-dire la séquence d'instructions) peut être stocké dans ou sur un médium de stockage amovible (par exemple une carte SD, ou un stockage de masse tel que un disque dur ex. un SSD) ou non-amovible, volatile ou non-volatile, ce médium de stockage étant lisible partiellement ou totalement par un ordinateur ou un processeur. Le support lisible par ordinateur peut être transportable ou communicable ou mobile ou transmissible (i.c. par un réseau de télécommunication 2G, 3G, 4G, Wifi, BLE, fibre optique ou autre). [0037] Matériellement, le calculateur permettant d'opérer le procédé décrit, peut être im- plémenté sur tablette ou ordinateur portable. A titre d'exemple d'architecture matérielle adaptée à mettre en oeuvre l'invention, un dispositif peut comporter un bus de communication auquel sont reliés une unité centrale (le traitement ou microprocesseur (CPU, acronyme de « Central Processing Unit » en anglais), lequel processeur peut être "mufticore" ou "many-core"; une mémoire moite (ROM, acronyme de «Rend OnIy Memory » en anglais) pouvant comporter les programmes nécessaires à la mise en oeuvre de l'invention; une mémoire vive ou mémoire cache (RAM, acronyme de Random Access Memory » en anglais) comportant des registres adaptés à enregistrer des variables et paramètres créés et modifiés au cours de l'exécution des programmes précités ; et une interface de communication ou E/S (I/O acronyme de « Input/ouput » en anglais) adaptée à transmettre et à recevoir des données. [0038] La figure 2 illustre un enchainement d'étapes (200) pour déterminer l'état d'un système selon les principes de l'invention. Dans une première étape (202), le procédé permet de collecter une pluralité de données du ou des capteurs utilisés pour surveiller un système. Selon un mode de réalisation de l'invention, les données sont enregistrées par le module d'acquisition de données (102) du dispositif (100) de la figure 1. Les données acquises comprennent des données utiles et du bruit. [0039] Dans une étape suivante (204), le procédé permet de générer à partir de tout ou partie 7 des données acquises, un signal pour lequel la régularité doit être analysée. Les données collectées étant bruitées, le signal qui est généré est un signal bruité composé du signal d'intérêt et du bruit. [0040] L'homme de l'art tonnait différentes méthodes pour l'analyse de régularité d'un signal, qui se ramènent à deux grandes familles : - celle des estimateurs reposants sur l'analyse fréquentielle du signal où l'on trouve par exemple la méthode de Whittle locale, la méthode du périodogramme, la méthode fondée sur la transformée en ondelettes ou encore la méthode set-ni-paramétrique. Plus récemment, des solutions reposant sur « l'Empirical Mode Dccomposition » (EMD) ou la transformée de Fourier fractionnelle sont apparues. - celle des estimateurs reposants sur l'analyse temporelle du signal où l'on trouve par exemple l'analyse appelée « Rescaled Range », la méthode « Aggregated variance », la méthode de « rabsolute-value » ou encore la méthode de la variance du résidu. [0041] Le procédé de l'invention rentre dans la famille de l'analyse temporelle d'un signal, par l'estimation du coefficient de Hurst 'H'. Avant de détailler la méthode revendiquée, il est fait un rappel des approches connues afin d'estimer le coefficient de Hurst d'un processus mono-fractal. [0042] Une méthode nommée en anglais « Fluctuation Analysis » (FA) a été proposée au début des années 1990. Le principe est le suivant : après intégration d'un signal à analyser, conduisant à une nouvelle séquence `yin,', il est calculé pour différentes valeurs de '1' la grandeur suivante : [Math.1] FF-4(1) = 1 représente la moyenne temporelle. [0043] Puisque FFA( I) c, 1H avec représentant une relation de proportionnalité, log ( FFA (1)) peut être représenté en fonction de log (I ) afin d'estimer la valeur de 'H' qui est alors égale à la valeur de la pente de la droite de régression obtenue. [0044] Plus récemment, deux méthodes dérivées de la méthode (FA) sont largement utilisées pour l'estimation du coefficient de Hurst : la méthode dite « Detrended Fluctuation Analysis » (DFA) et la méthode dite « Detrending Moving Average » (DMA). [0045] La méthode DFA est décrite par C. K. Peng, S. V. Buldyrev, S. Havlin, M. Simons, H. E. Stanley et A. L. Goldberger dans «Mosaic organization of DNA nucleotides» Physical Review, n° %149, pp. 1685-1689, 1994. [0046] La méthode DMA est décrite par D. Osborne dans «Moving average detrending and the analysis of business cycles» Oxford Bull. Econom. Statist, n° %157, pp. 547-558, 1995. 8 [0047] Ces deux approches reposent sur un même principe qui est l'estimation de la tendance du signal intégré, relativement à une grandeur 'N. , en vue d'obtenir la puissance de résidu associé. Le résidu est défini comme la différence entre la version intégrée du signal et la tendance estimée. Si le principe est le même pour les deux méthodes, la manière d'estimer la tendance diffère selon l'approche, tout comme la définition de 'N'. [0048] Pour la méthode DFA, l'estimation de la tendance du signal intégré se fonde sur un découpage du signal en segments de taille dont on estime la tendance. Pour la méthode DMA, le paramètre 'N' correspond à l'ordre du filtre passe-bas qui est appliqué au signal intégré pour en obtenir la tendance. [0049] Les étapes de l'analyse temporelle d'un signal selon les approches DFA et DMA sont rappelées ci-après. [0050] En considérant échantillons successifs I y ( ) d'un signal , m 'y', chaque méthode comprend quatre étapes principales : [0051] Étape 1 : Le signal est intégré et sa moyenne lui est retranchée. Cette étape aboutit à un nouveau signal appelé 'le profil de y' et noté: [0052] [Math.2] Yint( v ) = M / = 1 eJ où [Math.2] v désigne la moyenne de y. [0053] Étape 2 : La tendance du profil est estimée. [0054] Dans le cas de la méthode DFA, le profil est découpé en 'I,' segments de taille 'N'. Puisque 'M' le nombre d'échantillons successifs n'est pas nécessairement un multiple de 'N', on ne prend en compte que L WN échantillons. La 1-ième tendance locale (1 = 1, ... , L) est estimée sur le segment I au sens des moindres carrés. La tendance globale est alors définie comme la concaténation des L tendances locales. Un exemple de tendances de profil est illustré sur la partie gauche de la figure 5 pour différentes valeurs de N= (10, 20, 40). [0055] Dans le cas de la méthode DMA, le signal est filtré par un filtre passe bas dépendant de 'N'. En effet, la réponse impulsionnelle du filtre est définie par 7'. = 1\-r pour ri = 0, ... N - 1. [0056] Étape 3 : La tendance globale obtenue est retranchée au profil, dans le but d'obtenir le résidu relatif à N. Unc grandeur F(1V) est ensuite calculée comme la racine carrée de la puissance du résidu. 9 [0057] Étape 4 : Les étapes 2 et 3 sont répétées pour différentes valeurs de N. Puisque FL ) ec Na avec a H + 1 dû à l'intégration du signal, log (F (N )) est représenté en fonction de log (N). La valeur de la pente de la droite de régression des points obtenus est alors une estimation de a. et donc de la régularité du signal. Une représentation d'une droite de régression obtenue à partir de trois tendances de profil est illustrée sur la partie droite de la figure 5. [0058] Ces deux approches DFA et DMA qui ont été largement appliquées dans le domaine biomédical pour l'étude de pathologies dans la voix, des signaux EEG et ECG, ont montré de bonnes performances dans l'estimation du coefficient de Hurst. Leur utilisation reste plus limitée dans d'autres domaines comme dans un contexte industriel. En météorologie, un exemple est la demande de brevet W02017/132225 Al «Weather-based industry analysis system». [0059] Il existe des variantes de ces approches dont le but est de leur apporter des correctifs. Pour la méthode DFA, on peut citer par exemple les DFA d'ordre supérieur à 1, tels que « l'Adaptative Fractal Analysis » (AFA) ou le « Regularized DFA », reposant sur une régularisation de de Tikhonov, ou citer J. W. Kantelhardt, E. Koscielny-Bundc, H. H. A. Rego, S. Havlin et A. Bundc dans «Detecting long-range corrélations with detrended fluctuation analysis,» Physica A: Statistical Mcchanics and its Applications, n° %1295, pp. 441-454, 2001 qui propose une étude qui vise à ajouter un terme correctif pour les petites valeurs de N. [0060] Pour la méthode DMA, on peut citer les variantes « Centered DMA », « Wcigted DMA d'ordre 1 » (WDMA-1), « Weigted Centered DMA d'ordre 1 » (WCDM A-1). [0061] Il ressort que dans les contextes d'application connus, les techniques de ré- haussement de signal ne permettent pas de conserver la régularité du signal. Avantageusement, le procédé de la présente invention ne procède pas à un réhaussement du signal pour supprimer le bruit, et ensuite estimer la régularité du signal d'intérêt, mais il estime directement la régularité du signal en prenant en compte l'influence du bruit additif de mesure, sans perdre l'information d'intérêt, ni altérer la régularité du signal. Le procédé de l'invention s'affranchit de l'étape de rehaussement des méthodes usuelles, de manière à compenser directement l'influence du bruit de mesure dans le calcul de la puissance du résidu du signal. [0062] Ainsi, revenant à la figure 2, l'analyse du signal bruité selon la présente invention consiste à estimer (206) la régularité du signal bruité avec une compensation directe de l'influence du bruit de mesure dans l'analyse. [0063] Dans une étape suivante (208), le procédé permet de déterminer l'état du système selon le résultat de l'analyse de régularité du signal bruité. [0064] La figure 3 qui détaille l'étape (206), illustre un enchainement d'étapes pour estimer la régularité d'un signal bruité avec compensation de l'influence du bruit, selon un mode de réalisation de l'invention basé sur la méthode DFA. Cependant, le procédé peut être transposé pour toute autre méthode d'estimation de régularité reposant sur un calcul de tendance du signal. [0065] Le procédé permet dans une première étape (302) de définir un jeu de valeurs de N, correspondant à un découpage du signal bruité en segments de taille 'N', puis va réitérer (304) plusieurs étapes (306 à 308) pour chaque valeur de N. [0066] Dans l'approche classique, F (N ) est défini selon l'équation suivante : [0067] [Math.3] rpîvs ( - [0068] avec x 1 ( j ) le 77-ieine échantillon de la l-ième tendance locale. [00691 11 est également possible d'exprimer F (N) de manière matricielle selon l'équation suivante: [00701 IMath.41 N 1. [00711 avec F7 (.1\T ) correspondant à la puissance du résidu relatif à la méthode choisie (DFA, DMA, etc.) et .UN, . une matrice appropriée. [0072] Dans l'étape (306), le procédé permet de calculer F2 jj,NT} , la valeur de la is puissance du résidu du signal bruité relatif à la valeur 'N' en cours. Cela consiste à exprimer F ( y ) en fonction des échantillons du signal bruité : z(*) = y (k) b(k) [0073] Etant donné le développement mathématique ci-dessus, il vient que la valeur de la puissance du résidu du signal bruité relatif à la valeur 'N' en cours, peut s'exprimer comme suit: [0074] [Math.5] 11 [0075] Cette expression peut être approximée selon la fonction d' autocorrélation du signal bruité (305). [0076] En effet, mathématiquement, il est possible d'écrire : [0077] [Math.6] [0078] avec EH l'opérateur espérance mathématique. De ce fait, dans le cas d'un processus stationnaire au sens large, cette équation devient : [0079] [Math.7] [0080] Finalement, pour un processus stationnaire et ergodique, il est possible d'approximer la valeur de la puissance du résidu du signal d'intérêt F.2 ( N ) par E[E2 ( N)]. L'homme de l'art comprend que cela revient à approximer le produit z(k)z(k+r) par l'autocorrélation instantanéeR z z (r) prise pour un retard , T1 avec r =0, , LN -1. Dans ce cas, il vient : [0081] [Math.8] qui peut être réécrit selon l'équation suivante, comme : [Math.9] [0082] Dans l'étape suivante (308), le procédé permet de calculer la valeur de la puissance du résidu du signal d'intérêt n 2 en compensant la contribution du bruit dans le F.(N) 12 calcul, selon la fonction d'autocorrélation du signal d'intérêt (307), comme détaillé par la description de la figure 4. [0083] Les étapes 306 à 308 sont répétées pour chaque valeur de N, et quand les N valeurs sont traitées, le procédé permet dans une étape suivante (310) d'estimer la régularité du signal bruité, puis de continuer avec l'étape de détermination de l'état du système (208). [0084] La figure 4 détaille l'étape (308) permettant de déduire la valeur de la puissance du résidu du signal d'intérêt A2 (TV) avec compensation de l'influence du bruit. Dans une étape (402), le procédé permet de d'estimer les propriétés du bruit, puis de déterminer (404) l'expression de la fonction d'autocorrélation du bruit, notée R h(r). [0085] L'homme du métier comprend que les propriétés du bruit à estimer peuvent porter par exemple sur la stationnarité et/ou la nature colorée ou blanche du bruit et/ou la modélisation paramétrique, etc. [0086] A partir de la fonction d'autocorrélation du bruit (404) et de la fonction d'autocorrélation du signal bruité (400), le procédé permet de déterminer (406) la fonction d'autocorrélation du signal d'intérêt. [0087] Dans le cas où le bruit et le signal sont non corrélés et stationnaires, la fonction d'autocorrélation du signal bruité Z s'exprime comme la somme de la fonction d'autocorrélation du signal d'intérêt Y et celle du bruit additif h selon l'équation suivante : [0088] [Matli.10] R y (r ) = Rzz(r ) - (r ) avec [Math.10] r = 0, , L N - 1 [0089] Puis dans une étape suivante (408), le procédé permet de compenser mathéma- tiquement l'influence du bruit dans le calcul de la valeur de la puissance du résidu du signal d'intérêt A2 F (N) [0090] Lc procédé continue à l'étape (310) pour estimer la régularité du signal. [0091] Un exemple est donné avec un bruit blanc : soit h un bruit blanc gaussicn de moyenne nulle et de variance Dans ce cas, le signal bruité -7 cst défini commeta suit : 100921 Z(n)= y(n)+ b(n). 13 [0093] La fonction d' autocorrélation de b étant égale à R ( ) = ot 5(r) avec la fonction 15 définie par : 100941 6(0_ {1 si r= 0 il vient : 0 sinon [0095] Ry,),(r)=Rz z(r)-01,-6(r)- 100961 Il est alors possible d'exprimer F2 \ la puissance du résidu associé à N du signal bruité, en fonction de celle du signal d'intérêt, selon l'équation : [0097] [Math.11] [0098] Et il est possible alors de déduire une estimation de la valeur du résidu dans le cas non bruité, selon l'équation suivante : [0099] [Math.121 [0100] avec Tt (r N, .) la tract de la matricerN, , ic la somme de ses éléments diagonaux. [0101] La figure 6 illustre sur un graphe log(F(N)) en fonction de log(N), l'évolution de dif- férentes valeurs de A 2 en fonction de N par la méthode DFA. Le signal a une ( ) régularité estimée à a = 1-60- En bruitant ce signal avec un bruit blanc de même puissance que le signal, on peut constater que le bruit modifie de manière significative la valeur de la puissance du résidu pour les petites valeurs de N. La valeur de a estimée dans le cas bruité en est donc altérée et est estimée à a noisy - 1.35. En utilisant la méthode de compensation du bruit de la présente invention, il est possible de supprimer cette contribution dans le calcul de A 2 , et la valeur de a estimée F.(is) est améliorée. Dans le cas décrit en exemple, elle vaut : a corn p e rt s ' = 1.61. [0102] La figure 7 illustre des comparaisons de résultats obtenus pour le graphe du haut par une approche classique de rehaussement de signal avant l'estimation de la régularité et pour le graphe du bas par la méthode de compensation du bruit de l'invention. L'exemple est montré pour 50 signaux synthétiques dont la régularité est connue, et allant de 0.1 à 0.9. Ces signaux sont ensuite bruités avec un bruit blanc dont la puissance est égale à celle du signal. 14 [0103] La courbe du haut montre les résultats quand les signaux sont d'abord rehaussés puis que la régularité est estimée sur le signal rehaussé ou débruité, par une approche classique. -11 peut alors être constaté que pour les hautes valeurs de a , le fait de rehausser un signal altère la valeur de régularité estimée. Ce phénomène s'accentue lorsque la régularité a priori du signal de hase augmente. [0104] La courbe du bas montre les résultats avec une compensation directe du bruit lors de l'estimation de la régularité selon la méthode de l'invention. Il peut alors être constaté que le phénomène d'altération de la valeur de régularité estimée est atténué, et que le procédé de l'invention permet ainsi de mieux estimer en moyenne la régularité d'un signal bruité. [0105] De plus, la variance d'estimation avec la méthode de compensation de bruit de l'invention est également plus faible que celle obtenue avec l'approche de rehaussement classique. Avantageusement, il peut alors être accordé plus de confiance à la valeur de régularité a estimée sur un signal réel par le procédé de l'invention que celle par les approches connues. Cette méthode permet ainsi une estimation de la régularité d'un signal plus robuste et plus précise. [Revendication 1] [Revendication 2] [Revendication 3] [Revendication 4] [Revendication 5] REVENDICATIONSProcédé (200) mis en oeuvre par ordinateur pour déterminer l'état d'un système, le procédé comprenant des étapes consistant à : -collecter (202) des données relatives à un système, lesdites données étant des données bruitées comprenant des données utiles et du bruit ; -générer (204) un signal à analyser à partir des données collectées, ledit signal étant un signal bruité comprenant un signal d'intérêt et du bruit ; -effectuer (206) une analyse de régularité du signal bruité en compensant l'influence du bruit dans le calcul de la puissance du résidu du signal bruité; et -déterminer (208) l'état dudit système en fonction du résultat de l'analyse de régularité du signal bruité. Le procédé selon la revendication 1 dans lequel l'étape (206) d'analyse de régularité du signal bruité comprend les étapes de : -définir (302) un jeu de valeurs de 'N' : -pour chaque valeur de 'N' (304): -calculer (306) la valeur de la puissance du résidu du signal bruité relatif à la valeur 'N' à partir d'une fonction d'autocorrélation du signal bruité (305) ; -calculer (308) la valeur de la puissance du résidu du signal d'intérêt à partir d'une fonction d'autocorrélation du signal d'intérêt (307); et -estimer (310) la régularité du signal bruité. Le procédé selon la revendication 1 ou 2 comprenant après l'étape (204) de générer un signal à analyser, une étape (400) de calculer une fonction d'autocon-élation du signal bruité. Le procédé selon la revendication 3 dans lequel l'étape (308) de calculer la valeur de la puissance du résidu du signal d'intérêt comprend les étapes de : -estimer (402) les propriétés du bruit ; -calculer (404) une fonction d' autocorrélation du bruit ; -calculer (406) une fonction d' autocorrélation du signal d'intérêt à partir de la fonction d'autocorrélation du bruit. et de la fonction d'autocorrélation du signal bruité ; et - compenser mathématiquement par les fonctions d'autocorrélation l'influence du bruit dans le calcul (408) de la valeur de la puissance du résidu du sienal d'intérêt. Le procédé selon la revendication 4, dans lequel l'étape d'estimation des 16 propriétés du bruit consiste à estimer la stationnarité et/ou la nature colorée ou blanche du brui et/ou la modélisation paramétrique. [Revendication 6] Le procédé selon l'une quelconque des revendications 2 à 5 dans lequel l'étape (302) de définir un jeu de valeurs de 'N' consiste à découper le signal bruité en segments de taille 'N' et l'étape (310) d'estimer la régularité du signal bruité consiste à estimer la tendance dudit signal pour le jeu de valeurs de 'N'. [Revendication 7] Le procédé selon l'une quelconque des revendications 2 à 5 dans lequel l'étape (302) de définir un jeu de valeurs de 'N' consiste à appliquer un filtre passe-bas d'ordre 'N' au signal bruité intégré et l'étape (310) d'estimer la régularité du signal d'intérêt consiste à estimer la tendance dudit signal. [Revendication 8] Le procédé selon l'une quelconque des revendications précédentes dans lequel l'étape (202) de collecter des données, consiste à collecter des données depuis des capteurs couplés au système. [Revendication 9] Le procédé selon la revendication 8 dans lequel les capteurs sont configurés pour relever des données relatives à la réponse électrique du cerveau (EEG) « ElectroEncephaloGramme » et/ou à la réponse électrique du coeur (ECG) « ElectroCardioGramme » et/ou à la réponse électrodermale et/ou à la pupillométric. [Revendication 10] Le procédé selon la revendication 8 ou 9 dans lequel les capteurs com- prennent au moins un capteur inertiel qui comprend des accéléromètres et des gyromètres. [Revendication 11] Le procédé selon l'une quelconque des revendications 1 à 10 comprenant de plus une étape consistant à afficher l'état dudit système en fonction du résultat de l'analyse de régularité. [Revendication 12] Produit programme d'ordinateur, ledit programme d'ordinateur comprenant des instructions de code non transitoire permettant d'effectuer les étapes du procédé selon l'une quelconque des revendications 1 à 11, lorsque ledit programme est exécuté sur un ordinateur. [Revendication 13] Dispositif (100) pour déterminer l'état d'un système, le dispositif comprenant : - des moyens (102) pour collecter des données relatives à un système (101), lesdites données étant des données bruitées comprenant des données utiles et du bruit ; - des moyens (104) pour générer un signal à analyser à partir des données collectées, ledit signal étant un signal bruité comprenant un signal d'intérêt et du bruit ; 17 [Revendication 14] - des moyens (106) pour effectuer une analyse de régularité du signal bruité en compensant l'influence du bruit dans le calcul de la puissance du résidu du signal bruité ; et [Revendication 15] - des moyens (108) pour déterminer l'état dudit système en fonction du résultat de l'analyse de régularité. Le dispositif selon la revendication 13 comprenant de plus une interface utilisateur (110) configurée pour fournir l'état du système. Le dispositif selon la revendication 13 ou 14 dans lequel le système est une personne ou un équipement.