FR 2928219 A1 20090904 FR 0851329 A 20080229 B07-2415FR û AxC/EVH Société par actions simplifiée dite : RENAULT s.a.s. Procédé et dispositif de diagnostic d'un système de commande utilisant un modèle dynamique Invention de : LORIMIER Lionel Procédé et dispositif de diagnostic d'un système de commande utilisant un modèle dynamique La présente invention est relative au diagnostic du fonctionnement d'un système de commande d'au moins un paramètre de conduite d'un véhicule automobile utilisant un modèle dynamique. Dans le but d'améliorer la sécurité active et le plaisir de conduite, certains véhicules automobiles sont équipés de dispositifs d'aide à la conduite tels que des systèmes anti-patinage, des systèmes de freinage automatique, des systèmes de braquage de roues, etc. De tels systèmes sont pilotés par des lois de commande activées par un superviseur selon des fonctionnements répondant à un certain nombre de conditions. Les lois de commande sont implantées dans un calculateur embarqué dans le véhicule et génèrent périodiquement, à une certaine fréquence d'échantillonnage, des signaux de commande que l'on appelle des requêtes, destinés à des actionneurs agissant sur certains organes du véhicule. Dans le cas d'un système actif de braquage de roues arrière dans un véhicule comprenant au moins trois roues directrices, le calculateur émettra des requêtes de braquage du système de direction des roues arrière. Afin de pouvoir procéder au diagnostic de tels systèmes, dans le cadre d'un service après-vente ou en cas d'accident, un certain nombre d'informations sont enregistrées dans une mémoire non volatile qui peut ensuite être exploitée. Dans la majorité des cas, et en raison d'une taille insuffisante de la mémoire disponible, l'enregistrement des données est sous-échantillonné, c'est-à-dire réalisé avec une fréquence d'échantillonnage inférieure à la fréquence d'échantillonnage produisant les requêtes de commande lors du fonctionnement du système. Les données enregistrées étant traitées par des calculateurs, ont une certaine précision dont il est en outre nécessaire de tenir compte lorsque l'on veut effectuer un diagnostic. La demande de brevet japonais JP 2000/181 742 (Fujitsu) décrit un dispositif permettant d'effectuer un diagnostic de faute en ligne dans un système redondant. Une reconstruction de données est effectuée en fonction du résultat du diagnostic. La demande de brevet US 2004/122 639 (Bosch) décrit une méthode d'acquisition de paramètres de conduite d'un véhicule automobile utilisant un modèle tridimensionnel de la cinématique du véhicule, de façon à reconstruire ultérieurement le mouvement du véhicule à partir des signaux de mesure représentatifs des dynamiques latérale et longitudinale. Le modèle utilisé est un modèle cinématique. La présente invention a pour objet de permettre un diagnostic de commande de paramètres de conduite ou de requêtes produites par l'intermédiaire d'une loi de commande représentée par un modèle dynamique, c'est-à-dire un modèle dans lequel les signaux de sortie sont définis par des équations algébriques différentielles en fonction des entrées. La présente invention a également pour objet de proposer des moyens pour l'initialisation d'un modèle dynamique, utilisés pour reconstruire des signaux à partir de données enregistrées en vue d'effectuer un diagnostic. Selon un aspect général, il est proposé un procédé de diagnostic du fonctionnement d'un système de commande d'un paramètre de conduite d'un véhicule automobile, utilisant un modèle dynamique, le diagnostic étant fait à partir de données d'entrée et de sortie du système qui ont été enregistrées pendant le fonctionnement selon une certaine fréquence d'échantillonnage. Le procédé comprend les étapes suivantes : enregistrement de données d'entrée et de sortie du système avec une fréquence d'échantillonnage inférieure à la fréquence d'échantillonnage du système ; stimulation du modèle dynamique avec les données d'entrée enregistrées de façon à déterminer des données de sortie reconstituées; comparaison des données de sortie reconstituées avec les données de sortie enregistrées en vue d'un diagnostic de cohérence. Compte tenu de la fréquence d'échantillonnage utilisée lors de l'enregistrement, le procédé comprend de préférence une étape préalable d'interpolation des données enregistrées au pas d'échantillonnage du système. Avantageusement, on procède ensuite à une étape de reconstitution des paramètres et coefficients de correction statique à partir de données d'entrée enregistrées. L'étape de comparaison est faite par exemple en comparant l'écart entre les données reconstituées et les données enregistrées avec une valeur de seuil pour chaque donnée. On en déduit une information d'alerte si ledit écart est supérieur à la valeur de seuil. Pour pouvoir reconstituer correctement les données de sortie à l'aide du modèle dynamique, il est important de connaître avec précision l'état initial du système au moment où l'enregistrement des données est activé. A cet effet, le procédé comprend de préférence, avant l'étape de stimulation du modèle dynamique, une étape de reconstruction du vecteur d'état initial à partir de données d'entrée et de sortie enregistrées. Généralement, le modèle dynamique utilise, pour chaque pas d'échantillonnage, des équations dynamiques discrétisées faisant intervenir des variables d'état du modèle. L'étape de reconstruction du vecteur d'état initial est alors effectuée par inversion d'un système d'équations comprenant des données initiales enregistrées et des équations dynamiques précitées correspondant à un nombre minimal de pas d'échantillonnage, à partir de l'état initial. Selon un autre aspect, il est également proposé un système de diagnostic du fonctionnement d'un système de commande d'un paramètre de conduite d'un véhicule automobile, utilisant un modèle dynamique, comprenant des moyens d'enregistrement sur une mémoire non volatile, de données d'entrée et de sortie du système pendant le fonctionnement. Les moyens d'enregistrement sont conçus pour enregistrer lesdites données avec une fréquence d'échantillonnage inférieure à la fréquence d'échantillonnage du système. Le système comprend un modèle dynamique capable d'être stimulé avec les données d'entrée enregistrées de façon à déterminer des données de sortie reconstituées. Des moyens de comparaison sont également prévus pour comparer des données de sortie reconstituées avec les données de sortie enregistrées en vue d'un diagnostic de cohérence. De préférence, le modèle dynamique comprend des équations dynamiques discrétisées faisant intervenir, pour chaque pas d'échantillonnage, des variables d'état du modèle. Le système comprend des moyens de reconstruction du vecteur d'état initial par inversion d'un système d'équations comprenant des données initiales enregistrées et les équations dynamiques précitées correspondant à un nombre minimal de pas d'échantillonnage à partir de l'état initial. On va maintenant expliquer de manière plus précise, la manière dont on peut reconstruire un vecteur d'état initial à partir de données d'entrée et de sortie enregistrées avec un nombre minimal de pas d'échantillonnage. Le vecteur d'entrée peut être défini sous la forme Ul [k] U[k]= (1) U; [k] Chacune des composantes de ce vecteur d'entrée de 1 à j correspond à un instant d'échantillonnage noté k pendant la période d'échantillonnage Te pendant laquelle l'enregistrement est activé. Les m données de sortie sur lesquelles doit porter le diagnostic peuvent être exprimées par un vecteur de sortie Y sous la forme : [k] [k] _ (2) ~m [k] Enfin, l'état relatif au système de dimension n, qui représente les valeurs de sortie sur la base de données d'entrée rentrant dans le modèle dynamique, s'exprime sous la forme :25 [k] xn [k] Le modèle dynamique utilise, pour chaque pas d'échantillonnage, des équations dynamiques discrétisées, sous la forme : x[k+l]=Ak .X[k]+Bk(U[k]) (4) X [O] où k est positif ou nul, et où Ak et Bk sont des paramètres exprimés sous forme matricielle. I1 résulte de la forme des équations (4) ci-dessus, que l'évolution des valeurs de sortie issues du modèle dynamique est 10 linéaire par rapport à lui-même comme le montre le premier terme du produit matriciel Ak.X[k]. En revanche, l'évolution peut être non linéaire par rapport aux données d'entrée, comme l'exprime la deuxième relation sous forme du produit matriciel Bk(U[k]). Le vecteur de sortie Y dépend linéairement de l'état X et 15 éventuellement, de manière non linéaire, des données d'entrée U selon la relation : Y[k]=Ck .X[k]+Dk(U[k]) (5) X[k] = (3) où Ck et Dk sont des paramètres exprimés sous forme 20 matricielle. A l'instant initial qui correspond à k=0, les données d'entrée U[O] et les données de sortie Y[O] sont connues car elles ont fait l'objet d'un enregistrement dans la mémoire non volatile du système. L'équation (4) comporte n inconnues (X(0)) et m équations 25 sous la forme : Y[o]= co . x[o]+Do(U[o]) (6) En supposant que toutes les relations sont bien indépendantes, on peut donc construire un système de m équations. Si m est supérieur ou égal à n, et que la matrice Co est inversible, le système d'équations 30 obtenu permet de déterminer X[0]. Dans le cas contraire, il est 5 nécessaire d'utiliser les relations existant à l'instant d'échantillonnage suivant pour lequel k=l, pour obtenir plus d'équations. On fait alors intervenir les entrées U[1] et les sorties Y[1] et on utilise les relations représentées par les équations (4) et (5) ci-dessus au pas d'échantillonnage k=l. I1 en résulte que l'on ajoute les d'inconnues, et les équations X [1] = Ao • X [0] + Bo (U[0]) et Y[1] = C~ • X[1]+ DI(U[1]) au système d'équations. D'une manière générale, on a donc 2.n équations, à condition bien entendu que les indépendantes. Si le système d'équations est alors le vecteur X[0] par inversion matricielle. Dans le cas contraire, il est nécessaire de réitérer le processus une nouvelle fois. On considère alors l'instant d'échantillonnage suivant pour k=2 à l'instant 2.Te. Cela conduit à un vecteur de 3.n inconnues sous la forme : inconnues X(1) au vecteur inconnues pour 2.m+n relations soient bien inversible, on obtient (9) avec en outre 3.m+2.n équations. Si l'on poursuit encore l'itération jusqu'à l'instant k=p à l'instant p.Te7 on obtient un vecteur de p.n inconnues sous la forme _x[0]_ (10) avec (p+1).m+p.n équations. 5 On poursuit les itérations jusqu'à ce que l'on obtienne plus d'équations que d'inconnues. En ne conservant que le nombre d'équations nécessaires à l'inversion du système d'équations, on obtient alors l'état initial X[0] par l'équation matricielle : Co 0 0 0 -x[o]- Y[0]ùDo(U[0]) -Ao I 0 0 Bo (U[0]) 0 C1 0 x[l] (11) 0 0 ùA_ 1 I X[p] Bp-1(U[p -1]) 0 0 0 Cp Y[p] ù Dp (U[p]) où I est la matrice unité. La valeur de p correspond à la n ratio --1.valeur supérieure entière du 10 Selon un exemple de mise en oeuvre avantageux, le paramètre de conduite qui fait l'objet du diagnostic peut être une requête de braquage d'une roue arrière d'un véhicule comportant au moins trois roues directrices. Les données initiales enregistrées utilisées dans le système d'équations précité peuvent alors comprendre la vitesse longitudinale du 15 véhicule, l'angle de braquage des roues avant, la partie dynamique de l'angle de braquage de roue arrière, la partie statique de l'angle de braquage de roue arrière et la valeur de consigne de l'angle de braquage de roue arrière. Les équations dynamiques précitées comprennent comme inconnues, la valeur modélisée de l'angle de braquage de roue arrière, la 20 vitesse de lacet, la dérive latérale et une valeur intermédiaire de réaction positive de l'angle de braquage de roue arrière. Si l'on ajoute la valeur de consigne de l'angle de braquage de roue arrière, on est donc en présence de quatre états. La valeur de consigne de l'angle de braquage de roue arrière est toutefois entièrement déterminée par la 25 connaissance de l'entrée et de la sortie à l'instant k. I1 suffit dans ce cas de prendre en compte les équations précédentes pour quatre pas d'échantillonnage enregistrés, c'est à dire jusqu'à l'instant d'échantillonnage 3.Te (de 0 à p=3). 9 Dans une première application, le vecteur d'état initial permettant l'initialisation du modèle dynamique n'a pas été enregistré. Les données reconstituées utilisées dans l'étape de comparaison sont alors des données reconstituées à partir d'un état initial reconstruit. Dans une deuxième application, le vecteur d'état initial permettant l'initialisation du modèle dynamique a au contraire été enregistré. Le procédé comprend alors une étape supplémentaire de vérification préalable de cohérence entre le vecteur d'état initial enregistré et le vecteur d'état initial reconstruit par comparaison avec des valeurs de seuil, des écarts entre les composants du vecteur d'état initial enregistré et les composants du vecteur d'état initial reconstruit. Si la vérification préalable de cohérence montre une cohérence, l'étape de stimulation du modèle dynamique avec les données d'entrée enregistrées, en vue de déterminer des données de sortie reconstituées, est effectuée à partir du vecteur d'état initial enregistré. Si la vérification préalable de cohérence montre une incohérence, on procède, à partir du vecteur d'état initial enregistré, à une première stimulation du modèle dynamique avec les données d'entrée enregistrées pour déterminer des premières données de sortie reconstituées puis, à partir du vecteur d'état initial reconstruit, à une deuxième stimulation du modèle dynamique avec les données d'entrée enregistrées pour déterminer des deuxièmes données de sortie reconstituées et on compare les premières données de sortie reconstituées, les deuxièmes données de sortie reconstituées et les données de sortie enregistrées en vue du diagnostic de cohérence final. L'invention sera mieux comprise à l'étude de quelques modes de réalisation et de mise en oeuvre pris à titre d'exemples nullement limitatifs, et illustrés par les dessins annexés sur lesquels : la figure 1 illustre schématiquement les principaux éléments inclus dans un calculateur comportant un modèle dynamique pour la commande de braquage d'une roue arrière de véhicule automobile selon une première variante ; la figure 2 illustre les principaux éléments d'un système de diagnostic permettant de vérifier le fonctionnement du système de commande comportant le modèle dynamique illustré sur la figure 1; la figure 3 illustre les différentes étapes d'un procédé de diagnostic selon l'invention mis en oeuvre avec un système tel qu'illustré sur la figure 2; la figure 4 illustre une deuxième variante de calculateur embarqué comportant un modèle dynamique pour la détermination de requêtes de commande pour le braquage d'une roue arrière de véhicule automobile, avec cette fois l'enregistrement d'un nombre plus important de données ; et - la figure 5 illustre les différentes étapes d'un procédé de diagnostic mis en oeuvre à l'aide d'un système tel qu'illustré sur la figure 2, associé à un calculateur tel qu'illustré sur la figure 4. Les différents exemples illustrés à titre non limitatif, s'appliquent à un système de commande du braquage de roues arrière directrices d'un véhicule automobile, tel que décrit en particulier dans la demande de brevet français n° 2 864 002 (Renault) qui utilise une 20 loi de commande par placement de pôles pour déterminer une valeur de consigne d'un angle de braquage des roues arrière. Un tel système de commande permet de générer, au moyen d'un modèle dynamique, des valeurs de requêtes d'angle de braquage d'au moins une roue arrière, ces requêtes étant fournies à un dispositif 25 actionneur capable d'effectuer le braquage requis des roues arrière. Le système comprend un modèle dynamique permettant notamment de modéliser la dynamique latérale du véhicule par l'évolution d'un certain nombre de grandeurs d'étapes qui caractérisent le mouvement du véhicule dans l'espace. Le système comprend en outre un module de 30 réaction positive capable d'élaborer une valeur de consigne d'angle de braquage de roue arrière à partir d'une commande et permettant d'agir sur la dynamique de réponse transitoire. Le module élabore également une valeur de commande statique. 10 15 11 Le procédé de mise en oeuvre d'un tel système tel que décrit dans cette demande de brevet comporte en outre, l'activation ou la désactivation sélective des différents modules du système de façon à tenir compte des différentes situations auxquelles le véhicule est confronté afin d'obtenir, dans certaines situations, une valeur de consigne d'angle de braquage de roues arrière améliorant le comportement du véhicule et le confort de conduite. Le système de diagnostic selon l'invention comprend des moyens d'enregistrement sur une mémoire non volatile embarquée dans le véhicule, d'un certain nombre de données d'entrée et de sortie du système de commande. Compte tenu de la taille limitée de la mémoire prévue dans un calculateur embarqué sur un véhicule automobile, l'enregistrement de ces données ne se fait de préférence qu'à certains moments particuliers pour lesquels l'enregistrement des données semble important. Ce sera le cas, par exemple, lors du déclenchement d'un système anti-dérapage ou d'un système de braquage de roues arrière, ces systèmes entrant en fonction dans des situations particulières de conduite du véhicule. De plus et toujours pour tenir compte de la taille limitée de la mémoire disponible, les données enregistrées ne seront enregistrées qu'avec une fréquence d'échantillonnage inférieure à celle du système de commande. Un signal d'activation de l'enregistrement passe à l'instant t=0 le signal passe de la valeur 0 à la valeur 1. Ce front de montée entraîne l'activation de l'enregistrement. Te est la valeur de chaque pas d'échantillonnage du système de commande de braquage. L'enregistrement se fait dans un rapport 5 par rapport à l'échantillonnage du système de commande. Le pas d'enregistrement étant Tr, on a Tr = STe. A chaque pas d'enregistrement, on enregistre les valeurs des données d'entrée que constitue l'angle de braquage des roues avant aav (en radians). Cet angle est mesuré ou estimé par exemple à partir de la mesure de l'angle de rotation du volant de direction du véhicule. On enregistre également la vitesse longitudinale du véhicule vX en m/s. Cette vitesse est mesurée ou estimée par exemple à partir de la connaissance des vitesses de rotation des roues ou encore à partir de la dérivée filtrée de la position délivrée par un système de positionnement géographique du véhicule (GPS, marque déposée). De la même manière, à chaque pas d'échantillonnage on enregistre deux valeurs de sortie à savoir la requête de braquage statique de roues arrière aq;°' (en radians) et la requête de braquage dynamique des roues arrière aâ@" (en radians). Ces quatre données d'entrée et de sortie sont enregistrées avec un sous échantillonnage Tr par rapport à l'échantillonnage Te des calculs du système de commande de braquage. A l'instant t=0 on enregistre en outre, à titre de donnée d'entrée supplémentaire, la distance parcourue Dp au temps t=0, c'est-à-dire au début de l'enregistrement (en km). Ce signal est issu du compteur de kilomètres parcourus depuis le début de la vie du véhicule. On se reportera maintenant à la figure 1 qui montre schématiquement les principaux organes d'un calculateur embarqué sur un véhicule automobile et capable d'assurer une commande de braquage des roues arrière, comme indiqué par exemple dans la demande de brevet français n° 2 864 002. Sur la figure 1, le calculateur, référencé 1 dans son ensemble, comprend un bloc d'entrée 2 qui reçoit à chaque pas d'échantillonnage, à la fréquence d'échantillonnage Te, les valeurs mesurées de l'angle de braquage des roues avant aa,, et de la vitesse longitudinale du véhicule vx. Un bloc 3 de calcul de braquage statique reçoit sur ses deux entrées, les valeurs mesurées de l'angle de braquage des roues avant aav et de la vitesse longitudinale du véhicule vX issues du bloc d'entrée 2. Le bloc 3 délivre le taux de gain statique Tgs qui est un paramètre de réglage fonction de la vitesse du véhicule et de l'angle de braquage des roues avant. Ce paramètre est défini au moment de la mise au point du véhicule. Le bloc 3 délivre également un signal de requête de braquage statique des roues arrière °`. Cette valeur est, par exemple, calculée à partir de l'angle de braquage des roues avant et du taux de gain statique par la formule aâ'.°` =1 (1 - Tgs(a~ä , v.r )) ' aa,, (12) Le calculateur 1 comprend également un bloc de calcul 4 qui comporte deux modèles qui ne sont pas identifiés de manière précise sur la figure et qui sont, l'un un modèle de la dynamique latérale du véhicule et l'autre un modèle de la dynamique de l'actionneur de braquage des roues arrière. Le modèle de la dynamique latérale du véhicule tient compte de l'évolution des grandeurs d'état que sont la vitesse de lacet y'f et la dérive latérale du véhicule 3. On peut numériser les équations différentielles qui décrivent l'évolution de ces variables selon la méthode d'Euler de façon à obtenir un modèle linéaire décrit par les équations aux différences suivantes: 8[k+l]= M Dak] aAä[k]+M, xxa[k] aû [k]+ z[ (1+Dar.IlùDar '2.v[k] M•vx[k] 2 I, aa,[k]ù Te • Dar • I2 aü;.[k]+T • Dar 12 ùDu, I, .8[k] 1._ 2 z +1ùT Du"11' vxCk; IZ / ~V[k] (14) Avec k>0 le k'ème instant d'échantillonnage, Da,, la rigidité de dérive du train avant (N/rad), 20 Dar celle du train arrière (N/rad), I,, l'inertie en rotation du véhicule autour de son axe de lacet (verticale ascendante) (kg.m2), M la masse du véhicule (en kg), l la distance entre le centre de gravité et l'axe du train avant (m), 25 12 la distance entre le centre de gravité et l'axe du train arrière (m), et L= l,+ 12 l'empattement du véhicule. Le modèle de la dynamique de l'actionneur de braquage des roues arrière donne une estimation de l'évolution du braquage des roues arrière en fonction des consignes de braquage. Ce modèle est ùTe • 1ùT .(D+DQ,.)\,~CIÇ] M v[k] (13) yi[k+l]= également décrit par une équation aux différences résultant de la numérisation de l'équation différentielle caractérisant une dynamique du premier ordre de l'actionneur selon la méthode d'Euler : a"' [k +11=(1ùT •a"[k]+T . a' [k] z~ z où • r est la constante de temps caractéristique du modèle dynamique du premier ordre, a",. est la valeur modélisée de l'angle de braquage de roues arrière et, aQ, est la valeur de consigne de l'angle de braquage de roues arrière. 10 On notera qu'à chaque instant on a : d}gin stol a = a"r + a(, Y Le calculateur 1 comprend en outre un bloc 5 permettant le calcul d'une loi de commande par placement de pôles comme décrit par exemple dans la demande de brevet français n° 2 864 002. Ce bloc 15 délivre une variable intermédiaire qui correspond à une réaction positive dans le système, comme décrit dans la demande de brevet français précité. Cette variable intermédiaire est obtenue par l'équation : aâ''`y[k]=ùK,[k]•ÿ[k]ùKz[k]•8[k]ùK3[k]•aÜ;.[k]+K[k] •a",,[k] (17) 20 Où les coefficients du correcteur KI, K2 et K3 sont obtenus comme indiqué dans la demande de brevet précitée. Le coefficient K est donné par l'équation : K[k]=K,[k]+Tgs[k]• ùK2[k]• +0ù Tgs[k]) .0+ K3 [k]) 25 où C_DFF et C_da sont des coefficients fonction des paramètres géométriques et des rigidités de dérive des trains avant et arrière du véhicule. En l'occurrence, � � C_DFF=2 Met C da=M•2 Da'.ùQ•D",' L • D",, (15) (16) vll k +C_DFF• vx[k] +K,[k] v.K[k] [ ] L+Cda•v,.[k] (18) i L•Dav . D Le bloc 6 illustré sur la figure 1 et symbolisé par la mention 1/z, provoque un retard d'un pas d'échantillonnage qui se traduit par la formule suivante : aâr[k+1]=aarr`q[k] (19) On notera que la valeur de consigne de l'angle de braquage de roues arrière apr est initialisée à l'instant du début de l'enregistrement de manière indépendante et sans respecter cette équation. Le bloc 7 est un bloc d'addition qui reçoit sur son entrée positive la valeur de consigne aQ,. issue du bloc 5 et sur son entrée négative, la requête de braquage statique aa'r' issue du bloc 3. Le bloc additionneur 7 délivre donc la requête de braquage dynamique des roues arrière selon la formule : aarn [k] = aar [k] ù aâY' [k] (20) Une mémoire non volatile référencée 8 équipe le calculateur 1 et permet l'enregistrement, comme indiqué précédemment, des données d'entrée et de sortie à la période d'échantillonnage Tr. L'enregistrement de ces informations démarre dès que le signal d'activation passe de la valeur 0 à la valeur 1. L'enregistrement s'arrête automatique lorsque le nombre de données enregistrées requis est atteint. Les données d'entrées aa, et vX sont amenées à la mémoire 8 par les connexions 9 et 10. La mémoire 8 reçoit également au début de l'enregistrement, la distance parcourue Dp par la connexion 11. L'enregistrement dans la mémoire 8 débute lors de la réception du signal d'activation Act par l'intermédiaire de la connexion 12. Les données de sortie constituées par la requête de braquage statique aâ;.°' amenée par la connexion 13 et la requête de braquage dynamique aârn amenée par la connexion 14 sont également enregistrées, comme indiqué précédemment, selon la période d'enregistrement Tr supérieure ou égale à la période d'échantillonnage Te de la stratégie de commande de braquage de roue arrière, afin de limiter le nombre de données enregistrées. Un bloc de commande 15 reçoit la requête de braquage statique aâ;.°' par la connexion 16 et la requête de braquage dynamique aârn par 16 la connexion 17 et agit directement sur les actionneurs pour le braquage de roues arrière. La figure 2 illustre les principaux organes d'un système de diagnostic permettant d'établir un diagnostic de cohérence des données enregistrées par le calculateur 1 lors du fonctionnement du système de braquage de roues arrière illustré sur la figure 1. On retrouve sur la figure 2 le calculateur 1 comprenant la mémoire non volatile 8. Les données enregistrées dans la mémoire 8 peuvent être récupérées dans un simulateur référencé 18 dans son ensemble, par des moyens de transmission 19 de type classique et non décrits ici. La récupération des données enregistrées qui permet de lire le contenu de la mémoire non volatile 8, inclut différents traitements non illustrés ici, nécessaires pour rendre lisibles les données par le simulateur et pouvant inclure, par exemple, des étapes de décodage. Les données enregistrées dans la mémoire 8 se retrouvent donc dans un bloc d'entrée 20 à l'intérieur du simulateur 18. Ces données ont été enregistrées comme indiqué précédemment, avec une fréquence d'échantillonnage Tr. Il convient tout d'abord de procéder à une interpolation des données enregistrées de façon à reconstituer pour toutes les entrées et sorties enregistrées, la valeur des données au pas d'échantillonnage Te. Cette opération est effectuée dans le bloc d'interpolation 21. Avant de procéder aux étapes suivantes, il est nécessaire de reconstituer les valeurs du paramètre Tgs et des coefficients des correcteurs à chaque pas d'échantillonnage, ces coefficients de correcteurs ayant été utilisés pour calculer la requête de braquage des roues arrière. Cette opération est effectuée dans le bloc 22 à partir des données de braquage des roues avant et de la vitesse du véhicule, ces données étant interpolées et fournies par la connexion 23 issue du bloc 21. Ces informations sont ainsi obtenues pour toute la durée de l'enregistrement et à la période d'échantillonnage Te. Si k=0 est l'instant initial à partir duquel l'enregistrement a commencé, les variables d'état du modèle dynamique, à savoir la vitesse de lacet ~jr, la dérive latérale 8 et la valeur modélisée de l'angle de braquage de roues arrière aâ;. ont des valeurs inconnues qui ne sont pas nécessairement nulles. Afin de pouvoir reconstituer les requêtes dynamiques en stimulant une copie du modèle dynamique implanté dans le calculateur 1 avec les entrées qui ont été enregistrées, il est nécessaire de reconstruire l'état initial du modèle. Cette reconstruction se fait dans le bloc 24 d'une manière qui sera explicitée par la suite. Cet état initial reconstitué est amené par la connexion 24a à l'entrée d'un bloc 25 qui est identique au bloc 4 du calculateur 1, et qui comprend un modèle dynamique identique. Le bloc 25 reçoit sur ses entrées la valeur des données interpolées issues du bloc 21 pour l'angle de braquage des roues avant aa,, et la vitesse longitudinale du véhicule vx. Le bloc 25 reçoit également sur son entrée la valeur de consigne de l'angle de braquage de roues arrière aâ, qui est calculée par un bloc 26 correspondant au bloc 5 du calculateur 1 et qui contient la même loi de commande par placement de pôles. Le bloc 26 reçoit sur ses différentes entrées les valeurs déterminées par le modèle dynamique du bloc 25 que constituent la vitesse de lacet tfr, la dérive latérale 8 et la valeur modélisée de la requête de braquage de roue arrière . Le bloc 26 reçoit également par la connexion 48 la valeur interpolée de l'angle de braquage de roue avant aav. Le bloc 26 délivre à sa sortie la variable intermédiaire de réaction positive agi'' pour la requête d'angle de braquage de roue arrière qui fait l'objet, par le bloc 27, d'un retard d'un pas d'échantillonnage de façon à produire la valeur de consigne d'angle de braquage de roue arrière aâ, qui est ramenée en réaction par la connexion 29 à l'entrée du bloc 25. Cette valeur est également amenée sur l'entrée positive du bloc additionneur 28, qui reçoit par ailleurs sur son entrée négative, par la connexion 30, la valeur interpolée de la requête de braquage statique des roues arrière aa'r'' On obtient finalement à la sortie du bloc additionneur 28, comme c'était le cas à la sortie du bloc additionneur 7 du calculateur 1, une valeur de requête de braquage dynamique de roue arrière qui est cette fois obtenue par reconstitution et qui est notée ad'"'-recon,"`é Cette valeur reconstituée, qui a été obtenue par interpolation dans le bloc 21 des données enregistrées dans la mémoire 8 du calculateur 1, est comparée, dans un bloc 31 de comparaison et de diagnostic avec la valeur correspondante enregistrée de la requête de braquage dynamique aar'.nùe"'•eg'`"'é qui est amenée par la connexion 32 issue du bloc 20 jusqu'au bloc 31 en vue d'un diagnostic de cohérence. En se reportant à la figure 3, on voit que la mise en oeuvre du système tel qu'illustré sur la figure 2 se fait par une succession d'étapes. La première étape 33 qui se fait en amont du bloc 20 illustré sur la figure 2, permet la récupération des données enregistrées. Elle consiste à lire le contenu de la mémoire non volatile 8 du calculateur 1, qui contient les données enregistrées sous-échantillonnées. La deuxième étape 34 qui s'effectue dans le bloc 21 permet l'interpolation des données enregistrées au pas d'échantillonnage du système Te. Cette interpolation peut se faire par exemple de manière linéaire. Si l'on note le ratio de sous-échantillonnage de l'enregistrement des données par n=Tr/Te, on obtient alors l'angle de braquage des roues avant à l'instant n.k où k est un entier positif ou nul, par la formule : a interpolé [n . k] ù a e regtviré [n • k] (21) Pour chaque instant m compris entre n.k et n.(k+l), il est nécessaire de reconstituer la donnée enregistrée. On pourra par exemple effectuer une interpolation linéaire en se basant sur les données connues, à savoir aevreg;atré [n.k] et a enregi c"'é [n . (k + 1)] . L'interpolation peut se faire par l'équation : enregistré [n (k + 1)ù a enregisn é [n k] interpolé enregistré *aav av aav [n k+m]=aav [n k] +m (22) n La même opération d'interpolation se fait sur toutes les autres données enregistrées, dans les mêmes conditions.30 19 L'étape suivante consiste à calculer la valeur du paramètre que constitue le taux de gain statique Tgs ainsi que les coefficients des correcteurs pour chaque pas d'échantillonnage. Cette étape est notée 35 sur la figure 3 et elle est mise en oeuvre par le bloc 22 illustré sur la figure 2. Les valeurs reconstituées des coefficients de correcteur et pour Tgs sont obtenues à partir des données de braquage des roues avant et de vitesse du véhicule, interpolées à l'étape précédente. On effectue ensuite à l'étape 36 le calcul de l'état initial du modèle dynamique à partir de la seule connaissance des entrées et des sorties enregistrées ayant fait l'objet de l'interpolation. Etant donné que l'on dispose d'un nombre de valeurs enregistrées restreint en ce qui concerne les sorties du modèle, c'est-à-dire dans l'exemple illustré, les valeurs statique et dynamique de la requête d'angle de braquage de roue arrière, il est important d'utiliser le minimum de points pour reconstituer l'état initial. Si l'on utilise un trop grand nombre de points, le diagnostic final risque d'être faussé. En effet, le diagnostic est basé sur l'interprétation des écarts constatés entre les sorties simulées reconstituées sur la base de l'état initial également reconstruit, et les sorties enregistrées. Le fait d'utiliser trop d'échantillons enregistrés pour les données de sortie provoquerait une diminution de l'écart potentiel en raison du fait que ce ne serait plus l'état initial réel du modèle à l'instant de l'enregistrement sur le véhicule qui serait reconstruit, mais un état fictif différent de cet état réel. On a vu plus haut qu'il était possible de reconstruire l'état initial d'une manière générale à partir d'un nombre minimal d'équations de façon à rendre inversible un système d'équations obtenus à partir des données de sortie calculées sur la base des entrées enregistrées. Dans l'exemple illustré, on procède de la manière suivante. Au début de l'enregistrement des données, on dispose des entrées enregistrées va[0] et aav,[0] ainsi que de la sortie enregistrée aa'.'n[0], qui est la requête dynamique d'angle de braquage de roue arrière. On connaît encore la sortie enregistrée que constitue aü;u`[0], c'est-à-dire la requête de braquage statique de roue arrière. Les variables d'état du modèle dynamique qui constituent les entrées du bloc 26 de la figure 2, sont en revanche inconnues. Il s'agit de la valeur modélisée de l'angle de braquage de roue arrière a";.[0], de la vitesse de lacet 40] et de la dérive latérale 8[0]. I1 en est de même de la variable intermédiaire de réaction positive aâ;''"'[0]. A ce stade, on dispose donc d'une équation de la forme de l'équation (17) ci-dessus pour quatre inconnues. Il est donc impossible de déterminer précisément l'état du modèle aû [0],ôr[0],8[0]. Afin de disposer de davantage d'équations, on utilise les informations relatives à l'instant de calcul suivant. Les données d'entrée vx[l],au,,[l], la donnée de sortie a '."[1] ainsi que la donnée de sortie a;`°`[1] sont connues. On peut en déduire aisément, par une équation du type de l'équation (20), la variable intermédiaire a,,,.[1]. Les équations (13), (14), (15), (17) et (19) apportent cinq nouvelles équations avec quatre nouvelles inconnues, à savoir aâ;.[1],vi[1],8[1],a ,`'i'[l], soit un total cumulé de six équations pour huit inconnues, ce qui reste insuffisant pour déterminer l'état initial car cela constitue un système d'équations indéterminé. On prendra alors les informations à l'instant de calcul T2 qui fournit cinq nouvelles équations et quatre nouvelles inconnues, à savoir, a;;.[2],y [2],8[2],aa�''"'[2], soit un total cumulé de onze équations et douze inconnues. L'utilisation de données de sortie et d'entrée à l'instant T3 ajoute une nouvelle équation en combinant les équations (19) et (20) sans ajouter de nouvelle inconnue puisqu'elle permet de déterminer a 1 aqr 7;Yeq [2] A ce stade le système d'équation est donc inversible, et permet de déterminer toutes les inconnues, et finalement l'état initial du modèle dynamique, à savoir aQr [0], [0], 8 [0]. La combinaison des différentes équations mentionnées ci-dessus donne le système K, [0] K2 [0] K3[0] 0 0 0 0 0 0 0 0 T`ù1 0 0 1 0 0 0 z -al1 -a12 -a13 0 1 0 0 0 0 ùa21 -a22 -a23 1 0 0 0 0 0 0 0 0 K, [1] K2[1] K3[1] 0 0 0 0 0 0 0 0 T`ù1 0 0 1 z- 0 0 0 ùbä ùb12 ùb13 0 1 0 0 0 0 ù b2 ù b22 ù b,3 1 0 0 0 0 0 0 0 0 K,[2] K2[2] K3[2] (23) Dans ce système, les variables intermédiaires, y compris les valeurs de la variable intermédiaire aâ;'"d aux instants de calcul 0, 1 et 2, ont été éliminées. La matrice M contient les données interpolées à chaque pas 10 d'échantillonnage Te et s'écrit [2] 510 K[O]' aav [O]ùaara [1]ùcarat [i] Te (aarn[O]+aarat[O]) i T , Dav a [0] M vx [0] av Te Dav 11 aav[0] K[1].aav[1]ùaarn[2]ùaarat[2] T '(aar[1]+aaat[1]) i T ' D e av av [1] M'vx[1] Te ' Dav ' 11 • av [1] Izz K[2].aav[2]ùaara [3]ùaarat [3] On a de plus pour les coefficients aii : a =ùT ' 1+Dav 'l1 ùDal'. 12 M'vx[O]2 Dar • 12 ù Dav • 11 IZZ T 'D ar ' 12 M'vx[O] Les coefficients b,, sont définis avec les mêmes expressions que ail mais avec vx [1] au lieu de vx [O]. a12 = a13 '1 ù l e ' (Dav + Dar ) M'vx[O] , T. D e ar M'vx[O] ( T Dav h +Dar . e l2 Izz ' vx [O] M= IZZ (24) ( 11 e i Si l'on note Ainit la matrice 9x9 de l'équation (23), M étant un vecteur ligne, c'est-à-dire une matrice 9x1, le vecteur ligne résultant a une dimension 3x1. La matrice Ainit est inversible et l'on peut obtenir l'état initial du modèle dynamique par l'équation : _ w[o] _ 1 0 0 0 0 0 0 0 0 8[0] 0 1 0 0 0 0 0 0 0 . Ai n i t • M (25) aa [0] 0 0 1 0 0 0 0 0 0 Après qu'on ait ainsi reconstitué le vecteur d'état initial, on procède à l'étape 37 indiquée sur la figure 3, qui consiste à reconstituer les requêtes de braquage dans le simulateur 18 de la figure 10 2, les calculs s'effectuant dans les différents blocs 25, 26, 27 et 28. Cette étape consiste à calculer, pour les instants k allant de 0, qui correspond à l'initialisation jusqu'à un temps tenregistrement, l'état du modèle dynamique et les requêtes de braquage produites sur la base des équations (13), (14), (15), (16), (17), (18) et (19). On obtient 15 finalement la valeur de la requête de braquage dynamique des roues arrière reconstituée par l'équation : a dYn _reconstitué [k] ù aa,. [k] ù a star _ interposé [k] (26) 20 Avec tenregistrement ~ k . Te % 0, et où k est le k`éme instant d'échantillonnage. Pour ces calculs, on a utilisé l'état initial du modèle dynamique reconstitué au cours de l'étape 36 et calculé dans le bloc 24, ainsi que les coefficients de correcteur et Tgs calculés pour l'ensemble de 25 l'enregistrement à l'étape 35, le calcul étant effectué dans le bloc 22. La dernière étape 38 indiquée à la figure 3 consiste à effectuer un diagnostic par une vérification de la cohérence des données enregistrées avec les données reconstituées. On pourra de préférence tracer sur un même graphique les valeurs aâyn [j] et aad [n • k] où j et k 30 sont des nombres entiers positifs ou nuls, tels que j et n.k ne dépassent pas le nombre d'échantillons disponibles. La comparaison de ces5 24 valeurs permet de vérifier la cohérence des requêtes de braquage aux différents moments d'enregistrement, ainsi que la tendance générale d'évolution pendant la durée de l'enregistrement. Si l'écart entre les valeurs au'''''s'""` [n k] et a`;'"[n.k] où k est un nombre entier positif ou nul tel que n.k ne dépasse pas le nombre d'échantillons disponibles, dépasse un seuil autorisé qui prend en compte les incertitudes liées à l'interpolation linéaire, à la précision des données etc., un message d'alerte de diagnostic est fourni pour avertir d'une incohérence entre les données enregistrées et les données reconstituées. Une telle incohérence permettra de rechercher la cause d'un dysfonctionnement du système de commande de braquage de roues arrière embarqué sur le véhicule. La figure 4 illustre un deuxième mode de réalisation également appliqué à titre d'exemple au diagnostic d'un système de commande de braquage de roues arrière. Les organes identiques illustrés sur la figure 4 portent les mêmes références que ceux de la figure 1. La seule différence porte sur les entrées de la mémoire non volatile 8. En effet, dans ce mode de réalisation, on enregistre également au moment de l'activation de l'enregistrement (cet instant étant noté ko) la valeur courante des signaux issus du modèle dynamique contenu dans le bloc 4 du calculateur embarqué sur le véhicule. Ces valeurs sont enregistrées par les connexions 39, 40 et 41. Les valeurs ÿr[ko],8[ko],au;[4] sont donc enregistrées dans la mémoire 8. Ces valeurs correspondent à l'état initial du modèle dynamique. Dans ce deuxième mode de réalisation, le procédé se déroule comme illustré sur la figure 5, les quatre premières étapes étant identiques aux quatre premières étapes de la figure 3. En particulier, à l'étape 36, on calcule comme indiqué précédemment, l'état initial du modèle dynamique à partir des seules données d'entrée et de sortie enregistrées et interpolées comme précédemment. Une nouvelle étape 42 permet d'effectuer un diagnostic préliminaire en vérifiant tout d'abord la cohérence entre l'état initial reconstruit et l'état initial enregistré pour le modèle dynamique. Cette comparaison est réalisée sur chaque série de données enregistrées à analyser. Pour analyser cette cohérence, on prend en compte l'incertitude liée à la reconstitution de l'état initial du modèle dynamique à l'étape 36, sur la base de données présentant certaines imprécisions liées au type de mémoire utilisée, à la précision de l'interpolation, etc. On estimera que les données sont cohérentes si les trois conditions suivantes sont réunies : ~L0]û[ko] Si l'on constate une bonne cohérence, c'est-à-dire un écart inférieur au seuil visé ci-dessus, le processus se poursuit par l'étape 43, dans laquelle on reconstitue les requêtes de braquage au moyen d'un simulateur similaire à celui de la figure 2, dans lequel les calculs se font cependant à partir de l'état initial du modèle dynamique tel 20 qu'enregistré dans la mémoire 8. Comme indiqué précédemment, on calcule pour les instants k allant de l'initialisation à tenregistrement, l'état du modèle dynamique ainsi que les valeurs des requêtes de braquage de roues arrière sur la base des équations contenues dans les blocs 25, 26, 27 et 28, à savoir 25 les équations précédentes (13), (14), (15), (16), (17), (18) et (19). On obtient également les valeurs de requête de braquage dynamique de roues arrière reconstituées adyn-ircnnst"ué[k] par l'équation 26. Toutefois, dans ces différents calculs, le modèle dynamique contenu dans le bloc 25 reçoit comme entrée le vecteur d'état initial enregistré 30 ÿr[ko],8[ko],aâr[ko] étant donné que la cohérence entre l'état initial enregistré et l'état initial reconstitué a été constatée au cours de l'étape précédente référencée 42 sur la figure 5. A partir de ces valeurs reconstituées, on procède alors à l'étape de diagnostic référencé 44 sur la figure 5, où l'on compare les valeurs enregistrées avec les valeurs reconstituées. Cette étape se fait de la même manière que l'étape 38 précédemment explicitée pour le premier mode de réalisation. Dans le cas où une incohérence est détectée à l'étape 42, on procède tout d'abord à l'étape 45 qui consiste à reconstituer les requêtes de braquage de roues arrière à partir de l'état initial enregistré. On obtient une requête de braquage de roues arrière notée dyn recon.stitue 1 aa Puis, au cours d'une étape 46, on procède de la même manière à la reconstitution des requêtes de braquage de roues arrière mais cette fois à partir de l'état initial reconstruit. On obtient une autre valeur notée adyn-reconstitué-2 Le processus se poursuit par l'étape 47, dans laquelle on compare les valeurs obtenues au cours des étapes 45 et 46, avec les valeurs enregistrées. On pourra par exemple tracer sur un même graphique les valeurs obtenues aâ;" '[j] a''rec"ns7""é_2[j] et les valeurs enregistrées ap,'."'[n.k}où j et k sont des nombres entiers positifs ou nuls tels que j et n.k ne dépassent pas le nombre d'échantillons disponibles. A partir d'un tel tracé, on vérifie la cohérence des requêtes à chaque instant d'enregistrement, ainsi que la tendance générale d'évolution sur l'intervalle de temps enregistré. Plusieurs cas se présentent alors : Si la différence en valeur absolue entre toutes les données reconstituées et les données enregistrées est inférieure à un seuil déterminé, qui prend en compte les incertitudes liées à l'interpolation linéaire et à la précision des données, on pourra en conclure que les informations enregistrées sont globalement cohérentes avec l'ensemble des informations reconstituées. II est alors impossible de conclure quant au diagnostic, puisqu'une incohérence a été constatée à l'étape 42 en ce qui concerne l'état initial, incohérence que l'on ne retrouve plus lorsque les données ont été reconstruites à partir d'une part, de l'état initial enregistré et d'autre part, de l'état initial reconstruit. Dans un autre cas, il existera pour un enregistrement correspondant à un nombre entier k positif ou nul, simultanément les deux conditions suivantes : a dynreconsttué_l [n . k] ar ùa adyn k] > Qgapautorisé ar et a econstitaé_2 [n k] [n . k] résulte d'un problème d'enregistrement de l'état initial du modèle dynamique ou d'un problème concernant le calcul des requêtes finales. Dans une troisième situation, on pourra constater pour un instant d'enregistrement au moins, noté k (nombre entier positif ou nul), que l'on a simultanément : dyn_reconstitué_1 [n k] ù a dyn [n k] a dyn_ reconstitué _2 [n ki ù ar[n . k] > aap_autorisé Dans ce cas, la reconstitution des requêtes de braquage à partir des données d'entrée et de sortie et de l'état initial enregistré est cohérente avec les données enregistrées. L'incohérence constatée à l'étape 42 résulte donc d'un problème sur les deux premières séries d'échantillons enregistrés qui ont été utilisés pour reconstituer l'état initial. Dans une autre situation, on constatera qu'il existe deux instants correspondant à ki et k2 qui sont deux nombres entiers positifs ou nuls, pour lesquels on a simultanément : econstitué_1 [n kl ] _ k1 I > gap_autorisé ar et a dyn _ recons artué-2 [n ' k 1ù a a dr' [n k 1 > gap autorisé 2 2 ar Dans ce cas, les incohérences constatées concernent un problème de calcul des requêtes finales. On voit ainsi, à l'étude de ces exemples, qu'il est possible, en mettant en oeuvre l'invention, qu'il s'agisse de son premier ou de son deuxième mode de réalisation, d'obtenir un diagnostic sur la cohérence des données enregistrées par rapport à des données reconstituées et d'en déduire un renseignement sur un éventuel dysfonctionnement d'un dispositif de commande d'un ou de plusieurs paramètres de conduite d'un véhicule automobile, par exemple d'une requête de braquage de roues arrière. REVENDICATIONS 1. Procédé de diagnostic du fonctionnement d'un système de commande d'au moins un paramètre de conduite d'un véhicule automobile, utilisant un modèle dynamique, le diagnostic étant fait à partir de données d'entrée et de sortie du système qui ont été enregistrées pendant le fonctionnement, caractérisé par le fait qu'il comprend les étapes suivantes : enregistrement de données d'entrée et de sortie du système avec une fréquence d'échantillonnage inférieure à la fréquence d'échantillonnage du système ; stimulation du modèle dynamique avec les données d'entrée enregistrées de façon à déterminer des données de sortie reconstituées; comparaison des données de sortie reconstituées avec les données de sortie enregistrées en vue d'un diagnostic de cohérence. 2. Procédé selon la revendication 1 comprenant au préalable une étape d'interpolation des données enregistrées au pas d'échantillonnage du système. 3. Procédé selon l'une des revendications précédentes comprenant au préalable une étape de reconstitution des paramètres et coefficients de correction statique à partir de données d'entrée enregistrées. 4. Procédé selon l'une des revendications précédentes dans lequel l'étape de comparaison est faite en comparant l'écart entre les données reconstituées et les données enregistrées avec une valeur de seuil pour chaque donnée et on émet une information d'alerte si ledit écart est supérieur à la valeur de seuil. 5. Procédé selon l'une des revendications précédentes, comprenant, avant l'étape de stimulation du modèle dynamique, une étape de reconstruction du vecteur d'état initial à partir de données d'entrée et de sortie enregistrées. 6. Procédé selon la revendication 5 dans lequel le modèle dynamique utilise, pour chaque pas d'échantillonnage, des équations dynamiques discrétisées faisant intervenir des variables d'état du modèle, l'étape de reconstruction du vecteur d'état initial étanteffectuée par inversion d'un système d'équations comprenant des données initiales enregistrées et des équations dynamiques précitées correspondant à un nombre minimal de pas d'échantillonnage à partir de l'état initial. 7. Procédé selon la revendication 6 dans lequel le paramètre de conduite qui fait l'objet du diagnostic est une requête de braquage d'une roue arrière d'un véhicule comportant au moins trois roues directrices, les données initiales enregistrées utilisées dans le système d'équations précité comprenant la vitesse longitudinale du véhicule, l'angle de braquage des roues avant, la partie dynamique de l'angle de braquage de roue arrière, la partie statique de l'angle de braquage de roue arrière, la valeur de commande de l'angle de braquage de roue arrière, et les équations dynamiques précitées comprennent comme variables, la valeur modélisée de l'angle de braquage de roue arrière, la vitesse de lacet, la dérive latérale et une valeur intermédiaire de réaction positive de l'angle de braquage de roue arrière. 8. Procédé selon l'une des revendications 6 ou 7 dans lequel le vecteur d'état initial permettant l'initialisation du modèle dynamique n'a pas été enregistré, les données reconstituées utilisées dans l'étape de comparaison étant des données reconstituées à partir d'un état initial reconstruit. 9. Procédé selon l'une des revendications 6 ou 7 dans lequel le vecteur d'état initial permettant l'initialisation du modèle dynamique a été enregistré, le procédé comprenant une étape supplémentaire de vérification préalable de cohérence entre le vecteur d'état initial enregistré et le vecteur d'état initial reconstruit par comparaison avec des valeurs de seuil des écarts entre les composants du vecteur d'état initial enregistré et les composants du vecteur d'état initial reconstruit. 10. Procédé selon la revendication 9 dans lequel : si la vérification préalable de cohérence montre une cohérence, l'étape de stimulation du modèle dynamique avec les données d'entrée enregistrées de en vue de déterminer des données de sortie reconstituées est effectuée à partir du vecteur d'état initial enregistré ; si la vérification préalable de cohérence montre une incohérence, on procède,à partir du vecteur d'état initial enregistré, à une première stimulation du modèle dynamique avec les données d'entrée enregistrées pour déterminer des premières données de sortie reconstituées puis, à partir du vecteur d'état initial reconstruit, à une deuxième stimulation du modèle dynamique avec les données d'entrée enregistrées pour déterminer des deuxièmes données de sortie reconstituées et on compare les premières données de sortie reconstituées, les deuxièmes données de sortie reconstituées et les données de sortie enregistrées en vue du diagnostic de cohérence. 11. Système de diagnostic du fonctionnement d'un système de commande d'un paramètre de conduite d'un véhicule automobile, utilisant un modèle dynamique, comprenant des moyens d'enregistrement sur une mémoire non volatile de données d'entrée et de sortie du système pendant le fonctionnement, conçus pour enregistrer lesdites données avec une fréquence d'échantillonnage inférieure à la fréquence d'échantillonnage du système, caractérisé par le fait qu'il comprend un modèle dynamique (25) capable d'être stimulé avec les données d'entrée enregistrées (20) de façon à déterminer des données de sortie reconstituées et des moyens de comparaison (31) pour comparer des données de sortie reconstituées avec les données de sortie enregistrées en vue d'un diagnostic de cohérence. 12. Système selon la revendication 11 dans lequel le modèle dynamique comprend des équations dynamiques discrétisées faisant intervenir, pour chaque pas d'échantillonnage, des variables d'état du modèle, le système comprenant des moyens de reconstruction d'un vecteur d'état initial par inversion d'un système d'équations comprenant des données initiales enregistrées et les équations dynamiques précitées correspondant à un nombre minimal de pas d'échantillonnage à partir de l'état initial.