0% ont trouvé ce document utile (0 vote)
87 vues7 pages

Résolution d'équations de Bessel et Frobenius

Ce document contient la résolution de plusieurs exercices mathématiques. L'exercice 1 résout une équation de Bessel. L'exercice 2 résout une équation différentielle du second ordre. L'exercice 3 présente une autre équation différentielle.

Transféré par

Nseye Davyd
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 DOCX, PDF, TXT ou lisez en ligne sur Scribd
0% ont trouvé ce document utile (0 vote)
87 vues7 pages

Résolution d'équations de Bessel et Frobenius

Ce document contient la résolution de plusieurs exercices mathématiques. L'exercice 1 résout une équation de Bessel. L'exercice 2 résout une équation différentielle du second ordre. L'exercice 3 présente une autre équation différentielle.

Transféré par

Nseye Davyd
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 DOCX, PDF, TXT ou lisez en ligne sur Scribd

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

Vous aimerez peut-être aussi