Discrétisation — Wikipédia
![](http://upload.wikimedia.org/wikipedia/commons/thumb/a/a9/Finite_element_solution.svg/220px-Finite_element_solution.svg.png)
En mathématiques appliquées, la discrétisation est la transposition d'un état continu (fonction, modèle, variable, équation) en un équivalent discret. Ce procédé constitue en général une étape préliminaire à la résolution numérique d'un problème ou sa programmation sur machine. Un cas particulier est la dichotomisation où le nombre de classes discrètes est 2, où on peut approcher une variable continue en une variable binaire.
La discrétisation est aussi reliée aux mathématiques discrètes, et compte parmi les composantes importantes de la programmation granulaire. Dans le contexte, la discrétisation peut renvoyer à la modification de la granularité, quand plusieurs variables discrètes sont réunies ou des catégories discrètes fusionnées.
Discrétiser des données continues engendre systématiquement une erreur de discrétisation (en). Un des objectifs est donc de concevoir un modèle discret qui minimise au mieux cette erreur.
Il ne faut pas confondre discrétisation et quantification.
On compte également la méthode d'Euler-Maruyama (en) et le bloqueur d'ordre 0 parmi les méthodes de discrétisation.
Discrétisation de modèles d'état linéaires
[modifier | modifier le code]La discrétisation apparait dans la transformation d'équations différentielles continues en équations aux différence discrètes.
On considère le modèle d'état en espace, continu en temps :
où v et w sont des sources de bruit blanc avec une densité spectrale de puissance
peuvent être discrétisées, en supposant que le signal u est un bloqueur d'ordre 0 et une intégration continue pour le bruit v, donnant
avec des covariances
où
- , si A est régulière
et T est le temps d'échantillonnage, et est la transposée de A.
Une astuce pour calculer Ad et Bd en une étape consiste à utiliser la propriété[1]:p. 215
et donc
Discrétisation de bruits
[modifier | modifier le code]L'évaluation numérique de Qd est rendue plus délicate avec l'intégrale d'une exponentielle de matrice. On peut la calculer en deux temps, d'abord la construction de la matrice, puis le calcul de son exponentielle[2]
Le bruit discrétisé est ensuite évalué en multipliant la transposée du bloc en bas à droite de G avec celui en haut à droite :
Dérivation
[modifier | modifier le code]En partant du modèle continu
on sait que l'exponentielle de matrice est
et en multipliant à gauche le modèle :
on reconnait
L'intégration donne ainsi
ce qui est une solution analytique du modèle continu.
On veut désormais discrétiser cette expression. On suppose u constante sur chaque pas de temps.
On reconnait l'expression entre crochets dans le premier terme comme x[k], et le second terme peut être simplifié en faisant la substitution v(τ) = kT + T – τ, ce qui permet d τ = – d v. On suppose aussi que u est constante dans l'intégrale, ce qui donne :
qui est une solution exacte du problème de discrétisation.
Approximations
[modifier | modifier le code]Une discrétisation exacte peut parfois être impossible à cause d'une exponentielle de matrice lourde et des étapes d'intégrations. Il devient alors plus simple de calculer un modèle discret approché, basé sur de petits pas de temps de sorte qu'on ait . La solution approchée devient alors :
D'autres approximations possibles sont et . Chacun a des propriétés de stabilité différentes. On peut également mentionner la transformation bilinéaire, ou transformation de Tustin, qui préserve les propriétés de stabilité du système continu en temps.
Discrétisation d'équations différentielles
[modifier | modifier le code]![](http://upload.wikimedia.org/wikipedia/commons/thumb/d/d1/Finite_element_triangulation.svg/220px-Finite_element_triangulation.svg.png)
La résolution numérique d'une équation différentielle (ordinaire ou aux dérivées partielles) nécessite une discrétisation du domaine de définition de la solution (espace ou temps, voire les deux). Ainsi, d'une fonction u(x , t) définie sur un domaine Ω et un intervalle de temps [0 ; T], on ne calculera que des valeurs (u(xi , tn)), où les xi sont des points de Ω et tn des instants de [0 ; T]. Pour cela, les opérateurs différentiels sont également approchés par des versions discrètes, comme la dérivée seconde discrète :
La méthode de résolution (différences finies, éléments finis ou volumes finis, pour citer les plus courantes) permet de construire un problème discret dont la solution est une approximation de la solution du problème continu. L'erreur commise a deux sources :
- l'erreur de projection : en passant d'un espace continu à un espace discret, l'espace dans lequel la solution existante est changée ;
- l'erreur d'interpolation : le choix du schéma numérique et la définition de la grille espace-temps choisie pour la résolution va influer sur la qualité de l'approximation.
Discrétisation de caractéristiques continues
[modifier | modifier le code]En statistique et apprentissage machine, la discrétisation renvoie à la conversion de variables ou caractéristiques continues en variables ou caractéristiques discrètes nominales. Ce procédé est utile pour créer des fonctions de densité de probabilités.
Références
[modifier | modifier le code]- (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Discretization » (voir la liste des auteurs).
- (en) Raymond A. DeCarlo, Linear systems : A state variable approach with numerical implementation, Prentice-Hall, Inc., 1989.
- (en) Charles Van Loan, « Computing integrals involving the matrix exponential », IEEE transactions on automatic control, vol. 23, no 3, , p. 395-404.
Voir aussi
[modifier | modifier le code]Articles connexes
[modifier | modifier le code]- Espace discret
- Algèbre temporelle
- Simulation à événements discrets
- Simulation stochastique
- Méthode des volumes finis pour des flux instables
- Temps discret et temps continu
Liens externes
[modifier | modifier le code]- Robert Grover Brown & Patrick Y. C. Hwang, Introduction to random signals and applied Kalman filtering : with MATLAB exercises and solutions, 3rd, , 484 p. (ISBN 978-0-471-12839-7)
- Chi-Tsong Chen, Linear System Theory and Design, Philadelphia, PA, USA, Saunders College Publishing, (ISBN 0-03-071691-8)
- (en) Charles Van Loan, « Computing integrals involving the matrix exponential », IEEE transactions on automatic control, vol. 23, no 3, , p. 395-404
- R.H. Middleton et G.C. Goodwin, Digital control and estimation : a unified approach, , 33 p. (ISBN 0-13-211665-0)