Approximation of random dynamical systems with discrete time by stochastic differential equations: I. Theory
|
|
- Annabel Sherman
- 5 years ago
- Views:
Transcription
1 Random Operators / Stochastic Eqs. 15 7, 5 c de Gruyter 7 DOI / ROSE.7.13 Approximation of random dynamical systems with discrete time by stochastic differential equations: I. Theory Yuri A. Godin and Stanislav Molchanov Communicated by V. L. Girko Abstract. Approximation of one-dimensional stochastic differential equations and their additive functionals by dynamical systems with piecewise-constant random coefficients is obtained. We calculate asymptotic expansion of solution in terms of the step of discretization. Key words. Itô s stochastic equations; discretization; Lyapunov exponent; density of states. AMS classification. 58J5, 35J1, 35P99, 8B44, 35P5. 1. Introduction This paper was motivated by the analysis of one-dimensional Schrödinger operator d + ξt, dt with the white noise potential ξ. In this case one can introduce the usual phaseamplitude representation ψ = r λ sin ϑ λ, ψ = r λ cosϑ λ or use the Prüfer phase 1] for the solution of the Cauchy problem Hψ = λψ and get a stochastic differential equation SDE for ϑ λ t. The white noise potential here can be understood in several senses, including the Dirichlet form or Stratonovich SDE. Physical applications, however, are always dealing not directly with the white noise but rather with its regularization, i.e., a stationary Gaussian process having almost δ- correlator and, equivalently, a wide band spectrum. As an approximation, stochastic differential equations are obtained in the limit. To estimate the quality of this approximation, we have to calculate the remainder terms either in SDE for the phase or in stochastic integrals for the associated additive functional related to the Lyapunov exponent or the density of states, see ] 4]. Similar to ODEs, there exists an extensive literature on numerical methods for SDEs see for example 5] and references therein. However, unlike ODEs, we start from a dynamical system with discrete randomness that is more adequate physically, and approximate it or its functional by an SDE and stochastic integral. In doing so, we
2 6 Y. A. Godin and S. Molchanov consider only one-dimensional case. All calculations are performed for general equations in order to make them more transparent. The obtained result can also be applied in Monte Carlo simulation of SDEs or dynamical systems with discrete randomness. The second part of the paper will be devoted to the transition from localization pure point spectrum to the absolutely continuous spectrum in the spirit of the classical paper 6] and recent publications on the fast oscillating piecewise constant potentials 7] 8].. Description of the model Let w s = ws be the standard Brownian motion whose increment w n on the intervals of length 1 is a small parameter are independent N, random variables, w n = wn + 1 wn, n, w =. The simplest regularization ξ t of the white noise ξt = ẇt has the form ξ t = wn, t n,n + 1, n..1 If the origin t = is uniformly distributed on, ], the process ξ t will be a stationary Gaussian process with the correlator 1 t b t = E ξ ξ + t =, b t δ t if. in the sense of distributions. In what follows we put ξ t = w, t,..3 Assume now that Ax,Bx,x R, are smooth bounded functions. We consider the dynamical system ẋt = ξ tax + Bx..4 On each interval t n,n + 1] equation.4 has the form of usual ODE ẋt = wn Ax + Bx..5 The sequence xn, n, is the main goal of our study. The r.v. x n form usual Markov chain on R and its analytical properties are not simple. We prefer to work with SDE which has symbolically the form for dxt = Axtdw t + Bxtdt..6
3 Dynamical system with discrete randomness 7 Our goal is the analysis of the remainder in equation.6 which must be understood in the Stratonovich form 9]. We also want to study approximation by the stochastic integral of the following additive functional of x. Let i.e., for t = n y n = y n = ẏ t = ξ tfx t + Gx t,.7 n 1 k= w k k+1 k Fx s ds + k+1 In the first part of the paper we will construct the asymptotic expansion for or δxn = xn + 1 xn δyn = yn + 1 yn k Gx s ds..8 up to the order o using standard agreement that w = O 1/, w = O 3/, etc., since Ew =. Proposition.1. Let t = n,xs, ] be the solution of equation.6 on the interval n,n + 1]. Then for w = w n δx t = x t + x t = w A + B + w + 1 AA + 1 w A B + AB w3 A A + AA + 1 BB w AA B + AA B + A B + A B w4 A 3 A + 4A A + AA 3 + o..9 The remainder R in this formula has the form ER = E δx t w A w4 A 3 A + 4A A + AA 3] O 5, i.e., in fact R O 5/ in probability..1 The proof is based on direct calculations. Since solution xs, s, ], is a smooth function of the initial condition x and time t, it can be expanded into a Taylor series xt + s = x + ẋts xs ẍts 6 4 x4 s 4 + os 4..11
4 8 Y. A. Godin and S. Molchanov But ẋt + s = w Axt + s + B xt + s,.1 and for s = we have ẋt = w Ax + B x..13 Similarly, ẍ = w A + B ẋ = w A + B w A + B.14 and so on. Then.9 is obtained by substitution of these formulas in.11 and truncating the terms of the order O α, α 5/. In the case of additive functional given by.7 calculations are similar: δy t = y n+1 y n = y n + 1 yn ] n+1 w n = Fx s + Gx s ds n = ẏn + 1 ÿn yn y4 n + os For t = n, x = x n, w = w n the increments of the derivatives can be calculated as follows The final result represents Proposition.. ẏ = w Fx + Gx ÿ = w F x + G x ẋ = w F x + G x w Ax + Bx = w AF + w F B + AG + BG. δy t = w Fx + Gx + 1 w AF + 1 w BF + AG + 1 BG w3 A F + AA F w ABF + A BF + AB F + AA G + A G w4 A 3 F + 3A A F + A A F + AA F + o..16
5 Dynamical system with discrete randomness 9 Observe that both δx n and δy n are polynomials of only one Gaussian random variable w n, that makes numerical simulation Monte Carlo of the dynamical system x and the additive functional y very simple in contrast to simulation of SDEs see below. Consider a general SDE on the probabilistic space of the same Brownian motion W s,s, as our dynamical system x t dx t = αx t dw t + βx t dt, y = x,.17 with smooth coefficients α and β. We will select them in the form from the condition that the expression E δx n δxn ] αx = α x + α 1 x + o,.18 βx = β x + β 1 x + o,.19 x = E n + 1 xn + 1 xn=xn=x]. attains minimum value. On the interval t,t + ] = n,n + 1], with condition xn = x equation.17 for xs, s t, t + ], has integral Itô s form xs = x + αx s dw s + which can be solved by successive approximations as follows: x s = x, x 1 s = x + αxw s + βxs, x s = x + αx + αxw u + βxudw u + = x + αxw s + βxs + α x + β x αx + αxβx + 1 β xα x αx βx s ds,.1 βx + αxw u + βxudu, w u dw u + βx w u du + 1 βxs ] + 1 α x uw u dw u ] w u du α xα 3 x α x udw u ] w u dw u w 3 u du + o..
6 1 Y. A. Godin and S. Molchanov The third approximation gives the estimate to the order O One can easily see x 3 s = x + x 3 s = x s + α x + + α xβx + αxα xβ x + α xα xα x αx udw u + αxα x u ] v dw v dw u u w v dw v du βx udu..3 u w v dw v dw u u w u w v dw v dw u + o..4 Arranging in.,.4 the terms in their order in, we obtain Proposition.3. Expansion of δxn for x = xn, s, ] has the form δxn = xn + 1 xn = αxw + βx + αxα x + α xβx sdw s + αxβ x w s ds + 1 α xα x + αxα x w u dw u dw s + 1 βxβ x + αxα xβx α3 xα x + αxα xβ x + α xα xα x sw s dw s + 1 α xβ x w 3 s dw s + α x βx w u dw u ds ws w u dw u w s sw s + o. w s dw s w s dw s ds.5 udw u dw s Stochastic integrals in.5 can be written in the form containing random variable w and uncorrelated terms or even independent on w. Proposition.4. Integrals in.5 can be transformed using the formulas below. Proofs are based on the Itˆo formula and integration by parts.
7 Dynamical system with discrete randomness 11 a b c d e f g h i w s dw s = 1 w sw s dw s = 1 w sdw s = w udw u w s ds dw s = w 4 3 ws dw s = 1 3 w3 w s ds w u dw u dw s = 1 6 w3 1 w w s ds w s ds ws 3 dw s = 1 4 w4 3 ws ds w u dw u ds = 1 ws ds 1 4 w u dw u dw s = 1 4 w4 3 4 w w s ds.6 Observe that along with and w, δxn contains additional random variables w s ds and w s ds. These random variables have a joint Gaussian distribution with covariance matrix cov = Hence, standard calculation gives w s dw s = 1 w + where w s is an independent copy of w s, and it easy to see that E w s ds 1 w w = w,.8 and var w s ds = var 1 w + var 3 6 w = 3 3.
8 1 Y. A. Godin and S. Molchanov Note also that w is uncorrelated with all terms of order O in formulas.5.6. Indeed, after the transformation w s w s the terms of order O preserve their sign. For example, in.6 Ew s Therefore, w 3 s dw s = E w s Ew s w s 3 d w s = Ew s w 3 s dw s =. w 3 s dw s..9 Now we are ready to find optimal expressions for the diffusion coefficients αx and βx in in terms of the coefficients Ax and Bx of the dynamical system.4. To this end, we have to compare expansions in Proposition.1 with those in Proposition.3, taking into account formulas.6. Since for t = n and x t + x t = δx = Axw + Bx + AxA x + w AxA x + o.3 xt + xt = δx = αxw + βx + w we have Thus, we put αx = Ax + O, βx = Bx + 1 AxA x + o. αxα x + o,.31 α x = Ax,.3 β x = Bx + 1 AxA x..33 Note that the drift coefficient βx contains the Stratonovich term 1 AxA x. That means that in the transition from the difference relation δx t = w Axt + Bxt, t = n,.34 when we obtain in the leading term not the Itô quation dxt = Axtdw t + Bxtdt,.35
9 Dynamical system with discrete randomness 13 but rather the Stratonovich one dxt = Axt dw t + Bxtdt = Axtdw t + Bxt + 1 AA xt dt..36 In order to find next corrections α 1 and β 1, we put αx = α x + α 1 x = Ax + α 1 x + o, βx = β x + β 1 x = Bx + 1 AxA x + β 1 x + o. Compare now the terms of order O 3/ in.9 and.5. In.9 they are while in.5 we have 1 w A B + AB w3 A A + AA.37 w α 1 x + α β + 1 α α sdw s + α β w s ds ws dw s + α α w u dw u dw s..38 Using formulas.6 we transform the last expression to the following 1 α α + α α + w α β 1 α α 6 w3 + α β α β 1 α α w s ds..39 The terms containing w 3 are the same in both expressions.37 and.39. Comparison of the coefficients of w with.8 gives or α 1 x+α β + 1 α α + 1 α β α β 1 α α = 1 A B + AB.4 α α β + α β 1 4 α α + α α = α A B + 1 AA + 1 B A + 1 AA + A 1 AA + A A 4 = α A B + AB from.37 = 1 A B + AB..41
10 14 Y. A. Godin and S. Molchanov Thus, α 1 x =..4 However, we have a discrepancy between.37 and.39. The expression 3 6 w α β α β 1 α α = 3 6 AB A B. It gives contribution to the diffusion coefficient of the order O. Indeed, Dx = 1 varaxw + w Dx = A x + D 1 x = A x + O. Now we can find correction β 1 x to the drift coefficient.43 βx = Bx + 1 AxA x + β 1 x = β x + β 1 x..44 Let us take expectations in.9 and. and compare the terms of order O. In.9 such terms are In. they are 1 BB AA B + AA B + A B + A B A 3 A + 4A A A + AA β 1 x + 1 β β α β = β 1 x + 1 B + 1 AA B + 1 AA + A A B + 1 ] A A + AA,.46 where we used the identities E Ew 4 =, w ds = 1. From one can find the following expression for β 1 : β 1 x = 1 AA B A B + AA B A B Thus, we have proved the following
11 Dynamical system with discrete randomness 15 Theorem.5. The best approximation to the order of the dynamical system ẋ = ξ tax t + Bx t, by the SDE is given by Here ξ t = wn, t n,n + 1],.48 dxt = αx t dw w + βx t dt.49 αx = Ax + o,.5 βx = Bx + 1 AA + 1 E δx n δxn ] AA B A B + AA B A B + o. = E x n + 1 x n xn + 1 xn ] = O 4,.51 that is, x t xt = O in probability. The differential operator associated with the diffusion process xt has the form A = 1 α x d dx + βx d ] 1 d dx = A x + o dx.5 + Bx + 1 AA x + AA B A B + AA B A B ] d + o 1 dx. 3. Approximation of the additive functional Let us apply the techniques developed in the previous section to approximation of the additive functional y t = y + t ξ sfx sds + t Gx sds, 3.1 δy t = y n + 1 y n, t = n. 3. We have to expand δyt = yn + 1 yn
12 16 Y. A. Godin and S. Molchanov with respect to the diffusion process xt δyt = fδxsdw s + gδxs ds, 3.3 and then find fx = f x + f 1 x + o, 3.4 gx = g x + g 1 x + o, 3.5 from the condition that E δy δy = min. 3.6 Calculation of δyt is straightforward δyn = fx + δxsdw s + gx + δssds, 3.7 where δxs = αxw s + βxs + αxα x + αxβ x w u du + 1 α xα x w u dw u + α xβx w u du udw u u + αxα x dw u w v dw v + o 3.8 and αx = Ax + o, βx = Bx + 1 AxA x + β 1 x. 3.9
13 Dynamical system with discrete randomness 17 Evaluation of the first integral in 3.7 gives = fx + δxsdw s fx + f xδxs + 1 f xδ xs + 1 ] 6 f xδ 3 xs + o dw s = fxw + f x αx + αxα x + αxβ x w s dw s + βx dw s w u dw u + α xβx dw s sdw s w u du + 1 α xα x dw s udw u dw s wu dwu + αxα x s u dw s dw u w v dw v +o ]+ 1 f α x + αxβx f xα 3 x sw s dw s +α xα x ] w s dw s w u dw u + o w s dw s w 3 s dw s + o. 3.1 The second integral in 3.7 is expanded as follows gx + δxsds = gx + g x αx + αxα x + 1 g x α x ds w s ds + 1 βx ] w u dw u +o ] ws ds + o Similar to the case of xt, we have δyt = fxw + gx 1 αxf x + 1 f xαxw + o. 3.1 Due to Proposition., δy t = Fxw + Gx + 1 AxF xw + o. 3.13
14 18 Y. A. Godin and S. Molchanov This means that fx = Fx + o = f x + o, gx = Gx + 1 AxF x + o = g x + o, 3.14 f x = Fx, g x = Gx + 1 AxF x. Let us find the next corrections fx = f x + f 1 x + o, 3.15 gx = g x + g 1 x + o The terms of order O 3/ in Proposition. are the following 1 w AxG x+bxf x w3 A xf x+axa xf x Using Proposition.3 we can find terms of the same order in w f 1 x + f xβ x sdw s + αxα xf x + 1 f xα x ws dw s + α xg x w s ds = w f 1 x + f β w + 1 w f xα 3 x 3 dw s w u dw u w w s ds + α xα xf x 3 6 w w s ds + α xg x w s ds
15 Dynamical system with discrete randomness 19 = w f x + f xβ x 1 α xα xf x w3 + α xα xf x + α xf x w s ds α xg xβ xf x 1 f xα x = w f 1 x + f xβ x 1 α xα xf x + 1 α xg x 1 β xf x 1 4 f xα x w3 + α xα xf x + α xf x 3 6 w α xg x β xf x f xα x The terms containing w 3 in are the same, f x = Fx, α x = Ax Thus, for determining f 1 x we obtain the equation f 1 x F x Bx + 1 AxA x + Ax Gx + 1 AxF x AxA xf x A xf x 4 = AxG x + BxF x Hence, as in the previous case, f 1 x =. 3. In order to find the correction to the drift g 1 x, we have to compute the mean values of δyt and δy t at t = n. Eδy t = Gx + 1 axf x + 1 BxG x + AxBxF x 6 + AxB xf xa xbxf xaxa xg x + A xg x + A 3 xf x + 3A xa xf x + AxA x F x 8 + A xa xf x + o 3.1
16 Y. A. Godin and S. Molchanov and Eδyt = g x + g 1 x + g xβ x + α xg x 4 = + + g 1 x Gx + AxA x G x + A xf x + AxF x Bx + AxA x + 4 From the relation one can find expression for g 1 x G + A xf x + A xf x + AxF x Eδy t = Eδyt. 3. g 1 x = 1 AxBxF + AxB xf x A xbxf x AxA xg x A xg x. 3.3 Theorem 3.1. The best approximation of the additive functional δy t by the stochastic integral of xt has the form for δyt = fx = Fx + o, fx s dw s + gx s ds 3.4 gx = Gx + 1 AxF x + g 1 x + o, 3.5 g 1 x = 1 AxBxF + AxB xf x A xbxf x AxA xg x A xg x. Similar to Theorem.5, E δy n δyn = O 4, 3.6 that is, δy n δyn = O in probability.
17 Dynamical system with discrete randomness 1 Let dx t = αx t dwt + βx t dt 3.7 be a stochastic equation and yt = y + t fx s dw s + t gx s ds 3.8 be a linear functional. Below we formulate a proposition regarding the expectation λ y + ut,x = E x e t fx s dw s + which is a combination of Kac Feynman and Girsanov s formulas. t ] gx s ds, 3.9 Proposition 3.. The function ut, x is a solution of the following parabolic problem ut, x t = α u,x = e λy. To prove 3.3 observe that λ y + ut +,x = E x e That means u x + β + λfα u x + λg + λ fα u, 3.3 t fx s dw s + t ] gx s ds ut,x = E x e λαxfxw + λgx + f α w ] + o ut,x = E x 1 + λαfw + λgx + f α w + 1 λ α f w + o ] ut,x + u x αw + β + 1 ] u xx α w + o ut +,x ut,x Passing here to the limit we obtain 3.3. = λgxu + u xβ + λα fu x + 1 u xxα + 1 λ α f + o. 3.3
18 Y. A. Godin and S. Molchanov References 1. P. Hartman, Ordinary Differential Equations. John Wiley, New York, R. Carmona, J. Lacroix, Spectral Theory of Random Schrödinger operators. Birkhäuser, Boston, L. Pastur, A. Figotin, Spectra of Random and Almost Periodic Operators, Springer-Verlag, New York, S. A. Molchanov, Ideas in the theory of random media. Acta Appl. Math. 1991, P. E. Kloeden, E. Platen, Numerical Solution of Stochastic Differential Equations. Springer-Verlag, New York, S. Kotani, N. Ushiroya, One-dimensional Schrödinger operators with random decaying potentials. Communications in Mathematical Physics , A. Laptev, S. Naboko, O. Safronov, Absolutely continuous spectrum of Schrödinger operators with slowly decaying and oscillating potentials. Communications in Mathematical Physics 53 5, J. Holt, S. Molchanov, O. Safronov, Cwikel and quasi-szegö type estimates for random operators, to appear in Communications in Partial Differential Equations. 9. I. I. Gihman, A. V. Skorohod, Stochastic Differential Equations. Springer-Verlag, New York, 197. Received 1 May, 6; revised 5 March, 7 Author information Yuri A. Godin, Department of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC 83, USA. ygodin@uncc.edu Stanislav Molchanov, Department of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC 83, USA. smolchan@uncc.edu
Propagation of longitudinal waves in a random binary rod
Downloaded By: [University of North Carolina, Charlotte At: 7:3 2 May 28 Waves in Random and Complex Media Vol. 6, No. 4, November 26, 49 46 Propagation of longitudinal waves in a random binary rod YURI
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 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 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 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 informationwhere r n = dn+1 x(t)
Random Variables Overview Probability Random variables Transforms of pdfs Moments and cumulants Useful distributions Random vectors Linear transformations of random vectors The multivariate normal distribution
More informationSession 1: Probability and Markov chains
Session 1: Probability and Markov chains 1. Probability distributions and densities. 2. Relevant distributions. 3. Change of variable. 4. Stochastic processes. 5. The Markov property. 6. Markov finite
More 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 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 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 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 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 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 informationNumerical methods for solving stochastic differential equations
Mathematical Communications 4(1999), 251-256 251 Numerical methods for solving stochastic differential equations Rózsa Horváth Bokor Abstract. This paper provides an introduction to stochastic calculus
More informationNumerical 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 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 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 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 informationLAN property for sde s with additive fractional noise and continuous time observation
LAN property for sde s with additive fractional noise and continuous time observation Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Samy Tindel (Purdue University) Vlad s 6th birthday,
More information3 Applications of partial differentiation
Advanced Calculus Chapter 3 Applications of partial differentiation 37 3 Applications of partial differentiation 3.1 Stationary points Higher derivatives Let U R 2 and f : U R. The partial derivatives
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 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 informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
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 informationMean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations
Mean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations Ram Sharan Adhikari Assistant Professor Of Mathematics Rogers State University Mathematical
More informationStochastic 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 informationIntertwinings for Markov processes
Intertwinings for Markov processes Aldéric Joulin - University of Toulouse Joint work with : Michel Bonnefont - Univ. Bordeaux Workshop 2 Piecewise Deterministic Markov Processes ennes - May 15-17, 2013
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 informationELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process
Department of Electrical Engineering University of Arkansas ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process Dr. Jingxian Wu wuj@uark.edu OUTLINE 2 Definition of stochastic process (random
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 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 informationExact Linearization Of Stochastic Dynamical Systems By State Space Coordinate Transformation And Feedback I g-linearization
Applied Mathematics E-Notes, 3(2003), 99-106 c ISSN 1607-2510 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Exact Linearization Of Stochastic Dynamical Systems By State Space Coordinate
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 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 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 informationADAPTIVE WEAK APPROXIMATION OF REFLECTED AND STOPPED DIFFUSIONS
ADAPTIVE WEAK APPROXIMATION OF REFLECTED AND STOPPED DIFFUSIONS CHRISTIAN BAYER, ANDERS SZEPESSY, AND RAÚL TEMPONE ABSTRACT. We study the weak approximation problem of diffusions, which are reflected at
More informationSome Tools From Stochastic Analysis
W H I T E Some Tools From Stochastic Analysis J. Potthoff Lehrstuhl für Mathematik V Universität Mannheim email: potthoff@math.uni-mannheim.de url: http://ls5.math.uni-mannheim.de To close the file, click
More 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 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 informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More 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 informationAn adaptive numerical scheme for fractional differential equations with explosions
An adaptive numerical scheme for fractional differential equations with explosions Johanna Garzón Departamento de Matemáticas, Universidad Nacional de Colombia Seminario de procesos estocásticos Jointly
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 informationSequential maximum likelihood estimation for reflected Ornstein-Uhlenbeck processes
Sequential maximum likelihood estimation for reflected Ornstein-Uhlenbeck processes Chihoon Lee Department of Statistics Colorado State University Fort Collins, CO 8523 Jaya P. N. Bishwal Department of
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 informationDerivation of Itô SDE and Relationship to ODE and CTMC Models
Derivation of Itô SDE and Relationship to ODE and CTMC Models Biomathematics II April 23, 2015 Linda J. S. Allen Texas Tech University TTU 1 Euler-Maruyama Method for Numerical Solution of an Itô SDE dx(t)
More 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 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 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 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 informationLie Symmetry of Ito Stochastic Differential Equation Driven by Poisson Process
American Review of Mathematics Statistics June 016, Vol. 4, No. 1, pp. 17-30 ISSN: 374-348 (Print), 374-356 (Online) Copyright The Author(s).All Rights Reserved. Published by American Research Institute
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 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 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 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 informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More 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 informationStochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno
Stochastic Processes M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno 1 Outline Stochastic (random) processes. Autocorrelation. Crosscorrelation. Spectral density function.
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 informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationA PRACTICAL SPLITTING METHOD FOR STIFF SDES WITH APPLICATIONS TO PROBLEMS WITH SMALL NOISE
A PRACTICAL SPLITTING METHOD FOR STIFF SDES WITH APPLICATIONS TO PROBLEMS WITH SMALL NOISE HECTOR D. CENICEROS AND GEORGE O. MOHLER Abstract. We present an easy to implement drift splitting numerical method
More informationStochastic Calculus Made Easy
Stochastic Calculus Made Easy Most of us know how standard Calculus works. We know how to differentiate, how to integrate etc. But stochastic calculus is a totally different beast to tackle; we are trying
More 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 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 informationEngg. Math. I. Unit-I. Differential Calculus
Dr. Satish Shukla 1 of 50 Engg. Math. I Unit-I Differential Calculus Syllabus: Limits of functions, continuous functions, uniform continuity, monotone and inverse functions. Differentiable functions, Rolle
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 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 informationExample 4.1 Let X be a random variable and f(t) a given function of time. Then. Y (t) = f(t)x. Y (t) = X sin(ωt + δ)
Chapter 4 Stochastic Processes 4. Definition In the previous chapter we studied random variables as functions on a sample space X(ω), ω Ω, without regard to how these might depend on parameters. We now
More informationAsymptotic Coupling of an SPDE, with Applications to Many-Server Queues
Asymptotic Coupling of an SPDE, with Applications to Many-Server Queues Mohammadreza Aghajani joint work with Kavita Ramanan Brown University March 2014 Mohammadreza Aghajanijoint work Asymptotic with
More informationEstimates for the density of functionals of SDE s with irregular drift
Estimates for the density of functionals of SDE s with irregular drift Arturo KOHATSU-HIGA a, Azmi MAKHLOUF a, a Ritsumeikan University and Japan Science and Technology Agency, Japan Abstract We obtain
More informationCALCULUS JIA-MING (FRANK) LIOU
CALCULUS JIA-MING (FRANK) LIOU Abstract. Contents. Power Series.. Polynomials and Formal Power Series.2. Radius of Convergence 2.3. Derivative and Antiderivative of Power Series 4.4. Power Series Expansion
More informationS. Lototsky and B.L. Rozovskii Center for Applied Mathematical Sciences University of Southern California, Los Angeles, CA
RECURSIVE MULTIPLE WIENER INTEGRAL EXPANSION FOR NONLINEAR FILTERING OF DIFFUSION PROCESSES Published in: J. A. Goldstein, N. E. Gretsky, and J. J. Uhl (editors), Stochastic Processes and Functional Analysis,
More informationFast-slow systems with chaotic noise
Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 12, 215 Averaging and homogenization workshop, Luminy. Fast-slow systems
More informationNonparametric Drift Estimation for Stochastic Differential Equations
Nonparametric Drift Estimation for Stochastic Differential Equations Gareth Roberts 1 Department of Statistics University of Warwick Brazilian Bayesian meeting, March 2010 Joint work with O. Papaspiliopoulos,
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 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 informationConvergence to equilibrium for rough differential equations
Convergence to equilibrium for rough differential equations Samy Tindel Purdue University Barcelona GSE Summer Forum 2017 Joint work with Aurélien Deya (Nancy) and Fabien Panloup (Angers) Samy T. (Purdue)
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 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 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 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 6: Multiple Model Filtering, Particle Filtering and Other Approximations
Lecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations Department of Biomedical Engineering and Computational Science Aalto University April 28, 2010 Contents 1 Multiple Model
More 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 informationAdaptive wavelet decompositions of stochastic processes and some applications
Adaptive wavelet decompositions of stochastic processes and some applications Vladas Pipiras University of North Carolina at Chapel Hill SCAM meeting, June 1, 2012 (joint work with G. Didier, P. Abry)
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 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 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 informationApplications of controlled paths
Applications of controlled paths Massimiliano Gubinelli CEREMADE Université Paris Dauphine OxPDE conference. Oxford. September 1th 212 ( 1 / 16 ) Outline I will exhibith various applications of the idea
More informationLecture on Parameter Estimation for Stochastic Differential Equations. Erik Lindström
Lecture on Parameter Estimation for Stochastic Differential Equations Erik Lindström Recap We are interested in the parameters θ in the Stochastic Integral Equations X(t) = X(0) + t 0 µ θ (s, X(s))ds +
More informationLecture 1. Finite difference and finite element methods. Partial differential equations (PDEs) Solving the heat equation numerically
Finite difference and finite element methods Lecture 1 Scope of the course Analysis and implementation of numerical methods for pricing options. Models: Black-Scholes, stochastic volatility, exponential
More informationBV functions in a Gelfand triple and the stochastic reflection problem on a convex set
BV functions in a Gelfand triple and the stochastic reflection problem on a convex set Xiangchan Zhu Joint work with Prof. Michael Röckner and Rongchan Zhu Xiangchan Zhu ( Joint work with Prof. Michael
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 informationFormulas for probability theory and linear models SF2941
Formulas for probability theory and linear models SF2941 These pages + Appendix 2 of Gut) are permitted as assistance at the exam. 11 maj 2008 Selected formulae of probability Bivariate probability Transforms
More informationStrong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term
1 Strong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term Enrico Priola Torino (Italy) Joint work with G. Da Prato, F. Flandoli and M. Röckner Stochastic Processes
More information2012 NCTS Workshop on Dynamical Systems
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,
More informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
More informationStochastic Lagrangian Transport and Generalized Relative Entropies
Stochastic Lagrangian Transport and Generalized Relative Entropies Peter Constantin Department of Mathematics, The University of Chicago 5734 S. University Avenue, Chicago, Illinois 6637 Gautam Iyer Department
More informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr
More informationClases 11-12: Integración estocástica.
Clases 11-12: Integración estocástica. Fórmula de Itô * 3 de octubre de 217 Índice 1. Introduction to Stochastic integrals 1 2. Stochastic integration 2 3. Simulation of stochastic integrals: Euler scheme
More information