Selon un cinquième aspect, l’invention concerne également un procédé, mis en œuvre dans un calculateur embarqué au sein d’un véhicule automobile, d’optimisation d’une combinaison d’au moins deux critères relatifs au véhicule et choisis parmi le groupe consistant en : un niveau de pollution, une consommation en énergie, un état de santé et une durée d’un trajet prédéterminé, le véhicule comprenant un réservoir de carburant ou d’hydrogène, une batterie électrique et/ou des super-capacités apte(s) à fournir de l’énergie électrique, un moteur thermique alimenté par le réservoir de carburant ou une pile à combustible alimentée par le réservoir d’hydrogène, au moins une machine électrique alimentée par de l’énergie électrique fournie par la batterie et/ou les super-capacités, au moins un dispositif relatif au moteur thermique ou à la pile à combustible, au moins un dispositif relatif à la machine électrique, et au moins un dispositif relatif à la batterie électrique ou aux super-capacités, le calculateur étant configuré pour gérer la chaine de traction du véhicule automobile sur ledit trajet prédéterminé et étant apte à commander le moteur thermique ou la pile à combustible, la machine électrique et/ou les dispositifs par l’émission d’un ensemble de consignes, le moteur thermique ou la pile à combustible, la machine électrique, et les dispositifs étant chacun caractérisés par au moins une variable d’état, chaque variable d’état permettant de décrire l’état de fonctionnement du dispositif qu’elle caractérise, l’ensemble du moteur thermique ou de la pile à combustible, de la machine électrique, et des dispositifs étant représenté par un système d’équations d’état modélisant la dynamique du véhicule, lesdites équations d’état étant fonctions au moins des valeurs de consignes instantanées et des variables d’état, la combinaison de critères à optimiser étant représentée par une équation de critère fonction au moins des valeurs de consignes instantanées et des variables d’état, le véhicule comprenant en outre des moyens de mesure du ou des état(s) décrit(s) par la ou chaque variable d’état et des moyens de mesure de chacun des critères de la combinaison de critères à optimiser, ledit procédé étant mis en œuvre pendant une durée divisée en instants d’échantillonnage constants, le procédé comprenant, à chaque instant d’échantillonnage, les étapes de : gg) détermination d’un domaine de consignes applicables comprenant un ensemble de valeurs pour chaque consigne, hh) calcul, pour chaque consigne et dans le domaine de consignes applicables déterminé, de tous les gradients d’état possibles pour ladite consigne, en utilisant au moins les équations d’état, ii) calcul, pour chaque consigne et dans le domaine de consignes applicables déterminé, de toutes les valeurs possibles pour la combinaison de critères à optimiser, en utilisant au moins l’équation de critère, jj) calcul, pour chaque consigne et dans le domaine de consignes applicables déterminé, de toutes les valeurs possibles d’une fonction Hamiltonienne de ladite consigne, en utilisant les valeurs de gradients d’état et les valeurs possibles pour la combinaison de critères à optimiser déterminées pour ladite consigne, kk) détermination de la valeur de chaque consigne, dans le domaine de consignes applicables déterminé, pour laquelle ladite fonction Hamiltonienne est la plus faible, ll) calcul, pour chacune des consignes, en fonction de ladite valeur déterminée de la consigne, de la valeur du gradient d’état pour ladite consigne, et/ou calcul, pour chacune des consignes, en fonction de ladite valeur déterminée de la consigne, de la valeur de chaque critère de la combinaison de critères pour ladite consigne, mm) intégration temporelle de la valeur du gradient d’état calculée, fournissant une valeur d’état estimée, et/ou intégration temporelle, pour chaque critère de la combinaison de critères, de la valeur calculée de ce critère, fournissant une valeur de critère globale estimée, nn) comparaison de la valeur d’état estimée à une valeur d’état mesurée, l’écart entre les deux valeurs fournissant un premier terme de correction, le premier terme de correction ainsi déterminé étant fourni en entrée de l’étape de calcul hh) à l’instant d’échantillonnage suivant et étant utilisé dans ladite étape de calcul hh), et/ou comparaison, pour chaque critère de la combinaison de critères, de la valeur de critère globale estimée à une valeur de critère mesurée, l’écart entre les deux valeurs fournissant un second terme de correction, le second terme de correction ainsi déterminé étant fourni en entrée de l’étape de calcul ii) à l’instant d’échantillonnage suivant et étant utilisé dans ladite étape de calcul ii).