Chap 29 : Equations différentielles linéaires
Chap 29 : Equations différentielles linéaires
I. EDL scalaire
J intervalle de , A, B, C C ( J , ) ( ) : A( x) y '( x) B( x) y( x) C ( x)
Un point x0 J tel que A( x0 ) 0 est une singularité de l'équation différentielle
On se place sur I J intervalle non trivial tel que x I , A( x) 0
( ) devient : y ' a( x) y( x) b( x) L'ED homogène associée est : ( 0 ) y '( x) a( x) y ( x) 0
y ' ay b est une forme ordinaire/résolue
Problème de Cauchy : étant donné ( x0 , y0 ) I , existence et unicité de y solution de ( ) tq y( x0 ) y0
0 {Solutions de ( 0 )} C'est un sev de C 1 ( I , )
primitive de a sur I . L'ensemble des solutions de ( 0 ) est {x Ce ( x ) / C }
En pratique : " y '/ y a d'où ln y ( x) c " etc
y p sol. particulière de ( ) S y p 0
Variation des constantes : on cherche une solution particulière de ( ) sous la forme :
x
y ( x) Ce ( x ) , C C 1 ( I , ) C est une primitive de x e ( x ) b( x) : C ( x) e (t ) b(t )dt K 0
x0
x
Solution générale : y ( x) e ( x ) e (t ) b(t )dt C0 e (t )
x0
Le problème de Cauchy est résolu : il y a existence et unicité
f C1 ( , ), f f ' l f
l (l 0, via ED : ( x) f ( x) f '( x), sol, comparaison)
(a, b) C (, ), a T périodique y ' ay b possède une unique sol T périodique ssi (T ) 0
Si singularité : on regarde à la limite pour tenter de prolonger
x x
x [a, b], u ( x) K uv x [a, b], u ( x) K exp v
( I ) // HP //
Gronwall : u, v C ([a, b], ). Supp K
a a
x w' x
K 0 w K uv 0, C 1 , w ' uv v ln w( x) ln w(a) v K 0 limite K 0
a w a
II. Systèmes différentiels linéaires
I intervalle de , ou , E ev de dim n 1
A C ( I , L ( E )) B C ( I , E ) (ou matrices)
x '( t ) a11 ( t ) a1 n ( t ) x1 ( t ) b1 ( t )
1
( ) X ' AX B ( 0 ) X ' AX (càd : )
n
x '( t ) an 1 ( t ) an n ( t ) xn ( t ) bn ( t )
0 {sol de ( 0 )} est un sev de C 1 ( I , E). Si {sol de ( )} , c'est un sea de C 1 ( I , E ) dirigé par 0
m m
Si X k est sol de X ' AX Bk , k 1, m , k X k est sol de X ' AX k Bk
k 1 k 1
Fiches de maths - MP* - [Link] - 1
Chap 29 : Equations différentielles linéaires
Cauchy linéaire : A C ( I , L ( E )), B C ( I , E ), t0 I , X 0 E
Il existe une unique solution de ( ) X ' AX B définie sur I , telle que X (t0 ) Y0
On prend t , X 0 (t ) Y0 , X n1 (t ) Y0 ( A( s) X n ( s) B( s))ds , on étudie ( X n1 X n ) , majorations, récurrences,
t
// HP / /
t0
(t t0 ) n1
sup A(t ) , sup B(t ) X n1 (t ) X n (t ) n ( Y0 ) CVN sur segments CVU X sol
tJ tJ (n 1)!
t t
Unicité : [Link], X et X sol, Y X X sol/Y (t0 ) 0 Y '(t ) A(t ) Y (t ) , w Y ' 0 w A w + Gronwall
t0 t0
A C( , Mn ( )), B C ( , n
), X sol sur n
de X ' AX B, X (0) 0 t 0, X (t ) 0 (meth. itérative)
exp A ' Aexp A
t
t0
t
t0
Vrai si A commute avec A
t
t0
E 0 E
est une bijection affine, est un isomph de ev dim 0 dim E
X X (t0 ) X X ( t 0 )
( X1... X p ) fam. de sol. de X ' AX ( X1 ,..., X p ) est libre t0 I , ( X1 (t0 )... X p (t0 )) est libre t I ,.....
Une base ( X1... X n ) est un système fondamental de solutions de ( 0 )
Si on a ( X1 (t0 )... X n (t0 )) I n , R ( X1... X n ) est la résolvante d'indice t0 et vérifie : t I , R '(t ) A(t ) R(t )
Variation des constantes Trouver une sol part. de ( ) connaissant un système fondamental de ( 0 )
n
On cherche la solution sous la forme X (t ) k (t ) X k (t ),
k 1
det( X 1 (t ),..., X i 1 (t ), B(t ), X i 1 (t ),..., X n (t ))
Cramer : i est primitive de t
det( X 1 (t )... X n (t ))
X1 ,..., X n famille de solutions de ( 0 ) X ' AX , A C ( I , Mn ( ))
I
Le Wronskien de ( X 1... X n ) est l'application , noté w( X1... X n ) w
t det( X1 (t )... X n (t ))
n
u L( n
), x1...xn n
, det( x1...u ( xi )...xn ) tr(u )det( x1...xn ) (app. n-lin alt det , eval sur (can))
i 1
w est de classe C 1 et t I , w '(t ) tr( A(t ))w(t ) t0 I , w(t ) w(t0 )exp tr A(s)ds
t
t0
On a : soit w 0, soit w ne s'annule pas sur I
III. EDL scalaires d’ordre n supérieur à 2
J intervalle de , An ,..., A0 , B C ( J , ) ( ) An y ( n ) ... A1 y ' A0 y B
On se place sur I J où An ne s'annule pas (= hors des singularités)
ED ordinaire : ( ) y ( n ) an1 y ( n1) ... a0 y b
(I ) y0 ' y1 , y1 ' y2 ....
On attache à ( ) le système différentiel (S)
yn1 ' (a0 y0 ... an1 yn1 )b
Fiches de maths - MP* - [Link] - 2
Chap 29 : Equations différentielles linéaires
0 1 0 ,
y y
0 0
A( x) B( x) (S) : ( x) A( x) ( x) B( x)
0 1 0
a ( x) an 1 ( x ) b( x) yn 1 yn 1
0
y C (I , ) y est de classe C k , sol de ( )
y est de classe C 1 et il existe y1... yn C 1 ( I , ) tels que t ( y, y1 ,..., yn1 ) soit sol. de (S)
Cauchy : x0 I , ( z0 ...zn1 ) n
, il existe une unique solution de ( ) tq y( x0 ) z0 ,..., y ( n1) ( x0 ) zn1
S0
n
Si b 0, 0 {sol. de ( )}, x0 I , x0 ( n 1)
est un isomph de ev dim 0 n
y ( y( x0 ),..., y ( x0 ))
Pour a0 ...an1 constantes et b 0, on note 1 ...r les racines de multiplicité 1 ... r , de l'équation
caractéristique X n an1 X n1 ... a0 0 0 Vect ( x k ei )i 1,n
k 0 ,i 1
IV. EDL scalaires d’ordre 2
Sur J : Ay " By ' Cy D Sur I où A ne s'annule pas :
EDO : ( ) y " ay ' by c EDH : ( 0 ) y " ay ' by 0
,
y 0 1 y 0
(I )
Système associé : (S)
z b a z c
y C (I , ) y est C 2 , sol de ( ) y est C 1 , et z C1 ( I , ) tq ( y, z) sol de (S)
Cauchy : x0 I , ( y0 , y0 ') 2
, ! y tq y( x0 ) x0 et y '( x0 ) y0 '
0 2
x0 est un isomph de ev
y ( y ( x0 ), y '( x0 ))
, ( , ) est une base de , est un sys. fond. de sol de (S0 )
' '
0 0
/ x I ,
( x) ( x)
'( x) '( x)
0 x I , w( x) w( x0 )exp a x
x_0
(I )
On a : soit w 0, soit, x I , w( x) 0 Si ( 0 ) est de la forme y " qy 0, w est constant
Variation des constantes : on l'applique sur le système associé.
y ( x) ( x)
On cherche , C 1 ( I , ) tq ( x) ( x) ( x) soit sol.
z '( x) '( x)
( x) ( x) '( x) 0 '( x)
(I )
x I , On résout en , on intègre
'( x) '( x) '( x) c( x) '( x)
W ( x ) wronskienne
Chercher directement la forme ne mène à rien
f C 2 ( , ) tq f " f 0 x, f ( x) f ( x ) 0 (passer par f " f 0, VDC, expression intégrale...)
Fiches de maths - MP* - [Link] - 3
Chap 29 : Equations différentielles linéaires
Si on connait une solution particulière z de l'EDL Ay " By ' Cy D on cherche une solution
affinement indépendante de z sous la forme y uz EDL du premier ordre en u , quadrature
(I )
On peut utiliser le wronskien : si w yz ' y ' z , w( x) C exp a EDL d'ordre 1 en y
x
x0
Equation d'Euler : x2 y " axy ' by 0, a, b ctes, sur *
, on pose x et , z (t ) y(et ) EDL à coeffs csts
Compléments : méthodes usuelles
Variation des constantes : choix de constantes d'intégration solutions bornées... (en général, ()...)
x0
x et
f C ( , ) bornée y " y f possède une unique sol bornée sur (base : e x , e x , y Ae x f (t )dt ...)
0 2
y " ay ' by 0, a, b C ( I , ), , solutions à valeurs réelles.
( x0 ) 0 et 0 0 tq x [ x0 , x0 ] I , ( x x0 ( x) 0) (Cauchy : ' 0)
[a, b] I , {x [a, b] / ( x) 0} est fini. (sinon, suite inj CV de zéros) || ( x0 ) 0 ( , ) est liée
x1 x0 , ( x0 ) ( x1 ) 0 x [ x0 , x1 ] / ( x) 0 (abs, f / , f ' w / 2 strictement monotone...)
q1 , q2 C ( , ), q1 q2 , (E1 y " q1 y 0), (E2 y " q2 y 0), sol 0 de (E1 ) / (a) (b) 0, sol de (E2 )
s'annule sur [a, b] (wronskien mixte : w ' ', w ' (q2 q1 ) monotonie, signe)
Problème de Sturm-Liouville : q C ( , ), ( ) y " qy 0 a b
Existence et unicité de y tq (a) , (b) ? ( problème de Cauchy)
Si q 0, ( , ) 2
, ! / (a) et (b) 1 ,2 base, S0
1 ( a ) 2 ( a ) 0
( b ) ( b ) 0 , unique 0 de S0
2 2
Développement en série entière Franchissement des singularités
Analyse : Supp sol DSE, valeur des ak ? Synthèse : on prend cette sol, on vérifie la CV
Etude asymptotique : VDC, solution intégrale...
Série de Fourier : regarder les coeff 1 à 1 pour l'analyse (ne pas développer tout de suite)
V. SDL à coefficients constants
( ) X '(t ) AX (t ) B(t ), A Mn ( ), B C ( I , ) ( 0 ) X ' AX
Pour tout x0 n
, la solution X de ( ) tq X (0) x0 est : t , X (t ) exp(tA) x0
La solution M C ( I , Mn ( )) de M ' AM telle que M (0) I n est : t , M (t ) exp(tA)
A, B Mn ( ) [ A, B] 0 e A B e AeB (( ) M ' ( A B)M , et ( A B ) est la sol, etAetB aussi avec M (0) I n )
Au programme
(s, t ) 2
, e( st ) A esAetA
C ( , Mn ( )) tq : ( s, t ) 2
, ( s t ) ( s) (t ) et (0) I n (t ) et '(0)
(Intégrer, DL pour inversibilité, dériver, SDL)
Structure : Si V est vecteur propre de A pour la vp , t et V est sol de ( 0 )
n
i e i Vi ,(1 ... n )
t
Si A est diag dans la base (V1 ...Vn ), avec AVi i Vi , alors ( 0 ) n
i 1
t
B C( , n
) Une solution particulière de ( ) est t eta e sA B( s )ds
0
En général, trois possibilités : DL, diagonalisation ou calcul direct de etA
x ", y " On pose u x ', v y ' pour se ramener à l'ordre 1
Fiches de maths - MP* - [Link] - 4