Stochastic Analysis An Introduction
Prof. Dr. Andreas Eberle
June 23, 2015 Contents
Contents 2
1 Brownian Motion 8 1.1 FromRandomWalkstoBrownianMotion ...... 9 CentralLimitTheorem ...... 10 BrownianmotionasaLévyprocess...... 14 BrownianmotionasaMarkovprocess...... 15 WienerMeasure...... 18 1.2 BrownianMotionasaGaussianProcess ...... 20 Multivariatenormals ...... 21 Gaussianprocesses ...... 24 1.3 TheWiener-LévyConstruction ...... 31 Afirstattempt...... 32 TheWiener-Lévy representation of Brownian motion ...... 34 Lévy’sconstructionofBrownianmotion ...... 40 1.4 TheBrownianSamplePaths ...... 45 Typical Brownian sample paths are nowhere differentiable ...... 45 Höldercontinuity ...... 47 Lawoftheiteratedlogarithm ...... 48 Passagetimes ...... 49 1.5 Strong Markov property and reflection principle ...... 52 MaximumofBrownianmotion...... 52 StrongMarkovpropertyforBrownianmotion ...... 55
2 CONTENTS 3
Arigorousreflectionprinciple ...... 57
2 Martingales 59 2.1 DefinitionsandExamples...... 59 MartingalesandSupermartingales ...... 60 Somefundamentalexamples ...... 62 2.2 DoobDecompositionandMartingaleProblem ...... 69 DoobDecomposition ...... 69 ConditionalVarianceProcess...... 71 Martingaleproblem...... 74 2.3 Gamblingstrategiesandstoppingtimes ...... 76 Martingaletransforms...... 76 StoppedMartingales ...... 81 OptionalStoppingTheorems ...... 84 Wald’sidentityforrandomsums ...... 89 2.4 Maximalinequalities ...... 90 Doob’sinequality ...... 90 Lp inequalities ...... 92 Hoeffding’sinequality ...... 94
3 Martingales of Brownian Motion 97 3.1 Some fundamentalmartingalesof BrownianMotion ...... 98 Filtrations ...... 98 BrownianMartingales...... 99 3.2 OptionalSamplingTheoremandDirichletproblem ...... 101 OptionalSamplingandOptionalStopping ...... 101 Ruinprobabilitiesandpassagetimesrevisited ...... 104 ExitdistributionsandDirichletproblem ...... 105 3.3 MaximalinequalitiesandtheLIL...... 108 Maximalinequalitiesincontinuoustime ...... 108 ApplicationtoLIL ...... 110
University of Bonn Winter Term 2010/2011 4 CONTENTS
4 Martingale Convergence Theorems 114 4.1 Convergence in L2 ...... 114 Martingales in L2 ...... 114 Theconvergencetheorem...... 115 Summabilityofsequenceswithrandomsigns ...... 117 L2 convergenceincontinuoustime ...... 118 4.2 Almostsureconvergenceofsupermartingales ...... 119 Doob’supcrossinginequality ...... 120 ProofofDoob’sConvergenceTheorem ...... 122 Examplesandfirstapplications ...... 123 Martingales with bounded increments and a Generalized Borel-Cantelli Lemma126 Upcrossing inequality and convergence theorem in continuoustime . . 128 4.3 Uniform integrability and L1 convergence ...... 129 Uniformintegrability ...... 130 Definitive version of Lebesgue’s Dominated Convergence Theorem . . 133 L1 convergenceofmartingales ...... 135 BackwardMartingaleConvergence...... 138 4.4 Local and global densities of probability measures ...... 139 AbsoluteContinuity...... 139 Fromlocaltoglobaldensities...... 141 Derivationsofmonotonefunctions ...... 146 Absolutecontinuityofinfiniteproductmeasures ...... 146 4.5 TranslationsofWienermeasure ...... 150 TheCameron-MartinTheorem ...... 152 Passage times for Brownian motionwith constant drift ...... 157
5 Stochastic Integral w.r.t. Brownian Motion 159 5.1 Definingstochasticintegrals ...... 161 Basicnotation...... 161 Riemannsumapproximations ...... 161 Itôintegralsforcontinuousboundedintegrands ...... 163 5.2 SimpleintegrandsandItôisometryforBM ...... 165
StochasticAnalysis–AnIntroduction Prof.AndreasEberle CONTENTS 5
Predictablestepfunctions...... 165 Itôisometry;Variant1 ...... 167 Definingstochasticintegrals: Asecondattempt ...... 170 2 The Hilbert space Mc ...... 171 2 Itô isometry into Mc ...... 174 5.3 Itôintegrals for square-integrableintegrands ...... 174 DefinitionofItôintegral ...... 175 ApproximationbyRiemann-Itôsums ...... 176 Identificationofadmissibleintegrands ...... 177 Localdependenceontheintegrand ...... 180 5.4 Localization...... 182 Itô integrals for locally square-integrable integrands ...... 182 Stochasticintegralsaslocalmartingales ...... 184 ApproximationbyRiemann-Itôsums ...... 186
6 Itô’s formula and pathwise integrals 188 6.1 Stieltjesintegralsandchainrule ...... 189 Lebesgue-Stieltjesintegrals...... 190 ThechainruleinStieltjescalculus ...... 193 6.2 Quadr. variation, Itô formula, pathwise Itô integrals ...... 195 Quadraticvariation ...... 195 Itô’s formula and pathwise integrals in R1 ...... 198 Thechainruleforanticipativeintegrals ...... 201 6.3 First applications to Brownian motion and martingales ...... 202 QuadraticvariationofBrownianmotion ...... 203 Itô’sformulaforBrownianmotion ...... 205 6.4 Multivariateandtime-dependentItôformula ...... 210 Covariation ...... 210 ItôtoStratonovichconversion ...... 212 Itô’s formula in Rd ...... 213 Productrule,integrationbyparts ...... 215 Time-dependentItôformula ...... 216
University of Bonn Winter Term 2010/2011 6 CONTENTS
7 Brownian Motion and PDE 221 7.1 Recurrence and transience of Brownian motion in Rd ...... 222 TheDirichletproblemrevisited...... 222 Recurrence and transience of Brownian motion in Rd ...... 223 7.2 Boundaryvalue problems,exit and occupationtimes ...... 227 ThestationaryFeynman-Kac-Poissonformula ...... 227 Poissonproblemandmeanexittime ...... 229 OccupationtimedensityandGreenfunction ...... 230 Stationary Feynman-Kac formula and exit time distributions ...... 232 Boundary value problems in Rd andtotaloccupationtime...... 233 7.3 Heatequationandtime-dependentFKformula...... 236 BrownianMotionwithAbsorption ...... 236 Time-dependentFeynman-Kacformula ...... 239 Occupationtimesandarc-sinelaw ...... 241
8 SDE: Explicit Computations 244 8.1 StochasticCalculusforItôprocesses ...... 247 Stochasticintegralsw.r.t.Itôprocesses ...... 248 CalculusforItôprocesses...... 251 The Itô-Doeblin formula in R1 ...... 253 MartingaleproblemforsolutionsofSDE ...... 255 8.2 Stochasticgrowth ...... 256 Scalefunctionandexitdistributions ...... 257 Recurrenceandasymptotics ...... 258 GeometricBrownianmotion ...... 259 Feller’sbranchingdiffusion...... 261 Cox-Ingersoll-Rossmodel ...... 263 8.3 LinearSDEwithadditivenoise...... 264 Variationofconstants ...... 264 SolutionsasGaussianprocesses ...... 266 TheOrnstein-Uhlenbeckprocess ...... 268 Changeoftime-scale ...... 271
StochasticAnalysis–AnIntroduction Prof.AndreasEberle CONTENTS 7
8.4 Brownianbridge ...... 273 Wiener-Lévyconstruction ...... 273 Finite-dimensionaldistributions ...... 274 SDEfortheBrownianbridge ...... 276 8.5 Stochastic differential equations in Rn ...... 278 Itô processes driven by several Brownian motions ...... 279 MultivariateItô-Doeblinformula ...... 279 GeneralOrnstein-Uhlenbeckprocesses ...... 281 Examples ...... 281
University of Bonn Winter Term 2010/2011 Chapter 1
Brownian Motion
This introduction to stochastic analysis starts with an introduction to Brownian motion.
Brownian Motion is a diffusion process, i.e. a continuous-time Markov process (Bt)t 0 ≥ with continuous sample paths t B (ω). In fact, it is the only nontrivial continuous- → t time process that is a Lévy process as well as a martingale and a Gaussian process. A rigorous construction of this process has been carried out first by N. Wiener in 1923. Already about 20 years earlier, related models had been introduced independently for financial markets by L. Bachelier [Théorie de la spéculation, Ann. Sci. École Norm. Sup. 17, 1900], and for the velocity of molecular motion by A. Einstein [Über die von der molekularkinetischen Theorie der Wärme geforderte Bewegung von in ruhenden Flüssigkeiten suspendierten Teilchen, Annalen der Physik 17, 1905]. It has been a groundbreaking approach of K. Itô to construct general diffusion processes from Brownian motion, cf. [. . . ]. In classical analysis, the solution of an ordinary dif- ferential equation x′(t) = f(t, x(t)) is a function, that can be approximated locally for t close to t by the linear function x(t )+ f(t , x(t )) (t t ). Similarly, Itô showed, 0 0 0 0 − 0 that a diffusion process behaves locally like a linear function of Brownian motion – the connection being described rigorously by a stochastic differential equation (SDE).
The fundamental rôle played by Brownian motion in stochastic analysis is due to the central limit Theorem. Similarly as the normal distribution arises as a universal scal- ing limit of standardized sums of independent, identically distributed, square integrable
8 1.1. FROM RANDOM WALKS TO BROWNIAN MOTION 9 random variables, Brownian motion shows up as a universal scaling limit of Random Walks with square integrable increments.
1.1 From Random Walks to Brownian Motion
To motivate the definition of Brownian motion below, we first briefly discuss discrete- time stochastic processes and possible continuous-time scaling limits on an informal level.
A standard approach to model stochastic dynamics in discrete time is to start from a se- quence of random variables η , η ,... defined on a common probability space (Ω, ,P ). 1 2 A The random variables ηn describe the stochastic influences (noise) on the system. Often they are assumed to be independent and identically distributed (i.i.d.). In this case the collection (ηn) is also called a white noise, where as a colored noise is given by depen- d dent random variables. A stochastic process Xn, n =0, 1, 2,..., taking values in R is then defined recursively on (Ω, ,P ) by A
Xn+1 = Xn + Φn+1(Xn, ηn+1), n =0, 1, 2,.... (1.1.1)
Here the Φn are measurable maps describing the random law of motion. If X0 and
η1, η2,... are independent random variables, then the process (Xn) is a Markov chain with respect to P .
Now let us assume that the random variables ηn are independent and identically dis- tributed taking values in R, or, more generally, Rd. The easiest type of a nontrivial n stochastic dynamics as described above is the Random Walk Sn = ηi which satisfies i=1 Sn+1 = Sn + ηn+1 for n =0, 1, 2,....
Since the noise random variables ηn are the increments of the Random Walk (Sn), the law of motion (1.1.1) in the general case can be rewritten as
X X = Φ (X ,S S ), n =0, 1, 2,.... (1.1.2) n+1 − n n+1 n n+1 − n
This equation is a difference equation for (Xn) driven by the stochastic process (Sn).
University of Bonn Winter Term 2010/2011 10 CHAPTER 1. BROWNIAN MOTION
Our aim is to carry out a similar construction as above for stochastic dynamics in con- tinuous time. The stochastic difference equation (1.1.2) will then eventually be replaced by a stochastic differential equation (SDE). However, before even being able to think about how to write down and make sense of such an equation, we have to identify a continuous-time stochastic process that takes over the rôle of the Random Walk. For this purpose, we first determine possible scaling limits of Random Walks when the time steps tend to 0. It will turn out that if the increments are square integrable and the size of the increments goes to 0 as the length of the time steps tends to 0, then by the Central Limit Theorem there is essentially only one possible limit process in continuous time: Brownian motion.
Central Limit Theorem
Suppose that Y : Ω Rd, 1 i n < , are identically distributed, square- n,i → ≤ ≤ ∞ integrable random variables on a probability space (Ω, ,P ) such that Y ,...,Y A n,1 n,n are independent for each n N. Then the rescaled sums ∈ 1 n (Y E[Y ]) √n n,i − n,i i=1 converge in distribution to a multivariate normal distribution N(0,C) with covariance matrix (k) (l) Ckl = Cov[Yn,i ,Yn,i ]. To see, how the CLT determines the possible scaling limits of Random Walks, let us consider a one-dimensional Random Walk n
Sn = ηi, n =0, 1, 2,..., i=1 on a probability space (Ω, ,P ) with independent increments η 2(Ω, ,P ) nor- A i ∈ L A malized such that
E[ηi] = 0 and Var[ηi] = 1. (1.1.3)
Plotting many steps of the Random Walk seems to indicate that there is a limit process with continuous sample paths after appropriate rescaling:
StochasticAnalysis–AnIntroduction Prof.AndreasEberle 1.1. FROMRANDOMWALKSTOBROWNIANMOTION 11
4 10
2 5
5 10 15 20 20 40 60 80 100 2 5 − −
50 100
25 50
500 1000 5000 10000 25 50 − −
To see what appropriate means, we fix a positive integer m, and try to define a rescaled (m) Random Walk St (t =0, 1/m, 2/m, . . .) with time steps of size 1/m by
S(m) = c S (k =0, 1, 2,...) k/m m k
for some constants cm > 0. If t is a multiple of 1/m, then
Var[S(m)] = c2 Var[S ] = c2 m t. t m mt m
Hence in order to achieve convergence of S(m) as m , we should choose c t → ∞ m 1/2 (m) proportional to m− . This leads us to define a continuous time process (St )t 0 by ≥
1 S(m)(ω) := S (ω) whenever t = k/m for some integer k, t √m mt
k 1 k and by linear interpolation for t − , . ∈ m m University of Bonn Winter Term 2010/2011 12 CHAPTER 1. BROWNIAN MOTION
(m) St
1
√m
1 2 t m Figure 1.1: Rescaling of a Random Walk.
Clearly, E[S(m)] = 0 for all t 0, t ≥ and 1 Var[S(m)] = Var[S ] = t t m mt whenever t is a multiple of 1/m. In particular, the expectation values and variances for a fixed time t do not depend on m. Moreover, if we fix a partition 0 t < t <...