Modèles de Taylor

En calcul scientifique, les modèles de Taylor sont une approche permettant d'approximer les fonctions mathématiques de manière rigoureuse à l'aide de polynômes de Taylor et d'intervalles d'encadrement.

Contexte

Au croisement des mathématiques et de l'informatique, le calcul certifié consiste à estimer de manière rigoureuse l'influence de différentes incertitudes numériques sur le calcul. Ces incertitudes proviennent surtout de deux sources :

• Des erreurs de calculs dues à la précision finie des ordinateurs.
• Des incertitudes sur les variables du modèle de calcul.

Les méthodes d'intervalles permettent de prendre en compte ces deux sources d'erreurs, mais elles entraînent un coût non négligeable en complexité, notamment du fait qu'il devient nécessaire de contrôler à tout instant ces intervalles de contrôle dont la taille peut exploser lors d'applications algorithmiques.

Définition

Si ${\displaystyle v}$ est un entier, soit ${\displaystyle f}$ une fonction supposée ${\displaystyle {\mathcal {C}}^{n 1}}$ sur ${\displaystyle D_{f}\subset R^{v}}$, et soit ${\displaystyle {\vec {B}}\subset [a_{1},b_{1}]\times \cdots \times [a_{v},b_{v}]\subset D_{f}}$ une boîte d'intervalles contenant le point ${\displaystyle {\vec {x_{0}}}}$.

Soit ${\displaystyle T}$ le polynôme de Taylor de ${\displaystyle f}$ en ${\displaystyle {\vec {x_{0}}}}$. L'intervalle ${\displaystyle I}$ est une borne du reste d'ordre ${\displaystyle n}$ pour ${\displaystyle f}$ sur ${\displaystyle {\vec {B}}}$ si

$${\displaystyle f({\vec {x}})-T({\vec {x}})\in I\quad {\text{pour tout }}{\vec {x}}\in {\vec {B}}}$$

La paire ${\displaystyle (T,I)}$ est un modèle de Taylor de ${\displaystyle f}$ d'ordre ${\displaystyle n}$. L'ensemble de toutes les bornes de restes est appelée la famille des restes.

Propriétés

Le but des modèles de Taylor est de servir d'objets de calculs. Il est donc nécessaire de définir les opérations élémentaires d'addition et de multiplication de modèles de Taylor, ainsi que des fonctions élémentaires comme l'exponentielle ou le logarithme. Par exemple, l'addition :

Si ${\displaystyle (T_{f},I_{f})}$ et ${\displaystyle (T_{g},I_{g})}$ sont des modèles de Taylor d'ordre ${\displaystyle n}$, le modèle de Taylor ${\displaystyle (T_{f} T_{g},I_{f} I_{g})}$ est un modèle de Taylor de ${\displaystyle f g}$ étant donné que

$${\displaystyle (f g)({\vec {x}})\in (T_{f} T_{g})({\vec {x}}) (I_{f} I_{g})\quad {\text{pour tout }}{\vec {x}}\in {\vec {B}}}$$

En appliquant ces règles de manière répétée, il est possible de calculer des modèles de Taylor de toutes les fonctions qui consistent en la répétition d'opérations d'addition, multiplication, et de fonctions élémentaires - ce qui inclut toutes les fonctions représentables par un ordinateur.

Applications

• Lors de la résolution d'équations différentielles ordinaires.
• Dans des problèmes d'optimisation sous contraintes.
• Dans la résolution d'équations différentielles algébriques (en).

Voir aussi

• Modèles de Tchebychev, une approche similaire avec des polynômes de Tchebychev
• Méthode de Newton
• Méthode formelle

Notes et références

Bibliographie

• M. Berz, From Taylor series to Taylor models, In AIP Conference Proceedings CONF-961208, vol. 405, n^o 1, avril 1997, p. 1-23)
• M. Berz & G. Hoffstätter, Computation and application of Taylor polynomials with interval remainder bounds, Reliable Computing, 4 (1), 1998 : 83-97.
• K. Makino & M. Berz, Efficient control of the dependency problem based on Taylor model methods, Reliable Computing, 5 (1), 1999 : 3-12.
• J. Hoefkens, Rigorous numerical analysis with high-order Taylor models, dissertation de doctorat, Michigan State University, Department of Mathematics and Department of Physics and Astronomy, 1997.
• K. Makino & M. Berz, Taylor models and other validated functional inclusion methods, International Journal of Pure and Applied Mathematics, 4 (4), 2003 : 379-456.
• N. Revol, K. Makino & M. Berz, Taylor models and floating-point arithmetic: proof that arithmetic operations are validated in Cosy, The Journal of Logic and Algebraic Programming, 64 (1), 2005 : 135-154.
• R. Zumkeller, Formal global optimisation with Taylor models, In International Joint Conference on Automated Reasoning, août 2006, p. 408-422.
• M. Neher, K. R. Jackson & N. S. Nedialkov, On Taylor model based integration of ODEs, SIAM Journal on Numerical Analysis, 45 (1), 2007 : 236-262.
• X. Zhang, S. Redon, M. Lee, & Y. J. Kim, Continuous collision detection for articulated models using taylor models and temporal culling, In ACM Transactions on Graphics (TOG), vol. 26, n^o 3, août 2007, p. 15).
• A. M. Sahlodin, & B. Chachuat, Convex/concave relaxations of parametric ODEs using Taylor models, Computers & Chemical Engineering, 35 (5), 2011 : 844-857.
• X. Chen, Reachability analysis of non-linear hybrid systems using taylor models, dissertation de doctorat, PhD thesis, RWTH Aachen University, 2005.

• Portail de l'analyse