0% ont trouvé ce document utile (0 vote)
4 vues3 pages

Examen Final (MNAO) : D U X DX U X e U U

Le document présente un examen final en mathématiques appliquées, axé sur la résolution d'un problème différentiel avec différentes méthodes numériques. Les étudiants doivent appliquer la méthode des différences finies, la méthode de Gauss-Seidel et la méthode de Runge-Kutta pour trouver des solutions approchées. Il est également demandé de comparer les méthodes et de justifier celle qui est la plus précise.

Transféré par

ouariamina21
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)
4 vues3 pages

Examen Final (MNAO) : D U X DX U X e U U

Le document présente un examen final en mathématiques appliquées, axé sur la résolution d'un problème différentiel avec différentes méthodes numériques. Les étudiants doivent appliquer la méthode des différences finies, la méthode de Gauss-Seidel et la méthode de Runge-Kutta pour trouver des solutions approchées. Il est également demandé de comparer les méthodes et de justifier celle qui est la plus précise.

Transféré par

ouariamina21
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

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)
 e2x  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)  xe2x  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)  4e0.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)  4e0.15  8 (0.5 pt)
On trouve :
 2  u(0.05)  2u(0.025)  0.0025u(0.025)  0.0025e0.05  0.005...(0.5 pt)

 u(0.025)  u(0.075)  2u(0.05)  0.0025u(0.05)  0.0025e0.1  0.005..(0.5 pt)
u(0.05)  2.081873075  2u(0.075)  0.0025u(0.075)  0.0025e0.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.0251 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 e2*0.025  2(u(0.025))  4  0.022591198

K 2  0.025 e2*0.0375  2(u(0.025)  K1/ 2)  4 
 

0.02143977 (15x0.5pt)
K3  0.025 e2*0.0375  2(u(0.025)  K 2 / 2)  4 
  0.021468556
 
K 4  0.025 e2*0.05  2(u(0.025)  K3)  4  0.02035847
  

 u(0.05)  2.045241786

 K1  0.025e2*0.05  2(u(0.05))  4  0.020358846


K 2  0.0252*0.0625
e2*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 e2*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).

Vous aimerez peut-être aussi