Bibliographie & annexes
Quiz & test
Présentation
En anglais
RÉSUMÉ
On rappelle les définitions d’une fonction f(A) d’une matrice carrée à coefficients réels ou complexes. On présente ensuite des méthodes numériques permettant de calculer les éléments de f(A), l’action sur un vecteur donné f(A)v ou des formes bilinéaires utf(A)v où interviennent deux vecteurs.
Lire cet article issu d'une ressource documentaire complète, actualisée et validée par des comités scientifiques.
Lire l’articleABSTRACT
We first recall what is a matrix function f(A) of a square matrix with real or complex entries.Then we describe numerical methods to compute all the entries of f(A), the action on a given vector f(A)v or bilinear forms utf(A)v with two given vectors.
Auteur(s)
-
Gérard MEURANT : Ancien directeur de recherche au CEA
INTRODUCTION
Cet article est consacré au calcul de fonctions de matrices. Avant de définir ce qu’elles sont, expliquons brièvement ce qu’elles ne sont pas. Supposons que l’on ait une fonction f suffisamment régulière et une matrice carrée A d’ordre n à coefficients a i,j, réels ou complexes. La matrice f (A) d’ordre n n’est pas la matrice dont les éléments sont f (a i,j), auquel cas le calcul serait trivial. Les définitions de f (A) rappelées ci-dessous visent à reproduire, pour une matrice, la plupart des propriétés des fonctions scalaires. Dans une première partie, on présentera les définitions et les principales méthodes de calcul de tous les éléments de f (A). Cette partie est inspirée du livre qui contient l’état de l’art concernant le calcul de f (A), encore que certaines des méthodes décrites aient été légèrement améliorées depuis la parution de ce livre. On pourra également consulter avec profit.
Les algorithmes pour f (A) visent à calculer les n 2 éléments de la matrice. On utilise souvent des méthodes basées sur des factorisations de la matrice A à l’aide de transformations orthogonales et/ou des approximations de la fonction f permettant un calcul plus facile, par exemple des polynômes ou des fractions rationnelles. Les algorithmes correspondants ont donc un coût proportionnel à n 3 . Il n’est donc pas faisable, même avec les ordinateurs puissants dont on dispose aujourd’hui, de calculer f (A) pour des matrices de très grande taille. Il se trouve que de nombreuses applications n’ont besoin que de calculer f (A)v où v est un vecteur donné. Ceci peut être fait, sans calculer explicitement tous les éléments de f (A), à l’aide de méthodes itératives de Krylov qui peuvent s’appliquer à de très grandes matrices creuses et que nous décrirons dans une deuxième partie.
Enfin, il existe d’autres applications pour lesquelles on n’a besoin que de calculer des scalaires uT f (A)v, u et v étant des vecteurs donnés. Les méthodes pour calculer efficacement des bornes ou des approximations de ces quantités seront présentées dans une troisième et dernière partie.
L'expertise technique et scientifique de référence
KEYWORDS
matrix functions | mathematical software | scientific computing | engineering mathematics
DOI (Digital Object Identifier)
Cet article fait partie de l’offre
Mathématiques
(166 articles en ce moment)
Cette offre vous donne accès à :
Une base complète d’articles
Actualisée et enrichie d’articles validés par nos comités scientifiques
Des services
Un ensemble d'outils exclusifs en complément des ressources
Un Parcours Pratique
Opérationnel et didactique, pour garantir l'acquisition des compétences transverses
Doc & Quiz
Des articles interactifs avec des quiz, pour une lecture constructive
Présentation
Conclusion
En anglais
4. Méthodes pour uTf (A)v
Dans cette section on décrit comment utiliser des méthodes de quadrature pour calculer des approximations de uTf (A)v (voir ). On suppose que la matrice A est symétrique réelle et on utilise la décomposition spectrale de A qui s'écrit A = ZDZT , où Z est une matrice orthonormale dont les colonnes sont les vecteurs propres normalisés de A et D est une matrice diagonale dont les éléments diagonaux sont les valeur propres λ i de A, que l'on ordonne
Comme on l'a vu, la définition de f (A) est :
La forme bilinéaire uTf (A)v peut s'écrire :
La dernière somme peut être considérée comme une intégrale de Riemann-Stieltjes :
où la mesure α est constante par morceaux et définie par :
Lorsque u = v on a γ = β et α est une fonction positive croissante ainsi que lorsque γjβj > 0....
L'expertise technique et scientifique de référence
TEST DE VALIDATION ET CERTIFICATION CerT.I. :
Cet article vous permet de préparer une certification CerT.I.
Le test de validation des connaissances pour obtenir cette certification de Techniques de l’Ingénieur est disponible dans le module CerT.I.
de Techniques de l’Ingénieur ! Acheter le module
Cet article fait partie de l’offre
Mathématiques
(166 articles en ce moment)
Cette offre vous donne accès à :
Une base complète d’articles
Actualisée et enrichie d’articles validés par nos comités scientifiques
Des services
Un ensemble d'outils exclusifs en complément des ressources
Un Parcours Pratique
Opérationnel et didactique, pour garantir l'acquisition des compétences transverses
Doc & Quiz
Des articles interactifs avec des quiz, pour une lecture constructive
Méthodes pour uTf (A)v
BIBLIOGRAPHIE
-
(1) - AFANASJEW (M.), EIERMANN (M.), ERNST (O.G.), GÜTTEL (S.) - Implementation of a restarted Krylov subspace method for the evaluation of matrix functions - Linear Algebra Appl., vol. 429, p. 2293-2314 (2008).
-
(2) - AL-MOHY (A.H.), HIGHAM (N.J.) - A new scaling and squaring algorithm for the matrix exponential - SIAM J. Matrix Anal. Appl., vol. 31 n° 3, p. 970-989 (2009).
-
(3) - AL-MOHY (A.H.), HIGHAM (N.J.) - Computing the action of the matrix exponential with an application to exponential integrators - SIAM J. Sci. Comput., vol. 33 n° 2, p. 488-511 (2011).
-
(4) - AL-MOHY (A.H.), HIGHAM (N.J.) - Improved inverse scaling and squaring algorithms for the matrix logarithm - SIAM J. Sci. Comput., vol. 34 n° 4, p. C153-C169 (2012).
-
(5) - BREZINSKI (C.), VAN ISEGHEM (J.) - Padé approximations - in « Handbook of Numerical Analysis », vol. III, P.G. Ciarlet and J.L. Lions eds., North-Holland, Amsterdam, p. 47-222 (1994).
- ...
DANS NOS BASES DOCUMENTAIRES
ANNEXES
CISIA juin 2000 Le Bayésien (version pour Windows Vista)
CISIA 1 avenue Herbillon, 94160 Saint-Mandé, France
Les outils logiciels disponibles sont décrits dans le texte
HAUT DE PAGE
Boîte à outils mft_toolbox pour Matlab :
http://www.maths.manchester.ac.uk/∼higham
Boîte à outils expohit pour Matlab :
http://www.maths.uq.edu.au/expokit/
Fonctions Matlab contenant des intégrateurs exponentiels pour les équations différentielles :
http://www.math.ntnu.no/num/expint/matlab.php
Fonctions Matlab utilisant la méthode d'Arnoldi avec redémarrage pour une fonction quelconque :
Boîte à outils mmq_toolbox pour Matlab :
http://gerard.meurant.pagesperso-orange.fr/
HAUT DE PAGE
L'expertise technique et scientifique de référence
Cet article fait partie de l’offre
Mathématiques
(166 articles en ce moment)
Cette offre vous donne accès à :
Une base complète d’articles
Actualisée et enrichie d’articles validés par nos comités scientifiques
Des services
Un ensemble d'outils exclusifs en complément des ressources
Un Parcours Pratique
Opérationnel et didactique, pour garantir l'acquisition des compétences transverses
Doc & Quiz
Des articles interactifs avec des quiz, pour une lecture constructive
QUIZ ET TEST DE VALIDATION PRÉSENTS DANS CET ARTICLE
1/ Quiz d'entraînement
Entraînez vous autant que vous le voulez avec les quiz d'entraînement.
2/ Test de validation
Lorsque vous êtes prêt, vous passez le test de validation. Vous avez deux passages possibles dans un laps de temps de 30 jours.
Entre les deux essais, vous pouvez consulter l’article et réutiliser les quiz d'entraînement pour progresser. L’attestation vous est délivrée pour un score minimum de 70 %.
L'expertise technique et scientifique de référence
Cet article fait partie de l’offre
Mathématiques
(166 articles en ce moment)
Cette offre vous donne accès à :
Une base complète d’articles
Actualisée et enrichie d’articles validés par nos comités scientifiques
Des services
Un ensemble d'outils exclusifs en complément des ressources
Un Parcours Pratique
Opérationnel et didactique, pour garantir l'acquisition des compétences transverses
Doc & Quiz
Des articles interactifs avec des quiz, pour une lecture constructive
