UniversitéA.
MiradeBEJAIA Durée:1H30mn
FacultédeTechnologie,DépartementGEM
odule UEF 1122
Année2020/2021
Examen Final (MNAO)
On cherche à résoudre le problème différentiel suivant :
d 2u(x) 2 x
4u(x) 4e 8
dx
2
u(0) 2.................................(*)
u(0.1) 2.081873075
1. Donner le type des conditions aux limites.
2. Retrouver par la méthode des différences finies le schéma numérique (**) centré en
prenant un pas de discrétisation h=0.025 :
1
u(0.025) 2.0025 u(0.05) 1.997378
1 u(0.025) u(0.075) 0.0027379......(**)
u(0.05)
2.0025
1
u(0.075) 2.0025 u(0.05) 2.076176615
3. Donner les trois premières itérations des valeurs de X1=u(0.025), X2=u(0.05) et
X3=u(0.075) solution de l’équation (**) avec la méthode de Gauss-Seidel relaxée pour
un terme w=0.5 et un vecteur initial X 1(0) X2 (0) X
3 (0) 2 .
4. La même fonction u(x) est définie avec sa première dérivée dans le système
différentiel (***) suivant :
du(x)
e2x 2u(x) 4
........(***)
dx
u(0) 2
(Interrogation) Retrouver les valeurs précédentes u(0.025), u(0.05) et u(0.075) avec
la méthode de Runge-Kutta d’ordre 4
5. La fonction u(x) est déterminée analytiquement par u(x) xe2x 2 , quelle est la
meilleure méthode parmi ce que vous avez utilisé dans cette épreuve ? justifier !
NB : La précision est recommandée !
La question 4 est considérée comme une 2ème interrogation
Bonne réussite
Université[Link] Durée:1H30mn
FacultédeTechnologie,DépartementGEM
odule UEF 1122
Année2020/2021
Solution
1. Les conditions aux limites sont de type : Conditions de Dirichlet. (1pt)
2.
U(0.1)=2.081873075
U(0)=2 (1) (2) (3)
0 0.025 0.05 0.075 0.1
Schéma (1pt)
La discrétisation de l’équation (*) donne :
u(0) u(0.05) 2u(0.025) (0.025)2 4u(0.025) 4e0.05 8 (0.5 pt)
u(0.025) u(0.075) 2u(0.05) (0.025) 4u(0.05) 4e
2
0.1
8..(0.5 pt)
u(0.05) u(0.1) 2u(0.075) (0.025)2 4u(0.075) 4e0.15 8 (0.5 pt)
On trouve :
2 u(0.05) 2u(0.025) 0.0025u(0.025) 0.0025e0.05 0.005...(0.5 pt)
u(0.025) u(0.075) 2u(0.05) 0.0025u(0.05) 0.0025e0.1 0.005..(0.5 pt)
u(0.05) 2.081873075 2u(0.075) 0.0025u(0.075) 0.0025e0.15 0.005...(0.5 pt)
On aura le système d’équations suivant :
1 u(0.05) 1.997378074
u(0.025)
2.0025
u(0.05) 2.0025 u(0.025) u(0.075) 0.002737906455
1
1 u(0.05) 2.076176615
u(0.075)
2.0025
3.
En remplaçant par x1, x2 , x3 , on aura :
1.997378074
1
x 1 2.0025
x2
x x 0.002737906455
1
x
2
2.0025 1 3
x3
1
2.076176615
2.0025 x2
Les trois premières itérations sont données comme suit :
Application de la méthode (1pt)
Université[Link] Durée:1H30mn
FacultédeTechnologie,DépartementGEM
odule UEF 1122
Année2020/2021
x1(1) 1.99809688 x1(2) 1.99654426 x1(3) 1.996441124
x2 (1) 1.997592754 x2 (2) 2.000288821 2
x (3) (9x0.5pt)
2.003925365
x (1) 2.01717088 x (2) 2.0264295 x (3) 2.031966805
3 3 3
4.
En appliquant la méthode de Runge-Kutta, on trouvera :
K1 0.0251 2* 2 4 0.025
K 2 0.025 e 2(2 K1/ 2) 4 0.023757747
2*0.0125
2*0.0125
K3 0.025 e 2(2 K 2 / 2) 4 0.023788804
2*0.025
K 4 0.025 e 2(2 K3) 4 0.022591295
u(0.025) 2.023780733
K1 0.025 e2*0.025 2(u(0.025)) 4 0.022591198
K 2 0.025 e2*0.0375 2(u(0.025) K1/ 2) 4
0.02143977 (15x0.5pt)
K3 0.025 e2*0.0375 2(u(0.025) K 2 / 2) 4
0.021468556
K 4 0.025 e2*0.05 2(u(0.025) K3) 4 0.02035847
u(0.05) 2.045241786
K1 0.025e2*0.05 2(u(0.05)) 4 0.020358846
K 2 0.0252*0.0625
e2*0.0625 2(u(0.05) K1/ 2) 4 0.019291362
K3 0.025 e 2(u(0.05) K 2 / 2) 4
0.019318049
K 4 0.025 e2*0.075 2(u(0.05) K3) 4 0.018289707
u(0.075) 2.064553015
5.
La meilleure méthode dans ce cas est R-K-4 (vous pouvez vérifier en remplaçant dans
la fonction donnée dans cette question) (1pt) car la méthode de Runge Kutta est
d’ordre de précision (ordre 4) plus élevé que celui de la méthode des différences finies
(ordre 2).(1pt).