Qualitative behaviour of numerical methods for SDEs and application to homogenization
|
|
- Beverly Clarke
- 6 years ago
- Views:
Transcription
1 Qualitative behaviour of numerical methods for SDEs and application to homogenization K. C. Zygalakis Oxford Centre For Collaborative Applied Mathematics, University of Oxford. Center for Nonlinear Analysis, Carnegie Mellon University, 20/10/2011 K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 1 / 45
2 Outline 1 Modified Equations ODE theory. Main idea for SDEs. Different numerical methods and Associated Modified Equations. Numerical examples. 2 Application to Homogenization Long time behaviour and homogenization. Numerical algorithms/results. From homogenization to averaging in cellular flows. 3 Higher order numerical methods based on modified equations. Key idea. One simple example. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 2 / 45
3 Motivating example Introduction K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 3 / 45
4 Introduction Interesting Question The two numerical methods have the same order of convergence but completely different qualitative behaviour. Is there a way to distinguish between these two methods? A very powerful tool for addressing this question is backward error analysis (modified equations). K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 4 / 45
5 Introduction Modified equations for ODEs Ordinary Differential Equations dx dt = f (x), and let x n be a numerical approximation of x of order p: x(nh) x n = O(h p ). Can I find X (t) satisfying another ODE (modified equation) such that: X (nh) x n = O(h p+q ). K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 5 / 45
6 Introduction Euler method-one dimension Ordinary Differential Equations x n+1 = x n + hf (x n ). Modified equation: since dx dt = f (X ) h 2 f (X )f (X ), X (nh) x n = O(h 2 ). K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 6 / 45
7 Introduction Ordinary Differential Equations Sketch proof dx dt = f (X ) + hg(x ). X (h) = X (0) + h 0 (f (X (s)) + hg(x (s))) ds = X (0) + hf (X (0)) + h 2 g(x (0)) + h2 2 f (X (0))f (X (0)) + O(h 3 ). Assume x 0 = X (0) then X (h) x 1 = h (g(x 2 (0)) + 1 ) 2 f (X (0))f (X (0)) + O(h 3 ), and thus g(x) = 1 2 f (x)f (x). K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 7 / 45
8 Introduction Stochastic Differential Equations Stochastic Differential Equations and Numerical Methods dx = u(x)dt + σ(x)dw t, Euler method: θ-milstein method: x n+1 = x n + hu(x n ) + hσ(x n )ξ n, x n+1 = x n +θhu(x n+1 )+(1 θ)hu(x n )+ hσ(x n )ξ n + h 2 σ(x n)σ (1) (x n )(ξ 2 n 1), where ξ n N (0, 1). K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 8 / 45
9 Introduction Stochastic Differential Equations Weak and Strong Convergence Weak convergence: We look at E(φ(x(nh))) E(φ(x n )). Strong convergence: We look at E x(nh) x n. In general the weak and strong order of convergence of a numerical method NEEDS NOT to be the same!!! K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 9 / 45
10 Introduction Statement of the Problem Stochastic Differential Equations Let x(t) satisfy the following SDE: dx = u 1 (x)dt + σ 1 (x)dw t, and x n be its numerical approximation at T = nh by a weak p-order method i.e E(φ(x(T ))) E(φ(x n )) = O(h p ), φ C. We want to develop a procedure that allows us to evaluate the properties of our weak numerical scheme. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 10 / 45
11 Introduction Statement of the Problem Stochastic Differential Equations Let x(t) satisfy the following SDE: dx = u 1 (x)dt + σ 1 (x)dw t, and x n be its numerical approximation at T = nh by a weak p-order method i.e E(φ(x(T ))) E(φ(x n )) = O(h p ), φ C. We want to develop a procedure that allows us to evaluate the properties of our weak numerical scheme. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 10 / 45
12 First Modified Equation Introduction Stochastic Differential Equations We want to find a modified SDE of the form (i.e., find v 2 and σ 2 ) d x = [u 1 ( x) + hu 2 ( x)] + [σ 1 ( x) + hσ 2 ( x)] dw t, for which E(φ( x(t ))) E(φ(x n )) = O(h p+1 ), φ C. For the rest of the talk we concentrate in the case where p = 1. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 11 / 45
13 Main Idea Generators for ODEs and SDEs ODE: dx = h(x)dt, Lu := h(x) x u. SDE: dx = h(x)dt + σ(x)dw t, Lu := h(x) x u σ(x)σt (x) : x x u. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 12 / 45
14 Main Idea Backward Kolmogorov Equation u t = Lu, u(x, 0) = φ(x). Then u(x, t) = E(φ(x(t)) x(0) = x). K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 13 / 45
15 Stochastic B-series Main Idea By integrating over time the backward Kolmogorov Equation and taking a Taylor expansion of u(x, s) around s = 0, we obtain, (assuming appropriate smoothness of the drift and diffusion term) u(x, h) φ(x) = k=0 h k+1 (k + 1)! Lk+1 φ(x). Note that in the case where φ(x) = x, σ(x) = 0, this expansion correspond to the B-series expansion of the ODE dx = v 1 (x)dt. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 14 / 45
16 Main Idea Local Error/Global Error A weak first order numerical method has the following expansion u num (x, h) φ(x) = hlφ(x) + h 2 L e φ(x) + O(h 3 ), and so ( ) 1 u(x, h) u num (x, h) = h 2 2 L2 φ(x) L e φ(x), Local Error which implies that u(x, T ) u num (x, T ) = O(h). Global Error K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 15 / 45
17 Main Idea Generator of the Modified Equation Remember that the 1-st modified equation is of the form d x = [u 1 ( x) + hu 2 ( x)] + [σ 1 ( x) + hσ 2 ( x)] dw t. Its generator L can be written as L = L 0 + hl 1 + h 2 L 2, where L 0 is the generator of the original SDE and L 1 φ := u 2 (x) dφ dx + σ 1(x)σ 2 (x) d 2 φ dx 2. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 16 / 45
18 Main Equation Main Idea If we now subtract the Taylor expansion of the numerical method from the stochastic B-series of the modified equation we see that in order for the local error to be O( t 3 ) we need L 1 φ = L e φ 1 2 L2 0φ, φ C. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 17 / 45
19 Different Numerical Methods Euler-Maryama Method In the case of Euler-Maryama method in the case of multiplicative noise it turns out that a modified equation does not exist since L 1 φ + σ3 1 (x) σ (1) 1 2 (x)φ(3) (x). as L 1 is a second order partial differential operator!!! K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 18 / 45
20 θ-milstein Method Different Numerical Methods u 2 (x) = σ 2 (x) = ( θ 1 ) ( v 1 (x)v (1) 2 1 (x) + σ2 1 (x) 2 ( θ 1 ) σ 1 (x)v (1) 1 2 (x) 1 2 v 1(x)σ (1) v (2) 1 (x) ), 1 (x) σ2 1 (x) 4 σ (2) 1 (x). K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 19 / 45
21 Numerical examples Geometric Brownian motion SDEs Driven by Multiplicative Noise dx = µxdt + σxdw t, d X = [(µ h2 ) ] µ2 X dt + σ X (1 hµ) dw t error error original SDE modified SDE 10 3 original SDE modified SDE t First moment t Second moment K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 20 / 45
22 Numerical examples Modified equations Linear SDEs with additive noise dx = Axdt + ΣdW t, Numerical Approximation: x(h) = A(h)x + f (h, ω). Example (Euler-Maryama): A(h) = (I + ha), f (h, ω) = Σ hξ. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 21 / 45
23 Numerical examples Modified equations Modified Equation and its coefficients dx = Ãxdt + ΣdW t, Ã = log(a(h)), h eãh Σ Σ T eãt h Σ Σ T = ÃJ + JÃ T, where J = E(ff T ). K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 22 / 45
24 Numerical examples Orstein Uhlenbeck Process Modified equations Forward Euler: Backward Euler: dx = γxdt + σdw t. log(1 γh) Ã =, h 2 log(1 γh) Σ = σ (1 γh) 2 1. log(1 + γh) Ã =, h 2 log(1 + γh) Σ = σ 1 (1 + γh) 2. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 23 / 45
25 Invariant Measure Numerical examples Modified equations lim t E(x2 (t)) = lim t E(x2 (t)) = σ 2 2γ γ 2, Forward Euler h σ 2 (1 + γh) 2γ + γ 2, Backward Euler. h numerical solution modified equation numerical solution modified equation Var(x) 0.58 Var(x) t Forward Euler t Backward Euler Figure: lim t E(x 2 (t)) as a function of h. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 24 / 45
26 An application to homogenization Passive Tracers, Effective Diffusivity Long Time Behaviour and Homogenization dx = v(x)dt + σdw t, where v(x) is a periodic function. It is possible to show using homogenization that E(x(t) x(t)) lim = K. t 2t We will refer to K as the effective diffusivity matrix K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 25 / 45
27 An application to homogenization Long Time Behaviour and Homogenization Velocity Field of Interest Example We are interested in the following 2-dimensional incompressible velocity field v(x) = Ψ(x), where Ψ(x) = sin x 1 sin x 2 Result In this case it is known that K = DI 2, where D R depending on σ for passive tracers and that the following result is true for passive tracers D(σ) σ, σ 1 K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 26 / 45
28 An application to homogenization Key Property of the Velocity Field Long Time Behaviour and Homogenization Our velocity field v(x) can be written as v(x) = ( ) ( 1/2 1/2 sin(x +1/2 1 + x 2 ) + 1/2 = 2 ( d j v j ej, x ), j=1 where e j, d j R 2 with the property e j, d j = 0. ) sin(x 1 x 2 ), This is a key property for the construction of our method which is a stochastic extension of a splitting method proposed by Quispel K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 27 / 45
29 An application to homogenization Numerical algorithm/results Description of the Method: The method in the case of passive tracers involves these 3 steps: ( Step 1: Solve ẋ 1 = d 1 v 1 e1, x 1 ), ( Step 2: Solve ẋ 2 = d 2 v 2 e2, x 2 ), Step 3: Solve ẋ 3 = σ β 1. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 28 / 45
30 An application to homogenization The Deterministic Case Numerical algorithm/results Splitting method for σ = 0. Euler method for σ = 0. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 29 / 45
31 The Case σ 1 An application to homogenization Numerical algorithm/results Splitting method for σ = Euler method for σ = K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 30 / 45
32 An application to homogenization Calculating Effective Diffusivities Numerical algorithm/results Volume preserving Euler!Maruyama 10!1 10! !3 K K 10!1 10!4 10!2 10!5 10!3 10!2 10! Comparison of the two methods.! 10!6 10!3 10!2 10!1 Splitting Method.! K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 31 / 45
33 Mean Hamiltonian An application to homogenization Numerical algorithm/results We apply Itô s formula to H = Ψ we obtain dψ dt = σ2 Ψ + M.T which implies that the mean Hamiltonian decays like e σ2 t K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 32 / 45
34 An application to homogenization Numerical algorithm/results Numerical calculation of the mean Hamiltonian with the two methods Euler Method Stochastic Splitting Method Theory Mean Hamiltonian ! t x 10 4 Figure: Mean value of the Hamiltonian as a function of time, for t = 10 1, σ = K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 33 / 45
35 An application to homogenization Numerical algorithm/results Modified Equations for the Euler Method dx = + σ ( v(x) t 2 ( v(x))v(x) σ2 t 4 ) dw t. ( 1 t 2 v T (x) ) v(x) dt dψ dt = t 2 (cos2 x 1 + cos 2 x 2 )Ψ σ 2 Ψ(1 + t cos x 1 cos x 2 ) + σ2 2 t (cos 2 x 1 cos 2 x 2 Ψ Ψ 3 ) + M t. 4 K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 34 / 45
36 An application to homogenization From Homogenization to averaging in cellular flows. Statement of the problem u(x) = Ψ(x), Ψ(x) = 1 π sin πx 1 sin πx 2. X t satisfies the following SDE dx t = Au(X t )dt + 2dW t. Exit time problem τ + Au τ = 1 in D, τ = 0 on D, where D = [ L/2, L/2] [ L/2, L/2]. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 35 / 45
37 An application to homogenization Known results and open questions From Homogenization to averaging in cellular flows. 1 A fixed, L, homogenization (τ ). 2 L fixed, A, averaging (τ 0). 3 If D = B L a disk of radius L, then τ(x) L 2 x 2. 4 A = L α and L α < 4, homogenization (τ L 2 α/2 ). α > 4, averaging (τ L 2 α/2 ). Question How about α = 4? K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 36 / 45
38 An application to homogenization Numerical investigations I From Homogenization to averaging in cellular flows. 2 Define Y t = X At, then dy t = u(y t )dt + A dw t and τ(y) = Aτ(x) Asymptotic behaviour of exit time (α = 1) Exit problem from a disk (L = 40) K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 37 / 45
39 An application to homogenization Numerical investigations II, Case α = 4 From Homogenization to averaging in cellular flows. Asymptotic behaviour of exit time (α = 4) Comparison of c.d.f K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 38 / 45
40 Higher order methods Key idea i) Choose a numerical method for the original SDE. ii) Write down a suitably chosen SDE different than the original one (this SDE depends on the choice of the method from step i). iii) Apply the numerical method from step i to the SDE from step ii. Example dx t = u 1 (X t)dt + σ 1 (X t) i) x n+1 = x n +θhu 1 (x n+1 )+(1 θ)hu 1 (x n)+ hσ 1 (x n)ξ n + h 2 σ 1(x n)σ (1) 1 (xn)(ξ2 n 1) ii) dx t = ũ(x t)dt + σ(x t), ũ = u 1 hu 2, σ = σ 1 hσ 2 iii) x n+1 = x n + θhũ(x n+1 ) + (1 θ)hũ(x n) + h σ(x n)ξ n + h 2 σ 1(x n)σ (1) 1 (xn)(ξ2 n 1) K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 39 / 45
41 Higher order methods Application to an economy model for asset prices dx 1 = β 1 X 1 X 2 dw 1, dx 2 = (X 2 X 3 )dt + β 2 X 2 dw 2, dx 3 = α(x 2 X 3 )dt, N. Hofmann, E. Platen, M. Schweizer. Option pricing under incompletness and stochastic volatility. Mathematical Finance, 2(3): , (1992). K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 40 / 45
42 Higher order methods Numerical Investigations Error for E(X1 2). Nonstiff case α = 1. Error for E(X 1 2 ). Very stiff case α = 100. Error for E(X1 2). Stiff case α = 25. Error for E(X 2 2 ). Stiff case α = 25. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 41 / 45
43 Conclusions Summary 1 It is not always possible to write down a modified Itô SDE for a given numerical method. 2 In the case of linear SDEs with additive noise it is possible to write down an -modified equation that the numerical method satisfy exactly in the weak sense. 3 It is possible to generalize ideas from the backward error analysis of ODEs to SDEs. 4 Modified equations can be used as a tool for constructing higher order methods. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 42 / 45
44 Summary Future work 1 Find modified equations for numerical methods with respect to strong convergence. 2 Give a rigorous explanation for failing to find a modified SDE for the Euler method in case of multiplicative noise. 3 Use modified equations to characterize the invariant measure approximated by different numerical schemes. 4 Compare exit times from a square for different starting points, with the ones of the effective Brownian motion. 5 Study exit times in case where the inertia is important (inertial particles). K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 43 / 45
45 Acknowledgements Thank for your attention! Collaborators: A. Abdulle (EPFL), D. Cohen (Basel), G. Iyer (CMU), G. Pavliotis (Imperial), A. M. Stuart (Warwick), G. Villmart (ENS Cachan). Funding: David Crighton Fellowship and award KUK-C , made by King Abdullah University of Science and Technology (KAUST). K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 44 / 45
46 References Ernst Hairer, Christian Lubich, and Gerhard Wanner. Geometric numerical integration, volume 31 of Springer Series in Computational Mathematics. Springer-Verlag, Berlin, (2002). G.R.W. Quispel and D.I. McLaren, Explicit volume-preserving and symplectic integrators for trigonometric polynomial flows, J. Comp. Phys, , (2003). T. Shardlow. Modified equations for stochastic differential equations.bit, 46(1): , (2006). I. Lenane, K. Burrage and G. Lythe. Numerical methods for second order stochastic equations. SIAM J.Sci.Comp, 29(1): , (2008). G. A. Pavliotis, A. M. Stuart, K. C. Zygalakis. Calculating Effective Diffusivities in Vanishing Molecular Diffusion, J. Comp. Phys. 228, (2009). K. C. Zygalakis. On the existence and the applications of modified equations for stochastic differential equations. SIAM J. Sci. Comput.. 33, (2011). A. Abdulle, D. Cohen, G. Villmart, K. C. Zygalakis. High order weak methods for stochastic differential equations based on modified equations. Submitted to SIAM J. Sci. Comput.. G. Iyer, T. Komorowski, A. Novikov, L. Ryzhik. From homogenization to averaging in cellular flows. (Submitted) A. Abdulle, G. Villmart, K. C. Zygalakis. Explicit higher order methods for stiff stochastic differential equations. In preparation. K. C. Zygalakis (University of Oxford) Modified Equations for SDEs 45 / 45
Computing ergodic limits for SDEs
Computing ergodic limits for SDEs M.V. Tretyakov School of Mathematical Sciences, University of Nottingham, UK Talk at the workshop Stochastic numerical algorithms, multiscale modelling and high-dimensional
More informationc 2012 Society for Industrial and Applied Mathematics
MULTISCALE MODEL. SIMUL. Vol., No. 3, pp. 46 58 c 22 Society for Industrial and Applied Mathematics NUMERICAL STUDIES OF HOMOGENIZATION UNDER A FAST CELLULAR FLOW GAUTAM IYER AND KONSTANTINOS C. ZYGALAKIS
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 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 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 informationBackward error analysis
Backward error analysis Brynjulf Owren July 28, 2015 Introduction. The main source for these notes is the monograph by Hairer, Lubich and Wanner [2]. Consider ODEs in R d of the form ẏ = f(y), y(0) = y
More informationExit times of diffusions with incompressible drift
Exit times of diffusions with incompressible drift Gautam Iyer, Carnegie Mellon University gautam@math.cmu.edu Collaborators: Alexei Novikov, Penn. State Lenya Ryzhik, Stanford University Andrej Zlatoš,
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 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 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 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 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 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 informationQuantization of stochastic processes with applications on Euler-Maruyama schemes
UPTEC F15 063 Examensarbete 30 hp Oktober 2015 Quantization of stochastic processes with applications on Euler-Maruyama schemes Viktor Edward Abstract Quantization of stochastic processes with applications
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationNUMERICAL STUDIES OF HOMOGENIZATION UNDER A FAST CELLULAR FLOW. (1.2) dx = Av(X) dt + dw t,
MULTISCLE MODEL. SIMUL. Vol., No. 3, pp. 58 c Society for Industrial and pplied Mathematics NUMERICS FOR HOMOGENIZTION IN FST CELLULR FLOW 7 NUMERICL STUDIES OF HOMOGENIZTION UNDER FST CELLULR FLOW GUTM
More informationExit times of diffusions with incompressible drifts
Exit times of diffusions with incompressible drifts Andrej Zlatoš University of Chicago Joint work with: Gautam Iyer (Carnegie Mellon University) Alexei Novikov (Pennylvania State University) Lenya Ryzhik
More informationMean-square A-stable diagonally drift-implicit integrators of weak second order for stiff Itô stochastic differential equations
Mean-square A-stable diagonally drift-implicit integrators of weak second order for stiff Itô stochastic differential equations Assyr Abdulle, Gilles Vilmart, Konstantinos Zygalakis To cite this version:
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 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 informationHandbook of Stochastic Methods
C. W. Gardiner Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences Third Edition With 30 Figures Springer Contents 1. A Historical Introduction 1 1.1 Motivation I 1.2 Some Historical
More informationMATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics
MATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics MATHICSE Technical Report Nr. 35.01 September 01 Mean-square A-stable diagonally
More informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
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 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 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 informationNumerical treatment of stochastic delay differential equations
Numerical treatment of stochastic delay differential equations Evelyn Buckwar Department of Mathematics Heriot-Watt University, Edinburgh Evelyn Buckwar Roma, 16.10.2007 Overview An example: machine chatter
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More informationNUMERICAL METHODS FOR SECOND-ORDER STOCHASTIC DIFFERENTIAL EQUATIONS
SIAM J SCI COMPUT Vol 9, No, pp 45 64 c 7 Society for Industrial and Applied Mathematics NUMERICAL METHODS FOR SECOND-ORDER STOCHASTIC DIFFERENTIAL EQUATIONS KEVIN BURRAGE, IAN LENANE, AND GRANT LYTHE
More informationNumerical solution of second-order stochastic differential equations with Gaussian random parameters
Journal of Linear and Topological Algebra Vol. 2, No. 4, 23, 223-235 Numerical solution of second-order stochastic differential equations with Gaussian random parameters R. Farnoosh a, H. Rezazadeh a,,
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 informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
More informationA stochastic particle system for the Burgers equation.
A stochastic particle system for the Burgers equation. Alexei Novikov Department of Mathematics Penn State University with Gautam Iyer (Carnegie Mellon) supported by NSF Burgers equation t u t + u x u
More informationSummer 2017 MATH Solution to Exercise 5
Summer 07 MATH00 Solution to Exercise 5. Find the partial derivatives of the following functions: (a (xy 5z/( + x, (b x/ x + y, (c arctan y/x, (d log((t + 3 + ts, (e sin(xy z 3, (f x α, x = (x,, x n. (a
More informationHandbook of Stochastic Methods
Springer Series in Synergetics 13 Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences von Crispin W Gardiner Neuausgabe Handbook of Stochastic Methods Gardiner schnell und portofrei
More informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
More informationConvergence of Particle Filtering Method for Nonlinear Estimation of Vortex Dynamics
Convergence of Particle Filtering Method for Nonlinear Estimation of Vortex Dynamics Meng Xu Department of Mathematics University of Wyoming February 20, 2010 Outline 1 Nonlinear Filtering Stochastic Vortex
More informationDynamical systems with Gaussian and Levy noise: analytical and stochastic approaches
Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Noise is often considered as some disturbing component of the system. In particular physical situations, noise becomes
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 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 informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
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 informationLYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR
LYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR Peter H. Baxendale Department of Mathematics University of Southern California Los Angeles, CA 90089-3 USA baxendal@math.usc.edu
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 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 informationOn the construction of discrete gradients
On the construction of discrete gradients Elizabeth L. Mansfield School of Mathematics, Statistics and Actuarial Science, University of Kent, Canterbury, CT2 7NF, U.K. G. Reinout W. Quispel Department
More informationIntroduction to asymptotic techniques for stochastic systems with multiple time-scales
Introduction to asymptotic techniques for stochastic systems with multiple time-scales Eric Vanden-Eijnden Courant Institute Motivating examples Consider the ODE {Ẋ = Y 3 + sin(πt) + cos( 2πt) X() = x
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 informationStochastic Gradient Descent in Continuous Time
Stochastic Gradient Descent in Continuous Time Justin Sirignano University of Illinois at Urbana Champaign with Konstantinos Spiliopoulos (Boston University) 1 / 27 We consider a diffusion X t X = R m
More informationRough Burgers-like equations with multiplicative noise
Rough Burgers-like equations with multiplicative noise Martin Hairer Hendrik Weber Mathematics Institute University of Warwick Bielefeld, 3.11.21 Burgers-like equation Aim: Existence/Uniqueness for du
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 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 informationHigher order weak approximations of stochastic differential equations with and without jumps
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
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
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 informationRandom and Deterministic perturbations of dynamical systems. Leonid Koralov
Random and Deterministic perturbations of dynamical systems Leonid Koralov - M. Freidlin, L. Koralov Metastability for Nonlinear Random Perturbations of Dynamical Systems, Stochastic Processes and Applications
More informationSolution of Stochastic Optimal Control Problems and Financial Applications
Journal of Mathematical Extension Vol. 11, No. 4, (2017), 27-44 ISSN: 1735-8299 URL: http://www.ijmex.com Solution of Stochastic Optimal Control Problems and Financial Applications 2 Mat B. Kafash 1 Faculty
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 informationNested 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 informationA new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009
A new approach for investment performance measurement 3rd WCMF, Santa Barbara November 2009 Thaleia Zariphopoulou University of Oxford, Oxford-Man Institute and The University of Texas at Austin 1 Performance
More informationBackward martingale representation and endogenous completeness in finance
Backward martingale representation and endogenous completeness in finance Dmitry Kramkov (with Silviu Predoiu) Carnegie Mellon University 1 / 19 Bibliography Robert M. Anderson and Roberto C. Raimondo.
More informationMSc Dissertation topics:
.... MSc Dissertation topics: Omar Lakkis Mathematics University of Sussex Brighton, England November 6, 2013 Office Location: Pevensey 3 5C2 Office hours: Autumn: Tue & Fri 11:30 12:30; Spring: TBA. O
More informationRegularization by noise in infinite dimensions
Regularization by noise in infinite dimensions Franco Flandoli, University of Pisa King s College 2017 Franco Flandoli, University of Pisa () Regularization by noise King s College 2017 1 / 33 Plan of
More informationApplied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.
Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the
More informationExercises. 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 informationLocal vs. Nonlocal Diffusions A Tale of Two Laplacians
Local vs. Nonlocal Diffusions A Tale of Two Laplacians Jinqiao Duan Dept of Applied Mathematics Illinois Institute of Technology Chicago duan@iit.edu Outline 1 Einstein & Wiener: The Local diffusion 2
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 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 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 informationA random perturbation approach to some stochastic approximation algorithms in optimization.
A random perturbation approach to some stochastic approximation algorithms in optimization. Wenqing Hu. 1 (Presentation based on joint works with Chris Junchi Li 2, Weijie Su 3, Haoyi Xiong 4.) 1. Department
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 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 informationSPDEs, criticality, and renormalisation
SPDEs, criticality, and renormalisation Hendrik Weber Mathematics Institute University of Warwick Potsdam, 06.11.2013 An interesting model from Physics I Ising model Spin configurations: Energy: Inverse
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 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 informationAdvanced computational methods X Selected Topics: SGD
Advanced computational methods X071521-Selected Topics: SGD. In this lecture, we look at the stochastic gradient descent (SGD) method 1 An illustrating example The MNIST is a simple dataset of variety
More informationTHE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 4, Number 1, Winter 1992 THE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS J.-P. KAUTHEN ABSTRACT. We present a method of lines
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 informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
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 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 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 informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More information1. Introduction Systems evolving on widely separated time-scales represent a challenge for numerical simulations. A generic example is
COMM. MATH. SCI. Vol. 1, No. 2, pp. 385 391 c 23 International Press FAST COMMUNICATIONS NUMERICAL TECHNIQUES FOR MULTI-SCALE DYNAMICAL SYSTEMS WITH STOCHASTIC EFFECTS ERIC VANDEN-EIJNDEN Abstract. Numerical
More informationDegree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m.
Degree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m. Candidates should submit answers to a maximum of four
More informationc 2003 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 24, No. 5, pp. 1809 1822 c 2003 Society for Industrial and Applied Mathematics EXPONENTIAL TIMESTEPPING WITH BOUNDARY TEST FOR STOCHASTIC DIFFERENTIAL EQUATIONS KALVIS M. JANSONS
More informationInterest Rate Models:
1/17 Interest Rate Models: from Parametric Statistics to Infinite Dimensional Stochastic Analysis René Carmona Bendheim Center for Finance ORFE & PACM, Princeton University email: rcarmna@princeton.edu
More informationElliptic Operators with Unbounded Coefficients
Elliptic Operators with Unbounded Coefficients Federica Gregorio Universitá degli Studi di Salerno 8th June 2018 joint work with S.E. Boutiah, A. Rhandi, C. Tacelli Motivation Consider the Stochastic Differential
More informationSimulation and Parametric Estimation of SDEs
Simulation and Parametric Estimation of SDEs Jhonny Gonzalez School of Mathematics The University of Manchester Magical books project August 23, 2012 Motivation The problem Simulation of SDEs SDEs driven
More informationStochastic contraction BACS Workshop Chamonix, January 14, 2008
Stochastic contraction BACS Workshop Chamonix, January 14, 2008 Q.-C. Pham N. Tabareau J.-J. Slotine Q.-C. Pham, N. Tabareau, J.-J. Slotine () Stochastic contraction 1 / 19 Why stochastic contraction?
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 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 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 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 informationF1.9AB2 1. r 2 θ2 + sin 2 α. and. p θ = mr 2 θ. p2 θ. (d) In light of the information in part (c) above, we can express the Hamiltonian in the form
F1.9AB2 1 Question 1 (20 Marks) A cone of semi-angle α has its axis vertical and vertex downwards, as in Figure 1 (overleaf). A point mass m slides without friction on the inside of the cone under the
More informationParallel-in-time integrators for Hamiltonian systems
Parallel-in-time integrators for Hamiltonian systems Claude Le Bris ENPC and INRIA Visiting Professor, The University of Chicago joint work with X. Dai (Paris 6 and Chinese Academy of Sciences), F. Legoll
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 informationA Class of Fractional Stochastic Differential Equations
Vietnam Journal of Mathematics 36:38) 71 79 Vietnam Journal of MATHEMATICS VAST 8 A Class of Fractional Stochastic Differential Equations Nguyen Tien Dung Department of Mathematics, Vietnam National University,
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 information