Équation de Fokker-Planck


Équation de Fokker-Planck

L'équation de Fokker-Planck est une équation aux dérivées partielles linéaire que doit satisfaire la densité de probabilité de transition d'un processus de Markov. A l'origine, une forme simplifiée de cette équation a permis d'étudier le mouvement brownien. Comme la plupart des équations aux dérivées partielles, elle ne donne des solutions explicites que dans des cas bien particuliers portant à la fois sur la forme de l'équation, sur la forme du domaine où elle est étudiée (conditions réfléchissante ou absorbante pour les particules browniennes et forme de l'espace dans lequel elles sont confinées par exemple).

Sommaire

Équation

Cette équation concerne un processus de Markov, processus stochastique qui possède la propriété markovienne : la probabilité d'apparition d'un état du système à un instant donné ne dépend que de son histoire la plus récente. Ainsi le processus est caractérisé par une probabilité de transition entre deux états assortie d'une condition initiale.

Pour éviter les incohérences, la probabilité de transition doit satisfaire une condition de compatibilité qui s'exprime par une équation intégrale appelée équation de Chapman - Kolmogorov - Smoluchowski. Celle-ci signifie qu'étant donnés trois états successifs, la transition de l'état 1 à l'état 3 est obtenue en effectuant la transition de 1 à 2, puis la transition de 2 à 3. Si le processus est stationnaire, la densité de probabilité du processus est obtenue en faisant tendre la durée de la transition vers l'infini.

L'équation intégrale se transforme en une équation aux dérivées partielles selon une méthode définie par Wang et Uhlenbeck. L'équation présentée ci-dessous, limitée au deuxième ordre, n'est pas la plus générale.

À une dimension, l'équation de Fokker Planck avec un coefficient de diffusion D2(x,t)et la tendance (un drift) D1(x,t) s'écrit :

\frac{\partial}{\partial t}P(x,t)=-\frac{\partial}{\partial x}\left[ D_{1}(x,t)P(x,t)\right] + {1 \over 2} \frac{\partial^2}{\partial x^2}\left[ D_{2}(x,t)P(x,t)\right].

P est la probabilité de trouver la particule au point x et à l'instant t.

Dans le cas plus général à N variables la densité de probabilité jointe vérifie l'équation :

\frac{\partial f}{\partial t} = -\sum_{i=1}^{N} \frac{\partial}{\partial x_i} \left[ D_i^1(x_1, \ldots, x_N) f \right] + {1 \over 2} \sum_{i=1}^{N} \sum_{j=1}^{N} \frac{\partial^2}{\partial x_i \, \partial x_j} \left[ D_{ij}^2(x_1, \ldots, x_N) f \right],

D1 est le vecteur de tendance (drift) et D2 le tenseur de diffusion

Cas du mouvement brownien

Dans le cas d'un mouvement d'une particule et dans le cadre de l'équation de Smoluchowski qui concerne les particules telles que  \gamma v\gg m  a , typiquement des molécules, ou objets de masse "négligeable" (molécules atmosphériques, protéines en biologie...) :

\dot{X}+\frac{F(X)}{\gamma}=\sigma B

où B est un bruit blanc, \textstyle \gamma le coefficient de viscosité et F(x) un champ de forces. Si p(x,t) est la probabilité de trouver la particule au point x à l'instant t, par application du lemme d'Itô on a alors :

 \frac{\partial p(x,t)}{\partial t}=\frac{\sigma^2}{2\gamma^2} \triangle p (x,t) +\nabla \left(\frac{F(x)}{\gamma} p(x,t)\right)

\textstyle \frac{\sigma^2}{2\gamma^2}=D le coefficient de diffusion.

Cette équation de Fokker-Planck particulière permet alors, avec des conditions aux bords et à l'origine adéquates, d'étudier le mouvement brownien d'une particule dans un champ de forces.

Cas des systèmes dynamiques non linéaires

Rappel sur les systèmes linéaires et linéarisables

Chercher la réponse d'un système donné à une excitation est un problème courant en physique et dans diverses techniques. Cette réponse est généralement définie par une équation différentielle ou plus généralement par un système différentiel à plusieurs variables. Pour la simplicité de l'exposé on s'en tiendra à une équation classique du second ordre, cette équation correspondant en particulier aux oscillations d'un système mécanique muni d'une force de rappel et d'un amortissement (voir Systèmes oscillants à un degré de liberté).

Le problème le plus courant concerne une excitation sinusoïdale. Si l'équation est linéaire la réponse est elle-même sinusoïdale. Au contraire, lorsqu'elle contient des termes non linéaires on ne peut en général trouver que des solutions approchées : solutions numériques, recherche d'une approximation sinusoïdale par la linéarisation équivalente ou développement en une suite de sinusoïdes par la méthode des perturbations.

On rencontre parfois aussi des équations différentielles stochastiques dans lesquelles l'excitation est représentée par un processus aléatoire qui est en gros un ensemble de fonctions possédant les mêmes propriétés statistiques. Dans ces conditions la solution est également un processus aléatoire caractérisé par une densité de probabilité qui concerne, pour le second ordre, les excursions et les vitesses. Si on se limite à une excitation gaussienne, la réponse donnée par une équation linéaire est également gaussienne et peut se déterminer au moyen des techniques de description spectrale.

Systèmes non-linéaires excités par un bruit blanc

Dans le cas général un système mécanique à un degré de liberté est décrit par une équation différentielle de la forme

m \ddot{X} + g(X,\dot{X}) = F(t)

Si l'excitation peut être considérée comme une réalisation d'un processus aléatoire F(t)\, on parle d'équation différentielle stochastique dont la solution qui représente l'excursion est elle-même un processus aléatoire X(t)\,. Dans ces conditions, il est commode d'introduire le processus vitesse U(t)\,, ce qui conduit à deux équations du premier ordre :

\dot{X} = U
m \dot{U} + g(X,U) = F(t)

En considérant une petite variation \Delta t\, du temps, les équations différentielles sont remplacés par des équations aux différences (voir Méthode des différences finies) :

X_t = X_{t-dt} + U_t \Delta t\,
U_t = U_{t-dt} + {1 \over m} (F_t -g(X_t,U_t)) \Delta t

Si l'excitation est un bruit blanc gaussien, succession d'impulsions indépendantes, l'état du système ne dépend donc que de ce qui se passait à l'instant précédent. Il possède la propriété markovienne et une équation de Fokker-Planck permet de déterminer la densité de probabilité jointe du mouvement et de la vitesse. L'approximation du bruit blanc est d'autant meilleure que le système est moins amorti (résonance plus aiguë).

Comme la plupart des équations aux dérivées partielles, elle ne donne des solutions explicites que dans des cas bien particuliers portant à la fois sur la forme de l'équation et sur la nature de l'excitation. La méthode s'applique en particulier à l'équation possédant une force de rappel non linéaire :

m \ddot{X} + b \dot{X} + h(X) = F(t)

dans laquelle F(t)\, représente un bruit blanc gaussien de densité spectrale S_0\, (densité exprimée en unités au carré par radian par seconde).

La densité de probabilité jointe de l'excursion et de la vitesse s'écrit

p_{X \dot{X}}(x,\dot{x}) = C e^{-{b \mathcal{E}}/{\pi S_0}}

Dans cette formule, \mathcal{E} = {1\over 2}\,m\,{\dot{x}^2}\,+\,\int_0^x h(\xi) d\xi est l'énergie totale du système non amorti. La densité de probabilité de la vitesse reste gaussienne tandis que celle de l'excursion ne l'est plus. Elle le redevient évidemment lorsque la fonction h\, est linéaire.

Il est remarquable que cette solution, relativement simple, soit exacte alors qu'il n'existe rien de tel pour une excitation sinusoïdale. D'autre part, une telle solution exacte n'existe plus si c'est l'amortissement qui est non linéaire ; dans ce cas il existe une solution exacte pour une équation plus abstraite qui fournit une approximation non linéaire meilleure que l'approximation linéaire.

Références

  1. Philippe A. Martin, Introduction aux processus stochastiques en physique, http://seminairetransalpin.in2p3.fr/2006/talks/martin.pdf
  2. Ch. Ancey, Simulations stochastiques - Applications aux écoulements géophysiques et la turbulence, http://www.toraval.fr/articlePDF/stochastique.pdf
  3. I. Kosztin, Non-Equilibrium Statistical Mechanics chap 4: Smoluchowski Diffusion Equation, http://www.ks.uiuc.edu/~kosztin/PHYCS498NSM/LectureNotes/chp4.pdf
  4. H. Risken, "The Fokker–Planck Equation: Methods of Solutions and Applications", 2nd edition, Springer Series in Synergetics, Springer, ISBN 3-540-61530-X.
  5. S. Redner: A guide to first passage processes. Cambridge University Press, (2001)
  6. (en) Y. K. Lin, Probabilistic Theory of Structural Dynamics, New York, Robert E. Krieger Publishing Company, juillet 1976, 368 p. (ISBN 978-0-88275-377-5) (LCCN 75042154) 
  1. Thomas L. Paez, Random Vibrations : Assessment of the Sate of the Art, http://www.osti.gov/bridge/servlets/purl/3882-OeSEmc/webviewable/3882.pdf
  2. Thomas K. Caughey, Derivation and Application of the Fokker-Planck Equation to Discrete Nonlinear Dynamic Systems Subjected to White Random Excitation, http://authors.library.caltech.edu/4087/01/CAUjasa63a.pdf

Voir aussi

  1. formule d'Itô
  2. mouvement brownien

Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Équation de Fokker-Planck de Wikipédia en français (auteurs)

Regardez d'autres dictionnaires:

  • Equation de Fokker-Planck — Équation de Fokker Planck L équation de Fokker Planck est une équation aux dérivées partielles linéaire que doit satisfaire la densité de probabilité de transition d un processus de Markov. A l origine, une forme simplifiée de cette équation a… …   Wikipédia en Français

  • Équation de fokker-planck — L équation de Fokker Planck est une équation aux dérivées partielles linéaire que doit satisfaire la densité de probabilité de transition d un processus de Markov. A l origine, une forme simplifiée de cette équation a permis d étudier le… …   Wikipédia en Français

  • Fokker–Planck equation — [ thumb|A solution to the one dimensional Fokker–Planck equation, with both the drift and the diffusion term. The initial condition is a Dirac delta function in x = 1, and the distribution drifts towards x = 0.] The Fokker–Planck equation… …   Wikipedia

  • Fokker-Planck-Gleichung — Lösung der 1D Fokker Planck Gleichung mit Drift und Diffusionsterm. Die Anfangsbedingung ist eine Deltafunktion bei x = 1 und die Verteilung driftet nach links. Die Fokker Planck Gleichung (FPG, nach Adriaan Daniël Fokker (1887–1972) und Max… …   Deutsch Wikipedia

  • Ecuación de Fokker-Planck — Saltar a navegación, búsqueda La ecuación de Fokker–Planck, denominada así por Adriaan Fokker y Max Planck, y también conocida como ecuación avanzada de Kolmogórov (por Andréi Kolmogórov), describe la evolución temporal de la función de densidad… …   Wikipedia Español

  • Equation differentielle stochastique — Équation différentielle stochastique Une équation différentielle stochastique (EDS) est une généralisation de la notion d équation différentielle prenant en compte un terme de bruit blanc. Les EDS permettent de modéliser des trajectoires… …   Wikipédia en Français

  • Equation maitresse — Équation maîtresse En physique, une équation maîtresse est une équation différentielle décrivant l évolution temporelle d un système. C est une équation de taux pour les états du système. L évolution de la probabilité Pk d être dans l état… …   Wikipédia en Français

  • Équation différentielle stochastique — Une équation différentielle stochastique (EDS) est une généralisation de la notion d équation différentielle prenant en compte un terme de bruit blanc. Les EDS permettent de modéliser des trajectoires aléatoires, tels des cours de bourse ou les… …   Wikipédia en Français

  • Équation de Chapman-Kolmogorov — En théorie des probabilités, et plus spécifiquement dans la théorie des processus stochastiques markoviens, l équation de Chapman Kolmogorov est une égalité qui met en relation les lois jointes de différents points de la trajectoire d un… …   Wikipédia en Français

  • Équation maîtresse — En physique, une équation maîtresse est une équation différentielle décrivant l évolution temporelle d un système. C est une équation de taux pour les états du système. L évolution de la probabilité Pk d être dans l état discret k suit une… …   Wikipédia en Français


We are using cookies for the best presentation of our site. Continuing to use this site, you agree with this.