Higher order weak approximations of stochastic differential equations with and without jumps
|
|
- Esmond Dixon
- 5 years ago
- Views:
Transcription
1 Higher order weak approximations of stochastic differential equations with and without jumps Hideyuki TANAKA Graduate School of Science and Engineering, Ritsumeikan University Rough Path Analysis and Related Topics, Nagoya University January 26, / 37
2 Outline 1 Introduction 2 Approximations of SDEs without jumps Approximation of SDEs by random ODEs (Extended) Cubature on Wiener space Operator Splitting method as a Cubature 3 Approximations of SDEs with jumps I (finite jump intesity) Operator Splitting Jump-adapted approximation 4 Approximations of SDEs with jumps II (general Lévy process) Overview Small jump approximation How to deal with compound Poisson term 5 Further Research 2 / 37
3 Introduction Introduction 1 Introduction 2 Approximations of SDEs without jumps Approximation of SDEs by random ODEs (Extended) Cubature on Wiener space Operator Splitting method as a Cubature 3 Approximations of SDEs with jumps I (finite jump intesity) Operator Splitting Jump-adapted approximation 4 Approximations of SDEs with jumps II (general Lévy process) Overview Small jump approximation How to deal with compound Poisson term 5 Further Research 3 / 37
4 Introduction Topics 1 We want to know the value of the expectation E[f (X T )] where X t is a solution of a stochastic differential equation (SDE) 2 Applications: finance, control, filtering, physics, etc 3 However, we do not know the exact distribution of X t in general Therefore, we can not simulate X t by MonteCarlo 4 Goal: find a higher order approximation scheme X (n) t E[f (X T )] E[f (X (n) T )] = O(n k ) st (especially, k 2) 4 / 37
5 Introduction What are merits of higher order methods? Reducing computational cost (eg k = 1, n = 1 k = 2, n =several tens) Quasi Monte Carlo: It (sometimes) holds that E[f (X (n) T )] = g(y)p(dy) R α(n) g : R α(n) R p : measure on R α(n) If α(n) is not so large ( 1, 1?), QMC works well 5 / 37
6 Approximations of SDEs without jumps Approximations of SDEs without jumps 1 Introduction 2 Approximations of SDEs without jumps Approximation of SDEs by random ODEs (Extended) Cubature on Wiener space Operator Splitting method as a Cubature 3 Approximations of SDEs with jumps I (finite jump intesity) Operator Splitting Jump-adapted approximation 4 Approximations of SDEs with jumps II (general Lévy process) Overview Small jump approximation How to deal with compound Poisson term 5 Further Research 6 / 37
7 Approximations of SDEs without jumps Approximation of SDEs by random ODEs Setting Consider a Stratonovich SDE (R N -valued) X x t = x + d i= t V i (X x s ) dw i (s) (1) (W i ) 1 i d is a d-dimensional Brownian motion W (s) = s V i C b (RN ; R N ) Our purpose: Compute P t f (x) := E[f (X x t )] at time t = T, through time discretization methods 7 / 37
8 Approximations of SDEs without jumps Approximation of SDEs by random ODEs Random Ordinary Differential Equation Let t > be a fixed time-scaling parameter We want to construct an approximation process for r t X x r = x + d i= r V i (X x s)dω t i (s) (2) a driving random path ω has bounded variation paths ω W in some sense ω depends on fixed (small) time t; 8 / 37
9 Approximations of SDEs without jumps Approximation of SDEs by random ODEs Approximation Operator Approximation operator: For a scaling parameter t, Q t f (x) := E[f (X x t )] Time partition: π := { = t < t 1 < < t n = T } Operator Q t1 Q t2 t 1 Q tn t n 1 Markov chain whose transition follows from the random ODE (2) with scaling t = t j+1 t j Random ODE, defined in [, T ], driven by (ωπ) i Here ωπ(s)( i s T ) is defined for each time interval [t j, t j+1 ] independently: X t j,x r = x + d i= r ((ω tj+1 t j ) j=,,n 1 are independent) t j V i (X tj,x s )dω i t j+1 t j (s), r [t j, t j+1 ] 9 / 37
10 Approximations of SDEs without jumps Approximation of SDEs by random ODEs From local error to global error (Rough sketch): t j t j 1 T /n 1 Assume (P t Q t )g = O(t k+1 ) for smooth g 2 Then (P T QT n /n )f = (Pn T /n Qn T /n )f n 1 = Q j t/n (P T /n Q T /n )P (1 (j+1)/n)t f j= = n O(n (k+1) ) = O(n k ) 3 f : smooth & polynomial growth; see eg Talay & Tubaro(199), Tanaka & Kohatsu-Higa (29) 4 f : Lipschitz or more general, under Hörmander type condition; see Bally & Talay, Kohatsu-Higa, Kusuoka, etc 1 / 37
11 Approximations of SDEs without jumps Approximation of SDEs by random ODEs Stochastic Taylor expansion I By using Itô formula, (V i act as vector fields) f (X x t ) = f (x) + = f (x) + + i,j d d i= t d (V i f )(x) i= t s (V i f )(X x s ) dw i (s) t Notation for the order of convergence: For index (i 1,, i l ), i j =, 1,, d, dw i (s) (V j V i f )(X x r ) dw j (r) dw i (s) (i 1,, i l ) := l + #{j : i j = } 11 / 37
12 Approximations of SDEs without jumps Approximation of SDEs by random ODEs Stochastic Taylor expansion I (continued) Stochastic Taylor expansion: f (Xt x ) = (V i1 V il f )(x) (i 1,,i l ) m t s2 + R X m(t, f ) dw i 1 (s 1 ) dw i l (s l ) with E[ R X m(t, f ) ] C(T, f ) t (m+1)/2 A similar expansion holds for X x t ( dw is replaced by dω) 12 / 37
13 Higher order weak approximations of stochastic differential equations with and without jumps Approximations of SDEs without jumps (Extended) Cubature on Wiener space Cubature formulas (extended) 1 Definition: ω C,bv ([, 1]; R d+1 ) defines a cubature formula with degree m if for any index (i 1,, i l ) with (i 1,, i l ) m, [ 1 E s2 [ 1 = E s2 and for (i 1,, i l ) = m + 1 or m + 2, 1 ] dw i 1 (s 1 ) dw i l (s l ) ] dω i 1 (s 1 ) dω i l (s l ), s2 d ω i 1 (s 1 ) d ω i l (s l ) L 2 (Ω) 2 Lyons & Victoir (24) assume in addition, (We do not assume here) ω (s) = s ω has a discrete probability distribution on C,bv ([, 1]; R d+1 ) 13 / 37
14 Approximations of SDEs without jumps (Extended) Cubature on Wiener space Time Scaling 1 Note that t (i 1,,i l ) 1 t law = s2 2 Scaling for ω: For s t, s2 dw i 1 (s 1 ) dw i l (s l ) dw i 1 (s 1 ) dw i l (s l ) ω t (s) := tω (s/t) ω i t(s) := tω i t(s/t), 1 i d 14 / 37
15 Approximations of SDEs without jumps (Extended) Cubature on Wiener space Error estimates (Sketch) Assume that ω defines a cubature formula with degree m f C b Then E[f (Xt x )] E[f (X x t )] = E[Rm(t, X f )] E[Rm(t, X f )] = O(t (m+1)/2 ) Order k = (m 1)/2 scheme! For second order schemes, we need m = 5 15 / 37
16 Approximations of SDEs without jumps Operator Splitting method as a Cubature Stochastic Taylor expansion II Formal Taylor series of t e tl By taking expectations of stochastic Taylor expansions, E[f (X x t )] = f (x) + tlf (x) + t 2 2 L2 f (x) + t 3 3! L3 f (x) + where L := d i= L i ( the generator of P t ), L := V L i = 1 V 2 2 i, 1 i d Approximation of P t Approximation of exponential map 16 / 37
17 Approximations of SDEs without jumps Operator Splitting method as a Cubature Splitting of exponential maps on noncommutative algebra 1 Now we can not solve directly e tl = e t(l +L 1 + +L d ) 2 However, each e tl i may be solvable Approximation of e tl by {e tl,, e tl d } {L i } : noncommutative 17 / 37
18 Approximations of SDEs without jumps Operator Splitting method as a Cubature Second order method via splitting e t(l i +L j ) = I + t(l i + L j ) + t 2 2 (L i + L j ) 2 + e tl i e tl j = I + t(l i + L j ) + t2 2 (L2 i + L 2 j + 2L i L j ) + Socond-order methods: e t(l i +L j ) e t(l i +L j ) = 1 2 etl i e tl j etl j e tl i + O(t 3 ) = e t/2l i e tl j e t/2l i + O(t 3 ) The idea can be found in the works of Strang (196s) 18 / 37
19 Approximations of SDEs without jumps Operator Splitting method as a Cubature ODE and coordinate SDE What s e tl i? Notation: exp(v )x is the solution of at time t = 1 dz t (x) dt = V (z t (x)), z (x) = x exp(w i t V i)x is the solution of the SDE dx x i,t = x + t V i (X x i,s) dw i (s) R N So, e tl i is solvable in the sense of e tl i f (x) = E[f (X x i,t)] = E[f (exp(w i t V i )x)] e tl i e tl j two step flow exp(w j t V j) exp(w i t V i)x 19 / 37
20 Approximations of SDEs without jumps Operator Splitting method as a Cubature Ninomiya-Victoir scheme Ninomiya-Victoir scheme exp(t/2v ) exp( tz 1 V 1 ) exp( tz d V d ) exp(t/2v )x, X x if Λ = 1, t = exp(t/2v ) exp( tz d V d ) exp( tz 1 V 1 ) exp(t/2v )x, if Λ = 1 Z = (Z i ) 1 i d and Λ are independent Z N(, I d ) P(Λ = ±1) = 1/2 If the ODE has no closed-form solution, then we can use Runge-Kutta methods 2 / 37
21 Approximations of SDEs without jumps Operator Splitting method as a Cubature Ninomiya-Victoir scheme as an extended cubature on Wiener space Ninomiya-Victoir as cubature ω i (r) = r ηi (s)ds, r 1: η i (s) = (d + 1), if i =, s [, 1 2(d+1) ) [ 2d+1 2(d+1), 1), (d + 1)Z i, if i 1, Λ = 1, s [ 2d 2i+1 2(d+1), 2d 2i+3 2(d+1) ), (d + 1)Z i, if i 1, Λ = 1, s [ 2i 1 2(d+1), 2i+1 2(d+1) ),, otherwise ω : degree 5 formula (ω (s) s) Theorem Let f C 6 pol (RN ), and ω is defined by the above Then P T f (x) QT n const /nf (x) n 2 21 / 37
22 Approximations of SDEs without jumps Operator Splitting method as a Cubature Remark : Cubature formula for Gaussian measure Cubature for Gaussian measure with degree m: A discrete-valued random variable Z such that E[Z α ] = x α 1 R d ( 2π) d e x for α m 2 /2 dx N-V + Gaussian cubature with degree 5 Cubature on Wiener space with degree 5, and with finite number of paths (but ω (s) s) The number of paths is "2 (the number of points of Gaussian cubature)" 22 / 37
23 Approximations of SDEs with jumps I (finite jump intesity) Approximations of SDEs with jumps I (finite jump intesity) 1 Introduction 2 Approximations of SDEs without jumps Approximation of SDEs by random ODEs (Extended) Cubature on Wiener space Operator Splitting method as a Cubature 3 Approximations of SDEs with jumps I (finite jump intesity) Operator Splitting Jump-adapted approximation 4 Approximations of SDEs with jumps II (general Lévy process) Overview Small jump approximation How to deal with compound Poisson term 5 Further Research 23 / 37
24 Approximations of SDEs with jumps I (finite jump intesity) Setting: jump-type SDE Let J(t) be a compound Poisson process Consider the following SDE with jumps X x t = x + d i= t t V i (Xs x ) dw i (s) + h(xs )dj(s) x Question: How can we construct a higher order scheme? Point 1 : W and J are independent Point 2 : Number of jumps = finite 24 / 37
25 Approximations of SDEs with jumps I (finite jump intesity) Operator Splitting Approach (OS): Operator Splitting 1 Consider the equation t dxd+1,t x = x + h(xd+1,s )dj(s) x easy to simulate (if we can easily simulate jump size) 2 Generator: L d+1 f (x) := (f (x + h(x)y) f (x))ν(dy) R d where ν is a finite Lévy measure 25 / 37
26 Higher order weak approximations of stochastic differential equations with and without jumps Approximations of SDEs with jumps I (finite jump intesity) Operator Splitting Approach (OS): Operator Splitting (continued) 3 As in the continuous diffusion case, we can consider approximations of e t d+1 i= L i by {e tl,, e tl d, e tl d+1 } 4 We can construct a second order scheme via 1 e tl 2 etl d etl d+1 e tl = e t d+1 i= L i + O(t 3 ) (This is first considered by Fujiwara (26, Master thesis)) 26 / 37
27 Approximations of SDEs with jumps I (finite jump intesity) Jump-adapted approximation Approach (JSAS): Jump (Size) Adapted Simulation Ref : Mordecki & Szepessy & Tempone & Zouraris (28) 1 First, simulate J(t): Jump time: < τ 1 < < τ k < T 2 Since τ j+1 τ j may be large, we use another partition { t j } so that {t i } := { t j } {τ j } satisfies t j+1 t j 3 (Jump adapted simulation) Continuous term approximation: X tj = Euler or cubature scheme in time interval [t j 1, t j ) for continuous part X t j 1,x t j Jump simulation: = x + d tj i= t j 1 V i (X t j 1,x s ) dw i (s) { X tj + h(x tj ) (jump size), if t j : jump time, X tj = X tj, otherwise 27 / 37
28 Approximations of SDEs with jumps I (finite jump intesity) Jump-adapted approximation Practical Problems Which is better, (OS) or (JSAS)? The best (adapted) choice of partition t 1 < < t n How to simulate by QMC (We need a restriction of the number of jumps) 28 / 37
29 Approximations of SDEs with jumps II (general Lévy process) Approximations of SDEs with jumps II (general Lévy process) 1 Introduction 2 Approximations of SDEs without jumps Approximation of SDEs by random ODEs (Extended) Cubature on Wiener space Operator Splitting method as a Cubature 3 Approximations of SDEs with jumps I (finite jump intesity) Operator Splitting Jump-adapted approximation 4 Approximations of SDEs with jumps II (general Lévy process) Overview Small jump approximation How to deal with compound Poisson term 5 Further Research 29 / 37
30 Approximations of SDEs with jumps II (general Lévy process) Overview Setting : Infinite Activity Lévy-driven SDE Genaral Lévy-driven SDE: X x t = x + d i= t t V i (Xs x ) dw i (s) + h(xs )dz x (s) (3) where Z (t) is a Lévy process without Brownian term, and Z t p>1 L p Generator: L = L + + L d + L d+1, L d+1 f (x) := f (x)h(x)b + (f (x + h(x)y) f (x) f (x)h(x)1 y 1 )ν(dy), R d where b R d, ν is a Lévy measure (Note: y 1 y 2 ν(dy) < ) 3 / 37
31 Approximations of SDEs with jumps II (general Lévy process) Overview Basic Strategy Extend the schemes (OS) & (JSAS) to the general case (3) 1 (OS): Tanaka & Kohatsu(29) Kohatsu & Ngo(211, submitted) 2 (JSAS): Kohatsu & Tankov(29) without Brownian term ie V i, 1 i d Approach (OS): Decompose L d+1 = L (1) d+1 + L(2) d+1 + L(3) d+1 where L (1) d+1 f (x) = b ϵ(hf )(x), L (2) d+1 f (x) = (f (x + h(x)y) f (x) f (x)h(x)y)ν(dy), L (3) d+1 f (x) = (f (x + h(x)y) f (x))ν(dy), y ϵ y >ϵ where ϵ < 1, b ϵ := b ϵ< y 1 yν(dy) 31 / 37
32 Approximations of SDEs with jumps II (general Lévy process) Small jump approximation Asmussen-Rosinski approximation for L (2) d+1 L (1) d+1 and L(3) d+1 respectively correspond to "drift" & "compound Poisson" term ( ) What is L (2) d+1? : Define Σ ϵ := y ϵ y iy j ν(dy) 1 i,j d ( ) L (2) d+1 = h(x)σ ϵh (x)d 2 f (x) + O y 3 ν(dy) y ϵ So we replace L (2) d+1 by a new "small diffusion"(itô form) term L (2) d+1 := h(x)σ ϵh (x)d 2 f (x) as an approximation If necessary for simulation, we can modify (L (1) (2) (2) d+1, L d+1 ) so that L d+1 becomes the generator of the SDE of Stratonovich form Control (ϵ, n) as y ϵ y 3 ν(dy) O(n k ) 32 / 37
33 Higher order weak approximations of stochastic differential equations with and without jumps Approximations of SDEs with jumps II (general Lévy process) How to deal with compound Poisson term The case y ν(dy) < + y 1 1 The number of jumps wrt L (3) d+1 goes to +, as ϵ To avoid this, we need some restriction for jumps 2 If we assume that y ν(dy) < +, then y 1 Note: b ϵ < + So we can use (OS) approximations for (L + L (1) d+1 ) : (drift) (2) L 1,, L d, L d+1 : (2d-diffusion) L (3) d+1 : (Jumps) 3 For L (3) d+1, we can construct an approximation process which has single or double jumps in (fixed) small time interval See Tanaka & Kohatsu, or Kohatsu & Ngo 33 / 37
34 Higher order weak approximations of stochastic differential equations with and without jumps Approximations of SDEs with jumps II (general Lévy process) How to deal with compound Poisson term The case y ν(dy) = + y 1 1 Kohatsu & Ngo discuss a case study where Z is a subordinated Brownian motion They construct an algorithm under the case y ν(dy) = +, but yν(dy) < + y 1 y 1 2 If we assume that y ν(dy) = +, then in general, y 1 Note: b ϵ = + We should not decompose (L (1) d+1 + L(3) d+1 ), since L(1) d+1 includes the truncation function for L (3) d+1 Consider (OS) approximations for (2) L 1,, L d, L d+1 : (2d-diffusion) L + (L (1) d+1 + L(3) d+1 ) : (ODE with Jumps) What can we do for (L (1) d+1 + L(3) d+1 )? 34 / 37
35 Further Research Further Research 1 Introduction 2 Approximations of SDEs without jumps Approximation of SDEs by random ODEs (Extended) Cubature on Wiener space Operator Splitting method as a Cubature 3 Approximations of SDEs with jumps I (finite jump intesity) Operator Splitting Jump-adapted approximation 4 Approximations of SDEs with jumps II (general Lévy process) Overview Small jump approximation How to deal with compound Poisson term 5 Further Research 35 / 37
36 Further Research Problems in time discretization 1 Find higher-order methods with degree more than 7 (d 2) (We have known only "existence"!) 2 High order approximation of Lévy driven SDEs How to deal with small jumps (which occurs infinitely many times) How to simulate Lévy processes whose Lévy measure has high singularity ( y ν(dy) = + ) y 1 3 Problems in computing conditional expectations E[ F t ]: Pricing American/Bermudan options Simulating forward-backward SDEs (FBSDEs) Recombination techniques: eg Chevance(1997), Lyons-Litterer(forthcoming, AAP), Tanaka(211, submitted) 36 / 37
37 Further Research Thank you for your attention 37 / 37
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 informationWeak Approximations for SDE s Driven by Lévy Processes
Weak Approximations for SDE s Driven by Lévy Processes Arturo Kohatsu-Higa and Hoang-Long Ngo Abstract. In this article we briefly survey recent advances in some simulation methods for Lévy driven stochastic
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 informationOptimal simulation schemes for Lévy driven stochastic differential equations arxiv: v1 [math.pr] 22 Apr 2012
Optimal simulation schemes for Lévy driven stochastic differential equations arxiv:124.4877v1 [math.p] 22 Apr 212 Arturo Kohatsu-Higa Salvador Ortiz-Latorre Peter Tankov Abstract We consider a general
More information5 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 informationOptimal simulation schemes for Lévy driven stochastic differential equations
Optimal simulation schemes for Lévy driven stochastic differential equations Arturo Kohatsu-Higa Salvador Ortiz-Latorre Peter Tankov Abstract We consider a general class of high order weak approximation
More informationPROBABILITY: 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 informationControlled 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 informationPROBABILITY: 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 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 informationTHE GEOMETRY OF ITERATED STRATONOVICH INTEGRALS
THE GEOMETRY OF ITERATED STRATONOVICH INTEGRALS CHRISTIAN BAYER Abstract. We give a summary on the geometry of iterated Stratonovich integrals. For this exposition, we always have the connection to stochastic
More informationAn 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 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 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 informationSimulation methods for stochastic models in chemistry
Simulation methods for stochastic models in chemistry David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison SIAM: Barcelona June 4th, 21 Overview 1. Notation
More informationAccurate 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 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 informationWORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS. Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid
WORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid 1 I. OVERVIEW 2 The importance of splitting integrators
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 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 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 informationStochastic optimal control with rough paths
Stochastic optimal control with rough paths Paul Gassiat TU Berlin Stochastic processes and their statistics in Finance, Okinawa, October 28, 2013 Joint work with Joscha Diehl and Peter Friz Introduction
More informationEvaluation of the HJM equation by cubature methods for SPDE
Evaluation of the HJM equation by cubature methods for SPDEs TU Wien, Institute for mathematical methods in Economics Kyoto, September 2008 Motivation Arbitrage-free simulation of non-gaussian bond markets
More informationThe Hopf algebraic structure of stochastic expansions and efficient simulation
The Hopf algebraic structure of stochastic expansions and efficient simulation Kurusch Ebrahimi Fard, Alexander Lundervold, Simon J.A. Malham, Hans Munthe Kaas and Anke Wiese York: 15th October 2012 KEF,
More informationStochastic Volatility and Correction to the Heat Equation
Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century
More informationA 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 informationHardy-Stein identity and Square functions
Hardy-Stein identity and Square functions Daesung Kim (joint work with Rodrigo Bañuelos) Department of Mathematics Purdue University March 28, 217 Daesung Kim (Purdue) Hardy-Stein identity UIUC 217 1 /
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationc 2002 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 4 No. pp. 507 5 c 00 Society for Industrial and Applied Mathematics WEAK SECOND ORDER CONDITIONS FOR STOCHASTIC RUNGE KUTTA METHODS A. TOCINO AND J. VIGO-AGUIAR Abstract. A general
More informationNumerical Methods with Lévy Processes
Numerical Methods with Lévy Processes 1 Objective: i) Find models of asset returns, etc ii) Get numbers out of them. Why? VaR and risk management Valuing and hedging derivatives Why not? Usual assumption:
More informationOn the higher-order weak approximationof SDEs
On the higher-order weak approximationof SDEs Syoiti Ninomiya Tokyo Institute of Technology The 3rd Workshop on numerical methods for solving the filtering problem and high order methods for solving parabolic
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 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 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 informationEstimates for the density of functionals of SDE s with irregular drift
Estimates for the density of functionals of SDE s with irregular drift Arturo KOHATSU-HIGA a, Azmi MAKHLOUF a, a Ritsumeikan University and Japan Science and Technology Agency, Japan Abstract We obtain
More informationLecture 6: Bayesian Inference in SDE Models
Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs
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 informationThe Wiener Itô Chaos Expansion
1 The Wiener Itô Chaos Expansion The celebrated Wiener Itô chaos expansion is fundamental in stochastic analysis. In particular, it plays a crucial role in the Malliavin calculus as it is presented in
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 informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
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 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 informationMulti-Factor Lévy Models I: Symmetric alpha-stable (SαS) Lévy Processes
Multi-Factor Lévy Models I: Symmetric alpha-stable (SαS) Lévy Processes Anatoliy Swishchuk Department of Mathematics and Statistics University of Calgary Calgary, Alberta, Canada Lunch at the Lab Talk
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 informationOn Optimal Stopping Problems with Power Function of Lévy Processes
On Optimal Stopping Problems with Power Function of Lévy Processes Budhi Arta Surya Department of Mathematics University of Utrecht 31 August 2006 This talk is based on the joint paper with A.E. Kyprianou:
More informationGeometric 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 informationWEAK VERSIONS OF STOCHASTIC ADAMS-BASHFORTH AND SEMI-IMPLICIT LEAPFROG SCHEMES FOR SDES. 1. Introduction
WEAK VERSIONS OF STOCHASTIC ADAMS-BASHFORTH AND SEMI-IMPLICIT LEAPFROG SCHEMES FOR SDES BRIAN D. EWALD 1 Abstract. We consider the weak analogues of certain strong stochastic numerical schemes considered
More informationFast-slow systems with chaotic noise
Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 12, 215 Averaging and homogenization workshop, Luminy. Fast-slow systems
More informationSome 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 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 informationLAN property for ergodic jump-diffusion processes with discrete observations
LAN property for ergodic jump-diffusion processes with discrete observations Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Arturo Kohatsu-Higa (Ritsumeikan University, Japan) &
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 informationPoisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation
Poisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation Jingyi Zhu Department of Mathematics University of Utah zhu@math.utah.edu Collaborator: Marco Avellaneda (Courant
More informationSDE Coefficients. March 4, 2008
SDE Coefficients March 4, 2008 The following is a summary of GARD sections 3.3 and 6., mainly as an overview of the two main approaches to creating a SDE model. Stochastic Differential Equations (SDE)
More informationFast-slow systems with chaotic noise
Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 1, 216 Statistical properties of dynamical systems, ESI Vienna. David
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 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 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 informationLecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations
Lecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations Department of Biomedical Engineering and Computational Science Aalto University April 28, 2010 Contents 1 Multiple Model
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 information[Ahmed*, 5(3): March, 2016] ISSN: (I2OR), Publication Impact Factor: 3.785
IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY DENSITIES OF DISTRIBUTIONS OF SOLUTIONS TO DELAY STOCHASTIC DIFFERENTIAL EQUATIONS WITH DISCONTINUOUS INITIAL DATA ( PART II)
More informationErnesto Mordecki 1. Lecture III. PASI - Guanajuato - June 2010
Optimal stopping for Hunt and Lévy processes Ernesto Mordecki 1 Lecture III. PASI - Guanajuato - June 2010 1Joint work with Paavo Salminen (Åbo, Finland) 1 Plan of the talk 1. Motivation: from Finance
More informationJump-type Levy Processes
Jump-type Levy Processes Ernst Eberlein Handbook of Financial Time Series Outline Table of contents Probabilistic Structure of Levy Processes Levy process Levy-Ito decomposition Jump part Probabilistic
More informationStrong Predictor-Corrector Euler Methods for Stochastic Differential Equations
Strong Predictor-Corrector Euler Methods for Stochastic Differential Equations Nicola Bruti-Liberati 1 and Eckhard Platen July 14, 8 Dedicated to the 7th Birthday of Ludwig Arnold. Abstract. This paper
More informationIntroduction to Computational Stochastic Differential Equations
Introduction to Computational Stochastic Differential Equations Gabriel J. Lord Catherine E. Powell Tony Shardlow Preface Techniques for solving many of the differential equations traditionally used by
More informationFE610 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 informationADAPTIVE WEAK APPROXIMATION OF REFLECTED AND STOPPED DIFFUSIONS
ADAPTIVE WEAK APPROXIMATION OF REFLECTED AND STOPPED DIFFUSIONS CHRISTIAN BAYER, ANDERS SZEPESSY, AND RAÚL TEMPONE ABSTRACT. We study the weak approximation problem of diffusions, which are reflected at
More informationSTOCHASTIC PROCESSES FOR PHYSICISTS. Understanding Noisy Systems
STOCHASTIC PROCESSES FOR PHYSICISTS Understanding Noisy Systems Stochastic processes are an essential part of numerous branches of physics, as well as biology, chemistry, and finance. This textbook provides
More informationFiltrations, Markov Processes and Martingales. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition
Filtrations, Markov Processes and Martingales Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition David pplebaum Probability and Statistics Department,
More informationSimulation Method for Solving Stochastic Differential Equations with Constant Diffusion Coefficients
Journal of mathematics and computer Science 8 (2014) 28-32 Simulation Method for Solving Stochastic Differential Equations with Constant Diffusion Coefficients Behrouz Fathi Vajargah Department of statistics,
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 informationPrecise Asymptotics for Infinite Dimensional Itô-Lyons Maps of Brownian Rough Paths
Precise Asymptotics for Infinite Dimensional Itô-Lyons Maps of Brownian Rough Paths Jointwork with Yuzuru INAHAMA (Nagoya Univ.) Hiroshi KAWABI (Department of Mathematics, Faculty of Science Okayama University,
More informationJump-adapted discretization schemes for Lévy-driven SDEs
Jump-adapted discretization schemes for Lévy-driven SDEs Arturo Kohatsu-Higa Osaka University and Japan Science and Technology Agency Graduate School of Engineering Sciences Machikaneyama cho -3, Osaka
More informationDivergence theorems in path space II: degenerate diffusions
Divergence theorems in path space II: degenerate diffusions Denis Bell 1 Department of Mathematics, University of North Florida, 4567 St. Johns Bluff Road South, Jacksonville, FL 32224, U. S. A. email:
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 informationTyler 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 informationQuasi-invariant Measures on Path Space. Denis Bell University of North Florida
Quasi-invariant Measures on Path Space Denis Bell University of North Florida Transformation of measure under the flow of a vector field Let E be a vector space (or a manifold), equipped with a finite
More informationCUBATURE ON WIENER SPACE: PATHWISE CONVERGENCE
CUBATUR ON WINR SPAC: PATHWIS CONVRGNC CHRISTIAN BAYR AND PTR K. FRIZ Abstract. Cubature on Wiener space [Lyons, T.; Victoir, N.; Proc. R. Soc. Lond. A 8 January 2004 vol. 460 no. 2041 169-198] provides
More informationHairer /Gubinelli-Imkeller-Perkowski
Hairer /Gubinelli-Imkeller-Perkowski Φ 4 3 I The 3D dynamic Φ 4 -model driven by space-time white noise Let us study the following real-valued stochastic PDE on (0, ) T 3, where ξ is the space-time white
More informationWeak Convergence of Numerical Methods for Dynamical Systems and Optimal Control, and a relation with Large Deviations for Stochastic Equations
Weak Convergence of Numerical Methods for Dynamical Systems and, and a relation with Large Deviations for Stochastic Equations Mattias Sandberg KTH CSC 2010-10-21 Outline The error representation for weak
More informationNOTES III SAN DIEGO 4
SAN DIEGO 1 NOTES I The main original paper http://dmle.cindoc.csic.es/pdf/ MATEMATICAIBEROAMERICANA_1998_14_02_01.pdf is not so bad, and contains a proof of the "neo-classical" factorial estimate and
More informationOn rough PDEs. Samy Tindel. University of Nancy. Rough Paths Analysis and Related Topics - Nagoya 2012
On rough PDEs Samy Tindel University of Nancy Rough Paths Analysis and Related Topics - Nagoya 2012 Joint work with: Aurélien Deya and Massimiliano Gubinelli Samy T. (Nancy) On rough PDEs Nagoya 2012 1
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 informationRandom attractors and the preservation of synchronization in the presence of noise
Random attractors and the preservation of synchronization in the presence of noise Peter Kloeden Institut für Mathematik Johann Wolfgang Goethe Universität Frankfurt am Main Deterministic case Consider
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Centre for Computational Statistics and Machine Learning University College London c.archambeau@cs.ucl.ac.uk CSML
More informationBrownian 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 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 informationImprecise Filtering for Spacecraft Navigation
Imprecise Filtering for Spacecraft Navigation Tathagata Basu Cristian Greco Thomas Krak Durham University Strathclyde University Ghent University Filtering for Spacecraft Navigation The General Problem
More informationSession 1: Probability and Markov chains
Session 1: Probability and Markov chains 1. Probability distributions and densities. 2. Relevant distributions. 3. Change of variable. 4. Stochastic processes. 5. The Markov property. 6. Markov finite
More informationStochastic 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 informationStatistics 992 Continuous-time Markov Chains Spring 2004
Summary Continuous-time finite-state-space Markov chains are stochastic processes that are widely used to model the process of nucleotide substitution. This chapter aims to present much of the mathematics
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 informationCIMPA SCHOOL, 2007 Jump Processes and Applications to Finance Monique Jeanblanc
CIMPA SCHOOL, 27 Jump Processes and Applications to Finance Monique Jeanblanc 1 Jump Processes I. Poisson Processes II. Lévy Processes III. Jump-Diffusion Processes IV. Point Processes 2 I. Poisson Processes
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 informationWhat s more chaotic than chaos itself? Brownian Motion - before, after, and beyond.
Include Only If Paper Has a Subtitle Department of Mathematics and Statistics What s more chaotic than chaos itself? Brownian Motion - before, after, and beyond. Math Graduate Seminar March 2, 2011 Outline
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
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 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 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 informationResearch Article A Necessary Characteristic Equation of Diffusion Processes Having Gaussian Marginals
Abstract and Applied Analysis Volume 01, Article ID 598590, 9 pages doi:10.1155/01/598590 Research Article A Necessary Characteristic Equation of Diffusion Processes Having Gaussian Marginals Syeda Rabab
More information