BROWNIAN MOTION Antonella Basso, Martina Nardon

Similar documents
An Introduction to Stochastic Calculus

An Introduction to Stochastic Calculus

Modeling via Stochastic Processes in Finance

S t d with probability (1 p), where

Continuous Processes. Brownian motion Stochastic calculus Ito calculus

Homework Assignments

Drunken Birds, Brownian Motion, and Other Random Fun

A No-Arbitrage Theorem for Uncertain Stock Model

AMH4 - ADVANCED OPTION PRICING. Contents

Beyond the Black-Scholes-Merton model

Stochastic Dynamical Systems and SDE s. An Informal Introduction

3.1 Itô s Lemma for Continuous Stochastic Variables

1 Mathematics in a Pill 1.1 PROBABILITY SPACE AND RANDOM VARIABLES. A probability triple P consists of the following components:

Non-semimartingales in finance

1 IEOR 4701: Notes on Brownian Motion

IEOR E4703: Monte-Carlo Simulation

BROWNIAN MOTION II. D.Majumdar

IEOR 3106: Introduction to OR: Stochastic Models. Fall 2013, Professor Whitt. Class Lecture Notes: Tuesday, September 10.

Monte Carlo Methods in Financial Engineering

Continuous random variables

Stochastic calculus Introduction I. Stochastic Finance. C. Azizieh VUB 1/91. C. Azizieh VUB Stochastic Finance

STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL

Stochastic Calculus, Application of Real Analysis in Finance

Chapter 3: Black-Scholes Equation and Its Numerical Evaluation

The Black-Scholes Model

Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models

The Black-Scholes Model

Practical example of an Economic Scenario Generator

Continuous Time Finance. Tomas Björk

Introduction to Stochastic Calculus With Applications

Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.

Optimal stopping problems for a Brownian motion with a disorder on a finite interval

1.1 Basic Financial Derivatives: Forward Contracts and Options

Geometric Brownian Motions

1 The continuous time limit

IEOR E4703: Monte-Carlo Simulation

Dr. Maddah ENMG 625 Financial Eng g II 10/16/06

Are stylized facts irrelevant in option-pricing?

MAS3904/MAS8904 Stochastic Financial Modelling

Lévy models in finance

An Introduction to Stochastic Calculus

Reading: You should read Hull chapter 12 and perhaps the very first part of chapter 13.

Optimum Thresholding for Semimartingales with Lévy Jumps under the mean-square error

Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies

Quasi-Monte Carlo for Finance

The stochastic calculus

Equivalence between Semimartingales and Itô Processes

Pricing Dynamic Solvency Insurance and Investment Fund Protection

No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate

Lecture 1: Lévy processes

Simulating Stochastic Differential Equations

Randomness and Fractals

EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS

Math 416/516: Stochastic Simulation

Numerical schemes for SDEs

Functional vs Banach space stochastic calculus & strong-viscosity solutions to semilinear parabolic path-dependent PDEs.

23 Stochastic Ordinary Differential Equations with Examples from Finance

King s College London

Partial differential approach for continuous models. Closed form pricing formulas for discretely monitored models

A new approach for scenario generation in risk management

Lattice (Binomial Trees) Version 1.2

Ornstein-Uhlenbeck Theory

MATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS

The rth moment of a real-valued random variable X with density f(x) is. x r f(x) dx

NORMAL APPROXIMATION. In the last chapter we discovered that, when sampling from almost any distribution, e r2 2 rdrdϕ = 2π e u du =2π.

The value of foresight

UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions.

The Mathematics of Currency Hedging

Business Statistics 41000: Probability 3

Stochastic Calculus - An Introduction

TEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING

VOLATILITY FORECASTING IN A TICK-DATA MODEL L. C. G. Rogers University of Bath

12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006.

I Preliminary Material 1

Asymptotic Methods in Financial Mathematics

Logarithmic derivatives of densities for jump processes

M5MF6. Advanced Methods in Derivatives Pricing

Mgr. Jakub Petrásek 1. May 4, 2009

Short-time asymptotics for ATM option prices under tempered stable processes

10. Monte Carlo Methods

Market Risk Analysis Volume I

Universität Regensburg Mathematik

Math 239 Homework 1 solutions

Stochastic Processes and Financial Mathematics (part two) Dr Nic Freeman

Deterministic Income under a Stochastic Interest Rate

CHOICE THEORY, UTILITY FUNCTIONS AND RISK AVERSION

Module 10:Application of stochastic processes in areas like finance Lecture 36:Black-Scholes Model. Stochastic Differential Equation.

Convergence Analysis of Monte Carlo Calibration of Financial Market Models

A Numerical Approach to the Estimation of Search Effort in a Search for a Moving Object

Stochastic Differential equations as applied to pricing of options

MSc Financial Engineering CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL. To be handed in by monday January 28, 2013

Chapter 3 Common Families of Distributions. Definition 3.4.1: A family of pmfs or pdfs is called exponential family if it can be expressed as

Interest rate models in continuous time

Lecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing

Martingales, Part II, with Exercise Due 9/21

Forecasting Life Expectancy in an International Context

Weak Convergence to Stochastic Integrals

CRRAO Advanced Institute of Mathematics, Statistics and Computer Science (AIMSCS) Research Report. B. L. S. Prakasa Rao

TN 2 - Basic Calculus with Financial Applications

Monte Carlo Simulation of Stochastic Processes

Transcription:

BROWNIAN MOTION Antonella Basso, Martina Nardon basso@unive.it, mnardon@unive.it Department of Applied Mathematics University Ca Foscari Venice Brownian motion p. 1

Brownian motion Brownian motion plays a central role in probability theory, the theory of stochastic processes, physics, economics and finance,... Brownian motion p. 2

Definition of Brownian motion A stochastic process W = (W t ) t 0 is called (standard) Brownian motion or a Wiener process if the following conditions are satisfied: W 0 = 0 (the process starts at zero); it has stationary, independent increments; for every t > 0, W t has a Gaussian N(0,t) distribution; it has continuous sample paths (no jumps ). Brownian motion p. 3

Figure 1.3.1 Insert figure Brownian motion p. 4

A brief history of Brownian motion Brownian motion is named after the biologist Robert Brown whose research dates to the 1820s. Early in the last century, Louis Bachelier (1900), Albert Einstein (1905) and Norbert Wiener (1923) began developing the mathematical theory of Brownian motion. The construction of Bachelier (1900) was erroneous but captured many of the essential properties of the process. Wiener (1923) was the first to put Brownian motion on a firm mathematical basis. Brownian motion p. 5

Distribution a of Brownian motion The finite-dimensional distributions of Brownian motion are multivariate Gaussian; W is a Gaussian process. Brownian motion has independent Gaussian increments. Brownian motion p. 6

Distribution a of Brownian motion The random variables W t W s and W t s have a Gaussian N(0, t s) distribution for s < t. This follows from the stationarity of the increments. W t W s has the same distribution as W t s W 0 = W t s which is normal with zero mean and variance t s. Brownian motion p. 7

Distribution a of Brownian motion The variance of W t W s is proportional to the length of the interval [s, t]. The larger the interval, the larger the fluctuations of Brownian motion on this interval. Brownian motion p. 8

Figure 1.3.2 Insert figure Brownian motion p. 9

Remark The distribution identity W t W s = d W t s does not imply pathwise identity: in general, W t (ω) W s (ω) W t s (ω). Brownian motion p. 10

Brownian motion vs. Poisson process The definitions of Brownian motion and the Poisson process coincide insofar that they are processes with stationary, independent increments. The crucial difference is the kind of distribution of the increments. The requirement of the Poisson distribution makes the sample path pure jumps functions, whereas the Gaussian assumption makes the sample paths continuous. Brownian motion p. 11

Expectation and covariance functions of Brownian motion The Brownian motion has expectation function µ W (t) = E(W t ) = 0 t 0. Since the increments W s W 0 = W s and W t W s are independent for s < t, it has covariance function cov W (s,t) = E [ [(W t W s ) + W s ]W s ] = E [ (W t W s )W s ] + E [ (W s ) 2] 0 s < t = E(W t W s ) E(W s ) + s = 0 + s = s. Brownian motion p. 12

An alternative definition of Brownian motion Since a Gaussian process is characterized by its expectation and covariance functions, we can give an alternative definition. Brownian motion is a Gaussian process with and µ W (t) = E(W t ) = 0 cov W (s,t) = min(s,t). What can we say about variance? Brownian motion p. 13

Path properties In what follows, we fix one sample path W t (ω), t 0 (which is a function of t), and consider its properties. From the definition of Brownian motion, we know that its sample paths are continuous. Brownian paths are extremely irregular: they oscillate wildly. The main reason is that the increments of W are independent. In particular, increments of Brownian motion on adjacent intervals are independent whatever the length of the intervals. The continuity of the paths is rather surprisingly! Brownian motion p. 14

Question How irregular is a Brownian sample path? In order to answer this question, we need to introduce a class of stochastic processes which contains Brownian motion as a special case. Brownian motion p. 15

Self-similar stochastic processes A stochastic process X = (X t ) t 0 is H-self-similar, for some H > 0, if its finite-dimensional distributions satisfy the condition (T H X t1,..., T H X tn ) = d (X T t1,..., X T tn ), (1) for every T > 0, any choice of t i 0, i = 1,...,n, and n 1. Brownian motion p. 16

Remark Note that self-similarity is a distributional, not a pathwise property. In (1), one must not replace = d with =. Roughly speaking, self-similarity means that the properly scaled patterns of a sample path in any small or large time interval have similar shape, but they are not identical. Brownian motion p. 17

Non-Differentiability The sample paths of a self-similar process are nowhere differentiable. Brownian motion p. 18

Non-Differentiability Brownian motion is 1 2 -self-similar: (T 1/2 W t1,..., T 1/2 W tn ) = d (W T t1,..., W T tn ), (2) for every T > 0, any choice of t i 0, i = 1,...,n, and n 1. Its sample paths are nowhere differentiable. Brownian motion p. 19

Self-similarity The left-hand side of (2) are Gaussian vectors, therefore in order to check the distributional identity it suffices to verify that they have the same expectation and covariance matrix. Brownian motion p. 20

Differentiability of a function The graph of a differentiable function is smooth. If the limit f (x 0 ) = lim x 0 f(x 0 + x) f(x 0 ) x exists and is finite for some x 0 (0,t), say, then we may write for small x f(x 0 + x) = f(x 0 ) + f (x 0 ) x + h(x 0, x) x, where h(x 0, x) 0 as x 0. Brownian motion p. 21

Differentiability of a function Hence, in a small neighborhood of x 0, the function f is roughly linear (as a function of x). This explains its smoothness. Alternatively, differentiability of f at x 0 implies that there exists a unique tangent to the curve of the function at this point. Brownian motion p. 22

Figure 1.3.3 Insert figure Brownian motion p. 23

Question How is the graph of a nowhere differentiable function? The graph of such a function changes its shape in the neighborhood of any point in a completely non-predictable manner. It is very difficult to imagine a nowhere differentiable function, nevertheless, Brownian motion is considered as a very good approximation to many real-life phenomena. Brownian motion p. 24

Self-Similarity The self-similarity property of Brownian motion has an interesting consequence for the simulation of its sample paths. In order to simulate a path on [0, T] it suffices to simulate one path on [0, 1], then scale the time interval by the factor T and the sample path by the factor T 1/2. Brownian motion p. 25

Non-Differentiability The existence of nowhere differentiable continuous functions was discovered in the 19th century. One such function was constructed by Weierstrass. It was considered as a curiosity, far away from any practical application. Brownian motion is a process with nowhere differentiable sample paths. Currently it is considered as one of those processes which have a multitude of applications in very different fields. One of them is stochastic calculus and it is used, in particular, in mathematical finance. Brownian motion p. 26

Unbounded Variation A further indication of the irregularity of Brownian sample paths is the concept of unbounded variation. Brownian sample paths do not have bounded variation on any finite interval [0, T]. This means that sup τ n i=1 Wti (ω) W ti 1 (ω) = where the supremum is taken over all possible partitions τ : 0 = t 0 < t 1 < < t n = T of [0, T]. Brownian motion p. 27

Supremum The notation sup n a n = a means the supremum of a sequence of real numbers a n. If a R, a a n for all n and for every ε > 0 there exists a k such that a ε < a k. If a =, then for every M > 0 there exists a k such that a k > M. For a finite index set I, sup n I a n = max n I a n. Brownian motion p. 28

Observation The unbounded variation and non-differentiability of Brownian motion sample paths are major reasons for the failure of classical integration methods, when applied to these paths, and for the introduction of stochastic calculus. Brownian motion p. 29

Processes derived from Brownian motion Various Gaussian and non-gaussian stochastic processes of practical relevance can be derived from Brownian motion. We introduce some of those processes which will find further interesting applications in finance. As before, W = (W t ) t 0 denotes Brownian motion. Brownian motion p. 30

Brownian bridge Consider the process X t = W t t W 1 0 t 1. Obviously, and X 0 = W 0 0 W 1 = 0 X 1 = W 1 1 W 1 = 0. For this simple reason, the process X is called (standard) Brownian bridge or tied down Brownian motion. Brownian motion p. 31

Figure 1.3.4 Insert figure Brownian motion p. 32

Brownian bridge Using the formula for linear transformations of Gaussian random vectors, one can show that the finite-dimensional distributions of X are Gaussian. (Exercise) Brownian motion p. 33

Linear transformation of Gaussian random vectors Let X = (X 1,...,X n ) have a Gaussian N(µ, Σ) distribution and A be an m n matrix. Then AX has a Gaussian N(A µ,aσa ). Note that µ R n and Σ is an n n matrix. Brownian motion p. 34

Distribution, Expectation and Covariance of the Brownian bridge The Brownian bridge X is a Gaussian process. One can easily calculate the expectation and covariance functions of the Brownian bridge: and µ X (t) = E(X t ) = 0 t [0, 1] cov X (s,t) = min(s,t) s t s,t [0, 1]. Since X is Gaussian, the Brownian bridge is characterized by these two functions. Brownian motion p. 35

Brownian motion with drift Let us consider the process X t = µ t + σ W t t 0, for constants σ > 0 and µ R. X is a Gaussian process, with expectation and covariance and respectively. µ X (t) = E(X t ) = µ t t 0 cov X (s,t) = σ 2 min(s,t) s,t 0, Brownian motion p. 36

The drift The expectation function µ X (t) = µ t essentially determines the characteristic shape of the sample paths. The term µ t is the deterministic drift of the process X. Therefore X is called Brownian motion with (linear) drift. Brownian motion p. 37

Figure 1.3.6 Insert figure Brownian motion p. 38

Prices of risky assets With the fundamental discovery of Bachelier (1900) that prices of risky assets (stock indices, exchange rates, share prices, etc.) can be well described by Brownian motion, a new area of applications of stochastic processes was born. However, Brownian motion, as a Gaussian process, may assume negative values, which is not a very desirable property of a price. Brownian motion p. 39

Geometric Brownian motion Consider the process X t = e µ t+σw t t 0, i.e. it is the exponential of Brownian motion with drift. The process X is called geometric Brownian motion. X is not a Gaussian process. Brownian motion p. 40

Expectation of a lognormal random variable Let Z N(0, 1), then E [ e λ Z] = e λ2 /2 λ R. (3) Brownian motion p. 41

Proof E [ e λ Z] = 1 (2π) 1/2 + = e λ2 /2 1 (2π) 1/2 = e λ2 /2. + e λ z e z2 /2 dz e (z λ)2 /2 dz Here we used the fact that 1 (2 π) 1/2 e (z λ)2 is the density of a Gaussian N(λ, 1) random variable. Brownian motion p. 42

Expectation function of geometric Brownian motion From (3) and the self-similarity of Brownian motion, it follows that µ X (t) = e µ t E [ e ] σ W t ] = e µ t E [e σt1/2 W 1 (4) = e (µ+1/2σ2 )t. Brownian motion p. 43

Covariance function of geometric Brownian motion For s t, W t W s and W s are independent, and W t W s = d W t s. Hence cov X (s,t) = e (µ+1/2 σ2 )(s+t) ( ) e σ2s 1. (5) Brownian motion p. 44

Proof cov X (s,t) = E [X s X t ] E [X s ] E [X t ] = e µ(s+t) E [ e σ(w s+w t ) ] e (µ+1/2 σ2 )(s+t) = e µ(s+t) E [ e σ[(w t W s )+2W s ] ] e (µ+1/2 σ2 )(s+t) = e µ(s+t) E [ e ] σ(w t W s ) E [ e ] 2σW s e (µ+1/2 σ2 )(s+t) ( ) = e (µ+1/2 σ2 )(s+t) e σ2s 1. Brownian motion p. 45

Variance function of geometric Brownian motion In particular, geometric Brownian motion has variance function ( ) σx(t) 2 = cov X (t,t) = e (2µ+σ2 )t e σ2t 1. (6) Brownian motion p. 46

Figure 1.3.9 Insert figure Brownian motion p. 47