12.07.2015 Views

Commande boucle fermée multivariable pour le vol en ... - ISAE

Commande boucle fermée multivariable pour le vol en ... - ISAE

Commande boucle fermée multivariable pour le vol en ... - ISAE

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

250 5. MÉTHODOLOGIE POUR PILOTAGE EN ATTITUDE/TRANSLATIONLes Figs. 5.37 et 5.38 montr<strong>en</strong>t <strong>le</strong>s résultats d’une simulation utilisant la stratégie 2b. Ces résultatsne sont pas fondam<strong>en</strong>ta<strong>le</strong>m<strong>en</strong>t différ<strong>en</strong>ts des résultats obt<strong>en</strong>us avec la stratégie 2a. Cep<strong>en</strong>dant, il estbi<strong>en</strong> visib<strong>le</strong> que <strong>le</strong> temps passé dans <strong>le</strong> mode 2 est beaucoup plus é<strong>le</strong>vé <strong>en</strong> utilisant la stratégie 2b. Leschiffres exactes sont indiquées dans <strong>le</strong> Tab. 5.23 (page 254).états et états estimés3210−1−20 200 400 600 800 10004mode3210 200 400 600 800 1000temps [s]Figure 5.37 – Simulation dans laquel<strong>le</strong> la commutation est effectuée selon la stratégie 2b, c’est-à-diredès que dès que la fonction ˜Ξ T P ˜Ξ devi<strong>en</strong>t plus petite que 1. La partie supérieure trace l’é<strong>vol</strong>utiondes états du système (lignes continues) et du correcteur (lignes tiretées, à peine visib<strong>le</strong>s). La partieinférieure indique <strong>le</strong> mode actif <strong>en</strong> fonction du temps. Les lignes vertica<strong>le</strong>s noires sont <strong>le</strong>s instants decommutation.La Fig. 5.38 montre la fonction ˜Ξ T P ˜Ξ <strong>pour</strong> <strong>le</strong>s différ<strong>en</strong>ts modes parcourus et la compare avec lafonction de Lyapunov X T P X utilisée dans la stratégie 2a. D’abord, il est bi<strong>en</strong> visib<strong>le</strong> que la fonctionde Lyapunov est toujours inférieure. Ceci s’explique par <strong>le</strong> fait que la surface ˜Ξ T P ˜Ξ = 1 se trouvetoujours à l’intérieur de la surface X T P X = 1, cf. Fig. 5.36.L’augm<strong>en</strong>tation du temps passé dans <strong>le</strong> mode 2 s’explique par <strong>le</strong> fait que la fonction de Lyapunovdevi<strong>en</strong>t plus petite que 1 rapidem<strong>en</strong>t, tandis que ˜Ξ T P ˜Ξ arrête de décroître et recomm<strong>en</strong>ce à s’accroître<strong>pour</strong> dev<strong>en</strong>ir plus petit que 1 quelque temps après. Ce phénomène est lié au fait que la fonction deLyapunov n’est pas une fonction de Lyapunov du mode courant, mais du mode suivant.Stratégie itérative et exacteNous avons déjà m<strong>en</strong>tionné que <strong>le</strong> critère ˜Ξ T P ˜Ξ ≤ 1 est seu<strong>le</strong>m<strong>en</strong>t suffisant, mais pas nécessaire.Le critère nécessaire et suffisant est si la distance d <strong>en</strong>tre <strong>le</strong> point Ξ = M T 3 M T 2 M T 1 X et <strong>le</strong> point <strong>le</strong> plusproche de l’hyper-ellipsoïde ξ T Pξ = 1 est supérieure à 1.Le point <strong>le</strong> plus proche d’un hyper-ellipsoïde n’est pas trivial à calcu<strong>le</strong>r. En effet, nous n’avons pastrouvé de méthode directe dans la littérature. Cep<strong>en</strong>dant, nous avons pu développer deux algorithmesitératifs capab<strong>le</strong>s de déterminer si la distance d est supérieure ou inférieure à 1. Les détails sont donnésdans l’Annexe I. Notamm<strong>en</strong>t, nous avons montré que ces algorithmes converg<strong>en</strong>t et qu’ils n’ont besoin<strong>Commande</strong> <strong>bouc<strong>le</strong></strong> <strong>fermée</strong> <strong>multivariab<strong>le</strong></strong> <strong>pour</strong> <strong>le</strong> <strong>vol</strong> <strong>en</strong> formation de vaisseaux spatiaux

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!