09.01.2015 Views

ANALYSE NUMERIQUE Chapitre 3 Intégration numérique ... - lmpt

ANALYSE NUMERIQUE Chapitre 3 Intégration numérique ... - lmpt

ANALYSE NUMERIQUE Chapitre 3 Intégration numérique ... - lmpt

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Polynôme d’interpolation de Lagrange Lorsque f est connue au moins en<br />

n + 1 points(distincts) de subdivision et est de classe C n+1 , on écrit facilement<br />

le polynôme d’interpolation de Lagrange P n de f et l’erreur e,<br />

i=n<br />

∑<br />

1<br />

f(x) = P n (x) + e(x) = L i (x)f(x i ) +<br />

(n + 1)! L(x)f (n+1) (ζ x )<br />

avec L i (x) =<br />

j=n ∏<br />

j = 0<br />

j ≠ i<br />

i=0<br />

x − x j<br />

x i − x j<br />

et L(x) = ∏<br />

0≤i≤n<br />

Ainsi l’intégrale du polynôme de Lagrange donne :<br />

∫ b<br />

a<br />

i=n<br />

∑<br />

P n (x)dx = (<br />

i=0<br />

∫ b<br />

a<br />

L i (x)dx)f(x i ) = J n (f)<br />

(x − x i ) oú ζ x ∈ [a, b]<br />

Définition .2 La formule J n (f) est une formule de quadrature de type interpolation<br />

pour (1)<br />

On mesure la précision des quadratures numériques à l’aide des notions :<br />

Définition .3<br />

• L’erreur d’intégration est le réel E(f) = I − J n (f)<br />

• Une formule de quadrature est dite exacte sur F s.e.v. de C 0 [a, b], si<br />

∀g ∈ F E(g) = 0<br />

• Une formule de quadrature a un degré de précision p ∈ IN si elle<br />

est exacte sur l’espace vectoriel des polynômes IP p (de degré ≤ p), avec<br />

E(x p+1 ) ≠ 0.<br />

Ces définitions sont reliées entre elles !<br />

Proposition .1 L’erreur d’intégration, pour une formule de type interpolation<br />

peut se calculer en<br />

E(f) = f (n+1) ∫<br />

(θ) b<br />

L(x)dx<br />

(n + 1)! a<br />

pour un certain θ, lorsque le polynôme L reste de signe constant sur<br />

[a, b].<br />

Théorème .1 Une formule de quadrature à n + 1 points est exacte sur IP n si<br />

et seulement si elle est de type interpolation.<br />

REMARQUE : Les énoncés ci-dessus peuvent aussi s’écrire avec un ”poids<br />

d’intégration” µ ≥ 0 sur [a, b], pour des valeurs approchées de<br />

∫ b<br />

a<br />

µ(x)f(x)dx.<br />

2

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!