Brevet FR2889314A1 — publié le 20070202

La présente invention concerne un procédé et un système d'estimation de l'état de charge de moyens de stockage d'énergie électrique de véhicule automobile. Les moyens de stockage d'énergie électrique sont un élément central d'un véhicule automobile, intervenant aussi bien dans le fonctionnement du moteur, thermique ou électrique de celui-ci, que dans l'alimentation électrique d'organes comme l'éclairage, l'électronique embarquée, etc... Avec l'intégration toujours plus importante dans le véhicule d'éléments consommateurs d'électricité et le développement des véhicules utilisant l'électricité comme source d'énergie pour la propulsion tels que les véhicules hybrides ou tout électrique par exemple, il est aujourd'hui nécessaire de connaître en temps réel l'état de charge des moyens de stockage d'énergie électrique afin de gérer au mieux l'ensemble des fonctionnalités du véhicule. On connaît dans l'état de la technique des procédés d'estimation de l'état de charge de moyens de stockage d'énergie électrique de véhicule automobile, comme une batterie d'accumulateurs électrochimiques par exemple, qui consistent à acquérir le courant de la batterie puis à déterminer l'état de charge de la batterie en évaluant une cartographie prédéterminée d'états de charge en fonction de la valeur du courant acquise. De manière typique, ces procédés nécessitent régulièrement des recalages afin de corriger les dérives de leur estimation de l'état de charge, induites par une modification des caractéristiques de la batterie au cours du temps, comme son rendement de charge par exemple. Or de telles phases de recalage sont difficiles à exécuter car elles nécessitent beaucoup de temps ou une batterie totalement au repos. Le but de la présente invention est de résoudre le problème susmentionné en proposant un procédé et un système d'estimation de l'état de charge de moyens de stockage d'énergie électrique qui permettent une estimation de l'état de charge de ces moyens de stockage rapide, précise et qui, notamment, ne requièrent pas de phase de recalage. A cet effet, l'invention a pour objet un procédé d'estimation de l'état de charge de moyens de stockage d'énergie électrique de véhicule automobile, caractérisé en ce qu'il comprend les étapes consistant à : - appliquer à une entrée d'excitation des moyens de stockage un signal d'excitation prédéterminé et acquérir à une sortie de réponse de ceux-ci un signal de réponse correspondant; - identifier, en fonction du signal d'excitation et du signal de réponse correspondant, chaque modèle d'un ensemble fini de modèles du transfert entre l'entrée d'excitation et la sortie de réponse des moyens de stockage, ces modèles étant obtenus par linéarisation paramétrique d'un même modèle d'ordre non entier de ce transfert; sélectionner le modèle identifié le plus proche du transfert réel des 10 moyens de stockage; et - déterminer l'état de charge des moyens de stockage en fonction du modèle identifié sélectionné. Selon des modes particuliers de l'invention, le procédé comprend au moins l'une des caractéristiques suivantes: - le signal d'excitation est un courant, le signal de réponse est une tension, et le transfert est la transmittance des moyens de stockage; - le modèle d'ordre non entier de la transmittance des moyens de stockage est un modèle H selon la relation: 1+ 1+ H(p)=K co r où p est la variable de Laplace, K, cob' wr et n sont des paramètres de ce modèle, avec n un nombre positif, et l'ensemble des modèles est obtenu en attribuant au paramètre n des valeurs distinctes comprises dans une plage prédéterminée; - la plage prédéterminée est égale à [0,2; 0,8] pour des moyens de 25 stockage électrochimiques; - la détermination de l'état de charge des moyens de stockage est réalisée en fonction d'un vecteur des paramètres estimés du modèle d'ordre non entier selon la relation: ek =0k-1 +YX(Pk ek-1) où âk est le vecteur des paramètres estimés du modèle d'ordre non entier k, 0k_1 est un vecteur des paramètres estimés du modèle d'ordre non entier précédemment déterminé k_1, y est un nombre prédéterminé et cpk est le vecteur des paramètres du modèle d'ordre non entier correspondant au modèle identifié sélectionné ; - l'étape de détermination de l'état de charge des moyens de stockage comprend une étape d'interrogation d'une cartographie prédéterminée de valeurs d'état de charge en fonction de valeurs des paramètres du modèle d'ordre non entier; - l'étape d'interrogation de la cartographie de valeurs d'état de charge consiste à interroger celle-ci avec le vecteur Âk des paramètres estimés du modèle d'ordre non entier; - le signal d'excitation est un signal du type à séquence binaire pseudo-aléatoire; L'invention a également pour objet un système d'estimation de l'état de charge de moyens de stockage d'énergie électrique de véhicule automobile, caractérisé en ce qu'il comprend: - des moyens d'application à une entrée d'excitation des moyens de stockage d'un signal d'excitation et des moyens d'acquisition à une sortie de réponse de ceux-ci un signal de la réponse correspondant; - des moyens d'identification, en fonction du signal d'excitation et du signal de réponse correspondant, de chaque modèle d'un ensemble fini de modèles du transfert entre l'entrée d'excitation et la sortie de réponse des moyens de stockage, ces modèles étant obtenus par linéarisation paramétrique d'un modèle d'ordre non entier de ce transfert; - des moyens de sélection du modèle identifié le plus proche du transfert réel des moyens de stockage; et - des moyens de détermination de l'état de charge des moyens de stockage en fonction du modèle identifié sélectionné. Selon d'autres caractéristiques, ce système est propre à mettre en oeuvre un procédé du type susmentionné. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple, et faite en relation avec les dessins annexés dans lesquels: - la figure 1 est une vue schématique d'un système selon l'invention associé à une batterie d'accumulateurs de véhicule automobile; et - la figure 2 est un organigramme du fonctionnement du système de la figure 1 illustrant le procédé d'estimation selon l'invention mis en oeuvre par celui-ci. Sur la figure 1, on a illustré sous la référence générale 10 un système d'estimation de l'état de charge de moyens 12 de stockage d'énergie électrique de véhicule automobile, comme par exemple une batterie d'accumulateurs électrochimiques à électrolyte liquide ou étanche (plomb, lithium, ...). Cette batterie 12 est raccordée à un ensemble 14 d'organes du véhicule pour leur alimentation en énergie électrique, comme par exemple l'éclairage, l'électronique embarquée, le démarreur dans le cas d'un véhicule à moteur thermique, un moteur électrique dans le cas d'un véhicule électrique, etc..., ainsi qu'à d'autres organes pour la recharge de la batterie 12, comme par exemple l'alternateur du moteur d'un véhicule à moteur thermique, comme cela est connu en soi. Le système 10 d'estimation comprend, raccordés aux bornes de la batterie 12, des moyens 16 de sollicitation commandables de celle-ci, par exemple un générateur de courant alimenté par la batterie 12 elle-même. Ces moyens 16 sont adaptés pour délivrer à la batterie 12 un courant prédéterminé dont les caractéristiques sont choisies pour l'identification d'un nombre fini prédéterminé, par exemple égal à quatre, de modèles M1, M2, M3, M4 de la transmittance de la batterie, comme cela sera expliqué plus en détail par la suite. La transmittance est définie comme le rapport de la tension aux bornes de la batterie sur le courant qui lui est appliqué. De préférence, le courant de sollicitation des moyens 16 est un courant du type à séquence binaire pseudo-aléatoire dont le balayage de fréquences correspond à une plage de fréquences dans laquelle on souhaite identifier les modèles M1, M2, M3, M4 de la transmittance, comme cela est connu dans le domaine de l'identification. Le système 10 comprend également un capteur de tension 18 pour mesurer la tension u aux bornes de la batterie 12 et un capteur de courant 20 pour mesurer le courant i échangé par la batterie avec son environnement, c'est-à-dire les moyens 16 de sollicitation et l'ensemble 14 d'organes. Il comprend de plus une unité 22 de traitement du signal raccordée aux capteurs 18, 20 de tension et de courant et propre à estimer l'état de charge de la batterie 12 en fonction de la tension u et du courant i mesurés de la batterie 12. Enfin, le système 10 d'estimation comprend des moyens 24 de supervision raccordés aux différents éléments et organes venant d'être décrits pour leur commande, leur activation et leur désactivation, comme cela sera expliqué plus en détail par la suite. L'unité 22 de traitement du signal, par exemple mise en oeuvre sous la forme d'une unité à microprocesseur, est propre à estimer les paramètres d'un modèle H d'ordre non entier de la transmittance de la batterie 12, dont on sait qu'ils dépendent de l'état de charge de celle-ci, et à déterminer l'état de charge de la batterie 12 en fonction des paramètres estimés. Cette estimation des paramètres du modèle H d'ordre non entier se fonde sur l'identification d'un ensemble de modèles MI, M2, M3, M4 de la transmittance de la batterie 12, ces modèles M1, M2, M3, M4 étant obtenus à la suite de la linéarisation paramétrique du modèle H d'ordre non entier, comme cela est expliqué plus en détail par la suite. L'unité 22 de traitement du signal comprend un module 26 de prétraitement raccordé aux capteurs 18 et 28 et propre à traiter les signaux de mesure U, I délivrés par ceux-ci pour obtenir des signaux de forme requise pour une identification satisfaisante des modèles M1, M2, M3, M4 de la transmittance de la batterie 12. Le traitement mis en oeuvre par les moyens 26 consiste notamment à numériser, c'est-à-dire échantillonner en amplitude et en temps, les signaux de mesure u et i des capteurs 18 et 20 de façon appropriée et à réaliser un filtrage fréquentiel du type passebande pour éliminer les composantes continues du courant i et de la tension u mesurés et pour sélectionner la bande utile de fréquences dans laquelle les modèles M1, M2, M3, M4 de la transmittance sont identifiés, comme cela est connu en soi dans le domaine de l'identification. L'unité 22 comprend une mémoire tampon 28, raccordée au module 26, pour mémoriser le courant numérisé et filtré noté ik et la tension numérisée et filtrée notée uk délivrés par celui-ci. Cette mémoire tampon 28 est raccordée à un module 30 d'identification recevant de celle-ci une première partie du courant ik et une première partie correspondante de la tension uk à des fins d'identification des modèles M1, M2, M3, M4, les parties du courant ik et de la tension uk restantes, ou secondes parties, étant utilisées pour la validation des modèles MI, M2, M3, M4 identifiés, comme cela est connu en soi. Par exemple, si le courant ik et la tension uk filtrés et numérisés comprennent chacun N échantillons, {il, i2, ..., iN} et {u1, u2, ..., uN} respectivement, la première partie du courant ik consiste en les échantillons {il, i2, ..., i1} et la première partie de la tension uk correspond en les échantillons {ut u2, ..., ui}, où j est un nombre prédéterminé choisi préférentiellement égal à la partie entière de N/2. Les modèles M1, M2, M3, M4 sont obtenus par la linéarisation paramétrique du modèle H d'ordre non entier par exemple selon la relation: -n p wb) \n p \ où p est la variable de Laplace, K, wb, cor et n sont des paramètres de ce modèle, avec n positif non nécessairement entier. Le modèle H non entier selon la relation (1) modélise de manière efficace la transmittance de la batterie 12 pour une plage de fréquences étendue et ses paramètres dépendent de l'état de charge de la batterie 12. Toutefois, le modèle selon la relation (1) est non linéaire en ses paramètres. Aussi, l'identification de ceux-ci est typiquement réalisée en mettant en oeuvre un algorithme d'optimisation non linéaire. Par exemple, en utilisant une modélisation du type à erreur de sortie, l'identification des paramètres du modèle est réalisée au moyen d'un algorithme de régression 1+ U(10) = H(p) = K I(p) 1+ (1) non-linéaire, typiquement un algorithme de résolution du type Levenberg Marquardt. Ce type d'identification nécessite un nombre d'échantillons de courant et de tension important ainsi qu'une puissance de calcul importante. De plus, rien n'assure que les valeurs calculées par ce type d'algorithme soient optimales. En fixant le paramètre n à une valeur donnée, on obtient un modèle linéaire en ses paramètres (on parle alors de linéarisation paramétrique). Les paramètres restants, à savoir les paramètres K, 0b et cor sont alors identifiés de manière simplifiée. En effet, le modèle selon la relation (1) appartient à la classe de modèles dans le domaine temporel selon la relation: Dna u(t)+a1.Dnalu(t)+. ...+aq.Dnaqu(t)+...+aQ.DnaQu(t)= (2) b0.Dnb0 i(t) +b1.Dnbl i(t) +... +bm. Dnb,,, i(t) +... +bM.DnbM i(t) où Dnf(t) est la dérivée n1ème de la fonction f(t), aq et bq,0 On peut monter qu'une discrétisation de la relation (2), obtenue en approximant les dérivées continues à partir de la définition de Grunwald selon la relation: n Dnf(Kh) n E( 1)k hf((K k)h) (3) où h est un pas d'échantillonnage, K est le Kième instant in n(n -1)...(n k + 1) , , d'échantillonnage et = , s ecrit selon les relations: k! L M u(Kh)=- aqÉuq(Kh)+ bm.lm(Kh) (4) q=0 m=0 où : aq nQ aq= Q aq avec 0_ K uq(Kh)= 1( 1)k k=1 u((K k)h) ; et K (n s' im(Kh) = ( 1)k. bm k=O ski Comme on peut le noter la relation (4) définit un modèle linéaire en ses paramètres a'q et bm qu'il est possible d'identifier au moyen d'algorithme d'identification classique. Il s'agit donc d'une régression linéaire de sorte que l'identification des paramètres est par exemple réalisée au moyen d'une simple inversion de matrice, c'est-à-dire l'inversion d'un régresseur, ou bien les paramètres sont identifiés en utilisant un algorithme du type moindres carrés. Les valeurs des paramètres de la relation (2) sont alors calculées en fonction des valeurs des paramètres a'q et bm identifiées. L'identification d'un modèle linéaire M1, M2, M3, M4 en ses paramètres est donc plus rapide que l'identification globale du modèle H et nécessite en outre un nombre d'échantillons de mesure moins grand. De manière avantageuse, une plage prédéterminée de valeurs, dans laquelle le paramètre n est compris, est échantillonnée. Chaque échantillon sélectionné dans cette plage correspond alors à un modèle M1, M2, M3, M4 linéaire en ses paramètres que le module 30 d'identification identifie. Cette plage est déterminée lors d'une étude antérieure, par exemple sur un banc d'essai. i((K k)h) A titre d'exemple numérique, pour une batterie d'accumulateurs électrochimiques, n est compris dans la plage [0,2; 0,8] et l'ensemble des modèles de la transmittance à identifier correspond aux modèles M1, M2, M3 et M4 selon la relation (1) avec la valeur du paramètre n fixé respectivement à 0,2, 0,4, 0,6 et 0,8. Bien entendu, d'autres types de modèles que celui de la relation (1) sont possibles. Par exemple, il est possible d'utiliser un modèle H d'ordre non entier selon la relation: ( \-n1 ( -n2 1+p 1+p col) cor 1+p 1+p cor) coh H(p) = K (2) où K, cob, wr et coh n1, n2 sont les paramètres du modèle. Tout comme le modèle selon la relation (1), des modèles linéaires en leurs paramètres sont obtenus en fixant la valeur des paramètres n1 et n2 et des modèles linéaires en leurs paramètres sont sélectionnés en fixant les paramètres n1 et n2 à des valeurs comprises dans des plages de valeurs respectives prédéterminées. Le module 30 est donc propre à identifier chacun des modèles M1, M2, M3 et M4, c'est-à-dire déterminer la valeur des paramètres de celui-ci, en utilisant la première partie des signaux ik et Uk et en mettant en oeuvre un algorithme d'identification décrit précédemment par exemple. Le module 30 est raccordé à une mémoire 32 propre à mémoriser chacun des modèles identifiés. L'unité 22 de traitement du signal est par ailleurs propre à sélectionner le modèle identifié qui est le plus proche de ou semblable à la transmittance réelle de la batterie 12, la notion de proche ou semblable se comprenant vis-à-vis d'un critère de comparaison donné. L'unité 22 comprend à cet effet un module 34 de simulation et un module 36 de comparaison et de sélection. Le module 34 détermine, ou simule, la réponse en tension notée um de chacun des modèles M1, M2, M3 et M4 identifié par le module 30 à la seconde partie du courant ik d'une manière connue en soi. Le module 34 est raccordé au module 36 qui est propre à comparer la réponse en tension um de la batterie 12 de chacun des modèles M1, M2, M3, M4 à la seconde partie de la tension Uk, par exemple en calculant l'erreur quadratique entre ces deux signaux. Le module 34 est également propre à sélectionner le modèle qui donne la réponse en tension la plus proche de celle de la batterie 12, c'est-à- dire celui correspondant au minimum des erreurs quadratiques calculées. D'autres critères de comparaison sont possibles, comme par exemple la somme des valeurs absolues de la différence entre ces deux signaux, l'intercorrélation de ces deux signaux ou l'autocorrélation de l'erreur de sortie ou autres. L'unité 22 de traitement du signal comprend également un module 38 propre à calculer un vecteur âk de paramètres estimés du modèle H selon la relation (1), avec âk = (K wb n, où T est le symbole de la transposée, v représente l'estimée du paramètre v et k est le pas courant d'estimation de l'état de charge. Ce module 38 est raccordé au module 36 de comparaison ainsi qu'à une mémoire tampon 40 mémorisant la valeur âk_1 du vecteur des paramètres estimées calculé au pas d'estimation précédent k1 et est adapté pour calculer le vecteur âk courant selon la relation: 0k =0k-1+Yx((Pk -0k-1) (3) où y est un nombre prédéterminé, de préférence compris entre 0 et 1, et cpk=(K wb wr n)k est le vecteur des paramètres du modèle d'ordre non entier correspondant au modèle qui a été sélectionné par le module 36 de comparaison et de sélection au pas courant k. Enfin, l'unité 22 comprend un module 42 d'estimation de l'état de charge raccordé au module 38 et propre à estimer l'état de charge de la batterie 12 en fonction du vecteur âk des paramètres estimés par le module 38. Plus particulièrement, le module 42 comprend une cartographie prédéterminée de valeurs d'états de charge en fonction de valeurs de vecteur des paramètres 8k, par exemple déterminée expérimentalement lors de l'étude antérieure. Les moyens 42 évaluent cette cartographie en fonction du vecteur âk pour estimer l'état de charge courant de la batterie 12. Le fonctionnement du système susmentionné est maintenant décrit en relation avec l'organigramme de la figure 2 qui illustre le procédé d'estimation de l'état de charge selon l'invention mis en oeuvre par ce système. Suite à l'activation du système selon l'invention, dans une première étape 50, les moyens 24 de supervision scrutent l'occurrence de la satisfaction d'une condition de déclenchement de l'estimation de l'état de charge de la batterie 12. Par exemple, cette condition est une condition de périodicité, l'estimation étant déclenchée toutes les X minutes, où X est un nombre prédéterminé, ou de régularité, par exemple tous les Y kilomètres parcourus par le véhicule. En variante, l'estimation de l'état de charge est déclenchée lors d'une phase de repos du véhicule. Lorsque cette condition de déclenchement est satisfaite, les moyens 24 de supervision initialisent, dans une étape 52, l'estimation de l'état de charge. Cette initialisation consiste par exemple à stopper la régulation de la tension aux bornes de la batterie classiquement exécutée, en pilotant de manière appropriée l'alternateur. Une fois l'initialisation de l'étape 52 terminée, les moyens 24 de supervision commandent, en 54, les moyens 16 de sollicitation pour que ceux-ci délivrent à la batterie 12 le courant prédéterminé de sollicitation. Un courant d'excitation, somme de ce courant de sollicitation et du courant provenant de l'ensemble 14 d'organes, est alors appliqué en 54 à la batterie 12. Toujours en 52, ce courant d'excitation et la réponse en tension de la batterie 12 à celui-ci sont mesurés par le capteur de courant 20 et le capteur 18 de tension respectivement et délivrés au module 26. A la suite de l'étape 54 d'application du courant d'excitation et de la mesure de celui-ci et de la réponse en tension correspondante de la batterie 12, l'étape suivante du fonctionnement du système d'estimation est une étape 56 d'identification des paramètres de chacun des modèles M1, M2, M3 et M4 en fonction des mesures délivrées par les capteurs 18 et 20. A cet effet, le module 26 de prétraitement filtre et échantillonne en 58 les mesures du courant d'excitation I et la réponse en tension U correspondante et le courant ik et la tension uk numérisés et filtrés sont mémorisés dans la mémoire tampon 28. Dans une étape 60 suivante, les deux parties des signaux sont formées dans la mémoire 28, la première pour l'identification et la seconde pour la simulation des modèles, puis, dans une étape 62 suivante, les paramètres de chacun des modèles M1, M2, M3 et M4 sont identifiés par le module 30 d'identification en fonction du premier ensemble de signaux pour l'identification, ces paramètres identifiés étant mémorisés dans la mémoire 32. L'étape 56 d'identification des paramètres de chacun des modèles M1, M2, M3 et M4 est alors suivie d'une étape 64 de sélection du modèle identifié le plus proche de la transmittance réelle de la batterie 12. A cet effet, le module 34 détermine dans une étape 66 la réponse um de chacun des modèles identifiés à la seconde partie du courant ik. Dans une étape 68 suivante, le module 36 compare le signal um de ce modèle avec la seconde partie de la tension uk, puis sélectionne en 70, le modèle identifié présentant l'erreur quadratique la plus faible. L'étape 64 de sélection du modèle est alors suivi par une étape 72 d'estimation de l'état de charge, dans laquelle, en 74, le module 38 calculent le vecteur 8k des paramètres estimés du modèle H selon la relation (1) conformément à la relation (3), la valeur précédente de ce vecteur étant alors mémorisé dans le tampon 40. Dans une étape 76 suivante, l'état de charge de la batterie 12 est calculé par le module 42 d'estimation qui évalue sa cartographie pour la valeur du vecteur venant d'être calculée en 74. L'étape 76 est alors suivie, si nécessaire, d'une étape 78 de basculement du véhicule dans le mode dans lequel il fonctionnait avant l'exécution de l'étape d'initialisation 52, puis l'étape 78 boucle sur l'étape 50 pour un nouveau cycle d'estimation de l'état de charge de la batterie. Comme il est possible de le constater, le système et le procédé selon l'invention que met en oeuvre ce système permettent d'estimer rapidement un modèle d'ordre non entier de la transmittance de la batterie et d'en déduire de manière précise l'état de charge de celle-ci. En outre, l'estimation de l'état de charge de la batterie peut être réalisé alors que le véhicule est en fonctionnement, c'est-à-dire que la batterie est sollicitée par l'ensemble des organes qui lui sont raccordés. Dans le mode de réalisation décrit ci-dessus, le courant appliqué à la batterie est mesuré. En variante, l'estimation de l'état de charge de la batterie est réalisée lorsque le véhicule est à l'arrêt, la sollicitation de la batterie par les organes étant de ce fait minimale. Le capteur de courant n'est alors pas nécessaire puisque le courant appliqué à la batterie est connu et sensiblement égal à celui délivré par les moyens de sollicitation. Egalement en variante, d'autres types de signaux de sollicitations sont utilisés, comme par exemple des créneaux. Il a été décrit un système et un procédé d'estimation de l'état de charge de moyens de stockage d'énergie électrique se fondant sur l'identification d'un modèle non entier de la transmittance de ceux-ci. En variante, un autre transfert est utilisé pour cette estimation, par exemple l'inverse de la transmittance. Le système d'estimation est alors structuré et fonctionne d'une manière analogue à celui décrit précédemment. Il comprend des moyens de sollicitation des moyens de stockage selon une tension prédéterminée appropriée. Les moyens de mesure de celle-ci et du courant en réponse à cette sollicitation en tension. Ce système comprend de plus une unité de traitement d'informations propre à calculer les paramètres d'un modèle d'ordre non entier du transfert entre la tension et le courant d'une manière analogue à celle décrite précédemment et à estimer l'état de charge en fonction des paramètres calculés. Le nombre de modèles à identifier, c'est-à-dire le nombre d'échantillons du ou des paramètres étant la cause de la non linéarité vis-à-vis des paramètres du modèle de la transmittance sélectionné, est choisi pour satisfaire un compromis entre le temps de calcul total et la précision souhaitée sur l'état de charge. REVENDICATIONS 1. Procédé d'estimation de l'état de charge de moyens (12) de stockage d'énergie électrique de véhicule automobile, caractérisé en ce qu'il comprend les étapes consistant à : - appliquer (en 54) à une entrée d'excitation des moyens de stockage un signal d'excitation prédéterminé et acquérir à une sortie de réponse de ceux-ci un signal de réponse correspondant; - identifier (en 56), en fonction du signal d'excitation et du signal de réponse correspondant, chaque modèle d'un ensemble fini de modèles du transfert entre l'entrée d'excitation et la sortie de réponse des moyens de stockage, ces modèles étant obtenus par linéarisation paramétrique d'un même modèle d'ordre non entier de ce transfert; -sélectionner (en 64) le modèle identifié le plus proche du transfert réel des moyens de stockage; et - déterminer (en 72) l'état de charge des moyens de stockage en fonction du modèle identifié sélectionné. 2. Procédé selon la revendication 1, caractérisé en ce que le signal d'excitation est un courant, le signal de réponse est une tension, et le transfert est la transmittance des moyens de stockage. 3. Procédé selon la revendication 2, caractérisé en ce que le modèle d'ordre non entier de la transmittance des moyens de stockage est un modèle H selon la relation: où p est la variable de Laplace, K, wb, cor et n sont des paramètres de ce modèle, avec n un nombre positif, et en ce que l'ensemble des modèles est obtenu en attribuant au paramètre n des valeurs distinctes comprises dans une plage prédéterminée. -n P 4. Procédé selon la revendication 3, caractérisé en ce que la plage prédéterminée est égale à [0,2; 0,8] pour des moyens de stockage électrochimiques. 5. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la détermination (en 72) de l'état de charge des moyens de stockage est réalisée en fonction d'un vecteur des paramètres estimés du modèle d'ordre non entier selon la relation: ek -0k-1 +Yx((Pk ek-1) où âk est le vecteur des paramètres estimés du modèle d'ordre non entier k, ek_1 est un vecteur des paramètres estimés du modèle d'ordre non entier précédemment déterminé 1(4, y est un nombre prédéterminé et cpk est le vecteur des paramètres du modèle d'ordre non entier correspondant au modèle identifié sélectionné. 6. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'étape (52) de détermination de l'état de charge des moyens de stockage comprend une étape (76) d'interrogation d'une cartographie prédéterminée de valeurs d'état de charge en fonction de valeurs des paramètres du modèle d'ordre non entier. 7. Procédé selon les revendications 5 et 6, caractérisé en ce que 20 l'étape d'interrogation de la cartographie de valeurs d'état de charge consiste à interroger celle-ci avec le vecteur âk des paramètres estimés du modèle d'ordre non entier. 8. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le signal d'excitation est un signal du type à séquence 25 binaire pseudo-aléatoire. 9. Système d'estimation de l'état de charge de moyens (12) de stockage d'énergie électrique de véhicule automobile, caractérisé en ce qu'il comprend: - des moyens (14, 16) d'application à une entrée d'excitation des 30 moyens de stockage d'un signal d'excitation et des moyens (20) d'acquisition à une sortie de réponse de ceux-ci un signal de la réponse correspondant; - des moyens (30) d'identification, en fonction du signal d'excitation et du signal de réponse correspondant, de chaque modèle d'un ensemble fini de modèles du transfert entre l'entrée d'excitation et la sortie de réponse des moyens de stockage, ces modèles étant obtenus par linéarisation paramétrique d'un modèle d'ordre non entier de ce transfert; - des moyens (34, 36) de sélection du modèle identifié le plus proche du transfert réel des moyens de stockage; et -des moyens (38, 42) de détermination de l'état de charge des moyens de stockage en fonction du modèle identifié sélectionné. 10. Système selon la revendication 9, caractérisé en ce qu'il est adapté pour mettre en oeuvre un procédé conforme à l'une quelconque des revendication 2 à 8.