Résolution de l’exercice 1
(i) Equation de Bessel :
x 2 y ' ' + x y ' + ( x 2−υ 2) y=0 ; υ>0 (1)
(ii) Forme générale dans la méthode de Frobenius :
x 2 ( α 0 + α 1 x+ α 2 x 2 ) y +x( {β} rsub {0} + {β} rsub {1} x+ {β} rsub {2} {x} ^ {2} )y'+( {γ} rsub {0}
Par identification (1) = (2), on obtient :
α 0=1 ; α 1=0 ; α 2=0
β 0=1 ; β 1=0 ; β 2=0
2
γ 0=−υ ; γ 1=0 ; γ 2=1
(iii) Polynômes indiciels :
pi ( r )=α i r ( r −1 ) + β i r + γ i ; i=1,2,3
{
2 2 2
p 0 ( r )=r ( r−1 ) +r −υ =r −υ
p1 ( r )=0
p2 ( r )=1
(iv) Domaine de convergence :
I conv =¿ 0 ,+∞ ¿
(v) Equation indicielle :
2 2
p0 ( r )=0 ⇔ r −υ =0
⇔ 1
{ r =υ> 0
r 2 =−υ< 0
(vi) Coefficients a n :
a 0=1
− p1 ( r )
a 1= ⇒ a1=0
p0 ( r +1 )
−1
a n=
p 0 ( n+ r )
[ p1 ( n+r−1 ) a n−1 + p2 ( n+r −2 ) an−2 ] ; n ≥ 2
−1
⇒ an= an−2
( n+r )2−υ2
−1 −1
⇒ an= a =
2 n−2
a
2
( n+r ) −υ ( n+r +υ ) ( n+r−υ ) n−2
{
−1
a n= a ; pour r=υ
⇒ n(n+ 2υ) n−2
−1
an = an−2 ; pour r =−υ
n(n−2υ )
(vii) Généralisation :
{ {
−1 −1 −1
a2=
a0 = 2 a3 = a1=0
2(2+2 υ) 2 (1+υ) 2(2+2 υ)
n=2⇒ ; n=3 ⇒
−1 −1 −1
a2= a0 = 2 a3 = a =0
2(2−2 υ) 2 (1−υ) 2(2−2 υ) 1
{ {
1 1 −1
a4 =
a2 = 4 a5= a3=0
4 (4 +2υ ) 2 ∗2∗(2+υ )(1+ υ) 5(5+ 2υ)
n=4 ⇒ ; n=5 ⇒
1 1 −1
a4 = a2 = 4 a 5= a =0
2∗2∗2(4−2υ)(2−2 υ) 2 ∗2∗(2−υ)( 1−υ) 5(5−2υ ) 3
{ {
−1 −1 −1
a 6= a4 = 6 a7 = a 5=0
6(6+ 2υ) 2 ∗2∗3∗(3+υ)( 2+ υ)(1+υ) 7 (7+2 υ)
n=6 ⇒ ; n=7 ⇒
−1 −1 −1
a6 = a= a7= a =0
6 (6−2υ ) 4 26∗2∗3(3−υ)(2−υ)(1−υ ) 7(7−2 υ) 5
…..
n=2 n⇒ ¿; n=2 n+1⇒ a2 n+1=0
(viii) Système fondamental des solutions (SFS) :
r 1=υ et r 2=−υ ⇒r 1 −r 2=2 υ ∈ N , car υ ∈ N . Nous avons :
n
SFS={ y 1 } ⇒ y 1= ∑ ¿ ¿ ¿
k=0
Résolution de l’exercice 2
(i) Expression :
2
Ly=x y +2x {y} ^ {'} -3xy (1
(ii) Forme de la solution :
∞
y= y ( x )=∑ an (x−2)n (2)
n=0
Posons z=x−2 ⇒ x=z +2
Les relations (1) et (2) deviennent :
( z +2 )2 y +2(z+2) {y} ^ {'} -3(z+2)y (1'
∞
y=∑ an z (2')
n
n=0
(iii) Domaine de convergence
I conv =¿ 0 ; 4 ¿
(iv) Résolution :
∞
y ' =∑ an n z n−1 (3)
n=0
y = sum from {n=0} to {∞} {{a} rsub {n} n left (n-1 right ) {z} ^ {n-2}} (4
(2’), (3) et (4) dans (1’) :
∞ ∞ ∞
Ly=( z+ 2 )
2
∑ a n n ( n−1 ) z n−2 +2( z +2)∑ an n z n −1−3(z+ 2) ∑ a n z n
n=0 n=0 n=0
∞ ∞ ∞
¿(z + 4 z+ 4) ∑ a n n ( n−1 ) z +2(z +2) ∑ an n z −3 ( z+ 2)∑ a n z
2 n−2 n−1 n
n=0 n=0 n=0
∞ ∞ ∞ ∞ ∞ ∞ ∞
¿ ∑ a n n ( n−1 ) z n + ∑ 4 a n n ( n−1 ) z n−1 + ∑ 4 a n n ( n−1 ) z n−2 + ∑ 2 an n z n+ ∑ 4 an n z n−1 −∑ 3 an z n+1−∑ −6 an
n=0 n=0 n=0 n=0 n=0 n=0 n=0
∞ ∞ ∞ ∞
¿ ∑ 4 an n ( n−1 ) z n−2 + ∑ 4 an n2 z n−1 + ∑ an [n ( n−1 ) −6] z n−∑ 3 a n z n +1
n=0 n=0 n=0 n=0
∞ ∞ ∞ ∞
¿ 0 a 0+ 0 a1 +8 a2 + ∑ a n n ( n−1 ) z +0 a 0+ 4 a1+ ∑ 4 an n z −6 a0 + ∑ a n [ n ( n−1 )−6 ] z −∑ 3 an z
n −2 2 n−1 n n +1
n=3 n=2 n=1 n=0
∞ ∞ ∞ ∞
¿ 8 a2 + 4 a1 −6 a0 + ∑ an n ( n−1 ) z +∑ 4 a n n z + ∑ an [ n ( n−1 )−6 ] z −∑ 3 a n z
n−2 2 n−1 n n+1
n=3 n=2 n=1 n=0
∞ ∞ ∞ ∞
¿ 8 a2 + 4 a1 −6 a0 + ∑ an+2 (n+ 2) ( n+1 ) z + ∑ 4 an+1 (n+ 1) z + ∑ an [ n ( n−1 )−6 ] z −∑ 3 an−1 zn
n 2 n n
n=1 n=1 n =1 n=1
∞
¿ 8 a2 + 4 a1 −6 a0 + ∑ [ an +2 ( n+2 ) ( n+1 ) + 4 an +1 ( n+ 1 ) + an (n2−n−6)−3 an−1 ] z n
2
n=1
D’où :
∞
Ly=8 a2 +4 a1−6 a0 + ∑ [ an+2 ( n+2 ) ( n+1 ) +4 an +1 ( n+1 ) +a n (n2−n−6)−3 an−1 ] (x−2)n
2
n=1
Résolution de l’exercice 3
Equation :
9 x 2 y ' ' + 3 x (5+3 x−2 x 2) y ' + ( 1+ 12 x−14 x 2 ) y=0 (1)
(i) Forme générale dans la méthode de Frobenius :
x 2 ( α 0 + α 1 x+ α 2 x 2 ) y +x( {β} rsub {0} + {β} rsub {1} x+ {β} rsub {2} {x} ^ {2} )y'+( {γ} rsub {0}
Par identification (1) = (2), on obtient :
α 0=9 ; α 1=0 ; α 2=0
β 0=15 ; β 1=9 ; β 2=−6
γ 0=1 ; γ 1=12 ; γ 2=−14
(ii) Polynômes indiciels :
pi ( r )=α i r ( r −1 ) + β i r + γ i ; i=1,2,3
{
p 0 ( r )=9 r ( r−1 ) +15 r +1=(3 r +1)2
p 1 ( r )=9 r +12=3 (3 r +4)
p 2 ( r ) =−6 r−14=−2(3 r+ 7)
(iii) Domaine de convergence :
I conv =¿ 0 ,+∞ ¿
(iv) Equation indicielle :
2
p0 ( r )=0 ⇔(3 r +1) =0
−1
r 1=r 2=r =
3
(v) Forme des solutions :
∞ ∞
y 1=x r
∑ a n (r )x n
et y 2= y 1 ln ( x )+ x
r
∑ a' n (r ) x n
n=0 n=0
(vi) Coefficients a n :
a 0=1
−p 1 ( r ) −3 (3 r + 4 ) −3
a 1 ( r )= = =
p0 ( r +1 ) ( 3 r + 4 )2 3r+4
−1
a n (r )= [ p ( n+r−1 ) a n−1 + p2 ( n+r −2 ) an−2 ]
p 0 ( n+ r ) 1
1
¿
3 n+3 r +1
[−3 a n−1 ( r ) +2 an−2 ( r ) ] ; n≥ 2
9
a '1(r )=
( 3 r + 4 )2
1 1
'
a n ( r )=
3 n+3 r +1
[ −3 a' n−1 ( r )+ 2a ' n−2 ( r ) ]+ 2[
9 an−1 ( r )−6 an −2 ( r ) ] ; n ≥2
(3 n+ 3 r+ 1)
On a :
{
a 0=1
⇒ a1 (−1/ 3 )=−1
1
an (−1/ 3 )=
3n
[−3 an−1 (−1 /3 )+ 2a n−2 (−1/3 ) ]
{
a ' 1 (−1/ 3 )=1
Et 1 1
a' n (−1/3 )=
3n
[ −3 a ' n−1 (−1/3 ) +2 a ' n−2 (−1/3 ) ] + 2 [ 3 an−1 (−1/3 )−2 an−2 (−1/3 ) ]
3n
(vii) Programme de calcul de récurrence (sur MATLAB R2013a)
(viii) Valeurs obtenues après calcul
(ix) Solutions
−1 /3 5 2 1 3
y 1=x (1−x+ x − x +…)
6 2
11 25
y 2= y 1 ln ( x )+ x 2 /3 (1− x + x 2 +…)
12 36