Brevet FR2931003A1 — publié le 20091113

La présente invention concerne un procédé ainsi qu'un dispositif de rejet robuste de perturbations périodiques dans une boucle de contrôle de position d'un axe. Elle a des applications dans le domaine des automatismes dans lesquels des moteurs électriques commandés sont mis en oeuvre et en particulier au sein de boucles de rétrocontrôle ou d'asservissement. Elle s'applique en particulier dans les simulateurs de mouvements. Dans le domaine du contrôle d'axes, les moteurs en courant alternatif sans balai (� AC brushless �) sont aujourd'hui très largement utilisés. Il s'agit en fait de moteurs synchrones auxquels on adjoint une électronique assurant un autopilotage du moteur. Cet autopilotage a pour effet de rendre le couple développé par le moteur proportionnel à une grandeur d'entrée de commande de l'électronique. Cette linéarité entre le couple et la commande est une propriété très appréciable car elle permet l'utilisation d'algorithmes linéaires pour le contrôle de la vitesse ou de la position de l'axe du moteur. De plus, l'utilisation d'un moteur sans balais garantit un entretien minimal de la machine. Le mouvement de l'axe asservi peut être soumis à des perturbations périodiques dépendantes de la position dudit axe. Ces perturbations peuvent notamment provenir des ondulations de couple du moteur autopiloté qui ne sont jamais totalement réductibles, malgré les efforts des constructeurs. L'origine des ondulations de couple est multiple. Celles-ci peuvent provenir du couple de � cogging � qui est du à l'interaction entre les aimants permanents et les dents du stator, mais aussi provenir du couple de réluctance, du couple mutuel, ou encore des décalages (� offsets �) en courant continu des capteurs de courant dans l'électronique de commande. En pratique, l'origine prépondérante dépend du type de moteur et du constructeur. Quoi qu'il en soit, l'ondulation de couple est un phénomène périodique en fonction de la position angulaire de l'axe et sa fréquence est proportionnelle à la vitesse angulaire dudit axe. Une autre source de perturbation périodique sur un axe se trouve au niveau du capteur mesurant la position dudit axe. Dans le cas où l'on utilise un codeur incrémentai, des méthodes à base de réseaux optiques permettent de connaître la position de l'axe entre deux traits du codeur. Cette méthode dite d'interpolation de la position n'est pas parfaite et induit une non linéarité sur la mesure de la position angulaire quasiment reproductible d'un trait de codeur à l'autre. Dans ce cas aussi, la non linéarité induite par ce phénomène est périodique et proportionnelle à la vitesse angulaire de l'axe. Dans une boucle d'asservissement classique, où la variable contrôlée est soit la position, soit la vitesse angulaire de l'axe, ces perturbations périodiques entraînent des ondulations sur la vitesse angulaire de l'axe qui peuvent se révéler gênantes sur certaines applications. C'est le cas notamment des simulateurs mécaniques de mouvements qui sont des systèmes destinés en général au test de centrales inertielles et où la stabilité de vitesse est un élément déterminant de la qualité du système. Afin de limiter l'effet de ces perturbations périodiques, il a été proposé plusieurs méthodes permettant leur atténuation. Elles peuvent être subdivisées en trois principales catégories : Dans la première catégorie, c'est la conception de la machine qui est optimisée. Par exemple en ce qui concerne le moteur on élimine par construction notamment le couple de � cogging � et le couple réluctant. Ce travail ne peut être fait que par le constructeur du moteur. Dans la seconde catégorie, des méthodes de commande adaptative sont utilisées. Dans ces commandes adaptatives, l'amplitude et la phase des ondulations sont estimées au moyen d'un algorithme approprié, et elles sont utilisées dans un contrôleur de moteur pour compenser les perturbations. Il est à noter que ces méthodes ne sont efficaces que si les amplitudes des ondulations sont constantes sur un tour, ce qui n'est pas forcément le cas. Dans le cas contraire, on aboutit seulement à une compensation partielle. Dans la troisième catégorie, on range les algorithmes d'asservissement de l'axe en vitesse ou en position angulaire qui incluent un rejet de perturbation sinusoïdale à la fréquence des ondulations. Comme ladite fréquence est variable en fonction de la vitesse angulaire, le correcteur utilise forcément des gains séquencés en fonction de cette vitesse angulaire. Néanmoins, dans ce cas, le rejet des ondulations se fait de façon robuste, c'est-à-dire qu'il n'y a pas besoin de connaître la phase et l'amplitude des perturbations sur le couple afin de les éliminer. Ainsi, la/les fréquences de perturbation étant proportionnelles à la vitesse de l'axe, pour estimer cette/ces fréquences, on peut se baser sur la vitesse de consigne. Le rapport de proportionnalité entre la fréquence et la vitesse est constant et dépend de la construction du moteur. La présente invention se range dans cette troisième catégorie et permet notamment le rejet des perturbations périodiques dont la fréquence dépend de la vitesse de rotation d'un axe moteur asservi. Appartenant à cette troisième catégorie on peut signaler le document antérieur publié dans la revue � IEEE/ASME Transactions on mechatronics, vol.9 n°2, june 2004, sous le titre � torque and velocity ripple elimination of AC permanent magnet motor control systems using the internai mode) principle � et dont les auteurs sont Wai-Chuen et Li Qui. Ces auteurs proposent une méthode de contrôle de la vitesse d'un axe basée sur un correcteur à deux degrés de liberté de type RST, dont les blocs RST sont adaptés en fonction de ladite vitesse, de façon à rejeter les perturbations statiques et les perturbations sinusoïdales, à fréquence variable, des perturbations périodiques. Toutefois, ce document ne peut directement conduire à une application industrielle pour diverses raisons : - Les auteurs ne présentent pas une méthodologie de réglage du correcteur permettant d'assurer qu'à toutes les vitesses angulaires, le correcteur garantit un niveau de marges statiques et dynamiques nécessaires pour assurer la robustesse de la boucle. - Le montage des auteurs s'applique à une boucle de vitesse. Si l'on transpose le problème à une boucle de position, en ayant toujours en tête le contrôle de la vitesse, l'ordre du correcteur augmente. Or tous les paramètres des blocs RST devant être calculés en ligne, le volume de calculs à effectuer à chaque période d'échantillonnage risque de devenir rapidement rédhibitoire surtout si la période d'échantillonnage du dispositif est petite, comme l'on en trouve souvent dans le domaine du contrôle d'axes. - Le rejet des perturbations est systématique c'est-à-dire qu'il n'est pas laissé le choix d'activer ou non le rejet de la perturbation sinusoïdale. Afin de mieux faire comprendre l'invention, une description d'un dispositif correcteur à deux degrés de liberté de type RST classique avec filtre de consigne (préfiltre) est maintenant effectuée ainsi qu'une présentation de principes généraux en relation avec les problèmes et solutions en rapport avec l'état de la technique et l'invention. Le dispositif correspondant de correcteur RST avec filtre de consigne est symbolisé sur la Figure 1 où l'on a les paramètres suivants : c(t): consigne en position, w(t): sortie du filtre de consigne, 8(t) : position mesurée de l'axe, u(t) : commande en entrée du système, v(t) : perturbation de charge, p(t) : perturbation sur la sortie. Dans le cadre général des correcteurs RST, on peut signaler la demande de brevet FR07/52807 pour � Optimisation de la réponse fréquentielle d'un simulateur de mouvement par suivi adaptatif de consignes sinusoïdales � dans laquelle un filtre de consigne est mis en oeuvre dans une structure RST. Sur la Figure 1, la fonction de transfert Rr(z-1)/Fr(z-l) constitue le filtre de consigne destiné, comme son nom l'indique, à filtrer la consigne, dans le cas où celle-ci serait par exemple trop abrupte pour le correcteur. Soit un axe à commander relié à un moteur synchrone autopiloté et à un codeur de position mesurant la position dudit axe. Compte tenu de l'autopilotage, la commande u(t) en entrée de l'électronique de commande est proportionnelle au couple développé par le moteur. La position de l'axe du moteur est notée 8(t) . On considère que l'ensemble des frottements secs et visqueux sont des perturbations agissant sur le couple. Dans ces conditions, en appliquant le principe fondamental de la dynamique, la fonction de transfert entre la commande u(t) et la position 8(t) s'écrit simplement : 8(t) b u(t) s2 avec b proportionnel à lm et J l'inertie du système. Le système est échantillonné à la période Te. Soit z l'opérateur avance d'une période d'échantillonnage. La fonction de transfert du système discrétisé muni d'un bloqueur d'ordre 0 s'écrit : 0(t) _ 0.5•b•Te2 z 1 +0.5•b•Te2 z-2 u(t) 1ù 2z-1 + z-2 La fonction de transfert entre w(t) et 8(t) s'écrit : e(t) _ z-1T(z-1)B(z-1) w(t) - A(z-1)S(z-1)+z-1B(z-1)R(z-1) En présence de perturbations sur la sortie et des perturbations de charge, on définit les fonctions de sensibilité suivantes du système : Sep : Fonction de transfert entre la perturbation sur la sortie p(t) et la sortie 8(t) Sev : Fonction de transfert entre la perturbation de charge v(t) et la sortie 8(t). De la Figure 1, on déduit les relations suivantes : e(t) A(z-1)S(z-1) S û p(t) û A(z-1)S(z-1)+z-1B(z-1)R(z-1) S(z1) _ S~ v(t) A(z-1)S(z-1)+z-1B(z-1)R(z-1) Soit w(t) la vitesse angulaire de l'axe, c'est-à-dire la dérivée de la position angulaire 8(t), on a la relation : w(t) 1 • 9(t) Te On peut ainsi définir une fonction de transfert supplémentaire entre la perturbation de charge v(t) et la vitesse angulaire de l'axe w(t), nommée S,,. On peut écrire : w(t) û S(z-1)•(1ûz-1)/Te S _ v(t) A(z-1)S(z-1)+z-1B(z-1)R(z-1) Par ailleurs, on peut montrer que la fonction de transfert entre la A perturbation de charge v(t) et les variations relatives de vitesse angulaire (t)) de l'axe est égale à s,,,, Le réglage du correcteur se fait classiquement par la spécification d'un polynôme de stabilité arbitraire D(z-1) puis on calcule S(z-1) et R(z-1) tels que l'on ait : A(z-1)S(z-1)+z-1B(z-1)R(z-1) = D(z-1) Cette dernière équation étant une équation polynomiale dite équation de Bézout dont la résolution passe par la résolution d'un système linéaire dit système de Sylvester. On peut noter que le degré de D(z-1) est forcément au moins égal à 2n, où n est le degré de A(z-1) . Afin d'assurer le rejet robuste des perturbations sur la commande ou sur la sortie, il est nécessaire, dans l'équation de Bézout de pré-spécifier le bloc S. C'est-à-dire qu'il faut que le bloc S puisse s'écrire sous la forme : S(z-1) = Hs(z-1) • S'(z 1 ) On appelle cette méthode de pré-spécification, la méthode du � modèle interne �. L'équation de Bézout devient alors : A(z-1)Hs(z-1)S'(z-1) + z-1B(z-1 )R(z-1) = D(z-1) Les inconnues en sont maintenant S'(z1) et R(z 1). Si l'on souhaite rejeter d'une façon robuste les perturbations statiques il faut prendre : Hs(z-1) =1ûz-1 Si l'on souhaite rejeter une perturbation sinusoïdale à la fréquence f, il faut prendre Hs(z-� =1ù2cos(2TC• f •Te)•z-1 +z-2 Si l'on souhaite rejeter une perturbation statique et une perturbation sinusoïdale, il faut que Hs s'écrive : Hs(z-1)=(1ùz-1) .(1ù2cos(2z• f .Te)•z-1 +z 2) On voit donc que l'on peut rejeter ainsi des perturbations de nature diverses. En particulier ce rejet se traduit sur la courbe du module de la fonction de sensibilité Sep (d'une perturbation sur la sortie vers la sortie) par un module égal à zéro à la fréquence spécifiée du rejet comme on peut le voir sur les Figures 2 et 3. En particulier, on peut montrer que si l'on considère que le système est un double intégrateur et que HS contient un zéro en +1, alors, le module de la fonction de sensibilité Sep présente une pente de +60dB/décade dans les basses fréquence et une pente de OdB/décade en haute fréquence. Sous ces mêmes conditions, les pentes du module de Sev sont de +20 dB/décade en basse fréquence et -40 dB/décade en haute fréquence. L'allure du module de Sev en fonction de la fréquence qui en est déduite est représentée Figure 4. La fonction correspondante admet un maximum et on voit qu'il existe une zone fréquentielle critique où les perturbations de charge sont particulièrement amplifiées. On pourrait penser à première vue que l'on peut ainsi rejeter les perturbations à toutes les fréquences que l'on désire. Mais il n'en malheureusement rien, à cause d'un théorème fondamental de l'automatique qui est le théorème de Bode-Freudenberg-Looze. Ce théorème de Bode-Freudenberg- Looze s'écrit : 1n (j.2,z • f )1). df = 0 0 Ce théorème est valable pour les systèmes dont le degré relatif (différence entre le degré du numérateur et du dénominateur de la fonction de transfert) est au moins égal à 2. Il conduit à ce que sur la fonction du module de la sensibilité Sep, qui est représentée schématiquement sur la Figure 5, les aires I et Il sont forcément égales. En conséquence, toute � encoche � pour rejection créée dans la fonction de sensibilité Sep entraîne nécessairement une augmentation du module de cette fonction à une autre fréquence. On ne peut donc pas multiplier le rejet de perturbations sinusoïdales indéfiniment. La méthode de contrôle exposée dans l'article de Wai-Chuen et al. consiste à supposer que les perturbations périodiques sont strictement sinusoïdales. De plus la connaissance de la vitesse angulaire w de l'axe permet d'inférer la fréquence f de ces perturbations. A partir de là, Wai-Chuen et al. proposent de calculer en ligne les coefficients des blocs RST du correcteur suivant le principe du modèle interne par résolution à chaque période d'échantillonnage de l'équation de Bézout décrite précédemment : A(z-� Hs(z-� +z-lB(z-�R(z-� = D(z-� Avec : Hs(z-� = (1-z-�.(1-2cos(2TC• f •Te)z-1 +z-2) Dans cette équation de Bézout, A(z-'), B(z-') sont fixes (paramètres de la fonction de transfert du système à commander). Seul le polynôme Hs(z-')est variable dans le temps, en fonction de la fréquence de perturbation sinusoïdale. Enfin, D(z-'), polynôme des pôles de la boucle fermée, est également choisi fixe. Le polynôme R(z-') est ainsi calculé en ligne. On peut remarquer que les auteurs précités réalisent le contrôle en vitesse de l'axe. C'est-à-dire que dans le cas qui les concerne, le système à commander est d'ordre 1, il en résulte que le degré du polynôme R(z-1) à calculer en ligne est relativement simple puisqu'il est d'ordre 3. Par ailleurs, le nombre de pôles placés au moyen du polynôme D(z-1) est de 4. Or, dans le cas où l'on cherche à faire un asservissement de position, le système à commander est d'ordre 2. Il en résulte que les degrés de R(z-'), S(z-l), T(z-1) sont alors plus élevés, si bien que le volume de calculs requis pour la résolution de Bézout augmente sensiblement, ce qui peut poser des problèmes au niveau du calculateur si la période d'échantillonnage est courte. Par ailleurs, le choix de Wai-Chuen et al. de maintenir fixe le polynôme D(z-l) quelque soit la fréquence de perturbation à rejeter est limitant, dans la mesure où l'on optimise pas la robustesse de la boucle en fonction de cette fréquence de perturbation. Or, il est important dans un contexte industriel de garantir qu'un niveau minimal de marge de retard et de marge de module de la boucle fermée soit assuré quelque soit ladite fréquence de perturbation. En particulier, il est souhaitable, dans le cas où le rejet d'une perturbation entraînerait inévitablement des marges inférieures à celles spécifiées, que ce rejet ne soit pas réalisé. Cette dernière possibilité conduit à ce que le rejet d'une perturbation sinusoïdale puisse éventuellement ne pas avoir lieu à certaines fréquences, si la nature de ladite fréquence entraîne inévitablement une dégradation inadmissible des marges de robustesse. Il se pose ainsi également un problème de conservation d'une marge de retard admissible et on va maintenant présenter des principes généraux en relation avec la robustesse d'un système bouclé. Afin de mesurer la robustesse d'un système bouclé, on a recours à la notion de marges. Les marges usuelles utilisées sur les systèmes monovariables sont : - la marge de module, et - la marge de retard. La marge de module est la distance du lieu fréquentiel de la boucle ouverte par rapport au point -1, dans le plan de Nyquist. C'est un indicateur global de la robustesse d'une loi de commande linéaire. En général, on spécifie une marge de module supérieure à 0,5. La marge de retard représente la valeur maximale du retard parasite que l'on pourrait introduire dans la boucle avant que celle-ci ne devienne instable. A la différence de la marge de module, il n'y a pas de spécification générale sur la marge de retard. Son choix relève d'un compromis entre la robustesse et les performances attendues en matière de rejet de perturbation. On montre, en appelant MM la marge de module que : C'est-à-dire que la marge de module est égale à l'inverse du maximum du module de la fonction de sensibilité Sep. Aussi, lorsque à partir d'un système central bouclé RST d'un type classique rejetant les perturbations statiques, on cherche de surcroît à éliminer une 25 perturbation sinusoïdale (tout en gardant globalement les mêmes performances de rejet statique), � l'encoche � pour réjection créée à la fréquence correspondante va entraîner une remontée du module de Sep aux autres fréquences et en particulier à celles dont le module est supérieure à 1 (cf. théorème de Bode-Freudenberg-Looze). En conséquence, on peut penser que dans le cas général, la marge de 30 module va être diminuée par l'adjonction de ce rejet de perturbation sinusoïdale (toutes choses étant égales par ailleurs). On constate donc déjà que l'ajout d'un rejet de perturbation sinusoïdale sur un système bouclé rejetant initialement les perturbations statiques, n'est pas sans effet sur la robustesse de la boucle. Les conséquences sont également importantes sur la marge de retard du 35 système bouclé. Dans le cas où l'on suppose que le rejet d'une perturbation MM = 1Sepl_ 1 sinusoïdale s'effectue en totalité à la fréquence fo, on peut montrer que la marge de retard du système bouclé est contrainte par ce rejet. En particulier, si on nomme MR la marge de retard d'un système bouclé, on montre que le rejet parfait d'une perturbation sinusoïdale à une fréquence fo (en hertz) entraîne l'inégalité : MR Il est donc souhaitable de disposer d'une stratégie de contrôle différente de ce qui est connu, également basée sur une structure de correcteur fixe de type RST, mais qui évite notamment les problèmes mentionnés et permette une exploitation effective. De plus, dans le cadre d'une loi de commande à paramètres variant, il est nécessaire de pouvoir désactiver les moyens de rejet (par exemple par un interrupteur en sortie d'un bloc modificateur ou par forçage à une valeur nulle de la sortie d'un tel bloc avec ou non arrêt des calculs pour ledit bloc), en particulier dans le cas où la fréquence des perturbations est trop importante par rapport à un critère tel que celui mentionné précédemment. La présente invention y parvient en ajoutant à une structure de correcteur fixe de type RST, un bloc d'estimation de la perturbation de charge du système, ainsi qu'un bloc de transfert additif dont les paramètres sont calculés en ligne et varient en fonction de la vitesse de l'axe. La sortie du bloc de transfert additif est injectée en point particulier du correcteur RST fixe. Il en résulte qu'il est également possible de décider ou non de l'injection de cette sortie du bloc de transfert additif, c'est-à-dire de rejeter ou non les perturbations sinusoïdales. Ce bloc de transfert additif a pour effet de créer une ou plusieurs atténuations dans la fonction de 1 sensibilité du système entre une perturbation de charge et la sortie du système, aux fréquences des perturbations périodiques, tout en modifiant le moins possible ladite fonction de sensibilité aux autres fréquences. Dans le cadre de l'invention, il est également proposé une méthode de réglage hors ligne des paramètres R S et T associée à une méthode de calcul des paramètres du bloc additif de façon à assurer qu'à toute fréquence de perturbation, un niveau de marge de retard et un niveau de marge de module prédéfinis soient garantis, de façon à assurer la robustesse du système. Dans des variantes, il est possible de mettre en oeuvre en amont du dispositif RST avec blocs d'estimation et de transfert, un préfiltre de consigne. L'invention concerne donc un procédé de rejet robuste des perturbations périodiques dans une structure de commande de moteur en position, les perturbations périodiques étant assimilées à au moins une oscillation sinusoïdale à une fréquence déterminable proportionnelle à la vitesse de l'axe, la structure de commande recevant un signal de consigne en position c(t) mettant en oeuvre une boucle de contrôle de position de l'axe dudit moteur, ladite structure de commande étant du type correcteur RST et étant échantillonnée afin de permettre des calculs numériques dans un calculateur programmable pour produire des résultats destinés à la commande du moteur, le correcteur RST mettant en oeuvre des moyens de calcul correspondant, d'une part, à une fonction 1/So(z-l) en amont du moteur et produisant un signal de commande u(t) et, d'autre part, une fonction Ro(z-l) dans sa boucle qui reçoit en entrée un signal de position 8(t) d'un capteur de position sur l'axe dudit moteur, l'ensemble regroupant le moteur, le capteur de position, les moyens de calculs des fonctions 1/So(z-l) et Ro(z-l) avec bouclage étant appelé correcteur central. Selon l'invention : - on ajoute au correcteur central deux moyens de calcul supplémentaires, le premier étant un bloc estimateur E(z-l) et le second un bloc modificateur Q(z-l), ainsi qu'une boucle supplémentaire, - on modélise le moteur sous forme d'un système modélisé soumis à une perturbation de charge v(t), ledit système modélisé recevant en entrée le signal de commande u(t) et produisant en sortie le signal de position 8(t), le modèle du système appliquant sur le signal de commande u(t) une fonction z 1B(z') pour produire un signal additionné à la perturbation de charge v(t), la résultante étant soumise à une fonction 1/A(z-l) pour produire le signal de position 8(t), - on calcule un signal d'estimation de la perturbation de charge v(t) dans le bloc estimateur E(z-l) recevant en première entrée le signal de commande u(t) et en 11 2931003 seconde entrée le signal de position 0(t), dans ledit bloc estimateur le signal de commande u(t) étant soumis à la fonction z 1 B(z 1) pour produire un signal qui est soustrait au résultat de l'application de la fonction A(z 1) sur le signal de position 0(t), la soustraction donnant le signal d'estimation de la perturbation de charge s v(t) - on applique le signal d'estimation de la perturbation de charge v(t) au bloc modificateur Q(z-') pour produire un signal modificateur appliqué à une entrée négative (il en résulte que le signal modificateur est soustrait de l'autre signal arrivant sur l'additionneur) d'un additionneur supplémentaire disposé en amont du 10 moyen de calcul de la fonction 1/So(z-1) et formant ainsi la boucle supplémentaire, le bloc modificateur ayant une fonction de transfert variable en fonction de la au moins une fréquence déterminable et qui est calculée en temps réel ou quasi-réel sous la contrainte que la fonction de sensibilité de position Sep de la structure de commande comporte un creux pour la au moins une fréquence déterminable. 15 Dans divers modes de mise en oeuvre de l'invention, les moyens suivants pouvant être utilisés seuls ou selon toutes les combinaisons techniquement possibles, sont employés : - la boucle supplémentaire est intérieure à la boucle en Ro(z-') du correcteur central, l'additionneur supplémentaire étant disposé juste en amont du moyen de 20 calcul de la fonction 1/So(z-1), - le correcteur central est calculé (synthétisé) préalablement, hors ligne, pour rejeter des perturbations statiques et les fonctions So et Ro sont mémorisées pour être directement utilisées ultérieurement lors des calculs en temps réel, - lors du calcul du correcteur central dans lequel on met en oeuvre une étape de 25 détermination d'un horizon de filtrage To et d'un horizon de commande Tc tels que Tc A la base, l'invention consiste en l'ajout de deux blocs additifs sur une structure d'asservissement classique de type RST, cette dernière étant représentée Figure 6. Dans la description qui suit, un tel correcteur RST est nommé correcteur central. Le correcteur central, en lui-même, rejette uniquement les perturbations statiques, ses blocs polynomiaux fonctionnels constitutifs sont nommés Ro(z-') et So(z-�. De même, dans cette description, on considère la position de consigne ou la vitesse de consigne selon les cas mais ces valeurs se déduisent l'une de l'autre d'une manière simple car la vitesse de consigne est calculée par dérivation de la position de consigne. Les blocs fonctionnels additifs ajoutés au correcteur central sont au nombre de deux. Le premier bloc fonctionnel additif, E(z-�, produit un estimé v(t) des perturbations de charge v(t) et il est appelé bloc estimateur. Le second bloc fonctionnel additif, Q(z-l), a pour entrée v(t) et produit un signal qui peut être injecté en amont du bloc So(z-l) et ses coefficients sont calculés afin de modifier la fonction de sensibilité Sev du correcteur central de telle sorte à créer une, voire deux, encoches dans Sei,, pour une, voire deux fréquences (le cas de plus de deux fréquence s'avère plus difficile compte tenu du théorème de Bode-Freudenberg-Looze), tout en la modifiant le moins possible aux autres fréquences. Ce second bloc est appelé bloc modificateur. Ce bloc modificateur Q(z-l) possède des coefficients variables dans le temps. Comme représenté Figure 7, le bloc estimateur à ses deux entrées respectivement en amont et aval de la partie système du dispositif RST et sa sortie v(t) est envoyée en entrée du bloc modificateur dont la sortie est envoyée à un soustracteur (ou additionneur avec une entrée correspondante inverseuse ou négative) (la soustraction est effectuée ou non à volonté, du fait de la présence d'un symbole d'interrupteur sur la sortie du bloc modificateur), à l'entrée de la structure RST. On peut remarquer que cette Figure 7, la sortie du bloc modificateur est représentée comme venant se soustraire en aval de l'additionneur (à entrée inverseuse ou négative qui soustrait en fait la sortie de Ro(z-l) à l'entrée de la structure RST recevant la consigne pour constituer la rétroaction). Dans une variante non représentée, la sortie du bloc modificateur vient se soustraire en amont de l'additionneur à entrée inverseuse ou négative recevant la rétroaction Ro(z-l). Ainsi, selon le cas, la boucle supplémentaire est intérieure à la boucle en Ro(z-l) du correcteur central, l'additionneur supplémentaire (à entrée négative) étant disposé juste en amont du moyen de calcul de la fonction 1/So(z-l), ou, la boucle supplémentaire est extérieure à la boucle en Ro(z-l) du correcteur central, l'additionneur supplémentaire (à entrée négative) étant disposé en amont de l'additionneur fermant la boucle en Ro(z-1). En revenant à la Figure 7, on peut montrer que les pôles de la boucle fermée du correcteur central seul sont tels que : D(z-1)=A(z-1)So(z-1)+z-1B(z-1)Ro(z 1) Le bloc modificateur Q(Z 11 peut s'écrire : Q(z-1) = fl(z-1) a(z-1) et sa sortie est injecté dans le correcteur RST. En considérant que le système global bouclé ainsi obtenu, qui est muni de 10 ses deux blocs additifs supplémentaires (blocs estimateur et modificateur), est équivalent à un correcteur RST de blocs R(z-') et S(z-'), on a l'équivalence suivante : R(z-1) = Ro(z-1) + A(z-1) Q(z-1) S(z-1) = So(z-1) - z-' B(z-1) . Q(z-1) Avec: S(z-1) u(t) = ûR(z-1) 9(t) Ce qui aboutit à: So(z-1) • u(t) = ûRo(z-1) • 0(t) -Q(z-1) [A(z-1) • �t)- z-1B(z-1) • u(t)] = ûRo(z-1) .0(t) û Q(z-1) • v(t) Où : v(t) = A(z-1 ) • 9(t) - z 113(z -1 ) u(t) R(z-1) Ro(z-1).a(z-1)+A(z-1). fl(z 1) 20 On a alors : S(z-1) = So(z-1).a(z-1)-z-1B(z-1). fl(z 1) Le polynôme caractéristique de la boucle fermée s'écrit ainsi: Da =A(z-1)•(So(z-1)•a(z-1)_z-1B(z-1)fl(z-1))+z-1B(z-1)•(Ro (z-1)•a(z-1)+A(z-1). (z-1)) Da = D(z-1) • a(z-1) On voit donc que le dénominateur î(z-1) introduit des pôles supplémentaires dans la boucle fermée. 25 Par ailleurs, on peut montrer que la fonction de sensibilité Sepo du système bouclé constitué uniquement du correcteur central peut s'écrire : 15 A(z-�. So(z-� S = epo D(z ' ) La fonction de sensibilité correspondante Sep pour le système bouclé muni des deux blocs additifs supplémentaires (blocs estimateur et modificateur) s'écrit : A(z-�.So(z-� z-'B(z-�. Q(z-� Sep = - D(z-� D(z ) S (z1)=S (z~)_zIB(zI)Q(zl) ~~ D(z-l) La fonction de sensibilité Sep est donc affectée par ces blocs additifs supplémentaires (blocs estimateur et modificateur). On peut montrer qu'il en est de même pour toutes les fonctions de sensibilité de la boucle fermée. 10 Il en résulte que l'on peut utiliser ces deux blocs additifs supplémentaires (blocs estimateur et modificateur) afin de modifier les fonctions de sensibilité de l'asservissement et d'ajouter des pôles au système bouclé. Ces deux propriétés sont utilisées conjointement dans le cadre de l'invention. A partir d'un correcteur central conçu pour rejeter uniquement les 15 perturbations statiques, on adjoint donc ces deux blocs additifs supplémentaires (blocs estimateur et modificateur) pour modifier la fonction de sensibilité Sep (et donc Sei,) afin de rejeter en plus une ou deux perturbations sinusoïdales à une fréquence donnée (outre les perturbations statiques), tout en modifiant le moins possible les fonctions de sensibilité du correcteur central aux fréquences éloignées 20 de celles du rejet. La possibilité de rajouter des pôles à la boucle fermée est utile afin d'optimiser la robustesse de la boucle fermée. On a vu précédemment que l'on a la relation : S(z-1) = So(z-� ùz-1 . B(z-� . Q(z-� Afin que le système bouclé, tel que muni de ses deux blocs additifs 25 supplémentaires (blocs estimateur et modificateur), puisse rejeter les perturbations statiques, mais aussi les perturbations sinusoïdales à au moins une fréquence f1, il est nécessaire que, suivant le principe du modèle interne, le polynôme S(z-l) se factorise en : S(z-1) = Hs(z-1 ).S' (z-l) 30 Le bloc Hs se factorise lui même par : 1ù Z 1 afin qu'il y ait rejet des perturbations statiques. Soit : Hs(z-l) = (1ù z-1) . Hsl(z-l)5 Afin de déterminer les paramètres à utiliser dans le dispositif, le polynôme Hs1(z-l) et les pôles de Q(z-l) (zéros de a(z-l) ) doivent être calculés dans le but de créer une encoche dans les fonctions de sensibilité pour au moins la fréquence f1 de Sep (et donc Sei,), tout en produisant un effet négligeable aux autres 5 fréquences. Une façon de procéder à la détermination de Hsl et a consiste à calculer la fonction de transfert : Hsl(z � a(z-1) par discrétisation d'une fonction de transfert continue qui est dans le cas 10 d'une seule fréquence f1 : ll F(s)= s2+2 nnm •27z.f1 •s+(27z.f1)2 s2 +2 den 27t- . f1 .s+(27t- .f1)2 la discrétisation se faisant par la méthode de Tustin avec � prewarping �. Si l'on souhaite avoir plus de détails sur le choix de de Hsl et a, on peut se référer à l'ouvrage de loan Doré Landau : � Commande des systèmes � éditions 15 Hermès 2002 . L'atténuation Mt de Sep à la fréquence f1 est donnée par la relation : Mt = 20 .log num \ den / Il faut donc que num qui est une équation de Bézout et dont les inconnues sont S�z ) et ,8(z-1) . La résolution de cette équation permet d'obtenir fl(z-1) . Les paramètres a(z-') et (3(z-') étant déterminés, le bloc modificateur est également déterminé. On va maintenant décrire les moyens mis en oeuvre pour déterminer le correcteur central puis, ultérieurement, les moyens de rejet de perturbations sinusoïdales qui ont été vus dans leurs principes juste précédemment. Le procédé de placement de pôles développé par Philippe de Larminat est maintenant succinctement présenté. Ce procédé est applicable dans le cadre de l'invention moyennant quelques adaptations. La finalité du placement de pôles selon Philippe de Larminat consiste à faire rapprocher au maximum la commande RST d'une commande LQG monovariable. Cela suppose de choisir le nombre de pôles de la boucle fermée et de les placer approximativement là ou les placerait une commande LQG. Soit n=deg(A) et ns=deg(Hs). Philippe de Larminat propose que le nombre de pôles de la boucle fermée (qui est de toute façon supérieur ou égal à 2.n) soit égal à 2.n+ns. Les pôles sont alors regroupés en deux groupes : -le groupe des pôles dits de filtrage constitué de n+ns pôles, -le groupe des pôles dits de commande constitué de n pôles. Les termes � filtrage � et � commande � sont utilisés par analogie à une commande par retour d'état reconstruit. La détermination des pôles de filtrage résulte de l'analyse des zéros du polynôme A(z-').Hs(z-� et de leurs modifications ou non en fonction de leurs positions respectives dans le plan complexe. Ainsi, les zéros du polynôme A(z-').Hs(z-') sont projetés sur un cercle de centre zéro et de rayon exp(-Te/To) s'ils sont extérieurs à ce cercle et ils sont laissés inchangés s'ils sont à l'intérieurs ou sur le cercle. On appelle To � horizon de filtrage �, il s'agit d'une variable temporelle, et Te est la période d'échantillonnage. A noter que ce traitement est une transposition au cas discret de la méthodologie de Philippe de Larminat qui a été initialement conçue pour les systèmes continus et dont on peut trouver les détails dans � Automatique, commande des systèmes linéaires �, hermés,1996. Ce traitement s'applique pour les pôles stables ou en limite de stabilité. On peut noter que dans le cas présent, le système est du type double intégrateur et il n'y a donc pas de pôle strictement instable, les pôles peuvent seulement être en limite d'instabilité. La Figure 9 schématise cette détermination dans le cas où il y a projection d'un zéro du polynôme A(z-').Hs(z-� du fait qu'il est initialement à l'extérieur du cercle de rayon exp(-Te/To). La détermination des pôles de commande résulte de l'analyse des zéros du polynôme A(z-') et de leurs modifications ou non en fonction de leurs positions respectives. Le traitement appliqué à ces zéros est semblable au précédent sauf que l'horizon Tc dit � horizon de commande �, est cette fois plus petit que l'horizon de filtrage To. Le traitement consiste donc à prendre les zéros de A(z-') et à les projeter de la même façon sur un cercle de centre 0 et de rayon exp(-Te/Tc) s'ils sont en dehors dudit cercle et de les laisser inchangés dans le cas contraire. A noter encore une fois que ce traitement est une transposition au cas discret de la méthodologie de Philippe de Larminat qui a été initialement conçue pour les systèmes continus et dont on peut trouver les détails dans � Automatique, commande des systèmes linéaires �, hermés,1996. Comme indiqué, Tc qui est l'horizon de commande, est plus petit que To. De préférence, on prend : To/5 Le deuxième terme de l'équation de Bézout du calcul du correcteur RST D(z-� s'écrit : D(z-1) = F(z1)C(z1) En s'inspirant d'une commande par retour d'état reconstruit où l'observateur d'état est inobservable depuis la sortie comme dans la méthodologie de Philippe de Larminat, le polynôme T(z-l) du bloc RST est alors choisi tel que : T(z-� = F(z-� La fonction de transfert entre w(t) et 8(t) s'écrit alors : 8(t) z-lB(z-1) w(t) C(z-� On voit ainsi que les pôles de commande ont une influence claire sur le suivi de consigne. Les pôles de filtrage ont, quand à eux, une influence sur le rejet des perturbations. Si l'on considère à titre d'exemple le cas où le système se comporte en double intégrateur, A(z 1) possède deux zéros en +1. Pour ce qui concerne le rejet des perturbations statiques Hs(z-1) =1û z-1 possède un zéros en +1. Il en résulte lors du traitement des zéros que les trois zéros de A(z-1).Hs(z-l) sont projetés en exp(-Te/To) pour constituer le polynôme de filtrage F(z 1) . Afin de constituer le polynôme de commande C(z 1) , on va projeter les deux zéros de A(z 1) en exp(-Te/Tc). Les polynômes F(z 1) et C(z-l) étant ainsi déterminés, il reste à résoudre l'équation de Bézout afin de déterminer Ro(z 1) et So(z 1). Le procédé de synthèse de correcteur de Philippe de Larminat qui vient d'être présenté et appliqué à système particulier est utilisé dans le cadre de l'invention. Il sert de base à la synthèse du correcteur central et il est complété par un procédé destiné au rejet d'une/de perturbations sinusoïdales. On a vu précédemment que dans le cas où l'on cherche à éliminer une perturbation sinusoïdale à la fréquence f1 en plus des perturbations statiques, il faut que le polynôme Hs puisse s'écrire : Hs(z-1) _ (1û z-1) Hsl(z-1) Par ailleurs, il est nécessaire de rajouter des pôles au système bouclé pour des raisons de robustesse, c'est-à-dire au polynôme D(z 1) , afin que le polynôme de la boucle fermée puisse s'écrire : D(z-1).a(z-1) Hsl et a peuvent être choisis par le calcul indiqué précédemment. Comme on l'a vu, on arrive à la conclusion qu'il faut que num La valeur de den détermine la largeur de � l'encoche � qui apparaît dans la fonction de sensibilité. Une valeur de 0,1 peut être choisie et on peut la fixer une fois pour toutes. On va maintenant décrire à titre de nouvel exemple, le cas où l'on cherche à rejeter simultanément deux perturbations sinusoïdales aux fréquences f1 et f2 (en particulier où f2=2*f1, c'est-à-dire une fréquence fondamentale et un harmonique). Dans ce cas, il faut que la fonction de transfert discrète : Hsl(z-1) a(z-1) soit obtenue à partir de la fonction de transfert continue : F(S) = S2 + 2 num1 • 27t- . f1 • s + (27z • fi )2 s2 + 2 num2 • 27t- • f2 • s + (27z • f2 )2 2 + 2den1 •2?i • f1 • S + (27i • f1)2 S2 s + 2den2 •2?i • f2 • S + (27i • f2 )2 dans laquelle num1 et dent se rapportent à la fréquence f1 et num2 et dent se rapportent à la fréquence f2. 10 Dans ce cas, on peut prendre par exemple dent = den2 , 0,1 et choisir la valeur de num1 et num2 en fonction de f1 et f2. D'une façon générale, on peut dire que si les marges de retard et de module du correcteur central sont MRO et MMO, les marges du système bouclé rejetant en plus une voire deux perturbations sinusoïdales aux fréquences f1 et f2 seront 15 inférieures à MRO et MMO. Si l'on note MR(fl,f2) la marge de module du système bouclé rejetant deux perturbations sinusoïdales aux fréquences f1 et f2, on a la relation : MR(f1,f2) La méthode de réglage du correcteur central consiste en pratique à calculer To et Tc sous la contrainte que quelque soit la ou les fréquences de perturbations sinusoïdales, le respect de marges minimales de retard et de module soit effectif. Soit mrc : la marge de retard � cible �, c'est-à-dire la marge de retard à respecter impérativement quelque soit la fréquence de perturbation à rejeter. Soit mmc : la marge de module � cible �, c'est-à-dire la marge de module à respecter impérativement quelque soit la fréquence de perturbation à rejeter. On cherche à déterminer To et Tc de telle sorte que les performances dynamiques soient les meilleures possible, en tenant compte des contraintes sur mrc et mmc, marges fixées préalablement. L'algorithme présenté ci-après est destiné à être mis en oeuvre au moyen d'un automate programmé (par exemple câblé ou préprogrammé par mémoire morte) ou programmable avec le programme adapté et notamment dans un ordinateur ou calculateur programmable par programme. Pour simplifier, dans cet exemple, on suppose que den-1 et dent sont fixés une fois pour toutes et qu'il en est de même pour les valeurs de ipi et 4%2 qui sont également fixées. De plus, dans cet exemple simplifié, on prend f2=2*f1, ce qui fait que les fonctions qui vont suivre ne sont dépendantes que de f1. On utilise en particulier les fonctions suivantes : -[MR,MM]=marges(To,Tc) Cette fonction détermine la marge de retard MR et la marge de module MM du système bouclé, pour un jeu de paramètres To et Tc, par calcul du minimum des marges de retard sur toute la plage de fréquences de perturbations fi, en fonction de To et Tc donnés. Le calcul de MR se fait par appel à une sous-fonction � marge_ retard_ min �, le calcul de MM se fait par appel à une sous-fonction � marge_ module_ min �. Ces sous-fonctions sont basées sur des calculs d'optimisation qui sont plus efficaces si l'on peut partir de valeurs de départ pas trop éloignées de celles en rapport avec l'optimum recherché. Or on peut montrer que certaines de ces valeurs de départ peuvent être déduites de paramètres fixés ou connus/calculés à priori. En effet, on peut calculer la courbe de la marge de retard en fonction de fi (fréquence du fondamental de la perturbation et dans le cas où f2=2*f1) ainsi que la courbe de la marge de module en fonction de fi, pour un niveau d'atténuation donné du fondamental fi, 4)1 de la perturbation à rejeter et pour un niveau d'atténuation donné de l'harmonique f2, La2, de la perturbation à rejeter. Les courbes correspondantes possèdent chacune des minima locaux et un minimum global. De nombreux essais empiriques ont montré que pour un jeu fixé de paramètres 4)1 et La, la fréquence du minimum global sur chacune des courbes correspondantes varie seulement en fonction de To et Tc, et peut être grossièrement approximée à partir de l'indication de la fréquence du maximum de S6po qui est connue. On a ainsi pu montrer que la fréquence du minimum de MR(fl), pour une configuration donnée de 4)1 et La est grossièrement égale à la moitié de la fréquence du maximum de S6po et des essais extensifs ont également montré que cette propriété est conservé quelque soient To et Tc pourvu que 4)1 et y~2 ne varient pas. Il en est de même pour la courbe MM(fl), le minimum local de la courbe se trouve grossièrement à une fréquence qui est le double de la fréquence du maximum de S6po. Cette dernière propriété est conservée quelque soient To et Tc pourvu que 4)1 et y~2 ne varient pas. Ces indications grossières sur la localisation des fréquences des minima globaux de chacune des deux courbes fournissent un point de départ de l'algorithme d'optimisation (par exemple Newton), dont la convergence vers le minimum permet est assuré si la fréquence d'initialisation de l'algorithme est suffisamment proche du minimum global. On peut ainsi déterminer efficacement par optimisation le minimum de la marge retard ainsi que le minimum de la marge de module lorsque l'on assure le rejet de perturbations sinusoïdales. -[MR]=marge_ retard_min(f1,To,Tc) : Cette sous-fonction calcule le minimum de MR(fl) pour des valeurs données de To et Tc, sur toute la plage des valeurs de f1. Ce calcul peut se faire par optimisation avec un algorithme de descente par exemple l'algorithme de la sécante, à condition d'initialiser l'algorithme avec une valeur de f1 proche de l'optimum qui peut être grossièrement connu en fonction de To, ipi et 1P2. L'algorithme correspondant de calcul du minimum de la marge de retard pour un jeu de paramètre To, Tc donné (yi1 et La ayant été fixés une fois pour toutes préalablement aux calculs) est donné sur la Figure 12. Au départ on utilise une valeur f1(0) fonction de To, ipi et 4%2 et qui est proche de l'optimum comme expliqué ci-dessus. Par des étapes itératives dans lesquelles des variations à divers ordres sont considérées, on détermine le minimum de MR(To,Tc) en prenant comme critère d'arrêt le fait que la variation de f1 d'une itération à la suivante devient, en valeur absolue, inférieure à un seuil e prédéfini. -[MM]=marge_ module_min(f1,To,Tc) : Cette sous-fonction calcule le minimum de MM(f1) pour des valeurs données de To et Tc, sur toute la plage des valeurs de f1. Ce calcul peut se faire par optimisation avec un algorithme de descente par exemple l'algorithme de la sécante, à condition d'initialiser l'algorithme avec une valeur de f1 proche de l'optimum qui peut être grossièrement connue en fonction de To, ipi et 1P2. L'algorithme correspondant de calcul du minimum de la marge de module pour un jeu de paramètre To, Tc donné (yi1 et La ayant été fixés une fois pour toutes préalablement aux calculs) est donné sur la Figure 13. Au départ on utilise une valeur f1(0) fonction de To, ip1 et 1P2 et qui est proche de l'optimum comme expliqué ci-dessus. Par des étapes itératives dans lesquelles des variations à divers ordres sont considérées, on détermine le minimum de MM(To,Tc) en prenant comme critère d'arrêt le fait que la variation de f1 d'une itération à la suivante devient, en valeur absolue, inférieure à un seuil e prédéfini. Dans la fonction [MR,MM]=marges(To,Tc), une optimisation est faite de façon à déterminer les valeurs To et Tc qui permettent que MR=mrc et MM=mmc. Cette optimisation peut en particulier être faite en utilisant l'algorithme de � Newton projeté � de façon à respecter la contrainte To/5 On comprend que le terme � signal � utilisé dans cette description doit être considéré dans son acception large et que selon l'endroit où il se situe et les moyens matériels mis en oeuvre, il peut correspondre à des valeurs numériques (signal numérique échantillonné) ou analogique (courant moteur par exemple) et que des dispositifs de conversion adaptés sont mis en oeuvre afin de rendre compatible le signal avec les entités avec lesquelles il est en relation. Par exemple, un convertisseur numérique-analogique peut être mis en oeuvre dans le correcteur central pour l'attaque du moteur. De même, un convertisseur analogique-numérique ou fréquence-numérique peut être mis en ouvre en sortie du capteur. En pratique, les signaux qui sont soumis à des calculs sont échantillonnés et sont des signaux numériques. Les exemples donnés sont indicatifs et permettent de comprendre que les procédés utilisés sont applicables à d'autres configurations plus ou moins évoluées notamment pour ce qui concerne la précision du modèle du système ou son degré. On comprend que des moyens équivalents à ceux décrits entrent dans le cadre de l'invention. Ainsi, des données ou signaux peuvent se déduirent les uns des autres par des opérations simples comme par exemple la position évolutive et la vitesse et il est alors possible de suivre l'un ou l'autre paramètre pour en déduire l'autre ou l'un. De même les opérations/fonctions peuvent être décomposées ou réparties d'une manière différente. Par exemple le bloc modificateur peut produire un signal modificateur inversé qui sera alors envoyé à une entrée positive (il y a addition au lieu d'une soustraction contrairement au cas d'une entrée négative ou inversée comme représentée Figure 7) de l'additionneur amont de la fonction 1/So(z-l). REVENDICATIONS1. Procédé de rejet robuste des perturbations périodiques dans une structure de commande de moteur en position, les perturbations périodiques étant assimilées à au moins une oscillation sinusoïdale à une fréquence déterminable proportionnelle à la vitesse angulaire, la structure de commande recevant un signal de consigne en position c(t) mettant en oeuvre une boucle de contrôle de position de l'axe dudit moteur, ladite structure de commande étant du type correcteur RST et étant échantillonnée afin de permettre des calculs numériques dans un calculateur programmable pour produire des résultats destinés à la commande du moteur, le correcteur RST mettant en oeuvre des moyens de calcul correspondant, d'une part, à une fonction 1lSo(z 1) en amont du moteur et produisant un signal de commande u(t) et, d'autre part, une fonction Ro(z-') dans sa boucle qui reçoit en entrée un signal de position 0(t) d'un capteur de position sur l'axe dudit moteur, l'ensemble regroupant le moteur, le capteur de position, les moyens de calculs des fonctions 1/So(z-1) et Ro(z-1) avec bouclage étant appelé correcteur central, caractérisé en ce que : - on ajoute au correcteur central deux moyens de calcul supplémentaires, le premier étant un bloc estimateur E(z-') et le second un bloc modificateur Q(z-'), 20 ainsi qu'une boucle supplémentaire, - on modélise le moteur sous forme d'un système modélisé soumis à une perturbation de charge v(t), ledit système modélisé-recevant en entrée le signal de commande u(t) et produisant en sortie le signal de position 0(t), le modèle du système appliquant sur le signal de commande u(t) une fonction z-'B(z1) pour 25 produire un signal additionné à la perturbation de charge v(t), la résultante étant soumise à une fonction 1/A(z-1) pour produire le signal de position 0(t), - on calcule un signal d'estimation de la perturbation de charge v(t) dans le bloc estimateur E(z"1) recevant en première entrée le signal de commande u(t) et en seconde entrée le signal de position 0(t), dans ledit bloc estimateur le signal de 30 commande u(t) étant soumis à la fonction z-1B(z') pour produire un signal qui est soustrait au résultat de l'application de la fonction A(z-1) sur le signal de position 0(t), la soustraction donnant le signal d'estimation de la perturbation de charge v(t) - on applique le signal d'estimation de la perturbation de charge v(t) au bloc 35 modificateur Q(z-1) pour produire un signal modificateur appliqué à une entrée négative d'un additionneur supplémentaire disposé en amont du moyen de calculde la fonction 1/So(z-l) et formant ainsi la boucle supplémentaire, le bloc modificateur ayant une fonction de transfert variable en fonction de la au moins une fréquence déterminable et qui est calculée en temps réel ou quasi-réel sous la contrainte que la fonction de sensibilité de position Sep de la structure de commande comporte un creux pour la au moins une fréquence déterminable. 2. Procédé selon la revendication 1, caractérisé en ce que la boucle supplémentaire est intérieure à la boucle en Ro(z-l) du correcteur central, l'additionneur supplémentaire étant disposé juste en amont du moyen de calcul de la fonction 1/So(z-l). 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que le correcteur central est calculé préalablement, hors ligne, pour rejeter des perturbations statiques et les fonctions So et Ro sont mémorisées pour être directement utilisées ultérieurement lors des calculs en temps réel. 4. Procédé selon la revendication 3, caractérisé en ce que lors du calcul du correcteur central dans lequel on met en oeuvre une étape de détermination d'un horizon de filtrage To et d'un horizon de commande Tc tels que Tc 5. Procédé selon la revendication 4, caractérisé en ce que dans l'étape de détermination de l'horizon de filtrage To et de l'horizon de commande Tc sous contrainte, on calcule un ensemble de minima des marges de retard et des marges de module, les minima étant fonction des valeurs possibles de la/des fréquences à rejeter et étant également obtenus en faisant varier les valeurs de To et Tc et en ce qu'à partir de ensemble de minima des marges de retard et des marges de module on remonte aux valeurs To et Tc qui assurent le respect de la contrainte que la marge de retard soit supérieure à une marge de retard seuil prédéterminée mrc et la marge de module soit supérieure à une marge de module seuil prédéterminée mmc et, en plus, de préférence, que To/5 6. Procédé selon lune quelconque des revendications précédentes, caractérisé en ce qu'on calcule la fonction de transfert Q(z-l) du bloc modificateur en considérant le rapport de polynômes Q(z-1) = fl(z-1) a(z 1)à paramètres variables et on assimile par équivalence la structure de commande globale à un correcteur RST en R et S de blocs de fonctions polynomiales R(z-1) = Ro(z-1) + A(z-1) Q(z-1) S(z-1) = So(z-1) _ z-1 B(z-1) . Q(z-1) on factorise la fonction polynomiale S(z-1) telle que pré-spécifiée en S(z-1)=Hs(Z- 1).S'(z-1) pour faire apparaître un terme en (1-z-1) correspondant au rejet de perturbations statiques et des termes S'(z-1) et Hsl (z-1), S(z-1) s'exprimant S(z 1)= (1 -z 1).Hs1(z-1).S'(z-1) et en ce que l'on détermine les polynômes Hsl et a par discrétisation d'une fonction de transfert continue dépendant du niveau de rejet Mt souhaité à la fréquence considérée, le polynôme R étant ensuite obtenu par la résolution de l'équation de Bézout suivante S,(z-1).Hs(z-1)+z-1B(z-1). fl(z-1) = So(z-1).a(z-1) 7. Procédé selon la revendication 6, caractérisé en ce que la fonction de transfert continue dépendant du niveau de rejet souhaité à la fréquence considérée dans le cas d'une seule fréquence f1 est : F(s)= s2+2 nnm.2~z.f1 •s+(27c.f1)2 s2 +2 den ..f1 .s+(27c.f1)2 8. Procédé selon la revendication 6, caractérisé en ce que la fonction de transfert continue dépendant du niveau de rejet souhaité à la fréquence considérée 20 dans le cas de deux fréquences f1 et f2 est : s2+2 den1•2,r•f1's+(2,r.f1)2 s2+2 den2•2,r•f2s+(2,r•f2)2 9. Procédé selon la revendication 6, 7 ou 8 caractérisé en ce que le niveau de rejet Mt pour une fréquence f; donnée est calculé par Mt;=20.1og(num;/éden;) et, de préférence on prend den;=0,1. 25 10. Dispositif pour rejet robuste des perturbations périodiques dans une structure de commande de moteur en position, les perturbations périodiques étant assimilées à au moins une oscillation sinusoïdale à une fréquence déterminable proportionnelle à la vitesse angulaire de l'axe, la structure de commande recevant un signal de consigne en position c(t) mettant en oeuvre une boucle de contrôle de 30 position de l'axe dudit moteur, ladite structure de commande étant du type correcteur RST et étant échantillonnée afin de permettre des calculs numériques F(s) = s2 + 2&,m1 •2,r • fi s + (2,r . fi )2 s2 + 2 nnm2 •2,r • f2 s + (2,r • f2 )2dans un calculateur programmable pour produire des résultats destinés à la commande du moteur, le correcteur RST mettant en oeuvre des moyens de calcul correspondant, d'une part, à une fonction 1/So(z-l) en amont du moteur et produisant un signal de commande u(t) et, d'autre part, une fonction Ro(z-l) dans sa boucle qui reçoit en entrée un signal de position 8(t) d'un capteur de position sur l'axe dudit moteur, l'ensemble regroupant le moteur, le capteur de position, les moyens de calculs des fonctions 1/So(z-l) et Ro(z-l) avec bouclage étant appelé correcteur central, caractérisé en ce qu'il comporte, par ajout au correcteur central, des moyens pour l'exécution du procédé de l'une quelconque des revendications précédentes, lesdits moyens étant deux moyens de calcul supplémentaires, le premier étant un bloc estimateur E(z-l) et le second un bloc modificateur Q(z-l), ainsi qu'une boucle supplémentaire, le bloc estimateur E(z-l) permettant le calcul d'un signal d'estimation de la perturbation de charge v(t) , ledit bloc estimateur recevant en première entrée le signal de commande u(t) et en seconde entrée le signal de position 8(t), le bloc modificateur Q(z-l) recevant en entrée le signal d'estimation de la perturbation de charge v(t) pour produire un signal modificateur appliqué à une entrée négative d'un additionneur supplémentaire disposé en amont du moyen de calcul de la fonction 1/So(z-l) et formant ainsi la boucle supplémentaire, le bloc modificateur ayant une fonction de transfert variable en fonction de la au moins une fréquence déterminable et qui est calculée en temps réel ou quasi-réel sous la contrainte que la fonction de sensibilité de position Sep de la structure de commande comporte un creux pour la au moins une fréquence déterminable.