FR 2877097 A1 20060428 FR 0411477 A 20041027 La présente invention concerne un système d'estimation de l'état de charge de moyens de stockage d'énergie électrique de véhicule et un procédé de calibrage d'un tel système. Les moyens de stockage d'énergie électrique, constitués par exemple d'une batterie, sont un élément central pour le fonctionnement d'un véhicule, par exemple du type automobile, motocycle, ou autre. Ils participent en effet à l'alimentation électrique de divers organes fonctionnels du véhicule, tels que par exemple le démarrage, l'électronique embarquée, l'éclairage,... Ils peuvent également être la ou une source d'énergie pour la traction d'un véhicule électrique ou hybride. Le rôle central des moyens de stockage d'énergie électrique dans le fonctionnement d'un véhicule nécessite de connaître leur état de charge, par exemple pour des besoins de sécurité liés au système de freinage ou de direction, de gestion du choix de la motorisation dans un véhicule hybride ou d'avertissement de leur recharge nécessaire adressé à un utilisateur d'un véhicule électrique. On connaît dans l'état de la technique des systèmes d'estimation de l'état de charge de moyens de stockage d'énergie électrique de véhicule automobile qui mesurent le courant échangé entre ceux-ci et leur environnement, et qui estiment leur état de charge en évaluant une cartographie prédéterminée d'états de charge en fonction du courant échangé mesuré. Toutefois ces systèmes présentent une robustesse faible vis-à-vis de variations des caractéristiques des moyens de stockage d'énergie électrique, comme par exemple la variation de leur rendement de charge. De tels systèmes nécessitent alors de manière typique des recalages fréquents et complexes afin de compenser leur erreur d'estimation. Le but de l'invention est de résoudre le problème susmentionné en proposant un système alternatif d'estimation de l'état de charge de moyens de stockage d'énergie électrique de véhicule. A cet effet, l'invention a pour objet un système d'estimation de l'état de charge de moyens de stockage d'énergie électrique de véhicule, comprenant: - des moyens de sollicitation des moyens de stockage d'énergie électrique selon un courant prédéterminé d'identification; - des moyens de mesure de la réponse en tension des moyens de stockage d'énergie électrique à ce courant d'identification; - des moyens d'indentification de paramètres d'un modèle prédéterminé de la transmittance des moyens de stockage d'énergie électrique en fonction du courant d'identification et de la réponse en tension des moyens de stockage d'énergie électrique à celui-ci; et - des moyens d'estimation de l'état de charge des moyens de stockage d'énergie électrique en évaluant, pour les paramètres identifiés, une fonction prédéterminée d'estimation de l'état de charge dont les variables sont les 10 paramètres du modèle. Selon d'autres caractéristiques, le système se caractérise en ce que: les moyens de sollicitation sont adaptés pour solliciter les moyens de stockage d'énergie électrique selon une séquence binaire pseudo-aléatoire de courant prédéterminée, - le modèle de la transmittance des moyens de stockage d'énergie électrique est un modèle non entier prédéterminé, - le modèle non entier est un modèle H non entier selon la relation: ( -nl ( \-n2 1+ p l+ p wb wr H(p) = K 1+P 1+p wr) wh i où p est la variable de Laplace, et n1, n2, K, wb, wr et wh sont les 20 paramètres du modèle, les moyens d'identification des paramètres du modèle de la transmittance sont adaptés pour identifier les paramètres du modèle prédéterminé de la transmittance en mettant en oeuvre un algorithme d'identification du type erreur de sortie, notamment un algorithme du type Marquardt. - il comprend en outre des moyens de mesure de la température des moyens de stockage d'énergie électrique et des moyens de sélection de la fonction prédéterminée d'estimation parmi un ensemble prédéterminé de fonctions d'estimation de l'état de charge, en fonction de la température mesurée. Un autre objet de la présente invention concerne un procédé de calibrage du système d'estimation de l'état de charge des moyens de stockage d'énergie électrique de véhicule susmentionné. Le procédé est caractérisé en ce qu'il comprend les étapes consistant à: sélectionner un nombre N prédéterminé de valeurs de l'état de charge SOCS, j=1,...,N, des moyens de stockage d'énergie électrique et identifier, pour chacune de ces valeurs SOCS, la valeur Qi du vecteur Q des paramètres du modèle prédéterminé de la transmittance des moyens de stockage d'énergie électrique; et - déterminer la fonction d'estimation en interpolant les valeurs sélectionnées de l'état de charge SOC1 en fonction des vecteurs Qi des paramètres identifiés du modèle. Selon une autre caractéristique, le procédé est caractérisé en ce que l'étape de détermination de la fonction d'estimation consiste à : déterminer M fonctions F;, i=1,...,M, des paramètres du modèle, M étant un nombre prédéterminé supérieur ou égal à 1; - déterminer, pour chacune des M fonctions F; , un indice de confiance 8i dans les résultats renvoyés par la fonction F; ; et - déterminer la fonction d'estimation de l'état de charge comme étant la combinaison des M fonctions F; des paramètres selon la relation: M soc = Ml E Si Fi (Q) Es1 =1 i=1 où SOC est l'état de charge des moyens de stockage d'énergie électrique. La présente 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 du système selon l'invention associé à des moyens de stockage d'énergie électrique de véhicule; - la figure 2 est un organigramme du procédé de calibrage selon l'invention; - les figures 3a à 3d sont des graphiques de fonctions des paramètres du modèle de la transmittance des moyens de stockage d'énergie électrique; et - la figure 4 est un graphique illustrant les variations de pondération en de fonction de la moyenne du courant d'identification appliqué aux moyens de 5 stockage. Sur la figure 1, on a illustré de manière schématique un système d'estimation de l'état de charge (SOC) de moyens 10 de stockage d'énergie électrique, par exemple une batterie du type plomb/acide, de véhicule, par exemple du type automobile, motocycle ou autres. Un ensemble 12 d'organes fonctionnels du véhicule est raccordé aux bornes A et B des moyens 10 de stockage d'énergie électrique, et comprend par exemple des moyens de motorisation, l'éclairage, le système de freinage, le système de direction du véhicule, etc... Cet ensemble 12 d'organes fonctionnels échange du courant électrique avec les moyens 10 de stockage pour son alimentation en énergie électrique ou pour la recharge des moyens 10 de stockage au moyen d'un alternateur des moyens de motorisation, comme cela est connu en soi dans l'état de la technique. Le système d'estimation de l'état de charge des moyens 10 de stockage électrique comprend des moyens 14 de sollicitation de ceux-ci selon un courant prédéterminé d'identification, et de manière avantageuse selon une séquence binaire pseudo-aléatoire de courant prédéterminée. Cette séquence pseudo-aléatoire de courant est déterminée pour permettre une identification efficace d'un modèle de la transmittance des moyens 10 de stockage d'énergie électrique, comme cela sera expliqué plus en détail par la suite. Le système d'estimation de l'état de charge des moyens 10 de stockage d'énergie électrique comprend par ailleurs des moyens 16 de mesure de la réponse en tension des moyens 10 de stockage d'énergie électrique à un tel courant d'identification. Le courant d'identification et la réponse en tension mesurée correspondante sont alors délivrés à des moyens 18 de prétraitement qui réalisent un prétraitement de ceux-ci à des fins d'identification. Ce prétraitement consiste par exemple à soustraire, du courant d'identification et de la réponse en tension mesurée correspondante, les composantes continues qu'ils contiennent. Le courant d'identification prétraité et la réponse en tension mesurée correspondante prétraitée sont délivrés à des moyens 20 d'identification. Les moyens 20 identifient alors, en fonction de ceux-ci, des paramètres d'un modèle H non entier prédéterminé de la transmittance des moyens 10 de stockage d'énergie électrique selon la relation: i., ,\-n, / ., \-n2 H(p) = K 1+ cor 1+ p cor i êM) (1) où p est la variable de Laplace, et n1, n2, K, wb, cor et coh sont les paramètres du modèle. De manière avantageuse, le modèle non entier selon la relation (1) modélise de manière efficace la transmittance des moyens 10 de stockage pour une plage de fréquences étendue, et permet par ailleurs de satisfaire un bon compromis entre simplicité et précision. L'identification des paramètres du modèle est réalisée en mettant en oeuvre un algorithme non linéaire d'identification du type erreur de sortie, par 15 exemple un algorithme du type Marquardt. Les paramètres identifiés du modèle sont ensuite délivrés à des moyens 22 d'estimation de l'état de charge SOC. Les moyens 22 estiment l'état de charge SOC des moyens 10 de stockage en évaluant, pour ces paramètres identifiés, une fonction prédéterminée d'estimation de l'état de charge dont les variables sont les paramètres du modèle. Cette fonction d'estimation est par exemple mémorisée dans les moyens 22 d'estimation sous la forme d'une table de correspondance ou les moyens 22 d'estimation comprennent des moyens d'interpolation adaptés pour déterminer l'état de charge SOC des moyens 10 de stockage d'énergie électrique en interpolant des valeurs mémorisées de cette fonction pour des valeurs prédéterminées des paramètres du modèle, comme cela est connu en soi dans l'état de la technique. De manière avantageuse, l'estimation de l'état de charge des moyens 10 de stockage est en outre réalisée en fonction de leur température de fonctionnement pour prendre en compte la variabilité de leur état de charge en fonction de leur température de fonctionnement. A cet effet, le système d'estimation comprend en outre des moyens 24 de mesure de la température des moyens 10 de stockage. La mesure de la température est délivrée à des moyens 26 de sélection qui sélectionnent, en fonction de celle-ci, une fonction d'estimation pour 5 les moyens 22 d'estimation parmi un ensemble prédéterminé de fonctions d'estimation, représentatives chacune d'une température prédéterminée de fonctionnement des moyens 10 de stockage. Les moyens 22 d'estimation peuvent également être adaptés pour estimer l'état de charge des moyens 10 de stockage en évaluant une fonction d'estimation dépendant, non seulement des paramètres du modèle, mais également de la température des moyens 10 de stockage. De manière avantageuse, la fonction prédéterminée d'estimation de l'état de charge est une combinaison de fonctions selon la relation: M SOC= M, (Q) (2) 1sr i=l i=l où M est un nombre prédéterminé supérieur ou égal à 1, Q est le vecteur des paramètres du modèle de la transmittance des moyens 10 de stockage, c'est-à-dire le vecteur (K, wb, wr, wh, n1, n2) T pour le modèle selon la relation (1), F;, i=1, ..., M, est une fonction des paramètres du modèle et 8i, i=1, M, est un indice de confiance prédéterminé associé à la fonction Fi. L'état de charge SOC des moyens 10 de stockage est alors estimé en évaluant la fonction selon la relation (2) pour le vecteur Q des paramètres égal au vecteur des paramètres identifiés par les moyens 22 d'identification. II va être maintenant décrit ci-dessous, en relation avec la figure 2, un procédé de calibrage du système d'estimation de l'état de charge des moyens 10 25 de stockage selon l'invention. Ce procédé consiste, dans une première étape 40, à sélectionner un nombre prédéterminé N de valeurs de l'état de charge des moyens 10 de stockage, par exemple SOCS, j=1,...,N. Le procédé consiste ensuite, dans une étape 42, pour chacune des valeurs de l'état de charge SOCi susmentionnées, à identifier les paramètres du modèle de la transmittance des moyens 10 de stockage, c'est-à-dire le vecteur (K, wb, wr, wh, n1, n2)T pour le modèle selon la relation (1). Il est ainsi obtenu un ensemble de valeurs de l'état de charge SOCi en fonction d'un ensemble de vecteurs Qi, j=1,...,N, respectifs des paramètres du modèle. L'identification mise en oeuvre dans cette étape 42 est du même type que celle mise en oeuvre par le système d'estimation selon l'invention. Elle consiste, pour la valeur de l'état de charge SOCi, à solliciter les moyens 10 de stockage selon un courant d'identification prédéterminé, de préférence la séquence binaire pseudo-aléatoire des moyens 14 de sollicitation, à mesurer la réponse en tension correspondante des moyens 10 de stockage, à prétraiter ce courant d'identification et la réponse en tension mesurée correspondante aux fins d'identification, puis à mettre en oeuvre l'algorithme non linéaire du type à erreur de sortie, par exemple du type Marquardt, pour obtenir le vecteur Qi correspondant des paramètres du modèle. Une fois l'identification des paramètres du modèle réalisée pour chacune des valeurs sélectionnées de l'état de charge, l'étape 44 suivante du procédé consiste à déterminer la fonction d'estimation de l'état de charge utilisée par les moyens 22 d'estimation. Cette détermination est réalisée en interpolant les valeurs sélectionnées de l'état de charge SOCS en fonction des vecteurs Qi correspondants des paramètres identifiés du modèle. Plus particulièrement, pour déterminer la fonction d'estimation selon la relation (2), l'étape 44 de détermination comprend une première étape 44a de détermination des M fonctions Fi. Pour la détermination de chacune de celles-ci, le procédé consiste tout d'abord à déterminer une première fonction Hi (Q) dépendant des paramètres du modèle. Cette première fonction Hi (Q) est déterminée telle que la variation des valeurs sélectionnées de l'état de charge SOCi en fonction des valeurs de la première fonction Hi(Q) pour les vecteurs Qi des paramètres identifiés correspondants du modèle est monotone. Un exemple d'une telle fonction Hi (Q) est illustré sur la figure 3a pour le modèle de transmittance des moyens 10 de stockage selon la relation (1) . La figure 3a représente la variation de 4 valeurs sélectionnées de l'état de charge SOCi, ici 40%, 60%, 80% et 100%, en fonction des valeurs correspondantes de la fonction Hi (Q) = K x con pour les vecteurs Qi = (K, wb, wr, 2877097 8 wn, n1, n2)iT des paramètres du modèle selon la relation (1) identifiés pour 40%, 60%, 80% et 100% respectivement. Comme il est possible de le constater, cette variation est monotone. D'autres exemples d'une telle première fonction sont illustrés sur les figures 3b à 3d, où la fonction H; (Q) est égale à wb, wr et n1 respectivement. Une fois déterminée la première fonction H; (Q) , le procédé consiste, toujours en 44a, à interpoler les valeurs sélectionnées de l'état de charge SOC; en fonction des valeurs de la première fonction H; (Q) pour les vecteurs Q; des paramètres identifiés correspondants afin d'obtenir une seconde fonction Gi (Hi (Q)). De manière préférentielle, cette interpolation est une interpolation monotone, par exemple polynomiale du type spline , ou par morceau comme cela est illustré sur les figures 3a à 3d, ou autre. Enfin, en remplaçant la variable de la seconde fonction par l'expression littérale de la fonction H; (Q) , il est ainsi obtenu la fonction F; = Gi(Hi(Q)) , c'est- à-dire la composée de la seconde fonction par la première fonction. L'étape 44b suivante de l'étape 44 de détermination de la fonction d'estimation des moyens 22 d'estimation consiste à déterminer, pour chacune des M fonctions Fi, un indice de confiance 8; . Un premier type de détermination de la pondération, ou indice de confiance 8; , de la fonction F; consiste à calculer une valeur représentative de la confiance en les résultats renvoyés par la fonction Fi. Par exemple, la pondération est calculée comme étant égale au rapport selon la relation: Qmax Qmin (3) Qmax +Qmin où désigne une norme quelconque, Qmin et Qmax sont des valeurs du vecteur des paramètres Q du modèle tels que la plage de valeurs [F; (Q min) F; (Q max)] de la fonction F; est égale à une plage prédéterminée d'états de charge [SOCmin SOCmax] Un second type de détermination de la pondération Si associée à la fonction Fi consiste à calculer celle-ci comme étant une valeur représentative de la sensibilité de la fonction Fi à des perturbations, comme par exemple des variations du courant moyen i de sollicitation utilisé lors de l'identification du modèle, des variations des caractéristiques des moyens de stockage d'énergie électrique, des variations des caractéristiques des sollicitations qui perturbent les moyens de stockage avant et/ou pendant l'estimation de l'état de charge de ceux-ci, ou autres. Par exemple, la pondération Si est déterminée comme étant dépendante de la moyenne i du courant d'identification, comme cela est illustré sur la figure 4. La figure 4 représente l'évolution de pondérations Si, associées respectivement à m fonction Fi, par exemple M=4, en fonction de la valeur de la moyenne i du courant utilisé pour l'identification des paramètres du modèle. La détermination consiste tout d'abord à déterminer une plage [i0 iM prédéterminée de courants moyens Ï dans laquelle chaque fonction fi présente un intervalle de stabilité maximale, c'est-à-dire un intervalle de courants dans lequel la sensibilité de la fonction fi à une variation du courant moyen i est minimale, et donc un intervalle où le résultat renvoyé par Fi est fiable. La détermination se poursuit ensuite par la division la plage [I0 IM] de courants moyens en M intervalles, chacun des M intervalles étant associé à la fonction Fi qui présente la stabilité maximale sur celui-ci. Par exemple, dans le cas où M=4 comme cela est illustré sur la figure 4, la plage [io iM] est divisée en quatre intervalles [i0 [i1 i2], [i2 i3] et [13 Iq a, la fonction FI étant la fonction parmi les fonctions Fi qui représente la stabilité maximale sur l'intervalle [ï0 il], la fonction F2 étant la fonction parmi les fonctions Fi qui présente la stabilité maximale sur l'intervalle [il 12], etc. La pondération Si associée à la fonction Fi, i=1,..., M, est alors une pondération selon la relation: Si =1 si i e [Ii-1 11 =0 sinon. En variante, comme cela est illustré sur la figure 4, la pondération 8; est une fonction affine par morceaux présentant un maximum au centre de l'intervalle [ii_1 il] et étant nulle en dehors de celui-ci, par exemple une 5 pondération selon les relations: - si i =1 alors 8i =1 si i Ti Ii-1 Ii-1 Ii+1 =A;xi+B; si IE 2 2-Ii =0 sinon - si i= M alors 8i= Ci xi+D; si iE Ii-1 - Ii-2 I; - Ii-1 2 2 si 1E =1 I; - Ii-1 i - sinon = 0 sinon 8i=CixI+D; si IE Ii-1 - Ii-2 Ti - Ii-1 2 2 i; Ii+1 - Ii =Ai xi+B; sil E -1i-1 2 2 = 0 sinon où A; = 1 avec = Ii+1 - 2Ii +'i-1 - Ali 2 B. = b' avec b; ; 0I; 2 I; -24_1+ I;-2 C;= i avec ü-1= Di1 2 - cl- Di; avec d. = = Ii-1 - Ii-2 Di; 2 Comme on peut le constater, des fonctions F; affines par morceaux 20 telles qu'illustrées sur la figure 4 ont pour effet de lisser les transitions entre chaque fonction F; de sorte que les discontinuités dans l'estimation de l'état de charge des moyens de stockage sont minimisées. et Bien entendu, si les pondérations 8i dépendent d'un paramètre comme le courant moyen 1 par exemple, les moyens 22 d'estimation du système décrit en relation avec la figure 1 reçoivent en outre la valeur de ce paramètre de moyens d'acquisition ou de détermination appropriés Par exemple, dans le cas où les pondérations dépendent du courant moyen 1 utilisé pour l'identification des paramètres du modèle, les moyens 22 d'estimation reçoivent la valeur de ce courant moyen 1 des moyens 20 d'identification qui calculent cette valeur pour leurs propres besoins, et les moyens 22 sont propres à calculer les valeurs des pondérations en fonction de la valeur reçue de courant. En variante, chaque pondération Si est déterminée comme étant à la fois représentative de la confiance et de la sensibilité de la fonction Fi qui lui est associée, par exemple en calculant cette pondération 8i comme étant la multiplication de la valeur de pondération calculée selon le premier type de 15 détermination décrit ci-dessus par la valeur de pondération calculée selon le second type de détermination décrit cidessus. La combinaison selon la relation (2) des fonctions fi permet de limiter la dispersion de l'estimation mise en oeuvre par le système d'estimation selon l'invention, et présente donc une fiabilité accrue dans les résultats d'estimation 20 qu'elle renvoie. Enfin, en se référant à nouveau à la figure 2, le procédé consiste ensuite dans une étape 44c à combiner les M fonctions Fi susmentionnées avec leur indice de confiance pour obtenir la fonction d'estimation des moyens 22 d'estimation selon la relation (2). Les différentes étapes du procédé qui viennent d'être décrites sont mises en oeuvre pour une température prédéterminée de fonctionnement des moyens 10 de stockage, afin d'obtenir une fonction d'estimation représentative de cette température. Bien entendu, une fonction d'estimation est déterminée par le procédé 30 susmentionné pour chaque température de l'ensemble prédéterminé de températures décrit précédemment en relation avec la figure 1. Bien entendu d'autres modes de réalisation sont possibles, comme par exemple des moyens d'estimation utilisant d'autres types de modèles de transmittance des moyens de stockage d'énergie électrique, ou des modèles non entier plus complexes ou plus simples. REVENDICATIONS 1. Système d'estimation de l'état de charge de moyens (10) de stockage d'énergie électrique de véhicule, caractérisé en ce qu'il comprend: - des moyens (14) de sollicitation des moyens (10) de stockage d'énergie électrique selon un courant prédéterminé d'identification; - des moyens (16) de mesure de la réponse en tension des moyens (10) de stockage d'énergie électrique à ce courant d'identification; - des moyens (20) d'indentification de paramètres d'un modèle prédéterminé de la transmittance des moyens (10) de stockage d'énergie électrique en fonction du courant d'identification et de la réponse en tension des moyens (10) de stockage d'énergie électrique à celui-ci; et - des moyens (22) d'estimation de l'état de charge des moyens (10) de stockage d'énergie électrique en évaluant, pour les paramètres identifiés, une fonction prédéterminée d'estimation de l'état de charge dont les variables sont les paramètres du modèle. 2. Système selon la revendication 1, caractérisé en ce que les moyens (14) de sollicitation sont adaptés pour solliciter les moyens (10) de stockage d'énergie électrique selon une séquence binaire pseudo-aléatoire de courant prédéterminée. 3. Système selon la revendication 1 ou 2, caractérisé en ce que le modèle de la transmittance des moyens (10) de stockage d'énergie électrique est un modèle non entier prédéterminé. 4. Système selon la revendication 3, caractérisé en ce que le modèle non entier est un modèle H non entier selon la relation: ( \-n1 / \-nz 1+p 1+p col) car 1+P 1+ p cors o)h où p est la variable de Laplace, et n l, n2, K, wb, wr et wh sont les paramètres du modèle. 5. Système selon la revendication 1, 2 ou 3, caractérisé en ce que les moyens (20) d'identification des paramètres du modèle de la transmittance sont 30 adaptés pour identifier les paramètres du modèle prédéterminé de la H(p) = K transmittance en mettant en oeuvre un algorithme d'identification du type erreur de sortie, notamment un algorithme du type Marquardt. 6. Système selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend en outre des moyens (24) de mesure de la température des moyens (10) de stockage d'énergie électrique et des moyens (26) de sélection de la fonction prédéterminée d'estimation parmi un ensemble prédéterminé de fonctions d'estimation de l'état de charge, en fonction de la température mesurée. 7. Procédé de calibrage d'un système d'estimation de l'état de charge de moyens (10) de stockage d'énergie électrique de véhicule conforme à l'une quelconque des revendications 1 à 6, caractérisé en ce qu'il comprend les étapes consistant à : - sélectionner un nombre N prédéterminé de valeurs de l'état de charge SOCS, j=1,...,N, des moyens (10) de stockage d'énergie électrique et identifier, pour chacune de ces valeurs SOCi, la valeur Qi du vecteur Q des paramètres du modèle prédéterminé de la transmittance des moyens (10) de stockage d'énergie électrique; et - déterminer la fonction d'estimation en interpolant les valeurs sélectionnées de l'état de charge SOCi en fonction des vecteurs Qi des 20 paramètres identifiés du modèle. 8. Procédé selon la revendication 7, caractérisé en ce que l'étape de détermination de la fonction d'estimation consiste à : - déterminer M fonctions Fi, i=1,...,M, des paramètres du modèle, M étant un nombre prédéterminé supérieur ou égal à 1; - déterminer, pour chacune des M fonctions Fi, un indice de confiance 8; dans les résultats renvoyés par la fonction fi; et - déterminer la fonction d'estimation de l'état de charge comme étant la combinaison des M fonctions Fi des paramètres selon la relation: M SOC = M1 L, 8i Fi (Q) Esi i=1 i=1 où SOC est l'état de charge des moyens de stockage d'énergie électrique. 9. Procédé selon la revendication 8, caractérisé en ce que l'étape de détermination des M fonctions F; des paramètres au modèle consiste à déterminer chacune des M fonctions F; en: - déterminant une première fonction H; (Q) des paramètres du modèle telle que la variation des valeurs sélectionnées de l'état de charge SOCi en fonction des valeurs de la première fonction H; (Q) pour que les vecteurs Qj des paramètres du modèle identifiés correspondants est monotone; - déterminant une seconde fonction Gi(Hi(Q)) en interpolant de manière monotone les valeurs sélectionnées de l'état de charge SOCS en fonction 10 des valeurs de la première fonction H; (Q) pour les vecteurs Qj des paramètres identifiés correspondants; et - déterminant la fonction des paramètres du modèle comme étant la composée de la seconde fonction par la première fonction. 10. Procédé selon la revendication 7 ou 9, caractérisé en ce qu'il est adapté pour déterminer la fonction d'estimation en fonction de la température.