4.l'integration Numérique

Télécharger au format pptx, pdf ou txt
Télécharger au format pptx, pdf ou txt
Vous êtes sur la page 1sur 15

Chapitre 9

L’intégration numérique
PRÉSENTÉ PAR DR, MEHIAOUI M,Y,

1
9,1, La méthode des trapèzes
F(x)
Fonction à intégrer

Trapèzes

x
Intervalle d’intégration

9,1, La méthode des trapèzes 2


 Pour une intégrale bornée

La fonction étant connue en points , dans l’intervalle


Si l’échantillonage est régulier c à d si

F(x)
Fonction à intégrer

Trapèzes

x
Intervalle d’intégration

9,1, La méthode des trapèzes 3


z 1 3 5 7 9 11
F(z) 0,8 0,9 1,2 0,7 1,0 0,8

Exemple
 Soit une fonction F(z) connue  L’échantillonage est constant, soit , d’où
en des points déterminés,
Donner une approximation de
son l’intégral sur le domaine
[1 , 11] en utilisant la méthode
des trapèzes

9,1, La méthode des trapèzes 4


9,2, La méthode de quadrature de Gauss
Remplaçons
  l’intégrale bornée d’une fonction polynomiale par une combinaison linéaire de ses
valeurs aux points d’intégrations

Déterminons les coefficients et les points d’intégration

En remplaçant l’équation (4) dans l’intégrale (3), celle-ci devient

9,2, La méthode de Gauss 5


2 3 2 𝑟 −1
 𝑃 ( 𝜉 )=𝑎 1+ 𝑎2 𝜉❑
❑ + 𝑎 𝜉
3 ❑ + 𝑎 𝜉
4 ❑ +, , , ,++ 𝑎 𝜉
2𝑟 ❑

+1
  𝑟

∫ 𝑃(𝜉 )𝑑 𝜉 =𝜔1 𝑃 ( 𝜉 1 )+𝜔2 𝑃 ( 𝜉 2 )+, , ,, , ,+𝜔𝑟 𝑃 ( 𝜉 𝑟 ) =∑ 𝜔𝑖 𝑃 ( 𝜉 𝑖 )


−1 𝑖=1
 

5)
D’où

(+++,,,+

(6)

9,2, La méthode de Gauss 6


 L’équation (6) doit être vérifiée pour n’importe quel polynôme, i,e,

(7)
La solution devient

Soit
(8)

(9)

 Sous les conditions

9,2, La méthode de Gauss 7


 
La méthode peut intégrer un polynôme d’ordre soit si

Et donc
Exemple Pour que cette relation soit satisfaite, il faut que
 Calculdes r coefficients et Éqt (9)
positions de la méthode de
On peut démontrer ce résultats en remplaçant les
Gauss, pour intégrer le
expressions de et , dans la relation
polynôme entre
-1 et +1 La solution du système est
et

9,2, La méthode de Gauss 8


Position des points d’intégration et coefficients de
pondération
9,2, La méthode de Gauss 9
 Pour les intégrales de surface, une double sommation est utilisée, ainsi

Pour les intégrales de volumes, l’intégrale triple est transformée en une triple sommation

Remarque:
En pratique , le nombre de points de gauss est réduit afin de gagner du temps machine, Par expérience, la méthode
converge,

9,2, La méthode de Gauss 10


Quadrature de
Gauss pour une
intégrale 2D
02 points dans chaque
direction

9,2, La méthode de Gauss 11


 
 Solution exacte:

 Méthode de Gauss: Il faut d’abord faire un changement de


Exemple base entre et
d’application y 
Evaluer l’intégrale suivante
avec la méthode de gauss et
la méthode des trapèzes, 0 +6 x -1 +1 
Comparer avec la solution
exacte

La transformation géométrique est

9,2, La méthode de Gauss 12


 1 pnt de Gauss erreur = -70 %

2 pnts de Gauss erreur = -12,7 %

3 pnts de Gauss erreur = -0,9 %

4 pnts de Gauss 402,286901 erreur = -0,09 %

5 pnts de Gauss 402,425336 erreur = -0,001 %

A partir de 3 pnts, l’erreur devient inférieure à 1 %

• Méthode des trapèzes: En utilisant 05 pnts d’intégration

erreur = 18,1 %

9,2, La méthode de Gauss 13


9,3, La méthode de Hammer
 𝜂

Méthode de Méthode  ❑ +1 1 −𝜉
Eléments
quadrature de
triangulaires
Gauss inadaptée
de
Hammer
+1 ∫ 𝑑𝑉 =∫ ∫ 𝑑 𝜉 𝑑 𝜂
𝑉𝑒 0 0

dépendance entre
les bornes
Solution
d’intégrations de
et de  𝜉
0 +1

9,3, La méthode de Hammer 14


 La méthode permet d’intégrer un polynôme d’ordre m dont le terme le plus élevé est tel que

Ordre m

1 1 1/3 1/3 1/2

1/6 1/6
2 3 2/3 1/6 1/6
1/6 2/3

1/3 1/3 -27/96


1/5 1/5 25/96
3 4
3/5 1/5 25/96
1/5 3/5 25/96

9,3, La méthode de Hammer 15

Vous aimerez peut-être aussi