Brevet FR3071084A1 — publié le 20190315

DOMAINE DE L'INVENTION L'invention concerne une méthode de configuration d'un classifieur servant à identifier le régime de fonctionnement d'un système, une méthode et un dispositif de détermination du régime de fonctionnement d'un système, en particulier dans le cas où le système est un système électrochimique tel qu'une pile à combustible. Par 'classifieur', on désigne ici un dispositif de calcul, par exemple un ordinateur, apte, à partir de données représentant un état et/ou un régime de fonctionnement d'un système, à identifier l'état ou le régime de fonctionnement dans lequel le système se trouvait lorsque lesdites données ont été acquises. ARRIERE-PLAN DE L'INVENTION Les systèmes électrochimiques, qui sont les systèmes pour lesquels l'invention a été initialement conçue, sont des systèmes comprenant un milieu dans lequel une réaction ou un phénomène chimique qui se produit est couplé à des échanges d'énergie électrique. Les systèmes électrochimiques incluent notamment les systèmes d'électrolyse ou d'électrodéposition, les piles électriques, les piles à combustibles, les accumulateurs et les supercondensateurs. La durée de vie et la fiabilité de ces systèmes sont importantes pour les utilisateurs de ces systèmes. Notamment, dans le cas des piles à combustible, la durabilité et la fiabilité des systèmes à pile à combustible sont l'un des enjeux majeurs des développements actuels, tant pour les applications de transport (au moins 5 000 h de fonctionnement en sollicitations dynamiques pour l'automobile) que pour les applications stationnaires (objectif : 40 000 h à 60 000 h de fonctionnement). Il existe donc un besoin pour des méthodes et des dispositifs de diagnostic permettant de déterminer le régime de fonctionnement de systèmes électrochimiques, dans le but d'identifier au plus tôt les défaillances affectant ces 30 systèmes. Différentes méthodes ont été développées pour répondre à ce besoin. Des méthodes de caractérisation électrochimique telles que la méthode par courbe de polarisation, la méthode d'interruption de courant et la Spectroscopie d'Impédance Electrochimique (SIE) ont ainsi été utilisées pour le 35 diagnostic des accumulateurs et convertisseurs d'énergie. Cependant, ces méthodes présentent des limitations pour les applications en temps réel. Elles nécessitent en effet des contrôles et des sollicitations précises sur le courant ou la tension, ainsi que des équipements supplémentaires souvent coûteux et volumineux. D'autres méthodes de diagnostic, parfois appelées « méthodes internes », 5 se fondent sur une modélisation du système électrochimique étudié. Cependant, cette modélisation nécessite un travail préalable important. De plus, souvent, la complexité des phénomènes physico-chimiques rend leur modélisation difficile, ce qui entraîne une faible précision ou fiabilité de ces méthodes. Enfin, des méthodes exploitant des outils de traitement du signal ont 10 également été utilisées pour extraire des indices sur l'état de santé de piles à combustible. Cependant, ces méthodes exploitent les signaux acquis par un nombre élevé de capteurs. Ces méthodes nécessitent donc le traitement de quantités d'informations considérables. Elles sont complexes à mettre en oeuvre et nécessitent des ressources en calcul importantes. 15 OBJET ET RESUME DE L'INVENTION Par suite, un premier objectif de l'invention est de remédier aux inconvénients des méthodes de détermination du régime de fonctionnement de systèmes électrochimiques connues et de proposer une méthode de 20 détermination du régime de fonctionnement de systèmes, notamment de systèmes électrochimiques, qui permette d'identifier le régime de fonctionnement du système parmi une pluralité de régimes de fonctionnement prédéterminés dudit système, qui soit relativement simple, tant en ce qui concerne l'acquisition des informations nécessaires au diagnostic qu'en ce qui concerne le traitement 25 de ces informations. Pour atteindre cet objectif, selon la présente invention il est tout d'abord proposé une méthode de détermination d'un régime de fonctionnement d'un système. Cette méthode peut notamment être utilisée comme méthode de diagnostic, pour détecter que le système est dans un mode de fonctionnement anormal ou du moins non souhaité, afin de déclencher les actions correctives nécessaires. Mais de plus, comme cette méthode de détermination fait appel à un classifieur pour identifier le régime de fonctionnement du système parmi un certain nombre de régimes de fonctionnement prédéterminés, selon l'invention, il est proposé une méthode de configuration du classifieur utilisé dans la méthode de diagnostic indiquée précédemment. Selon l'invention, la méthode de détermination d'un régime de fonctionnement d'un système comporte les étapes suivantes : - on acquiert un vecteur de valeurs à N instants successifs d'une variable dont la valeur dépend d'un régime de fonctionnement du système ; 5 - on calcule un vecteur d'enveloppes instantanées et un vecteur de fréquences instantanées à partir du vecteur de valeurs acquises ; - pour chacun des vecteurs de valeurs acquises, d'enveloppes instantanées, et de fréquences instantanées, on calcule un vecteur de valeurs de descripteurs statistiques ; 10 - on détermine alors le régime de fonctionnement dans lequel se trouvait le système lorsque le vecteur de valeurs acquises a été acquis, en mettant en oeuvre un classifieur, configuré pour déterminer un régime de fonctionnement du système parmi une pluralité de régimes de fonctionnement prédéterminés de celui-ci, en fonction desdits trois vecteurs de descripteurs statistiques calculés à 15 l'étape de calcul de vecteurs de descripteurs statistiques indiquée ci-dessus. Le classifieur peut notamment exécuter un algorithme de classification de type SVM (de l'anglais 'Support Vector Machine' ou en français machines à vecteurs de support ou séparateurs à vaste marge) de type multi-classes. La variable peut être une variable scalaire ou vectorielle. Il peut s'agir par 20 exemple d'une seule grandeur, par exemple une mesure de tension, de pression, d'intensité de courant électrique, ou autre. Avantageusement, la méthode permet donc d'identifier le régime de fonctionnement du système simplement à partir des valeurs d'un vecteur de valeurs acquises, c'est-à-dire à partir d'une série de valeurs de la variable 25 étudiée. Il s'est avéré en effet que sur la base d'un vecteur de descripteurs statistiques de la variable étudiée, d'un vecteur de descripteurs statistiques des enveloppes instantanées de cette variable, et d'un vecteur de descripteurs statistiques des fréquences instantanées de cette variable, il est possible, en 30 utilisant un classifieur configuré de manière approprié, d'identifier le régime de fonctionnement du système considéré avec une fiabilité tout-à-fait acceptable. La méthode proposée est particulièrement adaptée pour déterminer le régime de fonctionnement d'un système, notamment un système électrochimique comprenant une pile à combustible. La variable indiquée précédemment peut 35 alors être par exemple la tension mesurée aux bornes de la pile à combustible. La méthode peut être appliquée à des systèmes de toute nature, et notamment à des systèmes électrochimiques. Par exemple, le système peut être une partie d'un être vivant, dont on cherche à déterminer un dysfonctionnement éventuel. On peut par exemple 5 chercher à détecter un dysfonctionnement du coeur, ou du cerveau. Dans le cas du coeur, la variable mesurée peut par exemple être la fréquence d'un signal électrique représentant l'activité du coeur ou du cerveau (ECG ou EEG) Dans la méthode proposée, le calcul des vecteurs d'enveloppes instantanées et de fréquences instantanées peut se faire de différentes manières. 10 Dans un mode de mise en oeuvre, l'étape de calcul de vecteurs d'enveloppes instantanées et de fréquences instantanées comporte les étapes suivantes : - on calcule la transformée de Hilbert H(x(n)) du vecteur de valeurs acquises x(n). 15 - on calcule le signal analytique complexe Zx(n) à partir de la transformée de Hilbert H(x(n)) ainsi calculée, en utilisant la formule : Zx(n) = x(n) + jH(x(n)) - on calcule alors le vecteur d'enveloppes instantanées (ex(n)) et le vecteur de fréquences instantanées (fx(n)) à partir du signal analytique complexe 20 Zx(n) ainsi calculé. Une fois calculés le vecteur d'enveloppes instantanées et le vecteur de fréquences instantanées, le calcul de descripteurs statistiques, respectivement pour ces trois vecteurs, peut être réalisé. Les descripteurs statistiques utilisés dans le vecteur de valeurs acquises, le vecteur d'enveloppes instantanées et le 25 vecteur de fréquences instantanées ne sont pas nécessairement les mêmes mais ils peuvent être les mêmes. Dans un mode de mise en oeuvre, les descripteurs statistiques d'un des vecteurs de descripteurs statistiques calculés, et éventuellement de ces trois vecteurs, sont un ou plusieurs descripteurs statistiques choisis parmi la moyenne, 30 la variance, la médiane, le kurtosis, le maximum, le minimum, l'écart-type et le coefficient de variation. (Une définition de ces descripteurs est donnée plus loin). En particulier les descripteurs statistiques d'un des vecteurs de descripteurs statistiques calculés, et éventuellement de ces trois vecteurs, peuvent être précisément la moyenne, la variance, la médiane, le kurtosis, le 35 maximum, le minimum, l'écart-type et le coefficient de variation. Une méthode pour faire et pour optimiser le choix des descripteurs statistiques utilisés est indiquée plus loin. Une fois calculés les trois vecteurs de descripteurs statistiques respectivement pour le vecteur de valeurs acquises, le vecteur d'enveloppes 5 instantanées, et le vecteur de fréquences instantanées, le régime de fonctionnement du système peut être identifié à l'aide d'un classifieur approprié. Tout classifieur permettant, à partir des trois vecteurs de descripteurs statistiques respectivement pour le vecteur de valeurs acquises, le vecteur d'enveloppes instantanées, et le vecteur de fréquences instantanées, d'identifier 10 le régime de fonctionnement du système, peut être utilisé. Avantageusement, le classifieur en général n'est pas fondé sur un modèle, et d'autre part, il peut fournir sa réponse sans qu'il soit nécessaire de soumettre le système à une sollicitation particulière. Dans un mode de réalisation, le classifieur est un classifieur à 15 apprentissage automatique, et l'opération de configuration réalisée est une opération d'apprentissage automatique ou apprentissage machine. Une telle opération consiste à fixer les valeurs de différents paramètres utilisés par le programme informatique de calcul. Cette opération d'apprentissage machine est normalement réalisée à partir d'une base de données associant 20 d'une part, des valeurs de triplets vecteurs de descripteurs statistiques respectivement pour le vecteur de valeurs acquises, le vecteur d'enveloppes instantanées et le vecteur de fréquences instantanées au régime de fonctionnement dans lequel était le système lorsque les valeurs de la variable considérée ont été acquises. 25 Dans un mode de réalisation, le classifieur est un classifieur apte à mettre en oeuvre un algorithme de type 'Supports Vectors Machines' multiclasses, ou apte à mettre en oeuvre la méthode 'des k plus proches voisins', ou la méthode `Gaussian mixture model' (GMM). Les algorithmes de type 'Supports Vectors Machines' sont des algorithmes 30 fondés sur la recherche de l'hyperplan de marge optimale. Le principe de la maximisation de la marge est formulé de sorte à choisir un hyperplan optimal qui doit maximiser la distance entre la frontière de séparation et les points de chaque classe qui lui sont les plus proches. La méthode des `k plus proches voisins' consiste à étiqueter un nouvel 35 échantillon par l'étiquette majoritaire dans ses k-plus proches voisins. La méthode 'Gaussian Mixture Model', ou GMM, est une approche probabiliste, utilisée pour estimer les paramètres d'une distribution de variables aléatoires en les modélisant par une densité de mélange. Chaque composant du mélange caractérise une classe. 5 D'autre part, le classifieur est de préférence un classifieur configuré par la méthode de configuration de classifieur selon l'invention qui va maintenant être présentée. La méthode de configuration de classifieur selon l'invention est une méthode de 10 configuration d'un classifieur pour rendre celui-ci apte à déterminer un régime de fonctionnement d'un système parmi une pluralité de régimes de fonctionnement prédéterminés dudit système. Cette méthode comprend d'une part, au moins une fois pour chaque régime parmi ladite pluralité de régimes de fonctionnement prédéterminés, les 15 étapes suivantes : - (S30) on acquiert un vecteur de valeurs à N instants successifs d'une variable dont la valeur dépend d'un régime de fonctionnement du système ; - (S40-S50-S60) on calcule un vecteur d'enveloppes instantanées et un vecteur de fréquences instantanées à partir du vecteur de valeurs acquises ; 20 - (S70) pour chacun des vecteurs de valeurs acquises, d'enveloppes instantanées, et de fréquences instantanées, on calcule un vecteur de valeurs de descripteurs statistiques. D'autre part, la méthode comprend en outre l'étape suivante : - (S100) on configure un premier classifieur de telle sorte qu'il associe à chacun 25 des vecteurs de valeurs de descripteurs statistiques calculés à l'étape de calcul de vecteurs de descripteurs statistiques et définie ci-dessus le régime de fonctionnement dans lequel se trouvait le système lorsque le vecteur de valeurs acquises qui a permis de calculer le vecteur de descripteurs statistiques considéré a été acquis. 30 Lors de l'étape S30, les valeurs acquises sont acquises lorsque le système est dans le régime de fonctionnement considéré. Pour la mise en oeuvre des méthodes selon l'invention, les descripteurs statistiques utilisés peuvent être choisis de différentes manières. Par défaut, on peut utiliser un nombre de descripteurs statistiques 35 relativement grand. Par exemple, on peut éventuellement calculer des vecteurs de valeurs de descripteurs statistiques pour l'ensemble des descripteurs statistiques cités dans la table ci-dessous. 5 8 Paramètres statistique Equation Moyenne N x(n) x m = n=1 N Ecart type N n=1(x(n)-xm)2 0"= N Variance N Z n=1(x(n)-xm)2 x var - N Valeur efficace (RMS) N n=1 (x(n))2 xeff = N Amplitude maximale xmax= maxix(n)1 Amplitude minimale X min = minix(n)1 Biais N Ln x,,i,e_ (N-1)o-' Kurtosis ..i(x(n)-xm )4 Z N xkur - (N-1)u4 Facteur de crête CF CF x max - X eff Facteur de clairance CLF xmax CLF -i 7 Z_, N pF(7)1)2 \ n=1 Facteur de forme SF X eff = SF N - lx(n)1 N n=1 Facteur d'impulsion IF X rnax IF = N lx(n)I 171 n=1 Déviation absolue N N (x(n)-Xm )1 moyenne MAD MAD = I n=1 n=1 N MOMENT Central MC N N (x(n)-xm)P .n--1 MC = >n n=i N Etendue Xrange = I X max - X mini Médiane Coefficient de variation CV CV =- Z Cependant, il est préférable d'optimiser le choix des descripteurs statistiques. Dans ce but, parmi tous les descripteurs statistiques envisageables, on commence alors par identifier un sous-ensemble de descripteurs statistiques 5 particulièrement pertinents pour identifier un régime de fonctionnement du système considéré parmi les régimes de fonctionnement prédéterminés. Ce sous-ensemble comporte les descripteurs statistiques particulièrement pertinents pour le vecteur des variables acquises, les descripteurs statistiques particulièrement pertinents pour le vecteur des enveloppes instantanées des variables acquises, et 10 les descripteurs statistiques particulièrement pertinents pour le vecteur des fréquences instantanées des variables acquises. Ensuite, une fois que ce sous-ensemble de descripteurs statistiques particulièrement pertinents a été identifié, on utilise dès lors uniquement les descripteurs statistiques de ce sous-ensemble de descripteurs statistiques dans le 15 cadre de l'invention : on calcule des vecteurs de descripteurs statistiques uniquement pour ces sous-ensembles de descripteurs statistiques ; ensuite, on met en oeuvre comme classifieur de données un classifieur de données prévu pour déterminer le régime de fonctionnement du système parmi les régimes de fonctionnement prédéterminés du système à partir de valeurs issues seulement 20 du sous-ensemble de descripteurs statistiques ainsi identifié. Dans la définition ci-dessus, le « sous-ensemble de descripteurs statistiques particulièrement pertinents pour identifier un régime de fonctionnement d'un système parmi les régimes de fonctionnement prédéterminés » désigne un ensemble de descripteurs statistiques contenant 25 moins de descripteurs que le nombre initial de descripteurs envisagés et tel qu'il soit possible en fonction uniquement des descripteurs statistiques de ce sous-ensemble d'identifier le régime de fonctionnement du système considéré parmi ladite pluralité de régimes de fonctionnement prédéterminés, avec un taux d'erreur considéré comme acceptable. Cette opération d'identification peut être 30 réalisée en mettant en oeuvre un algorithme de type mRMR ou `Wrapper'. La méthode mRMR est une méthode de filtrage fondée sur le calcul de mesures statistiques classiques comme l'information mutuelle, la corrélation, etc..., qui permet de mesurer la dépendance statistique d'une variable par rapport à une autre. A partir de cette mesure, on cherche à satisfaire deux 35 critères : minimiser la redondance (mR) et maximiser la pertinence (MR) entre les variables. 10 La stratégie 'Wrapper' quant à elle cherche à optimiser un critère de performance optimisé par le calcul du score de la classification. Son principe est de générer des sous-ensembles de candidats et de les évaluer grâce à un algorithme de classification. 5 Dans un mode particulier de réalisation, les différentes étapes de la méthode de détermination de régime de fonctionnement d'un système ou de la méthode de configuration de classifieur selon l'invention sont déterminées par des instructions de programmes d'ordinateurs. En conséquence, l'invention vise aussi un programme d'ordinateur sur un 10 support d'informations, ce programme étant susceptible d'être mis en oeuvre dans un ordinateur, ce programme comportant des instructions de code de programme pour l'exécution des étapes de la méthode de détermination de régime de fonctionnement d'un système ou de la méthode de configuration de classifieur définie précédemment lorsque ledit programme est exécuté sur un 15 ordinateur. Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code source, code objet, ou de code intermédiaire entre code source et code objet, tel que dans une forme partiellement compilée, ou dans n'importe quelle autre forme souhaitable. Le programme selon l'invention 20 peut être en particulier téléchargé sur un réseau de type Internet. L'invention vise aussi un support d'informations non volatiles lisible par un ordinateur, sur lequel est enregistré un programme d'ordinateur tel que défini ci-dessus. Le support d'informations peut être n'importe quelle entité ou dispositif 25 capable de stocker le programme. Par exemple, le support peut comporter un moyen de stockage, tel qu'une ROM, par exemple un CD ROM ou une ROM de circuit microélectronique, ou encore un moyen d'enregistrement magnétique, par exemple une clef USB ou un disque dur. Alternativement, le support d'informations peut être un circuit intégré 30 dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour être utilisé dans l'exécution de la méthode de diagnostic ou de la méthode de configuration de classifieur décrites précédemment. L'invention vise aussi un classifieur configuré par la méthode de configuration de classifieur définie précédemment. 35 L'invention vise également un dispositif de détermination d'un régime de fonctionnement d'un système, le dispositif comportant 11 - un module d'acquisition, configuré pour acquérir des vecteurs de valeurs à N instants successifs d'une variable dont la valeur dépend d'un régime de fonctionnement du système ; - un module de calcul, configuré : 5 pour calculer un vecteur d'enveloppes instantanées et un vecteur de fréquences instantanées, à partir d'un vecteur de valeurs acquises ; et, pour chacun des vecteurs de valeurs acquises, d'enveloppes instantanées, et de fréquences instantanées, pour calculer un vecteur de valeurs de descripteurs statistiques ; 10 - un classifieur, configuré pour déterminer un régime de fonctionnement dans lequel se trouvait le système lorsqu'un vecteur de valeurs acquises a été acquis, parmi une pluralité de régimes de fonctionnement prédéterminés, en fonction des trois vecteurs de ladite pluralité de descripteurs statistiques calculés par le module de calcul respectivement pour les vecteurs de valeurs acquises, 15 d'enveloppes instantanées, et de fréquences instantanées. Le dispositif de détermination peut être intégré au sein d'un ordinateur. Cet ordinateur peut notamment être embarqué et constituer ainsi un système électronique embarqué. Il peut s'agir par exemple d'un système électronique embarqué comportant un microcontrôleur ou à base de microcontrôleur. 20 Par suite, l'invention vise enfin un système équipé d'un dispositif de détermination de régime de fonctionnement tel que défini précédemment. Ce système peut être notamment un système électrochimique. Il peut comprendre notamment une pile à combustible, ainsi que les organes auxiliaires de celle-ci. 25 BREVE DESCRIPTION DES DESSINS L'invention sera bien comprise et ses avantages apparaîtront mieux à la lecture de la description détaillée qui suit, de modes de réalisation représentés à titre d'exemples non limitatifs. La description se réfère aux dessins annexés, sur 30 lesquels : - la figure 1 est une vue schématique d'une pile à combustible équipée d'un dispositif de détermination du régime de fonctionnement de la pile conforme à l'invention ; - la figure 2 est un logigramme présentant de manière schématique les 35 étapes d'une méthode selon l'invention permettant de configurer un classifieur 12 pour rendre celui-ci apte à déterminer le régime de fonctionnement de la pile de la figure 1 ; - la figure 3 est un logigramme présentant de manière schématique les étapes d'une méthode mise en oeuvre par le dispositif de détermination de la pile 5 de la figure 1 et permettant, conformément à l'invention, de déterminer le régime de fonctionnement de la pile de la figure 1 ; et - la figure 4 est une vue schématique de l'architecture fonctionnelle du dispositif de détermination de la pile représentée sur la figure 1. 10 DESCRIPTION DETAILLEE DE L'INVENTION En faisant référence aux différentes figures, un dispositif 10 de détermination de régime de fonctionnement d'un système, en l'occurrence un système électrochimique, et un système électrochimique 100 comportant le dispositif 10, vont maintenant être décrits. 15 Le système électrochimique 100 est constitué principalement par une pile à combustible également référencée 100. Dans l'exemple présenté, une charge 102 est reliée en série avec les bornes A et B de la pile 100. La pile 100 alimente donc la charge 102 en électricité. Cette charge peut être un moteur électrique, une lampe électrique, etc. 20 De manière connue en soi, la pile 100 est alimentée en hydrogène par un réservoir 104 via une conduite 106. L'hydrogène non consommé est rejeté via une conduite 108. L'air, dont l'oxygène est utilisé comme comburant, pénètre dans la pile via une conduite 110. L'air appauvri en oxygène mais enrichi en eau est rejeté par une conduite 112. 25 La tension x aux bornes A et B de la pile 100 est mesurée par un voltmètre 120. Le voltmètre 120 transmet les valeurs de la tension x mesurées au dispositif 10. La pile 100 peut être exploitée dans différents régimes de 30 fonctionnement, en fonction des conditions dans lesquelles elle est mise en oeuvre. Le fonctionnement de la pile 100 dépend ainsi notamment : - du courant de charge de la pile en mode galvanostatique ; - de la température de la pile ; - des débits et pressions d'entrée et de sortie dans les conduites respectivement 35 106,110 et 108,112; 13 - du taux d'humidité (ou teneur en eau) des réactifs d'alimentation de la pile circulant dans les conduites 106 et 110. Ces paramètres sont les principaux paramètres de fonctionnement de la pile. A chaque régime de fonctionnement de la pile 100 correspondent donc 5 certaines plages de valeur pour les différents paramètres de fonctionnement de la pile. Le but principal de l'invention est de pouvoir déterminer, dans quel régime de fonctionnement se trouve (ou se trouvait) la pile, à partir de mesures simples. 10 Dans le cas présent, le dispositif 10 permet de déterminer le régime de fonctionnement de la pile simplement à partir des valeurs de tension x mesurées aux bornes de la pile. Les régimes de fonctionnement que le dispositif 10 va permettre d'identifier sont donc différents régimes de fonctionnement correspondant à différentes valeurs ou plages de valeur des paramètres de 15 fonctionnement de la pile. Sur le plan matériel, le dispositif 10 est un dispositif de calcul, en l'occurrence un ordinateur. Le dispositif 10 peut être disposé à proximité de la pile 100, ou il peut être placé à distance de celle-ci. Il peut par exemple être 20 constitué principalement par un ou plusieurs serveurs accessibles via internet ou autre, par liaison filaire ou sans fil. L'architecture matérielle du dispositif 10 est illustrée schématiquement sur la figure 1. Il comprend notamment un processeur 15, une mémoire vive 12, une mémoire morte 14, ainsi que des moyens de communication 16 avec le voltmètre 25 120. Ces composants sont éventuellement partagés avec d'autres éléments de contrôle de la pile 100. La mémoire morte 14 du dispositif 10 constitue un support d'enregistrement conforme à l'invention, lisible par le processeur 15 et sur lequel est enregistré un programme d'ordinateur C, conforme à l'invention, comportant 30 des instructions pour l'exécution des étapes d'une méthode de détermination de régime de fonctionnement de la pile conforme à l'invention. Ces étapes sont décrites ci-dessous en référence à la figure 3. Le programme C comporte notamment un programme de classification. Pour que le dispositif 10 puisse déterminer le régime de fonctionnement 35 de la pile 100, il faut d'abord configurer ce programme de classification. Le 14 programme C comporte donc également des instructions pour l'exécution d'une méthode de configuration de classifieur conforme à l'invention. Le programme de classification intégré au programme C doit être potentiellement capable, une fois configuré, d'identifier le régime de 5 fonctionnement de la pile parmi les différents régimes de fonctionnement possibles de la pile, en fonction des valeurs de trois vecteurs de descripteurs statistiques Vx, Vex et Vfx. Chacun de ces vecteurs Vx, Vex et Vfx comprend les descripteurs statistiques suivants : la valeur maximale (max), la valeur minimale (min), la 10 moyenne (p), l'écart-type (a), le coefficient de variation (CV), la médiane (Med), le biais (Skew), le kurtosis (Kur). Ces descripteurs statistiques sont obtenus à partir des données suivantes : - le vecteur Vx comprend les descripteurs statistiques calculés à partir du vecteur 15 des tensions x acquises ; - le vecteur Vex comprend les descripteurs statistiques calculés à partir du vecteur ex des enveloppes instantanées calculées à partir du vecteur des tensions x acquises ; et - le vecteur Vfx comprend les descripteurs statistiques calculés à partir du 20 vecteur fx des fréquences instantanées calculées à partir du vecteur des tensions x acquises. Le programme de classification, lorsqu'il est exécuté, exécute un algorithme de classification de type SVM multi-classes. Pour pouvoir être mis en oeuvre, ce programme de classification doit 25 préalablement être configuré. Sa configuration est réalisée de la manière suivante : Tout d'abord, on définit les différents régimes de fonctionnement que l'on cherche à identifier ; chacun de ces régimes de fonctionnement est défini par des valeurs ou des plages de valeurs des paramètres de fonctionnement de la pile 30 100. Dans l'exemple présenté, on identifie dix régimes de fonctionnement, indexés par un index R variant de 1 à 10. Ensuite, pour chacun de ces régimes de fonctionnement, on procède aux opérations suivantes : - Pour le premier des dix régimes, à l'étape S10, l'index R est initialisé à la valeur 35 1, et le système électrochimique est placé dans le premier des dix régimes de fonctionnement identifiés. Cela consiste donc à faire fonctionner la pile en 15 imposant que ses paramètres de fonctionnement aient des valeurs correspondant au premier régime de fonctionnement identifié. - Pour les régimes suivants, (pour chaque nouvelle valeur de R), dans une étape SlOa, on place le système électrochimique dans un nouveau régime parmi les dix 5 régimes de fonctionnement identifiés, et on incrémente la valeur de R. Pour chacun des dix régimes de fonctionnement, on réalise un certain nombre de fois, en l'occurrence cent fois de suite, une même séquence d'acquisition. Les séquences d'acquisition sont indexées par un index S. A l'étape S20, l'index S est initialisé à la valeur 1. Le fait d'obtenir un relativement grand 10 nombre de données permet de disposer d'une base de données riche, ce qui est important pour pouvoir configurer le programme de classification de manière fiable. Pour les itérations suivantes, (pour chaque nouvelle valeur de S), dans une étape S20a, on incrémente la valeur de S, tout en laissant le système 15 électrochimique fonctionner dans le régime de fonctionnement en cours. Pour chaque séquence d'acquisition, on réalise les opérations suivantes : Le système électrochimique reste placé dans l'un des dix régimes de fonctionnement identifiés. 20 Pendant que la pile est dans ce régime de fonctionnement, on acquiert un vecteur x(n) de 1000 valeurs de la tension x (la tension x est donc dans cet exemple la variable représentant l'état du système 100) à 1000 instants ou pas de temps successifs (étape S30). On calcule le vecteur ex(n) contenant les 1000 valeurs de l'enveloppe 25 instantanée ex à partir du vecteur x(n) acquis à l'étape S30. On calcule le vecteur fx(n) contenant les 1000 valeurs de la fréquence instantanée fi( à partir du vecteur x acquis à l'étape 530. Les vecteurs ex(n) des enveloppes instantanées et fx(n) des fréquences instantanées peuvent être obtenus de différentes manières. 30 Dans le cas présent, ils sont calculés de la manière suivante : Dans une étape S40, on calcule la transformée de Hilbert H(x(n)) du vecteur de tension x(n). La transformée de Hilbert mentionnée ici est naturellement la transformée de Hilbert discrète. Dans une étape S50, on calcule le signal analytique complexe Zx(n) à 35 partir de la transformée de Hilbert H(x(n)) ainsi calculée, en utilisant la formule : Zx(n) = x(n) + jH(x(n)) 16 Dans une étape S60, on calcule alors le vecteur d'enveloppes instantanées ex(n) et le vecteur de fréquences instantanées fx(n) à partir du signal analytique complexe Zx(n) ainsi calculé. L'enveloppe ou l'amplitude instantanée, qui est égale au module du signal 5 analytique Zx(n), est définie par la relation : ex(n) = I Zx(n) La fréquence instantanée fx, qui est la dérivée de la phase du signal analytique Zx(n), est définie par la relation : fx (n) = 1 d(arg(Zx ))(n) 2n dn 10 Grâce aux formules ci-dessus, pour chaque valeur de tension x acquise à l'instant n, on peut calculer une enveloppe instantanée ex et une fréquence fx, ce qui permet de calculer le vecteur ex(n) des enveloppes instantanées et le vecteur fx(n) des fréquences instantanées. 15 Les étapes S40, S50 et S60 permettent donc d'obtenir, pour chacun des vecteurs x(n) de tensions x acquis, le vecteur ex(n) d'enveloppes instantanées et le vecteur fx(n) de fréquences instantanées. Ainsi, pour chaque acquisition de valeurs de tension réalisée, on obtient un triplet de vecteurs (x(n),ex(n),fx(n)), qui correspond à un régime de 20 fonctionnement de la pile bien identifié (le régime dans lequel se trouvait la pile lorsque les valeurs de tension x(n) ont été acquises). Ensuite, dans une étape S70 : - on calcule le vecteur Vx des valeurs des descripteurs statistiques choisis pour le vecteur de tension x(n) ; 25 - on calcule le vecteur Vex des valeurs des descripteurs statistiques choisis pour le vecteur de fréquences instantanées ex(n) ; et - on calcule le vecteur Vfx des valeurs des descripteurs statistiques choisis pour le vecteur d'enveloppes instantanées fx(n). On obtient donc des triplets de vecteurs (Vx, Vex, Vfx). 30 On vérifie alors la valeur de l'indice S : Si S reste inférieur à 100 on poursuit ces séquences d'acquisitions, en incrémentant S (étape S20a). 17 Dans le cas contraire, on vérifie alors la valeur de l'indice R : Si R reste inférieur à 10 on poursuit les séquences d'acquisitions, en incrémentant R (étape SlOa). A l'issue de ces acquisitions, on obtient une base de données comprenant 5 tous les triplets de vecteurs (Vx,Vex,Vfx) ainsi calculés. Comme chacun de ces triplets est associé à un régime de fonctionnement connu, dans une étape 5100, cette base permet alors de configurer le classifieur 10 ou plus précisément, le programme de classification contenu dans le programme C de détermination de régime de fonctionnement. Grâce à cette 10 étape de configuration, à l'issue de l'étape 5100, on dispose d'un classifieur apte, sur la base de vecteurs Vx, Vex et Vfx de descripteurs statistiques calculés à partir d'une séquence d'acquisition de la tension de la pile, à déterminer avec une probabilité élevée le régime de fonctionnement de la pile parmi les dix régimes de fonctionnement initialement identifiés. 15 Les étapes 5110 et 5120 présentées sur la figure 2 sont optionnelles et visent à obtenir un classifieur plus performant. Dans ce but, par la méthode mRMR, on sélectionne dans une étape S110 un sous-ensemble de descripteurs statistiques particulièrement pertinents pour la tension, un sous-ensemble de descripteurs statistiques particulièrement 20 pertinents pour l'enveloppe instantanée, et un sous-ensemble de descripteurs statistiques particulièrement pertinents pour la fréquence instantanée. On charge dans le dispositif de calcul 10 un programme de classification apte, à partir d'une séquence d'acquisition de la tension de la pile, à déterminer avec un taux de reconnaissance ou un taux d'identification élevé le régime de 25 fonctionnement de la pile parmi les dix régimes de fonctionnement identifiés, sur la base seulement d'un vecteur contenant ce sous-ensemble de descripteurs statistiques. Dans une étape S120, on configure alors le programme de classification de telle sorte que pour chacun des dix régimes de fonctionnement choisis, et 30 pour chaque triplet de vecteurs de descripteurs statistiques calculé à l'étape S110, le programme associe les valeurs des descripteurs statistiques calculées à l'étape S110 au régime de fonctionnement de la pile au moment où ces valeurs ont été enregistrées (de manière analogue à ce qui est fait à l'étape S100). Le programme de classification qui est configuré à l'étape S100, et/ou 35 celui qui est configuré à l'étape 5120 peut être un classifieur de type SVM, ou un programme de classification contenant des instructions pour mettre en oeuvre la 18 méthode 'des k plus proches voisins', ou la méthode GMM : 'Gaussian mixture model'. Une fois ces opérations réalisées, le dispositif 10 constitue un classifieur 5 au sens de l'invention, apte à partir des valeurs des descripteurs statistiques choisis pour la tension, l'enveloppe instantanée et la fréquence instantanée, à identifier le régime de fonctionnement dans lequel se trouvait la pile lorsque les valeurs de tension (ayant servi à déterminer les valeurs d'enveloppe instantanée et de fréquence instantanée) ont été enregistrées. 10 Le dispositif 10 peut être alors utilisé pour réaliser un diagnostic de la pile 100 en réalisant les opérations suivantes, présentées sur la figure 3 : S30) on acquiert un vecteur contenant les N valeurs à N instants successifs de la tension aux bornes A et B de la pile 100 (N=1000). S40-S50-S60) on calcule un vecteur d'enveloppes instantanées (ex(n)) et 15 un vecteur de fréquences instantanées (fx(n)) à partir du vecteur x(n) de valeurs de tension ainsi acquis. S70) alors, pour chacun des vecteurs de valeurs de tension, d'enveloppes instantanées, et de fréquences instantanées (x(n),ex(n),fx(n)), on calcule des valeurs des descripteurs statistiques (Vx,Vex,Vfx) choisis précédemment. 20 Les étapes S30, S40-S50-S60 et S80 citées ci-dessus sont identiques aux étapes correspondantes réalisées pour configurer le classifieur du dispositif 10, présentées précédemment. Dans une étape S200, on identifie alors le régime de fonctionnement dans lequel se trouvait la pile lorsque le vecteur de tension ayant servi à calculer les 25 vecteurs d'enveloppe instantanée et de fréquence instantanée a été enregistrée, en utilisant le classifieur 10 configuré de la manière indiquée précédemment. De manière plus générale, la méthode de diagnostic peut être mise en oeuvre en utilisant, à l'étape S200 d'identification du régime de fonctionnement, n'importe quel classifieur apte à identifier le régime de fonctionnement de la pile 30 en fonction des valeurs des descripteurs statistiques choisis, déterminées à partir du vecteur de tension acquis pendant la période considérée et des vecteurs d'enveloppes instantanées et de fréquences instantanées calculés à partir de ce vecteur de tension. La structure fonctionnelle du dispositif 10, en tant que dispositif de 35 diagnostic de système électrochimique, va maintenant être présentée en relation avec la figure 4. 19 Le dispositif 10 comporte : - un module d'acquisition (M11), configuré pour acquérir des vecteurs de valeurs à N instants successifs d'une variable (x(n)) dont la valeur dépend d'un régime de fonctionnement du système ; 5 - un module de calcul (M12), configuré : pour calculer un vecteur d'enveloppes instantanées (ex(n)) et un vecteur de fréquences instantanées (fx(n)), à partir d'un vecteur x(n) de valeurs acquises ; et, pour chacun des vecteurs de valeurs acquises, d'enveloppes instantanées, 10 et de fréquences instantanées, pour calculer un vecteur de valeurs de descripteurs statistiques (Vx,Vex,Vfx) ; - un classifieur (M13), configuré pour déterminer un régime de fonctionnement dans lequel se trouvait le système électrochimique lorsqu'un vecteur de valeurs acquises (x(n)) a été acquis, parmi une pluralité de régimes de fonctionnement 15 prédéterminés, en fonction des trois vecteurs Vx, Vfx et Vrx des descripteurs statistiques calculés à l'étape S70 respectivement pour les vecteurs de valeurs acquises, d'enveloppes instantanées, et de fréquences instantanées. Dans le mode de réalisation décrit ici, les modules fonctionnels d'acquisition (M11), de calcul (M12), et le classifieur (M13) décrits précédemment 20 sont des modules logiciels implémentés par le dispositif 10 pour réaliser la configuration du classifieur M13 ou la détermination du régime de fonctionnement de la pile 100, selon le cas. La méthode présentée ci-dessus a été testée et a permis d'identifier des régimes de fonctionnement de la pile à combustible avec un taux de réussite très 25 élevé (98%). Le temps de calcul étant très court (de l'ordre de quelques millisecondes), la méthode peut être utilisée de manière prédictive. Le système 10 peut donc être utilisé comme système de diagnostic, pour identifier en temps réel des dysfonctionnements de la pile. 30 Ainsi, dans le cas où la pile se met à fonctionner dans un régime de fonctionnement qui pourrait conduire à son endommagement, dès que ce régime de fonctionnement a été identifié, on applique à la pile ou plutôt au système à pile à combustible, et notamment aux organes auxiliaires de la pile comme son système d'alimentation en air et en hydrogène, une commande entraînant une 35 modification du mode de fonctionnement de la pile, de manière à éviter que celle-ci ne soit endommagée. 20 Avantageusement, cette commande est appliquée avant même que les conséquences négatives ne surviennent, c'est-à-dire avant même que la pile ne soit endommagée. Quoique la présente invention ait été décrite en se référant à des 5 exemples de réalisation spécifiques, il est évident que différentes modifications et changements peuvent être effectués sur ces exemples sans sortir de la portée générale de l'invention telle que définie par les revendications. En outre, des caractéristiques individuelles des différents modes de réalisation évoqués peuvent être combinées dans des modes de réalisation additionnels. Par 10 conséquent, REVENDICATIONS1. Méthode de configuration d'un classifieur pour rendre celui-ci apte à déterminer un régime de fonctionnement d'un système parmi une pluralité de régimes de fonctionnement prédéterminés dudit système, la méthode comprenant d'une part, au moins une fois pour chaque régime parmi ladite pluralité de régimes de fonctionnement prédéterminés, les étapes suivantes : (S30) on acquiert un vecteur de valeurs à N instants successifs d'une variable (x) dont la valeur dépend d'un régime de fonctionnement du système ; (S40,S50,S60) on calcule un vecteur d'enveloppes instantanées (ex(n)) et un vecteur de fréquences instantanées (fx(n)) à partir du vecteur de valeurs acquises ; (S70) pour chacun des vecteurs de valeurs acquises, d'enveloppes instantanées, et de fréquences instantanées, on calcule un vecteur de valeurs de descripteurs statistiques ; la méthode comprenant en outre l'étape suivante : (S100) on configure un premier classifieur de telle sorte qu'il associe à chacun des vecteurs de valeurs de descripteurs statistiques calculés à l'étape 570, le régime de fonctionnement dans lequel se trouvait le système lorsque le vecteur de valeurs acquises qui a permis de calculer le vecteur de descripteurs statistiques considéré a été acquis. 2. Méthode de configuration d'un classifieur selon la revendication 1, dans laquelle le classifieur configuré à l'étape 5100 est un classifieur à apprentissage automatique, et l'étape 5100 est réalisée par apprentissage automatique. 3. Méthode de configuration d'un classifieur selon la revendication 1 ou 2, dans laquelle le classifieur configuré à l'étape S100 est un classifieur de type SVM, ou met en oeuvre la méthode des k plus proches voisins, ou la méthode Gaussian mixture model (GMM). 4. Méthode de détermination d'un régime de fonctionnement d'un système, comportant les étapes suivantes : 22 530) on acquiert un vecteur de valeurs à N instants successifs d'une variable (x) dont la valeur dépend d'un régime de fonctionnement du système ; 540-550-560) on calcule un vecteur d'enveloppes instantanées (ex(n)) et un vecteur de fréquences instantanées (fx(n)) à partir du vecteur de valeurs 5 acquises; 570) pour chacun des vecteurs de valeurs acquises, d'enveloppes instantanées, et de fréquences instantanées, on calcule un vecteur de descripteurs statistiques ; S200) on détermine alors le régime de fonctionnement dans lequel se 10 trouvait le système lorsque le vecteur de valeurs acquises (x(n)) a été acquis, en mettant en oeuvre un classifieur, configuré pour déterminer un régime de fonctionnement du système parmi une pluralité de régimes de fonctionnement prédéterminés de celui-ci, en fonction desdits trois vecteurs de descripteurs statistiques calculés respectivement pour les vecteurs de valeurs acquises, 15 d'enveloppes instantanées, et de fréquences instantanées. 5. Méthode de détermination selon la revendication 4, dans laquelle l'étape de calcul de vecteurs d'enveloppes instantanées et de fréquences instantanées comporte les étapes suivantes 20 (S40) on calcule la transformée de Hilbert H(x(n)) du vecteur de valeurs acquises x(n). (550) on calcule le signal analytique complexe Zx(n) à partir de la transformée de Hilbert H(x(n)) ainsi calculée, en utilisant la formule : Zx(n) = x(n) + jH(x(n)) 25 (S60) on calcule alors le vecteur d'enveloppes instantanées (ex(n)) et le vecteur de fréquences instantanées (fx(n)) à partir du signal analytique complexe Zx(n) ainsi calculé. 6. Méthode de détermination selon la revendication 4 ou 5, dans laquelle les 30 descripteurs statistiques d'un vecteur de descripteurs statistiques calculé à l'étape S70 sont la moyenne, la variance, la médiane, le kurtosis, le maximum, le minimum, l'écart-type, le coefficient de variation. 7. Méthode de détermination selon l'une quelconque des revendications 4 à 35 6, dans laquelle le classifieur mis en oeuvre à l'étape 5200 met en oeuvre un algorithme de type Supports Vectors Machines multiclasses. 23 8. Méthode de détermination selon l'une quelconque des revendications 4 à 7, dans laquelle le classifieur mis en oeuvre à l'étape S200 est un classifieur créé par la méthode selon l'une quelconque des revendications 1 à 3. 5 9. Méthode de détermination selon l'une quelconque des revendications 4 à 8, dans laquelle le système comprend une pile à combustible, et la variable est une tension (x) mesurée à des bornes (A,B) de la pile à combustible (100). 10 10. Programme d'ordinateur sur un support d'informations, ce programme étant susceptible d'être mis en oeuvre dans un ordinateur, ce programme comportant des instructions de code de programme pour l'exécution des étapes de la méthode de détermination selon l'une quelconque des revendications 4 à 9 ou de la méthode de configuration de classifieur selon l'une quelconque des 15 revendications 1 à 3 lorsque ledit programme est exécuté sur un ordinateur. 11. Support d'informations (14) non volatiles lisible par un ordinateur sur lequel est enregistré un programme d'ordinateur selon la revendication 10. 20 12. Dispositif de détermination (10) d'un régime de fonctionnement d'un système , le dispositif comportant - un module d'acquisition (M11), configuré pour acquérir des vecteurs de valeurs à N instants successifs d'une variable (x) dont la valeur dépend d'un régime de fonctionnement du système ; 25 - un module de calcul (M12), configuré : pour calculer un vecteur d'enveloppes instantanées (ex(n)) et un vecteur de fréquences instantanées (fx(n)), à partir d'un vecteur x(n) de valeurs acquises ; et, pour chacun des vecteurs de valeurs acquises, d'enveloppes instantanées, 30 et de fréquences instantanées, pour calculer un vecteur de valeurs de descripteurs statistiques (Vx,Vex,Vfx) ; - un classifieur (M13), configuré pour déterminer un régime de fonctionnement dans lequel se trouvait le système lorsqu'un vecteur de valeurs acquises (x(n)) a été acquis, parmi une pluralité de régimes de fonctionnement prédéterminés, en 35 fonction des trois vecteurs (Vx,Vex,Vfx) de ladite pluralité de descripteurs statistiques calculés par le module de calcul respectivement pour les vecteurs 24 (x(n)) de valeurs acquises, (ex(n)) d'enveloppes instantanées, et (fx(n)) de fréquences instantanées. 13. Système (100), notamment un système électrochimique, notamment 5 comprenant une pile à combustible, équipé d'un dispositif de détermination (10) selon la revendication 12.