Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Analyse Numerique Google

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

Pour n+1 points de collocations, on aura

Et la formule de Lagrange

Formule pour 2 points:

Formule pour 3 points:

U. Laval Dept. Math & Stat MAT-2910


1
Chapitre 5
Méthode de Newton

Comment déterminer les bn?

Et par récurrence, on peut déterminer tous les coefficients

U. Laval Dept. Math & Stat MAT-2910


2
Chapitre 5
Différences divisées

On introduit la notation

On appellera première différence divisée, notée le rapport

On appellera deuxième différence divisée, notée le rapport

Et plus généralement la pième différence divisée, notée

U. Laval Dept. Math & Stat MAT-2910


3
Chapitre 5
Revenant aux coefficients du polynôme on a

Pour (x0,y0), (x1,y1)

Pour 3 points:

Et en général, pour n points:

Comment calculer de manière efficace les différences divisées?

U. Laval Dept. Math & Stat MAT-2910


4
Chapitre 5
Tableau des différences divisées
On exploite la structure particulière des diff. divisées. Soit (xi,f(xi)) i = 0,…,n, on
forme une table dont chaque colonne correspond à toutes les différences divisées d’un
même ordre. En augmentant l’ordre de la gauche vers la droite, chaque colonne se
construit en fonction de la colonne précédente. Pour 4 points,

U. Laval Dept. Math & Stat MAT-2910


5
Chapitre 5
Remarques:
• Presque sous la forme Horner, donc facilement transformable pour ramener à Horner.
• On peut facilement enlever/ajouter des points de collocations.
• Dans le cas d’une fonction correspondant à un polynôme de degré n, alors les différences
divisées d’ordre supérieur à n seront toutes nulles (unicité du polynôme).
• Si on n’utilise pas tous les points alors leur ordre d’apparition aura de l’importance: on
n’aura pas le même polynôme de collocation puisqu’on ne prendra pas les mêmes points!

Erreur d’interpolation

Soit (xi,f(xi)) i = 0,…,n. On défini l’erreur d’interpolation de f par

Où pn(x) est l’interpolant de f aux nœuds x0,x1,…,xn.

On s’intéresse à l’erreur « entre » les nœuds, on voudrait évaluer l’erreur pour .

Pour cela on utilisera Taylor et le théorème de Rolle.

U. Laval Dept. Math & Stat MAT-2910


6
Chapitre5
Théorème (sur l’erreur d’interpolation)

Soit x0 < x1 < x2 …< xn, les abscisses des points d’interpolation. Si f est définie sur [x0,xn]
et n+1 fois dérivable sur ]x0,xn[ alors pour tout , il existe tel que

Remarques
• Analogue au terme d’erreur de Lagrange, est inconnu et dépend de x. Pour estimer
l’erreur, on le fera disparaitre en faisant une majoration de la dérivée de f :

• En(xi) = 0 i=0,1,…,n: En est un polynôme de degré n+1 dont les racines


sont x0,…,xn.
• Plus il y aura de points de collocation plus l’erreur oscillera autour de l’axe des x.
L’augmentation du nombre de points pouvant mener à une augmentation importante de
l’erreur (effet de Runge)

U. Laval Dept. Math & Stat MAT-2910


7
Chapitre 5
• Si on peut choisir les points de collocations, en se basant sur l’erreur d’interpolation, il
est préférable de choisir les points le plus près des points d’intérêts.

• Si on ne connait pas f, on ne peut pas en calculer la dérivée, cependant puisque

Alors dans le cas où varie peu on aura

Dans ce cas:

Ce qui revient à estimer l’erreur depn(x) en prenant comme valeur de référence un


polynôme de degré n+1.

Attention, cet estimation n’est pas toujours très fiable.

U. Laval Dept. Math & Stat MAT-2910


8
Chapitre 5

Vous aimerez peut-être aussi