FR 3032919 A1 20160826 FR 1551400 A 20150219 1 Procédé et dispositif de détection d'un changement de comportement de conducteur d'un véhicule automobile. L'invention concerne un procédé et un dispositif de détection d'un changement de comportement de conducteur d'un véhicule automobile. Le comportement du conducteur désigne ici son état comportemental dont peut résulter, une façon de se comporter, en d'autres termes de réagir aux événements extérieurs. Un état de vigilance forte se traduit généralement par des facultés notamment d'anticiper et d'agir promptement, progressivement et efficacement. Un état de fatigue, de distraction ou d'énervement se traduit par une baisse de facultés du conducteur qui peuvent avoir des répercussions en termes anxiogènes pour les autres usagers, voire en termes d'accidents. La détection d'un changement de comportement du conducteur est donc utile pour alerter le conducteur, voire pour agir directement sur le véhicule. Les constructeurs automobile ont pris conscience très tôt du besoin et ont proposé des procédés et dispositifs pour répondre à ce besoin, mais sans apporter entière satisfaction. Les procédés de détection de la baisse de vigilance du conducteur peuvent être classés en deux grandes familles, la première par approche basée sur des mesures directes et la seconde par approche indirecte. Dans l'approche basée sur des mesures directes, on peut citer : - Les mesures de l'état physiologique du conducteur : fréquence cardiaque (électrocardiogramme), capteurs de l'activité cérébrale (électroencéphalogramme entre autres), l'électro-oculogramme, comme par exemple dans le document EP0699413. - Les mesures de l'attitude du conducteur basées sur l'observation du conducteur : clignement des yeux (Perclos qui caractérise le degré d'ouverture, fermeture des yeux), orientation et mouvement de la tête, posture du conducteur assis sur son siège, comme par exemple dans le document EP2368496, - Les mesures sur une réponse du conducteur à un événement préprogrammé comme par exemple dans les documents EP1957309 ou EP2512340. L'approche basée sur des mesures directes présente de nombreux inconvénients comme par exemple ceux de la nécessité d'installer un appareillage supplémentaire dans le 3032919 véhicule, parfois complexe et coûteux, ou encore ceux de se présenter particulièrement intrusifs sur la personne du conducteur. L'approche basée sur des mesures indirectes permet d'éviter les inconvénients de la première. On peut citer : 5 - Les mesures d'actions du conducteur sur le véhicule : angle volant, couple volant, pédale d'accélérateur, de frein, d'embrayage et changement de rapport, actions secondaires sur les clignotants, IHM audio, navigation, téléphonie main libre, climatisation, sans provocation particulière du conducteur par le procédé ou par le dispositif comme par exemple dans les documents US8593288 ou US6061610. 10 Mesure de la dynamique du véhicule : distance entre les véhicules, vitesse du véhicule, position du véhicule dans la voie, angle de cap, comme par exemple dans le document EP0952039. C'est sur cette seconde approche que l'invention porte sa préférence. Dans cette seconde approche, il est généralement nécessaire d'utiliser un jeu de 15 paramètres qui est réglé une fois pour toute dans le véhicule, lors de la phase de mise au point du véhicule, et ce, quelques soient les conditions d'utilisation et le type de conducteur. La difficulté de mise en oeuvre des procédés antérieurement connus dans le cadre de la seconde approche, consiste alors à trouver le bon jeu de réglage des paramètres 20 permettant de détecter une baisse de vigilance du conducteur dans des situations multiples et variées, comme celles liées au type de conducteur et à son type de conduite, aux modifications de configurations routières (autoroute / route départementale, trafic dense fluide) et aux modifications des caractéristiques du véhicule. Le compromis nécessaire pour satisfaire un maximum de situations possibles fait surgir l'inconvénient de ne pas avoir le bon 25 jeu de réglage des paramètres dans chaque situation prévue, lequel est amplifié pour des situations qui peuvent surgir de façon non prévues ou pour lesquelles on ne possède pas nécessairement tous les éléments au moment de la conception puis des essais. Pour remédier aux inconvénients de l'état antérieur de la technique, l'invention a pour objet un procédé de détection d'un changement de comportement de conducteur d'un véhicule automobile au moyen d'un modèle cybernétique défini par un premier 3032919 vecteur de mesure en relation avec au moins une entrée connectée à un instrument de mesure pour quantifier une donnée sensorirnotrice en relation avec le conducteur, par un vecteur de sortie comportant au moins une sortie pour estimer une donnée actionnelle produite par le conducteur, par un vecteur de paramétrage comportant au moins un 5 paramètre pour quantifier une donnée physiologique du conducteur, et par un vecteur d'état temporel comportant au moins un état temporel dudit modèle, caractérisé en ce qu'il comprend des étapes exécutées en temps réel pour : calculer au moins une valeur courante optimale d'un vecteur d'état augmenté à un instant présent en fonction d'au moins une valeur précédemment estimée dudit 10 vecteur d'état augmenté et d'un deuxième vecteur de mesure en relation avec ledit vecteur de sortie, ledit vecteur d'état augmenté comportant au moins une composante reproduisant une composante du vecteur d'état temporel et au moins une composante reproduisant une composante du vecteur de paramétrage ; évaluer un degré de stabilité d'au moins une composante de valeur courante du 15 vecteur d'état augmenté qui reproduit ladite au moins une composante du vecteur de paramétrage ; signaler un changement de comportement de conducteur lorsque ledit degré de stabilité est significatif d'une instabilité. Le modèle cybernétique du conducteur combiné à une identification en temps réel 20 des paramètres du modèle conducteur permet un diagnostic de la vigilance conducteur (distraction notamment) adapté à chaque conducteur et à chaque configuration de véhicule qu'il conduit. La variation au cours du temps de la valeur de ces paramètres, ou vitesse de variation des valeurs des paramètres du modèle, renseignent sur la modification du comportement du 25 conducteur (distraction), et dans une moindre mesure, de changement de l'état de la route. Par rapport à une technique d'identification classique telle que celles déjà divulguées dans la littérature, comme la méthode de prédiction d'erreur PEM (Prediction Error Method en anglais), la technique enseignée par l'invention, permet d'identifier la valeur instantanée du paramètre. Le diagnostic peut donc se faire sur un pas de temps très court (temps réel), 30 contrairement à la méthode PEM qui nécessite d'observer les signaux sur un intervalle de 3032919 temps plus long. L'avantage est d'obtenir des résultats plus proches de la réalité puisque le comportement du conducteur associé à la valeur des paramètres du modèle, varie sur des pas de temps très court. L'identification en temps réel permet de faire un diagnostic plus fidèle à la 5 dynamique rapide de changement du comportement du conducteur Particulièrement, le procédé comprend des étapes exécutées en temps réel pour calculer un gain de filtre bayésien à partir de la valeur précédemment estimée dudit vecteur d'état augmenté et dudit premier vecteur de mesure, le gain de filtre bayésien constituant un moyen technique de calcul de valeur courante du vecteur d'état augmenté. 10 De préférence, le procédé comprend des étapes exécutées en temps réel pour : générer un premier ensemble prédictif de sigma points associés au vecteur d'état augmenté, centrés sur la valeur précédemment estimée dudit vecteur d'état augmenté au moyen d'une matrice prédictive de covariance d'états statistiquement possibles; 15 calculer un deuxième ensemble prédictif de sigma points associés au vecteur de sortie à l'instant présent en fonction du premier ensemble prédictif de sigma points et dudit vecteur de mesure; calculer le gain de filtre bayésien à partir du deuxième ensemble prédictif de sigma points. 20 Particulièrement aussi, le procédé comprend des étapes exécutées en temps réel pour calculer une valeur future estimée dudit vecteur d'état augmenté à un instant suivant à partir de la valeur courante optimale dudit vecteur d'état augmenté et dudit premier vecteur de mesure. De préférence, le procédé comprend des étapes exécutées en temps réel pour : 25 générer un premier ensemble optimal de sigma points associés au vecteur d'état augmenté, centrés sur la valeur courante optimale dudit vecteur d'état augmenté au moyen d'une matrice optimale de covariance d'états statistiquement possibles; calculer un deuxième ensemble optimal de sigma points associés au vecteur d'état augmenté à l'instant suivant en fonction du premier ensemble optimal de sigma 30 points et dudit vecteur de mesure; 3032919 calculer la valeur future estimée dudit vecteur à partir du deuxième ensemble optimal de sigma points. Avantageusement, pour évaluer ledit degré de stabilité, le procédé comprend des étapes consistant à : 5 obtenir une valeur du vecteur de paramétrage à partir du vecteur d'état augmenté à l'instant présent ; calculer une évolution de paramètres à partir de la valeur à l'instant présent et d'une valeur à l'instant précédent du vecteur de paramétrage ; mémoriser dans un registre à décalage ladite évolution de paramètres de 10 manière à obtenir une quantité d'évolutions antérieures de paramètres ; calculer un écart type sur ladite quantité d'évolutions ; signifier une instabilité dudit degré de stabilité lorsque ledit écart type n'est pas inférieur à un seuil prédéterminé. L'invention a aussi pour objet un programme d'ordinateur comprenant des 15 instructions de code de programme pour l'exécution des étapes du procédé lorsque ledit programme est exécuté sur un ordinateur. L'invention a encore pour objet un dispositif de détection d'un changement de comportement de conducteur d'un véhicule automobile comprenant une unité de post traitement connectée à un ou plusieurs instruments de mesure pour quantifier une ou 20 plusieurs données sensorimotrices en relation avec le conducteur, une unité de calcul connectée à l'unité de post traitement et hébergée dans un ordinateur embarqué qui comporte en mémoire un modèle cybernétique défini par un vecteur de mesure comportant au moins une entrée, par un vecteur de sortie comportant au moins une sortie pour estimer une donnée actionnelle produite par le conducteur, par un vecteur de 25 paramétrage comportant au moins un paramètre pour quantifier une donnée physiologique du conducteur, et par un vecteur d'état temporel comportant au moins un état temporel dudit modèle, caractérisé en ce qu'il comprend en mémoire un programme selon l'invention et un processeur pour exécuter en temps réel ledit programme. L'invention a enfin pour objet un véhicule comprenant un dispositif selon 30 l'invention. 3032919 L'invention sera mieux comprise à l'aide d'exemples de réalisation d'un procédé et d'un dispositif conformes à l'invention en référence aux dessins annexés, dans lesquels : la figure 1 est une vue d'un véhicule auquel l'invention est applicable ; la figure 2 est une vue schématique du véhicule de la figure 1 évoluant sur une route ; la figure 3 est une vue schématique d'un dispositif particulier de mise oeuvre de l'invention ; la figure 4 est une vue schématique d'étapes de procédé conforme à l'invention ; 10 - la figure 5 illustre un modèle cybernétique de conducteur ; les figures 6 à 8 montrent des étapes détaillées d'une réalisation possible de procédé conforme à l'invention. La figure 1 montre un véhicule 12 équipé de différents capteurs et actionneurs non nécessairement représentés pour conduire le véhicule de manière connue en soi par 15 ailleurs. Une ou plusieurs caméras connectées à un calculateur embarqué 5, permettent de visualiser l'environnement du véhicule de manière comparable à ce que peut visualiser le conducteur à travers les vitres du véhicule ou pour aider le conducteur en relation avec des images qu'il ne peut pas visualiser directement. Le calculateur embarqué 5 est connecté à d'autres calculateurs au moyen d'un bus 6 de type CAN, LIN, Ethernet automobile ou de 20 tout autre type connu dans le domaine des réseaux embarqués. Un calculateur embarqué 1 contrôle et commande différents actionneurs tels que par exemple le groupe motopropulseur et les freins en liaison avec les pédales d'accélération et de freinage à disposition du conducteur et/ou des capteurs de vitesse 2 disposés sur tout ou partie des roues 21, 22 du véhicule. Un calculateur embarqué 7, souvent nommé ordinateur de bord, 25 est connecté au calculateur embarqué 1 et à d'autres calculateurs et capteurs au moyen d'un bus 8, souvent nommé bus véhicule, de type CAN, LIN, Ethernet automobile ou de tout autre type connu dans le domaine des réseaux embarqués. Parmi les signaux de capteurs transmis sur le bus 8, on peut citer une vitesse v du véhicule mesuré à partir du 3032919 ou des capteurs de vitesse 2, un angle de rotation Sd et un couple rd appliqués par le conducteur à un volant de direction 4. Le dispositif de détection d'un changement de comportement de conducteur du véhicule 12 selon l'invention concerne notamment un véhicule automobile dans le sens où le comportement du conducteur ou du pilote a une influence directe sur celui du véhicule. Ce peut être aussi bien un véhicule nautique ou aéronautique qu'un véhicule terrestre. Lorsque le véhicule 12 se déplace par exemple sur une route, le conducteur apprécie généralement un angle de visée prochee',' et un angle de visée lointainefar, comme illustré par la figure 2. L'angle de visée proche est l'angle du champ de vision sous lequel le 10 conducteur perçoit généralement la ligne médiane 14 de la chaussée. L'angle de visée lointain est l'angle fait entre la tangente à la ligne de direction du véhicule et la tangente à la courbure interne 13 de la chaussée. Lorsque la chaussée est rectiligne, la tangente à la courbure interne 13 est naturellement reportée à l'infini sur le point de fuite des bords de la chaussée. De manière plus générale, l'angle de visée lointainOfar correspond à 15 l'anticipation de la courbure de la route par le conducteur et l'angle de visée procheOnear correspond au contrôle de position latérale effectué par le conducteur en se repérant sur un point proche. L'exemple de réalisation du dispositif de détection d'un changement de comportement de conducteur selon l'invention illustré par la figure 3, est hébergé de 20 manière centralisée dans le calculateur embarqué 7 ou de manière répartie dans plusieurs calculateurs 1, 5, 7. Le dispositif comprend notamment une unité 10 de post-traitement connectée aux bus 6, 8 pour y lire les signaux en provenance des capteurs 2, 3 et de la ou des caméras pour élaborer un vecteur U de mesures d'entrées et un vecteur Z de mesures de sorties en 25 relation avec un modèle cybernétique expliqué dans la suite de la description en référence à la figure 5. On rappelle qu'au sens du traitement informatique, un vecteur de mesures est une liste ordonnée de variables dont la valeur de chacune est positionnée à partir du bus de données ou directement à partir d'un capteur. Le dispositif comprend aussi une unité 11 de calcul connectée à l'unité 10 de post-traitement pour recevoir le vecteur U de mesures 30 d'entrées et le vecteur Z de mesures de sorties, de manière à générer une alarme 9 sous 3032919 forme sonore, visuelle ou sou recq me de message pour un signaler un changement de comportement du conducteur. L'unité 11 de calcul est par exemple mais non nécessairement hébergée dans le calculateur ou ordinateur embarqué 7 qui comprend un processeur et, en mémoire un programme et des données pour l'exécution des étapes de procédé décrites ci-après en référence notamment à la figure 4. Des données variables d'entrée-sortie y sont réservées pour reproduire en temps réel le vecteur U de mesures d'entrées et le vecteur Z de mesures de sorties. Des données variables internes au programme y sont réservées pour mémoriser temporairement à chaque instant k ou sur une plage de temps déterminée, les 10 valeurs nécessaires à l'exécution des étapes. Des données constantes y sont réservées pour mémoriser de manière permanente, les valeurs constantes nécessaires à l'exécution des étapes, notamment d'une étape d'initialisation 100. La figure 5 illustre, à titre de pur exemple non limitatif, un modèle cybernétique possible du conducteur basé sur les connaissances des fonctions sensorimotrices de l'être 15 humain. Le modèle cybernétique ici représenté comprend quatre entrées et une sortie principale. En entrées, le couple d'autoalignement I-, qui traduit le couple exercé par le train avant sur la colonne de direction, l'angle de rotation du volant Sd , Ofar de visée à l'avant lointain du véhicule, et l'angle Onear de visée à l'avant proche du véhicule, sont des données perceptibles par le conducteur, être humain pour agir sur le véhicule en 20 fonction de son comportement. Le dispositif de détection de modification de comportement du conducteur, notamment l'unité de calcul 11 reçoit ces mêmes données à partir de mesures en temps réel comme expliqué en référence à la figure 3. En sortie principale, le couple conducteur rd est le couple que le conducteur est sensé appliquer sur le volant de direction en fonction des données d'entrée qu'il est sensé percevoir. Le 25 dispositif de détection de modification de comportement du conducteur peut recevoir par ailleurs en temps réel une valeur de mesure du couple que le conducteur exerce effectivement sur le volant de direction, comme expliqué en référence à la figure 3. Dans une variante de réalisation du dispositif de détection de distraction du conducteur, on peut aussi rajouter en sortie secondaire l'intention angulaire du 3032919 conducteur de tourner le volant de direction, dans le but de Vouver des paramètres stabilisant lors d'une phase d'identification. Sur la figure 5, la structure interne du modèle cybernétique comprend essentiellement quatre modules. Le module 20 modélise une anticipation visuelle au point 5 tangent de la trajectoire en multipliant l'angle Ofat. par un facteur de gain proportionnel Kp. Le module 21 modélise une compensation visuelle au point de visé devant le véhicule en divisant l'angle 0''' par la vitesse instantanée du véhicule et en le multipliant par un facteur de gain proportionnel Kc. à constante de tempsTi. Un module 22 introduit un délai à la somme des résultats obtenus en sortie des modules 20 et 21 pour intégrer le retard de 10 traitement Tp par le cerveau humain de façon à obtenir l'intention angulaire 85 du conducteur de tourner le volant de direction. Le module 23 modélise le système neuromusculaire humain qui transforme la consigne issue des traitements précédents en signal du couple f'd en tenant compte des retours proprioceptifs que sont le couple d'autoalignement I-, et l'angle volant 6d. 15 Plus précisément, le module 23 obtient l'image d'un couple utile pour atteindre la valeur d'intention angulaire 45'' du conducteur en sommant deux termes. Un premier terme correspond à la valeur d'intention angulaire 6, du conducteur multipliée par la vitesse v du véhicule et par un coefficient de proportionnalité K. Un deuxième terme correspond à l'écart entre la valeur d'intention angulaire 85' du conducteur et l'angle 20 volant 6d scruté à chaque instant en temps réel, multiplié par un coefficient de proportionnalité Kt. Le couple d'autoalignement r, est alors retranché au couple utile de façon à obtenir une différence de couple qui est égale au signal du couple rd à une constante de temps Tp près. Les paramètres Kp K, T1 Tp Kr. Kt Tt, du modèle sont regroupés en mémoire dans un 25 vecteur n Kp K, Ti Tp Kr Kt Tp qui est représentatif du comportement physiologique du conducteur. Les valeurs des paramètres qui viennent d'être mentionnés, sont généralement constantes pour un conducteur considéré dont elles quantifient le comportement. Toute modification constatée d'une ou plusieurs valeurs de ces paramètres, traduit un changement de comportement du conducteur, en d'autres termes 30 un changement de façon de réagir aux évènements. Le changement de comportement 3032919 10 peut résulter d'une baisse de vigilance due par exemple à la fatigue, d'une distraction causée par un panneau publicitaire, un appel téléphonique ou par une recherche de station radio. Le changement de comportement peut résulter aussi d'une adaptation à des conditions particulières de circulation, par exemple dues à un mauvais état de la chaussée. Une particularité de l'invention réside dans la définition d'un vecteur à valeur variable n, de paramétrage à chaque instant k d'exécution du procédé. L'unité de calcul 11 peut comprendre en mémoire une valeur constante no par défaut du vecteur de paramétrage du modèle, correspondant par exemple à des valeurs moyennes de paramètres obtenues préalablement sur un échantillon de conducteurs représentatifs. 10 L'objet du dispositif de détection de distraction du conducteur selon l'invention est de détecter cet éventuel changement de comportement du conducteur. Pour cela le dispositif de détection simule le modèle expliqué ci-dessus au moyen d'un vecteur d'état x = [x1 x2 x3]T. Chaque composante du vecteur d'état x(t) à l'instant t, a la valeur d'une variable d'état temporel du modèle : 1 + xi 15X3 Kt8d Krv Kt) Tr,s-Ft On rappelle qu'une description mathématique du modèle de la figure 5, ou de tout autre modèle cybernétique, peut se faire en considérant une représentation linéaire de la forme : fk f[x(t), u(t), H] , x(to) y = g[x(t), u(t), , z(k) = y(k) [i (k), k = 1, N Avec x = [x1 x2 x3]T le vecteur d'état de valeur initiale xo Les vecteurs 20 d'entrée et respectivement de sortie du modèle u = Mar 81ear rs et y = [rd SS,,,1T représentent respectivement le vecteur des entrées et le vecteur des sorties, TI = [Kp Kr T1 Tp Kr Kt est le vecteur des paramètres du modèle considérés ici inconnus ou en d'autres termes des paramètres dont on cherche à identifier les valeurs. Les fonctions réelles f et g donnent respectivement d'une part une dérivée k du vecteur d'état par 25 rapport au temps lorsque les paramètres sont considérés constants et d'autre part une 3032919 11 valeur y estimée en sortie du modèle. Le vecteur z des sorties mesurées à chaque période d'échantillonnage pour laquelle k est l'index du temps discret, prend en compte un bruit de mesure qui a pour matrice de variance covariance R. Dans un souci d'alléger la description on nomme matrice de covariance toute matrice de variance covariance sachant que sur la diagonale de la matrice, la covariance exprime la variance. La matrice R ou Ra est à coefficients invariants dans le temps. En considérant les variables d'états relatives au modèle de la figure 5 : Kc 1 s + 1 "car 1+ Tns+1 1 - Kt5d + + K ) La dérivée des variables d'état permet d'obtenir une expression analytique de la 10 fonction f xl +vTi near K f[x(t), u(t), H 2 2 -X1-TA X2 + 2 - Tp Tp -(Krv+Kt) (Krv-I-Kt) 1 2 x2 X Tn T/1 A partir de la structure du modèle de la figure 5 et des variables d'état définies ci-dessus, les sorties du modèles sont estimées comme suit : et 5, + 2x2 - KpOf' 15 L'expression analytique de la fonction g est déduite des équations des sorties : x3 + 2x2 KpOfm. Dans les rappels mathématiques présentés ci-dessus pour une meilleure compréhension de la suite de la description, le vecteur n des paramètres, est supposé connu et constant. Un système de base, linéaire en ses paramètres, vient ainsi d'être défini 20 à partir du modèle de la figure 5. Il sert à estimer le vecteur d'état x = ixi x2 x3]T duquel sont déduites des estimations des mesures. Le système ici élaboré à partir d'un modèle cybernétique particulier qui, par ailleurs ne fait pas partie en soi de l'invention, sera bien entendu différent pour un autre modèle cybernétique avec un nombre différent d'entrées, un nombre différent de paramètres, un nombre différent d'états et un nombre différent de g[x (0, u (t), 3032919 12 sorties. L'invention s'applique tout aussi bien à un simple régulateur PID (Proportionnel, Intégral, Dérivé). La discrétisation à la période d'échantillonnage T est simplement effectuée par la méthode d'Euler ; en posant xk = x(kT): xk+1 = fd [xk, uk, Yk = g[xk, uk, Où fd (kk, uk, rd) k T f(kk, uk, Tri) Dans le but d'estimer aussi le vecteur n des paramètres, l'invention transforme le système de base, comme par exemple celui évoqué ci-dessus en un système augmenté. Le système augmenté associé à la mise en oeuvre du modèle, est obtenu en augmentant la dimension du vecteur d'état X du système de base avec de nouvelles 10 variables d'état extraites du vecteur paramètre FI = [ Ke T1 Tp Kr Kt T. ni ; Avec x = X2 T (variables d'état du système initial) x2 x K K, Tp Kt TidT En considérant que le vecteur paramètre n varie très peu, la dynamique du système augmenté est représentée par les équations suivantes: 15 fa [ (t), u(t)i + a (t) , xao = xa(ta) y = ga[xa (t), u(0], z(k) = y(k) + µa (k), k = 1, . , N 0(7,1) est une façon symbolique de représenter un vecteur à sept composantes nulles a(0, u(t)1 0(7,1) Avec pa bruit de mesure du système augmenté de covariance Ra Nous venons ainsi de définir les éléments du modèle utiles pour mettre en oeuvre du procédé expliqué maintenant en référence à la figure 4. 0(7,1) fa [xa (0,11(01 Et ga = g 3032919 13 L'étape 100 d'initialisation consiste à affecter pour un instant initial d'indice k=0, des valeurs constantes contenues en mémoire aux variables internes du procédé amenées à évoluer à chaque instant k ultérieur. Une valeur X.,k+i du vecteur d'état augmenté, prédite pour l'instant suivant, est 5 initialisée par la formule x,,,k+1 xa_o où x est un vecteur à L composantes par exemple toutes égales à 0,5. A titre purement illustratif et non limitatif pour le modèle exposé ci-dessus avec un vecteur d'état de base à trois composantes et un vecteur de paramètres à sept composantes, L.=-10. On peut donner aussi aux sept dernières composantes initiales du vecteur des paramètres, des valeurs préalablement définies expérimentalement. 10 La valeur Pa k+1 de la matrice de covariance augmentée de la prédiction de l'état du système à l'instant k+1, notamment en statistique d'ordre 2, est initialisé par la formule Pa k+1 Pau OÙ Pan est une matrice à LxL composantes par exemple toutes égales à 0,1 sur la diagonale à une ou deux exceptions près basées sur l'expérience, et toutes nulles en dehors de la diagonale. 15 La valeur Ra k de la matrice de covariance augmentée du bruit de mesure est initialisé par la formule Rajk W Rai) où, pour un vecteur de sortie à deux composantes, Ra_o est une matrice à 2x2 composantes par exemple toutes inférieures à 0,001 sur la diagonale en se basant sur l'expérience et/ou les spécifications des organes de mesure, et toutes nulles en dehors de la diagonale. On peut prendre aussi une valeur constante Ra pour la matrice de 20 covariance augmentée du bruit de mesure qui est alors accédée directement en mémoire au cours de l'exécution des étapes sans besoin de la positionner dans l'étape d'initialisation. Les autres grandeurs qui gardent une valeur constante au cours du temps sont simplement stockées en mémoire permanente sans avoir à être initialisées dans l'étape 25 100. Le nombre L de composantes du vecteur d'état augmenté est égal à la somme du nombre d'états temporels du modèle cybernétique et du nombre de paramètres du modèle ajoutés. Un facteur d'échelle À a été déterminé de manière heuristique pour optimiser les temps de calcul, notamment en termes de convergence, la valeur retenue est À=3-L, d'autres valeurs de À peuvent être retenues sans sortir du cadre de l'invention. Le 3032919 14 facteur d'échelle X sert entre autre mais non uniquement comme nous le verrons par la suite, à calculer des coefficients de pondération gliP, variant de 0 à 21, au moyen des formules suivantes : Pour i=0, Wo À/(L-1-À) et pour i variant de 1 à 2L, TV,P .1/2(L-R). Les coefficients de 5 pondération dont la somme est égale à l'unité, sont des constantes qui peuvent être stockées en mémoire permanente ou calculées dans l'étape d'initialisation à partir du facteur d'échelle À. Le facteur d'oubli M est une constante positive qui traduit un nombre d'itérations antérieures prises en compte et qui est utilisé pour renforcer la confiance sur les 10 prédictions. L'ajout des paramètres dans le vecteur d'état augmenté peut avoir pour effet d'introduire, sauf à de rares exceptions près, des non linéarités dans le système augmenté, notamment dans celui évoqué ci-dessus du fait des multiplications qui en résultent. L'étape d'expansion 110 résout le problème de non linéarité en utilisant une 15 méthode de calcul associée à la statistique en relation avec une variable aléatoire soumise à des évolutions non linéaires, connue sous le nom de transformée «Unscented» ou transformée sans odeur en français, peut-être en raison de son caractère éphémère. On peut se dispenser de l'étape 110 lorsque l'ajout des paramètres du modèle aux états temporels du modèle pour obtenir le vecteur d'état augmenté, conserve un traitement 20 Le vecteur d'état augmenté est vu dans l'étape 110 comme une variable aléatoire Xa de dimension L, de moyenne Xa et de covariance Pa, qui subit une transformation non linéaire de la forme Y = g(Xa) . La transformée unscented permet d'approcher la moyenne Y, la covariance Pyy de Y, ainsi que la covariance Pxy croisée de X8 et de Y. 25 On peut détailler l'étape 110 en trois étapes 111, 112, 113. L'étape 111 consiste essentiellement à affecter aux valeurs courantes à l'instant k, .."1',,{ du vecteur d'état augmenté prédit et P de la matrice de covariance de la prédiction d'état, les valeurs .Y,,,k+1 du vecteur d'état augmenté précédemment prédit et Pa,k+1 de la matrice de covariance de la prédiction d'état précédemment considérées comme étant 3032919 15 celles de l'instant k+1 suivant, soit dans l'étape initiale 100 lors de la première itération de l'étape 110, soit à la suite d'une itération précédente de l'étape 110 pour les itérations suivantes. L'étape 112 consiste essentiellement à déterminer une racine de 1') au sens de Il, de vecteur de sortie permettent de calculer une matrice P', k de covariance croisée entre les états donnés par le premier ensemble prédictif de sigma points et les valeurs de sortie YA du modèle dans une étape 124, et de calculer une matrice P13, Â de covariance sur les valeurs de sortie )1 du modèle dans une étape 125. En relation avec l'une quelconque des étapes 124, 125, pour chaque point du premier ensemble prédictif de sigma points, on génère un vecteur d'écart d'état relatif au vecteur d'état augmenté xa,A prédit à L composantes et un vecteur d'écart de sortie [11 - yk] relatif à la valeur estimée ) à Q composantes. 10 Dans l'étape 124, chaque produit du vecteur d'écart d'état par la transposée du vecteur d'écart de sortie correspondant donne une matrice à L lignes et Q colonnes. La somme pondérée par les coefficients w( des 21+1 matrices ainsi obtenues donne la matrice Pt,, k de covariance croisée entre états et sorties. Dans l'étape 125, chaque produit du vecteur d'écart de sorties par la transposée du 15 vecteur d'écart de sortie correspondant donne une matrice à Q lignes et Q colonnes. On ajoute alors la matrice Ra de covariance de bruit de mesure des sorties à la somme pondérée par les coefficients I1" des 21+1 matrices ainsi obtenues pour donner la matrice de covariance croisée sur les sorties. L'étape 126 consiste alors à générer une valeur de gain de Kalman Kk à l'instant k en 20 multipliant la matrice P k de covariance croisée entre états et sorties par l'inverse de la matrice Pvi, t de covariance sur les sorties. Notons à ce sujet que le gain de Kalman est ici une matrice à L lignes et à Q colonnes. L'étape 130 consiste essentiellement à recadrer le vecteur d'état augmenté sur une valeur optimisée pour conduire aux valeurs de sortie effectivement mesurées, et ceci au 25 moyen du gain de Kalman généré par l'étape 120. On peut détailler l'étape 130 en deux étapes 131, 132. 3032919 18 L'étape 131 consiste es.5entieliement à lire le vecteur zk des sorties du modèle telles qu'elles sont effectivement mesurées à l'instant courant k. Dans l'étape 132, la différence à l'instant k entre le vecteur zk des sorties mesurées et la valeur estimée yx du vecteur des sorties à Q composantes, est multiplié par le gain de Kalman Kk à l'instant k dont le résultat qui est un vecteur à L composantes, est ajouté au vecteur d'état augmenté x,,,e prédit pour l'instant k, pour obtenir un vecteur d'état opt augmenté xa,k optimal, compte tenu du vecteur 11Â des entrées mesurées et du vecteur zk des sorties mesurées. Le produit de la matrice .13', de covariance croisée par la transposée du gain de Kalman Kk, est retranché de ia matrice Pa k de covariance de la prédiction d'état, 10 pour obtenir une matrice optimale le de covariance de la prédiction d'état. Une étape 140 a pour but de vérifier en permanence la stabilité des paramètres du p modèle cybernétique compris dans la partie augmentée du vecteur d'état augmenté .yot, optimal, en d'autres termes des composantes du vecteur d'état augmenté optimal qui ne sont pas des états temporels du modèle cybernétique. Dans une première phase après 15 démarrage du véhicule, une variation des valeurs du vecteur n de paramètres traduit une adaptation des paramètres du modèle cybernétique au comportement du conducteur, de nature avantageusement automatique. Une stabilisation des valeurs des paramètres est ensuite représentative du comportement du conducteur ou en d'autres termes de son état de vigilance, de fatigue ou de distraction qui influent sur sa façon de réagir aux 20 événements extérieur. Pour un comportement, c'est-à-dire un état donné du conducteur, les paramètres du modèle sont fixés et ce sont uniquement les variations de l'environnement perçu qui font varier les actions en sorties. Dans une deuxième phase, après stabilisation des valeurs du vecteur n de paramètres, une variation des valeurs du vecteur n de paramètres traduit une adaptation des paramètres du modèle cybernétique à 25 un comportement différent du conducteur, en d'autres termes à une modification de comportement, passage de son état de vigilance normal à un état d'énervement, de fatigue ou de distraction et réciproquement. On peut prévoir de lever un drapeau numérique associé à une première détection de stabilité selon un critère Ck dans l'étape 3032919 19 140, de sorte que seules les détections suivantes d'instabilité déclarent un changement de comportement du conducteur après avoir levé le drapeau numérique. Une variante possible de réalisation de l'étape 140, peut par exemple se détailler en sept étapes 141, 142, 143, 144, 145, 146, 147. L'étape 141 consiste essentiellement à affecter la valeur du vecteur de paramètre flk, contenue précédemment en mémoire pour l'instant courant k, au vecteur de paramètre nk_l, pour l'instant précédent k-1. L'étape 142 consiste essentiellement à générer une nouvelle valeur du vecteur de -opt paramètre flk, à partir du vecteur d'état augmenté optimal. L'étape 142 illustrée sur la , -op/ 10 figure 8 présente un vecteur de paramètre flk extrait du vecteur d'état augmenté xa,fr optimal au moyen d'une matrice à coefficients majoritairement nul sauf à valoir 1 au croisement de la colonne correspondant à l'ordre d'une composante du vecteur de paramètre flk et de la ligne correspondant à l'ordre de la composante du vecteur d'état -op! augmenté xn,i, optimal qui reproduit la valeur de paramètre de ladite composante du 15 vecteur de paramètre flk. On comprendra que d'autres méthodes d'extraction sont possibles comme par exemple celle consistant à utiliser un double pointeur se déplaçant d'une part sur les composantes du vecteur de paramètre flk et d'autre part sur les composantes du vecteur d'état augmenté x',z optimal qui reproduisent une valeur de paramètre et non pas une valeur d'état temporel du modèle cybernétique. 20 L'étape 143 consiste essentiellement à calculer une variation courante ànk entre la valeur du vecteur de paramètre flk à l'instant courant et la valeur du vecteur de paramètre rik_i à l'instant précédent. De préférence, la variation courante s'applique à une composante à la fois du vecteur de paramètres. En d'autres termes on observe pour chacun d'un ou plusieurs paramètres, individuellement sa variation. 25 L'étape 144 consiste essentiellement à calculer une moyenne arithmétique AII4 de la variation courante Atlk et de N-1 variations Aflk_i précédentes. On effectue ainsi une moyenne sur une fenêtre glissante à N éléments. 3032919 20 L'étape 145 consiste essentiellement à calculer un écart type Cui(k) (standard deviation en anglais) de la variation courante LIflk et des N-1 variations ànkj précédentes. La dispersion donnée par l'écart type CAn(k), est représentative d'un niveau de fluctuation des variations du vecteur de paramètre fl, d'autant plus élevé que le comportement du 5 conducteur est instable. L'étape 147 consiste alors à déclarer stable un degré de stabilité Ck lorsque l'écart type Con(k) est inférieur à un seuil prédéterminé. La valeur du seuil est prédéterminée par exemple expérimentalement sur un échantillon de conducteurs statistiquement représentatif, en mettant ces conducteurs dans divers états de distraction possibles. Le 10 critère de stabilité Ck est déclaré instable lorsque l'écart type Con(k) est supérieur ou égal au seuil prédéterminé. Lorsqu'on observe un seul paramètre, un seuil suffit. Lorsqu'on observe plusieurs paramètres, en d'autres termes plusieurs composantes du vecteur de paramètre, il existe un seuil pour chaque paramètre, identique ou non à un autre seuil. La comparaison est faite pour chaque paramètre observé. 15 Préalablement à l'étape 147 qui termine l'étape 140, une étape 146 est utile pour décaler d'un rang vers les instants antérieurs, les variations Ankj du vecteur de paramètre. Lorsque le degré de stabilité Ci, est celui d'une instabilité, une étape 150 a pour rôle de signaler un changement de comportement du conducteur, par exemple pour émettre une alarme en vue d'alerter le conducteur, voire pour déclencher des actions de 20 sauvegarde en cas de changement de comportement source de danger telles qu'un ralentissement du véhicule ou allumage des feux de détresse. A la suite de l'étape 150 ou directement à la suite de l'étape 140 si le degré de stabilité Ck est déclaré stable, une étape 160 utilise à nouveau la méthode de calcul, connue sous le nom de transformée «Unscented». , 25 Le vecteur d'état augmente optimal est vu dans l'étape 160 comme une variable -op! aléatoire Xe de dimension L, de moyenne :vo, et de covariance e, qui subit une transformation non linéaire de la forme X"'l fa( , Uk). La transformée unscented a,k+1 permet maintenant de prédire une valeur moyenne du vecteur d'état augmenté à 3032919 21 l'instant ultérieur kit, ainsi qu'une valeur Pa , de la matrice de covariance de la prédiction d'état. On peut détailler l'étape 160 en deux étapes 161, 162. L'étape 161 consiste essentiellement à déterminer une racine I"' de P°'; c'est-à- 5 dire une matrice 1° t vérifiant la relation Z°Pt(ïopt)T e par exemple en utilisant une factorisation de Cholesky. On note E;"'` la ième colonne de la matrice°Ptqui a L lignes et L colonnes de composantes. L'étape 162 consiste essentiellement à générer un ensemble de 21+1 points appelés sigma-points, dans un espace de dimension L. 10 Un sigma point central X, a pour coordonnées les composantes du vecteur d'état augmentéi'P.;. A chaque colonne rie de la racine de P,,'Pk' correspond un couple de sigma points [ dii7k" X,c,P;j*L ] centrés sur le sigma point central Les coordonnées de sigma points sont calculées au moyen des formules : 15 L'ensemble optimal de sigma points ainsi généré en étape 160 est utilisé en étape 170 pour anticiper un vecteur d'état augmenté A-',+1 à l'instant suivant k+1 de la manière expliquée à présent. On peut détailler l'étape 170 en cinq étapes 171, 172, 173, 174, 175. 20 Le résultat de l'étape 121 qui consiste essentiellement à lire le vecteur Uk des entrées du modèle telles qu'elles sont effectivement mesurées à l'instant courant k, est réutilisé ici. Dans l'étape 172, chaque sigma point xP," généré dans l'étape 160 est assimilé à un vecteur d'état augmenté possible permettant de prédire une valeur de vecteur d'état augmenté futur du modèle correspondant aux valeurs mesurées du vecteur Uk 25 d'entrée du modèle au moyen de la formule : X7i't = Uk Dans l'exemple illustré par le modèle cybernétique de la figure 5 3032919 22 3,k op n,k p,k K c,k "" oint 1,k ar're' e P >k r,k pt 1,k T n ,k (X ) où -11:;1,1,1 o ' opi api v api x,- opl appartiennent aux coordonnées ,k 3,k p,k r,k -"- i,k du sigma point Xfkl.' , où @far, nearr f5, 15d, appartiennent au vecteur Uk des entrées du modèle et où T désigne la période d'échantillonnage. 5 Le nombre L de composantes du vecteur anticipé, est ici égal à 10 parce que dans l'exemple illustré par le modèle cybernétique de la figure 5, on observe trois états temporels et sept paramètres du modèle. Le nombre L de composantes peut tout à fait être différent de 10 sans sortir du cadre de l'invention. Le nombre L est inférieur à 3 si on n'observe qu'un seul état temporel et un seul paramètre du modèle. Le nombre L est 10 supérieur à 3, voire à 10 si on observe plus d'états internes et plus de paramètres du modèle. L'étape 173 consiste à générer une valeur anticipée ou prédite xn,k+1 de vecteur d'état augmenté, égale à une somme des valeurs de sigma points Xr, pondérée par les coefficients iv/' de pondération mentionnés lors des explications de l'étape 100 15 Xa,k a,k+I 1 De la sorte, le premier ensemble optimal X-eji2 2L+1 sigma points, obtenu par a k expansion autour du singleton {x"P.;-° qui donne une valeur moyenne correspondant au vecteur d'état augmenté x°'P,k optimal, induit un deuxième ensemble optimal v - 12 L :i,k+1 3032919 23 2L+1 sigma points, chacun correspondant à un vecteur d'état augmenté futur possible du modèle induit par un point du premier ensemble optimal de sigma points. La compression du deuxième ensemble optimal frl.k+, ro de sigma points donne alors la valeur moyenne xa,k+1, correspondant à une valeur la plus probable du vecteur d'état augmenté à l'instant -opl 5 suivant k-F1 qu'il convient d'estimer à partir du vecteur d'état augmenté optimal_x',k et du vecteur[[ des entrées mesurées à l'instant courant k. Les valeurs possibles d'états X",k, du modèle et la valeur estimée x,,x+1 de vecteur de d'état augmenté à l'instant suivant k+1 permettent de calculer une matrice estimée de covariance sur les valeurs possibles d'états du modèle dans une étape 174. 10 En relation avec l'étape 174, pour chaque point du premier ensemble optimal de sigma points, on génère un vecteur d'écart d'état [x`',,, relatif au vecteur d'état augmenté xn,k+, prédit à L composantes. Dans l'étape 174, chaque produit du vecteur d'écart de sorties par la transposée du vecteur d'écart de sortie correspondant donne une matrice à L lignes et L colonnes. La 15 somme pondérée par les coefficients W,P des 2L+1 matrices ainsi obtenues donne alors la matrice ï-"a,t+1 de covariance sur les états anticipés ou prédits pour l'instant suivant k+1. L'étape 175 consiste alors à générer une valeur de matrice de covariance P',i,1 pour réitérer à l'instant k+1 suivant les étapes ci-dessus exposées. La matrice de covariance est obtenue en multipliant la matrice 76,k+, par un facteur d'oubli M/M-1 avec M>1, ou ce 20 qui revient au même par M+1/M avec M>0. Notons à ce sujet que la matrice de covariance /),,41 est ici une matrice symétrique à L lignes et à L colonnes. L'étape 170 reboucle ensuite sur l'étape 110 pour une nouvelle itération du procédé à l'instant suivant. L'invention ne se limite pas au modèle particulier exposé ci-dessus à titre purement 25 illustratif et qui ne constitue pas en lui-même un élément de l'invention. 3032919 24 L'exposé du procédé, du dispositif et du véhicule ci-dessus exposé a le mérite de montrer l'avantage de l'invention par rapport aux procédés, dispositifs et véhicules antérieurement connus. Par rapport aux systèmes industriels antérieurs basés sur une mesure indirecte 5 comme celle de l'angle volant, le diagnostic repose ici davantage sur les lois physiques liées au comportement du conducteur que sur des aprioris physiques qui ne sont pas vrais pour l'ensemble des conducteurs ou pour toutes les situations d'utilisation. Par rapport aux approches qui reposent sur une mesure directe, la réalisation de l'invention ne nécessite pas de rajout de capteurs supplémentaires et spécifiques. Les 10 capteurs qui sont utilisés par d'autres systèmes d'aide à la conduite (ADAS pour Advanced Driver Assistance System en anglais) peuvent être utilisés. REVENDICATIONS1. L Procédé de détection d'un changement de comportement de conducteur d'un véhicule automobile au moyen d'un modèle cybernétique défini par un premier vecteur (U) de mesure en relation avec au moins une entrée connectée à un instrument de mesure pour quantifier une donnée sensorimotrice en relation avec le conducteur, par un vecteur (Y) de sortie comportant au moins une sortie pour estimer une donnée actionnelle produite par le conducteur, par un vecteur (ri) de paramétrage comportant au moins un paramètre pour quantifier une donnée physiologique du conducteur, et par un vecteur (X) d'état temporel comportant au moins un état temporel dudit modèle, caractérisé en ce qu'il comprend des étapes exécutées en temps réel pour - calculer (130) au moins une valeur courante (ZP,; ) optimale d'un vecteur (X,) d'état augmenté à un instant présent (k) en fonction d'au moins une valeur précédemment estimée ( ) dudit vecteur (Xa) d'état augmenté et d'un deuxième vecteur (zk) de mesure en relation avec ledit vecteur (Y) de sortie, ledit vecteur (Xa) d'état augmenté comportant au moins une composante reproduisant une composante du vecteur (X) d'état temporel et au moins une composante reproduisant une composante du vecteur (n) de paramétrage ; - évaluer (140) un degré de stabilité C(k) d'au moins une composante de valeur courante ( xa,k du vecteur (X,) d'état augmenté qui reproduit ladite au moins une composante du vecteur (Il) de paramétrage ; - signaler (150) un changement de comportement de conducteur lorsque ledit degré de stabilité C(k) est significatif d'une instabilité. 2. Procédé selon la revendication 1, caractérisé en ce qu'il comprend des étapes exécutées en temps réel pour calculer (120) un gain (Kk) de filtre bayésien à partir de la valeur précédemment estimée ( x',k ) dudit vecteur (Xa) d'état augmenté et dudit premier vecteur (U) de mesure, le gain (Kk) de filtre bayésien constituant un moyen technique de calcul de valeur courante (e;LI"),e,ix ) du vecteur (Xa). 3032919 26 3. Procédé selon la revendication 2, caractérisé en ce qu'il comprend des étapes exécutées en temps réel pour générer (110) un premier ensemble prédictif de sigma points (X, x) associés au vecteur (X0) d'état augmenté, centrés sur la valeur précédemment estimée ( x3,0 dudit 5 vecteur (Xa) d'état augmenté au moyen d'une matrice prédictive (130,k) de covariance d'états statistiquement possibles; calculer (122) un deuxième ensemble prédictif de sigma points ( Yki ) associés au vecteur (Y) de sortie à l'instant présent (k) en fonction du premier ensemble prédictif de sigma points ( rio ) et dudit vecteur (Uk) de mesure; 10 calculer (126) le gain (Kk) de filtre bayésien à partir du deuxième ensemble prédictif de sigma points (1'. ). 4. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'il comprend des étapes exécutées en temps réel pour : 15 calculer (170) une valeur future estimée (x.,,,(-4.1) dudit vecteur (Xa) d'état augmenté à un instant suivant (k+1) à partir de la valeur courante (x',k) optimale dudit vecteur (Xa) d'état augmenté et dudit premier vecteur (U) de mesure. 5. Procédé selon la revendication précédente, caractérisé en ce qu'il comprend des 20 étapes exécutées en temps réel pour : générer (160) un premier ensemble optimal de sigma points ( X:Pk''') associés au l vecteur (X,) d'état augmenté, centrés sur la valeur courante x'Op,k optimale dudit vecteur (X3) d'état augmenté au moyen d'une matrice optimale ( /11,°:`;' ) de covariance d'états statistiquement possibles; 25 calculer (172) un deuxième ensemble optimal de sigma points ( ) associés au vecteur (X3) d'état augmenté à l'instant suivant (k+1) en fonction du premier ensemble optimal de sigma points ( et dudit vecteur (Uk) de mesure; 3032919 27 calculer (173) ia valeur future estimée x,,,k+1 ) dudit vecteur (X,) à partir du deuxième ensemble optimal de sigma pointsX 1 - k+1 - ' 6. Procédé selon l'une des revendications précédentes caractérisé en ce que 5 l'étape pour évaluer (140) ledit degré de stabilité C(k) comprend des étapes consistant à obtenir (142) une valeur (I%) du vecteur de paramétrage à partir du vecteur (Xa) d'état augmenté à l'instant présent ; calculer (144) une évolution de paramètres (Ank) à partir de la valeur (1-1k) à l'instant présent et d'une valeur (1-1k-1) à l'instant précédent du vecteur de paramétrage ; 10 mémoriser (146) dans un registre à décalage ladite évolution de paramètres (Ank) de manière à obtenir une quantité (N) d'évolutions antérieures de paramètres ; calculer (145) un écart type (Cui) sur ladite quantité (N) d'évolutions ; signifier (147) une instabilité dudit degré de stabilité C(k) lorsque ledit écart type (Ca') n'est pas inférieur à un seuil prédéterminé. 15 7. Programme d'ordinateur comprenant des instructions de code de programme pour l'exécution des étapes du procédé selon l'une des revendications 1 à 6 lorsque ledit programme est exécuté sur un ordinateur. 20 8. Dispositif de détection d'un changement de comportement de conducteur d'un véhicule automobile comprenant une unité de post traitement (10) connectée à un ou plusieurs instruments de mesure (6, 8) pour quantifier une ou plusieurs données sensorimotrices en relation avec le conducteur, une unité de calcul (11) connectée à l'unité de post traitement (10) et hébergée dans un ordinateur embarqué qui comporte en 25 mémoire un modèle cybernétique défini par un vecteur (U) de mesure comportant au moins une entrée, par un vecteur (Y) de sortie comportant au moins une sortie pour estimer une donnée actionnelle produite par le conducteur, par un vecteur (r1) de paramétrage comportant au moins un paramètre pour quantifier une donnée physiologique du conducteur, et par un vecteur (X) d'état temporel comportant au moins 3032919 28 un état temporel dudit modèle, caractérisé en ce qu'il comprend en mémoire un programme selon la revendication 7 et un processeur pour exécuter en temps réel ledit programme. 5 9. Véhicule comprenant un dispositif selon la revendication 8.