Introduction To Stochastic Process: Julia Hinkel
Introduction To Stochastic Process: Julia Hinkel
Introduction To Stochastic Process: Julia Hinkel
Outline
Basic concepts of probability theory; Markov process; Random walks; Master equation; FokkerPlank equation; From binomial distribution to Gaussian one;
Ficks law; Einsteins (Smoluchowskis) derivation of diffusion equation. Langevin formalizm; Wiener process; From microscopic description of stochastic differential equations to mesoscopic one; Fluctuationdissipation theorem.
Motivations
The study of complex systems coming from chemical, biological and social sciences; Investigations in this eld have been determined as interdisciplinary research. For example: Trafc ow; Granular matter; Ant colony behaviour; Transactions in nancial markets.
These systems are interesting not only as objects of natural sciences but also from the physical points of view for fundamental understanding and detailed analysis of such exotic phenomena. What does it mean to understand into details??? The main goal of such investigations is to invent a theoretical model which can describe the general features of the phenomenon under consideration; To provide the detailed theoretical analysis and computer simulation of these models.
What the difference between these two approaches? The way of the studying: To forecast the behaviour of the described object for the specic time (or OBSERVATION time) tobs .
Deterministic process = one its realization and everything clear; Stochastic process = one realization, two, ..., ten, ..., hundred, may be more??? Ensemble of realization; Probability distributions.
The solution for the initial condition m = n = 0 is the binomial distribution. Master Equation: Continuous time and discrete space process; Transformation t = n t with the xed time step t;
t 0 = Master equation p q p+q P (m, t) = P (m 1, t) + P (m + 1, t) P (m, t) t t t t
(3)
DriftDiffusion Equation
Continuous space and time process; Transformation x = m x with the xed space step x;
x 0; P (m, t) p(x, t) dx;
FokkerPlanck equation
p(x, t) p(x, t) 2 p(x, t) = v +D t x x2
(4)
Ficks Law
Diffusion Equation
2 p(x, t) p(x, t) =D t x2
(5)
Normalization condition:
+ p(x, t)
= 1;
continuity equation for certain position x: change of probability at a xed position in time equals to the spatial derivative of a probability ux:
p(x, t) jdif f (x, t) =0 t x
(6)
Brownian Motion
The Brownian motion which observed originally by Robert Brown in 1827 has initiated the development of a new branch in physics; By modelling of the nature phenomena, it has been concluded, that the wellinvestigated deterministic approach based on the Newton theory could not give the correct description of such processes; It has been found that the systems under study do not have the deterministic behaviour. The detailed experimental investigations of the Brownian motion have shown that the observed motion is irregular and unpredictable
First explanations of the Brownian motion have been given by Einstein and Smoluchowski and were provided in terms of probability density distribution; Langevin has used another approach based on Newtons law There are two forces having the inuence on the particle dynamics: deterministic force given by friction v ; stochastic force which uctuates rapidly (t).
The onedimensional stochastic differential equation in terms of the Langevin notation has the form
dx = (t) dt
(7)
The stochastic force (t) presented by Gaussian white noise; There is the problem that the Langevin equation is impossible to integrate; The integral of the random term is not differentiable; As result, another approach has to be used = stochastic integration; The stochastic integral is dened as a limit of the partial sums which include the increment of the Wiener process.
Wiener Process
The Wiener process is the simplest continuous in time onedimensional stochastic process and its dynamics obeys the equation
dx(t) = dW (t)
(8)
and variance
(W (t))2 W (t)
= t.
(10)
All increments
W (t) = W (t + t) W (t)
(11)
and
(W (t))2 = t .
(12)
For the same realization the Box-Mller method proposes the necessary transformation after the generation U1 , U2 (0, 1) only by one step
Z= 2 ln U1 cos (2 U2 ) .
(14)
The difference between both described methods is insignicant and connected only with the simulation time.
The uniform distribution at the interval [0, 1] 8 > 0 > < 1 > > : The mean value = 1/2, because of = Z
+
(15)
x dx =
1 2
(16)
1 0
1 x 2
dx =
1 . 12
(17)