Scribd
Importer
6 vues3 pages

Schéma aux différences finies en Python

Transféré par

Xan Dutheil
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
6 vues3 pages

Schéma aux différences finies en Python

Transféré par

Xan Dutheil
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

Projet Python

[Link]
Décembre 2022 -Janvier 2023

1 Construction du shémas aux différences finies


Soit N ∈ N , on pose h = 1/N .
On pose ∀i ∈ [0, N ] : xi = ih de telle sorte que x0 = 0, xN = 1.
On pose ∀i ∈ [0, N ], ui = u(xi )
On a alors grâce aux formulles de Taylor:

h2 ′′ h3 (3)
ui+1 = ui + hu′i + u + u + O(h4 ) (1)
2 i 6 i
h2 h3 (3)
ui−1 = ui − hu′i + u′′i − u + O(h4 ) (2)
2 6 i
Il vient alors en faisant (1) + (2):
ui−1 − 2ui + ui+1
u′′i = + O(h2 ) (3)
h2
En notant f (x) = 23 ex le problème continue 5.3 devient le problème discret
suivant:  ui−1 −2ui +ui+1
 h2 − ui = fi ∀i ∈ [1, N − 1]
u0 = 0 (4)
uN = e/3

On note alors:
 
1 0 0 ··· ··· 0
1 .. 
 2
h − h22 − 1 1
h2 0 ··· . 
1 2 1

0
h2 − h2 − 1 h2 0 · · ·
A=
 .

 .. .. .. .. 
 0 . . . 

1
 0 h2 − h22 − 1 1 
h2
··· 0 1

1
et  
0
  f1
..
 
 



 .
B= f
 i 

 . 
 .. 
 
fN −1 
e/3
Le problème discret (4) devient alors:

A(ui )i∈[0,N ] = B

2 Calcul de l’ordre de convergence


2.1 Solution exact
Par des calculs classiques de résolution d’équation differentiels, on obtient que
x 7→ 13 xex est la solution exact au problème (5.3).

3 Implémentation
3.1 Simulation
On peut voir ici une simulation pour n = 10 et la solution exacte.

2
On peut voir ici une simulation pour n = 10 et la solution exacte.

3.2 Convergence
On peut voir ici l’erreur en fonction de n. En echelle log-log.

Vous aimerez peut-être aussi