Exercise Sheet 8 - Solutions
|
|
- Theodore Johnson
- 5 years ago
- Views:
Transcription
1 Exercise Sheet 8 - Solutions Alessandro Gnoatto June 16, Exercise 1 Let X t = x + µ t + σ W t denote a scaled Brownian motion with constant diffusion coefficient σ and constant drift µ. 1. Prove that Y t = e Xt where Y (t = Y (0 + t 0 dy with Y (0 = ex and dy t = (µ + 1 σ Y t dt + σy t dw t. We shall assume that the SDE for dy has a unique 1 solution in [0, t]. Since the solution to the SDE for dy is unique, we will prove the claim by obtaining d(e Xt and verifying that it is the same SDE as dy. Let us then start with dx t. By definition of Ito processes, dx t = µdt + σdw t and therefore, using Ito s lemma to obtain df(x t with f(z = e z. We have f (z = f (z = e z and d(e Xt =e Xt dx t + 1 σ e Xt dt d(e Xt =(µ + σ ext dt + σe Xt dw t. Note that f(x 0 = e X0 = e x = Y (0. Hence f(x t = e Xt and Y t are both solutions to the same SDE, i.e., Y t = e Xt. Another way is to take d(log Y t and show that its solution equals X t. We have g(y = log y, g (y = 1 y, and g (y = 1 y. Thus, hence also i.e. d(log Y t =(µ + σ Y t dt Y t Y t Y =µdt + σdw t, log Y t = log Y (0 + t σ dt + σ Y t Y t dw t t 0 µds + t log Y t = x + µt + σw t = X t. 0 σdw s 1 For more details about the necessary conditions for uniqueness of SDE, see Section 3, Chapter 5 of the lecture notes Lawrence Evans, An Introduction to Stochastic Differential Equations, ver 1., UC Berkeley. 1
2 . An example Java program of the solution can be found in 4 package de. math. lmu. exercise08 ; import java. text. DecimalFormat ; import net. finmath. montecarlo. BrownianMotion ; 6 import net. finmath. montecarlo. RandomVariable ; import net. finmath. stochastic. RandomVariableInterface ; 8 import net. finmath. time. TimeDiscretization ; 10 public class LogNormalTest { 1 s t a t i c f i n a l DecimalFormat formattersci6 = new DecimalFormat (" E00 ; E00 "; 14 public s t a t i c void main ( String [] args { int seed = 134; 16 int numberofpaths = 1000; //double dt = 0.01; 18 double sigma = 0.1; double mu = 0.1; 0 double initialvalueybar = Math. exp (1.0 ; double initialvaluex = 1.0; //define time grid 4 int numberoftimes = 00; double t0 =0.0; 6 double T =.0; double[] times = new double [00]; 8 for ( int countertime =0; countertime < numberoftimes ; countertime ++ { times [ countertime ] = t0 + countertime *( T- t0/(( double numberoftimes ; 30 3 //brownian motion TimeDiscretization mytimes = new TimeDiscretization ( times ; BrownianMotion W = new BrownianMotion ( mytimes,1, numberofpaths, seed ; RandomVariableInterface previousvalueybar = new RandomVariable ( t0, initialvalueybar ; 38 RandomVariableInterface initialx = new RandomVariable ( t0, initialvaluex ; System. out. println (" time "+ " "+" mean value Y"+ " "+ " mean value Ybar " +" "+ " theor. mean "+ " "+ " error Y"+ " "+" error Ybar " + " "+ " Y - Ybar " ; 40 for ( int countertime =0; countertime < numberoftimes -1; countertime ++ { 4 double dt = times [ countertime +1] - times [ countertime ]; RandomVariableInterface dw = W. getbrownianincrement ( countertime, 0;//.getIncrementAtIndex( countertime; 44 RandomVariableInterface nextvalueybar = previousvalueybar. mult (1+( mu +0.5* sigma * sigma *dt. add (dw. mult ( previousvalueybar. mult ( sigma ; double empiricalmeanybar = nextvalueybar. getaverage
3 46 (; RandomVariableInterface nextvaluex = initialx. add (( mu* times [ countertime ]. add (W. getbrownianincrement ( countertime, 0. mult ( sigma ; 48 RandomVariableInterface nextvaluey = nextvaluex. exp (; double empiricalmeany = nextvaluey. getaverage (; 50 double theoreticalmean = initialvalueybar * Math. exp (( mu +0.5* sigma * sigma * times [ countertime ]; //continuous model 5 double absdifference = Math. abs ( empiricalmeanybar - empiricalmeany ; 54 double errory = Math. abs ( empiricalmeanybar - theoreticalmean ; double errorybar = Math. abs ( empiricalmeany - theoreticalmean ; 56 System. out. println ( formattersci6. format ( times [ countertime ] + " " + formattersci6. format ( empiricalmeany + " " + formattersci6. format ( empiricalmeanybar 58 + " " + formattersci6. format ( theoreticalmean + " "+ formattersci6. format ( errory + " "+ formattersci6. format ( errorybar +" " + formattersci6. format ( absdifference ; previousvalueybar = new RandomVariable ( nextvalueybar ; Partial solution: We can express the approximations Ȳ in the following Euler-Maruyama recursion Ȳ n = t Ȳ +σȳ W for n =, 3,... with Ȳ0 = e x. From the solution to the previous exercise sheet we have that ( V ar[ȳn] = e x ( t l + σ t l 3 t l
4 and if we take t l, l, ( V ar[ȳn] =e x ( + σ ( ( n =e x + σ ( ( n n =e x (1 + (µ + σ σ + (µ + + σ ( ( n n =e x 1 + µ + σ + (µ + σ + σ ( ( n n =e x 1 + (µ + σ + (µ + σ n In order to compute V ar[y n ] = V ar[e Xn ], with X n = X + µ t + σ W we need to know the first and second moments of e Xn. Note that X n = x + µt n + σ W l therefore E [ e Xn] = e x+µtn E [e σ ] W l = e x+µtn E [ e Z] where Z := σ W l is just an N (0, σ n = N (0, σ t n distributed r.v. hence E [ e Z] = e σ tn and E [ e Xn] σ x+(µ+ = e tn. Now, computing E[(e Xn ] = E[e Xn ], we have Note that hence and finally E [ e Xn] = e x+µtn E [e σ ] W l = e x+µtn E [ e Z]. V ar[y n ] =V ar[e Xn ] E [ e Z] = e σ t n, E [ e Xn] = e x+(µ+σ t n =E[e Xn ] (E[e Xn ] =e x+(µ+σ σ t n x+(µ+ (e tn =e x+(µ+σ σ t n x+(µ+ e tn σ x+(µ+ =e tn (e σ t n 1 4
5 Now, we can compare V ar[ȳn] and V ar[y n ]. We note that, since e x is a common factor of both, it is enough to compare V ar[yn] e V ar[ȳn] x e with 0. We have that x n e (µ+σ t n (1 + (µ + σ + (µ + σ > e (µ+σ t n ( 1 + (µ + σ + (µ + σ n = e (µ+σ t n (1 + (µ + σ n n n ( = e (µ+σ t n 1 + (µ + σ t n n n > 0. Note that the right hand term is a lower bound for the left hand term (and its limit as n is the left hand term. Similarly n σ (µ+ e tn < 0. It remains to show whether n t n (1 + (µ + σ + (µ + σ n e(µ+σ < σ (µ+ e tn or not. If it is true, we conclude that V ar[y n ] e x V ar[ȳn] e x < 0, i.e. V ar[y n ] < V ar[ȳn], otherwise V ar[y n ] > V ar[ȳn]. Exercise The example Java program containing the solution is in provided in the following files package de. math. lmu. exercise08 ; 4 import java. text. DecimalFormat ; import net. finmath. montecarlo. BrownianMotion ; 6 import net. finmath. montecarlo. RandomVariable ; import net. finmath. stochastic. RandomVariableInterface ; 8 import net. finmath. time. TimeDiscretization ; 10 public class OptionsTest { 1 14 s t a t i c f i n a l DecimalFormat formattersci6 = new DecimalFormat (" E00 ; E00 "; public s t a t i c void main ( String [] args { 16 int numberoftimes = 800; 5
6 double t0 =0.0; 18 double T = 1.0; double initialvalue = 1.0; 0 double mu = 0.1; double r = 0.0; double sigma = 1.0; double strikes [] = {0.8,0.9,1.0,1.1,1.; 4 double price [] = new double[ strikes. length ]; double riskyvalue [] = new double[ strikes. length ]; 6 double pricebs [] = new double[ strikes. length ]; 8 double[] times = new double[ numberoftimes ]; for ( int countertime =0; countertime < numberoftimes ; countertime ++ { 30 times [ countertime ] = t0 + countertime *( T- t0/(( double numberoftimes ; 3 int seed = 34; 34 int numberofpaths = 10000; 36 TimeDiscretization mytimes = new TimeDiscretization ( times ; 38 BrownianMotion brownianw = new BrownianMotion ( mytimes,1, numberofpaths, seed ; //RandomVariableInterface initialxbar = new RandomVariable (t0,0.0; 40 RandomVariableInterface previousvaluep = new RandomVariable ( t0, initialvalue ; 4 RandomVariableInterface previousvalueq = new RandomVariable ( t0, initialvalue ; System. out. println (" strike "+ " "+" risky value "+ " "+"MC price " + " "+ "BS price "; for ( int countertime =0; countertime < numberoftimes -1; countertime ++ { 48 double dt = times [ countertime +1] - times [ countertime ]; RandomVariableInterface dw = brownianw. getbrownianincrement ( countertime, 0; 50 RandomVariableInterface nextvaluep = previousvaluep. mult (1+ mu*dt. add (dw. mult ( previousvaluep. mult ( sigma ; 5 RandomVariableInterface nextvalueq = previousvalueq. mult (1+ r*dt. add (dw. mult ( previousvalueq. mult ( sigma ; //continuous model // double MCerror = Math.abs(empiricalMean theoreticalmean; 58 previousvaluep = new RandomVariable ( nextvaluep ; previousvalueq = new RandomVariable ( nextvalueq ; 60 //System.out.println(formatterSci6.format(times[ countertime]+ " "+formattersci6.format( empiricalmean+ " "+formattersci6.format( 6
7 6 theoreticalmean+ " ; "+formattersci6.format(mcerror for ( int counterstrike = 0; counterstrike < strikes. length ; counterstrike ++ { 70 price [ counterstrike ] = Math. exp (-r*(t-t0* previousvalueq. sub ( strikes [ counterstrike ]. floor (0. getaverage (; riskyvalue [ counterstrike ] = Math. exp (-r*(t-t0* previousvaluep. sub ( strikes [ counterstrike ]. floor (0. getaverage (; 7 74 pricebs [ counterstrike ] = net. finmath. functions. AnalyticFormulas. blackscholesoptionvalue ( initialvalue, 76 r, sigma, 78 T, strikes [ counterstrike ]; System. out. println ( formattersci6. format ( strikes [ counterstrike ] +" " + formattersci6. format ( riskyvalue [ counterstrike ] + " " + formattersci6. format ( price [ counterstrike ] + " " + formattersci6. format ( pricebs [ counterstrike ]; 86 package de. math. lmu. exercise08 ; 4 import java. text. DecimalFormat ; import net. finmath. montecarlo. BrownianMotion ; 6 import net. finmath. montecarlo. RandomVariable ; import net. finmath. stochastic. RandomVariableInterface ; 8 import net. finmath. time. TimeDiscretization ; 10 public class OptionsTest { 1 14 s t a t i c f i n a l DecimalFormat formattersci6 = new DecimalFormat (" E00 ; E00 "; public s t a t i c void main ( String [] args { 16 int numberoftimes = 800; 18 double t0 =0.0; double T = 1.0; 0 double initialvalue = 1.0; double mu = 0.1; double r = 0.01; double sigma = 1.0; 7
8 4 double strikes [] = {0.8,0.9,1.0,1.1,1.; double price [] = new double[ strikes. length ]; 6 double riskyvalue [] = new double[ strikes. length ]; double pricebs [] = new double[ strikes. length ]; 8 double[] times = new double[ numberoftimes ]; for ( int countertime = 0; countertime < numberoftimes ; countertime ++ { 30 times [ countertime ] = t0 + countertime *( T- t0/(( double numberoftimes ; 3 int seed = 34; 34 int numberofpaths = 10000; 36 TimeDiscretization mytimes = new TimeDiscretization ( times ; 38 System. out. println ( mytimes. getnumberoftimes (; 40 BrownianMotion brownianw = new BrownianMotion ( mytimes,1, numberofpaths, seed ; RandomVariableInterface initialxbar = new RandomVariable ( t0,0.0 ; 4 44 RandomVariableInterface myw = new RandomVariable (0.0 ; for ( int i = 0; i < mytimes. getnumberoftimesteps (; i ++ { 46 myw = myw. add ( brownianw. getbrownianincrement (i, 0; System. out. println (" strike "+ " "+" risky value "+ " "+"MC price " + " "+ "BS price "; 5 RandomVariableInterface valuelogxbarp = initialxbar. add ((( mu -0.5* sigma * sigma * times [ numberoftimes -1]. add ( myw. mult ( sigma ; 54 RandomVariableInterface valuexbarp = valuelogxbarp. exp (; 56 RandomVariableInterface valuelogxbarq = initialxbar. add ((( r -0.5* sigma * sigma * times [ numberoftimes -1]. add ( myw. mult ( sigma ; RandomVariableInterface valuexbarq = valuelogxbarq. exp (; for ( int counterstrike = 0; counterstrike < strikes. length ; counterstrike ++ { 6 price [ counterstrike ] = Math. exp (-r*(t-t0* valuexbarq. sub ( strikes [ counterstrike ]. floor (0. getaverage (; riskyvalue [ counterstrike ] = Math. exp (-r*(t-t0* valuexbarp. sub ( strikes [ counterstrike ]. floor (0. getaverage (; 64 pricebs [ counterstrike ] = net. finmath. functions. AnalyticFormulas. blackscholesoptionvalue ( 66 initialvalue, r, 68 sigma, T, 70 strikes [ counterstrike ]; 8
9 7 System. out. println ( formattersci6. format ( strikes [ counterstrike ] +" " + formattersci6. format ( riskyvalue [ counterstrike ] + " " + formattersci6. format ( price [ counterstrike ] + " " + formattersci6. format ( pricebs [ counterstrike ]; Each uses a different approach to approximate S, the first uses the Euler- Maruyama scheme on ds and the second computes the exact value of S using the Euler-Maruyama scheme on logs. (The second approach is faster. The correct measure to compute is the pricing measure Q. In order to implement this measure by Monte-Carlo, the empirical mean must be computed on the MC paths simulated under the Q-model. The Q-model for S is ds t = r t S t dt + σs t dw Q t and since r t 0, we have the following Q-model for S ds t = σs t dw Q t. 9
Exercise Sheet 7 - Solutions
Exercise Sheet 7 - Solutions Alessandro Gnoatto June 10, 015 1 Exercise 1 1. Let Y t = log X t, then dy t = 1 dx t 1 X t Xt Xt σ dt by Ito s Lemma. Hence This implies that dy t = (µ σ )dt + σdw t. log
More informationThe 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 informationStochastic 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 informationStochastic 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 informationMA8109 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 informationStochastic 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 informationLecture 12: Diffusion Processes and Stochastic Differential Equations
Lecture 12: Diffusion Processes and Stochastic Differential Equations 1. Diffusion Processes 1.1 Definition of a diffusion process 1.2 Examples 2. Stochastic Differential Equations SDE) 2.1 Stochastic
More informationQuestion 1. The correct answers are: (a) (2) (b) (1) (c) (2) (d) (3) (e) (2) (f) (1) (g) (2) (h) (1)
Question 1 The correct answers are: a 2 b 1 c 2 d 3 e 2 f 1 g 2 h 1 Question 2 a Any probability measure Q equivalent to P on F 2 can be described by Q[{x 1, x 2 }] := q x1 q x1,x 2, 1 where q x1, q x1,x
More informationMalliavin 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 informationIntroduction 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 informationDUBLIN CITY UNIVERSITY
DUBLIN CITY UNIVERSITY SAMPLE EXAMINATIONS 2017/2018 MODULE: QUALIFICATIONS: Simulation for Finance MS455 B.Sc. Actuarial Mathematics ACM B.Sc. Financial Mathematics FIM YEAR OF STUDY: 4 EXAMINERS: Mr
More informationPartial Differential Equations with Applications to Finance Seminar 1: Proving and applying Dynkin s formula
Partial Differential Equations with Applications to Finance Seminar 1: Proving and applying Dynkin s formula Group 4: Bertan Yilmaz, Richard Oti-Aboagye and Di Liu May, 15 Chapter 1 Proving Dynkin s formula
More informationI 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 informationStochastic Differential Equations
Chapter 5 Stochastic Differential Equations We would like to introduce stochastic ODE s without going first through the machinery of stochastic integrals. 5.1 Itô Integrals and Itô Differential Equations
More informationBrownian 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 informationShort-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility
Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,
More informationNumerical 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 informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationThe moment-generating function of the log-normal distribution using the star probability measure
Noname manuscript No. (will be inserted by the editor) The moment-generating function of the log-normal distribution using the star probability measure Yuri Heymann Received: date / Accepted: date Abstract
More informationBernardo 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 informationStochastic 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 informationChapter 4: Monte-Carlo Methods
Chapter 4: Monte-Carlo Methods A Monte-Carlo method is a technique for the numerical realization of a stochastic process by means of normally distributed random variables. In financial mathematics, it
More information1 Brownian Local Time
1 Brownian Local Time We first begin by defining the space and variables for Brownian local time. Let W t be a standard 1-D Wiener process. We know that for the set, {t : W t = } P (µ{t : W t = } = ) =
More informationLecture 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 informationLECTURE 2: LOCAL TIME FOR BROWNIAN MOTION
LECTURE 2: LOCAL TIME FOR BROWNIAN MOTION We will define local time for one-dimensional Brownian motion, and deduce some of its properties. We will then use the generalized Ray-Knight theorem proved in
More informationKolmogorov 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 informationExact Simulation of Diffusions and Jump Diffusions
Exact Simulation of Diffusions and Jump Diffusions A work by: Prof. Gareth O. Roberts Dr. Alexandros Beskos Dr. Omiros Papaspiliopoulos Dr. Bruno Casella 28 th May, 2008 Content 1 Exact Algorithm Construction
More informationlim n C1/n n := ρ. [f(y) f(x)], y x =1 [f(x) f(y)] [g(x) g(y)]. (x,y) E A E(f, f),
1 Part I Exercise 1.1. Let C n denote the number of self-avoiding random walks starting at the origin in Z of length n. 1. Show that (Hint: Use C n+m C n C m.) lim n C1/n n = inf n C1/n n := ρ.. Show that
More informationDISCRETE-TIME STOCHASTIC MODELS, SDEs, AND NUMERICAL METHODS. Ed Allen. NIMBioS Tutorial: Stochastic Models With Biological Applications
DISCRETE-TIME STOCHASTIC MODELS, SDEs, AND NUMERICAL METHODS Ed Allen NIMBioS Tutorial: Stochastic Models With Biological Applications University of Tennessee, Knoxville March, 2011 ACKNOWLEDGEMENT I thank
More informationBrownian Motion and Stochastic Calculus
ETHZ, Spring 17 D-MATH Prof Dr Martin Larsson Coordinator A Sepúlveda Brownian Motion and Stochastic Calculus Exercise sheet 6 Please hand in your solutions during exercise class or in your assistant s
More informationAffine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009
Affine Processes Econometric specifications Eduardo Rossi University of Pavia March 17, 2009 Eduardo Rossi (University of Pavia) Affine Processes March 17, 2009 1 / 40 Outline 1 Affine Processes 2 Affine
More informationTheoretical Tutorial Session 2
1 / 36 Theoretical Tutorial Session 2 Xiaoming Song Department of Mathematics Drexel University July 27, 216 Outline 2 / 36 Itô s formula Martingale representation theorem Stochastic differential equations
More informationExpectation, variance and moments
Expectation, variance and moments John Appleby Contents Expectation and variance Examples 3 Moments and the moment generating function 4 4 Examples of moment generating functions 5 5 Concluding remarks
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing
More informationNumerical 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 informationNumerical Simulation of Stochastic Differential Equations: Lecture 2, Part 1. Recap: SDE. Euler Maruyama. Lecture 2, Part 1: Euler Maruyama
Numerical Simulation of Stochastic Differential Equations: Lecture, Part 1 Des Higham Department of Mathematics Universit of Strathclde Lecture, Part 1: Euler Maruama Definition of Euler Maruama Method
More informationStochastic 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 informationAdaptive timestepping for SDEs with non-globally Lipschitz drift
Adaptive timestepping for SDEs with non-globally Lipschitz drift Mike Giles Wei Fang Mathematical Institute, University of Oxford Workshop on Numerical Schemes for SDEs and SPDEs Université Lille June
More informationStochastic 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 informationIntroduction to the Numerical Solution of SDEs Selected topics
0 / 23 Introduction to the Numerical Solution of SDEs Selected topics Andreas Rößler Summer School on Numerical Methods for Stochastic Differential Equations at Vienna University of Technology, Austria
More informationOptimal portfolio strategies under partial information with expert opinions
1 / 35 Optimal portfolio strategies under partial information with expert opinions Ralf Wunderlich Brandenburg University of Technology Cottbus, Germany Joint work with Rüdiger Frey Research Seminar WU
More informationExercises in stochastic analysis
Exercises in stochastic analysis Franco Flandoli, Mario Maurelli, Dario Trevisan The exercises with a P are those which have been done totally or partially) in the previous lectures; the exercises with
More informationDerivation of Itô SDE and Relationship to ODE and CTMC Models
Derivation of Itô SDE and Relationship to ODE and CTMC Models Biomathematics II April 23, 2015 Linda J. S. Allen Texas Tech University TTU 1 Euler-Maruyama Method for Numerical Solution of an Itô SDE dx(t)
More informationIntroduction. Stochastic Processes. Will Penny. Stochastic Differential Equations. Stochastic Chain Rule. Expectations.
19th May 2011 Chain Introduction We will Show the relation between stochastic differential equations, Gaussian processes and methods This gives us a formal way of deriving equations for the activity of
More informationStochastic 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 informationMoment Generating Functions
MATH 382 Moment Generating Functions Dr. Neal, WKU Definition. Let X be a random variable. The moment generating function (mgf) of X is the function M X : R R given by M X (t ) = E[e X t ], defined for
More informationWeak solutions of mean-field stochastic differential equations
Weak solutions of mean-field stochastic differential equations Juan Li School of Mathematics and Statistics, Shandong University (Weihai), Weihai 26429, China. Email: juanli@sdu.edu.cn Based on joint works
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationIntroduction 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 informationB8.3 Mathematical Models for Financial Derivatives. Hilary Term Solution Sheet 2
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
More informationStochastic 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 informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations Simo Särkkä Aalto University, Finland November 18, 2014 Simo Särkkä (Aalto) Lecture 4: Numerical Solution of SDEs November
More informationPopulations Models: Part I
Populations Models: Part I Phil. Pollett UQ School of Maths and Physics Mathematics Enrichment Seminars 26 August 2015 Phil. Pollett (UQ School of Maths and Physics) Populations Models: Part I 1 / 48 Growth
More informationOn the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem
On the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem Koichiro TAKAOKA Dept of Applied Physics, Tokyo Institute of Technology Abstract M Yor constructed a family
More informationExact Simulation of Multivariate Itô Diffusions
Exact Simulation of Multivariate Itô Diffusions Jose Blanchet Joint work with Fan Zhang Columbia and Stanford July 7, 2017 Jose Blanchet (Columbia/Stanford) Exact Simulation of Diffusions July 7, 2017
More informationIEOR 4701: Stochastic Models in Financial Engineering. Summer 2007, Professor Whitt. SOLUTIONS to Homework Assignment 9: Brownian motion
IEOR 471: Stochastic Models in Financial Engineering Summer 27, Professor Whitt SOLUTIONS to Homework Assignment 9: Brownian motion In Ross, read Sections 1.1-1.3 and 1.6. (The total required reading there
More informationThe multidimensional Ito Integral and the multidimensional Ito Formula. Eric Mu ller June 1, 2015 Seminar on Stochastic Geometry and its applications
The multidimensional Ito Integral and the multidimensional Ito Formula Eric Mu ller June 1, 215 Seminar on Stochastic Geometry and its applications page 2 Seminar on Stochastic Geometry and its applications
More informationLecture on Parameter Estimation for Stochastic Differential Equations. Erik Lindström
Lecture on Parameter Estimation for Stochastic Differential Equations Erik Lindström Recap We are interested in the parameters θ in the Stochastic Integral Equations X(t) = X(0) + t 0 µ θ (s, X(s))ds +
More informationHOMEWORK 4 1. P45. # 1.
HOMEWORK 4 SHUANGLIN SHAO P45 # Proof By the maximum principle, u(x, t x kt attains the maximum at the bottom or on the two sides When t, x kt x attains the maximum at x, ie, x When x, x kt kt attains
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations Simo Särkkä Aalto University Tampere University of Technology Lappeenranta University of Technology Finland November
More informationA 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 informationMATH4210 Financial Mathematics ( ) Tutorial 7
MATH40 Financial Mathematics (05-06) Tutorial 7 Review of some basic Probability: The triple (Ω, F, P) is called a probability space, where Ω denotes the sample space and F is the set of event (σ algebra
More informationA MODEL FOR THE LONG-TERM OPTIMAL CAPACITY LEVEL OF AN INVESTMENT PROJECT
A MODEL FOR HE LONG-ERM OPIMAL CAPACIY LEVEL OF AN INVESMEN PROJEC ARNE LØKKA AND MIHAIL ZERVOS Abstract. We consider an investment project that produces a single commodity. he project s operation yields
More informationUniversity Of Calgary Department of Mathematics and Statistics
University Of Calgary Department of Mathematics and Statistics Hawkes Seminar May 23, 2018 1 / 46 Some Problems in Insurance and Reinsurance Mohamed Badaoui Department of Electrical Engineering National
More informationSome 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 informationIntroduction to Diffusion Processes.
Introduction to Diffusion Processes. Arka P. Ghosh Department of Statistics Iowa State University Ames, IA 511-121 apghosh@iastate.edu (515) 294-7851. February 1, 21 Abstract In this section we describe
More informationMoment 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 information1. 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 informationBernardo 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 informationBridging the Gap between Center and Tail for Multiscale Processes
Bridging the Gap between Center and Tail for Multiscale Processes Matthew R. Morse Department of Mathematics and Statistics Boston University BU-Keio 2016, August 16 Matthew R. Morse (BU) Moderate Deviations
More informationTMS165/MSA350 Stochastic Calculus, Lecture on Applications
TMS165/MSA35 Stochastic Calculus, Lecture on Applications In this lecture we demonstrate how statistical methods such as the maximum likelihood method likelihood ratio estimation can be applied to the
More informationThe Multivariate Normal Distribution 1
The Multivariate Normal Distribution 1 STA 302 Fall 2017 1 See last slide for copyright information. 1 / 40 Overview 1 Moment-generating Functions 2 Definition 3 Properties 4 χ 2 and t distributions 2
More informationLecture 21: Stochastic Differential Equations
: Stochastic Differential Equations In this lecture, we study stochastic differential equations. See Chapter 9 of [3] for a thorough treatment of the materials in this section. 1. Stochastic differential
More informationHJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011
Department of Probability and Mathematical Statistics Faculty of Mathematics and Physics, Charles University in Prague petrasek@karlin.mff.cuni.cz Seminar in Stochastic Modelling in Economics and Finance
More informationPoisson random measure: motivation
: motivation The Lévy measure provides the expected number of jumps by time unit, i.e. in a time interval of the form: [t, t + 1], and of a certain size Example: ν([1, )) is the expected number of jumps
More informationMultilevel Monte Carlo for Lévy Driven SDEs
Multilevel Monte Carlo for Lévy Driven SDEs Felix Heidenreich TU Kaiserslautern AG Computational Stochastics August 2011 joint work with Steffen Dereich Philipps-Universität Marburg supported within DFG-SPP
More informationExact and high order discretization schemes. Wishart processes and their affine extensions
for Wishart processes and their affine extensions CERMICS, Ecole des Ponts Paris Tech - PRES Université Paris Est Modeling and Managing Financial Risks -2011 Plan 1 Motivation and notations 2 Splitting
More informationStochastic Numerical Analysis
Stochastic Numerical Analysis Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Stoch. NA, Lecture 3 p. 1 Multi-dimensional SDEs So far we have considered scalar SDEs
More informationUses of Asymptotic Distributions: In order to get distribution theory, we need to norm the random variable; we usually look at n 1=2 ( X n ).
1 Economics 620, Lecture 8a: Asymptotics II Uses of Asymptotic Distributions: Suppose X n! 0 in probability. (What can be said about the distribution of X n?) In order to get distribution theory, we need
More informationMean-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 informationStochastic Models (Lecture #4)
Stochastic Models (Lecture #4) Thomas Verdebout Université libre de Bruxelles (ULB) Today Today, our goal will be to discuss limits of sequences of rv, and to study famous limiting results. Convergence
More informationSTAT/MATH 395 A - PROBABILITY II UW Winter Quarter Moment functions. x r p X (x) (1) E[X r ] = x r f X (x) dx (2) (x E[X]) r p X (x) (3)
STAT/MATH 395 A - PROBABILITY II UW Winter Quarter 07 Néhémy Lim Moment functions Moments of a random variable Definition.. Let X be a rrv on probability space (Ω, A, P). For a given r N, E[X r ], if it
More informationDerivation of SPDEs for Correlated Random Walk Transport Models in One and Two Dimensions
This article was downloaded by: [Texas Technology University] On: 23 April 2013, At: 07:52 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationStochastic 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 informationDiscretization 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 information1 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 informationStochastic Computation
Stochastic Computation A Brief Introduction Klaus Ritter Computational Stochastics TU Kaiserslautern 1/1 Objectives Objectives of Computational Stochastics: Construction and analysis of (A) algorithms
More informationGeneralised Fractional-Black-Scholes Equation: pricing and hedging
Generalised Fractional-Black-Scholes Equation: pricing and hedging Álvaro Cartea Birkbeck College, University of London April 2004 Outline Lévy processes Fractional calculus Fractional-Black-Scholes 1
More informationVariance 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 informationStochastic 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 informationMonte-Carlo MMD-MA, Université Paris-Dauphine. Xiaolu Tan
Monte-Carlo MMD-MA, Université Paris-Dauphine Xiaolu Tan tan@ceremade.dauphine.fr Septembre 2015 Contents 1 Introduction 1 1.1 The principle.................................. 1 1.2 The error analysis
More informationON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER
ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER GERARDO HERNANDEZ-DEL-VALLE arxiv:1209.2411v1 [math.pr] 10 Sep 2012 Abstract. This work deals with first hitting time densities of Ito processes whose
More informationJUST THE MATHS UNIT NUMBER LAPLACE TRANSFORMS 3 (Differential equations) A.J.Hobson
JUST THE MATHS UNIT NUMBER 16.3 LAPLACE TRANSFORMS 3 (Differential equations) by A.J.Hobson 16.3.1 Examples of solving differential equations 16.3.2 The general solution of a differential equation 16.3.3
More informationLinear 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 informationGARCH processes continuous counterparts (Part 2)
GARCH processes continuous counterparts (Part 2) Alexander Lindner Centre of Mathematical Sciences Technical University of Munich D 85747 Garching Germany lindner@ma.tum.de http://www-m1.ma.tum.de/m4/pers/lindner/
More informationLeast Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises
Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Hongwei Long* Department of Mathematical Sciences, Florida Atlantic University, Boca Raton Florida 33431-991,
More informationData Analysis and Machine Learning Lecture 10: Applications of Kalman Filtering
Data Analysis and Machine Learning Lecture 10: Applications of Kalman Filtering Lecturer, Imperial College London Founder, CEO, Thalesians Ltd 2017.11.08 The Newtonian system I Let r be the position of
More informationLecture 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 informationMultilevel Monte Carlo for Stochastic McKean-Vlasov Equations
Multilevel Monte Carlo for Stochastic McKean-Vlasov Equations Lukasz Szpruch School of Mathemtics University of Edinburgh joint work with Shuren Tan and Alvin Tse (Edinburgh) Lukasz Szpruch (University
More informationParameter Estimation for Random Differential Equation Models
Parameter Estimation for Random Differential Equation Models H.T. Banks and Michele L. Joyner Center for Research in Scientific Computation North Carolina State University Raleigh, NC, United States and
More information