Approximation de Fonction (Cours Math 4 Usthb)
Approximation de Fonction (Cours Math 4 Usthb)
Approximation de Fonction (Cours Math 4 Usthb)
Problèmes Pratiques
Approximation de fonctions
Introduction 1/2
1
Introduction 2/2
2
Approximation de fonctions 2/2
-1
-2
-3
-4
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8
❚ Interpolation d'Hermite
❙ informations sur les dérivées de la fonction à approcher
❚ …
Ph. Leray Analyse Numérique 5
3
Interpolation polynomiale : Lagrange
❚ Théorème
❙ Soient n+1 points distincts xi réels et n+1 réels yi,
il existe un unique polynôme p ∈ Pn tel que
p(xi) = yi pour i = 0 à n
❚ Démonstration
n
par construction
❙ p( x ) = ∑ y i Li ( x )
i =0
n (x − x )
avec Li polynôme de Lagrange Li ( x ) = ∏
j
j =0 (x − x )
i j
j ≠i
❙ Propriétés de Li
❘ Li polynôme de degré n
❘ Li(xi)=1 Li(xj)=0 (j ≠ i)
Ph. Leray Analyse Numérique 6
4
Lagrange : exemple n°2 1/2
L0 ( x ) =
(x − 2 )(x − 4 ) L1 ( x ) =
x(x − 4 )
L2 ( x ) =
x (x − 2 )
8 −4 8
0 2 4 0 2 4 0 2 4
25
15
10
0
-1 0 1 2 3 4 5
5
Lagrange : exemple n°3 1/2
50
40
30
20
10
-10
-1 0 1 2 3 4 5
60
50
40
30
20
10
-10
-20
-1 0 1 2 3 4 5
6
Lagrange : erreur d'interpolation
❚ Théorème :
❙ si f est n+1 dérivable sur [a,b], ∀ x ∈ [a,b], notons :
❘ I le plus petit intervalle fermé contenant x et les xi
❘ φ(x)=(x-x0)(x-x1)…(x-xn)
(ξ )φ (x )
( n +1 )
❚ Utilité
❙ contrôle de l’erreur d’approximation
❙ « réglage » de la qualité de l’approximation
Ph. Leray Analyse Numérique 12
Lagrange : choix de n
❚ Supposons que l'on possède un nb élevé de points
pour approcher f … faut-il tous les utiliser ?
❙ (calculs lourds)
❚ Méthode de Neville :
❙ on augmente progressivement n
❙ on calcule les Li de manière récursive
❙ on arrête dès que l'erreur est inférieure à un seuil
(d’où l’utilité du calcul de l’erreur)
7
Lagrange : l’algorithme 1/2
p( x ) =
(x − x j )Q j (x )− (x − xk )Q k (x )
❘ on a
(xk − x j )
❘ avec Qj(x) (resp. Qk(x)) polynôme d’interpolation obtenu avec
tous les points sauf xj (resp. xk)
❚ Interpolation en 1 point x :
Pour i allant de 1 à n
Pour j allant de 1 à i
Q ij ( x ) =
(x − xi− j )Q i , j −1 (x ) − (x − xi )Q i−1, j −1 (x )
(xi − xi− j )
FinPour
Approximation de degré i = Qii(x)
FinPour
❚ Complexité du calcul : n2
Ph. Leray Analyse Numérique 15
8
Interpolation polynomiale : Newton
❚ Polynômes de Newton :
❙ base = {1, (x-x0), (x-x0)(x-x1), …, (x-x0)(x-x1)…(x-xn-1)}
❚ Définition :
❙ Soit une fonction f dont on connaît les valeurs en des
points distincts a, b, c, …
❙ On appelle différence divisée d’ordre 0, 1, 2,...,n
les expressions définies par récurrence sur l’ordre k :
❙ k=0 f [a] = f (a)
❙ k=1 f [a,b] = ( f [b] - f [a] ) / ( b - a )
❙ k=2 f [a,b,c] = ( f [a,c] - f [a,b] ) / ( c - b )
9
Newton : différences divisées 2/3
❚ Théorèmes :
❙ détermination des coefficients de p(x) dans la base de
Newton :
❙ erreur d'interpolation :
10
Newton : exemple
❚ (ex. n°2) : n=2 (0,1), (2,5) et (4,17)
a0
0 f [x0]=1
Newton : l’algorithme
❚ Interpolation en 1 point x :
Pour i allant de 1 à n
Pour j allant de 1 à i
Fi , j −1 − Fi −1, j −1
Fi , j =
xi − xi − j
FinPour
ai = Fi,i
FinPour
❚ Complexité du calcul : n2
Ph. Leray Analyse Numérique 21
11
A bas les polynômes
❚ Ex : y=2(1+tanh(x)) - x/10 avec 9 points
❙ entre les points, le polynôme fait ce qu'il veut…
et plus son degré est élevé plus il est apte à osciller !
6
-2
-4
-6
-6 -4 -2 0 2 4 6
12
Splines cubiques : définition
❚ Définition :
❙ On appelle spline cubique d’interpolation une fonction
notée g, qui vérifie les propriétés suivantes :
❘ g ∈ C2[a;b] (g est deux fois continûment dérivable),
❘ g coïncide sur chaque intervalle [xi; xi+1] avec un polynôme de
degré inférieur ou égal à 3,
❘ g(xi) = yi pour i = 0 … n
❚ Remarque :
❙ Il faut des conditions supplémentaires pour définir la
spline d’interpolation de façon unique
❙ Ex. de conditions supplémentaires :
❘ g"(a) = g"(b) = 0 spline naturelle.
Ph. Leray Analyse Numérique 24
❙ Notations :
❘ hi = xi+1 - xi pour i = 0 … n-1
❘ δi= [xi; xi+1]
❘ gi(x) le polynôme de degré 3 qui coïncide avec g sur l’intervalle δi
13
Splines cubiques : détermination 2/4
mi hi 2
❙ gi(xi) = yi yi = + bi 1
6 mi +1hi 2
❙ gi(xi+1) = yi+1 y i +1 = + ai hi + bi 2
6
❙ 1 et 2 ai =
1
(yi +1 − yi ) − hi (mi +1 − mi )
hi 6
14
Splines cubiques : détermination 4/4
1
hi −1 mi −1 + 2(hi + hi −1 )mi + hi mi +1 = 6 (yi +1 − yi ) − (yi − yi −1 )
1
4
hi hi −1
❘ T inversible (tridiagonale)
❚ Complexité du calcul : n
15
Splines cubiques : exemple
❚ Ex : y=2(1+tanh(x)) - x/10 avec 9 points
6
4
spline
2
-2
-4
polynôme de degré 8
-6
-6 -4 -2 0 2 4 6
Conclusion
❚ Interpolation polynomiale
❙ évaluer la fonction en un point :
polynôme de Lagrange ➨ méthode de Neville
❙ compiler la fonction (calcul des coefficients du polynôme) :
polynôme de Newton
❚ Interpolation polynomiale par morceau : splines
❙ spline cubique d’interpolation
❙ spline cubique d’approximation (on régularise)
❙ b spline
❙ spline généralisée : splines gaussiennes
16