A cette fin et selon un premier aspect, l’invention concerne un procédé, mis en œuvre dans un calculateur embarqué au sein d’un véhicule automobile, d’optimisation d’un niveau de pollution du véhicule, 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 un 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, le niveau de pollution à optimiser étant représenté 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 du niveau de pollution du véhicule, 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 : a) détermination d’un domaine de consignes applicables comprenant un ensemble de valeurs pour chaque consigne, b) 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, c) calcul, pour chaque consigne et dans le domaine de consignes applicables déterminé, de toutes les valeurs possibles pour le niveau de pollution à optimiser, en utilisant au moins l’équation de critère, d) 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 le niveau de pollution à optimiser déterminées pour ladite consigne, e) 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, f) 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 du niveau de pollution pour ladite consigne, g) intégration temporelle de la valeur du gradient d’état calculée, fournissant une valeur d’état estimée, et/ou intégration temporelle de la valeur du niveau de pollution calculée, fournissant une valeur de niveau de pollution globale estimée, h) 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 b) à l’instant d’échantillonnage suivant et étant utilisé dans ladite étape de calcul b), et/ou comparaison de la valeur de niveau de pollution globale estimée à une valeur de niveau de pollution 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 c) à l’instant d’échantillonnage suivant et étant utilisé dans ladite étape de calcul c).