Procédé d'identification du modèle d'une enceinte thermique et application à la régulation d'une enceinte thermique. La présente invention concerne un procédé d'identification du modèle d'une enceinte à réguler et son application à la régulation d'une enceinte thermique. Il est connu dans les régulations d'utiliser des régulateurs analogiques notamment proportionnel intégral dérivé P I D. qui ont l'inconvénient de devoir être réglés chaque fois que l'on change de point de fonctionnement Ces réglages sont de plus en général longs, peu précis et doivent faire appel à des spécialistes. Il est également connu par la demande de brevet FR 2 372 466 déposée au nom d'ALSTHOM ATLANTIQUE de réaliser une régula- tion utilisant un modèle de représentation du processus à commander, mais ici ce modèle est réalisé analogiquement et ne peut être qu'une approximation de la représentation du processus De plus ce modèle est utilisé avec un régulateur P.I analogique. Un premier but de l'invention est de réaliser un procédé d'identification d'un modèle numérique de représentation d'une enceinte thermique, qui soit valable quel que soit le point de fonctionnement de l'enceinte. Un deuxième but de l'invention est la réalisation d'un régulateur numérique utilisant le modèle de représentation de l'enceinte et admettant des désadaptations du modèle par -2- rapport à la réaction de l'enceinte aux actions appliquées par le régulateur. Un troisième but de l'invention est d'automatiser l'identi- fication et le réglage optimal du régulateur. Suivant l'invention le procédé d'identification du modèle comprend les phases suivantes: détermination approximative du temps de réponse de l'enceinte par mesure et enregistrement de la réponse en température de l'enceinte à un échelon de puissance, élaboration d'un signal d'identification constitué d'une multiplicité de créneaux aléatoires de hauteur et de largeur, comprise entre une limite inférieure et une limite supérieure en rapport avec le temps de réponse du four, variant de façon aléatoire, mesure par échantillonnage, calcul et mise en mémoire de la réponse impulsionnelle de l'enceinte au signal d'identi- fication pour chaque créneau, calcul du modèle de représentation de l'enceinte à partir des réponses impulsionnelles mémorisées. Suivant une deuxième caractéristique de l'invention, le procédé identifie un modèle constitué: par l'enregistrement dans une première série de mémoires d'une première multiplicité de premières valeurs détermi- nées après échantillonnage à des intervalles de temps constants, lesdites valeurs étant représentatives de la réponse impulsionnelle de l'enceinte en fonction du temps pour une température de fonctionnement donnée, lesdites valeurs étant rangées à des adresses correspondant au rang de l'échantillonnage, 12985 -3- par l'enregistrement dans une deuxième série de mémoires d'une deuxième multiplicité de valeurs représentatives de la variation du gain des réponses impulsionnelles en fonction de la température de fonctionnement, lesdites deuxièmes valeurs étant rangées à des adresses correspondant aux dif- férentes températures de fonctionnement, par l'enregistrement dans une troisième série de mémoires d'une troisième multiplicité de valeurs représentatives de la variation de la période d'échantillonnage en fonction des différentes températures de fonctionnement, à des adresses correspondant aux températures. Suivant une troisième caractéristique de l'invention, le procédé est appliqué à la régulation d'une enceinte comprenant de façon connue un élément de mesure 4 de la variable à réguler, un élément de régulation 5 élaborant à partir des valeurs de variables et du point de consigne l'action à imposer à l'enceinte, un organe de commande 11 agissant sur l'enceinte en réponse au signal d'action carac- térisé en ce que l'élément de régulation élabore le signal de commande par un premier ensemble d'algorithme qui tient compte de la mesure de la sortie de l'enceinte, calcule le comportement futur-souhaité en fonction du temps de réponse souhaité du système bouclé, estime grâce au modèle l'évolu- tion future de la sortie et ajuste les commandes futures, applique la première commande, mesure et compare avec les valeurs estimées la mesure pour appliquer la commande suivante en cas d'adaptation ou effectue un nouveau calcul, une nouvelle estimation et un nouvel ajustement avant d'appliquer la nouvelle première commande. L'invention sera mieux comprise à la lecture de la descrip- tion faite en référence aux dessins annexés dans lesquels v La figure 1 représente une régulation constituant l'état de la technique. 4 - La figure 2 est une vue schématique du dispositif de mise en oeuvre de l'invention. La figure 3 est un graphique des courbes représenta- tives des réponses impulsionnelles permettant de constituer le modèle. Les figures 4, 5, 6 représentent les courbes repré- sentatives du modèle de représentation de l'enceinte. La figure 7 est une variante de réalisation de la régulation. La figure 1 représente un procédé de régulation d'une enceinte thermique utilisant de façon connue en soi, un élément de mesure 4 constitué par un thermocouple ou capteur de température, un élément de calcul 5 élaborant à partir des valeurs de la température d'enceinte 9 e et du point de consigne Oc le signal de commande W à appliquer à un organe de commande 11 agissant sur le processus à réguler. L'enceinte thermique, comporte une enceinte calorifugée 1, un actionneur ou une résistance électrique 3 et une charge 2 constituée par les objets que l'on veut chauffer Dans ce type de réalisation classique l'élément de calcul 5 comprenait un comparateur 7 qui élaborait un signal d'écart entre la valeur de consigne et la valeur mesurée et un régu- lateur autonome P I D à action proportionnelle intégrale et dérivée qui élaborait un signal de commande W analogique à partir du signal d'écart. Toutefois ces régulateurs qui fonctionnent correctement pour un point de fonctionnement donné, nécessitent des réglages fastidieux et délicats chaque fois que l'on veut changer de point de fonctionnement. De plus, lorsque l'on applique ce type de régulateur à une 251,2985 - enceinte dont les caractéristiques évoluent dans le temps, non seulement en fonction de son entrée mais aussi en fonction d'autres facteurs tels que par exemple pour un four, les fuites thermiques, la puissance rayonnée par les résistances, les coefficients d'absorption des parois; il est nécessaire d'effectuer des réglages fréquents des para- mètres Ces inconvénients et bien d'autres sont évités par le dispositif de l'invention comme nous le verrons dans la suite de la description. Le dispositif de régulation qui va être décrit à titre d'exemple en se référant à la figure 2 sert à réguler la température d'un four ou enceinte thermique Ce dispositif comprend comme précédemment une enceinte 1, une résistance électrique de chauffage branchée à un dispositif de commande fournissant la puissance, ce dispositif pouvant être un organe de commande connu en soi, une charge 2 et un thermo- couple 4 pour mesurer la température de l'enceinte La sortie de ce thermocouple 61 est reliée à l'entrée d'un convertisseur analogiquenumérique 51 qui est relié à une ligne bus d'entrée 55 d'un microprocesseur 53 Le micropro- cesseur est relié par une première ligne bus d'entrée-sortie 57 à des mémoires externes 54 et par une deuxième ligne bus 58 à une face avant 81 comportant un clavier de touches 82, un élément d'affichage 83 et une prise de dialogue 84 qui peut être reliée par un câble de liaison entréesortie 91 représenté en pointillé à une console 10 Le rôle de la liaison 91 avec la console sera expliqué par la suite Le microprocesseur ou microsystème est relié par une ligne bus de sortie 56 à l'entrée d'un convertisseur numérique analo- gique 52 dont la sortie est branchée par la ligne 62 à l'entrée du gradateur 11 Le microprocesseur peut également comprendre au moins une entrée tout ou rien 59 reliée par un contact 60, servant à signaler par exemple, l'ouverture de la porte du four, à une source de tension V. Le dispositif est utilisé en deux phases, une première dite -6- d'identification et une deuxième dite de commande. Le dispositif est utilisé dans la deuxième phase dite de commande, sans la liaison 91 avec la console 10 Au début de cette phase les mémoires 54 sont chargées par un modèle interne de représentation de l'enceinte constitué par un premier ensemble de valeurs stockées dans un premier bloc de mémoires, cet ensemble de valeurs étant les coefficients représentatifs de la courbe de réponse impulsionnelle de l'enceinte, pris à des intervalles de temps donnés, définis- sant les adresses de stockage de ces premières valeurs; un deuxième ensemble de valeurs représentatives de la courbe de variation du gain G des réponses impulsionnelles pour diffé- rentes températures de fonctionnement 9 j définissant les adresses de ces deuxièmes valeurs; un troisième ensemble de valeurs représentatives de la courbe de variation des périodes d'échantillonnages T permettant le même nombre de points de définition que pour les premières valeurs, pour différentes températures 9 j définissant les adresses de ces troisièmes valeurs Ces deux dernières courbes permettent par deux homothéties ou par d'autres transformations mathé- matiques appliquées à la première courbe, de reconstituer par calcul toutes les courbes de réponse impulsionnelle de l'enceinte en fonction des différentes températures de fonctionnement comprises entre une température de fonction- nement maximum et une température minimum, et d'obtenir ainsi une représentation du modèle mathématique du processus en boucle ouverte. On utilise ensuite ce modèle stocké dans les mémoires 54 avec un modèle de référence donnant le comportement souhaité du système bouclé pour déterminer par calcul la succession de signaux de commande à appliquer au gradateur, en fonction du temps de réponse souhaité pour le système bouclé de la température de consigne à atteindre, et de la température de départ qui est celle mesurée dans l'enceinte par le thermo- couple 4 ou un capteur de température Le microprocesseur ou -7microsystème contient un programme qui élabore en fonction des contraintes que l'on veut imposer au système (temps de réponse bouclé, gradient de température à ne pas dépasser, etc) une stratégie de commande qu'il va appliquer au modèle en boucle ouverte et à l'enceinte Du modèle en boucle ouverte il recevra la température que doit atteindre l'enceinte et de l'enceinte il recevra la température mesurée La comparaison de ces deux températures permet au microprocesseur de modifier si besoin sa stratégie de commande dans le cas d'une désadaptation. Les algorithmes utilisés pour élaborer la stratégie de commande sont bien connus des automaticiens et seront expli- cités plus loin. Les signaux de commande W sont envoyés sous forme numérique par le bus 56 au convertisseur numérique analogique qui transmet un signal de commande analogique à l'organe de commande 11 La température de l'enceinte évoluant en réponse au signal de commande est mesurée par un thermocou- ple et transmise sous forme analogique à un convertisseur analogiquenumérique qui transmet au bus d'entrée 55 un signal numérique Les touches 82 et l'affichage 83 permet- tent en entrée de faire prendre en compte par le micropro- cesseur la température de consigne à atteindre et le temps de réponse souhaité pour le système bouclé, et les contrain- tes à imposer à l'enceinte qui peuvent être de niveau (valeur à ne pas dépasser), ou de gradient (contrainte sur les incréments) L'affichage permret entre autre, de visualiser en sortie, la température du four. Pour pouvoir utiliser le dispositif dans la phase dite de commande il est nécessaire de procéder lors de la mise en route du régulateur à une première phase dite d'identifica- tion Cette phase d'identification consiste à élaborer le modèle mathématique de représentation de l'enceinte du processus en boucle ouverte. 8 - On choisit de représenter le modèle par la réponse impul- sionnelle de l'enceinte Pour obtenir cette réponse impulsionnelle on procède en quatre phases. La première phase consiste en la détermination approximative du temps de réponse de l'enceinte Pour cela la console envoie un signal en forme d'échelon à l'enceinte et enregistre et stocke en mémoire la courbe de réponse de l'enceinte aux stimulis envoyés par la console Puis la console calcule, par exemple, par la méthode de Naslin ou la méthode de la pente à l'origine, le temps de réponse du processus en boucle ouverte et passe à la deuxième phase. Dans la deuxième phase, la console élabore un protocole d'identification Ce protocole consiste en une succession de créneaux de largeur variable et de-hauteur variable dont les largeurs restent comprises entre une largeur maximale et une largeur minimale qui sont en relation avec le temps de réponse de l'enceinte Ceci afin d'éviter d'envoyer à l'enceinte des signaux d'excitation provoquant une réponse non significative pour la modélisation de l'enceinte La console détermine également la période d'échantillonnage (T) en fonction du temps de réponse ( lr) de l'enceinte et du nombre (N) de points souhaités pour la réponse impulsionnel- le La période d'échantillonnage T est calculée de façon-que N x T > T tout en respectant le théorème de Shannon pour le choix de T. Une fois le protocole élaboré, la console passe à la troisième phase qui permet la constitution du modèle Pour cela la console applique le protocole d'identification qui consiste en une succession de créneaux à l'enceinte qui, par échantillonnage, enregistre en mémoire les réponses au protocole et effectue les calculs pour déterminer la réponse impulsionnelle de l'enceinte. Pour effectuer ces calculs la console utilise un algorithme 9 d'identification tel que: h (n+l)= h (n) _ (n) e(n-i) ( N + l = h 1 E En) -(n) avec e(n) = SM (n) -(n) =n M(n) en-i men- i= 1 h(n) = ime coefficient de la réponse impulsionnelle de l'enceinte SM: sortie modèle SO: sortie objet O ET = len' e e e n-2 nlenî 11 n, en_ 2, en N et En e n-i en-N E est le vecteur entrée constitué par les valeurs des entrées aux instants d'échantillonnage n-l, n-2 n-N. c (n) est égal à la différence entre la sortie du modle à l'instant d'échantillonnage N et la sortie de l'enceintedite sortie objet. est un coefficient dont le choix permet d'obtenir une convergence plus ou moins rapide et donc d'arriver plus ou moins rapidement à la détermination de la valeur h M (n) représentant la ieme composante du vecteur HM à l'instant d'échantillonnage n Le vecteur HM est la représentation mathématique du modèle constitué par les coefficients de la réponse impulsionnelle du modèle de l'enceinte HMT = hi h MM pour un point de fonctionnement = 1h N donné. On a donc E(n) donné à un instant d'échantillonnage par la console et SO(n) mesurée à l'instant d'échantillonnage n. Connaissant par les calculs précédents h(n), on en déduit C (n) et h 1 (n + 1) et on détermine ainsi de proche en proche par une succession de mesures effectuées par échantillonnage sur le créneau appliqué à l'entrée et sur le signal récupéré - en sortie de l'enceinte et par une succession de calculs les composantes définissant la réponse impulsionnelle de l'enceinte. On opère de la même manière pour une succession de créneaux dont la valeur moyenne se situe à différents niveaux corres- pondant à différents points de fonctionnement On obtient ainsi une famille de vecteurs représentant une famille de courbes telles que celles représentées à la figure 3 Cette famille de courbe Gj Nl hij On procède de même pour la période d'échantil- lonnage en calculant les valeurs de la période d'échantil- lonnage T pour que compte tenu du temps de réponse aux différentes températures de fonctionnement T Qj on ait toujours le même nombre de points d'échantillonnage N pour une courbe donnée On obtient ainsi la variation de la période d'échantillonnage en fonction de la température du point de fonctionnement La période d'échantillonnage est déterminée par T = T Et aui permettrait de représenter la N - courbe de la figure 4 Cette succession de coefficients représentant une réponse impulsionnelle échantillonnée pour une température de fonctionnement donnée, les variations du gain en fonction de la température et les variations de la période d'échantillonnage en fonction de la température sont stockées dans la mémoire de la console jusqu'à la fin de l'identification et transférées dans la mémoire 54 du régulateur à la fin de la phase d'identification Cette représentation paramétrée du modèle permet de reconstituer par homothétie, ou par toute autre transformation mathémati- que, les courbes de réponses impulsionnelles aux différentes températures de fonctionnement. L'étape d'identification se termine par une quatrième phase de validation qui consiste à faire générer par la console un autre protocole différent du premier, à calculer par le modèle les valeurs des sorties et à vérifier que les valeurs mesurées correspondent aux valeurs calculées A la fin de cette phase, la console transfère le modèle dans le régulateur et le régulateur peut fonctionner seul pour commander l'enceinte sans que la console soit branchée. Pour la régulation le microprocesseur utilise donc le modèle interne représentant le système en boucle ouverte, lun modèle de référence qui représente le comportement -désiré du système bouclé et un algorithme de commande permettant de calculer les commandes à appliquer à l'enceinte compte tenu des commandes présentes et passées de façon qu'il y ait coïncidence entre le comportement souhaité calculé et le comportement réel du système. Dans l'exemple considéré on a choisi comme représentation du modèle de référence: SMR (n+l) = a SO(n) + (lc-a) CONS 3 ( 2 ' et SMR (n+h) = SMR(n+ h-l) + (L-é) CONS;h a 2 -T avec CL = e h T période d'échantillonnage, Tr temps de réponse souhaité du système, SMR sortie du modèle de référence, SO sortie objet mesurée, CONS consigne. L'algorithme de commande est l'algorithme dual de l'identi- fication. 12 - La sortie du processus étant calculée par la formule faisant M intervenir les coefficients h du modèle interne SM In) = eni (n) hi M Donc connaissant par la mesure SO(n) à l'instant d'échantil- lonnage n, on calcule les sorties du modèle de référence souhaitée aux instants N + 1, N + 2, N + h par les formules ( 2) dans lesquelles on connait la sortie objet à l'instant N et la consigne à atteindre CONS. On estime l'évolution future de la sortie en calculant: T * SMIF(n+h 3 = HM Eh * j Eh est le vecteur commande compose de h commandes futures et N-h commandes passées. Connaissant la sortie du modèle interne futur on calcule E(n+h) = SMIF(n+ h) SMR(n+h) Les commandes futures sont ajustées de façon itérative, chaque itération se décomposant en h ajustements pour satisfaire les h points de colncidence souhaités, en commençant par le plus éloigné dans le futur. Cet ajustement se fait, pour chacune des composantes du vecteur commande futur, en chacun des points de coïncidence suivant l'algorithme *j = j X(n + h i + 1) h. ej =j HM HM avec j variant de 1 à h i + 1 i variant de 1 à h coefficient de relaxation ej = valeur de la jeme commande future obtenue soit à l'itération précédente, soit par décalage pour la lère itération * ej prenant la place de ej à l'itération suivante h M = coefficient du modèle interne j= E(n + h i + 1)= écart prédit à l'instant N au h i + 1 point de coïncidence. On applique la première commande à l'enceinte, cette 12985 13 - commande provoquant une sortie SQ (n+l) Connaissant SQ (n+l) on comparera avec la sortie prédite et on réitérera le calcul dans le cas d'un écart, dans le cas de la coinci- dence on applique la commande future calculée Les opérations se sont donc déroulées de la façon suivante: mesure de la sortie du processus, calcul du comportement futur souhaité, estimation de l'évolution future de la sortie, ajustement des commandes futures, application de la première commande, mesure, comparaison avec les valeurs estimées et nouvelle phase de calcul en cas de désadapta- tion, sinon on applique les valeurs du premier calcul. La figure 7 représente une variante de régulation utilisant le principe connu de régulation en cascade des résistances de chauffage et de l'enceinte Cette variante ayant pour but d'améliorer les performances thermiques de l'enceinte et de protéger les résistances de chauffage. Dans la figure 7 on retrouve les éléments du dispositif de la figure 2 avec les mêmes références; mais, dans le micro- processeur 53, on a rajouté une deuxième entrée constituée par le bus d'entrée 123 qui est relié par un convertisseur analogique numérique 122 à un thermocouple 12 mesurant la température de la résistance de chauffage du four électri- que Cette deuxième entrée est utilisée par le micro-processeur pour calculer les commandes W à appliquer au gradateur en fonction de la commande e déterminée de la façon décrite ci-dessus et de la température de la résistan- ce de chauffage TRA. Pour cela le microprocesseur utilise l'algorithme suivant, mais tout autre algorithme peut être utilisé. (n) (n-l) + G Pre(n) -(nl) + GIT (n-1) avec E(n) =e(n) TRA(n) e(n) consigne température résistance 14 - Dans cette variante chaque fois que le microprocesseur a déterminé la valeur en à appliquer comme précédemment on calcule la valeur W(n) qui sera appliquée à la commande de la résistance chauffante Cet algorithme est la réalisation numérique d'un régulateur P I (proportionnel, intégral) dont le gain proportionnel est GP et le gain intégral GI; T représente la période d'échantillonnage. Cette variante permet de tenir compte du décalage existant entre la montée en température des résistances chauffantes et la montée en température de l'enceinte du four, et d'éviter de soumettre les résistances chauffantes à des contraintes pouvant entraîner leur destruction. Il va de soi que toute variante ou modification apportée aux dispositifs et procédés décrits, à la portée de l'homme de l'art, font également partie de l'esprit de l'invention. - Revendications de brevet. 1 Procédé d'identification du modèle mathématique d'une enceinte thermique, caractérisé en ce qu'il comprend les phases suivantes: détermination approximative du temps de réponse de l'enceinte par mesure et enregistrement de la réponse en température de l'enceinte à un échelon de puissance, élaboration d'un signal d'identification constitué d'une multiplicité de créneaux aléatoires de hauteur et de largeur, comprise entre une limite inférieure et une limite supérieure en rapport avec le temps de réponse du four, variant de façon aléatoire, mesure par échantillonnage, calcul et mise en mémoire de la réponse impulsionnelle de l'enceinte au signal d'identi- fication pour chaque créneau, calcul du modèle de représentation de l'enceinte à partir des réponses impulsionnelles mémorisées. 2 Procédé selon la revendication 1, caractérisé en ce que le modèle de représentation de l'enceinte est constitué: par l'enregistrement dans une première série de mémoires d'une première multiplicité de premières valeurs détermi- nées après échantillonnage à des intervalles de temps constants, lesdites valeurs étant représentatives de la réponse impulsionnelle de l'enceinte en fonction du temps, lesdites valeurs étant rangées à des adresses correspondant au rang de l'échantillonnage, par l'enregistrement dans une deuxième série de mémoires d'une deuxième multiplicité de valeurs représentatives de la 16 - variation du gain des réponses impulsionnelles en fonction de la température de fonctionnement, lesdites deuxièmes valeurs étant rangées à des adresses correspondant aux différentes températures de fonctionnement, par l'enregistrement dans une troisième série de mémoires d'une troisième multiplicité de valeurs représentatives de la variation de la période d'échantillonnage en fonction des différentes températures de fonctionnement, à des adresses correspondant aux températures. 3 Procédé selon la revendication 1 ou 2, caractérisé en ce qu'il comprend une phase de validation et de contrôle du modèle consistant à appliquer des signaux d'action à l'enceinte, à calculer au moyen du modèle la réponse aux signaux de l'enceinte, et à comparer la réponse calculée avec la réponse mesurée. 4 Procédé selon l'une des revendications 1 à 3, caractérisé en ce que le calcul de la réponse impulsionnelle et du modèle est élaboré par une console qui utilise un algorithme d'identification tel que h.(n+l) = h (n) (n)T(n (n) (n) Application du procédé d'identification du modèle d'une enceinte à la régulation de l'enceinte comprenant de façon connue un élément de mesure 4 de la variable à réguler, un élément de régulation 5 élaborant a partir des valeurs de variables et du point de consigne l'action à imposer à l'enceinte, un organe de commande 11 agissant sur l'enceinte en réponse au signal d'action, caractérisée en ce que l'élément de régulation élabore le signal de commande par un premier ensemble d'algorithme qui tient compte de la mesure de la sortie de l'enceinte, calcule le comportement futur souhaité en fonction du temps de réponse souhaité du système bouclé, estime grâce au 17 - modèle l'évolution future de la sortie et ajuste les com- mandes futures, applique la première commande, mesure et compare avec les valeurs estimées la mesure pour appliquer la commande suivante en cas d'adaptation ou effectue un nou- veau calcul, une nouvelle estimation et un nouvel ajustement avant d'appliquer la nouvelle première commande. 6 Application du procédé selon la revendication 5, caractérisée en ce que l'élément de régulation est constitué par un microprocesseur ou un microsystème relié pour les entrées à un convertisseur analogique digital recevant le signal d'au moins un capteur de température, pour la sortie du signal de commande ainsi qu'à des mémoires de stockage du modèle mathématique. 7 Application du procédé d'identification du modèle d'un processus à la régulation selon la revendication 6, caractérisée en ce que le microprocesseur ou microsystème est relié lors de l'identification du modèle à une console incorporée ou non élaborant le modèle et venant charger les mémoires de stockage du régulateur avec le modèle. 8 Application du procédé selon la revendication 6 à la régulation avec double boucle, caractérisée en ce qu'un deuxième capteur mesure un paramètre de l'actionneur et fournit un signal au micropro- cesseur ou microsystème utilisé par un deuxième algorithme qui calcule les commandes à appliquer à l'actionneur en fonction du paramètre mesuré de cet actionneur et des commandes calculées par le premier ensemble d'algorithmes. 9 Application du procédé selon la revendication 5, caractérisée en ce que le premier ensemble d'algorithmes est constitué par SMR(n+h) = ct SMR(n+hl) + (l -ea) CONS 18 - -T avec a= e Tr et SMR(n+l) =O SO + ( 1 a) CONS * x-M s ke =e -_À ú (n+hi+l) hj ej =l e HMT HM avec SM avec (n+h i+l) SMIF(n+h i+l) SMR(n+h i+l) SMIF =HM Eh (n+h) h Application du procédé selon la revendication 8, caractérisée en ce que le deuxième algorithme peut être un algorithme P I numérique de la forme W(n) = W(n_) + GP (n) (n-l)l + GI T (n-l) avec ú (n) = e(n) TRA(n) GP: Action proportionnelle GI: Action intégrale