0% ont trouvé ce document utile (0 vote)
45 vues6 pages

Corrigé TD3 - Calcul Variationnel

Le document présente des exercices de mathématiques appliquées à l'optimisation et à la théorie des fonctions. Chaque exercice aborde des problèmes d'infimum de fonctionnelles avec des conditions aux limites variées, en utilisant des méthodes comme l'équation d'Euler-Lagrange. Les solutions détaillées incluent des calculs d'intégrales et des démonstrations de convergence vers des valeurs spécifiques.

Transféré par

kenza.chebil.ch06
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
0% ont trouvé ce document utile (0 vote)
45 vues6 pages

Corrigé TD3 - Calcul Variationnel

Le document présente des exercices de mathématiques appliquées à l'optimisation et à la théorie des fonctions. Chaque exercice aborde des problèmes d'infimum de fonctionnelles avec des conditions aux limites variées, en utilisant des méthodes comme l'équation d'Euler-Lagrange. Les solutions détaillées incluent des calculs d'intégrales et des démonstrations de convergence vers des valeurs spécifiques.

Transféré par

kenza.chebil.ch06
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

École Supérieure de la statistique et de l’Analyse de l’Informatque 2éme Année

Période 4 2019-2020

Corrigé de TD3

Exercice 1
On s’intéresse au problème suivant :
Z 1
inf J(y) = xy 02 (x)dx, y(0) = 1, y(1) = 0.
0

1. Calculer J(yn ) avec xn la fonction valant 1 sur [0, n1 [


log(x)
avec yn (x) = − log(n) pour x ∈ [ n1 , 1].
2. Monter que l’infinium du problème est 0.

Corrigé
Z 1
inf J(y) = xy 02 (x)dx
0

1. Calculer J(yn )
1
Z 1 Z Z 1  
n 1
J(yn ) = xyn02 dx = 0dx + x dx
0 0 n
1 x2 log(n)2
Z 1
1 1
= 2
dx
log(n) n1 x
1
J(yn ) = , ∀n ∈ N.
log(n)

2. Rappel : un élément a d’un ensemble A est dit une borne inférieure si et seulement si :
— a est un minorant de A.
— Il existe une suite (xn ) ⊂ A telle que xn →n→+∞ a
On a 0 est un minorant de J(y) car J(y) > 0.
Il faut trouver une suite (zn )n ⊂ A telle que zn → 0 quand n tend vers +∞ avec
A = {z ∈ R| z = J(y), ∀y ∈ C 1 ([0, 1])}, (zn )n ⊂ A c’est à dire il existe yn tel que
zn = J(yn ) donc on doit trouver yn telle que J(yn ) → 0 quand n → +∞.
1
Or d’après la première question on a J(yn ) = log(n) →n→+∞ 0. Alors 0 est un infinimum
du problème.

Exercice 2
Résoudre le problème suivant :
Z 1
inf J(y) = (y 02 (x) + y 2 (x))dx
0

dans les cas suivants :

1
1. Sans conditions aux limites.
2. Avec les conditions y(0) = 0, y(1) = 1.
3. Avec la condition y(0) = 1.

Corrigé
F (x, y, y 0 ) = y 2 + y 02
L’équation d’Euler-Lagrange :
d
Fy − Fy0 = 0.
dx
Donc on a
d 0
2y − 2y = 0
dx
Si et seulement si
y 00 − y = 0. (1)
L’équation caractéristique de (1) est r2 − 1 = 0 qui a pour solution r1 = 1 et r2 = −1. Donc la
solution générale de (1) est y(x) = Aex + Be−x .
1. Sans condition aux bords : ona Fy0 |x=0 = 0 et Fy0 |x=1 = 0. Donc on a

y 0 (0) = 0 et y 0 (1) = 0

Donc on obtient le système suivant :



A−B =0
Ae1 − Be−1 = 0.

Donc on obient A = B = 0. D’où y(x) = 0


2. Avec les conditions y(0) = 0 et y(1) = 1 : on a y(x) = Aex + Be−x donc on obtient le
système suivant : 
A+B =0
Ae + Be−1 = 1.
1

2sh(x)
Donc on obtient y(x) = 1
e1 −e−1
(ex − e−x ) = e1 −e−1
.
x −x
3. Avec la condition y(0) = 1 : on a y(x) = Ae + Be . Donc on a :

y(0) = 1
Fy0 |x=1 = 0.

Donc on obtient 
A+B =1
Ae − Be−1 = 0.
1

e−1 e1
D’où y(x) = e1 +e−1
ex + e1 +e−1
e−x .

Exercice 3
Résoudre le problème suivant :
Z 1
inf J(y) = (y 02 (x) + xy(x))dx + y 2 (1)
0

dans les cas suivants :

2
1. Sans conditions aux limites.
2. Avec la condition y(0) = 1.

Corrigé
On a : Fy = x, Fy0 = 2y 0 donc l’équation d’Euler-Lagrange est

x − 2y 00 = 0

si et seulement si
1
y 00 (x) = x
2
si et seulement si
1
y 0 (x) = x2 + C1
4
si et seulement si
1 3
y(x) = x + C1 x + C2 .
12
1. on a 
Fy0 |x=0 = 0
(Fy0 + fy )|x=1 = 0.
si et seulement si
y 0 (0) = 0

0
2y (1) + 2y(1) = 0.
Donc on trouve C1 = 0 et C2 = − 31 . Donc y(x) = 1 3
12
x − 13 .
1 3
2. On a y(x) = 12
x + C1 x + C2 et y(0) = 1 donc on a

y(0) = C2 = 1
2y 0 (1) + 2y(1) = 0.

Donc C1 = − 32 et C2 = 1
1 2
y(x) = 12 x − 23 x + 1

Exercice 4
Déterminer les extrémales de la fonctionnelle
Z 1 02
y
J(y) = ( + yy 0 + y 0 + y)dx
0 2
lorsqu’aucune condition au bord n’est fixée en x = 0 ou en x = 1.

Corrigé

On a Fy = y 0 + 1, Fy0 = y 0 + y + 1 donc l’équation d’Euler-Lagrange est :

y 00 (x) = 1

si et seulement si
y 0 (x) = x + C1
si et seulement si
x2
y(x) = + C1 x + C2
2
3
On a les deux condtions aux bords sont libres, donc on a

Fy0 |x=0 = 0
Fy0 |x=1 = 0.

si et sulement si 
C1 + C2 + 1 = 0
3
2
+ 2C1 + C2 = 0.
Donc C1 = − 25 et C2 = 72 .

Exercice 5
Soit Z π
q1 (q1 + 3q20 ) + q2 (q2 − q10 ) − 2x(q1 + q2 ) dx

J(q) =
0
Donner la forme générale des éxtremales de J.
On fixe les conditions au bord en x = 0 aux valeurs q1 (0) = 0 et q2 (0) = 2.
Déterminer q1 (π) et q2 (π).

corrigé :
On a F (x, q1 , q2 , q10 , q20 ) = q1 (q1 + 3q20 ) + q2 (q2 − q10 ) − 2x(q1 + q2 )
On note y = (q1 , q2 ) et y 0 = (q10 , q20 ) et par suite, on peut écrire

F (x, y, y 0 ) = q1 (q1 + 3q20 ) + q2 (q2 − q10 ) − 2x(q1 + q2 )

l’équation d’Euler Lagrange :



d 0 Fq1 − d Fq0 = 0

Fy − Fy = 0 ⇔ dx 1
dx d
Fq2 − Fq0 = 0

dx 2
( (
2q1 + 4q20 = 2x (1) (2) ⇔ q2 = x + 2q10
⇔ ⇔
2q2 − 4q10 = 2x (2) (1) ⇔ q1 − 4q100 = x + 2
de même (
(1) ⇔ q1 = x − 2q20
(2) ⇔ q2 + 4q200 = x − 2

alors on cherche à trouver q1 et q2 tels que


(
q1 + 4q100 = x − 2
q2 − 4q200 = x + 2

On commence par résoudre le système avec les équation homogénes


( (
q1 + 4q100 = 0 q2 (x) = Ae−x/2 + Bex/2

q2 − 4q200 = 0 q1 (x) = Ccos(x/2) + Dsin(x/2)

4
il est claire que q1p (x) = x − 2 et q2p (x) = x + 2 Donc
(
q2 (x) = Ae−x/2 + Bex/2 + x + 2
q1 (x) = Ccos(x/2) + Dsin(x/2) + x − 2
(
q1 (0) = C − 2 = 0, ⇒ C = 2
on a ⇔
q2 (0) = A + B + 2 = 2 ⇒ A + B = 0 ⇒ A = −B
(
C = 2 ⇒ q1 (x) = 2cos(x/2) + Dsin(x/2) + x − 2
A+B =0

Pour déterminer q1 (π) et q2 (π) on utilise les conditions naturelles suivantes


(
Fq10 |x=π = −q2 (π) = 0
Fq20 |x=π = 3q1 (π) = 0
alors q2 (π) = q1 (π) = 0 par suite, on peut déterminer q1 (x) et q2 (x)

q (π) = 0 ⇒ A(e π2 − e−π/2 ) + π + 2 = 0 ⇒ A = − π + 2 sh(π/2)
2
2
q (π) = 2cos(π/2) + Dsin(π/2) + π − 2 = 0 ⇒ D = −π + 2
1

Finalement,
q (x) = − (π + 2)sh(x) + x + 2

2
sh(π/2)
q1 (x) = 2cos(x/2) − (π − 2)sin(x/2) + x − 2.

Exercice 6
Soit Z 1
J(y) = y 02 dx
0
sous la contrainte Z 1
K1 (y) = y 2 dx = 2.
0

Donner les extrémales du problème y(0) = 0 et la valeur x = 1 est libre puis pour
Z 1
K2 (y) = xydx = 1.
0

Corrigé :
R1 R1
J(y) = 0 y 02 (x)dx sous la contrainte K1 (y) = 0 y 2 dx = 2.

1. L’équation d’Euler-Lagrange :
∂ 02 d ∂ 02
[y − λy 2 ] − [y − λy 2 ] = 0
∂y dx ∂y 0
si et seulement
λy + y 00 = 0
On distingue trois cas :
(a) Si λ = 0 on a y 00 = 0 donc y(x) = C1 x + C2 .

5
√ √
−λx
(b) λ < 0 on a y(x) = C1 e + C2 e− −λx .
√ √
(c) Si λ > 0 on a y(x) = C1 cos( λx) + C2 sin( λx).
Pour le premier cas on a 
y(0) = 0
Fy0 |x=1 = 0.
si et seulement si 
C2 = 0
y 0 (1) = C1 = 0 = 0.
donc y(x) = 0. De même pour les deux autres cas on trouve que y(x) = 0.
R1
2. La contrainte K2 (y) = 0 xydx = 1. L’équation d’Euler-Lagrange est :

∂ 02 d ∂ 02
[y − λxy] − [y − λxy] = 0
∂y dx ∂y 0
donc on a
λ
y 00 (x) = − x
2
si et seulement si
λ
y 0 (x) = − x2 + C1
4
si et suelement si
λ 3
y(x) = − x + C1 x + C3
12

y(0) = 0 ⇒ C3 = 0
0
y (1) = − λ4 + C1 = 0 ⇒ C1 = λ4 .
Il faut remplacer dans la contrainte pour déterminer la valeur de λ et par suite la valeur
de C1
On a
Z 1 Z 1  
λ 3 λ
K2 (y) = xydx = x x − x dx
0 0 12 4
 1  1
λ 4 λ 2
= x − x
48 0 8 0
=1
−48
Donc λ = 5
.
Corrigé de l’exercice 9 :
L’équation d’Euler-Lagrange est :

y 00 (x) = x
si et seulement si
2 3
y 0 (x) = x 2 + C1
3
si et seulement si
4 5
y(x) = x 2 + C1 x + C2
15
On a y(0) = 0 et la condition en 1 est libre donc on a y(0) = C2 = 0 et Fy0 |x = 1 = 0 donc on
4 52
obtient C1 = − 32 et C2 = 0. D’où y(x) = 15 x − 23 x.

Vous aimerez peut-être aussi