0% ont trouvé ce document utile (0 vote)
58 vues11 pages

Intégration et Équations Différentielles

Ce document présente deux méthodes numériques pour résoudre des intégrales et équations différentielles: la méthode des trapèzes, de Simpson et d'Euler. Il décrit les formules mathématiques utilisées et propose de comparer les performances de ces méthodes sur différents exemples.

Transféré par

hupselle bitassi
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)
58 vues11 pages

Intégration et Équations Différentielles

Ce document présente deux méthodes numériques pour résoudre des intégrales et équations différentielles: la méthode des trapèzes, de Simpson et d'Euler. Il décrit les formules mathématiques utilisées et propose de comparer les performances de ces méthodes sur différents exemples.

Transféré par

hupselle bitassi
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

Année 2023–2024

LIE–OAN: Outils d’analyse numérique pour l’EEA

Travaux Pratiques
2

Partie I: Méthode d’intégration


Méthodes d’intégration 3

On cherche à calculer les intégrales


Z 4 Z 4
5x + 1 dx et x 3 + 2x + 1 dx
0 0

avec les deux méthodes:


Méthode des Trapèzes
Méthode de Simpson
Méthode des Trapèzes 4

" N−1
#
h X
Formule générale : S = f (a) + f (b) + 2 f (xi )
2 i=1

Définir de la fonction f (x) à intégrer


b−a
Définir le nombre d’intervalle, et calculer la taille d’intervalle h =
N
Boucle pour calculer la somme
Méthode de Simpson 5

( N−1 N−2
)
h X X
Formule générale : S = f (a) + f (b) + 4 f (xi ) + 2 f (xi )
3 i=1, impair i=2, pair

Boucle pour calculer la somme


Objective 6

Objective:
Tester les performances des deux méthodes en fonction de nombre
d’intervalle
deux différentes fonctions à intégrer
comparer avec la solution exacte
donner les conclusions
7

Partie II: Résolution des équations


différentielles
Résolution des équations différentielles 8

Objective: On cherche une solution approchée d’une équation


différentielle dans un intervalle [a, b]. Pour un système differentiel du 1er
ordre, on cherche une fonction y (t) définie et dérivable dans l’intervalle
tel que
dy (t)
= f (t, y (t)), ∀ t ∈ [a, b]
dt

y (a) = η

Euler explicite:
y (ti+1 ) = y (ti ) + ∆t · f (ti , y (ti ))
Euler implicite:

y (ti+1 ) = y (ti ) + ∆t · f (ti+1 , y (ti+1 ))

Runge-Kutta (d’ordre 2)

1
y (ti+1 ) = y (ti ) + (K1 + K2 )
2
avec K1 = ∆t · f (ti , y (ti )) et K2 = ∆t · f (ti+1 , y (ti ) + K1 )
Résolution des équations différentielles 9

R = 5 Ohms
L = 0.1 H

V (t) = 2Veff sin(ωt + φ)
Veff = 220 Volts
f = 50 Hz
Utiliser les trois méthodes pour calculer le courant i(t) entre 0 et 70 ms avec
un pas de temps de 1 ms
Résolution des équations différentielles 10

ODE à résoudre:

di(t) 1
= (V (t) − Ri(t))
dt L

= f (t, i(t))
Résolution des équations différentielles 11

Definition de la function V (t)


Definition de la function f (t, i(t))
Definition: pas de temps (∆t), temps de simulation (T = 70 ms),
T
nombre d’intervalle (n= ∆t ),
Boucle pour les méthodes d’Euler

Vous aimerez peut-être aussi