Stochastic calculus

Stochastic calculusMain

Community hub

Stochastic calculus

7 pages, 0 posts

0 subscribers

Recent from talks

Be the first to start a discussion here.

Recent from talks

Be the first to start a discussion here.

Contribute something

About hubMembersContent overviewUpdatesRules

Main reference articles

Stochastic calculus

View on Wikipedia

from Wikipedia

This article includes a list of references, related reading, or external links, but its sources remain unclear because it lacks inline citations. Please help improve this article by introducing more precise citations. (August 2011) (Learn how and when to remove this message)

Differential

Definitions
Derivative (generalizations) Differential infinitesimal of a function total
Concepts
Differentiation notation Second derivative Implicit differentiation Logarithmic differentiation Related rates Taylor's theorem
Rules and identities
Sum Product Chain Power Quotient L'Hôpital's rule Inverse General Leibniz Faà di Bruno's formula Reynolds

Integral

Definitions
Lists of integrals Integral transform Leibniz integral rule
Antiderivative Integral (improper) Riemann integral Lebesgue integration Contour integration Integral of inverse functions
Integration by
Parts Discs Cylindrical shells Substitution (trigonometric, tangent half-angle, Euler) Euler's formula Partial fractions (Heaviside's method) Changing order Reduction formulae Differentiating under the integral sign Risch algorithm

Series

Convergence tests
Geometric (arithmetico-geometric) Harmonic Alternating Power Binomial Taylor
Summand limit (term test) Ratio Root Integral Direct comparison Limit comparison Alternating series Cauchy condensation Dirichlet Abel

Stochastic calculus is a branch of mathematics that operates on stochastic processes. It allows a consistent theory of integration to be defined for integrals of stochastic processes with respect to stochastic processes. This field was created and started by the Japanese mathematician Kiyosi Itô during World War II.

The best-known stochastic process to which stochastic calculus is applied is the Wiener process (named in honor of Norbert Wiener), which is used for modeling Brownian motion as described by Louis Bachelier in 1900 and by Albert Einstein in 1905 and other physical diffusion processes in space of particles subject to random forces. Since the 1970s, the Wiener process has been widely applied in financial mathematics and economics to model the evolution in time of stock prices and bond interest rates.

The main flavours of stochastic calculus are the Itô calculus and its variational relative the Malliavin calculus. For technical reasons the Itô integral is the most useful for general classes of processes, but the related Stratonovich integral is frequently useful in problem formulation (particularly in engineering disciplines). The Stratonovich integral can readily be expressed in terms of the Itô integral, and vice versa. The main benefit of the Stratonovich integral is that it obeys the usual chain rule and therefore does not require Itô's lemma. This enables problems to be expressed in a coordinate system invariant form, which is invaluable when developing stochastic calculus on manifolds other than Rⁿ. The dominated convergence theorem does not hold for the Stratonovich integral; consequently it is very difficult to prove results without re-expressing the integrals in Itô form.

Itô integral

[edit]

The Itô integral is central to the study of stochastic calculus. The integral $\int H\,dX$ is defined for a semimartingale X and locally bounded predictable process H. ^{[citation needed]}

Stratonovich integral

[edit]

The Stratonovich integral or Fisk–Stratonovich integral of a semimartingale $X$ against another semimartingale Y can be defined in terms of the Itô integral as

\int _{0}^{t}X_{s-}\circ dY_{s}:=\int _{0}^{t}X_{s-}dY_{s}+{\frac {1}{2}}\left[X,Y\right]_{t}^{c},

where [X, Y]_t^c denotes the optional quadratic covariation of the continuous parts of X and Y, which is the optional quadratic covariation minus the jumps of the processes $X$ and $Y$ , i.e.

\left[X,Y\right]_{t}^{c}:=[X,Y]_{t}-\sum \limits _{s\leq t}\Delta X_{s}\Delta Y_{s}

The alternative notation

\int _{0}^{t}X_{s}\,\partial Y_{s}

is also used to denote the Stratonovich integral.

Applications

[edit]

An important application of stochastic calculus is in mathematical finance, in which asset prices are often assumed to follow stochastic differential equations. For example, the Black–Scholes model prices options as if they follow a geometric Brownian motion, illustrating the opportunities and risks from applying stochastic calculus.

Stochastic integrals

[edit]

Besides the classical Itô and Fisk–Stratonovich integrals, many other notions of stochastic integrals exist, such as the Hitsuda–Skorokhod integral, the Marcus integral, and the Ogawa integral.

References

[edit]

Thomas Mikosch, 1998, Elementary Stochastic Calculus, World Scientific, ISBN 981-023543-7
Fima C Klebaner, 2012, Introduction to Stochastic Calculus with Application (3rd Edition). World Scientific Publishing, ISBN 9781848168312
Szabados, T.S.; Székely, B.Z. (2008). "Stochastic Integration Based on Simple, Symmetric Random Walks". Journal of Theoretical Probability. 22: 203–219. arXiv:0712.3908. doi:10.1007/s10959-007-0140-8. Preprint

Industrial and applied mathematics

Computational

Mathematical software	Arbitrary-precision arithmetic Finite element analysis Tensor software Interactive geometry software Optimization software Statistical software Numerical-analysis software Numerical libraries Solvers

Algebraic structures	Algebra of physical space Particle physics and representation theory Feynman integral Poisson algebra Quantum group Renormalization group Spacetime algebra Superalgebra Supersymmetry algebra

Decision sciences

Other applications

Mathematics

Organizations

Category
Mathematics portal / outline / topics list

Revisions and contributors Edit on Wikipedia Read on Wikipedia

View on Grokipedia

from Grokipedia

Stochastic calculus is a branch of mathematics that extends classical calculus to handle stochastic processes, particularly through the development of integrals and differential equations involving random phenomena such as Brownian motion.^[1] It provides essential tools for modeling and analyzing systems influenced by uncertainty, with foundational concepts including the Itô integral and stochastic differential equations (SDEs).^[2] The origins of stochastic calculus trace back to the late 19th and early 20th centuries, beginning with Thorvald Nicolai Thiele's 1880 modeling of Brownian motion for time series analysis.^[3] Louis Bachelier's 1900 thesis applied Brownian motion to stock price fluctuations, introducing the idea of independent, normally distributed increments in financial markets.^[3] Albert Einstein's 1905 physical interpretation of Brownian motion further solidified its theoretical basis, while Norbert Wiener's 1923 rigorous construction using measure theory formalized the Wiener process.^[3] Andrey Kolmogorov's 1931 work on Markov processes connected diffusions to partial differential equations, laying groundwork for later developments.^[3] The pivotal advancement came in 1944 with Kiyosi Itô's introduction of stochastic integration, followed by his 1951 formulation of Itô's lemma, which enables differentiation of stochastic processes and is central to solving SDEs.^[3] By the 1960s and 1970s, contributions from Paul-André Meyer and others, including the Doob-Meyer decomposition in 1962 and the concept of semimartingales in 1970, broadened the theory beyond Markov processes, establishing stochastic calculus as a robust framework.^[3] At its core, stochastic calculus revolves around stochastic processes, which are collections of random variables evolving over time, defined on a probability space

(\Omega, \mathcal{F}, P)

with a filtration

\{\mathcal{F}_t\}

representing accumulating information.^[4] The Wiener process, or Brownian motion, serves as the canonical example: a continuous-time process with independent, normally distributed increments

W_t - W_s \sim \mathcal{N}(0, t-s)

for

t > s

, exhibiting quadratic variation

[W, W]_t = t

.^[2] The Itô integral,

\int_0^t h_s \, dW_s

, extends integration to non-deterministic integrands

h

adapted to the filtration, defined as an

L^2

limit for simple processes and possessing properties like zero mean (

\mathbb{E}[\int h \, dW] = 0

) and Itô isometry (

\mathbb{E}[(\int h \, dW)^2] = \int \mathbb{E}[h^2] \, ds

).^[2] This integral underpins stochastic differential equations of the form

dX_t = b(t, X_t) \, dt + \sigma(t, X_t) \, dW_t

, where

b

is the drift and

\sigma

the diffusion coefficient, solved using Itô's lemma—a chain rule analogue that accounts for the quadratic variation of the Wiener process: for

f(t, X_t)

df = \frac{\partial f}{\partial t} dt + \frac{\partial f}{\partial x} dX + \frac{1}{2} \frac{\partial^2 f}{\partial x^2} (dX)^2

, with

(dX)^2 = \sigma^2 dt

.^[5] Applications of stochastic calculus span multiple fields, most notably mathematical finance, where it derives the Black-Scholes equation for option pricing under the geometric Brownian motion model

dS_t = \mu S_t \, dt + \sigma S_t \, dW_t

, enabling risk-neutral valuation and hedging strategies.^[1] In physics and engineering, it models particle diffusion, noise in signal processing, and optimal control problems, such as minimizing costs in SDEs via the Hamilton-Jacobi-Bellman equation.^[4] Filtering theory, including the Kalman filter extension to nonlinear cases, uses stochastic calculus for state estimation in noisy environments, as in tracking or navigation systems.^[4] Additionally, it connects to partial differential equations through Feynman-Kac representations, linking SDEs to solutions of parabolic PDEs like the heat equation.^[2] These tools have transformed quantitative modeling, emphasizing martingales—processes with constant conditional expectation—for fair pricing and no-arbitrage principles.^[5]

Introduction

Overview

Stochastic calculus is the branch of mathematics that extends the methods of calculus to stochastic processes, particularly emphasizing integration and differentiation in environments characterized by randomness.^[1] It provides tools for analyzing systems where outcomes are probabilistic rather than deterministic, enabling the rigorous treatment of uncertainty in continuous time.^[5] The core components of stochastic calculus include stochastic integrals, which define integration with respect to random processes; stochastic differential equations (SDEs), which model dynamics driven by noise; and Itô's lemma, a fundamental theorem analogous to the chain rule but adapted for stochastic settings.^[1] These elements allow for the manipulation of expressions involving randomness, such as computing expectations or solving equations under uncertainty.^[5] In contrast to deterministic calculus, which assumes smooth and differentiable paths, stochastic calculus addresses the irregularities of random paths, such as those with infinite variation but finite quadratic variation, exemplified by Brownian motion.^[1] This distinction necessitates new definitions and rules to handle the non-differentiability inherent in noise-driven evolutions.^[5] Stochastic calculus is essential for modeling phenomena with intrinsic randomness, including financial markets where asset prices fluctuate unpredictably and physical processes like particle diffusion.^[1] A representative example is the basic SDE

dX_t = \mu \, dt + \sigma \, dW_t,

where

\mu \, dt

captures the deterministic drift,

\sigma \, dW_t

the random volatility term, and

W_t

the Wiener process representing the noise source.^[1]

Historical development

The foundations of stochastic calculus emerged in the early 20th century, building on probabilistic models of random phenomena. In 1900, Louis Bachelier presented his doctoral thesis Théorie de la spéculation, which modeled stock price fluctuations as an arithmetic Brownian motion, introducing the concept of a continuous-time random walk to financial mathematics for the first time. This work laid an early groundwork for applying stochastic processes to economic systems, though it initially received limited attention. Five years later, Albert Einstein's seminal paper on the Brownian motion of suspended particles provided a physical interpretation through the lens of molecular diffusion, rigorously deriving the mean squared displacement proportional to time and connecting microscopic chaos to macroscopic randomness. The 1920s and 1930s saw mathematical formalization that enabled a rigorous framework for stochastic analysis. Norbert Wiener constructed the Wiener process in 1923, defining Brownian motion as a continuous but nowhere differentiable path in a probabilistic space, which became the canonical model for random fluctuations. This was complemented by Andrey Kolmogorov's 1933 axiomatization of probability theory, which provided the measure-theoretic foundations necessary for handling infinite-dimensional path spaces and ensuring the consistency of stochastic integrals. During World War II, practical needs in electronics and signal processing, such as modeling thermal noise in circuits and radar interference, accelerated research into stochastic processes, influencing the development of tools for noisy dynamical systems. A pivotal breakthrough occurred in 1944 when Kiyosi Itô invented the Itô stochastic integral, motivated by his efforts to solve stochastic differential equations describing physical systems perturbed by Brownian noise, such as turbulence and electronic fluctuations. Itô expanded on this in his 1951 memoir On Stochastic Differential Equations, establishing the calculus for non-anticipating integrands with respect to Brownian motion. In the 1950s, Joseph Doob's Stochastic Processes formalized martingale theory, offering a probabilistic structure essential for convergence and optional sampling in stochastic settings. The 1960s introduced Ruslan Stratonovich's integral, developed around 1961 for applications in physics where ordinary chain rules apply, facilitating modeling in quantum mechanics and control theory. From the 1970s onward, stochastic calculus exploded in applications, particularly in finance and beyond. The 1973 Black-Scholes model for option pricing relied on Itô calculus to derive a partial differential equation for asset prices under geometric Brownian motion, revolutionizing quantitative finance. Later extensions addressed limitations with rougher paths; in the 1990s, Terry Lyons developed rough path theory to generalize stochastic integrals to signals with finite p-variation where p > 2, enabling solutions to equations driven by paths beyond semimartingales. These advancements continue to underpin modern probability and its interdisciplinary uses.

Prerequisites

Stochastic processes

A stochastic process is formally defined as a family of random variables

\{X_t : t \in T\}

, where

T

is an index set typically representing time, and each

X_t

is defined on a common probability space

(\Omega, \mathcal{F}, P)

.^[6] This collection describes the evolution of a random phenomenon over time, with realizations known as sample paths or trajectories. In the context of stochastic calculus, continuous-time processes—where

T = [0, \infty)

or a similar interval—are of primary interest, as they model phenomena like asset prices or physical systems with smooth temporal progression.^[7] Stochastic processes are classified based on key properties, such as the Markov property, which states that the future state depends only on the current state and not on the history, leading to Markov processes.^[8] Another important class is Lévy processes, characterized by stationary and independent increments, starting at zero almost surely, and having right-continuous paths with left limits (càdlàg).^[9] Stationarity refers to the invariance of the process's statistical properties over time shifts; strict stationarity requires the joint distribution of any finite collection of variables to remain unchanged under time translation, while weak (or second-order) stationarity assumes constant mean and autocovariance depending only on the time lag, provided second moments exist.^[10] Independent increments mean that the differences

X_t - X_s

for disjoint intervals

(s, t]

are independent random variables, a property central to processes like Lévy.^[11] Examples illustrate these concepts: the Poisson process, a counting process with independent increments and jumps at random times, models events like arrivals in a queue, where the number of events in an interval follows a Poisson distribution with mean proportional to the interval length.^[12] Gaussian processes, where every finite-dimensional distribution is multivariate normal, provide smooth examples with continuous paths, serving as a bridge to processes without jumps.^[13] Path properties are crucial for analysis; continuity implies no jumps, while càdlàg paths—right-continuous with left limits—accommodate jumps common in financial modeling, ensuring well-defined limits for stochastic integrals.^[14] A fundamental result connecting discrete to continuous processes is that the central limit theorem implies sums of independent random variables, suitably scaled and centered, approximate a Brownian motion in distribution, justifying the use of continuous paths for large-scale random walks.^[15] Brownian motion stands out as a key example of a continuous-time stochastic process with these properties, underpinning much of stochastic calculus.^[16]

Brownian motion

Brownian motion, also known as the Wiener process, originates from the empirical observation of erratic particle movement in fluids, first systematically documented by the Scottish botanist Robert Brown in 1827 while examining pollen grains under a microscope.^[17] This phenomenon was later mathematically formalized by Norbert Wiener in 1923, who provided the first rigorous construction of the process as a continuous-time stochastic model. The standard Brownian motion

W = (W_t)_{t \geq 0}

is defined on a probability space as a stochastic process starting at

W_0 = 0

almost surely, with continuous sample paths with probability 1, independent increments, and normally distributed increments such that for

0 \leq s < t

W_t - W_s \sim \mathcal{N}(0, t - s)

.^[18] The covariance function of the process is given by

\mathbb{E}[W_s W_t] = \min(s, t)

for

s, t \geq 0

, which encapsulates its Gaussian nature and the stationary variance of increments.^[19] Key path properties distinguish Brownian motion: almost every sample path is continuous but nowhere differentiable, reflecting its infinite variation despite bounded quadratic variation defined as

[W]_t = t

.^[20] Additionally, the process exhibits scaling invariance, where for any

c > 0

, the rescaled process satisfies

W_{ct} \stackrel{d}{=} \sqrt{c} \, W_t

Wct=dc

Wt in distribution.^[19] One standard construction of Brownian motion proceeds as the scaling limit of symmetric random walks on the integers, where the position after $n$ steps approximates $W_t$ for $t = n / N$ as the step size $1/\sqrt{N}$

tends to zero. Path regularity, including continuity, follows from the Kolmogorov-Chentsov theorem applied to the moments of increments, ensuring the existence of a continuous modification.^[21] The process is unique up to its version, meaning any two Brownian motions agreeing on finite-dimensional distributions possess indistinguishable laws; the canonical version lives on the space of continuous functions $C[0, \infty)$ , while cadlag versions may be considered on the Skorokhod space $D[0, \infty)$ for broader applications.^[20] Brownian motion serves as the foundational driving noise for stochastic integrals and underlies diffusion processes in probabilistic modeling.^[18]

Martingales and filtrations

In stochastic processes, a filtration

\{\mathcal{F}_t\}_{t \geq 0}

is defined as an increasing family of

\sigma

-algebras on a probability space

(\Omega, \mathcal{F}, P)

, where

\mathcal{F}_s \subseteq \mathcal{F}_t

for all

0 \leq s \leq t

, representing the accumulation of information available up to time

t

.^[22] This structure models the progressive revelation of events in a random system, with

\mathcal{F}_0

containing the initial information and

\mathcal{F}_t

incorporating all observable events by time

t

. A stochastic process

X = (X_t)_{t \geq 0}

is said to be adapted to the filtration if, for each

t

, the random variable

X_t

\mathcal{F}_t

-measurable, meaning its value at time

t

depends only on the information up to that time.^[22] Central to the theory are martingales, which formalize the notion of a "fair game" in probabilistic terms. A stochastic process

M = (M_t)_{t \geq 0}

is a martingale with respect to a filtration

\{\mathcal{F}_t\}_{t \geq 0}

if it is adapted to the filtration,

\mathbb{E}[|M_t|] < \infty

for all

t

, and satisfies the conditional expectation property

\mathbb{E}[M_t \mid \mathcal{F}_s] = M_s

almost surely for all

0 \leq s \leq t

. This property implies that the expected future value of the process, given the current information, equals its current value, capturing predictability in an average sense. The optional sampling theorem extends this by stating that, under suitable conditions such as bounded stopping times

\tau

, the stopped process

M_{\tau \wedge t}

remains a martingale, allowing evaluation at random times without altering the fairness property. Martingales possess several key properties that underpin their utility in analysis. Doob's maximal inequalities provide bounds on the supremum of the process: for a nonnegative submartingale

X

and

p > 1

\mathbb{E}[\sup_{0 \leq s \leq t} X_s^p]^{1/p} \leq \frac{p}{p-1} \mathbb{E}[X_t^p]^{1/p}

, controlling the likelihood of large deviations. The martingale convergence theorem states that if

(M_t)

is a martingale bounded in

L^1

(i.e.,

\sup_t \mathbb{E}[|M_t|] < \infty

), then

M_t

converges almost surely to some

M_\infty \in L^1

t \to \infty

. For convergence in

L^1

, uniform integrability is required: a family

\{|M_t| : t \geq 0\}

is uniformly integrable if

\sup_t \mathbb{E}[|M_t| \mathbf{1}_{\{|M_t| > K\}}] \to 0

K \to \infty

, ensuring the limit preserves the

L^1

norm.^[23] Submartingales and supermartingales generalize martingales through inequalities in the conditional expectation. A process

X

is a submartingale if

\mathbb{E}[X_t \mid \mathcal{F}_s] \geq X_s

a.s. for

s \leq t

, modeling processes with a nonnegative drift, while a supermartingale satisfies

\mathbb{E}[X_t \mid \mathcal{F}_s] \leq X_s

a.s., indicating a nonpositive drift.^[22] These definitions extend the convergence and inequality results from martingales, with submartingales converging almost surely under uniform integrability.^[23] Classic examples illustrate these concepts. Standard Brownian motion

(B_t)_{t \geq 0}

, adapted to its natural filtration, is a martingale because

\mathbb{E}[B_t \mid \mathcal{F}_s] = B_s

for

s \leq t

, reflecting its zero-drift property.^[22] Similarly, the compensated Poisson process

M_t = N_t - \lambda t

, where

N

is a Poisson process with rate

\lambda

and natural filtration, forms a martingale since the compensator

\lambda t

subtracts the expected increments, yielding

\mathbb{E}[M_t \mid \mathcal{F}_s] = M_s

.^[23] In stochastic calculus, filtrations and martingales provide the foundational framework for constructing integrals and ensuring their well-definedness, as adapted integrands and martingale properties allow the resulting processes to maintain predictability and avoid pathological behaviors.^[23] This structure is essential for preserving key probabilistic features in more advanced developments.^[22]

Stochastic Integrals

Itô integral

The Itô integral provides a framework for integrating adapted stochastic processes with respect to Brownian motion, addressing the challenges posed by the irregular paths of the latter. For a standard Brownian motion

W

on a probability space

(\Omega, \mathcal{F}, (\mathcal{F}_t), P)

with the natural filtration

(\mathcal{F}_t)

, consider progressively measurable processes

\phi = \{\phi_t\}_{t \geq 0}

that are square-integrable in the sense that

\mathbb{E}\left[\int_0^T \phi_t^2 \, dt\right] < \infty

for each

T > 0

. The Itô integral

\int_0^t \phi_s \, dW_s

is constructed as the

\mathbb{L}^2(P)

-limit of approximating sums using left-endpoint evaluations to ensure non-anticipating behavior. The construction begins with simple processes, which are linear combinations of indicator functions of the form

\phi_t = \sum_{i=1}^n c_i(\omega) \mathbf{1}_{[t_{i-1}, t_i)}(t)

, where

0 = t_0 < t_1 < \cdots < t_n = t

, the

c_i

are

\mathcal{F}_{t_{i-1}}

-measurable, and

\phi

satisfies the square-integrability condition. For such processes, the Itô integral over

[0, t]

is defined as the stochastic sum

\int_0^t \phi_s \, dW_s = \sum_{i=1}^n c_i (W_{t_i} - W_{t_{i-1}}),

evaluated at the left endpoint of each subinterval to preserve the adaptedness and avoid anticipation of future Brownian increments. For indicator functions specifically, such as

\phi_s = \mathbf{1}_{[0, u]}(s)

with

u \leq t

, the integral simplifies to

W_t - W_u

almost surely. This left-endpoint choice distinguishes the Itô integral from classical Riemann sums, as the quadratic variation of Brownian motion,

[W, W]_t = t

, introduces non-commutativity: interchanging the order of integration with respect to

W

and a smooth function yields an extra term proportional to the quadratic variation, preventing convergence in the Riemann-Stieltjes sense.^[24] The definition extends to the full class of square-integrable progressively measurable processes by completing the space of simple processes in the norm

\|\phi\|^2 = \mathbb{E}\left[\int_0^t \phi_s^2 \, ds\right]

. For a Cauchy sequence

\{\phi^{(n)}\}

of simple processes converging to

\phi

in this norm, the corresponding Itô integrals

\int_0^t \phi^{(n)}_s \, dW_s

form a Cauchy sequence in

\mathbb{L}^2(P)

, converging in

\mathbb{L}^2(P)

to a limit denoted

\int_0^t \phi_s \, dW_s

. This extension is possible due to the Itô isometry,

\mathbb{E}\left[\left( \int_0^t \phi_s \, dW_s \right)^2 \right] = \mathbb{E}\left[ \int_0^t \phi_s^2 \, ds \right],

which holds for simple processes by independence of Brownian increments and extends continuously to the closure, ensuring the space of Itô integrals is complete (a Hilbert space). The Itô integral has zero expectation,

\mathbb{E}\left[\int_0^t \phi_s \, dW_s\right] = 0

, as each approximating sum has mean zero. Key properties include the martingale property: if

\mathbb{E}\left[\int_0^\infty \phi_s^2 \, ds\right] < \infty

, then

M_t = \int_0^t \phi_s \, dW_s

is a square-integrable martingale with respect to

(\mathcal{F}_t)

, as the increments are orthogonal to the past filtration. More generally, it is a local martingale. The map

\phi \mapsto \int_0^t \phi_s \, dW_s

is continuous from the space of integrands to

\mathbb{L}^2(P)

, and the resulting process has continuous paths almost surely. While the construction is specific to Brownian motion, the Itô integral extends to integration with respect to semimartingales via a decomposition into a continuous local martingale part (integrated against like Brownian motion) and a finite-variation part (via Stieltjes integration), though the latter is not derived here.

Stratonovich integral

The Stratonovich integral, introduced by Ruslan Stratonovich in the context of stochastic equations driven by random noise, provides a symmetric interpretation of stochastic integration that contrasts with the forward-looking Itô integral.^[25] It is particularly suited for applications where classical calculus rules are desirable, and is defined for progressively measurable integrand processes

\phi

with respect to a semimartingale integrator

X

, such as Brownian motion

W

. The integral

\int_0^t \phi_s \circ dX_s

is constructed as the limit in probability (or pathwise under suitable conditions) of Riemann-Stieltjes sums using evaluations at midpoints of partition intervals.^[25]^[26] Specifically, for a partition

0 = t_0 < t_1 < \cdots < t_n = t

[0, t]

with mesh size approaching zero, the approximating sums are

\sum_{i=0}^{n-1} \phi_{\frac{t_i + t_{i+1}}{2}} (X_{t_{i+1}} - X_{t_i}),

where

\phi

is evaluated at the midpoint

\frac{t_i + t_{i+1}}{2}

.^[25] This midpoint evaluation yields the Stratonovich integral as the limit, first defined for simple step functions and then extended by density arguments to square-integrable or more general adapted processes via approximation in the sup norm or quadratic variation.^[26] Unlike the Itô integral, the Stratonovich integral does not generally produce an

L^2

-martingale but instead defines a semimartingale, allowing integration against processes of finite variation or with controlled quadratic variation.^[27] A key property of the Stratonovich integral is that it satisfies the ordinary chain rule of calculus, enabling straightforward application of transformation formulas without additional correction terms, much like in deterministic integration.^[28] For instance, the integral

\int_0^t W_s \circ dW_s = \frac{1}{2} W_t^2

, mirroring the classical result.^[28] This feature makes it valuable in contexts requiring intuitive calculus manipulations, such as deriving stochastic differential equations (SDEs) with physical interpretability.^[25] The Stratonovich integral arises naturally in physical and engineering models through the Wong–Zakai approximation theorem, which shows that when Brownian motion is approximated by smooth paths (e.g., piecewise linear or mollified versions), the corresponding ordinary Riemann–Stieltjes integrals converge to the Stratonovich integral as the approximation granularity increases.^[26] This convergence holds under mild conditions on the driving noise and integrand, justifying its use in stochastic averaging and systems where noise is modeled as a limit of correlated, smooth perturbations, common in mechanics and signal processing.^[26]^[25] In relation to the Itô integral, the Stratonovich integral can be expressed as

\int_0^t \phi_s \circ dW_s = \int_0^t \phi_s \, dW_s + \frac{1}{2} \langle \phi, W \rangle_t

, where

\langle \phi, W \rangle_t

denotes the quadratic covariation process between

\phi

and

W

.^[28] For continuous semimartingales, this integral is uniquely defined pathwise, independent of the choice of approximating sequence, provided the quadratic variation remains finite.^[27]

Properties and relations

Both the Itô integral and the Stratonovich integral share several fundamental properties as stochastic integrals with respect to Brownian motion. They are both semimartingales, ensuring that they can serve as integrators in further stochastic integrations, and they exhibit continuity in probability, meaning that the processes converge in probability as the time partition refines. Additionally, both integrals are defined for progressively measurable integrand processes, which allows adaptation to the underlying filtration generated by the Brownian motion.^[29] A key property specific to the Itô integral is its quadratic variation. For an Itô integral process

I_t = \int_0^t \phi_s \, dW_s

, where

\phi

is a progressively measurable integrand satisfying

\mathbb{E}\left[\int_0^t \phi_s^2 \, ds\right] < \infty

and

W

is a standard Brownian motion, the quadratic variation is given by

[I]_t = \int_0^t \phi_s^2 \, ds.

This follows from the definition of quadratic variation for martingales and the Itô isometry, which equates the expected squared increment of the integral to the integral of the squared integrand.^[30] The Itô and Stratonovich integrals are related through a precise conversion formula derived from their differing Riemann sum approximations. The Stratonovich integral

\int_0^t \phi_s \circ dW_s

is the limit of midpoint Riemann sums, while the Itô integral

\int_0^t \phi_s \, dW_s

uses left-endpoint sums. The relation is

\int_0^t \phi_s \circ dW_s = \int_0^t \phi_s \, dW_s + \frac{1}{2} [\phi, W]_t,

where

[\phi, W]_t

is the quadratic covariation between

\phi

and

W

. This covariation captures the second-order interaction due to the non-zero quadratic variation of Brownian motion,

[W]_t = t

. For

\phi

an Itô process

d\phi_t = \mu_t \, dt + \sigma_t \, dW_t

, the covariation is

[\phi, W]_t = \int_0^t \sigma_s \, ds

. In specific cases, such as when the diffusion coefficient

\sigma_t = \phi_t

(e.g., in geometric Brownian motion), this becomes

\frac{1}{2} \int_0^t \phi_s \, ds

. The formula follows from the definitions and properties of quadratic covariation, with the correction arising as the mesh of the partition approaches zero.^[31]^[32] The choice between Itô and Stratonovich interpretations often depends on the application domain. In mathematical finance, the Itô integral is preferred because it preserves the martingale property essential for arbitrage-free pricing and risk-neutral valuation, aligning with the non-anticipating nature of market information. In contrast, the Stratonovich integral is favored in physics and engineering for its geometric interpretation, where it satisfies the ordinary chain rule, facilitating modeling of physical systems like diffusion processes with multiplicative noise.^[33] These integrals generalize to semimartingales beyond Brownian motion. The Stratonovich integral extends using the continuous part of the predictable quadratic covariation

\langle X, Y \rangle^c_t

, defined as

\int_0^t \phi_s \circ dX_s = \int_0^t \phi_s \, dX_s + \frac{1}{2} \langle \phi, X \rangle^c_t

, where the conversion mirrors the Brownian case but accounts for the continuous finite variation component. This framework, introduced by Meyer, unifies the theory for processes with jumps while preserving chain rule properties for the continuous part.^[34] A notable relation is provided by the Wong-Zakai theorem, which states that approximating Brownian motion by smooth paths (e.g., piecewise linear interpolations) and computing ordinary Riemann-Stieltjes integrals yields convergence to the Stratonovich integral in the limit of refinement. This justifies the Stratonovich interpretation as a natural extension of deterministic calculus to noisy paths. Finally, the Itô and Stratonovich integrals coincide when the integrand

\phi

is deterministic, as the covariation

[\phi, W]_t = 0

, eliminating the correction term. In stochastic differential equations, this equivalence holds only in such cases; otherwise, the interpretations differ by a diffusion-induced drift adjustment of

\frac{1}{2} \sigma \sigma'

in the Itô form, affecting long-term behavior like stationary distributions.^[29]

Stochastic Differential Equations

Definition and solutions

Stochastic differential equations (SDEs) provide a framework for modeling systems influenced by random noise, extending ordinary differential equations by incorporating stochastic integrals. In the Itô form, an SDE is expressed as

dX_t = b(t, X_t) \, dt + \sigma(t, X_t) \, dW_t,

where

X_t

is the state process,

b: \mathbb{R}_+ \times \mathbb{R}^d \to \mathbb{R}^d

is the drift coefficient,

\sigma: \mathbb{R}_+ \times \mathbb{R}^d \to \mathbb{R}^{d \times m}

is the diffusion coefficient, and

W_t

is an

m

-dimensional Brownian motion.^[35] This differential notation corresponds to the integral equation

X_t = X_0 + \int_0^t b(s, X_s) \, ds + \int_0^t \sigma(s, X_s) \, dW_s,

where the second integral is an Itô stochastic integral.^[35] In contrast, the Stratonovich form uses symmetric integrals, denoted by

\circ dW_t

, which arise naturally in some physical derivations but require conversion to Itô form for standard stochastic calculus tools. Solutions to SDEs are classified as strong or weak. A strong solution is a process

X = (X_t)_{t \geq 0}

adapted to the filtration generated by the driving Brownian motion

W

, satisfying the integral equation pathwise almost surely on a given probability space.^[36] A weak solution, however, exists on some probability space with a Brownian motion

\tilde{W}

(possibly different from

W

) such that

X

and

\tilde{W}

satisfy the integral equation, with the law of

X

matching that required by the equation; every strong solution is weak, but not conversely.^[35] For linear SDEs of the form

dX_t = (a(t) X_t + c(t)) \, dt + (B(t) X_t + D(t)) \, dW_t

, explicit solutions can be obtained using an integrating factor analogous to the deterministic case, yielding

X_t = \Phi(t) \left( X_0 + \int_0^t \Phi(s)^{-1} (c(s) \, ds + D(s) \, dW_s) \right)

, where

\Phi

solves the associated homogeneous equation.^[35] A prominent example is geometric Brownian motion,

dS_t = \mu S_t \, dt + \sigma S_t \, dW_t

, whose solution is

S_t = S_0 \exp\left( \left( \mu - \frac{\sigma^2}{2} \right) t + \sigma W_t \right),

widely used in modeling asset prices.^[35] Existence and uniqueness of solutions rely on conditions on the coefficients. Under global Lipschitz continuity—i.e.,

|b(t,x) - b(t,y)| + |\sigma(t,x) - \sigma(t,y)| \leq K |x - y|

for some

K > 0

and all

t, x, y

—along with linear growth bounds, Picard iteration establishes the existence of a unique strong solution via successive approximations converging in appropriate norms.^[37] These Lipschitz conditions ensure pathwise uniqueness, implying strong existence from weak existence via Yamada-Watanabe theorems.^[36] Solutions to SDEs with time-homogeneous coefficients

b(t,x) = b(x)

and

\sigma(t,x) = \sigma(x)

possess the Markov property: the future distribution of

X_t

given the past depends only on the current state

X_s

for

s < t

.^[35] For numerical approximation, the Euler-Maruyama scheme discretizes the SDE on a grid

0 = t_0 < t_1 < \cdots < t_N = T

Y_{n+1} = Y_n + b(t_n, Y_n) \Delta t_n + \sigma(t_n, Y_n) \Delta W_n

, where

\Delta t_n = t_{n+1} - t_n

and

\Delta W_n = W_{t_{n+1}} - W_{t_n}

, converging weakly to the true solution under Lipschitz and growth conditions.^[38]

Itô's lemma

Itô's lemma provides the stochastic chain rule for functions of Itô processes, accounting for the quadratic variation inherent in stochastic differentials, unlike the classical chain rule for deterministic functions. Consider a twice continuously differentiable function

f(t, x)

and an Itô process

X_t

satisfying the stochastic differential equation (SDE)

dX_t = \mu(t, X_t) \, dt + \sigma(t, X_t) \, dW_t

, where

W_t

is a standard Brownian motion. Then, the differential of

Y_t = f(t, X_t)

is given by

dY_t = \left( \frac{\partial f}{\partial t}(t, X_t) + \mu(t, X_t) \frac{\partial f}{\partial x}(t, X_t) + \frac{1}{2} \sigma^2(t, X_t) \frac{\partial^2 f}{\partial x^2}(t, X_t) \right) dt + \sigma(t, X_t) \frac{\partial f}{\partial x}(t, X_t) \, dW_t.

This formula arises from a second-order Taylor expansion of

f

, truncated at higher-order terms that vanish in the limit. Specifically,

df(t, X_t) = f_t \, dt + f_x \, dX_t + \frac{1}{2} f_{xx} (dX_t)^2 + o(dt)

. Substituting

dX_t = \mu \, dt + \sigma \, dW_t

yields

(dX_t)^2 = \sigma^2 (dW_t)^2 + 2\mu\sigma \, dt \, dW_t + \mu^2 (dt)^2

. The multiplication rules for stochastic differentials simplify this to

(dX_t)^2 = \sigma^2 \, dt

, since

(dW_t)^2 = dt

dt \, dW_t = 0

, and

(dt)^2 = 0

, with the latter following from the quadratic variation of Brownian motion being

\langle W \rangle_t = t

. A formal proof relies on the definition of the Itô integral and stochastic integration by parts. For simple processes, the result holds by direct computation; extension to general cases uses approximation by smooth functions and limits in probability. In the multidimensional setting, let

X_t = (X_t^1, \dots, X_t^d)

be a vector Itô process with

dX_t^i = \mu^i(t, X_t) \, dt + \sum_{k=1}^m \sigma^{i k}(t, X_t) \, dW_t^k

, where

W_t^k

are independent Brownian motions. For

f(t, x)

with continuous second partial derivatives, Itô's lemma states

dY_t = \left( f_t + \sum_{i=1}^d \mu^i f_{x_i} + \frac{1}{2} \sum_{i,j=1}^d \sum_{k=1}^m \sigma^{i k} \sigma^{j k} f_{x_i x_j} \right) dt + \sum_{i=1}^d \sum_{k=1}^m \sigma^{i k} f_{x_i} \, dW_t^k,

where the second-order term incorporates the covariation matrix of the diffusion components. For general semimartingales

X

and

Y

, possibly discontinuous, the full Itô formula for

f \in C^{1,2}

f(t, X_t) - f(0, X_0) = \int_0^t f_s(s, X_{s-}) \, ds + \int_0^t f_x(s, X_{s-}) \, dX_s^c + \frac{1}{2} \int_0^t f_{xx}(s, X_{s-}) \, d\langle X^c \rangle_s + \sum_{0 < s \leq t} \left[ f(s, X_s) - f(s, X_{s-}) - f_x(s, X_{s-}) \Delta X_s \right],

with the continuous parts using covariation

\langle X^c, Y^c \rangle

for vector processes, and the jump term capturing discontinuities via

\Delta X_s = X_s - X_{s-}

. This extends the diffusion case to include finite variation and jump components. In contrast, the Stratonovich integral, defined as

\int H \circ dX = \int H \, dX + \frac{1}{2} \langle H, X \rangle

, yields the classical chain rule:

df(t, X_t) = f_t \, dt + f_x \circ dX_t

. The Itô-Stratonovich conversion adjusts for the quadratic covariation term. Itô's lemma is essential for deriving solutions to SDEs and applications such as derivative pricing in mathematical finance.

Existence and uniqueness

The existence and uniqueness of solutions to stochastic differential equations (SDEs) of Itô type are established under conditions analogous to those in the deterministic Picard-Lindelöf theorem, adapted to account for the stochastic nature of the driving noise. Specifically, suppose the drift coefficient

b(t, x)

and diffusion coefficient

\sigma(t, x)

are locally Lipschitz continuous in

x

uniformly in

t

and satisfy the linear growth condition

|b(t, x)| + |\sigma(t, x)| \leq K(1 + |x|)

for some constant

K > 0

. Then, for a given initial condition

X_0

independent of the Brownian motion

W

, there exists a unique strong solution to the SDE

dX_t = b(t, X_t) dt + \sigma(t, X_t) dW_t

on the interval

[0, \tau)

, where

\tau

is the explosion time, which is almost surely positive. This result, originally due to Itô, relies on Picard iteration applied to the integral form of the SDE, leveraging the contraction properties in appropriate Banach spaces of stochastic processes. If the coefficients

b

and

\sigma

are globally Lipschitz continuous in

x

(i.e.,

|b(t, x) - b(t, y)| + |\sigma(t, x) - \sigma(t, y)| \leq K |x - y|

) and satisfy the linear growth condition, the solution exists and is unique on the entire interval

[0, \infty)

with probability 1, and the explosion time

\tau = \infty

almost surely. These global conditions prevent finite-time explosions, ensuring the solution remains well-defined without reaching infinity in finite time. In contrast, violations of linear growth can lead to explosions; for instance, SDEs with superlinear growth in the coefficients may have solutions that explode in finite time with positive probability.^[36] For cases where strong existence fails, the Yamada-Watanabe theorem provides a bridge between weak and strong solutions: if weak existence holds (i.e., there exists a probability space with a Brownian motion and a process satisfying the SDE in distribution) and pathwise uniqueness holds (i.e., any two solutions driven by the same Brownian motion coincide almost surely), then there exists a unique strong solution. This result, which generalizes earlier work on martingale problems, underscores the role of pathwise uniqueness in upgrading weak solutions to strong ones. A classic counterexample illustrating the necessity of suitable conditions for strong uniqueness is Tanaka's SDE,

dX_t = \operatorname{sign}(X_t) \, dW_t, \quad X_0 = 0,

which admits weak solutions (e.g.,

X_t = |W_t|

in law under a suitable measure change) but no strong solution, as pathwise uniqueness fails—different versions of the absolute value process can be constructed on the same probability space without coinciding. Weak existence can often be established even without Lipschitz conditions on the diffusion coefficient using the Girsanov theorem, which allows construction of a weak solution by starting with a Brownian motion (satisfying a simpler SDE with zero drift) and changing the probability measure via an exponential martingale to incorporate the desired drift, provided the Novikov condition or similar integrability holds. Under the standard global Lipschitz and linear growth assumptions, numerical approximations like the Euler-Maruyama scheme converge strongly to the true solution with order

1/2

, meaning the expected error satisfies

\mathbb{E}[|X_T - X_T^n|] \leq C h^{1/2}

for time step

h

, where

X_T^n

is the numerical approximation at final time

T

. For Stratonovich SDEs, which arise naturally in physical applications due to their chain rule properties, existence and uniqueness theorems follow similar lines but are stated in terms of the Stratonovich coefficients. If the Stratonovich drift

b

and diffusion

\sigma

satisfy global Lipschitz continuity and linear growth (adjusted for the symmetric integral interpretation, often via conversion to an equivalent Itô SDE with an additional drift term

\frac{1}{2} \sum_i \partial_i \sigma_{i j} \sigma_{i j}

), a unique strong global solution exists on

[0, \infty)

. These conditions ensure the Stratonovich integral is well-defined and the resulting process behaves analogously to the Itô case under the specified regularity.

Applications

Mathematical finance

In mathematical finance, stochastic calculus provides the foundational tools for modeling asset prices and derivatives under uncertainty. The Black-Scholes model assumes that the stock price

S_t

follows the stochastic differential equation (SDE)

dS_t = \mu S_t \, dt + \sigma S_t \, dW_t,

where

\mu

is the drift,

\sigma

is the volatility, and

W_t

is a standard Brownian motion.^[39] This model facilitates the pricing of European call options through the risk-neutral measure, obtained via Girsanov's theorem, which changes the probability measure to eliminate the drift's risk premium.^[40] The explicit solution to this SDE is geometric Brownian motion.^[39] Applying Itô's lemma to the option value

V(S_t, t)

yields the Black-Scholes partial differential equation (PDE)

\frac{\partial V}{\partial t} + \frac{1}{2} \sigma^2 S^2 \frac{\partial^2 V}{\partial S^2} + r S \frac{\partial V}{\partial S} - r V = 0,

where

r

is the risk-free rate, under the risk-neutral dynamics.^[39] Solving this PDE gives the 1973 Black-Scholes formula for a European call option:

V = S N(d_1) - K e^{-rT} N(d_2),

with

d_1 = \frac{\ln(S/K) + (r + \sigma^2/2)T}{\sigma \sqrt{T}}

d1=σT

ln(S/K)+(r+σ2/2)T and $d_2 = d_1 - \sigma \sqrt{T}$

, where $N(\cdot)$ is the cumulative distribution function of the standard normal, $K$ is the strike price, and $T$ is the time to maturity.^[39] Hedging in this framework involves constructing a delta-hedging portfolio, where the hedge ratio is $\Delta = \partial V / \partial S$ , and the portfolio remains self-financing through stochastic integrals that replicate the option payoff without arbitrage.^[39] For exotic options, such as barrier options that activate or deactivate upon the underlying price hitting a barrier, pricing often relies on reflected Brownian motion to account for the boundary conditions in the risk-neutral framework.^[41] Extensions of the Black-Scholes model incorporate more realistic dynamics; for instance, Merton's 1976 jump-diffusion model adds Poisson jumps to the SDE to capture sudden price discontinuities.^[42] Similarly, the Heston 1993 model introduces stochastic volatility via an additional SDE for the variance process, allowing correlation between asset returns and volatility shocks.^[43] Risk management in these models employs Value-at-Risk (VaR), which quantifies potential losses at a given confidence level, computed via Monte Carlo simulations of the underlying SDEs to generate future price paths and estimate tail distributions.

Physics and engineering

Stochastic calculus provides essential tools for modeling systems in physics and engineering where random fluctuations, such as thermal noise, play a significant role. In physics, it enables the description of particle dynamics under irregular forces, while in engineering, it supports the design of robust systems for estimation and control amid uncertainty. These applications often involve stochastic differential equations (SDEs) to capture both deterministic evolution and diffusive noise. A foundational example is the Langevin equation, which models the motion of a Brownian particle subject to friction and random kicks from the surrounding medium. For a particle of mass

m

, velocity

v

, friction coefficient

\gamma

, Boltzmann constant

k

, temperature

T

, and Wiener process

W

, the equation is given by

m \, dv = -\gamma v \, dt + \sqrt{2 \gamma k T} \, dW,

mdv=−γvdt+2γkT

dW, where the noise term represents Gaussian white noise with variance tied to temperature. In physical contexts, the Stratonovich interpretation of this SDE is preferred because it preserves the chain rule from ordinary calculus and aligns with the continuous limit of microscopic collisions, ensuring consistency with equilibrium thermodynamics.^[44]^[45] From such SDEs, the Fokker-Planck equation describes the evolution of the probability density $p(x, t)$ of the system's state $x$ . For an Itô SDE $dx = \mu(x, t) \, dt + \sigma(x, t) \, dW$ , the corresponding Fokker-Planck equation is $\frac{\partial p}{\partial t} = -\frac{\partial}{\partial x} \left( \mu p \right) + \frac{1}{2} \frac{\partial^2}{\partial x^2} \left( \sigma^2 p \right),$ derived via Itô's lemma applied to the density evolution. This equation quantifies how noise diffuses probability mass while drift terms shift it, providing insights into stationary distributions and relaxation times in noisy physical systems.^[46] The fluctuation-dissipation theorem connects the strength of these noise fluctuations to dissipative processes and temperature, ensuring that equilibrium is reached with the correct thermal distribution. In the Langevin framework, the noise amplitude $\sqrt{2 \gamma k T}$

directly embodies this relation, where the diffusion coefficient scales with temperature to balance frictional damping, as derived from linear response theory.^[47] In quantum mechanics, stochastic unraveling techniques decompose the Schrödinger equation into ensembles of stochastic trajectories, facilitating simulations of open quantum systems interacting with noisy environments. These methods represent the density matrix evolution as an average over nonlinear stochastic Schrödinger equations driven by quantum noise processes, enabling efficient numerical treatment of decoherence and measurement effects.^[48]^[49] Engineering applications leverage stochastic calculus for state estimation in noisy systems, exemplified by the Kalman-Bucy filter. For a linear system with state $X$ , observation $Y$ , dynamics matrix $A$ , observation matrix $C$ , and gain $K$ , the filter update is $d\hat{X} = A \hat{X} \, dt + K \left( dY - C \hat{X} \, dt \right),$ where stochastic integrals handle the noisy measurements, minimizing estimation error variance in real-time applications like navigation and signal processing.^[50] In control theory, stochastic optimal control addresses decision-making under uncertainty, with the Hamilton-Jacobi-Bellman (HJB) equation providing the value function for optimality. For Itô SDEs with control $u$ , the HJB equation arises from the dynamic programming principle and Itô's lemma, balancing infinitesimal generators of the controlled process to minimize expected costs. This framework is crucial for designing controllers in stochastic environments, such as robotic systems with sensor noise.^[51] Representative examples include the stochastic logistic equation for population dynamics, $dN = r N (1 - N/K) \, dt + \sigma N \, dW$ , which models growth limited by carrying capacity $K$ amid environmental fluctuations, leading to quasi-stationary distributions around the deterministic equilibrium. In electrical engineering, thermal noise in circuits is captured by SDEs for resistor-inductor networks, where Johnson-Nyquist noise drives voltage fluctuations proportional to $\sqrt{4 k T R \Delta f}$

, analyzed via stochastic integrals to predict signal integrity.^[52]^[53]

Authority control databases
International	GND
National	United States Israel
Other	Yale LUX

History

Stochastic calculus

Recent from talks

Recent from talks

Contribute something

Contribute something

Media Pages

Timelines

Articles

Notes collections

Notes

Notes

Days in Chronicle

Stochastic calculus

Itô integral

Stratonovich integral

Applications

Stochastic integrals

See also

References