FR 3107367 A1 20210820 FR 2001504 A 20200214 Dispositif comprenant au moins un système électromécanique et une unité électronique de commande comportant un bouclage d’asservissement qui relie l’entrée et la sortie du système électromécanique et qui comporte un correcteur numérique. Le dispositif comprend un module logiciel ayant une entrée et une sortie reliées respectivement à la sortie et à l’entrée du système électromécanique. Le module logiciel est programmé pour, pour chaque fréquence d’un ensemble de fréquences à analyser: émettre un signal sinusoïdal d’excitation sur sa sortie et récupérer un signal de réponse sur son entrée; modéliser par une représentation de Fresnel le signal de réponse et le filtrer autour de la fréquence d’excitation par au moins un filtre de Kalman pour obtenir un vecteur d’état du signal de réponse filtré; obtenir au moins une fonction de transfert du système à partir de la modélisation et du signal d’excitation et du signal d’entrée. FIGURE DE L’ABREGE: Fig. 1 Mesure de fonction de transferts dans un système mécatronique ARRIERE PLAN DE L’INVENTION La présente invention concerne le domaine des systèmes asservis comme des systèmes mécatroniques ou électromécaniques. Dans un tel système, l’entrée et la sortie du système sont liées entre elles par une fonction de transfert. La sortie d’un système comprenant un actionneur pour asservir la sortie du système en tenant compte d’une mesure de cette sortie dépend de l’entrée du système mais également de perturbations internes ou externes, de défauts dans la chaîne de commande, et/ou dans la boucle de rétroaction, et/ou dans l’actionneur… de sorte que la fonction de transfert du système peut être difficile à déterminer. Or, la connaissance de cette fonction de transfert permet de piloter plus précisément le système ainsi bouclé par l’actionneur et le ou les capteurs mesurant la sortie du système. Il est connu de raccorder entre l’entrée et la sortie du système mécatronique ou électromécanique un analyseur de spectre ou un oscilloscope mettant en œuvre des fonctions mathématiques ou de traitement du signal pour en déterminer la fonction de transfert. Ceci impose cependant de pouvoir raccorder un matériel externe en utilisant éventuellement des interfaces telles que des convertisseurs analogiques/numériques ou numériques/analogiques. La détermination de la fonction de transfert nécessite donc d’interrompre le fonctionnement du système pour ouvrir la boucle de rétroaction et conduire une opération de maintenance en faisant intervenir un opérateur. OBJET DE L’INVENTION L’invention a notamment pour but de permettre de déterminer au moins une des fonctions de transfert d’un système électromécanique ou mécatronique avec un impact limité sur son fonctionnement et avec un équipement limité. A cet effet, on prévoit, selon l’invention un dispositif comprenant au moins un système électromécanique ayant une entrée et au moins une sortie, et une unité électronique de commande comportant un bouclage d’asservissement qui relie l’entrée et la sortie du système électromécanique et qui comporte un correcteur numérique. Le dispositif comprend un module logiciel ayant une entrée reliée à la sortie du système électromécanique et une sortie reliée à une entrée du système électromécanique, et le module logiciel est programmé pour, chaque fréquence d’un ensemble de fréquences à analyser: engendrer au moins un signal sinusoïdal d’excitation ; émettre le signal sinusoïdal d’excitation sur sa sortie et récupérer un signal de réponse sur son entrée; modéliser par une représentation de Fresnel le signal de réponse et le filtrer autour de la fréquence d’excitation par au moins un filtre de Kalman pour obtenir un vecteur d’état du signal de réponse filtré; multiplier le vecteur d’état par une matrice de rotation pour obtenir une partie réelle et une partie imaginaire; effectuer une division de nombre complexe pour obtenir au moins une fonction de transfert du système à partir de la modélisation et du signal d’excitation et du signal d’entrée. Le procédé de l’invention permet de déterminer les fonctions de transfert suivantes pour un même système: la fonction de transfert de la boucle fermée avant et après l’injection du signal sinusoïdal; la fonction de transfert de la boucle ouverte la fonction de transfert de la partie electromécanique; la fonction de transfert de la partie correctrice. Ces fonctions de transfert peuvent être déterminées par le procédé de l’invention sans qu’il soit besoin de déboucler l’asservissement et les calculs nécessaires à l’estimation sont relativement simples et ne nécessitent pas de ressources de calcul importantes du fait de l’utilisation: d’une représentation de Fresnel pour modéliser les signaux d’excitation et de mesure dans le système; d’un filtre de Kalman à gain convergé pour filtrer des mesures autour de la fréquence d’excitation; d’un simple produit matricielle pour calculer la fonction de transfert (sous forme complexe) entre l’excitation et un des mesures (ce qui donne la boucle fermée avant ou après la perturbation suivant l’emplacement du point de mesure); différents points de mesures de la boucle (pour une même excitation) et par simples divisions des fonctions transfert complexes entre elles, il est possible d’obtenir: la fonction de transfert de la boucle ouverte, la fonction de transfert de la partie électromécanique, la fonction de transfert de la partie correctrice. L’invention est par exemple utilisable dans un système d’orientation et de stabilisation de la ligne de visée d’un dispositif optronique en pilotant la motorisation d’une plateforme gyro-stabilisée sur laquelle est monté le dispositif optronique. Il est aussi possible d’utiliser l’invention dans d’autres systèmes d’actionnement tels que le système d’actionnement des gouvernes d’un aéronef. D’autres caractéristiques et avantages de l’invention ressortiront à la lecture de la description qui suit d’un mode de réalisation particulier et non limitatif de l’invention. Il sera fait référence aux dessins annexés, parmi lesquels: la figure 1 est une vue schématique en perspective d’un dispositif optronique mettant en œuvre l’invention; la figure 2 est une vue schématique de l’unité de commande de ce dispositif; la figure 3 est une vue schématique de détail de cette unité de commande de ce dispositif. la figure 4 est une représentation de Fresnel du signal de mesure et de la fonction de transfert du système; la figure 5 est une vue schématique en perspective d’un dispositif d’actionnement de gouvernes mettant en œuvre l’invention. DESCRIPTION DETAILLEE DE L’INVENTION En référence en particulier à la figure 1, l’invention est ici décrite en application à un dispositif optronique comprenant un capteur optique 1 monté sur une plateforme gyro-stabilisée 2. La plateforme gyro-stabilisée 2 comprend deux axes motorisés, à savoir: un axe de gisement z permettant de régler l’orientation d’un axe de visée v du capteur optique 1 autour de l’axe de gisement z et un axe de site y perpendiculaire à l’axe de gisement z et à l’axe de visée v pour régler l’orientation de l’axe de visée v autour de l’axe de site y. Chacun des axes z, y est pourvu d’un moteur 3, 4 agencé pour orienter la ligne de visée v du capteur 1 autour de l’axe z, y respectivement et associé chacun à un capteur angulaire 3bis, 4bis comme un codeur ou résolveur. Les moteurs 3, 4 et les capteurs angulaires associés 3bis, 4bis sont reliés à un dispositif de pilotage 10 qui transmet des commandes m1, m2 aux moteurs 3, 4 et qui reçoit des mesures c1, c2 des capteurs angulaires 3bis, 4bis. Le capteur optique 1 est pourvu d’au moins un capteur angulaire inertiel 5 pour détecter une vitesse angulaire c3 de l’axe de visée v par rapport à un premier axe de référence colinéaire à l’axe de site y et une vitesse angulaire c4 de l’axe de visée v par rapport à un deuxième axe de référence perpendiculaire à l’axe de visée v et à l’axe de site y. Le capteur angulaire inertiel est ici un gyromètre à deux axes sensibles comme celui produit sous la marque GSL par la demanderesse. En variante, tout autre dispositif de mesure angulaire est utilisable et par exemple deux capteurs inertiels ayant chacun un axe sensible sont utilisables, comme celui produit sous la marque QUAPASON par la demanderesse. Le capteur inertiel peut être un capteur mécanique à résonateur vibrant, par exemple un résonateur hémisphérique, ou utiliser d’autres technologies comme les fibres optiques (capteur dit FOG). Le capteur inertiel 5 est lui aussi relié au dispositif de pilotage 10 à qui il transmet les signaux de mesure c3, c4. En référence également aux figures 2 et 3, chaque ensemble formé, pour chaque axe, par un moteur 3, 4 et un des capteurs angulaires associés 5, 3bis, 4bis (relatifs et/ou inertiels) forme un système linéaire S ayant une entrée (la ligne d’alimentation fournissant un signal de tension au moteur 3, 4) et une sortie (la ligne de sortie des capteurs angulaires 5, 3bis, 4bis fournissant un signal de mesure). Le dispositif de pilotage 10 comprend une unité électronique de commande 11 comportant, pour chaque système linéaire noté H, une boucle d’asservissement qui relie l’entrée et la sortie du système linéaire par un correcteur numérique 12. L’unité électronique de commande 11 est ainsi agencée pour piloter via le correcteur numérique 12 les moteurs en fonction de commandes provenant d’une interface utilisateur 20 (visible sur la figure 1) et des signaux de mesure provenant soit des codeurs soit du ou des capteurs inertiels. On comprend que chaque signal de commande m en entrée d’un des moteurs 3, 4 provoque un mouvement correspondant du capteur optique 1, mouvement qui est détecté par le capteur angulaire relatif 3bis, 4bis ou inertiel 5 qui est sensible selon l’axe correspondant au moteur 3, 4 et qui émet un signal de mesure c correspondant à la convolution de la fonction de transfert H par l’écart du signal de mesure c avec la consigne, soit: m=K*(consigne-c) et c=H*m. Ce mode de fonctionnement est classique et ne sera pas décrit plus avant ici. Le dispositif de pilotage 10 échange des données avec des modules mémoires 30 et 40. Le module mémoire 30 comporte un tableau des fréquences f p à balayer, de l'amplitude P 0 du signal d'excitation associé, du vecteur des gains de Kalman K f associés. Le dispositif de pilotage 10 comprend un module logiciel 13 ayant une sortie p reliée à une sortie e3 du correcteur numérique 12, au moins une première entrée de mesure e1 reliée en amont d’une entrée du correcteur numérique 12 et deux entrées de mesure e2 et e3 reliées à la sortie du correcteur numérique 12 en aval et en amont de la liaison de ladite sortie du correcteur numérique 12 avec la sortie p du module logiciel 13. Ainsi, la sortie p additionnée à la sortie e3 du correcteur numérique 12 donne la commande e2 qui sera appliquée au moteur via le convertisseur numérique/analogique et le signal m Le module logiciel 13 est programmé pour: former un générateur d’excitation à base d’une matrice d’incrément de rotation et d’une modélisation de Fresnel, le générateur étant agencé pour émettre (vers une entrée p du système S) des signaux sinusoïdaux d’excitation P de différentes fréquences et/ou amplitudes sur la sortie du correcteur en utilisant une somme logicielle e2=e3+P; récupérer sur chacune de ses entrées un signal sinusoïdal de mesure e1 à e6 pour chacun des signaux sinusoïdaux d’excitation; modéliser les signaux de mesures e1 à e6 (figure 3) à partir d’un vecteur de Fresnel et d’une représentation d’état; filtrer les signaux de mesure e1 à e6 par des filtrages passe-bande autour de la fréquence d’excitation au moyen de plusieurs filtres de Kalman mis en parallèle pour obtenir des vecteurs d’état des signaux de mesures filtrés; multiplier les vecteurs d’état des signaux filtrés par une matrice de rotation «T» pour obtenir les parties réelle et imaginaire des fonctions de transfert. Les composantes de la matrice de rotation «T» sont obtenues à partir des composantes du vecteur d’excitation. mémoriser les parties réelle et imaginaire des fonctions de transfert (référence 40 sur la figure 1) et obtenir les fonctions de transfert du bouclage en utilisant des divisions de nombres complexes. Le module logiciel 13 est agencé pour balayer plusieurs fréquences en récupérant pour chacune d’elle l’amplitude de l’excitation et les gains des filtres de Kalman stockés dans une mémoire (30 -figure 1). Le module logiciel 13 est agencé par exemple pour pouvoir estimer plusieurs fonctions de transfert, identifiées FT1 à FT6 sur la figure 3, représentative chacune d’un ensemble de six mesures (par exemple courants, vitesses inertielles, positions relatives) appartenant à un système électromécanique S mis en œuvre dans un même dispositif. Le principe, selon l’invention, de la détermination d’une des fonctions de transfert va maintenant être décrit dans ses trois phases principales: génération d’un signal d’excitation sinusoïdal à la fréquence fp, filtrage passe-bande en utilisant un filtrage de Kalman autour de la fréquence fp, multiplication par une matrice de rotation pour obtenir la partie réelle et imaginaire de la fonction de transfert entre l’excitation et la mesure considérée à la fréquence fp. Le signal d’excitation (perturbateur) est un signal sinusoïdal, par exemple un cosinus d’amplitude P 0 à la fréquence , et peut se modéliser sous la forme d’un vecteur qui est représentatif d’un nombre complexe et qui tourne à une vitesse constante (représentation dite de Fresnel). A l’instant t k , le nombre complexe s’écrit: La perturbation en tension réellement injectée peut s’exprimer comme étant la partie réelle du nombre complexe, soit: La partie imaginaire du nombre complexe est, elle, laissée libre: Ceci s’écrit sous forme matricielle de la façon suivante: Pour une fréquence fixée, la matrice de transition entre deux périodes d’échantillonnage T e correspond à la rotation d’angle où T e est la période d’échantillonnage du calculateur (T e = t k+1 -t k ). La matrice de transition s’écrit: La relation matricielle entre deux instants consécutifs est donc: On remarquera pour la suite que la matrice de rotation à l’instant k est le produit de k incréments de transition: On remarque que la matrice de rotation inverse est simplement la transposée. Elle se calcule aussi avec les composantes du vecteur d’excitation: La notation désigne la «n ième » composante du vecteur représentatif de la perturbation. En considérant la perturbation P k comme la première composante du vecteur d’état, on obtient le système dynamique discret suivant: dans lequel B=D=0 et C=[1 0]. Suivant la figure 2, on applique le signal d’excitation (perturbation) sur la sortie du correcteur numérique en utilisant une somme logicielle e 2 = e 3 + P. La forme matricielle du système dynamique sera adaptée pour faire du filtrage de Kalman. On prend P 0 comme condition initiale d’amplitude du signal d’excitation (perturbateur) et une phase nulle de sorte que le vecteur initial est: Concernant la réponse fréquentielle du système électromécanique et produit matriciel, la formulation sous forme de système dynamique linéaire permet d’assimiler le système physique à un filtre linéaire qui reçoit en entrée la perturbation P k pour donner en sortie la mesure M k . Avec la transformée en z et en notant la fonction de transfert que l’on cherche à estimer, on obtient: Avec l’opérateur retard , on a: Pour un système linéaire excité avec un cosinus pur à la fréquence , la mesure est un cosinus d’amplitude M P et de phase ψ P et prend la forme matricielle suivante: La mesure directement accessible est alors que l’on cherche : Qui s’écrit sous forme matricielle: La mesure «matricielle» et la fonction de transfert sont unies par les relations trigonométriques suivantes (voir la représentation de Fresnel de la figure 4): Soit, en l’inversant: Et: Les composantes réelles et imaginaires de la fonction de transfert sont obtenues par un simple produit matriciel du signal de mesure idéal avec une simple matrice de rotation dont les composantes sont déjà calculées par la génération du signal d’excitation: avec i l’imaginaire pur tel que i 2 = -1. On aborde désormais le filtrage passe-bande par filtrage de Kalman. Pour estimer la partie réelle et la partie imaginaire de la mesure à partir de la mesure réelle M k , on utilise un filtre de Kalman à gain convergé qui fait office de filtre passe-bande sélectif autour de la fréquence d’excitation: Dans lequel: w k est le bruit du modèle scalaire, v k est le bruit de mesure scalaire, ∆A la matrice de transition, C=[1 0] est la matrice d’observation du générateur de perturbation, B T =[1 1] est la matrice d’injection des bruits d’état. La variance du bruit de mesure est . La matrice de variance du bruit de modèle est: Le filtre peut se représenter sous la forme de la représentation d’état suivante: La matrice K f est la matrice de gain convergé du filtre de Kalman telle que: Dans lesquelles Z est la solution de l’équation de Riccati. En pratique, pour déterminer en temps réel la fonction de transfert du système linéaire, on connaît la fréquence f P de chaque signal de perturbation, l’amplitude P 0 , la durée du filtrage T f et le gain du filtre de Kalman soit: Pour effectuer le balayage de toutes les fréquences dans un intervalle donné, on met en œuvre l’algorithme suivant: On notera qu’ici les équations du filtre de Kalman ont été mises sous la forme d’un filtre récursif avec un seul vecteur d’état comme dans l’équation 2.31 du document D. Alazard, Introduction au filtrage de Kalman. Il est cependant possible de les écrire aussi sous la forme prédiction/recalage avec deux vecteurs d’état (voir les équations 2.23 et 2.25 du même document) comme dans la suite de la présente description en relation avec la composante continue. On obtient donc ainsi une estimation d’au moins trois fonctions de transfert du système bouclé, entre le signal d’excitation et au moins trois mesures, telles que sans avoir eu besoin de l’intervention d’un opérateur pour séparer la boucle d’asservissement de l’actionneur. Les différentes fonctions de transfert du bouclage sont ensuite déterminées de la manière suivante. Avec les notations de la figure 2, et en utilisant des divisions de nombres complexes, on réaliser les opérations suivantes: a/ Pour la fréquence courante, mémoriser pour la fréquence courante les parties réelles d'au moins les trois boucles fermée suivantes: b/ Calculer la boucle ouverteavec : c/ Calculer la fonction de transfert de la partieélectromécanique : On aborde ensuite un mode de réalisation préféré dans lequel on procède à une reformulation et une extension pour prendre en compte la composante continue des mesures. Selon le mode de réalisation préféré, on injecte, non plus un cosinus, mais un sinus d’excitation ayant une fréquence f P prédéterminée (en Hz). Lors de l’injection du signal d’excitation de type sinus dans un système linéaire, le système répond par une raie sinus amplifiée d’un gain et déphasée d’une phase La représentation d’état de la raie sinusoïdale de la mesure sans composante continue est donc de dimension 2 et donnée par: Le vecteur d’état représente l’abscisse (cosinus) et l’ordonnée (sinus) d’un vecteur tournant sur un cercle de rayon à la fréquence f p depuis la position initiale déterminée par l’angle . Ainsi, on a les matrices de la représentation d’état du vecteur de Fresnel: où, de manière classique en représentation d’état, est la matrice de transition, B la matrice de commande, C la matrice d’observation et D la matrice d’action directe. Cependant, en pratique, il est probable que la réponse du système comprendra également une composante continue correspondant au point de fonctionnement autour duquel le système est excité. Cette composante continue constitue un biais dont la représentation d’état discrète est donnée ci-dessous: En ajoutant cette représentation à la représentation d’état de la raie sinus, on obtient la représentation d’état d’un sinus avec composante continue de dimension 3 pour les différentes mesures, à savoir: Le générateur d’excitation n’a pas besoin lui d’une composante continue. Pour concevoir le générateur d’excitation sinusoïdale, il suffit de simuler le modèle décrit par la représentation d’état de la raie sinus. Le choix de l’amplitude et de la phase est effectué en imposant les conditions initiales du modèle. On choisit une phase nulle pour le sinus à générer de manière à démarrer progressivement l’excitation du système (sin(t=0)=0). On obtient comme suit un générateur de sinus d’amplitude P 0 : de condition initiale avec la dynamique du signal d'excitation: et les matrices de la représentation d'état: Chaque mesure issue d’un capteur est une raie sinus présentant une composante continue, dont on estimera l’état du vecteur de Fresnel grâce à un filtre de Kalman formant un estimateur correspondant à un filtre passe bande autour de la fréquence fp pour les deux premières composantes du vecteur d’état et un filtre passe bas pour la partie continue, en se basant sur le modèle (1), (2). Cet état une fois estimé, on est capable d’estimer la fonction de transfert correspondante entre l’excitation (entrée) et la mesure (sortie). En sortie de l’estimateur de fonction transfert, on aura directement la partie réelle et la partie imaginaire de la fonction de transfert à la fréquence f p . Le modèle de l’estimateur de Kalman augmenté prend la forme générale suivante: Avec: comme précédemment, U k représente les commandes, v k le bruit de mesure, w k le bruit du modèle, et la matrice M choisie est la suivante Comme les matrices sont nulles, il est possible de simplifier le modèle: On note V[1x1] et W[3x3] les matrices de covariance du bruit de mesure v k et du bruit de modèle w k qui sont supposés gaussiens dans le cas d’un estimateur de Kalman linéaire. Le régime permanent de l’estimateur est déterminé uniquement par le choix des matrices de covariance V et W. Vu que l’estimateur est implémenté en gain convergé, on peut fixer arbitrairement V=1 sans perdre de généralité. On choisit W, la matrice de covariance du bruit du modèle, comme suit: Ainsi, le choix du couple de variances permet de définir complètement l’estimateur de Kalman en régime permanent. Une représentation d’état de l’estimateur est notée: On obtient alors: Dans cette représentation, Z est l’unique solution stabilisante de l’équation algébrique de Riccati: On notera que cette équation présente une singularité dans le cas où la fréquence f p est égale à la moitié de la fréquence d’échantillonnage F e =1/T e . Dans le cas général, cette équation peut être résolue par exemple avec la fonction «dare» du langage «Matlab» qui permet également de déterminer les valeurs propres de A k permettant de déterminer la constante de temps du filtre passe-bande pour les deux premières composantes du vecteur d’état et du filtre passe-bas de la composante continue. Comme indiqué ci-dessus, le réglage de l’estimateur (c’est-à-dire la largeur de bande autour de f p et de la fréquence de coupure pour le filtre passe-bas) est obtenu par le choix du couple de variances qui permet de définir complètement l’estimateur de Kalman en régime permanent. Considérant la variance fixée, une diminution de la variance induit une diminution de la largeur du filtrage passe-bande et donc une meilleure estimation de la raie du fait d’un meilleur filtrage du bruit et des éventuelles non-linéarités du système. Si la constante de temps du filtre est égale à celle de la raie, cette diminution de la variance s’accompagne d’une augmentation de la constante de temps du filtre. Il faut donc rechercher un compromis entre la qualité de l’estimation (étroitesse du filtre) et le temps de convergence de l’estimateur. Considérant la variance fixée, une diminution de la variance induit une dégradation du filtrage des basses fréquences pour les composantes cosinus et sinus. En outre la fréquence de coupure de la réponse correspondant à la composante continue (biais) diminue lorsque la variance diminue (cela revient à augmenter le bruit de mesure relativement au bruit du modèle), d’où un filtrage plus conservatif. A titre d’exemple, pour une fréquence d’échantillonnage F e donnée, des fréquences minimum et maximum f min et f max du signal d’excitation, et un nombre N points de fréquences à identifier dans cet intervalle de fréquences, on propose de régler l’estimateur de la manière suivante: détermination de la variance pour obtenir la largeur souhaitée du filtrage passe-bande. Ceci est effectué plutôt en moyenne/haute fréquence dans l’intervalle de fréquences de sorte que le filtre ne soit pas perturbé par la composante continue du modèle. On s’assurera que le filtre est bien symétrique autour de la fréquence sélectionnée; détermination de la variance sur la fréquence la plus basse à identifier (f min ) de sorte que la sélectivité du filtre reste acceptable à cette fréquence; calculer la durée totale de l’identification induite par ce paramétrage. Si la durée n’est pas satisfaisante, augmenter la variance ou réduire le nombre de fréquences à identifier. Le type de balayage a un fort impact sur la durée d’identification: pour réduire celle-ci, il est préférable de choisir un balayage linéaire plutôt que logarithmique. Dans le cas présent, on choisit pour une fréquence d’échantillonnage F e de 3 kHz. Ce réglage est adapté par exemple pour le pointage d’un dispositif optronique. Pour une autre application, les valeurs précitées peuvent servir de points de départ pour converger vers un réglage plus adapté à ladite application. En particulier, une telle application pourrait nécessiter une fréquence d’échantillonnage différente. Or, on sait que pour un réglage donné du couple , la largeur du filtre passe-bande augmente avec la fréquence d’échantillonnage: l’identification est donc plus rapide mais de qualité moindre. Ce problème peut être résolu en recalculant les valeurs du couples garantissant une constante de temps quasi-identique et une réponse fréquentielle équivalente lors du passage d’une fréquence d’échantillonnage à une autre. En variante, il est possible d’alléger les calculs en utilisant l’estimateur non plus sous la forme d’une représentation d’état comme précédemment (c’est-à-dire sous la forme d’un filtre récursif à un seul vecteur d’état) mais sous la forme prédiction/recalage avec deux vecteurs d’état en notant l’état prédit et l’état recalé (sortie de l’estimateur), les équations de récurrence du filtre à l’instant k+1 sont: Le gain de Kalman K f est le même que celui défini ci-dessus. Cette implémentation est mathématiquement équivalente à la précédente mais moins coûteuse en temps de calcul. L’estimateur de fonction de transfert calcule les parties réelle et imaginaire de la fonction de transfert entre le signal d’excitation sinusoïdal et la mesure obtenue, de l’une des trois boucles fermées suivantes : On cherche donc à estimer (i étant l’imaginaire pur tel que i 2 =-1): On note l’état du générateur d’excitation à l’instant k et l’état du processus estimé par le filtre de Kalman à partir de la mesure. On obtient: On obtient le couple à partir de l’état du processus par une matrice T construite à partir des termes et P 0 : Comme précédemment, est une matrice de rotation. On a donc: De là découle une estimation de la partie réelle R p et de la partie imaginaire I p telle que: Cette estimation ne doit être réalisée qu’une fois que l’estimateur de Kalman a suffisamment convergé. Pour s’en assurer, on laisse s’écouler une durée égale à plusieurs fois la constante de temps caractéristique de l’estimateur avant d’estimer la fonction de transfert. En prenant une durée égale à six fois la constante de temps caractéristique , on assure une convergence de l’erreur d’estimation de 0,25%. On notera que, pour s’assurer que les biais du système mécanique (par exemple liés à la présence d’un balourd ou d’un organe de rappel dans le cas d’une transmission de mouvement) ne viennent pas perturber la détermination de la fonction de transfert, il faut s’assurer que l’amplitude du signal sinusoïdal soit suffisante pour qu’on puisse négliger les effets qui ne sont pas engendrés par le signal sinusoïdal. En référence à la figure 5, pour un dispositif de gouverne d’avion, la ligne de visée V est remplacé par le plan de l’aileron. L’aileron 1 est relié à l’ensemble moteur 3 et codeur 3bis par un arbre de transmission 2 en liaison pivot avec le reste de l’avion. Bien entendu, l’invention n’est pas limitée au mode de réalisation décrit mais englobe toute variante entrant dans le champ de l’invention telle que définie par les revendications. En particulier, l’invention est applicable à tout type de système électromécanique modélisable par une fonction de transfert linéaire et n’est donc pas limitée à une application à un dispositif de pointage et de stabilisation d’un capteur optique ou à une gouverne d’avion. L’invention est notamment utilisable avec tout dispositif comprenant un actionneur commandé via un asservissement et des mesures accessibles par un calculateur. D’autres paramètres de modélisation peuvent être choisis en fonction de l’application envisagée. La fréquence d’échantillonnage, les variances, le nombre de fréquences à estimer… Par signal sinusoïdal, on entend un signal sinus comme un signal cosinus. Le module logiciel peut recevoir en entrée soit la mesure de la variable de sortie du système linéaire et/ou la différence entre la mesure de la variable de sortie du système linéaire d’une part et la consigne de commande d’autre part. Bien que dans les modes de réalisation décrits la matrice de rotation soit définie en fonction des composantes du vecteur représentatif du signal sinusoïdal d’excitation, ce qui est avantageux en terme de simplification des calculs, ceci n’est pas obligatoire. Par ailleurs, on notera que, dans le mode de réalisation décrit, on prend en compte les trois mesures e1, e2, e3: e2 et e3 pour la «boucle ouverte» (K*H) soit une en amont du sommateur (e3) et une en aval du sommateur (e2), e1 et e2 pour la fonction de transfert du système électromécanique (H) Par contre, si l’on veut mesurer à la fois les deux vitesses inertielles, la position des deux codeurs et la mesures de deux courants moteurs, on peut utiliser e4, e5, e6 en plus de e1. La matrice T peut avoir différentes formes : Dispositif comprenant au moins un système électromécanique ayant une entrée et au moins une sortie, et une unité électronique de commande comportant un bouclage d’asservissement qui relie l’entrée et la sortie du système électromécanique et qui comporte un correcteur numérique, caractérisé en ce que le dispositif comprend un module logiciel ayant une entrée reliée à la sortie du système électromécanique et une sortie reliée à une entrée du système électromécanique, et en ce que le module logiciel est programmé pour, pour chaque fréquence d’un ensemble de fréquences à analyser: -engendrer au moins un signal sinusoïdal d’excitation (P); -émettre le signal sinusoïdal d’excitation sur sa sortie et récupérer un signal de réponse sur son entrée; -modéliser par une représentation de Fresnel le signal de réponse et le filtrer autour de la fréquence d’excitation par au moins un filtre de Kalman pour obtenir un vecteur d’état du signal de réponse filtré; -multiplier le vecteur d’état par une matrice de rotation pour obtenir une partie réelle et une partie imaginaire; -effectuer une division de nombre complexe pour obtenir au moins une fonction de transfert du système à partir de la modélisation et du signal d’excitation et du signal d’entrée. Dispositif selon la revendication 1, dans lequel le générateur de signaux est basé sur une représentation d’un vecteur de Fresnel en rotation suivant le modèle suivant: -condition initiale: -la dynamique du signal d'excitation : -signal à appliquer en sortie du correcteur: -avec les matrices de la représentation d'état: dans lesquelles . Dispositif selon la revendication 1 ou 2, dans lequel la dynamique des signaux d’excitation et la dynamique des signaux de mesure sont définies par les équations suivantes: dans lesquelles comporte au moins comme composantes le cosinus et le sinus du vecteur d’état, G p est le gain, f p est la fréquence d’excitation, est la phase, dt est la période d’échantillonnage, représente un bruit du modèle et v k un bruit de mesure, b 0 est une composante continue du système, sont des matrices d’état, et la matrice M est telle que Dispositif selon la revendication 3, dans lequel les matrices sont déterminées de la manière suivante: Avec la période d'échantillonage du calculateur et la fréquence courante du balayage. Dispositif selon la revendication 3, dans lequel les matrices sont déterminées de la manière suivante: Avec la période d'échantillonnage du calculateur et la fréquence courante du balayage. Dispositif selon la revendication 3, dans lequel les matrices ont déterminées de la manière suivante: Avec la période d'échantillonnage du calculateur et la fréquence courante du balayage. Dispositif selon la revendication 3, dans lequel les matrices sont déterminées de la manière suivante: Avec la période d'échantillonnage du calculateur et la fréquence courante du balayage. Dispositif selon la revendication 3, dans lequel le filtre de Kalman est exploité sous la forme simplifiée suivante: Dispositif selon l’une quelconque des revendications 2 à 8, dans lequel la matrice de covariance du bruit de mesure est choisie égale à 1 et la matrice de covariance W du bruit de modèle est choisie comme suit: Avec la variance des biais et la variance du signal en entrée du module logiciel. Dispositif selon la revendication 9, dans lequel la variance est fixée 10 -2 et la variance est fixée 10 -4 pour une fréquence d’échantillonnage égale à 3 kHz. Dispositif selon l’une quelconque des revendications 2 à 10, dans lequel la matrice de rotation est définie en fonction des composantes du vecteur représentatif du signal sinusoïdal d’excitation. Dispositif selon la revendication 11, dans lequel la matrice de rotation a la forme suivante: