Roberto Andrés Martínez Merchán Méthodes numériques

2161881

EXTRAPOLATION DE RICHARDSON

Introduction.
Dans de nombreux problèmes d'application, il est nécessaire d'obtenir des résultats.
numériques, qui en plus d'être de bons estimations de valeurs réelles, n'impliquent pas un
haut effort computationnel.
Le processus d'obtention d'une estimation améliorée pour la valeur des dérivées, des équations
différentiels, entre autres, sur la base de l'application d'une ou plusieurs formules, en utilisant
Différentes longueurs d'intervalle se désignent par extrapolation.

La différenciation numérique.
C'est une technique d'analyse numérique qui nous permet de calculer une approximation à la
dérivée d'une fonction en un point en utilisant les valeurs et propriétés de celle-ci. La
dérivée d'une fonction f en xo, il est défini de la manière suivante :

Pour de petites valeurs de h, nous pouvons approcher la dérivée de f en xode la suivante

manière

Nous savons que cette équation entraîne une erreur associée qui est déterminée
en nous servant de la série de Taylor, cela amène avec lui une détermination où naissent les
proposiciones ‘Diferenciación hacia delante ’,‘Diferenciación hacia atrás ’ y
‘Différenciation centrée’. Où chacun se comprend plus facilement ci-dessous
comme :
 a) Différenciation en amont.

Où ∆fi est connu comme la première différence

vers l'avant et à h on l'appelle la taille du pas ou
increment; c'est-à-dire, la longueur de l'intervalle sur
lequel se réalise l'approche. O(h) est l'erreur de
truncation.

b) Différenciation en amont.

Tronquant l'équation après la première

en dérivant et en réarrangeant les termes, on obtient :

Où l'erreur est O(h).

c) Différenciation centrée.
Ici, les développements de Taylor sont soustraits.
auparavant, obtenant ainsi :

Notez que l'erreur de troncature est de l'ordre

de h2en contraste avec les approches vers
en avant et en arrière. Par conséquent, ceci est un
représentation plus exacte de la première dérivée.
Extrapolation de Richardson.
Cette technique améliore la précision du calcul numérique de la dérivée d'une fonction,
partant de la base de la série de Taylor. Cela survient en raison du désir de réduire l'erreur
existent dans les méthodes utilisées auparavant.
Cette méthode est basée sur l'utilisation de deux estimations de la dérivée pour calculer un
troisième approximation, qui sera beaucoup plus précise.
Pour atteindre cette précision, nous utiliserons la formule :

1
≅ ( ℎ2)+ ℎ [ ( ℎ2)− ℎ( 1] ) (1.01)
( ℎ 1 ) 2−1
2

Où D (ℎ1 ) y D (ℎ2 ) ses estimations de la dérivée obtenues en utilisant deux tailles

de paso, h1y h2.
⁄ de maintenant en avant est ce qui sera recherché dans la
ℎ1
Quand ils conviennent queℎ2= (que
2
la majorité des exercices à résoudre) l'équation (1.01) se transformera en :
Partant de la dérivée centrale, qui a déjà été vue précédemment, elle génère
meilleures approches et ayant les deux étapesℎ, ℎ/2.

Maintenant, pour le pas égal à h.

(1)
Pour le pas égal à h/2.

(2)

Multipliquons (1) par ¼ et soustrayons (2) alors,

1 ′(
1 ( ) − (
+ℎ − )ℎ 1
) = + (ℎ 2 )
4 4 2ℎ 4

′(
( )
+ ℎ/2− ( − ℎ/2) 1
− ) = + (ℎ 2 )
ℎ 4
ℎ ℎ
3 1 ( +ℎ )− ( − ℎ) ( +2 )− ( − 2)
− ′( ) = −4 + (3)
4 2ℎ ℎ
 ( + ) −( −) ( + )− ( −)
Maintenant, soit ( ) = y ( )= ,remplacé
en (3) nous avons :

3 ′(
1 ℎ
− )= (ℎ ) − ( )
4 4 2
Où :

′( ) = ( ) − ( ) [1.02]
NOTES : Lorsque les étapes prises correspondent à ce qu'elles sont
exactement la moitié de l'ancien sera utilisé l'équation [1.02] et
quand non, la générale [1.01]

Ceci étant l'une des formules que nous utiliserons pour l'estimation de l'extrapolation
de Richardson, où les étapes prises coïncident être :ℎ2= ℎ1⁄
2
Pour la détermination de l'erreur, il faut connaître la valeur réelle exacte de la
la dérivée au point évalué, cela dépendra des pas h qui peuvent être
prendre.
− ℎ
| = | ∗ 100

En résumé, lorsqu'on a la relation mentionnée entre les étapes h1 et h2,

on a ce qui suit.

′(
4 1
) = ( ℎ) 2− (ℎ1)
3 3
( +ℎ) − ( )
−ℎ
Où : ℎ ( =) ;
2ℎ

Exercice.
Calculer la première dérivée de ( )= 0.5x+ 0.25x 4 3
+ 0,75 − 2+
3.2, au point xo =0.5 ayant des pas h1=0.5 et h2=0.25. Déterminer
l'erreur qui se produit à chaque estimation de la dérivée sachant
que la valeur exacte de la première dérivée évaluée en ce point est
0,1875.
Solution.

On a que les étapes sont liées par h2=h1/2, c'est pourquoi nous utilisons le second
situation présentée.

Alors, h2=0.5/2=0.25

Remplaçant.

Pour h1=0.5
( 0,5+ 0,5 −
) ( − 0,5 )
0.5 ( 1−
) 0( )
( )
0,5= =
( )
20,5 1
( 0.5 ( 0.5+) 40.25 0,5( + 0,750,5−
)3 0,5 ) 2 0,5 0+ 0.25
( +3.2− ) 0(+ 0,750−
( )4 0 +3.2( 3) ( )2 )
=
1
3.7- 3.2
( )
0.5= = 0,5
1
( 0.1875− 0,5 )
= ( ) ∗ 100166,6 %
0.1875

Pour h2=0.25
( 0.5+ 0.25−) 0.5 ( − 0.25 ) ( 0,75−
) (0.25))
(
0.25 )= =
2(0.25) 0,5
( 0,5+( 0.25) 40,75 + 0,750,75−
( ) 3 0,75 +3.2−
( 2
)0,5 0.25 +0.25) 0.25 ) 4 0.25− 0.25
( (+ 0,75 )3
( +3.2 ( )2 )
=
0,5
3.1355- 3.0027
(
D'po 0.)25 = = 0.2656
0,5
( 0.1875− 0.2656 )
= ∗ 100 = 41.65%
0.1875

Maintenant en appliquant la formule suivante :

′(
4 1 ′(
4 1
) =
= 0.5 ( )ℎ2− ( ℎ 1=
) 0,5 =) 0.2656(− 0,5 ) ( )
3 3 3 3
′( = 0.5) = 0.18746
| (0.1875− 0.18746| )
= ∗ 100 = 0,05%
0.1875

Nous voyons que l'erreur qui se produit est minimale, elle peut être attribuée à des chiffres qui se perdent.
dans les calculs. Pour les fonctions qui ne sont pas des polynômes, cette estimation également
fonctionnerait, mais avec un peu plus de probabilité que l'erreur soit plus grande.
Bibliographie.
• https://sites.google.com/site/procesosnumericoslrl20151/segunda-
livraison/méthodes-itératives
• Méthodes numériques pour ingénieurs. 5e éd. - Steven C. Chapra & Raymond P.
canal.
• Unable to access content from the provided link.
extrapolation-de-richardson-pour-derivation-numerique.html
Étapes pour la programmation en Matlab de la méthode :

1. S'assurer d'avoir toutes les fenêtres propres (Command Windows,

Espace de travail…).
2. Rendre symbolique la lettre x dans Matlab qui sera notre variable.
3. Demander à l'utilisateur la fonction à travailler.
4. Demander à l'utilisateur le point (xo) où il souhaite connaître son premier
dérivée dans la fonction.
5. Demander les étapes qui seront prises, il faut clarifier que celles-ci
Ils entrent par ordre décroissant.
6. Déterminer les approximations D(h1) et D(h2).
7. Vérifier si h2=h1/2.
8. Dans le cas où cela coïncide :
a. Évaluer les approximations calculées précédemment dans la fonction
[1.02].
b. Connaissant la valeur réelle de la fonction évaluée au point et celle de la
évaluation réalisée au point précédent, calculer l'erreur
pourcentage que cela peut générer.
9. À ces endroits ne pas coïncider :

a. Saisir les valeurs de h1 et h2 ainsi que celles des approximations dans le

fonction générale [1.01].
b. Connaissant la valeur réelle de la fonction évaluée au point et celle de la
évaluation réalisée au point précédent, calculer l'erreur pourcentage
que cela puisse générer.
10. Imprimer les résultats trouvés ; la valeur de xo dans la première dérivée
de la fonction et de l'erreur.