Mean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations
|
|
- Shanon Perkins
- 5 years ago
- Views:
Transcription
1 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 Finance, Probability, and Partial Differential Equation Conference Rutgers University, New Brunswick, May 17-19, 2017
2 Table of contents Introduction 1 Introduction
3 Introduction Consider the system of stochastic differential equations X(s) =x + R s 0 b (X(r)) dr + P R M s k=1 0 k (X(r) k dw k (r), X(0) =x 2 R d, (1.1) where b : R d! R d, k : R d! R, k 2 R d and W k (t) are independent, one dimensional Brownian motions. Assume b and k are such that a weak solution to (1.1) exists and is unique in probability law. We would like to construct an accurate approximations to solutions of (1.1) on fixed time intervals.
4 Introduction SDEs have wide range of applications: Population dynamics, protein kinetics and genetics, psychology and neuronal activity, Math finance, turbulent diffusion and radio-astronomy, Seismology and structural mechanics and so on. As with most ordinary differential equations we can not generally solve SDEs explicitly. Two ways of approximating solution to SDE: Strong approximation. Weak approximation.
5 Introduction We may be interested to approximate E(g(X T )), T > 0. For instance E e rt (X(T) K) +. Better approximation on the probability distribution of solutions is sufficient. Definition Weak convergence of order Ef(X(T)) Ef(Y(N)) apple Kh, N = T h.
6 The Numerical Method Literature review [Anderson and Mattingly, 2011],[Bruti-Liberati and Platen, 2008], [Gaines and Lyons, 1997],[Kloeden and Platen, 1992]. We propose a numerical method to approximate the solution of (1.1) with order of convergence one. Let := {t 0 = 0 < t 1 <... < t N = T} be a subdivision of [0, T]. Consider the same initial condition X(0) =Y 0 = x 0. (i) 1k, (i) : i 2 N, k 2{1, 2,...M} is a collection of mutually 2k independent normal random variables with mean zero and variance 1. Define [a] + = max{a, 0} = a _ 0 for all a 2 R.
7 The Numerical Method Fix 2 (0, 1) and define, 1 = 1 2 (1 ) and 2 = 1 1 = Discretization step h > 0. For i = 1, 2, 3,..., Step 1 : Y 1 = Y i 1 + b(y i 1 ) h + Step 2 : Y i = Y 1 + b(y i 1)(1 MX + k=1 MX k= (1 ) )h q 2 [ 1 k (Y 1 ) 2 2 k (Y i 1)] + k i 2k k (Y i 1 ) k i 1k p h, q (1 )h.
8 Motivation Behind The Algorithm Equation (1.1) is distributionally equivalent to Z t X(t) =X(0)+ + MX Z 1 k k=1 0 b (X(s)) ds 0 Z t 0 I [0, 2 k (X(s))(u) k (du ds). (2.1) k are independent space-time white noise processes. Must approximate k (A [0,h) ( 2 )) to approximate the diffusion k term in (2.1) over the interval [0, h). A [0,h) ( 2 k ) is the region under the curve 2 (X(t)) for 0 apple t apple h. k
9 Motivation Behind The Algorithm Figure: 1 A graphical illustration of Extended Euler-Maruyama Scheme for = 1 2 To approximate X(h), we focus on the double integral in (2.1) for a single k.
10 Motivation Behind The Algorithm Take = 1 2 for illustration. Approximate X( h 2 ) using the Euler-Maruyama method on the interval [0, h 2 ) Denote it by Y 1 To do so, we calculate k (Region 1), where Region 1 is the grey shaded region in the figure. k (Region 1) d = N(0, 2 k (X(0))h 2 ) d = k (X(0)) r h N(0, 1) 2 This is equivalent in distribution to step 1 of our algorithm and exactly carried by step 1 of our algorithm.
11 Motivation Behind The Algorithm Consider the whole region Region 5. k (Region 5) = d N 0, ( 2 k (X(0)) 2( 2 k (X(0)) 2 k (Y 1 ))h 2 d = N 0, (2 2 k (Y 1 ) 2 k (X(0)))h 2 Returning q to step 2 of our q algorithm: For i = 1, we have (2 2 k (Y 1 )) 2 k (X(0)) h 2 N(0, 1) in step 2 of our algorithm.
12 Functional Setting We define C k (R d ) in the following way: C k (R d )={f : R d! R s.t. D f(x) exists, bounded and continuous for all x 2 R d }, (3.1) where is such that i 2 Z + [{0} for each i = 1,...,d and applek. We define the norm kfk k = sup n D f(x) : x 2 R d, =( 1,..., d ), applek o. (3.2) B(R d ) is the family of real-valued, bounded and Borel measurable functions defined on R d.
13 Functional Setting Define the Markov semigroup P t : B(R d )!B(R d ) related to (1.1) by (P t f)(x) : def = E x f(x(t)), where X(0) =x (3.3) The Markov semigroup P h : B(R d )!B(R d ) by. (P h f)(y) : def = E y f(y 1 ), where Y 0 = y (3.4)
14 Assumptions Introduction (A1) The coefficients of (1.1) satisfy the Lipschitz and linear growth conditions: MX b(x) b(y) + k (x) k (y) appleapple x y, b(x) + k=1 MX k (x) appleapple(1 + x ), k=1 (3.5) for all k = 1,...,M and x, y 2 R d, where apple is a positive constant.
15 Assumptions Introduction (A2) For each k = 1,...,M, we have inf x2r d{ k (x)} > 0. In addition, there exists a positive constant 2 (0, 1] so that for any x, 2 R d we have 2 apple T a(x) apple 1 2, (3.6) where T denotes the transpose of and a(x) := P M k=1 (A3) For all multi-index with apple4, we have D b(x) + 2 k (x) k T k. MX D k (x) applek(1 + x p ), for all x 2 R d, (3.7) k=1 where K and p are positive numbers.
16 The Generators Introduction (Af)(x) =f 0 [b(x)](x)+ 1 2 (B 1 f)(x) =f 0 [b(x 0 )](x)+ 1 2 MX k=1 MX k=1 2 k (x)f 00 [ k, k ](x), (3.8) k (x 0 ) 2 f 00 [ k, k ](x), (3.9) (Bf)(x) =f 0 [b(x 0 )](x) (3.10) + 1 MX [ k (Y 1 ) 2 2 k (x 0)] + f 00 [ k, k ](x). k=1 Where f 0 [ ](z) is the derivative of f in the direction evaluated at the point z.
17 Main Theorems Introduction Theorem (Local Approximation ) Assume (A1) (A3). Then there exist a constant K so that P h P h 4!0 apple Kh 2 for all h > 0 sufficiently small. Theorem (Global Approximation ) Assume that b 2 C 4 and 8k, k 2 C 4 with inf x k (x) > 0. Then for any T > 0 there exists a constant C(T) such that sup P nh 0applenhappleT P n h 4!0 apple C(T)h
18 Main Theorems Introduction The proof of the Local Approximation Theorem depends on the following lemma Lemma Assume (A1) (A4). Then for all h > 0 sufficiently small and f 2 C 4, we have E[f(Y 1 )+(Bf)(Y 1 )(1 )h] =f(x 0)+(Af)(x 0 )h + O(h 2 ).
19 Example 1 Geometric Brownian Motion For one dimensional case we have the Black-Scholes model. dx(t) = X(t)dt + µx(t)dw(t), X(0) =x 2 R (4.1) This Stochastic differential equation (4.1) is used to model asset price in financial mathematics. The solution to this SDE is X(t) =X(0)e ( 1 2 µ2 )t+µw(t) (4.2) The expected value of the solution is E[X(t)] = E[X(0)]e t (4.3)
20 Example 1 Geometric Brownian Motion We test the weak convergence of our method over [0, 1] for = 2, µ = 0.1, and X(0) =1. We used step sizes h p = 2 p 10, for 1 apple p apple 5 to sample over discretized brownian paths of the equation (4.1). We then compute the error and the outcome of the numerical experiment is shown in Figure 2 where we have plotted the weak error against h on log-log scale.
21 Example 1 Geometric Brownian Motion Figure: 2 Log-Log plots of error versus step-size.
22 Example 1 Geometric Brownian Motion Detailed calculation shows that E[(X(t)) 2 ]=E[X(0) 2 ]e (2 +µ2 )t (4.4) We test the weak convergence of our method over [0, 1] for = 2, µ = 0.1, and X(0) =1. We used step sizes h p = 2 p 10, for 1 apple p apple 5 to sample over discretized brownian paths of the equation (4.1). The result of the numerical experiment is shown in the Figure 3 where we have plotted the weak error against h on log-log scale.
23 1 Geometric Brownian Motion Figure: 3 Log-Log plots of error versus step-size. It is observed that the slope of the best fit line is empirically one.
24 Background For Roundoff and truncation errors produced during simulation is natural. The spreading of such errors is absolutely necessary to understand in simulation. The usefulness of a numerical method depends upon its ability to control the propogation of errors. Numerical schemes with better numerical stability are demanded.
25 Background For The solutions of a differential equation are continuous in their initial values, at least over a finite time interval. The concept of stability is an extension of this idea to an infinite time interval, [Kloeden and Platen, 1992]. In some sense, the stability of a numerical scheme refers to the conditions under which the impact of an error vanishes asymptotically over time, [Bruti-Liberati and Platen, 2008]. Consider scalar Itô equation dx(t) =b(t, X(t))dt + (t, X(t))dW(t), X(t 0 )=x 0 2 R d (5.1) With a steady solution X(t) 0, so b(t, 0) =0 and (t, 0) =0.
26 Background For We recall the definitions for pth-mean stability and asymptotically stable in pth-mean from [Kloeden and Platen, 1992]. Definition Steady solution X t 0 is called stable in pth-mean if for every >0 and t = (t 0, ) > 0 such that E[ X t 0,x 0 t p ] < for all t t 0 and x 0 <, asymptotically stable in pth- mean if it is stable in pth-mean and there exists a 0 = 0 (t 0 ) > 0 such that lim t!1 E[ X t 0,x 0 t p ]=0 for all x 0 < 0.
27 Background For We give the following probabilistic definitions for stochastic stability from [Kloeden and Platen, 1992] Definition The steady solution X(t) 0 is called stochastically stable if for any >0 and t 0 0 lim P sup X t 0,x 0 x 0!0 t = 0, t t 0 and stochastically asymptotically stable if, in addition, lim x 0!0 P lim t!1 X t 0,x 0 t!0 = 1.
28 Background For Now it is natural to ask the following questions: 1 Do the numerical solutions of SDEs preserve the stability properties of the original SDEs? 2 and if yes for what stepsizes h does the numerical method reproduce the characteristics of the test equation? The above two questions have received quite a lot of attention and in fact stability analysis of numerical methods is motivated by them. We try to answer the above two questions. We will focus on the linear test equation X(t) =X(0)+ Z t 0 X(s)ds + Z t 0 µx(s)dw(s), t 0. (5.2)
29 Mean-square Stability Analysis For Stochastic ODEs The SDE (5.2) can be written as dx(t) = (X(t))dt + µ(x(t))dw(t), t 2 [0, T] X(0) =x 2 R (5.3) We assume that and µ are real constants. Assuming that X(0), 0 with probability 1, solutions of (5.3) have the following properties: lim t!+1 E( X(t) 2 )=0, 2 + µ 2 < 0 (5.4) lim t!+1 X(t) = 0, with probability 1, 1 2 µ2 < 0. (5.5) Here E(.) denotes the expected value.
30 Mean-square Stability Analysis For Stochastic ODEs For simplicity, let S P denotes the sets of order pairs of real problem parameters for which problem is stable. Clearly, S P = n, µ 2 R : 2 + µ 2 < 0 o. Applying the extended Euler-Maruyama method to (5.2) produces the following iterative sequence: Y n = h r i apple A + B (n) Yn (1 )h C + D (n) + E (n) 2 + (n) Y n 1, (5.6) for n = 1, 2,...,where { (n) 1, (n), n = 1, 2,...} are mutually 2 independent normal random variables with mean zero and variance one.
31 Mean-square Stability Analysis For Stochastic ODEs A, B, C, D and E in (5.6) are given by A := 1 + h, B := µ p h, C := µ 2 ( h h) D := 2 1 µ 3 ( h + 1) p h, E := 1 µ 4 h > 0. (5.7) Apparently C is positive when 0. When <0, we can find 0 < h < 1 positive. 1 + p 2 / 1 such that C is
32 Mean-square Stability Analysis For Stochastic ODEs The sequence (5.6) is mean-square stable if lim n!1 E( Y n 2 )=0) [Higham, 2000b]. It follows from (5.6) that E[ Y n 2 ]=E[ Y n 1 2 ] A 2 + B 2 +(1 appleh (n) )he C + D + E (n) 2 i (5.8) Lemma For 0 < h < p 2 / 1, we have apple E [C + D (n) + E( (n) 1 1 )2 ] + apple C + E + o(h 2 ). (5.9)
33 Mean-square Stability Analysis For Stochastic ODEs Putting (5.9) into (5.8) and using the expressions for A,...,E in (5.7), detailed computations reveal that E[ Y n 2 ] < E[ Y n 1 2 ] A 2 + B 2 +(1 )h(c + E + o(h 2 )) " = E[ Y n 1 2 ] 1 +(2 + µ 2 )h [ µ 2 + µ 4 ]h 2 # + o(h 2 ). (5.10) Furthermore, for the expression inside the brackets of the right-hand side of (5.10), we notice that µ 2 + µ 4 = µ2 + 4 µ4 > 0. (5.11)
34 Mean-square Stability Analysis For Stochastic ODEs Next we compute the discriminant =(2 + µ 2 ) ( µ 2 + µ 4 ) = µ 4 < 0. (5.12) Therefore, it follows that for any h > 0, 1 +(2 + µ 2 )h [ µ 2 + µ 4 ]h 2 > 0. The condition for mean square stability of the Extended Euler-Maruyama(EEM) method for (5.2) is 1 +(2 + µ 2 )h [ µ 2 + µ 4 ]h 2 < 1. (5.13)
35 Mean-square Stability Analysis For Stochastic ODEs We can rewrite equation (5.13) as 0 < h < 2(2 + µ 2 ) µ 2 + µ 4. (5.14) Sufficient condition for mean square stability of the weak Simpson method for (5.2) is ( 2(2 + µ 2 ) 0 < h < min µ 2 + µ, p ) 2 / 1. (5.15)
36 Mean-square Stability Analysis For Stochastic ODEs We have the following theorem: Theorem The following assertions are true: (a) Given (, µ) 2 S P, the extended Euler-Maruyama method is mean-square stable if the discretization stepsize h satisfies (5.15). Therefore the mean square stability of the process (5.2) implies the mean square stability of the extended Euler-Maruyama method if the discretization stepsize h satisfies (5.15). (b) Conversely, if the extended Euler-Maruyama method with discretization stepsize h > 0 is mean-square stable for (5.2), then the parameters and µ of (5.2) satisfies (, µ) 2 S P. In other words, the mean square stability of the extended Euler-Maruyama method implies that of the underlying stochastic process (5.2).
37 Mean-square Stability Analysis For Stochastic ODEs We can visualize the stability region for 2(2 + µ 2 ) µ 2 + µ < p 2 / 1. Using (5.13), we have S M := n (x, y) 2 R 2 : 2x 2 + 2xy + y 2 + 4x + 2y < 0 o. (5.16) The plot of mean-square stability domain for the EEM method S M is shown in Figure 9.
38 Mean-square Stability Analysis For Stochastic ODEs Figure: 9 Real mean-square stability domain for weak EEM and EM method
39 Thank You!
40 Anderson, D. F. and Mattingly, J. C. (2011). A weak trapezoidal method for a class of stochastic differential equations. Commun. Math. Sci., 9(1): Asmussen, S. and Glynn, P. W. (2007). Stochastic simulation: algorithms and analysis, volume 57 of Stochastic Modelling and Applied Probability. Springer, New York. Atkinson, K. E. (1989). An introduction to numerical analysis. John Wiley & Sons, Inc., New York, second edition.
41 Bruti-Liberati, N. and Platen, E. (2008). Strong predictor-corrector Euler methods for stochastic differential equations. Stoch. Dyn., 8(3): Ethier, S. N. and Kurtz, T. G. (1986). Markov processes: Characterization and Convergence. Wiley Series in Probability and Mathematical Statistics: Probability and Mathematical Statistics. John Wiley & Sons Inc., New York. Gaines, J. G. and Lyons, T. J. (1997). Variable step size control in the numerical solution of stochastic differential equations.
42 SIAM J. Appl. Math., 57(5): Higham, D. J. (2000a). A-stability and stochastic mean-square stability. BIT, 40(2): Higham, D. J. (2000b). Mean-square and asymptotic stability of the stochastic theta method. SIAM J. Numer. Anal., 38(3): (electronic). Higham, D. J. (2001). An algorithmic introduction to numerical simulation of stochastic differential equations. SIAM Rev., 43(3): (electronic).
43 Karatzas, I. and Shreve, S. E. (1991). Brownian motion and stochastic calculus, volume 113 of Graduate Texts in Mathematics. Springer-Verlag, New York, second edition. Kloeden, P. E. and Platen, E. (1992). Numerical solution of stochastic differential equations, volume 23 of Applications of Mathematics (New York). Springer-Verlag, Berlin. Øksendal, B. (2003). Stochastic differential equations, An introduction with applications. Universitext. Springer-Verlag, Berlin, sixth edition.
44 Saito, Y. and Mitsui, T. (1996). Stability analysis of numerical schemes for stochastic differential equations. SIAM J. Numer. Anal., 33(6): Saito, Y. and Mitsui, T. (2002). Mean-square stability of numerical schemes for stochastic differential systems. Vietnam Journal of Mathematics 30(SI): Sauer, T. (2012). Numerical solution of stochastic differential equations in finance.
45 In Handbook of computational finance, Springer Handb. Comput. Stat., pages Springer, Heidelberg. Talay, D. and Tubaro, L. (1990). Expansion of the global error for numerical schemes solving stochastic differential equations. Stochastic Anal. Appl., 8(4): (1991). Steele, J. (2010). Stochastic Calculus and Financial Applications. Springer-Verlag, Newyork. Khasminskii, R. (2012).
46 Stochastic stability of differential equations, volume 66 of Stochastic Modelling and Applied probability,with contributions by G.N. Milstein and M.B. Nevelson. Heidlberg,second edition. Mao, X. (1997). Stochastic differential equations and application. Horwood, Chichester.
A Weak Simpson Method for a Class of Stochastic Differential Equations and Numerical Stability Results
University of Wisconsin Milwaukee UWM Digital Commons Theses and Dissertations August 15 A Weak Simpson Method for a Class of Stochastic Differential Equations and Numerical Stability Results Ram Sharan
More informationOn Mean-Square and Asymptotic Stability for Numerical Approximations of Stochastic Ordinary Differential Equations
On Mean-Square and Asymptotic Stability for Numerical Approximations of Stochastic Ordinary Differential Equations Rózsa Horváth Bokor and Taketomo Mitsui Abstract This note tries to connect the stochastic
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 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 informationPeter E. Kloeden Eckhard Platen. Numerical Solution of Stochastic Differential Equations
Peter E. Kloeden Eckhard Platen Numerical Solution of Stochastic Differential Equations Peter E. Kloeden School of Computing and Mathematics, Deakin Universit y Geelong 3217, Victoria, Australia Eckhard
More informationIntroduction to Diffusion Processes.
Introduction to Diffusion Processes. Arka P. Ghosh Department of Statistics Iowa State University Ames, IA 511-121 apghosh@iastate.edu (515) 294-7851. February 1, 21 Abstract In this section we describe
More 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 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 informationNumerical Simulation of Stochastic Differential Equations: Lecture 2, Part 1. Recap: SDE. Euler Maruyama. Lecture 2, Part 1: Euler Maruyama
Numerical Simulation of Stochastic Differential Equations: Lecture, Part 1 Des Higham Department of Mathematics Universit of Strathclde Lecture, Part 1: Euler Maruama Definition of Euler Maruama Method
More informationEULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS
Qiao, H. Osaka J. Math. 51 (14), 47 66 EULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS HUIJIE QIAO (Received May 6, 11, revised May 1, 1) Abstract In this paper we show
More 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 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 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 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 informationJournal of Computational and Applied Mathematics. Higher-order semi-implicit Taylor schemes for Itô stochastic differential equations
Journal of Computational and Applied Mathematics 6 (0) 009 0 Contents lists available at SciVerse ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
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 informationSecond order weak Runge Kutta type methods for Itô equations
Mathematics and Computers in Simulation 57 001 9 34 Second order weak Runge Kutta type methods for Itô equations Vigirdas Mackevičius, Jurgis Navikas Department of Mathematics and Informatics, Vilnius
More informationfor all f satisfying E[ f(x) ] <.
. Let (Ω, F, P ) be a probability space and D be a sub-σ-algebra of F. An (H, H)-valued random variable X is independent of D if and only if P ({X Γ} D) = P {X Γ}P (D) for all Γ H and D D. Prove that if
More informationStationary distribution and pathwise estimation of n-species mutualism system with stochastic perturbation
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 6), 936 93 Research Article Stationary distribution and pathwise estimation of n-species mutualism system with stochastic perturbation Weiwei
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 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 informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
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 informationConvergence rate estimates for the gradient differential inclusion
Convergence rate estimates for the gradient differential inclusion Osman Güler November 23 Abstract Let f : H R { } be a proper, lower semi continuous, convex function in a Hilbert space H. The gradient
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 informationEXPONENTIAL MEAN-SQUARE STABILITY OF NUMERICAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL EQUATIONS
London Mathematical Society ISSN 1461 157 EXPONENTIAL MEAN-SQUARE STABILITY OF NUMERICAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL EQUATIONS DESMOND J. HIGHAM, XUERONG MAO and ANDREW M. STUART Abstract Positive
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 informationGillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde
Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation
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 informationSome Properties of NSFDEs
Chenggui Yuan (Swansea University) Some Properties of NSFDEs 1 / 41 Some Properties of NSFDEs Chenggui Yuan Swansea University Chenggui Yuan (Swansea University) Some Properties of NSFDEs 2 / 41 Outline
More informationA MODEL FOR THE LONG-TERM OPTIMAL CAPACITY LEVEL OF AN INVESTMENT PROJECT
A MODEL FOR HE LONG-ERM OPIMAL CAPACIY LEVEL OF AN INVESMEN PROJEC ARNE LØKKA AND MIHAIL ZERVOS Abstract. We consider an investment project that produces a single commodity. he project s operation yields
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationOn the quantiles of the Brownian motion and their hitting times.
On the quantiles of the Brownian motion and their hitting times. Angelos Dassios London School of Economics May 23 Abstract The distribution of the α-quantile of a Brownian motion on an interval [, t]
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations Simo Särkkä Aalto University, Finland November 18, 2014 Simo Särkkä (Aalto) Lecture 4: Numerical Solution of SDEs November
More 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 informationBernardo D Auria Stochastic Processes /12. Notes. March 29 th, 2012
1 Stochastic Calculus Notes March 9 th, 1 In 19, Bachelier proposed for the Paris stock exchange a model for the fluctuations affecting the price X(t) of an asset that was given by the Brownian motion.
More information(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS
(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS Svetlana Janković and Miljana Jovanović Faculty of Science, Department of Mathematics, University
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 informationAPPLICATION OF THE KALMAN-BUCY FILTER IN THE STOCHASTIC DIFFERENTIAL EQUATIONS FOR THE MODELING OF RL CIRCUIT
Int. J. Nonlinear Anal. Appl. 2 (211) No.1, 35 41 ISSN: 28-6822 (electronic) http://www.ijnaa.com APPICATION OF THE KAMAN-BUCY FITER IN THE STOCHASTIC DIFFERENTIA EQUATIONS FOR THE MODEING OF R CIRCUIT
More informationUNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE
Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (2010), 275 284 UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE Iuliana Carmen Bărbăcioru Abstract.
More 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 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 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 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 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 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 informationQualitative behaviour of numerical methods for SDEs and application to homogenization
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,
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 informationWeak convergence and large deviation theory
First Prev Next Go To Go Back Full Screen Close Quit 1 Weak convergence and large deviation theory Large deviation principle Convergence in distribution The Bryc-Varadhan theorem Tightness and Prohorov
More informationDUBLIN CITY UNIVERSITY
DUBLIN CITY UNIVERSITY SAMPLE EXAMINATIONS 2017/2018 MODULE: QUALIFICATIONS: Simulation for Finance MS455 B.Sc. Actuarial Mathematics ACM B.Sc. Financial Mathematics FIM YEAR OF STUDY: 4 EXAMINERS: Mr
More informationSquared Bessel Process with Delay
Southern Illinois University Carbondale OpenSIUC Articles and Preprints Department of Mathematics 216 Squared Bessel Process with Delay Harry Randolph Hughes Southern Illinois University Carbondale, hrhughes@siu.edu
More informationAsymptotic behaviour of the stochastic Lotka Volterra model
J. Math. Anal. Appl. 87 3 56 www.elsevier.com/locate/jmaa Asymptotic behaviour of the stochastic Lotka Volterra model Xuerong Mao, Sotirios Sabanis, and Eric Renshaw Department of Statistics and Modelling
More informationOPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS
APPLICATIONES MATHEMATICAE 29,4 (22), pp. 387 398 Mariusz Michta (Zielona Góra) OPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS Abstract. A martingale problem approach is used first to analyze
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 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 informationHow does it work? QM describes the microscopic world in a way analogous to how classical mechanics (CM) describes the macroscopic world.
Today we use Quantum Mechanics (QM) on a regular basis: Chemical bonds NMR spectroscopy The laser (blue beam in Blue-ray player; red beam in a DVD player for example) The fingerprint /signature of a molecule
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationOn the relation between the Smith-Wilson method and integrated Ornstein-Uhlenbeck processes
Mathematical Statistics Stockholm University On the relation between the Smith-Wilson method and integrated Ornstein-Uhlenbeck processes Håkan Andersson Mathias Lindholm Research Report 213:1 ISSN 165-377
More informationApproximating diffusions by piecewise constant parameters
Approximating diffusions by piecewise constant parameters Lothar Breuer Institute of Mathematics Statistics, University of Kent, Canterbury CT2 7NF, UK Abstract We approximate the resolvent of a one-dimensional
More informationB8.3 Mathematical Models for Financial Derivatives. Hilary Term Solution Sheet 2
B8.3 Mathematical Models for Financial Derivatives Hilary Term 18 Solution Sheet In the following W t ) t denotes a standard Brownian motion and t > denotes time. A partition π of the interval, t is a
More 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 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 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 informationSTAT 331. Martingale Central Limit Theorem and Related Results
STAT 331 Martingale Central Limit Theorem and Related Results In this unit we discuss a version of the martingale central limit theorem, which states that under certain conditions, a sum of orthogonal
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 informationBernardo D Auria Stochastic Processes /10. Notes. Abril 13 th, 2010
1 Stochastic Calculus Notes Abril 13 th, 1 As we have seen in previous lessons, the stochastic integral with respect to the Brownian motion shows a behavior different from the classical Riemann-Stieltjes
More informationStochastic Integration (Simple Version)
Stochastic Integration (Simple Version) Tuesday, March 17, 2015 2:03 PM Reading: Gardiner Secs. 4.1-4.3, 4.4.4 But there are boundary issues when (if ) so we can't apply the standard delta function integration
More informationDownloaded 12/15/15 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 40, No. 4, pp. 56 537 c 2002 Society for Industrial and Applied Mathematics PREDICTOR-CORRECTOR METHODS OF RUNGE KUTTA TYPE FOR STOCHASTIC DIFFERENTIAL EQUATIONS KEVIN BURRAGE
More informationThe Continuity of SDE With Respect to Initial Value in the Total Variation
Ξ44fflΞ5» ο ffi fi $ Vol.44, No.5 2015 9" ADVANCES IN MATHEMATICS(CHINA) Sep., 2015 doi: 10.11845/sxjz.2014024b The Continuity of SDE With Respect to Initial Value in the Total Variation PENG Xuhui (1.
More informationGaussian, Markov and stationary processes
Gaussian, Markov and stationary processes Gonzalo Mateos Dept. of ECE and Goergen Institute for Data Science University of Rochester gmateosb@ece.rochester.edu http://www.ece.rochester.edu/~gmateosb/ November
More information2008 Hotelling Lectures
First Prev Next Go To Go Back Full Screen Close Quit 1 28 Hotelling Lectures 1. Stochastic models for chemical reactions 2. Identifying separated time scales in stochastic models of reaction networks 3.
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 informationNUMERICAL ANALYSIS OF STOCHASTIC DIFFERENTIAL EQUATIONS WITH EXPLOSIONS. 1. Introduction
NUMERICAL ANALYSIS OF STOCHASTIC DIFFERENTIAL EQUATIONS WITH EXPLOSIONS JUAN DÁVILA, JULIAN FERNÁNDEZ BONDER, PABLO GROISMAN, JULIO D. ROSSI, AND MARIELA SUED Abstract. Stochastic ordinary differential
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 informationON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER
ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER GERARDO HERNANDEZ-DEL-VALLE arxiv:1209.2411v1 [math.pr] 10 Sep 2012 Abstract. This work deals with first hitting time densities of Ito processes whose
More informationApplied Mathematics Letters. Stationary distribution, ergodicity and extinction of a stochastic generalized logistic system
Applied Mathematics Letters 5 (1) 198 1985 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Stationary distribution, ergodicity
More informationAn Analysis of Katsuura s Continuous Nowhere Differentiable Function
An Analysis of Katsuura s Continuous Nowhere Differentiable Function Thomas M. Lewis Department of Mathematics Furman University tom.lewis@furman.edu Copyright c 2005 by Thomas M. Lewis October 14, 2005
More informationWeakly interacting particle systems on graphs: from dense to sparse
Weakly interacting particle systems on graphs: from dense to sparse Ruoyu Wu University of Michigan (Based on joint works with Daniel Lacker and Kavita Ramanan) USC Math Finance Colloquium October 29,
More informationScaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations
Scaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations Jan Wehr and Jack Xin Abstract We study waves in convex scalar conservation laws under noisy initial perturbations.
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 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 informationPathwise uniqueness for stochastic differential equations driven by pure jump processes
Pathwise uniqueness for stochastic differential equations driven by pure jump processes arxiv:73.995v [math.pr] 9 Mar 7 Jiayu Zheng and Jie Xiong Abstract Based on the weak existence and weak uniqueness,
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 informationStochastic Numerical Analysis
Stochastic Numerical Analysis Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Stoch. NA, Lecture 3 p. 1 Multi-dimensional SDEs So far we have considered scalar SDEs
More informationBROWNIAN BRIDGE HYPOTHESIS TESTING FOR THE INITIAL TRANSIENT PROBLEM
Proceedings of the 211 Winter Simulation Conference S. Jain, R. R. Creasey, J. Himmelspach, K. P. White, and M. Fu, eds. BROWNIAN BRIGE HYPOTHESIS TESTING FOR THE INITIAL TRANSIENT PROBLEM Peter W. Glynn
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 informationReaction-Diffusion Equations In Narrow Tubes and Wave Front P
Outlines Reaction-Diffusion Equations In Narrow Tubes and Wave Front Propagation University of Maryland, College Park USA Outline of Part I Outlines Real Life Examples Description of the Problem and Main
More informationFinite element approximation of the stochastic heat equation with additive noise
p. 1/32 Finite element approximation of the stochastic heat equation with additive noise Stig Larsson p. 2/32 Outline Stochastic heat equation with additive noise du u dt = dw, x D, t > u =, x D, t > u()
More informationWeak solutions of mean-field stochastic differential equations
Weak solutions of mean-field stochastic differential equations Juan Li School of Mathematics and Statistics, Shandong University (Weihai), Weihai 26429, China. Email: juanli@sdu.edu.cn Based on joint works
More informationScaling Limits of Waves in Convex Scalar Conservation Laws Under Random Initial Perturbations
Journal of Statistical Physics, Vol. 122, No. 2, January 2006 ( C 2006 ) DOI: 10.1007/s10955-005-8006-x Scaling Limits of Waves in Convex Scalar Conservation Laws Under Random Initial Perturbations Jan
More informationIntroduction to the Numerical Simulation of Stochastic Differential Equations with Examples
Introduction to the Numerical Simulation of with Examples Prof. Michael Mascagni Department of Computer Science Department of Mathematics Department of Scientific Computing Florida State University, Tallahassee,
More informationDerivation of SPDEs for Correlated Random Walk Transport Models in One and Two Dimensions
This article was downloaded by: [Texas Technology University] On: 23 April 2013, At: 07:52 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationStochastic perspective and parameter estimation for RC and RLC electrical circuits
Int. J. Nonlinear Anal. Appl. 6 (5 No., 53-6 ISSN: 8-68 (electronic http://www.ijnaa.semnan.ac.ir Stochastic perspective and parameter estimation for C and LC electrical circuits P. Nabati a,,. Farnoosh
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 informationFinite-time Ruin Probability of Renewal Model with Risky Investment and Subexponential Claims
Proceedings of the World Congress on Engineering 29 Vol II WCE 29, July 1-3, 29, London, U.K. Finite-time Ruin Probability of Renewal Model with Risky Investment and Subexponential Claims Tao Jiang Abstract
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 informationMan Kyu Im*, Un Cig Ji **, and Jae Hee Kim ***
JOURNAL OF THE CHUNGCHEONG MATHEMATICAL SOCIETY Volume 19, No. 4, December 26 GIRSANOV THEOREM FOR GAUSSIAN PROCESS WITH INDEPENDENT INCREMENTS Man Kyu Im*, Un Cig Ji **, and Jae Hee Kim *** Abstract.
More informationON THE HÖLDER CONTINUITY OF MATRIX FUNCTIONS FOR NORMAL MATRICES
Volume 10 (2009), Issue 4, Article 91, 5 pp. ON THE HÖLDER CONTINUITY O MATRIX UNCTIONS OR NORMAL MATRICES THOMAS P. WIHLER MATHEMATICS INSTITUTE UNIVERSITY O BERN SIDLERSTRASSE 5, CH-3012 BERN SWITZERLAND.
More informationSome Terminology and Concepts that We will Use, But Not Emphasize (Section 6.2)
Some Terminology and Concepts that We will Use, But Not Emphasize (Section 6.2) Statistical analysis is based on probability theory. The fundamental object in probability theory is a probability space,
More informationAn introduction to rough paths
An introduction to rough paths Antoine LEJAY INRIA, Nancy, France From works from T. Lyons, Z. Qian, P. Friz, N. Victoir, M. Gubinelli, D. Feyel, A. de la Pradelle, A.M. Davie,... SPDE semester Isaac Newton
More information