Résolution numérique de problèmes de contrôle optimal ... - Enseeiht
Résolution numérique de problèmes de contrôle optimal ... - Enseeiht
Résolution numérique de problèmes de contrôle optimal ... - Enseeiht
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Résolution numérique <strong>de</strong><br />
problèmes <strong>de</strong> contrôle <strong>optimal</strong> via<br />
la condition nécessaire, application<br />
au problème <strong>de</strong> transfert d’orbite à<br />
faible poussée<br />
Présentation TIPE 2003: contrôle <strong>optimal</strong><br />
14 janvier 2003 Lycée Fermat Toulouse<br />
J. Gergaud, J. Noailles<br />
gergaud@enseeiht.fr, jnoaille@enseeiht.fr<br />
ENSEEIHT-IRIT, UMR CNRS 5505<br />
Toulouse, FRANCE<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.1/40
Plan<br />
Principe du Maximum <strong>de</strong> Pontriaguine<br />
Exemples simples<br />
Application au problème <strong>de</strong> transfert d’orbite à<br />
poussée faible en temps minimum<br />
Application au problème <strong>de</strong> maximisation <strong>de</strong> la<br />
masse finale<br />
Les métho<strong>de</strong>s directes<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.2/40
¢<br />
¡<br />
<br />
£<br />
<br />
<br />
'<br />
¥¥£<br />
%%<br />
<br />
¦ §¨§¨§¨§¨§¨§¨§¨§¨© §§§¨§§§¨§§<br />
¥ ¤<br />
£<br />
)<br />
<br />
¡<br />
2<br />
<br />
*<br />
¢<br />
2<br />
<br />
*<br />
111<br />
<br />
<br />
Problème <strong>de</strong> contrôle <strong>optimal</strong><br />
fixés<br />
et<br />
¥"<br />
¥ <br />
<br />
¥£<br />
<br />
! £ <br />
£<br />
¢ ¥<br />
£<br />
¢<br />
<br />
¥<br />
<br />
<br />
£<br />
& ¡<br />
¢<br />
dans<br />
£ <br />
<br />
£ ¥ <br />
¥<br />
$<br />
#<br />
£<br />
¥ ( £<br />
/0 <br />
111<br />
(<br />
, ¡.-<br />
£ ¡ ¥ <br />
<br />
+*<br />
¢ -<br />
/0 <br />
(<br />
,<br />
£<br />
¢ ¥ <br />
<br />
+*<br />
où<br />
;<br />
43<br />
¥ ( £<br />
;<br />
53<br />
¥ ( £<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.3/40
3<br />
3¤£<br />
3<br />
§<br />
¦<br />
¥¥<br />
¥<br />
!<br />
¦<br />
¨<br />
0<br />
3<br />
3¤£<br />
3<br />
§<br />
¦<br />
¥¥<br />
¥<br />
!<br />
<br />
§<br />
¦<br />
¨<br />
%<br />
<br />
<br />
©£%<br />
$ £ <br />
:<br />
<br />
<br />
©£%<br />
$ £ <br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.4/40<br />
Hamiltonien<br />
43<br />
¥ £<br />
3<br />
4 £<br />
) £<br />
¢¡<br />
<br />
¥<br />
£ <br />
§<br />
%¡<br />
%¡<br />
%<br />
£ <br />
<br />
<br />
¥©1désigne le produit scalaire.<br />
où<br />
£ 1<br />
la fonction ne dépendra plus <strong>de</strong><br />
%¡<br />
Remarque . Si<br />
%¡<br />
43<br />
4 £<br />
) £<br />
¢¡<br />
¥<br />
£ <br />
<br />
£
$<br />
!<br />
$<br />
!<br />
¢¡<br />
$ <br />
!<br />
<br />
¡<br />
£<br />
¤<br />
<br />
)<br />
Principe du Maximum<br />
Hypothèses simplifiées<br />
(i)<br />
et<br />
sont continues<br />
(ii) Les dérivées partielles<br />
et<br />
existent et sont continues<br />
(iii)<br />
est <strong>de</strong> classe<br />
(iv)<br />
est compact<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.5/40
¤¤<br />
¢<br />
<br />
§<br />
£¢<br />
¤<br />
¤¤<br />
¡£¢ ¤£¥<br />
¦£¢<br />
¤¤<br />
<br />
¦<br />
¡<br />
<br />
<br />
<br />
¤¤<br />
¦<br />
¡<br />
¤<br />
<br />
<br />
¦<br />
¡<br />
¤¤<br />
¡<br />
¤¤<br />
<br />
¦<br />
¡<br />
<br />
<br />
<br />
¤¤<br />
<br />
¦<br />
¡<br />
<br />
¨&<br />
<br />
<br />
¤<br />
¤<br />
¡<br />
*<br />
<br />
<br />
¥<br />
non<br />
¥¢¢¢¤£¥<br />
§<br />
§¨ ¥<br />
§ <br />
¦ ¥<br />
¤£¥<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.6/40<br />
et<br />
§ ¢<br />
!<br />
¥<br />
¦<br />
"<br />
#<br />
§¨ ¥<br />
§ <br />
%"<br />
& '<br />
¤¤<br />
¤£¥<br />
©<br />
<br />
<br />
('<br />
PMP<br />
Théorème Si<br />
est solution alors il existe<br />
simultanément nuls tels que l’on ait<br />
§¨<br />
(i) l’équation adjointe<br />
§¨ ¥<br />
§ <br />
¤£¥<br />
¤£¥<br />
<br />
¥<br />
§ <br />
¤<br />
<br />
¤£¥<br />
¥<br />
<br />
<br />
§¨<br />
<br />
¤¤§<br />
¤£¥<br />
<br />
¥<br />
(ii) la minimisation <strong>de</strong> l’Hamiltonien<br />
§¨ ¥<br />
§ <br />
¤ ¥<br />
¤£¥<br />
¥<br />
<br />
<br />
¤<br />
<br />
¥<br />
<br />
¤<br />
¤ ¥<br />
(iii) les conditions <strong>de</strong> transversalité<br />
%"<br />
§$<br />
¨<br />
<br />
'<br />
)<br />
('<br />
§$
¦ <br />
¤§<br />
¥ ¦ ¨<br />
!<br />
¦<br />
¤<br />
¨<br />
¦ <br />
¤<br />
¤<br />
<br />
¤<br />
<br />
¨<br />
<br />
¡<br />
<br />
¤ <br />
¡<br />
'<br />
¡<br />
§<br />
¤<br />
<br />
<br />
¤<br />
¨<br />
¨<br />
¦<br />
<br />
¤¤<br />
¡ <br />
¡<br />
§<br />
<br />
©<br />
<br />
©<br />
<br />
§<br />
<br />
©¤<br />
<br />
©<br />
§<br />
¦<br />
¤¤<br />
<br />
¦<br />
¤ § <br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.7/40<br />
<br />
<br />
<br />
¢<br />
<br />
<br />
¤<br />
<br />
<br />
¡ <br />
§ <br />
¤<br />
<br />
Exemple 1<br />
¡<br />
¢¡£¡¤¡£¡£¡¤¡¤¡<br />
© ¦ <br />
©¤<br />
¡£¡¤¡£¡¤¡£¡£¡ <br />
Equation adjointe<br />
Minimisation <strong>de</strong> l’Hamiltonien<br />
¥<br />
!<br />
<br />
<br />
¡<br />
¥<br />
¤£¥<br />
<br />
¤<br />
si<br />
<br />
¤<br />
<br />
<br />
sinon<br />
*<br />
§ <br />
<br />
¤<br />
<br />
<br />
Conditions <strong>de</strong> transversalité<br />
<br />
¤
¥£<br />
<br />
<br />
¦<br />
¥¥£ ©<br />
£ <br />
¥£<br />
%<br />
¥ <br />
¤<br />
¦ §§¨§§§¨© §§§§¨§<br />
¥ ¡<br />
£<br />
¤<br />
<br />
¦¤1<br />
¥ <br />
<br />
2<br />
1<br />
0<br />
p<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.8/40<br />
Problème aux <strong>de</strong>ux bouts<br />
¥ £<br />
#<br />
£ ¥<br />
¨§<br />
¦ §§¨© §§<br />
#% £ ¥<br />
£¤<br />
¢¡<br />
£<br />
£¤<br />
¥<br />
£¥<br />
¦¤1<br />
£¥<br />
¥<br />
1.5<br />
0.5<br />
u<br />
−0.5<br />
−1<br />
−1.5<br />
−2<br />
−3 −2 −1 0 1 2 3<br />
#© £ ¥<br />
©<br />
¤<br />
©<br />
¤
3<br />
3<br />
§<br />
¡<br />
¦<br />
¥<br />
§<br />
¡<br />
¢<br />
¡<br />
¨<br />
¦<br />
¥<br />
¡<br />
©<br />
¥£<br />
<br />
¥£<br />
¢<br />
¡¢<br />
0<br />
z<br />
Fonction <strong>de</strong> tir<br />
¦<br />
¦¤1<br />
¤ <br />
£¥<br />
©<br />
¤<br />
¡<br />
¥ <br />
£<br />
est la solution du système à valeur initiale suivant.<br />
¢ ¤ 1<br />
£<br />
où<br />
¥ £<br />
¦ ©<br />
¤<br />
#©<br />
¤<br />
£ ¥<br />
© £<br />
£ ¥<br />
#©¤£<br />
¦ §¨§§§¨§© §§¨§§§¨<br />
£,<br />
£ ¡ ¥<br />
¥ <br />
¤<br />
£¤<br />
©<br />
¤<br />
¥ <br />
£¤<br />
©¥£<br />
0.4<br />
0.2<br />
−0.2<br />
−0.4<br />
S(z)<br />
−0.6<br />
−0.8<br />
−1<br />
−1.2<br />
−1.4<br />
−2 −1.5 −1 −0.5 0 0.5 1 1.5 2<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.9/40
¤<br />
¨<br />
¦ <br />
¤<br />
¤<br />
<br />
¤<br />
<br />
<br />
<br />
¤ <br />
¡<br />
<br />
<br />
¡<br />
<br />
<br />
<br />
¤<br />
¤¦¥<br />
¤<br />
<br />
¡<br />
¥<br />
<br />
<br />
¥ ¦ ¨<br />
!<br />
¤<br />
¥<br />
¤<br />
¨<br />
¦ <br />
¤<br />
¤<br />
<br />
¤<br />
¨<br />
¦ <br />
¤<br />
¤<br />
<br />
© ¦ <br />
©¤<br />
¤<br />
<br />
¤<br />
<br />
<br />
2<br />
<br />
¤ <br />
¡<br />
¢&<br />
£ ¤<br />
<br />
<br />
¤ <br />
¡<br />
<br />
1<br />
<br />
¡<br />
0<br />
§<br />
p<br />
Résoudre<br />
¢<br />
¥<br />
<br />
où est<br />
la solution <strong>de</strong><br />
¡ <br />
¥<br />
<br />
¡ <br />
¡<br />
¢¡¤¡¤¡£¡£¡<br />
<br />
¤<br />
<br />
§<br />
§<br />
¤<br />
¡ ¡¤¡¤¡ ¡ <br />
¥<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.10/40<br />
<br />
¡ <br />
¡<br />
Exemple 2<br />
¢¡<br />
£ ¤<br />
¤©§<br />
<br />
¡ <br />
¡£¡¤¡£¡¤¡ <br />
<br />
¤<br />
<br />
§<br />
§<br />
<br />
¤<br />
¡ ¡¤¡¤¡£¡ <br />
¢<br />
où<br />
1.5<br />
0.5<br />
−0.5<br />
−1<br />
−1.5<br />
−2<br />
−3 −2 −1 0 1 2 3<br />
u<br />
© ¦<br />
¡<br />
¢<br />
<br />
¤<br />
¢¡¤¡£¡£¡¤¡£¡¤¡ <br />
¡ ¡¤¡¤¡£¡£¡¤¡
Fonction <strong>de</strong> tir<br />
0.4<br />
0.2<br />
0<br />
−0.2<br />
−0.4<br />
S(z)<br />
−0.6<br />
−0.8<br />
−1<br />
−1.2<br />
−1.4<br />
−2 −1.5 −1 −0.5 0 0.5 1 1.5 2<br />
z<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.11/40
¤<br />
¤<br />
<br />
¤<br />
<br />
<br />
<br />
¤<br />
<br />
'<br />
¤<br />
<br />
¡<br />
¡<br />
¤¦¥<br />
¤<br />
<br />
<br />
¥ ¦ ¨<br />
!<br />
¥<br />
<br />
¡<br />
<br />
¤<br />
<br />
¤<br />
¡<br />
¡<br />
¥<br />
¤<br />
¤<br />
<br />
¨<br />
¤<br />
¡<br />
¤<br />
<br />
'<br />
¤<br />
¤<br />
<br />
¡<br />
<br />
©¤<br />
<br />
¤<br />
<br />
¤<br />
<br />
<br />
¤ <br />
<br />
<br />
¤<br />
<br />
2<br />
¢&<br />
£ ¤<br />
<br />
<br />
¤ <br />
¤<br />
<br />
1<br />
<br />
<br />
¨<br />
¤<br />
¡<br />
0<br />
<br />
<br />
§<br />
Résoudre<br />
¨ ¥<br />
<br />
est<br />
¨ ¥<br />
<br />
où<br />
la solution <strong>de</strong><br />
<br />
¡¦ <br />
<br />
¤<br />
<br />
¦<br />
¡ ¡¦ §<br />
¡£¡£¡¤¡£¡¤¡£¡£¡¤¡¤¡ <br />
<br />
¤<br />
<br />
§<br />
<br />
¡<br />
¤<br />
§¦<br />
¡ ¡£¡¤¡£¡¤¡£¡£¡¤¡¤¡<br />
¥<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.12/40<br />
<br />
¡¦ <br />
<br />
¤<br />
<br />
¦<br />
¡ ¡¦ §<br />
¡ ¡£¡¤¡£¡¤¡£¡£¡¤¡¤¡ <br />
<br />
¤<br />
<br />
§<br />
<br />
¡<br />
<br />
¤<br />
§¦<br />
¡¤¡ ¡¤¡£¡¤¡£¡£¡¤¡¤¡<br />
<br />
¡<br />
où<br />
1.5<br />
0.5<br />
−0.5<br />
−1<br />
−1.5<br />
−2<br />
p 2<br />
−3 −2 −1 0 1 2 3<br />
¢¡<br />
£ ¤<br />
u<br />
Exemple 3<br />
¤©§<br />
¦ ©<br />
<br />
¡¦ <br />
<br />
¤<br />
<br />
¦<br />
¢<br />
¡<br />
¡¦<br />
<br />
¤<br />
<br />
¤<br />
© ¦<br />
¡¤¡¤¡£¡£¡¤¡¤¡£¡£¡¤¡£¡¤¡ ¡¤¡¤¡ ¡ <br />
¡<br />
¡¦<br />
¡<br />
¡¦<br />
¡ ¡¤¡¤¡£¡£¡¤¡£¡¤¡£¡£¡¤¡¤¡ ¡¤¡¤¡
Fonction <strong>de</strong> tir<br />
2<br />
2<br />
1<br />
1.5<br />
1<br />
0<br />
0.5<br />
S 1<br />
(z)<br />
−1<br />
S 2<br />
(z)<br />
0<br />
−0.5<br />
−2<br />
−1<br />
−1.5<br />
−3<br />
−4<br />
−2<br />
−4<br />
−2<br />
−4<br />
−2<br />
−4<br />
0<br />
z2<br />
2<br />
4<br />
4<br />
2<br />
0<br />
z1<br />
−2<br />
0<br />
z2<br />
2<br />
4<br />
4<br />
2<br />
0<br />
z1<br />
−2<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.13/40
£<br />
¥<br />
¡<br />
¦<br />
£<br />
¡<br />
¥<br />
¡<br />
¡<br />
¥<br />
£ 0<br />
¥<br />
¨§<br />
¡<br />
¡<br />
¡<br />
¦<br />
¦<br />
¥ ¡<br />
¥<br />
<br />
¨§<br />
¡<br />
¡<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.14/40<br />
Algorithme <strong>de</strong> Newton<br />
Exemple<br />
¡<br />
¥ <br />
£<br />
¤ <br />
¥<br />
¤<br />
¥£¡<br />
¡£<br />
¡<br />
£<br />
¡<br />
¥<br />
0 ¥<br />
¤ <br />
25<br />
20<br />
15<br />
10<br />
S(z)<br />
5<br />
0<br />
z3z2 z1<br />
z0<br />
−5<br />
−10<br />
−1 0 1 2 3 4 5<br />
z
£<br />
¡<br />
¥<br />
£<br />
¡<br />
¥ <br />
¡£<br />
¡<br />
£<br />
¡<br />
¥<br />
¥£¡<br />
¦<br />
¤<br />
¥<br />
¡<br />
¡£<br />
¡<br />
¥<br />
¡£<br />
¡<br />
¨§<br />
¤<br />
¥£¡<br />
*<br />
¦<br />
¤ <br />
¡<br />
¡<br />
¡£¢<br />
£<br />
¡<br />
¥<br />
¥<br />
¥<br />
¢<br />
*<br />
¡<br />
¡<br />
¡<br />
¦<br />
¦<br />
¡<br />
¡<br />
¡£<br />
£<br />
43<br />
¡<br />
£<br />
¡<br />
¦<br />
¡<br />
¥¥¤<br />
¤<br />
¥<br />
43<br />
£<br />
¡<br />
¥<br />
Cas dimension<br />
Objectif: résoudre<br />
£<br />
¡<br />
¥ <br />
où<br />
est une fonction <strong>de</strong><br />
dans<br />
Par définition <strong>de</strong> la dérivée on a au voisinage <strong>de</strong><br />
avec ici<br />
L’algorithme <strong>de</strong> Newton s’écrit alors <strong>de</strong> la même façon:<br />
Remarque . On n’inverse jamais un système linéaire pour résoudre une<br />
équation linéaire !!!<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.15/40
¨<br />
¥<br />
<br />
¡<br />
<br />
<br />
¢¤<br />
¢¤ ¦¥ ¥<br />
¨<br />
¤<br />
£<br />
¥<br />
¤<br />
¤¦¥<br />
¤<br />
Algorithme<br />
Initialisation<br />
choisir<br />
choisir<br />
choisir nbitmax<br />
k=0<br />
Corps<br />
Répéter<br />
Résoudre<br />
k:=k+1<br />
Jusqu’à (k=nbitmax ou<br />
¤¦¥<br />
¢¤<br />
¢¤ ¤ ¥<br />
)<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.16/40
§<br />
¦<br />
¤<br />
¡<br />
¥<br />
¢<br />
¥<br />
¥<br />
§<br />
¢<br />
¦<br />
¡<br />
¥<br />
¡<br />
¡<br />
¡<br />
£<br />
¡ ¡<br />
¨<br />
¡<br />
¦<br />
¦<br />
Théorème<br />
. On suppose que<br />
43<br />
43<br />
¡<br />
-<br />
Soit<br />
est un ouvert convexe <strong>de</strong><br />
43<br />
sur<br />
£<br />
¤<br />
est<br />
inversible<br />
¡£<br />
¡ ¡<br />
¥<br />
et<br />
¡ ¡<br />
¥<br />
£<br />
tel que<br />
Il existe<br />
¡ ( ¡<br />
:<br />
est Lipschitzienne sur<br />
¡ ¡<br />
£<br />
¢<br />
tel que<br />
£ ¤<br />
¢<br />
Il existe<br />
¨ ¥<br />
¡£ © ¥<br />
©¦<br />
£<br />
¡£<br />
¡ ¡<br />
£<br />
¢<br />
£ <br />
© ¥(<br />
¤¦¥ <br />
l’algorithme <strong>de</strong> Newton<br />
alors il existe tel que pour tout dans<br />
converge et la convergence est quadratique:<br />
¡ ¡<br />
£<br />
¢<br />
¥<br />
¤<br />
¡ ¡<br />
©<br />
¡<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.17/40
'<br />
<br />
¢<br />
¡<br />
Problèmes <strong>de</strong> transfert orbital<br />
CNES<br />
Forte excentricité initiale (<br />
Poussées faibles (60 à 0.2 Newtons)<br />
Minimisation <strong>de</strong><br />
Maximisation <strong>de</strong> la masse finale<br />
)<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.18/40
¢<br />
¢<br />
¥£<br />
£<br />
§<br />
§¦5<br />
¥<br />
¦<br />
¡<br />
£ ¥<br />
£<br />
¤<br />
¦ §¨§§§¨§§§¨§§§¨§§§¨§© §§¨§§§¨§§§¨§§§¨§§§¨<br />
¥ ¡<br />
£<br />
¥£<br />
©<br />
<br />
¦<br />
¥£<br />
¨<br />
<br />
¡5<br />
¡<br />
¥<br />
¨<br />
¡ £<br />
¥<br />
¨<br />
£<br />
¡<br />
¡<br />
libre<br />
¨<br />
¥£<br />
¢<br />
¥<br />
¤<br />
¢<br />
£ £<br />
¥ ¢<br />
£<br />
<br />
<br />
¥£<br />
¡<br />
¥ <br />
¤ £ ¡<br />
¥£<br />
¦<br />
¥£<br />
<br />
¡<br />
£¤<br />
<br />
Min<br />
libre<br />
et<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.19/40<br />
£¤<br />
¡ <br />
£ #¢¡<br />
£¤<br />
¥<br />
# £ £ ¥<br />
#¢¨<br />
£ ¥<br />
¥ <br />
£ ¡<br />
avec
&¤ <br />
0 '<br />
(<br />
¨<br />
¢<br />
¡<br />
¡<br />
<br />
£0<br />
¥<br />
¢<br />
¥£¡<br />
£<br />
¡<br />
£<br />
¥<br />
¢<br />
¦£¢<br />
¥<br />
§<br />
¥ ¥£¡<br />
§¦£¢5<br />
©<br />
£ ¡<br />
¦<br />
¦ §§§¨§§§¨§§§¨§§§¨§§§¨§© §§¨§§§¨§§§¨§§§¨§§§¨§§<br />
¥ ¡<br />
£<br />
¥<br />
¤<br />
£ ¡<br />
¥<br />
¨<br />
¡5<br />
¡<br />
¥<br />
¨<br />
¡ £<br />
¥<br />
¨<br />
£<br />
¡<br />
¡<br />
libre<br />
¨<br />
¥£<br />
¢<br />
¥<br />
¤<br />
¥ <br />
£ £<br />
¢<br />
<br />
<br />
¥<br />
¢<br />
¥£¡<br />
£¤<br />
¥<br />
¢<br />
¥£¡<br />
£¤<br />
¡ <br />
£¤ £ ¡<br />
£<br />
¡<br />
£<br />
¡<br />
¡<br />
¥ <br />
£<br />
£<br />
¡<br />
£<br />
¢<br />
<br />
#¢¡<br />
# £<br />
#¢¨<br />
# ¢<br />
£¤<br />
¥<br />
¥<br />
¥ <br />
£ ¡<br />
fixé<br />
et<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.20/40
¤¤<br />
£<br />
¢<br />
<br />
¤ ¥<br />
<br />
¤£¥<br />
<br />
¤ ¥<br />
¢ <br />
<br />
¤£¥<br />
£<br />
¤ ¥<br />
<br />
<br />
¤<br />
¡<br />
<br />
¡<br />
¤<br />
¤<br />
§<br />
¨<br />
<br />
¦ <br />
§<br />
¨<br />
¦<br />
<br />
¥¤<br />
¨<br />
©<br />
¤£¥ ¨<br />
<br />
¥<br />
¢ ¥<br />
¢<br />
¤<br />
§<br />
$ ¦ ¤<br />
¦ <br />
<br />
<br />
<br />
¢<br />
¤<br />
¤¥<br />
¦<br />
§ ¥<br />
<br />
<br />
$ ¦<br />
¢ <br />
<br />
<br />
¨<br />
§<br />
¨<br />
£<br />
£<br />
¨<br />
<br />
Paramètres orbitaux<br />
Il faut, pour <strong>de</strong>s raisons <strong>de</strong> stabilité lors <strong>de</strong> l’intégration numérique travailler avec le système <strong>de</strong><br />
coordonnées <strong>de</strong> Gauss:<br />
où<br />
est le vecteur excentricité,<br />
.<br />
est le vecteur rotation du plan <strong>de</strong> l’orbite,<br />
est la longitu<strong>de</strong> vraie,<br />
Z<br />
satellite<br />
périgée<br />
plan équatorial<br />
w<br />
ω<br />
X<br />
orbite<br />
Ω<br />
Y<br />
i<br />
Figure : Coordonnées <strong>de</strong> Gauss<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.21/40
¡<br />
¡ ¡<br />
¤<br />
£<br />
¢<br />
¡<br />
¥<br />
£<br />
¡<br />
£<br />
¡<br />
¥<br />
£<br />
¡<br />
Repère ortho-radial<br />
k<br />
j<br />
r<br />
w<br />
S<br />
v<br />
s<br />
q<br />
O<br />
i<br />
Figure : Repère ortho-radial<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.22/40
¢<br />
¥£<br />
¥£<br />
<br />
¦<br />
¤<br />
¦ §¨§§§¨§§§¨§§§¨§§§¨§§§¨§§§¨© §¨§¨§¨§¨§¨§¨§¨§¨§¨§¨§¨§¨§¨§¨§¨§¨§¨§¨§¨<br />
¥ ¡<br />
£<br />
<br />
¡5<br />
¡<br />
¥¦ <br />
¤§¦<br />
¥<br />
¤<br />
1<br />
00<br />
¥¢<br />
£<br />
¡<br />
¤<br />
Problème en coordonnées <strong>de</strong><br />
Gauss<br />
&¤ <br />
0 '<br />
£0<br />
¥<br />
(<br />
<br />
'<br />
¢<br />
¥£<br />
£ <br />
¥£<br />
¥ <br />
¤ 5<br />
§¦ ¢<br />
¥<br />
&¡<br />
£ <br />
¥£<br />
#<br />
£ ¥<br />
¥£<br />
'<br />
¢<br />
&©<br />
# ¨<br />
£ ¥<br />
# ¢ £ ¥<br />
¥£<br />
¨<br />
¥<br />
¤<br />
£¤<br />
£¤<br />
£<br />
¦<br />
¨ <br />
¤ ¡<br />
£¤<br />
¥ <br />
¤¤¦<br />
£<br />
<br />
¥<br />
0<br />
£¤<br />
©<br />
<br />
¨<br />
£¤<br />
¥ <br />
¨<br />
¤ <br />
¦<br />
£¤<br />
¥<br />
¤¤1<br />
0¢¥ <br />
¡<br />
£¤<br />
¥ <br />
¥<br />
¤<br />
£0<br />
¥<br />
¤ <br />
£0<br />
¥0<br />
¦¢<br />
¤§¦<br />
¨ <br />
¤ ¡<br />
£0<br />
¥ <br />
¦<br />
£0<br />
¥<br />
¥<br />
¤ <br />
£0<br />
¦<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.23/40
§<br />
¡<br />
*<br />
¡<br />
<br />
¦<br />
¤<br />
£<br />
<br />
¤<br />
<br />
¤¤<br />
<br />
<br />
¡<br />
¨<br />
¤¤<br />
<br />
¡<br />
¨<br />
¡<br />
<br />
¤<br />
<br />
¡ <br />
<br />
¤<br />
© §¥¤<br />
<br />
¤<br />
§<br />
¡<br />
<br />
¡<br />
<br />
<br />
¤<br />
<br />
¤<br />
<br />
¡ <br />
<br />
¡<br />
¤<br />
<br />
¡ <br />
<br />
<br />
¡<br />
<br />
<br />
¦<br />
§<br />
¡<br />
<br />
<br />
<br />
<br />
"<br />
¤<br />
<br />
¡ <br />
<br />
<br />
¡<br />
<br />
§<br />
¡<br />
<br />
<br />
¦<br />
<br />
<br />
<br />
<br />
<br />
¤ <br />
<br />
<br />
<br />
Minimisation <strong>de</strong> l’Hamiltonien<br />
<br />
¤¤<br />
¤£¥<br />
¤£¥<br />
¤£¥<br />
¢ <br />
<br />
§¥¤<br />
£ ¥<br />
! ¢<br />
¥<br />
<br />
¤<br />
¦ ¤<br />
§<br />
<br />
¤<br />
¤¤§¤<br />
<br />
¤ '<br />
¤<br />
¦ <br />
©<br />
¤©§<br />
¤¤¦<br />
¢<br />
¤¤© §<br />
£¢ ¤<br />
<br />
<br />
¤¤<br />
¤¤§<br />
¢<br />
<br />
¤¤<br />
§¥¤<br />
¤ ¥<br />
¤£¥<br />
¤£¥<br />
¢ <br />
alors<br />
<br />
% <br />
¤¤§<br />
Si<br />
¤¤§<br />
<br />
<br />
<br />
¤<br />
alors<br />
<br />
¤¤<br />
% <br />
§¥¤<br />
¤£¥<br />
¤£¥<br />
¤£¥<br />
¢ <br />
¥<br />
<br />
<br />
<br />
<br />
Si<br />
¤¤§<br />
<br />
<br />
<br />
¨ ¥<br />
<br />
¤<br />
alors<br />
<br />
¤¤<br />
§¤<br />
¤£¥<br />
¤£¥<br />
¤£¥<br />
¢ <br />
¥<br />
<br />
<br />
<br />
<br />
Si<br />
Conditions <strong>de</strong> transversalité<br />
§<br />
<br />
¤<br />
<br />
¤<br />
§<br />
<br />
¥<br />
§ ¤<br />
¥<br />
§<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.24/40
Résultats numériques<br />
Le calcul <strong>de</strong> la fonction <strong>de</strong> tir nécessite l’intégration numérique d’un<br />
système différentiel ordinaire<br />
La dérivée <strong>de</strong> la fonction <strong>de</strong> tir est approchée numériquement par<br />
différences finies<br />
Il faut faire <strong>de</strong> la mise à l’échelle (scaling) <strong>de</strong>s données<br />
Logiciel tfmin <strong>de</strong> J.B. Caillau<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.25/40
Résultats numériques: 60N<br />
P<br />
40<br />
20<br />
0<br />
0 5 10<br />
1<br />
t<br />
L<br />
10<br />
5<br />
0<br />
0 5 10<br />
0.2<br />
t<br />
e x<br />
0<br />
−1<br />
0 5 10<br />
0.5<br />
t<br />
h x<br />
0.1<br />
0<br />
0 5 10<br />
0.02<br />
t<br />
r 3<br />
−0.50<br />
−1<br />
e y<br />
m<br />
0<br />
−0.5<br />
0 5 10<br />
1500<br />
t<br />
1400<br />
h y<br />
0<br />
−0.02<br />
0 5 10<br />
t<br />
1300<br />
0 2 4 6 8 10 12 14<br />
1000<br />
t<br />
20<br />
10<br />
0<br />
r 2<br />
−10<br />
−20<br />
−40<br />
−30<br />
−20<br />
r 1<br />
−10<br />
0<br />
10<br />
u 1<br />
u 2<br />
0<br />
−1000<br />
0 2 4 6 8 10 12 14<br />
1000<br />
t<br />
0<br />
−1000<br />
0 2 4 6 8 10 12 14<br />
500<br />
t<br />
r 2<br />
20<br />
10<br />
0<br />
−10<br />
r 3<br />
2<br />
1<br />
0<br />
−1<br />
u 3<br />
0<br />
−500<br />
0 2 4 6 8 10 12 14<br />
t<br />
−20<br />
−40 −30 −20 −10 0 10<br />
r 1<br />
−2<br />
−40 −20 0 20 40<br />
r 2<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.26/40
Résultats numériques: 60N<br />
0.5<br />
0<br />
p P<br />
0<br />
p L<br />
−5<br />
−0.5<br />
0 5 10<br />
40<br />
t<br />
−10<br />
0 5 10<br />
2.29<br />
t<br />
p ex<br />
20<br />
p hx<br />
2.285<br />
p ey<br />
0<br />
0 5 10<br />
10<br />
t<br />
0<br />
−10<br />
0 5 10<br />
t<br />
p hy<br />
0 5 10<br />
t<br />
−0.826<br />
−0.828<br />
−0.83<br />
−0.832<br />
−0.834<br />
0 5 10<br />
t<br />
4 x 10−5 t<br />
3<br />
µ<br />
2<br />
1<br />
0<br />
0 2 4 6 8 10 12 14<br />
0.5<br />
ψ 3<br />
0<br />
−0.5<br />
1<br />
0<br />
−1<br />
−2<br />
−0.5<br />
0<br />
0.5<br />
ψ 2<br />
ψ 1<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.27/40
Résultats numériques: 3N<br />
P<br />
40<br />
20<br />
0<br />
0 50 100 150 200 250<br />
1<br />
t<br />
L<br />
200<br />
100<br />
0<br />
0 50 100 150 200 250<br />
0.1<br />
t<br />
e x<br />
0.5<br />
h x<br />
0<br />
0<br />
0 50 100 150 200 250<br />
0.02<br />
t<br />
−0.1<br />
0 50 100 150 200 250<br />
t<br />
r 3<br />
40<br />
−2 02 r 1<br />
e y<br />
0<br />
−0.02<br />
0 50 100 150 200 250<br />
1500<br />
t<br />
h y<br />
1 x 10−3 0<br />
t<br />
−1<br />
0 50 100 150 200 250<br />
20<br />
0<br />
0<br />
20<br />
40<br />
m<br />
1400<br />
−20<br />
−20<br />
1300<br />
0 50 100 150 200 250<br />
50<br />
t<br />
u 1<br />
0<br />
40<br />
r 2<br />
−40<br />
−40<br />
2<br />
−50<br />
0 50 100 150 200 250<br />
50<br />
t<br />
20<br />
1<br />
u 2<br />
0<br />
r 2<br />
0<br />
r 3<br />
0<br />
−50<br />
0 50 100 150 200 250<br />
50<br />
t<br />
−20<br />
−1<br />
u 3<br />
0<br />
−50<br />
0 50 100 150 200 250<br />
t<br />
−40<br />
−60 −40 −20 0 20 40<br />
r 1<br />
−2<br />
−40 −20 0 20 40<br />
r 2<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.28/40
Résultats numériques: 0.2N<br />
P<br />
40<br />
20<br />
0<br />
0 1000 2000 3000 4000<br />
1<br />
t<br />
L<br />
2000<br />
1000<br />
0<br />
0 1000 2000 3000 4000<br />
0.1<br />
t<br />
e x<br />
0.5<br />
h x<br />
0<br />
1<br />
0<br />
0 1000 2000 3000 4000<br />
t<br />
0<br />
−1<br />
0 1000 2000 3000 4000<br />
1500<br />
t<br />
e y<br />
−0.1<br />
0 1000 2000 3000 4000<br />
t<br />
h y<br />
1 x 10−4 0<br />
t<br />
−1<br />
0 1000 2000 3000 4000<br />
r 3<br />
40<br />
20<br />
0<br />
0<br />
20<br />
40<br />
−2 02 r 1<br />
m<br />
1400<br />
−20<br />
−20<br />
1300<br />
0 500 1000 1500 2000 2500 3000 3500 4000<br />
5<br />
t<br />
u 1<br />
0<br />
40<br />
r 2<br />
−40<br />
−40<br />
2<br />
−5<br />
0 500 1000 1500 2000 2500 3000 3500 4000<br />
5<br />
t<br />
20<br />
1<br />
u 2<br />
0<br />
r 2<br />
0<br />
r 3<br />
0<br />
−5<br />
0 500 1000 1500 2000 2500 3000 3500 4000<br />
5<br />
t<br />
−20<br />
−1<br />
u 3<br />
0<br />
−5<br />
0 500 1000 1500 2000 2500 3000 3500 4000<br />
t<br />
−40<br />
−60 −40 −20 0 20 40<br />
r 1<br />
−2<br />
−40 −20 0 20 40<br />
r 2<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.29/40
¡<br />
¢<br />
¡5<br />
¡5<br />
¨<br />
<br />
¢<br />
£<br />
¨¢<br />
¤<br />
¤<br />
©<br />
¤<br />
Relation<br />
,<br />
60N 14.8 h<br />
3N 285.6 h<br />
0.2N 4260.6 h<br />
12 jours<br />
6 mois<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.30/40
¥£<br />
"<br />
<br />
<br />
<br />
¢<br />
¥£<br />
¡ £ ¥<br />
#<br />
£<br />
¥£<br />
¥<br />
<br />
§¦ 5 <br />
¤ £ ¡<br />
¥£<br />
# £ £ ¥<br />
#¢¨<br />
£ ¥<br />
¦ §¨§§§¨§§§¨§§§¨§§§¨§© §§¨§§§¨§§§¨§§§¨§§§¨<br />
¥ ¡<br />
£<br />
¥£<br />
©<br />
<br />
¦<br />
¡5<br />
¨ 0<br />
¥£<br />
<br />
¡<br />
¥ <br />
# ¡<br />
¡ <br />
¥<br />
¥ <br />
¨<br />
¨<br />
¡<br />
¡<br />
libre<br />
¨<br />
¥£<br />
¢<br />
¤<br />
¥¢<br />
£<br />
¢ ¥ £<br />
<br />
¥£<br />
¡<br />
¥ <br />
¤ £ ¡<br />
¥£<br />
¦<br />
¥£<br />
<br />
¡<br />
<br />
©<br />
"<br />
<br />
<br />
¥£<br />
¨<br />
¥£ ¡<br />
¨<br />
£<br />
¢ ¥<br />
¦<br />
¡5<br />
Max<br />
et<br />
fixé<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.31/40<br />
£¤<br />
#¢¡<br />
£¤<br />
# ¡<br />
¡ <br />
£¤<br />
£ ¡<br />
avec<br />
Remarque .
§<br />
¦<br />
¡<br />
¨<br />
<br />
¤<br />
<br />
§ §¤<br />
<br />
©<br />
¨<br />
£<br />
¦<br />
<br />
<br />
§<br />
©<br />
<br />
¤¤<br />
<br />
§<br />
¡<br />
<br />
<br />
<br />
§<br />
¢<br />
<br />
<br />
¥<br />
<br />
<br />
¦<br />
§<br />
¡<br />
<br />
¥<br />
<br />
<br />
"<br />
<br />
<br />
<br />
§<br />
¡<br />
<br />
<br />
¦<br />
<br />
<br />
¥<br />
<br />
<br />
<br />
¡<br />
<br />
¦<br />
§<br />
¡<br />
§<br />
<br />
¤<br />
<br />
¤<br />
©<br />
¤<br />
"<br />
¦<br />
<br />
¥<br />
§ §¥¤<br />
<br />
<br />
¤<br />
©<br />
¡<br />
"<br />
¦<br />
<br />
<br />
¤<br />
§ §¤<br />
<br />
<br />
<br />
<br />
¤<br />
©<br />
<br />
<br />
¡<br />
<br />
¦<br />
§ §¥¤<br />
<br />
<br />
Minimisation <strong>de</strong> l’Hamiltonien<br />
L’Hamiltonien est<br />
<br />
¤¤<br />
§¥¤<br />
¤£¥<br />
¤£¥<br />
¥ ¤<br />
¤£¥<br />
¢ <br />
¥<br />
<br />
¤¤<br />
<br />
¤©§¡£<br />
¤¤© §<br />
<br />
¤¤<br />
¢<br />
¡<br />
<br />
<br />
¤¤<br />
¨<br />
<br />
¤©§¡<br />
¤<br />
¦ <br />
¤¤<br />
¤¤<br />
¦ ¤<br />
<br />
<br />
¤<br />
% <br />
¦ ¤<br />
§<br />
<br />
¤¤<br />
§¤<br />
¤<br />
<br />
¤<br />
§ ¤<br />
¤ ¦<br />
§¤<br />
¤£¥<br />
¥ ¤<br />
¤ ¥<br />
¢ <br />
<br />
¤<br />
<br />
¤<br />
§¡<br />
<br />
alors<br />
<br />
¤<br />
<br />
¤¤<br />
§¤<br />
¤£¥<br />
¤£¥<br />
¤£¥<br />
¢ <br />
<br />
Si<br />
¨ ¥<br />
<br />
¤<br />
§¡<br />
<br />
alors<br />
<br />
¤<br />
<br />
¤¤<br />
§¤<br />
¤£¥<br />
¤£¥<br />
¤£¥<br />
¢ <br />
<br />
Si<br />
alors<br />
<br />
¤<br />
<br />
¤¤<br />
§¥¤<br />
¤£¥<br />
¤£¥<br />
¤£¥<br />
¢ <br />
Si<br />
<br />
<br />
¤<br />
alors<br />
<br />
¤<br />
¥<br />
¦ ¤<br />
Si<br />
alors<br />
<br />
¤<br />
¥<br />
¦ ¤<br />
Si<br />
alors<br />
<br />
¤<br />
¦ ¤<br />
Si<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.32/40
Difficultés<br />
Le contrôle Bang-Bang induit un second membre d’équation<br />
différentiel discontinu.<br />
Dérivabilité <strong>de</strong> la fonction <strong>de</strong> tir?<br />
Le tir simple diverge ici<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.33/40
Resultats pour 2.5N,<br />
40<br />
30<br />
20<br />
10<br />
0 200 400<br />
1500<br />
1450<br />
1400<br />
1350<br />
0 200 400<br />
6<br />
4<br />
2<br />
1<br />
0.5<br />
0<br />
−0.5<br />
−1<br />
0 200 400<br />
0.1<br />
0.05<br />
0<br />
−0.05<br />
−0.1<br />
0 200 400<br />
1<br />
0.5<br />
0<br />
−0.5<br />
1<br />
0.5<br />
0<br />
−0.5<br />
−1<br />
0 200 400<br />
1<br />
0.5<br />
0<br />
−0.5<br />
−1<br />
0 200 400<br />
0<br />
0 200 400<br />
−1<br />
0 50 100 150 200 250 300 350 400<br />
Figure : état et contrôle<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.34/40
Resultats pour 2.5N,<br />
1<br />
0.5<br />
0<br />
−0.5<br />
−1<br />
0 50 100 150 200 250 300 350 400 450<br />
1<br />
0.5<br />
0<br />
−0.5<br />
−1<br />
0 50 100 150 200 250 300 350 400 450<br />
2<br />
1<br />
0<br />
−1<br />
0 50 100 150 200 250 300 350 400 450<br />
Figure : contrôle<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.35/40
£<br />
<br />
£ $<br />
)<br />
<br />
<br />
<br />
<br />
4<br />
£<br />
<br />
£ !<br />
4<br />
¤<br />
<br />
4<br />
2<br />
(<br />
¡<br />
Métho<strong>de</strong>s directes<br />
Problème <strong>de</strong> Mayer<br />
¥¥£<br />
¢<br />
£<br />
¢<br />
<br />
¥<br />
<br />
<br />
£ <br />
¥ <br />
<br />
¥¥£<br />
#<br />
£ ¥<br />
¥ ( £<br />
£ ¡ ¥<br />
¡£ ¡<br />
<br />
¦ §¨§§§¨§§§¨§© §§¨§§§¨§§§¨<br />
¥£¤<br />
¢ 3<br />
¤(<br />
¢3<br />
¤(<br />
£<br />
¢ ¥<br />
¢£<br />
¢<br />
<br />
Remarque . On peut toujours se ramener à un problème <strong>de</strong> Mayer. Il<br />
suffit <strong>de</strong> définir un état tel que:<br />
¥<br />
¤<br />
¥ <br />
<br />
¥¥£<br />
#<br />
¤ ¥<br />
¥<br />
£<br />
¡¤<br />
¥<br />
¨ ¤<br />
<br />
£<br />
©£¢<br />
¥<br />
53<br />
0 <br />
111<br />
) <br />
/ <br />
On prendra<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.36/40
¤<br />
<br />
¨<br />
¤<br />
<br />
<br />
¢<br />
<br />
¤ ¢<br />
©¢ ©<br />
¨<br />
¤¤<br />
¢<br />
¨<br />
¤<br />
<br />
¤ ¢<br />
¢<br />
<br />
¤<br />
$<br />
¢<br />
¨<br />
$<br />
<br />
¢<br />
¨<br />
¤<br />
$<br />
$<br />
¢<br />
<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.37/40<br />
Discrétisation d’une EDO<br />
<br />
¤¤<br />
¡<br />
<br />
¢<br />
¥<br />
<br />
¢&<br />
£ ¤<br />
¨ <br />
fixé)<br />
(<br />
'<br />
'<br />
£<br />
¥<br />
¥<br />
¢ ¢¢<br />
<br />
¥<br />
¨<br />
(pas constant)<br />
$<br />
<br />
¤ <br />
$ <br />
Schéma d’Euler<br />
$ ¥<br />
$<br />
$<br />
$<br />
¨<br />
¤ <br />
¤<br />
$ <br />
On pose<br />
¤ <br />
$<br />
$ ¥<br />
¥<br />
¢<br />
L’erreur locale (c’est-à-dire l’erreur d’un pas) <strong>de</strong> ce schéma est en<br />
Métho<strong>de</strong> <strong>de</strong>s trapèzes<br />
On prend le schéma:<br />
¤ <br />
¤¤<br />
$<br />
¨<br />
¢<br />
<br />
¤ <br />
$<br />
¤ ¥<br />
$ <br />
$ ¥<br />
¢<br />
¥<br />
¢<br />
¦¤<br />
L’erreur locale est en<br />
¦¤ ¢
¥<br />
*<br />
¤<br />
¡<br />
0<br />
¤<br />
<br />
<br />
<br />
¦<br />
¦<br />
¦<br />
111<br />
¡<br />
0<br />
©£¢<br />
¥<br />
* £<br />
¡<br />
¦<br />
111<br />
111<br />
<br />
<br />
Métho<strong>de</strong> directe<br />
On discrétise l’équation d’état et on obtient le problème d’optimisation en<br />
dimension finie suivant:<br />
<br />
£ <br />
pour<br />
¥<br />
*<br />
£ *<br />
* <br />
<br />
$<br />
¥<br />
*<br />
<br />
¤ <br />
¢ <br />
0 <br />
<br />
¤ <br />
pour<br />
¦ §§§¨§§§¨§§© §¨§¨§¨§¨§¨§¨§¨§¨<br />
¨ ¤<br />
¡¥<br />
¤<br />
¡£ ¡<br />
<br />
¥ <br />
¤<br />
¢£<br />
¢<br />
<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.38/40
¥£<br />
&<br />
<br />
'<br />
¥£<br />
<br />
¡<br />
<br />
£<br />
¥<br />
£<br />
<br />
*<br />
<br />
<br />
¤<br />
¢ ¥<br />
£<br />
<br />
<br />
¤<br />
¡<br />
£<br />
£<br />
<br />
©£¢<br />
¥<br />
* £<br />
¡<br />
0<br />
¡<br />
¦<br />
111<br />
111<br />
<br />
<br />
Tir direct<br />
On réduit l’espace <strong>de</strong>s contrôles:<br />
¥<br />
¤<br />
constant sur<br />
& *<br />
*<br />
¥<br />
¤<br />
affine sur<br />
& *<br />
*<br />
...<br />
et<br />
Cas <strong>de</strong>s contrôles constants par morceaux On peut alors en fonction <strong>de</strong><br />
<strong>de</strong>s calculer , puis ,... ,et enfin<br />
¥£<br />
¤<br />
¤¥<br />
¡ 1 11<br />
<br />
£ <br />
¡ <br />
<br />
¥<br />
£<br />
D’où le problème d’optimisation en dimension finie<br />
¤¥¥<br />
¡ 1 11<br />
<br />
<br />
<br />
pour<br />
¨ ¤<br />
¢ <br />
0 <br />
<br />
¤ <br />
¡¥<br />
¤<br />
¡£ ¡<br />
<br />
¦ §¨§§§¨§© §§§¨§§<br />
¥<br />
¤<br />
¢£<br />
¢<br />
<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.39/40
Remarques<br />
Les métho<strong>de</strong>s directes sont plus robustes<br />
Il est plus facile d’introduire <strong>de</strong>s contraintes sur l’état avec les<br />
métho<strong>de</strong>s directes<br />
Les résultats <strong>de</strong>s métho<strong>de</strong>s directes sont moins précis, en particulier sur<br />
le contrôle<br />
L’interaction mathématique-numérique est fondamentale ici.<br />
Cette présentation peut aussi concerner le point discret, continu<br />
TIPE contrôle <strong>optimal</strong>, 2002-2003 – p.40/40