B8.3 Mathematical Models for Financial Derivatives. Hilary Term Solution Sheet 2

Size: px
Start display at page:

Download "B8.3 Mathematical Models for Financial Derivatives. Hilary Term Solution Sheet 2"

Transcription

1 B8.3 Mathematical Models for Financial Derivatives Hilary Term 18 Solution Sheet In the following W t ) t denotes a standard Brownian motion and t > denotes time. A partition π of the interval, t is a sequence of points = t < t 1 < t < < t n = t and π = max k t k+1 t k ). On a given partition W k W tk, δw k W k+1 W k, δt k t k+1 t k and if f is a function on, t, f k ft k ) and δf k f k+1 f k. 1. Show that if t, s then EW s W t = mins, t). If s = t then we have EW s W t = EW t = t. If not, assume s < t and write W t = W s + W t W s. Then EW s W t = E W s W s + W s Wt W s ) = EW s W s + E W s Wt W s ) = EW s = s = mins, t).. Suppose we define the following two stochastic integrals, the anti-itô integral fw s, s) dw s = lim and the Stratonovich integral π k= fw k+1, t k+1 )W k+1 W k ), fw s, s) dw s = lim π k= 1 fwk+1, t k+1 ) + fw k, t k ) ) W k+1 W k ). Show that W s dw s = Wt + t, E W s dw s = Wt, E W s dw s W s dw s = t, = t 1

2 W s dw s = lim π k= π k= ) W k+1 Wk+1 W k = lim W k+1 Wk ) ) + lim Wk+1 W k = W n W + lim π k= π k= ) Wk+1 W k = Wt + t. Since we know E Wt = t, we see that E W s dw s = E Wt + t = t. Similarly we have and so W s dw s = lim π k= = lim π k= = W n W = W t, ) ) Wk+1 + W k Wk+1 W k W k+1 Wk ) E W s dw s = E Wt = t. 3. Assuming that both the integral and its variance exist, show that var fw s, s) dw s = E fw s, s) ds. Note: if the integral and its variance exist then it is legitimate to interchange the order of expectation and integration. If the integral exists then its expected value is zero so ) var fw s, s) dw s = E fw s, s) dw s. Set f k = fw k, t k ) and consider ) f k δw k = k= = f j f k δw j δw k k= j= fk δw k) + f j f k δw j δw k. k= k=1 j<k

3 We can eliminate the double sum using the tower law, E f j f k δw j δw k = E E tk fj f k δw j δw k = E f j f k δw j E tk δwk since t j < t k, which means f j and δw j are known at time t k, as is f k = fw k, t k ). As E tk δwk =, this term vanishes, leaving ) f k δw k E k= = E fk δw k) = E = E k= E tk f k δw k ) k= fk E t k δwk ) k= = E fk δt k = E fk δtk k= k= and as we refine the partition, π, we see lim π k= E fk δtk E fw s, s) ds. 4. Use the differential version of Itô s lemma to show that a) W s ds = t W t We have which integrates to show Rearranging gives s dw s = t s) dw s, d t W t ) = Wt dt + t dw t t W t = W s ds = t W t W s ds + s dw s = s dw s. t s) dw s. 3

4 b) W s dw s = 1 3 W 3 t W s ds, This time d Wt 3 ) = 3 W t dw t + 3 W t dt which integrates to give W 3 t = 3 W s dw s + 3 Dividing by 3 and rearranging gives W s dw s = 1 3 W 3 t W s ds. W s ds. c) E e awt a t/ = 1. Set fw, t) = expaw 1 a t) and f t = fw t, t). Itô s lemma gives df t = 1 a f t dt + a f t dw t + 1 a f t dt = a f t dw t. This integrates to give f t f = a exp aw t 1 a t ) = 1 + a Taking expectations gives E exp aw t 1 a t ) = 1 + a E = 1. f s dw s or, in full, exp aw s 1 a s ) dw s exp aw s 1 a s ) dw s 5. Define X t to be the area under a Brownian motion, X = and X t = W u du for t >. Show that X t is normally distributed with From Question 4a) we have E X t =, E X t = 1 3 t3. W s ds = t s) dw s. 4

5 As t s does not depend on W s, the integral is normally distributed with E W s ds = E t s) dw s = and var W s ds = var t s) dw s = = 1 3 t3. t s) ds Note that X t is continuously differentiable for t >, with Ẋt = W t recall W t is continuous in t). Now define Y t as the average area under a Brownian motion, { if t =, Y t = X t /t if t >. Show that Y t has EY t =, E Yt = t/3 and that Yt is continuous for all t. Is 3 Y t a Brownian motion? Give reasons for you answer. For any t >, Y t is normal because X t is, with and E Y t = E Xt /t = var Y t = var Xt /t = var X t /t = 1 3 t. For t > we have Y t = X t /t which is the ratio of two differentiable functions and as the denominator is never zero it follows that Y t is differentiable for t >, which implies continuous. Moreover, it means we can use l Hopital s rule to show lim Y W t t = lim t + t + 1 =, which shows that Y t is continuous at t =. The function 3 Y t is not only continuous but differentiable for t > and therefore it can t be a Brownian motion. The hard way to do this part of the question is to show that the increments over disjoint intervals are not independent.) 6. Show that if ds t S t = µ dt + σ dw t, 5

6 and f t = fs t ) = S 3 t, then Write the SDE for S t as df t f t = 3 µ + σ ) dt + 3 σ dw t. ds t = µ S t dt + σ S t dw t so Itô s lemma for a general F t = F S t, t) is df t = F t S t, t) + 1 σ S t F ) S S t, t) dt + F S S t, t) ds t. With fs) = S 3, so f S) = 3 S and f S) = 6 S, this becomes df t = 3 σ St 3 dt + 3 St ds t = 3 σ f t dt + St ) ) µ St dt + σ S t dw t = 3 f t µ + σ ) ) dt + σ dw t which we could write as df t f t = 3 µ + σ ) dt + 3 σ dw t. 7. Find solutions of the Black-Scholes terminal value problem when + r y) S r V =, S >, t < T, S S V S, T ) = fs), S >, a) fs) = C, where C is a constant; In this case it makes sense to look for solutions which don t depend on S, so / S = and V/ S = identically. Then we have V = V t) and the partial differential equation reduces to the ODE dv dt r V =, V T ) = C which has the solution V = C e rt t). b) fs) = S α, where α is a constant. Hint: you don t need the Feynman-Kǎc formula to do either of these. Look for simple functional forms of the solution. 6

7 In this case, let s look for separable solutions We find that t = Sα ft), V S, t) = S α ft), ft ) = 1. S S = α Sα ft), S V S = αα 1) Sα ft) and the Black-Scholes equation becomes ) S α ft) + 1 σ αα 1) ft) + r y) α ft) r ft) =. As this has to hold for all S > we have an ordinary differential equation for ft) which we can write as where The solution is and ft) = λ ft), ft ) = 1, λ = r r y) α 1 σ αα 1). λt t) ft) = e V S, t) = e λt t) S α. 7

8 Optional questions 8. The Black-Scholes equation from a binomial method. One step of the Cox, Ross & Rubinstein binomial method can be written as V S, t) = e rδt q V u + 1 q) V d) where V u = V S u, t + δt ), V d = V S d, t + δt ), δt S u = S e σ δt, S d = S e σ δt, q = erδt e σ e σ δt e, σ δt σ > is the volatility, r is the risk-free rate and δt > is the length of the time-step. Supposing this is true for all S > and that V S, t) may be expanded in a Taylor series in both S and t, show that as δt q = 1 + r 1 σ δt + Oδt), σ V u = V + δt σ S S + δt t + 1 σ S V d = V δt σ S S + δt S + S V )) S + Oδt 3/ ), t + 1 σ S S + S V )) S + Oδt 3/ ), where V and all its partial derivatives are evaluated at S, t). Hence show that in the limit δt the option price satisfies the zero dividend-yield) Black-Scholes equation As δt we can expand S + r S S S r V =. e σ δt = 1 + σ δt + 1 σ δt + Oδt 3/ ), from which it follows that e σ δt = 1 σ δt + 1 σ δt + Oδt 3/ ), e rδt = 1 + rδt + Oδt ), q = σ δt + r 1 σ )δt + Oδt 3/ ) σ δt + Oδt 3/ ) = 1 + r 1 σ δt + Oδt), σ 8

9 and that V S e σ δt, t + δt) = V S, t) + δt t S, t) + σ δt + 1 σ δt ) S S, t) S + 1 σ δt S V S S, t) + Oδt3/ ) = V + δt σ S S + δt V S e σ δt, t + δt) Therefore t + 1 σ S )) S + S V S + Oδt 3/ ) = V δt σ S S + δt t + 1 σ S )) S + S V S + Oδt 3/ ). q V S e σ δt, t + δt) + 1 q) V S e σ δt, t + δt) = V + δt S + 1 σ S S = V + δt S + r S S ) + r 1 σ ) δt S S + Oδt3/ ) ) + Oδt 3/ ) and hence ) e rδt q V S e σ δt, t + δt) + 1 q) V S e σ δt, t + δt) = 1 rδt) V + δt S + r S S = V + δt S + r S S r V )) + Oδt 3/ ) ) + Oδt 3/ ). Substituting this back into the binomial equation gives V = V + δt S + r S ) S r V + Oδt 3/ ) and cancelling V and dividing by δt then gives S + r S ) S r V + Oδt 1/ ) =. Taking the limit δt gives the Black-Scholes equation, S + r S S r V =. 9

10 9. The variation of a function, or stochastic process, over, t, is f t = lim π k= f k+1 f k. If f t is finite on, t we say f has bounded variation on, t. Show that: a) if f is C 1, t then f t = f t) dt < ; If f is C 1 on, t then the intermediate value theorem asserts that f k+1 f k = f ξ k )t k+1 t k ) for some ξ k t k, t k+1. Since δt k = t k+1 t k ) >, f t = lim π k= f ξ k ) δt k = f s) ds. If f s) is continuous on, t then so too is f s) and as, t is a closed and bounded set, f s) must be bounded on, t. Therefore the integral must be finite and hence f has bounded variation on, t. b) if f is a continuous function with f t < then its quadratic variation is zero, f t = ; We have ) f t = lim fk+1 f k and we also have π k= ) ) fk+1 f k max fj+1 f j fk+1 f k ). k= j If f has bounded variation then the sum on the right-hand side is finite as π. If f is continuous then the maximum value of f j+1 f j as π. Thus the product of the two vanishes as π. c) Brownian motion does not have bounded variation; Brownian motion is continuous in t and has non-zero quadratic variation. Therefore it must have unbounded variation by the previous part of the question. d) the arc length of the graph of a Brownian motion is infinite for any t >. Hint: if y = X t has an arc length s then ds = dy + dt dy = dy. k= 1

11 Assume that we can define arc length for a Brownian motion. Then it is given by s = where if y = X t then ds is given by Now observe that ds ds = dt + dy, ds = dt + dy. ds = dt + dy dy = dy. If the integral of dy existed it would equal the variation of W t, but we know that W t doesn t have bounded variation. This implies that the length of the graph of a Brownian motion is infinite for any t > ). 11

MATH 56A SPRING 2008 STOCHASTIC PROCESSES 197

MATH 56A SPRING 2008 STOCHASTIC PROCESSES 197 MATH 56A SPRING 8 STOCHASTIC PROCESSES 197 9.3. Itô s formula. First I stated the theorem. Then I did a simple example to make sure we understand what it says. Then I proved it. The key point is Lévy s

More information

Stochastic Calculus Made Easy

Stochastic Calculus Made Easy Stochastic Calculus Made Easy Most of us know how standard Calculus works. We know how to differentiate, how to integrate etc. But stochastic calculus is a totally different beast to tackle; we are trying

More information

The concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.

The concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt. The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes

More information

Lecture 4: Ito s Stochastic Calculus and SDE. Seung Yeal Ha Dept of Mathematical Sciences Seoul National University

Lecture 4: Ito s Stochastic Calculus and SDE. Seung Yeal Ha Dept of Mathematical Sciences Seoul National University Lecture 4: Ito s Stochastic Calculus and SDE Seung Yeal Ha Dept of Mathematical Sciences Seoul National University 1 Preliminaries What is Calculus? Integral, Differentiation. Differentiation 2 Integral

More information

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539 Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory

More information

Stochastic Calculus. Kevin Sinclair. August 2, 2016

Stochastic Calculus. Kevin Sinclair. August 2, 2016 Stochastic Calculus Kevin Sinclair August, 16 1 Background Suppose we have a Brownian motion W. This is a process, and the value of W at a particular time T (which we write W T ) is a normally distributed

More information

Itô s formula. Samy Tindel. Purdue University. Probability Theory 2 - MA 539

Itô s formula. Samy Tindel. Purdue University. Probability Theory 2 - MA 539 Itô s formula Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Itô s formula Probability Theory

More information

Stochastic Calculus for Finance II - some Solutions to Chapter VII

Stochastic Calculus for Finance II - some Solutions to Chapter VII Stochastic Calculus for Finance II - some Solutions to Chapter VII Matthias hul Last Update: June 9, 25 Exercise 7 Black-Scholes-Merton Equation for the up-and-out Call) i) We have ii) We first compute

More information

n E(X t T n = lim X s Tn = X s

n E(X t T n = lim X s Tn = X s Stochastic Calculus Example sheet - Lent 15 Michael Tehranchi Problem 1. Let X be a local martingale. Prove that X is a uniformly integrable martingale if and only X is of class D. Solution 1. If If direction:

More information

5 December 2016 MAA136 Researcher presentation. Anatoliy Malyarenko. Topics for Bachelor and Master Theses. Anatoliy Malyarenko

5 December 2016 MAA136 Researcher presentation. Anatoliy Malyarenko. Topics for Bachelor and Master Theses. Anatoliy Malyarenko 5 December 216 MAA136 Researcher presentation 1 schemes The main problem of financial engineering: calculate E [f(x t (x))], where {X t (x): t T } is the solution of the system of X t (x) = x + Ṽ (X s

More information

Discretization of SDEs: Euler Methods and Beyond

Discretization of SDEs: Euler Methods and Beyond Discretization of SDEs: Euler Methods and Beyond 09-26-2006 / PRisMa 2006 Workshop Outline Introduction 1 Introduction Motivation Stochastic Differential Equations 2 The Time Discretization of SDEs Monte-Carlo

More information

FE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology

FE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology FE610 Stochastic Calculus for Financial Engineers Lecture 3. Calculaus in Deterministic and Stochastic Environments Steve Yang Stevens Institute of Technology 01/31/2012 Outline 1 Modeling Random Behavior

More information

Stochastic Integration and Stochastic Differential Equations: a gentle introduction

Stochastic Integration and Stochastic Differential Equations: a gentle introduction Stochastic Integration and Stochastic Differential Equations: a gentle introduction Oleg Makhnin New Mexico Tech Dept. of Mathematics October 26, 27 Intro: why Stochastic? Brownian Motion/ Wiener process

More information

Stochastic differential equation models in biology Susanne Ditlevsen

Stochastic differential equation models in biology Susanne Ditlevsen Stochastic differential equation models in biology Susanne Ditlevsen Introduction This chapter is concerned with continuous time processes, which are often modeled as a system of ordinary differential

More information

Stochastic Processes and Advanced Mathematical Finance

Stochastic Processes and Advanced Mathematical Finance Steven R. Dunbar Department of Mathematics 23 Avery Hall University of Nebraska-Lincoln Lincoln, NE 68588-13 http://www.math.unl.edu Voice: 42-472-3731 Fax: 42-472-8466 Stochastic Processes and Advanced

More information

Laplace Transforms Chapter 3

Laplace Transforms Chapter 3 Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important

More information

Accurate approximation of stochastic differential equations

Accurate approximation of stochastic differential equations Accurate approximation of stochastic differential equations Simon J.A. Malham and Anke Wiese (Heriot Watt University, Edinburgh) Birmingham: 6th February 29 Stochastic differential equations dy t = V (y

More information

Brownian Motion. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Brownian Motion

Brownian Motion. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Brownian Motion Brownian Motion An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Background We have already seen that the limiting behavior of a discrete random walk yields a derivation of

More information

Introduction to numerical simulations for Stochastic ODEs

Introduction to numerical simulations for Stochastic ODEs Introduction to numerical simulations for Stochastic ODEs Xingye Kan Illinois Institute of Technology Department of Applied Mathematics Chicago, IL 60616 August 9, 2010 Outline 1 Preliminaries 2 Numerical

More information

1. Stochastic Processes and filtrations

1. Stochastic Processes and filtrations 1. Stochastic Processes and 1. Stoch. pr., A stochastic process (X t ) t T is a collection of random variables on (Ω, F) with values in a measurable space (S, S), i.e., for all t, In our case X t : Ω S

More information

Stochastic Differential Equations.

Stochastic Differential Equations. Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)

More information

MFE6516 Stochastic Calculus for Finance

MFE6516 Stochastic Calculus for Finance MFE6516 Stochastic Calculus for Finance William C. H. Leon Nanyang Business School December 11, 2017 1 / 51 William C. H. Leon MFE6516 Stochastic Calculus for Finance 1 Symmetric Random Walks Scaled Symmetric

More information

An Uncertain Control Model with Application to. Production-Inventory System

An Uncertain Control Model with Application to. Production-Inventory System An Uncertain Control Model with Application to Production-Inventory System Kai Yao 1, Zhongfeng Qin 2 1 Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China 2 School of Economics

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

Numerical Integration of SDEs: A Short Tutorial

Numerical Integration of SDEs: A Short Tutorial Numerical Integration of SDEs: A Short Tutorial Thomas Schaffter January 19, 010 1 Introduction 1.1 Itô and Stratonovich SDEs 1-dimensional stochastic differentiable equation (SDE) is given by [6, 7] dx

More information

MA50251: Applied SDEs

MA50251: Applied SDEs MA5251: Applied SDEs Alexander Cox February 12, 218 Preliminaries These notes constitute the core theoretical content of the unit MA5251, which is a course on Applied SDEs. Roughly speaking, the aim of

More information

Stochastic Calculus February 11, / 33

Stochastic Calculus February 11, / 33 Martingale Transform M n martingale with respect to F n, n =, 1, 2,... σ n F n (σ M) n = n 1 i= σ i(m i+1 M i ) is a Martingale E[(σ M) n F n 1 ] n 1 = E[ σ i (M i+1 M i ) F n 1 ] i= n 2 = σ i (M i+1 M

More information

Stochastic Calculus and Black-Scholes Theory MTH772P Exercises Sheet 1

Stochastic Calculus and Black-Scholes Theory MTH772P Exercises Sheet 1 Stochastic Calculus and Black-Scholes Theory MTH772P Exercises Sheet. For ξ, ξ 2, i.i.d. with P(ξ i = ± = /2 define the discrete-time random walk W =, W n = ξ +... + ξ n. (i Formulate and prove the property

More information

Topic 5.1: Line Element and Scalar Line Integrals

Topic 5.1: Line Element and Scalar Line Integrals Math 275 Notes Topic 5.1: Line Element and Scalar Line Integrals Textbook Section: 16.2 More Details on Line Elements (vector dr, and scalar ds): http://www.math.oregonstate.edu/bridgebook/book/math/drvec

More information

Geometric projection of stochastic differential equations

Geometric projection of stochastic differential equations Geometric projection of stochastic differential equations John Armstrong (King s College London) Damiano Brigo (Imperial) August 9, 2018 Idea: Projection Idea: Projection Projection gives a method of systematically

More information

I forgot to mention last time: in the Ito formula for two standard processes, putting

I forgot to mention last time: in the Ito formula for two standard processes, putting I forgot to mention last time: in the Ito formula for two standard processes, putting dx t = a t dt + b t db t dy t = α t dt + β t db t, and taking f(x, y = xy, one has f x = y, f y = x, and f xx = f yy

More information

Gaussian, Markov and stationary processes

Gaussian, Markov and stationary processes Gaussian, Markov and stationary processes Gonzalo Mateos Dept. of ECE and Goergen Institute for Data Science University of Rochester gmateosb@ece.rochester.edu http://www.ece.rochester.edu/~gmateosb/ November

More information

BIRKBECK (University of London) MATHEMATICAL AND NUMERICAL METHODS. Thursday, 29 May 2008, 10.00am-13.15pm (includes 15 minutes reading time).

BIRKBECK (University of London) MATHEMATICAL AND NUMERICAL METHODS. Thursday, 29 May 2008, 10.00am-13.15pm (includes 15 minutes reading time). BIRKBECK (University of London) MSc EXAMINATION FOR INTERNAL STUDENTS MSc FINANCIAL ENGINEERING School of Economics, Mathematics, and Statistics MATHEMATICAL AND NUMERICAL METHODS ANSWERS EMMS011P 200804301248

More information

MA8109 Stochastic Processes in Systems Theory Autumn 2013

MA8109 Stochastic Processes in Systems Theory Autumn 2013 Norwegian University of Science and Technology Department of Mathematical Sciences MA819 Stochastic Processes in Systems Theory Autumn 213 1 MA819 Exam 23, problem 3b This is a linear equation of the form

More information

Bernardo D Auria Stochastic Processes /12. Notes. March 29 th, 2012

Bernardo D Auria Stochastic Processes /12. Notes. March 29 th, 2012 1 Stochastic Calculus Notes March 9 th, 1 In 19, Bachelier proposed for the Paris stock exchange a model for the fluctuations affecting the price X(t) of an asset that was given by the Brownian motion.

More information

Stochastic Modelling in Climate Science

Stochastic Modelling in Climate Science Stochastic Modelling in Climate Science David Kelly Mathematics Department UNC Chapel Hill dtbkelly@gmail.com November 16, 2013 David Kelly (UNC) Stochastic Climate November 16, 2013 1 / 36 Why use stochastic

More information

Economics 2010c: Lectures 9-10 Bellman Equation in Continuous Time

Economics 2010c: Lectures 9-10 Bellman Equation in Continuous Time Economics 2010c: Lectures 9-10 Bellman Equation in Continuous Time David Laibson 9/30/2014 Outline Lectures 9-10: 9.1 Continuous-time Bellman Equation 9.2 Application: Merton s Problem 9.3 Application:

More information

Clases 11-12: Integración estocástica.

Clases 11-12: Integración estocástica. Clases 11-12: Integración estocástica. Fórmula de Itô * 3 de octubre de 217 Índice 1. Introduction to Stochastic integrals 1 2. Stochastic integration 2 3. Simulation of stochastic integrals: Euler scheme

More information

Kolmogorov Equations and Markov Processes

Kolmogorov Equations and Markov Processes Kolmogorov Equations and Markov Processes May 3, 013 1 Transition measures and functions Consider a stochastic process {X(t)} t 0 whose state space is a product of intervals contained in R n. We define

More information

UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE

UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (2010), 275 284 UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE Iuliana Carmen Bărbăcioru Abstract.

More information

Malliavin Calculus in Finance

Malliavin Calculus in Finance Malliavin Calculus in Finance Peter K. Friz 1 Greeks and the logarithmic derivative trick Model an underlying assent by a Markov process with values in R m with dynamics described by the SDE dx t = b(x

More information

Lecture 7 - Separable Equations

Lecture 7 - Separable Equations Lecture 7 - Separable Equations Separable equations is a very special type of differential equations where you can separate the terms involving only y on one side of the equation and terms involving only

More information

(B(t i+1 ) B(t i )) 2

(B(t i+1 ) B(t i )) 2 ltcc5.tex Week 5 29 October 213 Ch. V. ITÔ (STOCHASTIC) CALCULUS. WEAK CONVERGENCE. 1. Quadratic Variation. A partition π n of [, t] is a finite set of points t ni such that = t n < t n1

More information

Verona Course April Lecture 1. Review of probability

Verona Course April Lecture 1. Review of probability Verona Course April 215. Lecture 1. Review of probability Viorel Barbu Al.I. Cuza University of Iaşi and the Romanian Academy A probability space is a triple (Ω, F, P) where Ω is an abstract set, F is

More information

A numerical method for solving uncertain differential equations

A numerical method for solving uncertain differential equations Journal of Intelligent & Fuzzy Systems 25 (213 825 832 DOI:1.3233/IFS-12688 IOS Press 825 A numerical method for solving uncertain differential equations Kai Yao a and Xiaowei Chen b, a Department of Mathematical

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

Nested Uncertain Differential Equations and Its Application to Multi-factor Term Structure Model

Nested Uncertain Differential Equations and Its Application to Multi-factor Term Structure Model Nested Uncertain Differential Equations and Its Application to Multi-factor Term Structure Model Xiaowei Chen International Business School, Nankai University, Tianjin 371, China School of Finance, Nankai

More information

Introduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0

Introduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0 Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods

More information

Exercises. T 2T. e ita φ(t)dt.

Exercises. T 2T. e ita φ(t)dt. Exercises. Set #. Construct an example of a sequence of probability measures P n on R which converge weakly to a probability measure P but so that the first moments m,n = xdp n do not converge to m = xdp.

More information

Example 4.1 Let X be a random variable and f(t) a given function of time. Then. Y (t) = f(t)x. Y (t) = X sin(ωt + δ)

Example 4.1 Let X be a random variable and f(t) a given function of time. Then. Y (t) = f(t)x. Y (t) = X sin(ωt + δ) Chapter 4 Stochastic Processes 4. Definition In the previous chapter we studied random variables as functions on a sample space X(ω), ω Ω, without regard to how these might depend on parameters. We now

More information

Solution for Problem 7.1. We argue by contradiction. If the limit were not infinite, then since τ M (ω) is nondecreasing we would have

Solution for Problem 7.1. We argue by contradiction. If the limit were not infinite, then since τ M (ω) is nondecreasing we would have 362 Problem Hints and Solutions sup g n (ω, t) g(ω, t) sup g(ω, s) g(ω, t) µ n (ω). t T s,t: s t 1/n By the uniform continuity of t g(ω, t) on [, T], one has for each ω that µ n (ω) as n. Two applications

More information

STATISTICS 385: STOCHASTIC CALCULUS HOMEWORK ASSIGNMENT 4 DUE NOVEMBER 23, = (2n 1)(2n 3) 3 1.

STATISTICS 385: STOCHASTIC CALCULUS HOMEWORK ASSIGNMENT 4 DUE NOVEMBER 23, = (2n 1)(2n 3) 3 1. STATISTICS 385: STOCHASTIC CALCULUS HOMEWORK ASSIGNMENT 4 DUE NOVEMBER 23, 26 Problem Normal Moments (A) Use the Itô formula and Brownian scaling to check that the even moments of the normal distribution

More information

Brownian Motion and Poisson Process

Brownian Motion and Poisson Process and Poisson Process She: What is white noise? He: It is the best model of a totally unpredictable process. She: Are you implying, I am white noise? He: No, it does not exist. Dialogue of an unknown couple.

More information

Stochastic Calculus (Lecture #3)

Stochastic Calculus (Lecture #3) Stochastic Calculus (Lecture #3) Siegfried Hörmann Université libre de Bruxelles (ULB) Spring 2014 Outline of the course 1. Stochastic processes in continuous time. 2. Brownian motion. 3. Itô integral:

More information

Laplace Transforms. Chapter 3. Pierre Simon Laplace Born: 23 March 1749 in Beaumont-en-Auge, Normandy, France Died: 5 March 1827 in Paris, France

Laplace Transforms. Chapter 3. Pierre Simon Laplace Born: 23 March 1749 in Beaumont-en-Auge, Normandy, France Died: 5 March 1827 in Paris, France Pierre Simon Laplace Born: 23 March 1749 in Beaumont-en-Auge, Normandy, France Died: 5 March 1827 in Paris, France Laplace Transforms Dr. M. A. A. Shoukat Choudhury 1 Laplace Transforms Important analytical

More information

Controlled Diffusions and Hamilton-Jacobi Bellman Equations

Controlled Diffusions and Hamilton-Jacobi Bellman Equations Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter

More information

A Short Introduction to Diffusion Processes and Ito Calculus

A Short Introduction to Diffusion Processes and Ito Calculus A Short Introduction to Diffusion Processes and Ito Calculus Cédric Archambeau University College, London Center for Computational Statistics and Machine Learning c.archambeau@cs.ucl.ac.uk January 24,

More information

Infinite-dimensional methods for path-dependent equations

Infinite-dimensional methods for path-dependent equations Infinite-dimensional methods for path-dependent equations (Università di Pisa) 7th General AMaMeF and Swissquote Conference EPFL, Lausanne, 8 September 215 Based on Flandoli F., Zanco G. - An infinite-dimensional

More information

GAUSSIAN PROCESSES; KOLMOGOROV-CHENTSOV THEOREM

GAUSSIAN PROCESSES; KOLMOGOROV-CHENTSOV THEOREM GAUSSIAN PROCESSES; KOLMOGOROV-CHENTSOV THEOREM STEVEN P. LALLEY 1. GAUSSIAN PROCESSES: DEFINITIONS AND EXAMPLES Definition 1.1. A standard (one-dimensional) Wiener process (also called Brownian motion)

More information

If g is also continuous and strictly increasing on J, we may apply the strictly increasing inverse function g 1 to this inequality to get

If g is also continuous and strictly increasing on J, we may apply the strictly increasing inverse function g 1 to this inequality to get 18:2 1/24/2 TOPIC. Inequalities; measures of spread. This lecture explores the implications of Jensen s inequality for g-means in general, and for harmonic, geometric, arithmetic, and related means in

More information

Optimal Trade Execution with Instantaneous Price Impact and Stochastic Resilience

Optimal Trade Execution with Instantaneous Price Impact and Stochastic Resilience Optimal Trade Execution with Instantaneous Price Impact and Stochastic Resilience Ulrich Horst 1 Humboldt-Universität zu Berlin Department of Mathematics and School of Business and Economics Vienna, Nov.

More information

The Cameron-Martin-Girsanov (CMG) Theorem

The Cameron-Martin-Girsanov (CMG) Theorem The Cameron-Martin-Girsanov (CMG) Theorem There are many versions of the CMG Theorem. In some sense, there are many CMG Theorems. The first version appeared in ] in 944. Here we present a standard version,

More information

Definition of Derivative

Definition of Derivative Definition of Derivative The derivative of the function f with respect to the variable x is the function ( ) fʹ x whose value at xis ( x) fʹ = lim provided the limit exists. h 0 ( + ) ( ) f x h f x h Slide

More information

Central European University

Central European University Central European University Department of Applied Mathematics thesis Pathwise approximation of the Feynman-Kac formula Author: Denis Nagovitcyn Supervisor: Tamás Szabados May 18, 15 Contents Introduction

More information

Math 162: Calculus IIA

Math 162: Calculus IIA Math 62: Calculus IIA Final Exam ANSWERS December 9, 26 Part A. (5 points) Evaluate the integral x 4 x 2 dx Substitute x 2 cos θ: x 8 cos dx θ ( 2 sin θ) dθ 4 x 2 2 sin θ 8 cos θ dθ 8 cos 2 θ cos θ dθ

More information

The method to find a basis of Euler-Cauchy equation by transforms

The method to find a basis of Euler-Cauchy equation by transforms Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 12, Number 5 (2016), pp. 4159 4165 Research India Publications http://www.ripublication.com/gjpam.htm The method to find a basis of

More information

Some Terminology and Concepts that We will Use, But Not Emphasize (Section 6.2)

Some Terminology and Concepts that We will Use, But Not Emphasize (Section 6.2) Some Terminology and Concepts that We will Use, But Not Emphasize (Section 6.2) Statistical analysis is based on probability theory. The fundamental object in probability theory is a probability space,

More information

Bernardo D Auria Stochastic Processes /10. Notes. Abril 13 th, 2010

Bernardo D Auria Stochastic Processes /10. Notes. Abril 13 th, 2010 1 Stochastic Calculus Notes Abril 13 th, 1 As we have seen in previous lessons, the stochastic integral with respect to the Brownian motion shows a behavior different from the classical Riemann-Stieltjes

More information

MARTINGALES: VARIANCE. EY = 0 Var (Y ) = σ 2

MARTINGALES: VARIANCE. EY = 0 Var (Y ) = σ 2 MARTINGALES: VARIANCE SIMPLEST CASE: X t = Y 1 +... + Y t IID SUM V t = X 2 t σ 2 t IS A MG: EY = 0 Var (Y ) = σ 2 V t+1 = (X t + Y t+1 ) 2 σ 2 (t + 1) = V t + 2X t Y t+1 + Y 2 t+1 σ 2 E(V t+1 F t ) =

More information

1. Stochastic Process

1. Stochastic Process HETERGENEITY IN QUANTITATIVE MACROECONOMICS @ TSE OCTOBER 17, 216 STOCHASTIC CALCULUS BASICS SANG YOON (TIM) LEE Very simple notes (need to add references). It is NOT meant to be a substitute for a real

More information

AMS 216 Stochastic Differential Equations Lecture 03 Copyright by Hongyun Wang, UCSC

AMS 216 Stochastic Differential Equations Lecture 03 Copyright by Hongyun Wang, UCSC Lecture 03 Copyright by Hongyun Wang, UCSC Review of probability theory (Continued) We show that Sum of independent normal random variable normal random variable Characteristic function (CF) of a random

More information

Bayesian Inference and the Symbolic Dynamics of Deterministic Chaos. Christopher C. Strelioff 1,2 Dr. James P. Crutchfield 2

Bayesian Inference and the Symbolic Dynamics of Deterministic Chaos. Christopher C. Strelioff 1,2 Dr. James P. Crutchfield 2 How Random Bayesian Inference and the Symbolic Dynamics of Deterministic Chaos Christopher C. Strelioff 1,2 Dr. James P. Crutchfield 2 1 Center for Complex Systems Research and Department of Physics University

More information

From Random Variables to Random Processes. From Random Variables to Random Processes

From Random Variables to Random Processes. From Random Variables to Random Processes Random Processes In probability theory we study spaces (Ω, F, P) where Ω is the space, F are all the sets to which we can measure its probability and P is the probability. Example: Toss a die twice. Ω

More information

MULTIDIMENSIONAL WICK-ITÔ FORMULA FOR GAUSSIAN PROCESSES

MULTIDIMENSIONAL WICK-ITÔ FORMULA FOR GAUSSIAN PROCESSES MULTIDIMENSIONAL WICK-ITÔ FORMULA FOR GAUSSIAN PROCESSES D. NUALART Department of Mathematics, University of Kansas Lawrence, KS 6645, USA E-mail: nualart@math.ku.edu S. ORTIZ-LATORRE Departament de Probabilitat,

More information

Some Tools From Stochastic Analysis

Some Tools From Stochastic Analysis W H I T E Some Tools From Stochastic Analysis J. Potthoff Lehrstuhl für Mathematik V Universität Mannheim email: potthoff@math.uni-mannheim.de url: http://ls5.math.uni-mannheim.de To close the file, click

More information

Mat104 Fall 2002, Improper Integrals From Old Exams

Mat104 Fall 2002, Improper Integrals From Old Exams Mat4 Fall 22, Improper Integrals From Old Eams For the following integrals, state whether they are convergent or divergent, and give your reasons. () (2) (3) (4) (5) converges. Break it up as 3 + 2 3 +

More information

High order Cubature on Wiener space applied to derivative pricing and sensitivity estimations

High order Cubature on Wiener space applied to derivative pricing and sensitivity estimations 1/25 High order applied to derivative pricing and sensitivity estimations University of Oxford 28 June 2010 2/25 Outline 1 2 3 4 3/25 Assume that a set of underlying processes satisfies the Stratonovich

More information

Math 212-Lecture 8. The chain rule with one independent variable

Math 212-Lecture 8. The chain rule with one independent variable Math 212-Lecture 8 137: The multivariable chain rule The chain rule with one independent variable w = f(x, y) If the particle is moving along a curve x = x(t), y = y(t), then the values that the particle

More information

Stochastic Differential Equations

Stochastic Differential Equations CHAPTER 1 Stochastic Differential Equations Consider a stochastic process X t satisfying dx t = bt, X t,w t dt + σt, X t,w t dw t. 1.1 Question. 1 Can we obtain the existence and uniqueness theorem for

More information

Multivariate Generalized Ornstein-Uhlenbeck Processes

Multivariate Generalized Ornstein-Uhlenbeck Processes Multivariate Generalized Ornstein-Uhlenbeck Processes Anita Behme TU München Alexander Lindner TU Braunschweig 7th International Conference on Lévy Processes: Theory and Applications Wroclaw, July 15 19,

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Introduction to Malliavin calculus and its applications Lecture 3: Clark-Ocone formula David Nualart Department of Mathematics Kansas University University of Wyoming Summer School 214 David Nualart

More information

Numerical methods for solving stochastic differential equations

Numerical methods for solving stochastic differential equations Mathematical Communications 4(1999), 251-256 251 Numerical methods for solving stochastic differential equations Rózsa Horváth Bokor Abstract. This paper provides an introduction to stochastic calculus

More information

Lecture 4: Introduction to stochastic processes and stochastic calculus

Lecture 4: Introduction to stochastic processes and stochastic calculus Lecture 4: Introduction to stochastic processes and stochastic calculus Cédric Archambeau Centre for Computational Statistics and Machine Learning Department of Computer Science University College London

More information

EULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS

EULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS Qiao, H. Osaka J. Math. 51 (14), 47 66 EULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS HUIJIE QIAO (Received May 6, 11, revised May 1, 1) Abstract In this paper we show

More information

Stochastic Integration.

Stochastic Integration. Chapter Stochastic Integration..1 Brownian Motion as a Martingale P is the Wiener measure on (Ω, B) where Ω = C, T B is the Borel σ-field on Ω. In addition we denote by B t the σ-field generated by x(s)

More information

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3 Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................

More information

Week 9 Generators, duality, change of measure

Week 9 Generators, duality, change of measure Week 9 Generators, duality, change of measure Jonathan Goodman November 18, 013 1 Generators This section describes a common abstract way to describe many of the differential equations related to Markov

More information

Sample Final Questions: Solutions Math 21B, Winter y ( y 1)(1 + y)) = A y + B

Sample Final Questions: Solutions Math 21B, Winter y ( y 1)(1 + y)) = A y + B Sample Final Questions: Solutions Math 2B, Winter 23. Evaluate the following integrals: tan a) y y dy; b) x dx; c) 3 x 2 + x dx. a) We use partial fractions: y y 3 = y y ) + y)) = A y + B y + C y +. Putting

More information

Tyler Hofmeister. University of Calgary Mathematical and Computational Finance Laboratory

Tyler Hofmeister. University of Calgary Mathematical and Computational Finance Laboratory JUMP PROCESSES GENERALIZING STOCHASTIC INTEGRALS WITH JUMPS Tyler Hofmeister University of Calgary Mathematical and Computational Finance Laboratory Overview 1. General Method 2. Poisson Processes 3. Diffusion

More information

Mean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations

Mean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations Mean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations Ram Sharan Adhikari Assistant Professor Of Mathematics Rogers State University Mathematical

More information

On pathwise stochastic integration

On pathwise stochastic integration On pathwise stochastic integration Rafa l Marcin Lochowski Afican Institute for Mathematical Sciences, Warsaw School of Economics UWC seminar Rafa l Marcin Lochowski (AIMS, WSE) On pathwise stochastic

More information

1 Solution to Problem 2.1

1 Solution to Problem 2.1 Solution to Problem 2. I incorrectly worked this exercise instead of 2.2, so I decided to include the solution anyway. a) We have X Y /3, which is a - function. It maps the interval, ) where X lives) onto

More information

Risk-Minimality and Orthogonality of Martingales

Risk-Minimality and Orthogonality of Martingales Risk-Minimality and Orthogonality of Martingales Martin Schweizer Universität Bonn Institut für Angewandte Mathematik Wegelerstraße 6 D 53 Bonn 1 (Stochastics and Stochastics Reports 3 (199, 123 131 2

More information

Time Response Analysis (Part II)

Time Response Analysis (Part II) Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary

More information

Homogeneous Equations with Constant Coefficients

Homogeneous Equations with Constant Coefficients Homogeneous Equations with Constant Coefficients MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 General Second Order ODE Second order ODEs have the form

More information

Moment Properties of Distributions Used in Stochastic Financial Models

Moment Properties of Distributions Used in Stochastic Financial Models Moment Properties of Distributions Used in Stochastic Financial Models Jordan Stoyanov Newcastle University (UK) & University of Ljubljana (Slovenia) e-mail: stoyanovj@gmail.com ETH Zürich, Department

More information

Practice Midterm Solutions

Practice Midterm Solutions Practice Midterm Solutions Math 4B: Ordinary Differential Equations Winter 20 University of California, Santa Barbara TA: Victoria Kala DO NOT LOOK AT THESE SOLUTIONS UNTIL YOU HAVE ATTEMPTED EVERY PROBLEM

More information

Variance Reduction Techniques for Monte Carlo Simulations with Stochastic Volatility Models

Variance Reduction Techniques for Monte Carlo Simulations with Stochastic Volatility Models Variance Reduction Techniques for Monte Carlo Simulations with Stochastic Volatility Models Jean-Pierre Fouque North Carolina State University SAMSI November 3, 5 1 References: Variance Reduction for Monte

More information

Topic 2-2: Derivatives of Vector Functions. Textbook: Section 13.2, 13.4

Topic 2-2: Derivatives of Vector Functions. Textbook: Section 13.2, 13.4 Topic 2-2: Derivatives of Vector Functions Textbook: Section 13.2, 13.4 Warm-Up: Parametrization of Circles Each of the following vector functions describe the position of an object traveling around the

More information

Linear Filtering of general Gaussian processes

Linear Filtering of general Gaussian processes Linear Filtering of general Gaussian processes Vít Kubelka Advisor: prof. RNDr. Bohdan Maslowski, DrSc. Robust 2018 Department of Probability and Statistics Faculty of Mathematics and Physics Charles University

More information