FR 3029297 A1 20160603 FR 1461616 A 20141128 PROCEDE AUTOMATIQUE D'ESTIMATION DE L'ETAT DE CHARGE D'UNE CELLULE D'UNE BATTERIE [001] L'invention concerne un procédé automatique d'estimation de l'état de charge 5 d'une cellule d'une batterie ainsi qu'un support d'enregistrement et un système de gestion de batterie pour mettre en oeuvre ce procédé. L'invention a également pour objet un véhicule automobile comportant ce système de gestion de batterie. [002] Des procédés connus d'estimation de l'état de charge d'une cellule d'une batterie comportent : 10 a) l'acquisition, à un instant k, d'une valeur mesurée yk de la tension entre des bornes de la cellule, et d'une intensité mesurée ik du courant de charge ou de décharge de la cellule, b) l'estimation de l'état de charge SOCk de la cellule, à l'instant k, à partir de la valeur mesurée yk, de l'intensité mesurée ik et d'une capacité Cn,k3 de la cellule, cette 15 capacité Cn,k3, exprimée en Ah, représentant la quantité maximale d'énergie électrique stockable dans la cellule à l'instant k3, l'instant k3 étant l'instant auquel la capacité Cn,k3 a été estimée le plus proche de l'instant k, et c) l'estimation de la capacité Cn,k3 à un instant k3, les instants k3 étant moins fréquent que les instants k, cette estimation comportant : 20 - le calcul d'une prédiction de la capacité Cn,k3 à l'aide d'un modèle d'état qui relie la capacité Cn,k3 à la capacité Cn,k3-1 de cette même cellule à un instant précédent k3-1, puis - la correction de la prédiction de la capacité Cn,k3 en fonction d'une intensité mesurée, à l'instant le plus proche de l'instant k3, du courant de charge ou de 25 décharge de la cellule. [003] Un tel procédé est par exemple décrit dans le paragraphe 4.2 de la partie 3 de l'article suivant : L. Plett, et al. : « Extended Kalman filtering for battery management systems of LiPB-based HEV battery packs », journal of Power Sources, 2004, page 252-292. Par la suite, cet article est désigné par l'abréviation « Plett 30 2004 ». [004] L'invention vise à améliorer le procédé décrit dans le paragraphe 4.2 de Plett 2004 afin d'augmenter la précision de l'estimation de l'état de charge de la cellule tout en diminuant la puissance de calcul requise pour mettre en oeuvre ce procédé. Elle a donc pour objet un tel procédé d'estimation dans lequel la correction de la prédiction 35 de la capacité Cn,k3 comporte : - l'acquisition, à l'instant k3, d'une grandeur physique mesurable zk3 définie par la relation suivante : 302 92 97 2 où k est l'instant le plus proche de l'instant k3 et N est un entier supérieur ou égal à un, zk3 étant égale à ik-1 lorsque N est égal à un, - le calcul d'une prédiction 2k3 de la grandeur physique mesurable zk3 à l'aide du modèle d'observation suivant : 600 . C où : - k est l'instant le plus proche de l'instant k3, - SOCk_N est l'état de charge de la cellule à l'instant antérieur k-N, - Cn,k3 est, dans ce modèle d'observation, la prédiction de la capacité de la cellule à l'instant k3 calculée juste avant d'être corrigée, et - la correction de la prédiction de la capacité Cn,k3 en fonction de l'écart entre la grandeur physique acquise zk3 et la prédiction 2k3 calculée à partir du modèle d'observation. [5] Dans le procédé revendiqué, le modèle d'observation pour calculer la prédiction 2k3 est une fonction linéaire de la capacité Cn,k3. Par conséquent, le calcul de la prédiction 2k3 n'est pas entaché d'erreurs liées à la linéarisation du modèle d'observation autour de la valeur Cn,k3. A l'inverse, le modèle d'observation décrit au chapitre 4.2.2 de Plett 2004 est non linéaire et doit donc nécessairement être linéarisé autour de la capacité Cn,k3. Une telle linéarisation se traduit nécessairement par des erreurs d'estimation sur la capacité corrigée Cn,k3. De plus, cette opération de linéarisation du modèle d'observation consomme de la puissance de calcul. [6] Les modes de réalisation de ce procédé automatique d'estimation peuvent comporter une ou plusieurs des caractéristiques suivantes : - la correction de la prédiction de la capacité Cn,k3 est inhibée tant que la valeur d'un paramètre n'a pas franchi un premier seuil prédéterminé, ce paramètre étant choisi dans le groupe composé de la valeur mesurée yk, de l'état de charge estimé SOCk et d'une quantité Qk de courant débité par la cellule entre l'instant k et un instant précédent, et la correction de la prédiction de la capacité Cn,k3 est déclenchée en réponse au fait que ce paramètre descend en dessous du premier seuil prédéterminé si le paramètre est la valeur mesurée yk ou l'estimation de l'état de charge SOCk, ou monte au-dessus du premier seuil prédéterminé si ce paramètre est la quantité Qk de courant débitée ; - le procédé comprend le comptage du nombre d'itération de l'étape a) entre un instant t1 où la valeur mesurée yk ou l'état de charge estimé SOCk descend en dessous d'un second seuil prédéterminé et un instant t2 suivant où la valeur mesurée yk ou l'état de charge estimé SOCk descend en dessous du premier seuil prédéterminé, puis lors de la correction de la prédiction de la capacité Cn,k3 l'entier N est choisi égal au nombre d'itérations comptées entre les instants t1 et t2 et strictement supérieur à un ; - le premier seuil prédéterminé correspond à un état de charge de la cellule inférieur à 20 % et le second seuil prédéterminé correspond à un état de charge de la cellule supérieur à 90 %; - le procédé comprend le réglage de la covariance du bruit sur la mesure de l'intensité ik puis la correction de la prédiction de la capacité Cn,k3 en fonction de cette covariance du bruit sur la mesure de l'intensité ik ; - le modèle d'état utilisé lors du calcul de la prédiction de la capacité Cn,k3 est le suivant : Cn,k3 = Cn,k3-1, et le procédé comporte : - le calcul d'une prédiction de la covariance P3,k3 de l'erreur d'estimation de la capacité Cn,k3 en fonction d'une covariance Q3,0 du bruit de l'erreur de modélisation introduite par le modèle d'état utilisé, et - le réglage automatique de la covariance Q3,0 à l'aide de la relation suivante : Qd,k3 = 20 RyCn')/(3NceoiNs)]2, où : - y est une constante prédéterminée dont la valeur est strictement inférieure à un et supérieure à zéro, - C'' est la capacité maximale nominale de la cellule, - Nced est le nombre de cycles de décharge de la cellule avant que celle-ci atteigne sa 25 fin de vie, ce nombre étant prédéterminé, - Ne st le nombre d'itérations de l'étape c) dans un cycle de décharge de la cellule ; - les étapes a) et b) sont exécutées au moins N-fois entre les instants k3 et k3-1, où N est un nombre entier supérieur à trois ; - le procédé comporte également le calcul de l'état de santé SOHk3 de la cellule 30 à l'aide de la relation suivante : SOHk3 = Cn,k3/ Cf', où : - Cn,k3 est la capacité de la cellule estimée à l'instant k3 le plus récent, et - C'' est la capacité maximale nominale de la cellule. [007] Ces modes de réalisation du procédé automatique d'estimation présentent en outre les avantages suivants : 35 - Le fait de déclencher l'estimation de la capacité Cn,k3 automatiquement en réponse au franchissement d'un seuil permet d'adapter automatiquement la fréquence d'exécution de cette estimation à l'utilisation actuelle de la cellule. Cela permet notamment d'éviter d'estimer la capacité Cn,k3 lorsque cela n'est pas nécessaire et donc de limiter la puissance de calcul requise pour exécuter ce procédé sans pour 40 autant détériorer la précision de l'estimation de cette capacité. De plus, déclencher 302 92 97 4 l'estimation de la capacité Cn,k3 seulement lorsqu'un seuil d'état de charge ou un seuil de tension ou un seuil de quantité de charges débitées est franchi permet de réaliser une seule estimation par cycle de charge et de décharge de la cellule, ce qui est largement suffisant pour estimer avec précision cette capacité. 5 - Utiliser N mesures de l'intensité ik réalisée entre deux estimations successives de la capacité Cri,k3 pour estimer cette capacité permet d'augmenter la précision de cette estimation. - Utiliser des premier et second seuils prédéterminés qui délimitent un cycle de décharges de la cellule presque complet permet d'estimer la capacité Cn,k3 en tenant 10 compte de l'ensemble des valeurs mesurées durant ce cycle presque complet de décharge de la cellule. Cela augmente la précision de l'estimation de la capacité Cn,k3. - Dans Plett 2004, la prédiction Cri,k3 est corrigée en fonction de la covariance du bruit sur la mesure d'une grandeur dk3 qui n'est pas mesurée. Dans ces conditions, il est très difficile de déterminer la covariance du bruit de mesure de cette grandeur dk3. À 15 l'inverse, dans le procédé revendiqué, c'est la covariance du bruit sur la mesure de l'intensité ik qui est utilisée. Cette covariance du bruit de mesure de l'intensité ik est facile à déterminer expérimentalement par des tests réalisés sur l'ampèremètre utilisé pour mesurer l'intensité ik OU à partir de données fournies par le constructeur de l'ampèremètre. Cela rend donc la mise en oeuvre du procédé revendiqué beaucoup plus facile que la mise en oeuvre du procédé décrit dans Plett 2004. [8] L'invention a également pour objet un support d'enregistrement d'informations comportant des instructions pour l'exécution du procédé automatique d'estimation ci-dessus lorsque ces instructions sont exécutées par un calculateur électronique. [9] L'invention a également pour objet un système de gestion d'une batterie 25 équipé d'au moins une cellule, ce système comportant un calculateur électronique programmé pour : a) acquérir, à un instant k, une valeur mesurée yk de la tension entre des bornes de la cellule, et une intensité mesurée ik du courant de charge ou de décharge de la cellule, b) estimer l'état de charge SOCk de la cellule, à l'instant k, à partir de la valeur 30 mesurée yk, de l'intensité mesurée ik et d'une capacité Cn,k3 de la cellule, cette capacité Cn,k3, exprimée en Ah, représentant la quantité maximale d'énergie électrique stockable dans la cellule à l'instant k3, l'instant k3 étant l'instant auquel la capacité Cn,k3 a été estimée le plus proche de l'instant k, et c) estimer la capacité Cn,k3 à un instant k3, les instants k3 étant moins fréquent que 35 les instants k, cette estimation comportant : - le calcul d'une prédiction de la capacité Cn,k3 à l'aide d'un modèle d'état qui relie la capacité Cn,k3 à la capacité Cn,k3-1 de cette même cellule à un instant précédent k3-1, puis - la correction de la prédiction de la capacité Cn,k3 en fonction d'une intensité 40 mesurée, à l'instant le plus proche de l'instant k3, du courant de charge ou de décharge de la cellule, dans lequel le calculateur électronique est également programmé pour réaliser la correction de la prédiction de la capacité Cn,k3 en exécutant les étapes suivantes : - l'acquisition, à l'instant k3, d'une grandeur physique mesurable zk3 définie par la relation suivante : ni = k-1 où k est l'instant le plus proche de l'instant k3 et N est un entier supérieur ou égal à un, zk3 étant égale à ik-1 lorsque N est égal à un, - le calcul d'une prédiction 2k3 de la grandeur physique mesurable zk3 à l'aide du modèle d'observation suivant : ' ,k3 où : - k est l'instant le plus proche de l'instant k3, - SOCk_N est l'état de charge de la cellule à l'instant antérieur k-N, - Cn,k3 est, dans ce modèle d'observation, la prédiction de la capacité de la cellule à l'instant k3 calculée juste avant d'être corrigée, et - la correction de la prédiction de la capacité Cn,k3 en fonction de l'écart entre la 15 grandeur physique acquise zk3 et la prédiction 2k3 calculée à partir du modèle d'observation. [0010] Enfin, l'invention a également pour objet un véhicule automobile comportant : - au moins une roue motrice, - un moteur électrique apte à entraîner en rotation cette roue motrice pour déplacer le 20 véhicule automobile, - une batterie comportant au moins une cellule apte à stocker de l'énergie électrique et, en alternance, à restituer de l'énergie électrique pour alimenter le moteur électrique, cette cellule comportant deux bornes par l'intermédiaire desquelles elle est électriquement raccordée au moteur électrique, 25 - un voltmètre raccordé électriquement entre les bornes de la cellule pour mesurer la tension entre ces bornes, - un ampèremètre raccordé en série avec la cellule électrique pour mesurer l'intensité du courant de charge ou de décharge de cette cellule, et - le système revendiqué de gestion de la batterie raccordé au voltmètre et à 30 l'ampèremètre, ce système de gestion comportant un calculateur électronique programmable apte à estimer l'état de charge de la cellule de la batterie à partir des mesures du voltmètre et de l'ampèremètre. 302 92 9 7 6 [0011] L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif et faite en se référant aux dessins sur lesquels : - la figure 1 est une illustration schématique partielle d'un véhicule automobile équipé 5 d'une batterie électrique, - la figure 2 est une illustration schématique d'un modèle électrique d'une cellule de la batterie du véhicule de la figure 1 ; - la figure 3 est une illustration schématique d'un agencement d'estimateurs utilisé pour estimer l'état de charge d'une cellule de la batterie du véhicule de la figure 1, 10 - les figures 4 à 9 représentent des équations de différents modèles d'état et d'observation utilisés par les estimateurs de la figure 3 ; - la figure 10 est un organigramme d'un procédé d'estimation de l'état de charge d'une cellule à l'aide des estimateurs de la figure 3 ; - la figure 11 est un organigramme d'un procédé de détermination de l'état de charge 15 de la batterie du véhicule de la figure 1 ; - la figure 12 est un organigramme d'un procédé de planification des instants de rafraîchissement des estimations des états de charge de différentes cellules d'une batterie ; - la figure 13 est un chronogramme illustrant différents instants de 20 rafraîchissement planifiés à l'aide du procédé de la figure 12; - la figure 14 est une illustration schématique d'un autre agencement d'estimateurs utilisé pour estimer l'état de charge d'une cellule de la batterie du véhicule de la figure 1; - les figures 15 et 16 représentent, respectivement, un modèle d'état et un modèle 25 d'observation utilisé par les estimateurs de la figure 14; - la figure 17 est un organigramme d'un procédé d'estimation de l'état de charge d'une cellule à l'aide des estimateurs de la figure 14 ; - la figure 18 est une illustration d'un autre modèle d'état possible pour prédire la capacité et la résistance interne d'une cellule d'une batterie. 30 [0012] Dans ces figures, les mêmes références sont utilisées pour désigner les mêmes éléments. Dans la suite de cette description, les caractéristiques et fonctions bien connues de l'homme du métier ne sont pas décrites en détail. [0013] Dans cette description par « puissance de calcul », on désigne le nombre d'opérations à exécuter par un calculateur électronique. Ainsi, diminuer la puissance 35 de calcul signifie diminuer le nombre d'opérations à réaliser pour atteindre le même résultat ou un résultat de même nature. [0014] La figure 1 représente un véhicule automobile 2 à tractions électriques plus connu sous le terme de « véhicule électrique ». Les véhicules électriques sont bien connus et seuls les éléments structurels nécessaires pour comprendre la suite de 40 cette description sont présentés. Le véhicule 2 comporte : - un moteur électrique 4, apte à entraîner en rotation des roues motrices 6 pour faire rouler le véhicule 2 sur une chaussée 8, et - une batterie 10 qui alimente en énergie électrique le moteur 4. [0015] La batterie 10 comporte deux bornes 12, 14 de raccordement électrique et plusieurs cellules électriques raccordées électriquement entre ces bornes 12 et 14. Les bornes 12 et 14 sont raccordées aux charges électriques à alimenter. Ici, elles 5 sont donc notamment raccordées au moteur électrique 4. [0016] Pour simplifier la figure 1, seules quatre cellules électriques 18 à 21 sont représentées. Typiquement, ces cellules électriques sont regroupées en plusieurs étages et ces étages sont raccordés en série entre les bornes 12 et 14. Ici, seuls deux étages sont représentés. Le premier étage comporte les cellules 18 et 19, et le 10 second étage comporte les cellules 20 et 21. Chaque étage comporte plusieurs branches raccordées en parallèle. Chaque branche d'un étage comporte une cellule électrique ou plusieurs cellules électriques en série. Ici, le premier étage comporte deux branches, et chaque branche comporte une seule cellule électrique. Le deuxième étage est structurellement identique au premier étage dans l'exemple 15 représenté sur la figure 1. [0017] Ici, toutes les cellules de la batterie 10 sont structurellement identiques aux tolérances de fabrication près. Par conséquent, seule la cellule 18 est maintenant décrite plus en détail. [0018] La cellule 18 comporte deux bornes de raccordement électrique 30, 32 qui la 20 raccorde électriquement aux autres cellules et aux bornes 12 et 14 de la batterie 10. La cellule 18 est aussi fixée mécaniquement, sans aucun degré de liberté aux autres cellules de la batterie 10 pour former ce que l'on appelle fréquemment un « pack » de cellules. La cellule 18 est capable de stocker de l'énergie électrique lorsqu'elle n'est pas utilisée. Cette énergie électrique stockée est ensuite utilisée pour alimenter le 25 moteur 4, ce qui décharge la cellule 18. En alternance, la cellule 18 peut aussi recevoir de l'énergie électrique ce qui la charge. La décharge complète d'une cellule suivie de sa recharge complète constitue ce que l'on appelle un cycle de charge/décharge appelé simplement par la suite « cycle d'une cellule ». [0019] La cellule 18 est une cellule de type connu, par exemple, il s'agit d'une cellule 30 LiPB (Lithium-ion Polymer Battery) ou autre. [0020] La cellule 18 se caractérise par une capacité nominale initiale Cn'n', une résistance interne initiale RO'', une intensité Imax de courant débité, une tension maximale Umax, une tension minimale Um,n et une fonction OCV(SOCk). La capacité Cn'n' est la capacité initiale de la cellule 18. La capacité d'une cellule représente la 35 quantité maximale d'énergie électrique qui peut être stockée dans cette cellule. Cette capacité est exprimée en Ah. Au fur et à mesure que la cellule 18 vieillit, c'est-à-dire au fur et à mesure que le nombre de cycles de charge et de décharge augmente, la capacité de la cellule diminue. À l'instant k, la capacité nominale de la cellule 18 est notée Cn,k par la suite. 40 [0021] La résistance interne initiale RO'' est la valeur de la résistance interne de la cellule 18 avant que celle-ci ne commence à vieillir. La résistance interne d'une cellule est une grandeur physique que l'on retrouve dans la plupart des modèles électriques d'une cellule électrique. Lorsque la cellule vieillit, typiquement, la résistance interne augmente. À l'instant k, la résistance interne de la cellule 18 est notée ROk. [0022] Imax est l'intensité maximale du courant qui peut être délivrée par la cellule 18 5 sans l'endommager. [0023] Umax est la tension maximale qui peut être présente en permanence entre les bornes 30 et 32 de la cellule sans l'endommager. [0024] La tension Um,, est la tension minimale entre les bornes 30 et 32 lorsque la cellule 18 est complètement déchargée. Par la suite, on considère que Imax, Umax, Um,n 10 sont des grandeurs physiques constantes qui ne varient pas au cours du temps. [0025] OCV(SOCk) est une fonction prédéterminée qui retourne la tension à vide de la cellule 18 en fonction de son état de charge SOCk. La tension à vide est la tension mesurable entre les bornes 30 et 32 après que la cellule 18 ait été isolée électriquement de toute charge électrique pendant deux heures. 15 [0026] L'état de charge à l'instant k de la cellule 18 est noté SOCk. L'état de charge représente le taux de remplissage de la cellule 18. Il est égal à 100 % lorsque la quantité d'énergie électrique stockée dans la cellule 18 est égale à sa capacité Cn,k. Il est égal à 0 % lorsque la quantité d'énergie stockée dans la cellule 18 est nulle, c'est-à-dire qu'il ne peut plus être extrait d'énergie électrique de la cellule 18 pour alimenter 20 une charge électrique. [0027] Les paramètres Cn'n', RO, Imax, Umax, Um'-, et la fonction OCV(SOCk) sont des paramètres connus de la cellule. Par exemple, ils sont donnés par le constructeur de la cellule ou sont déterminés expérimentalement à partir de mesures réalisées sur cette cellule. 25 [0028] La batterie 10 comporte également pour chaque cellule : - un voltmètre qui mesure la tension entre les bornes de cette cellule, et - un ampèremètre qui mesure l'intensité du courant de charge ou de décharge de cette cellule. [0029] Pour simplifier la figure 1, seuls un voltmètre 34 et un ampèremètre 36 de la 30 cellule 18 ont été représentés. [0030] Contrairement aux différents paramètres de la cellule 18 introduits précédemment, l'état de charge SOCk de la cellule 18 n'est pas mesurable. Il doit donc être estimé. À cet effet, le véhicule 2 comporte un système 40 de gestion de la batterie 10 plus connu sous l'acronyme BMS (« Battery Management System »). Ce 35 système 40 a notamment pour fonction de déterminer l'état de charge de la batterie 10 ainsi que l'état de santé de cette batterie. Pour déterminer cet état de charge et cet état de santé, le système 40 est capable d'estimer l'état de charge et l'état de santé de chaque cellule de la batterie 10. L'état de santé d'une cellule représente l'état d'avancement du vieillissement de cette cellule. Ici, l'état de santé d'une cellule, 40 à l'instant k, est noté SOHk. Par la suite, il est mesuré par le rapport C,,k/Cn'n'. Pour calculer l'état de santé d'une cellule, le système 40 est donc aussi capable d'estimer la capacité Cri,k de cette cellule à l'instant courant k. 302 92 9 7 9 [0031] Pour réaliser ces différentes estimations, le système 40 est électriquement raccordé à chaque voltmètre et chaque ampèremètre de la batterie 10 pour acquérir les mesures de la tension et de l'intensité du courant entre les bornes de chaque cellule. [0032] Ici, le système 40 comporte une mémoire 42 et un calculateur électronique 44 programmable, apte à exécuter des instructions enregistrées dans la mémoire 42. À cet effet, la mémoire 42 comporte les instructions nécessaires pour l'exécution des procédés des figures 10 à 12 et/ou de la figure 17. Cette mémoire 42 comporte aussi les valeurs initiales des différents paramètres nécessaires à l'exécution de ces procédés. La structure du système 40 est donc identiques ou similaire à celles des systèmes connus de gestion d'une batterie et n'est pas décrite plus en détail. [0033] La figure 2 représente un modèle électrique 50 de la cellule 18. Ce modèle est connu sous le terme de « Modèle de Thévenin du premier ordre » ou « Lumped parameter model ». Il comporte successivement raccordé en série en partant de la 15 borne 32 jusqu'à la borne 30: - un générateur 52 de la tension à vide OCV(SOCk), - un circuit RC parallèle 54, et - une résistance interne 56 appelée par la suite, à l'instant k, « résistance interne ROk ». 20 [0034] Le circuit 54 comporte un condensateur de capacité CD raccordé en parallèle avec une résistance de valeur RD. Par la suite, on considère que ces deux paramètres CD et RD du modèle 50 sont connus et constants au cours du temps. La tension à l'instant k aux bornes du circuit 54 est notée VD,k. La valeur à l'instant k de la tension entre les bornes 30 et 32 de la cellule 18 est notée yk et l'intensité, au même 25 instant, du courant de charge ou de décharge de la cellule 18 est notée ik. [0035] La figure 3 représente un premier mode de réalisation d'un agencement d'estimateurs implémenté dans le système 40 pour estimer l'état de charge et l'état de santé de la cellule 18. Chaque estimateur est implémenté sous la forme d'un algorithme d'estimation exécuté par le calculateur. Ainsi, on parlera pas la suite aussi 30 bien « d'exécution d'un estimateur » que « d'exécution d'un algorithme d'estimation ». Dans ce premier mode de réalisation, le système 40 comporte un estimateur 60 de l'état de charge SOCk et de la tension VD,k à partir de la valeur mesurée yk de la tension et de l'intensité mesurée ik. L'estimateur 60 est ici implémenté sous la forme d'un filtre de Kalman. Il utilise donc un modèle d'état 62 (figure 4) et un modèle 35 d'observation 64 (figure 5). Dans ces figures 4 et 5, les équations de ces modèles sont représentées en utilisant les notations précédemment définies. Les notations R0k2 et Cn,k3 représentent, respectivement, la capacité et la résistance interne de la cellule 18, respectivement, à des instants k2 et k3. Ces instants k2 et k3 sont définis par la suite. De plus, dans le modèle 62, xk désigne le vecteur d'état [sock, vp,k]r à 40 l'instant k. Dans cette description, le symbole « T » désigne l'opération transposée mathématique. L'opération multiplication est représentée par l'opérateur « . » ou « * ». [0036] Par la suite, on considère que l'origine des temps correspond à la valeur nulle de l'instant k. Dans ces conditions, l'instant k courant est égal à kTe, où T, est la période d'échantillonnage des mesures des ampèremètres et voltmètres de la batterie 10. Ainsi, Te est la période de temps qui sépare deux instants quelconques 5 successifs k et k-1 d'acquisition de la tension et de l'intensité du courant par le système 40. La période Te est typiquement une constante comprise entre 0,1 s et 10 s. Ici, la période Te est égale à 1 s à plus ou moins 20 % près. Par exemple, Te est égale à une seconde. [0037] Dans le modèle 62, wk est un vecteur de bruit d'état. Ici, le bruit wk est un bruit 10 blanc gaussien centré. Ce bruit représente l'incertitude sur le modèle utilisé. La matrice de covariance, à l'instant k, du bruit wk est notée Qk Elle est définie par la relation suivante : Qk = E(Wk * WkT), où E(...) est la fonction espérance mathématique. Le modèle 62 s'écrit également sous la forme Xk+1 = FkXk Bkik Wk, où - Fk est la matrice de transition d'état à l'instant k, 15 - Bk est le vecteur de commande à l'instant k. [0038] Le modèle 62 permet notamment de prédire l'état de charge SOCk.i à l'instant k+1 à partir de l'état de charge précédant SOCk. [0039] Le modèle 64 permet de prédire la valeur yk de la tension à l'instant k à partir de l'état de charge SOCk, de la tension VD,k et de l'intensité mesurée ik. Dans ce 20 modèle, vk est un bruit de mesure blanc gaussien centré. La matrice de covariance du bruit vk à l'instant k est notée Rk par la suite. Dans le cas particulier décrit ici, cette matrice Rk est une matrice à une seule colonne et à une seule ligne. Elle est définie par la relation Rk = E(Vk * VkT). Ce bruit vk est indépendant du bruit wk et du vecteur d'état initial xo. 25 [0040] On notera que le modèle 64 est non linéaire, car la fonction OCV(SOCk) est généralement non linéaire. À cause de cela, l'estimateur 60 implémente la version étendue du filtre de Kalman, plus connu sous l'acronyme EKF (Extended Kalman Filter). Dans cette version étendue, on se ramène à un modèle linéaire d'observation de la forme yk = Hkxk + ROk2.ik + vk en linéarisant le modèle 64 au voisinage du 30 vecteur xk. Typiquement, le modèle 64 est développé en série de Taylor au voisinage du vecteur xk. Puis on néglige les contributions des dérivées à partir du second ordre. Ici, la matrice Hk est donc égale à la dérivée première de la fonction OCV au voisinage de l'état de charge SOCk. Cette linéarisation du modèle 64 est typiquement réalisée pour chaque nouvelle valeur de l'état de charge SOCk. 35 [0041] L'estimateur 60 a besoin de connaître la capacité Cri,k3 et la résistance interne R0k2 pour être capable d'estimer l'état de charge SOCk+i. La capacité et la résistance interne de la cellule 18 varient au fur et à mesure qu'elle vieillit. Pour tenir compte de ce vieillissement, la capacité et la résistance interne de la cellule 18 sont estimées, respectivement, à des instants k3 et k2. Ici, un estimateur 66 estime la résistance 40 interne R0k2 à partir de la valeur mesurée yk2, de l'intensité mesurée ik2 et de l'état de charge SOCk2. Un autre estimateur 68 estime la capacité Cn,k3 à partir de l'intensité ik3 et de l'état de charge SOCk3. 3 02 92 9 7 11 [0042] La résistance interne et la capacité de la cellule 18 varient plus lentement que son état de charge. Ainsi, pour limiter la puissance de calcul requise pour estimer l'état de charge de la cellule sans pour autant dégrader la précision de cette estimation, les estimateurs 66 et 68 sont exécutés moins fréquemment que l'estimateur 60. Par la suite, les instants d'exécution des estimateurs 66 et 68 sont notés, respectivement, k2 et k3 pour les distinguer des instants k. Ici, l'ensemble des instants k2 et l'ensemble des instants k3 sont des sous-ensembles de l'ensemble des instants k. Entre deux instants successifs k2 et k2-1 et entre deux instants successifs k3 et k3-1, il s'écoule donc plusieurs périodes Te et plusieurs instants k. [0043] Ces estimateurs 66 et 68 sont aussi implémentés chacun sous la forme d'un filtre de Kalman. L'estimateur 66 utilise un modèle d'état 70 (figure 6) et un modèle d'observation 72 (figure 7). Dans ces modèles, les bruits w - 2,k2 - et - - v 2,k2 sont des bruits blancs gaussiens centrés. Les covariances des bruits W - 2,k2 et V2,k2 sont notées, respectivement, 0 -2,k2 et R2,k2 par la suite. Le modèle 72 d'observation permet de prédire la valeur d'une grandeur physique uk2 directement mesurable. La grandeur physique uk2 est ici la somme des N dernières valeurs mesurées yk. Elle est définie par la relation suivante : Uk2 Y k [0044] N est un entier strictement supérieur à un qui est compté comme cela sera décrit plus loin. Dans la relation ci-dessus et dans le modèle 72, l'instant k est égale à 20 l'instant k2. [0045] Le modèle 72 tient compte non seulement de l'état de charge SOCk, de la tension VD,k et de l'intensité mesurée ik à l'instant k = k2 mais également des N précédentes estimations de l'estimateur 60 et des N précédentes intensités mesurées, comprises entre les instants k2 et k2-1. Le fait de prendre en compte les 25 mesures et estimations intermédiaires entre les instants k2 et k2-1, permet d'augmenter la précision de l'estimation de la résistance interne R0k2. [0046] L'estimateur 68 utilise un modèle d'état 74 (figure 8) et un modèle d'observation 76. Dans les modèles 74 et 76, les bruits W 3,k3 et V3,k3 sont des bruits blancs gaussiens centrés. Les covariances des bruits W 3,k3 et V3,k3 sont notées 30 respectivement Q3,k3 et R3,k3 par la suite. On remarquera que le modèle 76 est un modèle linéaire de sorte qu'un simple filtre de Kalman peut être utilisé pour l'estimateur 68 au lieu d'un filtre de Kalman étendu. [0047] Le modèle 76 d'observation permet d'estimer une grandeur physique zk3 directement mesurable. La grandeur physique zk3 est ici la somme des N dernières 35 intensités mesurées ik. Elle est définie par la relation suivante : 302 92 97 12 [0048] Dans la relation ci-dessus et dans le modèle 76, l'instant k est égal à l'instant k3. Cette grandeur physique zk3 tient compte non seulement de l'intensité mesurée ik-1 à l'instant k-1 précédent l'instant k3 mais également de N précédentes intensités mesurées entre les instants k3 et k3-1. Ici, N est un entier strictement supérieur à un 5 qui est compté comme cela sera décrit plus loin. Il n'est pas nécessairement égal au N introduit dans le modèle 72. Le fait de prendre en compte des mesures et estimations intermédiaires entre les instants k3 et k3-1, permet d'augmenter la précision de l'estimation de la capacité Cn,k3. [0049] Le fonctionnement des estimateurs 60, 66 et 68 va maintenant être décrit à 10 l'aide du procédé de la figure 10 et dans le cas particulier de l'estimation de l'état de charge de la cellule 18. [0050] Le procédé débute par une phase 100 de réglage des différentes matrices de covariance nécessaires pour exécuter les estimateurs 60, 66 et 68. Plus précisément, lors d'une opération 102, les matrices de covariance Qk et Rk de l'estimateur 60 sont 15 automatiquement réglées à l'aide des relations suivantes : Qk = [NoGo,k(No)]-1 et Rk = I, où - No est un entier prédéterminé strictement supérieur à 1, - I est la matrice identité, et - Go,k(No) est défini par la relation suivante : 20 Go .k Fie =0 [0051] No est généralement choisi lors de la conception du système 40 puis figé une fois pour toute. Généralement, No est inférieur à 100. Par exemple, No est compris entre Set 15. Ici, No est choisi égal à 10. [0052] Le fait d'utiliser les relations précédentes simplifie considérablement le 25 réglage des matrices Qo et Ro ainsi que le réglage des matrices Qk et Rk comme on le verra pas la suite. En effet, le seul paramètre à choisir est la valeur de l'entier No. [0053] Lors d'une opération 104, les covariances Q2,0 et R2,0 sont aussi réglées. Par exemple, Q2,0 est choisi égal à [(13*ROini)/(3*NCeoms)]2, où - 13 est une constante choisie supérieure ou égale à 0,3 ou 0,5 et, de préférence, 30 supérieur à 0,8 et généralement inférieure à trois, - Nced est le nombre prévu de cycles de charge et de décharge de la cellule 18 avant que celle-ci atteigne sa fin de vie, et - Ns est le nombre de fois où la résistance interne est estimée par cycle de charge et de décharge de la cellule18. 35 [0054] La constante 13 représente, exprimée en pourcentage divisé par 100, l'écart entre la valeur de la résistance interne initiale RO'' et sa valeur en fin de vie. Typiquement, 13 est fixé par l'utilisateur ou mesuré expérimentalement. Nced est un 302 92 9 7 13 nombre de cycles qui peut être mesuré expérimentalement ou obtenu à partir des données du constructeur de la cellule 18. Ns est fixé par le procédé d'estimation de l'état de charge mis en oeuvre par le calculateur 44. Dans ce mode de réalisation, comme on le verra plus loin, la résistance interne est estimée une seule fois par 5 cycle. Par conséquent, Ns est pris égale à 1. [0055] A titre d'illustration, la covariance R2,0 est choisie égale à (2cmUmax/300)2, où Cm est l'erreur maximale du voltmètre 34 exprimée en pourcentage. [0056] Par la suite, les covariances 0 -2,k2 et R2,k2 sont considérées comme étant constantes et égales, respectivement, à Q2,0 et R2,0. 10 [0057] Lors d'une opération 106, les covariances Q3,0 et R3,0 sont réglées. Par exemple, la covariance Q3,0 est prise égale à [y*Cri'n'/(3*Nceol*Ns)]2, où y représente, exprimée en pourcentage divisé par 100, l'écart entre la capacité Cn'n' et la capacité de la cellule 18 en fin de vie. y est une constante choisie par l'utilisateur entre 0,05 et 0,8 et, de préférence, entre 0,05 et 0,3. Ici, y = 0,2. 15 [0058] La covariance R3,0 est par exemple choisie égale à [2*E,m*Imax/300]2, où En est l'erreur maximale de l'ampèremètre 36 exprimé en pourcentage. [0059] Par la suite, les covariances Q3,k3 et R3,k3 sont considérées comme étant constantes et prises égales, respectivement, à Q3,0 et R3,0. [0060] Une fois les matrices de covariance réglées, l'estimation de l'état de charge 20 de la cellule 18 peut commencer. [0061] Lors d'une phase 110, à chaque instant k, le voltmètre 34 et l'ampèremètre 36 mesurent, respectivement, la valeur yK et l'intensité ik et ces mesures sont immédiatement acquises par le système 40 et enregistrées dans la mémoire 42. La phase 110 est réitérée à chaque instant k. 25 [0062] En parallèle, l'estimateur 60 exécute une phase 114 d'estimation de l'état de charge à l'instant k de la cellule 18. [0063] Pour cela, lors d'une étape 116, l'estimateur 60 calcule une prédiction SÔCkik-i et une prédiction VD,k/k_i de, respectivement, l'état de charge de la cellule 18 et la tension VD aux bornes du circuit 54 à l'instant k. Dans les notations utilisées ici, 30 l'indice k/k-1 indique que cette prédiction est faite en prenant en compte seulement les mesures réalisées entre les instants 0 et k-1. On parle alors de prédiction a priori. L'indice k/k indique que la prédiction à l'instant k prend en compte l'ensemble des mesures réalisées entre les instants 0 et k. On parle alors de prédiction a posteriori. Les prédictions SÔCi et VD,kik-1 sont calculées à l'aide du modèle 62, de l'intensité 35 mesurée ik-1 et de la capacité Cn,k3. On notera que dans le modèle 62, la matrice F k_i de transition d'état est constante quel que soit k et n'a donc pas besoin d'être réévaluée à chaque instant k. [0064] Lors d'une étape 117, l'estimateur 60 calcule également la prédiction Pkik_i d'une matrice de covariance d'erreur d'estimation sur le vecteur d'état xk. 40 Typiquement, cela est réalisé à l'aide de la relation suivante : Pk/k-1 = Fk-1Pk-1/k-1Fk-1T+Qk-1 302 92 9 7 14 [0065] Ces différentes matrices Fk-1, Pk-1/k-1 et Qk-1 ont déjà été précédemment définies. [0066] Ensuite, lors d'une étape 118, l'estimateur 60 construit la matrice Hk en linéarisant le modèle 64 autour des prédictions SOCkik_i et VD,k/k-1. 5 [0067] Lors d'une étape 120, les matrices de covariance Qk et Rk sont automatiquement mises à jour. Ici, pour cela, l'étape 120 est identique à l'opération 102 en prenant en compte cette fois-ci la matrice Hk construite lors de l'étape 118. [0068] Après cela, lors d'une étape 122, l'estimateur 60 corrige les prédictions SOCkik-i et VD,k/k-1 en fonction d'un écart entre la valeur mesurée yk et une valeur j/k 10 prédite à partir du modèle 64. Cet écart est connu sous le nom de « innovation ». Cette étape 122 comprend typiquement : - une opération 124 de calcul de la prédiction j/k, puis - une opération 126 de correction des prédictions SOCkik-i et VD,k/k-1 et de la matrice Pk/k-1 pour obtenir les prédictions corrigées SQCk/k, VD,k/k et Pk/k. 15 [0069] Lors de l'opération 124, la prédiction j/k est calculée à l'aide du modèle 64 dans lequel la valeur de l'état de charge est prise égale à SOCkik_i et la valeur de la tension VD,k est prise égale à VD,k/k-i. L'écart entre la valeur mesurée yk et sa prédiction j/k est notée Ek par la suite. [0070] De nombreuses méthodes existent pour corriger les estimations à priori 20 SÔCkik-i et VD,k/k-1 à partir de l'innovation Ek. Par exemple, lors de l'opération 126, ces estimations sont corrigées à l'aide du gain de Kalman Kk. Le gain Kk est donné par la relation suivante Kk = Pk/k-1H-rk(HkPk/k-1H-rk+Rk)i. Ensuite, les prédictions a priori sont corrigées à l'aide de la relation suivante : xkik = KkEk. [0071] La matrice Pk/k_i est corrigée à l'aide de la relation suivante : 25 Pk/k = Pk/k-1-KkHkPk/k-1- [0072] Les étapes 116 à 122 sont réitérées à chaque instant k où une nouvelle estimations de l'état de charge de la cellule 18 doit être réalisée. Lors de chaque nouvelle itération, le vecteur d'état xk_i est initialisé avec les valeurs obtenues lors de la précédente itération de la phase 114 pour la cellule 18. 30 [0073] En parallèle, lors d'une étape 130, le calculateur 44 compare chaque nouvelle mesure de l'intensité ik à un seuil SH, de courant prédéterminé. Tant que l'intensité mesurée ne dépasse pas ce seuil SH' l'exécution de l'estimateur 66 et inhibée. À l'inverse, dès que l'intensité ik mesurée dépasse ce seuil SH' alors l'estimateur 66 est immédiatement exécuté. Le seuil SH, est généralement supérieur à Imax/2 et, 35 avantageusement supérieur à 0,8*Imax ou 0,9*Imax. [0074] L'estimateur 66 exécute une phase 140 d'estimation de la résistance interne R0k2 à l'instant k2. Ici, l'instant k2 est égal à l'instant k où l'intensité ik franchit le seuil SH,. [0075] Pour cela, lors d'une étape 142, l'estimateur 66 calcule la prédiction a priori 40 RQi de la résistance interne à partir du modèle 70. [0076] Ensuite, lors d'une étape 144, l'estimateur 66 calcule la prédiction P - 2,k2/k2-1 de la matrice de covariance de l'erreur d'estimation sur la résistance interne. Par 302 92 9 7 15 exemple, cette prédiction est calculée à l'aide de la relation suivante : P - 2,k2/k2-1 = P2,k2- 1/k2-1 ± Q2,0. On notera qu'ici, le modèle 72 est une fonction linéaire de la variable d'état. Il n'est donc pas nécessaire de la linéariser au voisinage de la prédiction R k2/k2- 1 pour obtenir la matrice H2,k2. Ici, cette matrice H2,k2 est égale à -N. 5 [0077] Lors d'une étape 148, l'estimateur 66 corrige la prédiction RÔ k2/k2-1 en fonction de l'écart entre la grandeur physique mesurée uk2 et une prédiction Ûk2 de cette même grandeur physique. Ici, N est une constante prédéterminée choisie strictement à un et, de préférence, supérieure à 10 ou 30. La grandeur uk2 est acquise par l'estimateur 66 au fur et à mesure que les valeurs yk sont mesurées et acquises. 10 [0078] Plus précisément, lors d'une opération 150, le calculateur 44 acquière la grandeur mesurée uk2 et calcule la prédiction Ûk2. L'acquisition de la grandeur uk2 est réalisée en sommant les N dernières mesures de la valeur mesurée yk. La prédiction Ûk2 est calculée à l'aide du modèle 72. Dans ce modèle 72, la valeur RO k2 est prise égale à la valeur ROk2/k2-1 précédemment calculée. 15 [0079] Ensuite, lors d'une opération 152, l'estimateur 66 corrige la prédiction RÔ k2/k2-1 en fonction de l'innovation E k2. L'innovation E k2 est égale à l'écart entre la grandeur mesurée uk2 et la grandeur prédite Ûk2. Par exemple, lors de l'opération 152, la même méthode que celle mise en oeuvre lors de l'opération 126 est employée. Ainsi, cette opération 152 n'est pas décrite ici plus en détail. Ensuite, la nouvelle estimation ROk2/k2 est utilisée lors des exécutions suivantes de l'estimateur 60 à la place de la précédente estimation R0k2-1 /k2-1. [0080] Le fait de déclencher l'exécution de l'estimateur 66 uniquement lorsque l'intensité mesurée ik est élevée permet d'accroître la précision de l'estimation de la résistance interne tout en diminuant en même temps la puissance de calcul nécessaire pour mettre en oeuvre ce procédé. En effet, la précision de la mesure de ampèremètre est plus élevée lorsque l'intensité ik est plus élevée. [0081] Également en parallèle des phases 110 et 114, le procédé comporte une étape 160 lors de laquelle à chaque instant k, l'estimation SOCk est comparée à un seuil haut prédéterminé SHs'. Si l'estimation SOCk descend en dessous de ce seuil SHs', alors le procédé se poursuit immédiatement par des étapes 162 et 164. Dans le cas contraire, l'étape 160 est réitérée au prochain instant k. Typiquement, le seuil SH0 est compris entre 90 % et 100 %. [0082] Lors de l'étape 162, le calculateur 44 commence par initialiser un compteur à zéro puis l'incrémente de 1 à chaque nouvelle mesure de l'intensité ik depuis le début de cette étape. De plus, à chaque instant k, l'intensité mesurée ik et l'estimation SOCk générée au même instant sont enregistrées, associée à cet instant k, dans une base de données. [0083] En parallèle de l'étape 162, lors de l'étape 164, le calculateur 44 compare chaque nouvelle estimation SOCk à un seuil prédéterminé SLs'. Le seuil SLs' est par exemple compris entre 0 % et 10 %. Tant que l'estimation SOCk reste supérieure à ce seuil SLs', l'étape 162 est réitérée à l'instant k suivant. Dans le cas contraire, dès que l'estimation SOCk pour la cellule 18 descend en dessous de ce seuil SLs' alors, le 302 92 9 7 16 calculateur 44 déclenche immédiatement l'exécution de l'estimateur 68 et arrête d'incrémenter le compteur. Ainsi, tant que ce seuil SLs' n'est pas franchi, l'exécution de l'estimateur 68 est inhibée. [0084] L'estimateur 68 estime, lors d'une phase 166, capacité C n,k3 à l'instant k3. 5 L'instant k3 est donc égal à l'instant k où l'exécution de l'estimateur 68 est déclenchée. [0085] Comme pour la phase 140, étant donné que l'estimateur 68 n'est pas exécuté à chaque instant k, l'instant k3-1 ne correspond pas à l'instant k-1. Au contraire, les instants k3 et k3-1 sont séparés par un intervalle de temps supérieur ou égal à NT, 10 où N est le nombre compté lors de l'étape 162. [0086] Les paramètres du filtre de Kalman de l'estimateur 68 sont initialisés avec les valeurs précédentes de ces paramètres obtenus à l'issue de la précédente itération à l'instant k3-1 de la phase 166. [0087] La phase 166 comprend : 15 - le calcul, lors d'une étape 170, de la prédiction Cn,k3/k3-1 à l'aide du modèle 74, - le calcul, lors d'une étape 172, de la prédiction P - 3,k3/k3-1 de la matrice de covariance de l'erreur d'estimation de la capacité, et - la correction, lors d'une étape 174, des prédictions Cn, k3/k3-1 et P - 3,k3/k3-1 - [0088] Lors des étapes 172 et 174, la matrice d'observabilité H 3,k3 est égale à [(SOCk 20 - SOCk_N)]*3600/(NT,). N est ici le nombre d'instants k écoulés entre l'instant où l'état de charge estimé est descendu en dessous du seuil SHs' et l'instant où l'état de charge estimé est descendu en dessous du seuil SLs'. La valeur N est égale à la valeur comptée lors de l'étape 162. [0089] L'étape 174 comprend une opération 176 d'acquisition de la grandeur 25 physique mesurée zk3 et de calcul de la prédiction 2k3 de la grandeur zk3. L'acquisition de la grandeur zk3 consiste ici à calculer la somme des N dernières intensités mesurées entre les instants k-1 et k-N. La prédiction 2k3 est obtenue à partir du modèle 76. [0090] Ensuite, lors d'une opération 178, l'estimateur 68 corrige la prédiction C n,k3/k3-1 30 en fonction de l'écart entre la grandeur mesurée zk3 et la grandeur prédite 2k3 pour obtenir l'estimation a posteriori de la capacité C n,k3/k3. Cette correction est par exemple réalisée comme décrit lors de l'étape 126. [0091] Ensuite, la capacité Cn,k3/k3 est transmise à l'estimateur 60 qui l'utilise pour estimer l'état de charge de la cellule 18 aux instants suivants. 35 [0092] Le fait de déclencher l'exécution de l'estimateur 68 uniquement après que la cellule 18 se soit en grande partie déchargée permet d'accroître la précision de l'estimation tout en diminuant en même temps la puissance de calcul nécessaire pour mettre en oeuvre ce procédé. [0093] À l'issue de la phase 166, lors d'une étape 180, le calculateur calcule l'état de 40 santé SOHk3 l'instant k3 à l'aide de la formule suivante : SOHk3= Cn,k3/Cn'n'. [0094] La figure 11 représente un procédé de détermination de l'état de charge de la batterie 10. À l'instant k, l'état de charge de la batterie 10 est déterminé à partir de l'état de charge de chacune des cellules de cette batterie. Par exemple, cela est réalisé de la façon suivante. Lors d'une étape 190, le calculateur 44 détermine l'état de charge de chaque étage de la batterie en additionnant l'état de charge de chaque cellule de cet étage. [0095] Ensuite, lors d'une étape 192, l'état de charge de la batterie est pris égal au plus petit des états de charge d'un étage déterminé lors de l'étape 190. [0096] Comme illustré par le procédé de la figure 11, la détermination de l'état de charge de charge de la batterie à chaque instant k nécessite de disposer d'une estimation de l'état de charge pour chacune des cellules à l'instant k. Une première solution consiste donc à exécuter en parallèle, pour chacune des cellules, le procédé d'estimation de la figure 10 en exécutant à chaque instant k la phase 114. Toutefois, pour limiter la puissance de calcul requise, sans dégrader la précision de l'état de charge déterminé pour la batterie, il est aussi possible de planifier l'exécution des estimations des états de charge des cellules comme décrit en référence au procédé de la figure 12. [0097] Le procédé de la figure 12 est décrit dans le cas simplifié où l'on utilise seulement trois niveaux de priorité, appelés respectivement niveau élevé, niveau moyen et niveau faible de priorité. De plus, on prend comme hypothèse que l'état de charge d'une cellule dont le niveau de priorité est élevé doit être estimé à chaque instant k et donc à une fréquence fe. L'état de charge d'une cellule dont le niveau de priorité est moyen doit seulement être estimé avec une fréquence trois fois plus faible et donc à une fréquence fe/3. Enfin, l'état de charge des cellules de niveau de priorité faible doit être estimé à une fréquence dix fois plus faible est donc à une fréquence fe/10. Dans cet exemple, pour les niveaux de priorité élevé et moyen, il existe un nombre de places limités connus à l'avance. En d'autres termes, le nombre de cellules affectées au niveau de priorité élevé est limité à un nombre maximal prédéterminé à l'avance. Il en est de même pour le nombre de cellules affectés au niveau de priorité moyen. [0098] Pour planifier les instants auxquels doivent être rafraîchies les estimations de 30 l'état de charge de chacune des cellules, le calculateur commence par affecter, lors d'une étape 198, un niveau de priorité à chaque cellule. [0099] L'étape 198 débute par une opération 200 lors de laquelle le système 40 acquiert la valeur mesurée yk de la tension entre les bornes de chacune des cellules. [00100] Ensuite, lors d'une opération 202, si la valeur mesurée yk est au-dessus d'un 35 seuil haut SHy ou, au contraire, en dessous d'un seuil bas SLy, alors le calculateur 44 affecte à cette cellule le niveau de priorité élevé tant qu'il reste assez de place dans ce niveau. Le seuil SHy est supérieur ou égal à 0,9*Umax et, de préférence, supérieur à 0,95*Umax. Le seuil SLy est quant à lui supérieur ou égal à Umin, et inférieur à 1,1*Unr' ou 1,05*Unn. Il est important de rafraîchir fréquemment l'estimation de l'état de charge 40 des cellules dont la tension est proche de Umax ou au contraire proche de Unn. En effet, une erreur sur l'estimation de l'état de charge d'une cellule dans une telle situation peut conduire à une dégradation des propriétés électriques et mécaniques de cette cellule. [00101] Ensuite, pour les autres cellules, lors d'une opération 204, le calculateur 44 calcule l'écart de tensions entre la valeur mesurée actuelle yk et une précédente 5 valeur yk-x, où X est un entier prédéterminé supérieur ou égal à un et généralement inférieur à 5 ou 10. Ici, X = 1. [00102] Lors d'une opération 205, le calculateur 44 identifie des cellules jumelles. Des cellules sont considérées comme étant « jumelles » si, au même instant k, elles ont le même écart de tensions et la même valeur mesurée yk. Pour cela, lors de l'étape 205, 10 le calculateur 44 compare l'écart de tensions et la valeur mesurée yk pour une cellule aux écarts de tensions et aux valeurs mesurées yk pour les autres cellules au même instant afin d'identifier, parmi ces autres cellules, la ou les cellules jumelles de cette cellule. L'identifiant de cette cellule et les identifiants de la ou des cellules identifiées comme étant ses jumelles sont alors regroupés dans un ensemble qui est enregistré 15 dans la mémoire 42. La comparaison ci-dessus est par exemple réalisée pour chacune des cellules de la batterie 10 dont l'identifiant n'a pas déjà été incorporé dans l'un des ensembles enregistrés de cellules jumelles. Par la suite, un niveau de priorité est affecté à une seule des cellules de chaque ensemble de cellules jumelles. Ainsi, l'opération 206 et les étapes 208 et 210 suivantes sont réalisées uniquement 20 pour les cellules qui n'ont pas de cellule jumelle et pour une seule cellule de chaque ensemble de cellules jumelles. [00103] Lors d'une opération 206, le calculateur classe les cellules par ordre décroissant de valeur absolue de l'écart calculé lors de l'opération 204. Puis, il affecte aux premières cellules de ce classement les places restantes associées à des 25 niveaux de priorité élevés. Ensuite, il affecte les places restantes associées à un niveau de priorité moyen aux cellules suivantes dans ce classement. Enfin, il affecte aux dernières cellules de ce classement le niveau de priorité faible. [00104] Une fois qu'un niveau de priorité a été affecté à chaque cellule, lors d'une étape 208, le calculateur 44 planifie les instants de rafraîchissement des estimations 30 des états de charge des cellules en fonction de leur niveau de priorité. L'étape 208 est réalisée de manière à respecter la fréquence de rafraîchissement des estimations associée à chacun des niveaux de priorité. Pour cela, par exemple, le calculateur 44 réserve d'abord les instants auxquels les estimations des cellules de niveau de priorité élevé doivent être rafraîchies. Ensuite, il réserve les instants auxquels les 35 estimations de l'état de charge des cellules de niveau de priorité moyen doivent être rafraîchies en tenant compte des instants de rafraîchissement déjà réservés. Enfin, il procède de même avec les cellules affectées d'un niveau de priorité faible. [00105] Pour illustrer cela, on suppose qu'un niveau de priorité élevé a été affecté à la cellule 18, qu'un niveau de priorité moyen a été affecté aux cellules 19 et 20 et qu'un 40 niveau de priorité faible a été affecté à la cellule 21. De plus, on impose que pendant une période T,, le calculateur exécute au plus deux fois la phase 114 du procédé de la figure 10. Le résultat obtenu avec ces hypothèses est représenté sur la figure 13. Sur cette figure, les instants k à k+11 ont été représentés en abscisse. Au-dessus de chacun de ces instants k, deux cases symbolisent le fait que le calculateur 44 peut, à chaque instant k, exécuter deux fois la phase 114 du procédé de la figure 10. Dans chacune de ces cases, on a indiqué le numéro de la cellule pour lequel la phase 114 est exécutée. Lorsque aucun numéro ne figure dans cette case, cela signifie que le procédé de la figure 10 n'est pas exécuté et donc que la puissance de calcul épargnée peut être utilisée à d'autres fins comme, par exemple, l'exécution des estimateurs 66 et 68. [00106] Enfin, lors d'une étape 210, pour chaque cellule affectée d'un niveau de priorité, le calculateur 44 exécute la phase 114 à l'instant planifié pour cette cellule. En dehors de ces instants planifiés, le calculateur inhibe l'exécution complète de la phase 114 pour cette cellule. De même, l'exécution de la phase 114 pour les cellules jumelles auxquelles aucun niveau de priorité n'a été affecté est également inhibé. [00107] En parallèle, lors d'une étape 212, pour chaque cellule jumelle à laquelle aucun niveau de priorité n'a été affecté, l'estimation de l'état de charge de cette cellule est prise égale à la dernière estimation calculée lors de l'étape 210 pour une cellule jumelle de cette cellule. Ainsi, la phase 114 est exécutée pour une seule des cellules jumelles. Cela permet de réduire la puissance de calcul requise pour déterminer l'état de charge de batterie sans dégrader la précision de cette détermination. [00108] Éventuellement, en parallèle de l'étape 210, à chaque instant k, le calculateur 44 exécute aussi une étape 214 de prédiction de l'état de charge pour chacune des cellules qui ne sont pas traitées lors de l'étape 210 à cet instant k. L'étape 214 consiste à exécuter seulement l'étape 116 de prédiction, sans exécuter l'étape 122 de correction, pour toutes les cellules pour lesquelles, au même instant, la phase 114 complète d'estimation n'est pas exécutée. En effet, l'étape 116 de prédiction est beaucoup moins gourmande en puissance de calcul que l'étape 122 et peut donc être exécutée, par exemple, à chaque instant k. Ainsi, lorsque l'étape 214 est mise en oeuvre, on dispose à chaque instant k d'une nouvelle estimation de l'état de charge pour chacune des cellules de la batterie. [00109] Les étapes 198 et 208 sont réitérées à intervalle régulier pour mettre à jour le niveau de priorité affecté à chacune de ces cellules et donc les fréquences de rafraîchissement de l'estimation de l'état de charge de ces cellules. Ce procédé de planification des instants de rafraîchissement des estimations des états de charge des cellules permet de limiter la puissance de calcul requise sans pour autant dégrader la précision de l'état de charge déterminé pour la batterie. En effet, le procédé de la figure 12 exploite le fait que les cellules dont les écarts de tensions sont faibles sont des cellules qui se déchargent ou se chargent peu et donc dont l'état de charge ne varie pas rapidement. On peut donc estimer l'état de charge de ces cellules à une fréquence moins élevée sans que cela dégrade la précision de l'état de charge déterminé pour la batterie. [00110] Au cours de l'exécution des procédés des figures 10 et 11, à chaque fois que l'état de charge SOCk d'une cellule à un instant donné doit être utilisé pour un calcul, l'état de charge SOCk est pris égal au dernier état de charge estimé ou prédit pour cette cellule. Autrement dit, on considère que l'état de charge reste constant entre deux instants successifs où il est estimé ou prédit. [00111] On notera également qu'à chaque fois que le calculateur 44 exécute la phase 114 d'estimation pour une cellule, il récupère les informations nécessaires pour cette exécution à partir des valeurs obtenues à l'issue de la précédente exécution de cette phase pour la même cellule. C'est notamment le cas par exemple pour les variables d'état. On notera cependant que l'instant d'exécution précédent n'est alors pas nécessairement l'instant k-1 mais peut être l'instant k-3 ou k-10 selon le niveau de priorité affecté à cette cellule. [00112] De nombreux autres modes de réalisation du procédé d'estimation de l'état de charge d'une cellule sont possibles. Par exemple, la figure 14 représente un autre agencement d'estimateurs. Cet autre agencement est identique à celui de la figure 3 sauf que les estimateurs 66 et 68 sont remplacés par un seul estimateur 230. L'estimateur 230 estime simultanément la capacité et la résistance interne de la cellule 18. L'estimateur 230 est exécute moins fréquemment que l'estimateur 60. Ici, on note k4 les instants d'exécution de l'estimateur 230 et donc C n,k4 et R0k4 la capacité et la résistance interne estimée à l'instant k4. L'ensemble des instants k4 est un sous-ensemble des instants k. [00113] L'estimateur 230 estime en même temps la capacité Cn,k4 et la résistance interne ROka. Cet estimateur 230 implémente un filtre de Kalman qui utilise un modèle d'état 232 (figure 15) et un modèle d'observation 234 (figure 16). [00114] Le fonctionnement de cet estimateur 230 va maintenant être décrit en référence au procédé de la figure 17 et dans le cas particulier de la cellule 18. Ce procédé de la figure 17 est identique au précédé de la figure 10 sauf que les étapes 130 à 174 sont remplacées par des étapes 240, 242, 244 et une phase 246 d'estimation de la capacité et de la résistance interne. [00115] Lors de l'étape 240, le calculateur 44 compare à chaque instant k la valeur mesurée yk à un seuil haut SHy2. Typiquement, ce seuil SHy2 est supérieur ou également à 0,8*Umax ou 0,9*Umax. Les étapes 242 et 244 sont exécutées uniquement si valeur mesurée yk descend en dessous de ce seuil SHy2. [00116] Lors de l'étape 242, la calculateur 44 commence par initialiser un compteur à zéro puis incrémente ce compteur de 1 à chaque nouvel instant k. De plus, à chacun de ces instants k, l'intensité mesurée ik, la valeur yk, l'état de charge SOCk et la tension estimée VD,k sont enregistrés, associés à cet instant k, dans une base de données. [00117] En parallèle de l'étape 242, lors de l'étape 244, le calculateur 44 compare, à 40 chaque instant k, la nouvelle valeur mesurée yk à un seuil bas SLy2 de tension. Ce seuil SLy2 est inférieur ou égal à 1,2*Um,, ou à 1,1*Um,n et supérieur ou égal à Umm 302 92 97 21 [00118] Dès que la valeur mesurée yk descend en dessous du seuil SLy2, l'incrémentation du compteur lors de l'étape 242 est arrêtée et l'exécution de l'estimateur 230 est déclenchée. À l'inverse, tant que la valeur mesurée yk reste supérieure à ce seuil SLy2, l'exécution de l'estimateur 230 est inhibée. 5 [00119] L'estimateur 230 exécute la phase 246. Comme précédemment, on note que les instants ka et k4-1 sont séparés par un intervalle de temps supérieur ou également à NTe, où N est la valeur du compteur incrémenté lors de l'étape 242. Le fonctionnement de l'estimateur 230 se déduit du fonctionnement décrit précédemment pour les estimateurs 66 et 68. Il ne sera donc pas décrit ici plus en 10 détail. [00120] D'autres modèles électriques et donc d'autres modèles d'état peuvent être utilisés pour estimer l'état de charge de la cellule 18. Par exemple, dans une variante simplifiée, le circuit 54 est omis. À l'inverse, un modèle électrique plus complexe peut contenir plusieurs circuits RC parallèles raccordés électriquement en série les uns 15 avec les autres. Le modèle d'état de la cellule 18 doit alors être modifié en conséquence pour correspondre à ce nouveau modèle électrique de la cellule. Toutefois, tout ce qui a été décrit précédemment s'applique sans difficulté avec un tel modèle d'état modifié. Pour des exemples de modèles d'état modifiés, le lecteur peut se référer à la demande W02006057468. 20 [00121] Les paramètres RD et CD du modèle 50 peuvent aussi être estimé au lieu d'être considérés comme des paramètres constants prédéterminés. A cet effet, ces deux paramètres RD et CD sont par exemple introduits dans le vecteur d'état xk qui devient alors [SOCk, VD,k, RD,k et CD,k]T. Par exemple, le modèle d'état est modifié pour incorporer les deux équations suivantes RD,k+1 = RD,k et CD,k+1 = CD, k. 25 [00122] Le vecteur d'état xk peut aussi être complété par la température de la cellule de manière à estimer cette température en même temps que l'état de charge de cette cellule. [00123] La cellule peut aussi être équipé de capteurs supplémentaires tel qu'un capteur de température. Dans ce cas, le modèle d'observation est modifié pour tenir 30 compte de ces grandeurs physiques mesurées supplémentaires. Le lecteur peut se référer à la demande W02006057468 pour des exemples de modèles d'observation modifiés. [00124] D'autres modèles électriques possibles pour modéliser la cellule électrique sont également présentés dans la partie 2 de Plett 2004, au chapitre 3.3. 35 [00125] L'ajustement automatique en continu des matrices de covariance Rk et Qk peut être réalisé différemment. Par exemple, la méthode dite de « Covariance Matching » décrite dans l'article suivant peut être appliquée : Mehra, R.K: « On the identification of variances and adaptative Kalman Filtering », Automatic Control, IEEE Transaction on, Volume 15, n°2, pages 175-184, avril 1970. Cette méthode s'applique après un réglage initial des matrices Ro et Qo, par exemple, comme décrit lors de l'opération 102. [00126] Dans une autre variante, les matrices Qo, Ro, Qk et Rk ne sont pas réglées comme décrit en référence aux opérations 102 et 120. Par exemple, ces matrices sont réglées en mettant en oeuvre une méthode conventionnelle. Dans un cas simplifié, elles sont constantes. Par exemple, la matrice Ro est alors réglée à partir des données fournies par le constructeur des capteurs ou à partir de tests réalisés sur ces capteurs et la matrice Qo par essais successifs. [00127] L'étape 122 de correction de la prédiction peut être réalisée différemment. Par exemple, dans une méthode préférée, la correction de la prédiction de l'état de charge et de la tension VD,k, est réalisée en minimisant une fonction de coût J 10 quadratique qui comporte deux termes : - un terme lié à l'erreur de prédiction de la valeur mesurée, et - un autre terme lié à l'erreur d'estimation du vecteur d'état. [00128] Cette méthode est décrite en détail au chapitre 10.5.2 du livre suivant : Y. Bar-Shalom, et Al. : « Estimation With Applications to Tracking and Navigation, 15 Theory Algorithms and Software », Wiley Inter-science, 2001. [00129] Dans une autre variante, l'estimateur 60 n'est pas implémenté sous la forme d'un filtre de Kalman. Par exemple, l'état de charge est estimée en modélisant son évolution au cours du temps sous la forme d'un filtre RII (Réponse Impulsionnelle Infinie) dont les coefficients sont estimés par la méthode RLS (Recursive Least 20 Square). [00130] D'autres modèles d'état peuvent être utilisés pour estimer la résistance interne et la capacité de la cellule. Par exemple, le modèle 232 peut être remplacé par un modèle 250, représenté sur la figure 18. Dans le modèle 250, a, 13 et y sont des constantes dont les valeurs sont obtenues à partir de données du constructeur de 25 la cellule ou mesurées expérimentalement. Typiquement : - a est égal à 1 à plus ou moins 30% près ou 10% près, -13 est également égal à 1 à plus ou moins 30 % ou 10 %, et - y est typiquement compris entre 0,1 et 0,5. Par exemple, y est égal à 0,2 à plus ou moins 30 % près ou 10 %. 30 [00131] Dans le modèle 250, Nck est égal au nombre de cycles de charge/décharge de la cellule réalisés avant l'instant k. Ce nombre de cycles est par exemple mesuré en comptant le nombre de fois où l'état de charge de la cellule descend en dessous du seuil haut SHs' puis en dessous du seuil bas SLs'. wad,k est un bruit blanc gaussien centré. y est l'écart, exprimé en pourcentage divisé par 100, entre la 35 capacité initiale Cn'n' de la cellule et sa capacité en fin de vie. Ce modèle prend en compte le fait : - que la résistance interne augmente au fur et à mesure que la cellule vieillit, et - que la capacité de la cellule diminue au fur et à mesure que cellule vieillit. [00132] De façon similaire, le modèle d'état 70 peut être remplacé par le modèle 40 d'état suivant : R0k2.1 = (a+13Nck2/N )RÛ w CEOL,-- k2+ - 2, k2, où les différents symboles de ce modèle ont déjà été précédemment définis. 302 92 9 7 23 [00133] Le modèle d'état 74 peut être remplacé par le modèle d'état suivant : Cn,k3+1 = (1 -YNCk3/NCEOL)Cn, k3+V3,k3 où les différents symboles de ce modèle ont déjà été précédemment définis. [00134] Selon le modèle d'observation utilisé par l'estimateur 68, la grandeur zk3 peut 5 être calculée différemment. Par exemple, la grandeur zk3 est égale à la somme des N dernières intensités mesurées entre les instants k et k-N+1. Dans ce cas, lorsque N est égal à 1, zk3 = 1k3. [00135] Ce qui a été décrit précédemment pour l'initialisation des matrices de covariance Qk et Rk peut aussi être appliqué pour l'initialisation des matrices de 10 covariance des estimateurs 68 et 230. [00136] En variante, l'estimateur 68 n'est pas implémenté sous la forme d'un filtre de Kalman. Par exemple, la capacité est estimée en modélisant son évolution au cours du temps sous la forme d'un filtre RII (Réponse Impulsionnelle Infinie) dont les coefficients sont estimés par la méthode RLS (Recursive Least Square). 15 [00137] Les procédés des figures 10 et 17 peuvent être simplifiés en prenant N égal à une constante prédéterminée. Dans ce cas, N n'est pas compté et les étapes 160, 162, 240 et 242 peuvent être omises. Par exemple, N est choisi égale à un ou, au contraire, strictement supérieur à 1 ou Sou 10. [00138] Dans une autre variante, à chaque instant k entre les instants k3 et k3-1, 20 seule l'étape 170 de calcul d'un prédiction Cn,k est exécutée mais l'étape 174 de correction de cette prédiction n'est pas exécutée. Ainsi, on obtient une nouvelle prédiction de la capacité de la cellule à chacun de ces instants k tout en limitant la puissance de calcul requise. De façon similaire, à chaque instant k entre les instants k4 et k4-1, seule l'étape de calcul des prédictions de la capacité et de la résistance 25 interne est exécutée sans exécuter l'étape de correction de ces prédictions. Ainsi, dans ces variantes, la capacité de la cellule est prédite à chaque instant k mais cette prédiction est corrigée uniquement aux instants k3 ou k4. L'algorithme d'estimation de cette capacité n'est donc que partiellement exécuté entre les instants k3 et k3-1 ou k4 et k4-1 et complètement exécuté uniquement à l'instant k3 ou k4. 30 [00139] A chaque instant k entre les instants k3 et k3-1 ou entre les instants k4 et k41, la capacité peut être estimée en exécutant un premier algorithme, puis à l'instant k3 ou k4, la capacité est estimée en exécutant un second algorithme différent du premier algorithme et nécessitant une puissance de calcul plus importante. Les premier et second algorithmes ne correspondent pas nécessairement, comme décrit 35 précédemment, respectivement, à la seule étape 170 et à la phase 166 ou 246 d'un filtre de Kalman. Il peut aussi s'agir de deux algorithmes d'estimation complètement différents. [00140] L'étape 166 ou 246 d'estimation de la capacité de la cellule peut être déclenchée en réponse au franchissement d'un seuil sur l'état de charge, comme 40 décrit en référence à la figure 10, ou en réponse au franchissement d'un seuil sur la tension comme décrit en référence à la figure 17. Ces étapes 166 et 246 peuvent aussi être déclenchées en réponse au franchissement d'un seuil de quantité de 302 92 9 7 24 courant débité. Pour cela, à partir du moment où la tension ou l'état de charge de la cellule est descendu en dessous d'un seuil haut prédéterminé, à chaque instant k, le calculateur 44 calcule la quantité QCk de courant débité à l'aide de la relation suivante : QCk = QCk-1 + kTe. Dès que QCk dépasse un seuil haut SHQ, alors la phase 5 166 ou 246 est exécutée. A l'inverse, tant que la quantité QCk reste supérieure au seuil SHQ, l'exécution des phases 166 ou 246 est inhibée. En variante, la quantité QCk peut aussi être calculée sur une fenêtre glissante contenant les N derniers instants k, où N est une constante prédéterminée. [00141] Dans un autre mode de réalisation, le déclenchement des estimations de la capacité et/ou de la résistance interne en réponse au franchissement d'un seuil est omis. Par exemple, ces estimations sont déclenchées à intervalle régulier. Cette intervalle régulier est égal à Te si la puissance de calcul disponible est suffisante pour estimer cette capacité et cette résistance interne à chaque instant k. [00142] De nombreux autres modes de réalisation du procédé de la figure 12 sont 15 possibles. Par exemple, l'opération 205 peut être omise. Dans ce cas, aucune cellule jumelle n'est identifiée et l'étape 212 est aussi omise. [00143] L'opération 202 peut être réalisée différemment. Par exemple, un seul des seuils haut et bas est utilisé. L'opération 202 peut aussi être omise. [00144] Le nombre de niveaux de priorité peut être quelconque et supérieur à au 20 moins deux ou trois. D'autres méthodes pour affecter un niveau de priorité aux cellules sont possibles. Par exemple, le niveau de priorité d'une cellule peut être calculé à l'aide d'une formule qui relie son niveau de priorité à son écart de tension et de sa tension. Dans ce dernier cas, les opérations de comparaison sont omises. [00145] La méthode décrite pour associer des instants de rafraîchissement aux 25 cellules en fonction de leurs niveaux de priorité n'est qu'un exemple. Tout autre procédé connu d'ordonnancement de tâches en fonction du niveau de priorité de ces tâches peut être adapté au cas décrit ici de l'ordonnancement des instants de rafraîchissement des estimations des états de charge des cellules. [00146] La planification des instants de rafraîchissement de l'estimation de l'état de 30 charge de chacune des cellules décrite en référence à la figure 12 peut être omise. Par exemple, cela sera le cas si la puissance de calcul nécessaire pour l'estimation de l'état de charge de chacune des cellules à chaque instant k est disponible. [00147] En variante, le calculateur 44 comprend plusieurs sous-calculateurs programmables capables d'exécuter chacun et en parallèle, le procédé d'estimation 35 de la figure 10 ou 17 pour des cellules respectives. [00148] L'état de santé d'une cellule peut aussi être calculé à l'aide de la relation suivante : SOHK = ROK/RO". [00149] La batterie 10 peut être remplacée par tout type de batterie, comme par exemple, une batterie au plomb, une super-capacité ou une pile à combustible. Dans 40 ce cas, le modèle d'état et/ou le modèle d'observation de l'estimateur 60 peuvent éventuellement être adaptés pour tenir compte de la technologie de la batterie. [00150] Ce qui a été décrit précédemment s'applique aussi au véhicule hybride, c'est-à-dire au véhicule dont l'entraînement des roues motrices est à la fois ou, en alternance, assuré par un moteur électrique et un moteur thermique à combustion interne. Le véhicule 2 peut aussi être un camion, une moto ou un tricycle et de façon générale, tout véhicule capable de se déplacer en entraînant des roues motrices à l'aide d'un moteur électrique alimenté par une batterie. Par exemple, il peut s'agir d'un monte-charge. [00151] La batterie 10 peut être rechargée par l'intermédiaire d'une prise électrique qui permet de la raccorder électriquement à un réseau de distribution d'électricité. La 10 batterie 10 peut aussi être rechargée par un moteur thermique à combustion interne. REVENDICATIONS1. Procédé automatique d'estimation de l'état de charge d'une cellule d'une batterie, ce procédé comportant : a) l'acquisition (110), à un instant k, d'une valeur mesurée yk de la tension entre des bornes de la cellule, et d'une intensité mesurée ik du courant de charge ou de décharge de la cellule, b) l'estimation (114) de l'état de charge SOCk de la cellule, à l'instant k, à partir de la valeur mesurée yk, de l'intensité mesurée ik et d'une capacité Cn,k3 de la cellule, cette capacité Cn,k3, exprimée en Ah, représentant la quantité maximale d'énergie électrique stockable dans la cellule à l'instant k3, l'instant k3 étant l'instant auquel la capacité Cn,k3 a été estimée le plus proche de l'instant k, et c) l'estimation (166; 246) de la capacité Cn,k3 à un instant k3, les instants k3 étant moins fréquent que les instants k, cette estimation comportant : - le calcul (170) d'une prédiction de la capacité Cn,k3 à l'aide d'un modèle d'état qui relie la capacité Cn,k3 à la capacité Cn,k3-1 de cette même cellule à un instant précédent k3-1, puis - la correction (174) de la prédiction de la capacité Cn,k3 en fonction d'une intensité mesurée, à l'instant le plus proche de l'instant k3, du courant de charge ou de décharge de la cellule, caractérisée en ce que la correction de la prédiction de la capacité Cn,k3 comporte : - l'acquisition (176), à l'instant k3, d'une grandeur physique mesurable zk3 définie par la relation suivante : ni -= k- I =k- N où k est l'instant le plus proche de l'instant k3 et N est un entier supérieur ou égal à 25 un, zk3 étant égale à ik-i lorsque N est égal à un, - le calcul (176) d'une prédiction 2k3 de la grandeur physique mesurable zk3 à l'aide du modèle d'observation suivant : SUC SOC .3600.0 k3 '7 où : - k est l'instant le plus proche de l'instant k3, 30 - SOCk_N est l'état de charge de la cellule à l'instant antérieur k-N, - Cn,k3 est, dans ce modèle d'observation, la prédiction de la capacité de la cellule à l'instant k3 calculée juste avant d'être corrigée, et- la correction (178) de la prédiction de la capacité Cn,k3 en fonction de l'écart entre la grandeur physique acquise zk3 et la prédiction 2k3 calculée à partir du modèle d'observation. 2. Procédé selon la revendication 1 dans lequel : - la correction de la prédiction de la capacité Cn,k3 est inhibée (164; 244) tant que la valeur d'un paramètre n'a pas franchi un premier seuil prédéterminé, ce paramètre étant choisi dans le groupe composé de la valeur mesurée yk, de l'état de charge estimé SOCk et d'une quantité Qk de courant débité par la cellule entre l'instant k et un instant précédent, et - la correction de la prédiction de la capacité Cn,k3 est déclenchée (164; 244) en réponse au fait que ce paramètre descend en dessous du premier seuil prédéterminé si le paramètre est la valeur mesurée yk ou l'estimation de l'état de charge SOCk, ou monte au-dessus du premier seuil prédéterminé si ce paramètre est la quantité Qk de courant débitée. 3. Procédé selon la revendication 2, dans lequel le procédé comprend le comptage (162 ; 242) du nombre d'itération de l'étape a) entre un instant t1 où la valeur mesurée yk ou l'état de charge estimé SOCk descend en dessous d'un second seuil prédéterminé et un instant t2 suivant où la valeur mesurée yk ou l'état de charge estimé SOCk descend en dessous du premier seuil prédéterminé, puis lors de la correction (178) de la prédiction de la capacité Cn,k3 l'entier N est choisi égal au nombre d'itérations comptées entre les instants t1 et t2 et strictement supérieur à un. 4. Procédé selon la revendication 3, dans lequel le premier seuil prédéterminé correspond à un état de charge de la cellule inférieur à 20 % et le second seuil prédéterminé correspond à un état de charge de la cellule supérieur à 90 %. 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel le 30 procédé comprend le réglage (106) de la covariance du bruit sur la mesure de l'intensité ik puis la correction (178) de la prédiction de la capacité Cn,k3 en fonction de cette covariance du bruit sur la mesure de l'intensité ik. 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel le 35 modèle d'état utilisé lors du calcul de la prédiction de la capacité Cn,k3 est le suivant : Cn,k3 = Cn,k3-1, et le procédé comporte : - le calcul (172) d'une prédiction de la covariance P3,k3 de l'erreur d'estimation de la capacité Cn,k3 en fonction d'une covariance Q3,0 du bruit de l'erreur de modélisation introduite par le modèle d'état utilisé, et 302 92 9 7 28 - le réglage (106) automatique de la covariance Q3,0 à l'aide de la relation suivante : Qd ,k3 = RyCn)/(3NceoiNs)]2, où : - y est une constante prédéterminée dont la valeur est strictement inférieure à un et supérieure à zéro, 5 - C'' est la capacité maximale nominale de la cellule, - Nced est le nombre de cycles de décharge de la cellule avant que celle-ci atteigne sa fin de vie, ce nombre étant prédéterminé, - Ns est le nombre d'itérations de l'étape c) dans un cycle de décharge de la cellule. 10 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel les étapes a) et b) sont exécutées au moins N-fois entre les instants k3 et k3-1, où N est un nombre entier supérieur à trois. 15 8. Procédé selon l'une quelconque des revendications précédentes, dans lequel le procédé comporte également le calcul (180) de l'état de santé SOHk3 de la cellule à l'aide de la relation suivante : SOHk3 = Cn,k3/ C'', où : - Cn,k3 est la capacité de la cellule estimée à l'instant k3 le plus récent, et - C'' est la capacité maximale nominale de la cellule. 20 9. Support (42) d'enregistrement d'informations, caractérisé en ce qu'il comporte des instructions pour l'exécution d'un procédé d'estimation conforme à l'une quelconque des revendications précédentes, lorsque ces instructions sont exécutées par un calculateur électronique. 25 10. Système de gestion d'une batterie équipé d'au moins une cellule, ce système comportant un calculateur électronique (44) programmé pour : a) acquérir, à un instant k, une valeur mesurée yk de la tension entre des bornes de la cellule, et une intensité mesurée ik du courant de charge ou de décharge de la cellule, b) estimer l'état de charge SOCk de la cellule, à l'instant k, à partir de la valeur mesurée yk, de l'intensité mesurée ik et d'une capacité Cn,k3 de la cellule, cette capacité Cn,k3, exprimée en Ah, représentant la quantité maximale d'énergie électrique stockable dans la cellule à l'instant k3, l'instant k3 étant l'instant auquel la capacité Cn,k3 a été estimée le plus proche de l'instant k, et c) estimer la capacité Cn,k3 à un instant k3, les instants k3 étant moins fréquent que les instants k, cette estimation comportant : - le calcul d'une prédiction de la capacité Cn,k3 à l'aide d'un modèle d'état qui relie la capacité Cn,k3 à la capacité Cn,k3-1 de cette même cellule à un instant précédent k3-1, puis- la correction de la prédiction de la capacité Cn,k3 en fonction d'une intensité mesurée, à l'instant le plus proche de l'instant k3, du courant de charge ou de décharge de la cellule, caractérisée en ce que le calculateur électronique (44) est également programmé 5 pour réaliser la correction de la prédiction de la capacité Cn,k3 en exécutant les étapes suivantes : - l'acquisition, à l'instant k3, d'une grandeur physique mesurable zk3 définie par la relation suivante : on =k - N où k est l'instant le plus proche de l'instant k3 et N est un entier supérieur ou égal à 10 un, zk3 étant égale à ik-i lorsque N est égal à un, - le calcul d'une prédiction 2k3 de la grandeur physique mesurable zk3 à l'aide du modèle d'observation suivant : 1 S C - 50 C 6 0 0 . où : - k est l'instant le plus proche de l'instant k3, 15 - SOCk_N est l'état de charge de la cellule à l'instant antérieur k-N, - Cn,k3 est, dans ce modèle d'observation, la prédiction de la capacité de la cellule à l'instant k3 calculée juste avant d'être corrigée, et - la correction de la prédiction de la capacité Cn,k3 en fonction de l'écart entre la grandeur physique acquise zk3 et la prédiction 2k3 calculée à partir du modèle 20 d'observation. 11. Véhicule automobile comportant : - au moins une roue motrice (6), - un moteur électrique (4) apte à entraîner en rotation cette roue motrice pour 25 déplacer le véhicule automobile, - une batterie (10) comportant au moins une cellule (18-21) apte à stocker de l'énergie électrique et, en alternance, à restituer de l'énergie électrique pour alimenter le moteur électrique, cette cellule comportant deux bornes (30, 32) par l'intermédiaire desquelles elle est électriquement raccordée au moteur électrique, 30 - un voltmètre (34) raccordé électriquement entre les bornes de la cellule pour mesurer la tension entre ces bornes,- un ampèremètre (36) raccordé en série avec la cellule électrique pour mesurer l'intensité du courant de charge ou de décharge de cette cellule, et - un système (40) de gestion de la batterie raccordé au voltmètre et à l'ampèremètre, ce système de gestion comportant un calculateur électronique (44) programmable 5 apte à estimer l'état de charge de la cellule de la batterie à partir des mesures du voltmètre et de l'ampèremètre, caractérisé en ce que le système (40) de gestion de la batterie est conforme à la revendication 10.