1. Introduction. Much effort has been invested in studying numerical schemes for stochastic differential equations of the form
|
|
- Molly Carter
- 5 years ago
- Views:
Transcription
1 SIAM J. NUMER. ANAL. Vol. 4, No. 6, pp c 005 Society for Industrial and Applied Mathematics NUMERICAL ANALYSIS OF STOCHASTIC SCHEMES IN GEOPHYSICS BRIAN D. EWALD AND ROGER TÉMAM Abstract. We present and study the stability, convergence, and order of convergence of a numerical scheme used in geophysics, namely, the stochastic version of a deterministic implicit leapfrog scheme which has been developed for the approximation of the so-called barotropic vorticity model. Two other schemes which might be useful in the context of geophysical applications are also introduced and discussed. Key words. numerical methods, stochastic differential equations, leapfrog scheme, Adams Bashforth scheme, geophysical fluid dynamics AMS subject classifications. 65C30, 60H35, 65N1, 86A10 DOI /S Introduction. Much effort has been invested in studying numerical schemes for stochastic differential equations of the form (1.1) du t = a(u t ) dt + b(u t ) dw t, where U t R d, a is a function from R d into itself, W is a Wiener process on R m, and b is a function from R d into R d m. For the so-called weak approximation of (1.1), in which the approximation of the expectation of functions of U is considered, extensive work is due, for example, to Talay and his collaborators, work relying on probabilistic methods more involved than those used in this article (see, e.g., 1,, 1 and the references therein). The question of strong approximation of (1.1), in which the approximation of sample paths of U is desired, has also been much studied. Mil shtein, in 8, introduced the scheme (1.) U k n+1 = U k n + + d b k,j (U n ) Wn j + a k (U n ) t j=1 d j 1,j,l=1 b l,j1 (U n ) bk,j x l (U n) tn+1 t n (W j1 s W j1 t n ) dw j s, which converges to U to the order of t in mean-square error. His method involved the consideration of a functional analytic Taylor series for the infinitesimal generator of a semigroup corresponding to U and W. Rümelin later investigated a stochastic analogue of Runge Kutta (RK) schemes in 10, in which he compared them to Mil shtein s scheme. The RK schemes which he derives can be arranged to converge Received by the editors November 0, 00; accepted for publication (in revised form) September 5, 003; published electronically March 31, 005. This work was ported in part by a grant from the National Science Foundation (NSF-DMS ) and by the Research Fund of Indiana University. Department of Mathematics, Texas A&M University, College Station, TX (ewald@math. tamu.edu). Department of Mathematics, Institute for Scientific Computing and Applied Mathematics, Indiana University, Bloomington, IN (temam@indiana.edu). 57
2 58 BRIAN D. EWALD AND ROGER TÉMAM to U when the stochastic integral is interpreted in the sense of Itô, Stratonovich, or in fact for any stochastic calculus whatsoever. However, the issue of the accuracy of these RK schemes (which is not the same as in the deterministic case) is not fully addressed in 10 and is mostly unresolved. In fact, there are indications that these stochastic RK schemes are of significantly lower orders of accuracy than their deterministic counterparts (see 3). The book by Kloeden and Platen 6 and the companion volume by Kloeden, Platen, and Schurz 7 offer a systematic investigation of numerical schemes for (1.1) in both the sense of Itô and of Stratonovich, the two stochastic calculi which in applications are by far the most useful. Their methods are analytic and are applicable to proving the convergence of a wide range of numerical schemes, and they derive a very general scheme (formula (1.6.) of 6) which, for various choices of parameters, includes stochastic analogues of such deterministic schemes as the explicit and implicit Euler schemes, the Crank Nicholson scheme, and the leapfrog scheme. In the geophysics community, an enormous amount of work has been spent in developing large, complex numerical models of the oceans and atmosphere. The questions therefore arise: Is it possible to add stochastic numerical noise to these already existing models in such a way that it is known to what the scheme converges (e.g., to the Itô or Stratonovich solution of some stochastic differential equation), to what order they may be expected to converge, etc.? While we certainly do not answer these complex questions here, we consider a simple implicit leapfrog scheme for a barotropic model (plied to us by Cecile Penland and Prashant Sardeshmukh) and demonstrate one way of adding stochastic noise to it so that these questions can be answered for the resulting stochastic scheme (section 4). This scheme and the scheme in section 5 have been applied in investigating El Niño (see 4). In this paper, the schemes were used for the numerical timestepping to determine if a linear inverse model of El Niño (see 9) can be reconciled with the observed skew toward warm events in the Pacific. It was found that the observed skew is well within the range predicted by the model, although the observed trend is not. We also propose a stochastic analogue for the deterministic Adams Bashforth scheme, using methods similar to those of 6, as an attempt to produce alternate schemes which are higher order in time (studied in section 3, following the preliminary results in section ). Last, we examine the derivatives of a and b which occur naturally in the above schemes, and which can prove to be troublesome in certain applications in which these functions, especially b, are given by physical parametrizations (i.e., by tables ) and not by analytic expressions. We consider how these derivatives can be replaced by finite differences derived from space-discretization while still maintaining the existing rate of convergence (section 5). We realize that the results of this article, while very useful in our opinion, are just some small contribution to an outstanding problem, namely, the numerical analysis of stochastic differential equations which raise with more difficulty the same issues as in the deterministic case: consistency, convergence, and accuracy. All of these issues partly due to the form of the stochastic Taylor formula are considerably more difficult than in the deterministic case; in particular, consistency includes here the issue of the type of stochastic calculus (Itô, Stratonovich, or otherwise) to which the scheme converges. In the case of the geoscience scheme, the scheme that we study in section 4 is the closest we could get, at this time, to a scheme actually used in the geo-
3 STOCHASTIC SCHEMES IN GEOPHYSICS 59 sciences, without any prior information on its consistency, convergence, and accuracy properties. The version of the Adams Bashforth scheme studied in section 3 has given very good numerical results in simulations for simple (one-dimensional) stochastic differential equations; see section 6 and forthcoming articles. The numerical example in section 6 is actually based on a mistake: one of the stochastic processes involved in the scheme was mistakenly believed to be (and treated as) a Gaussian process; this did not affect the accuracy of the scheme, which remains at order two as predicted by the theory. This raises perhaps an interesting probabilistic problem about the approximation of certain non-gaussian processes used in numerical schemes. Another issue of probabilistic nature is raised by the Adams Bashforth-type scheme studied in section 3: in several or in high dimensions, a large number of stochastic processes need to be simulated, which could make the cost prohibitive. It is not excluded that future probabilistic developments will improve this situation. In particular, the first author, using some ideas of Gaines and Lyons 5, is trying at this time to develop methods of generating the needed stochastic increments. As we have said, there are, of course, a great many mathematical difficulties which this paper does not address. However, methods involving stochastic noise are already in common use in numerical simulations for the geosciences and turbulence (and, no doubt, many other areas in science). As mathematicians, we can attempt to help these scientists develop the necessary numerical tools, or watch as they do it themselves.. Preliminary results. We consider a stochastic differential equation (.1) du t = a(t, U t )dt + b(t, U t )dw t for U =(u 1,...,u d ) R d, where a : R + R d R d, b : R + R d R d m, and W is a Wiener process in R m adapted to a filtration {F t } t 0. We then have the Itô formula, which states that if F : R + R d R ˆd, then F t = F (t, U t ) satisfies the stochastic differential equation F df t = t + ak (F t ) F u k + 1 bij (F t )b kj F (.) (F t ) u i u k dt + b ij (F t ) F u i dw j t ; here we use the Einstein convention for repeated indices. We use the following notation from 6: We call a row vector α =(j 1,j,...,j l ), where each j i {0, 1,...,m}, a multi-index of length l = l(α) {1,,...}. We also use ν to denote the multi-index of length 0, i.e., l(ν) =0. We define n(α) tobethe number of entries of α which are 0. For adapted, right-continuous functions f, and stopping times ρ, τ such that 0 ρ τ T almost surely, we define f(τ) if l(α) =0, τ (.3) I α f( ) ρ,τ = ρ I α f( ) ρ,s ds if l(α) 1, j l(α) =0, τ p I α f( ) ρ,s dw j l(α) s if l(α) 1, j l(α) 0. Here α is α with its final component removed. We define the spaces H α as follows. First, H ν is the space of adapted right-continuous stochastic processes f with left limits such that f(t) is almost surely finite for each t 0. Next, H (0) contains those elements of H ν such that (.4) t 0 f(s) ds <
4 60 BRIAN D. EWALD AND ROGER TÉMAM almost surely for each t 0; and H (j) for j 0 contains those elements of H ν such that (.5) t 0 f(s) ds < almost surely for each t 0. Finally, if l(α), we define H α recursively as those elements of H ν that satisfy (.6) I α f( ) 0,t H (jl(α) ) almost surely for all t 0. We also define the operators L 0 = t + ak u k + 1 bkj b lj (.7) u k u l, (.8) L j = b kj u k, and, if f C h (R + R d, R), where h l(α)+n(α), we set { f if l(α) =0, (.9) f α = L j1 f α if l(α) 1. Here α is α with its first component removed. We note that if f(t, u) u, then f (0) = a, f (j) = b j, etc. In what follows, unless explicitly stated otherwise, we will assume that f is this identity function. A set, A, of multi-indices is said to be a hierarchical set if A, α A l(α) <, and α Awhenever α A {ν}. We then define the remainder set B(A) of A by B(A) ={α α/ A and α A}. We can now provide a stochastic Taylor expansion for U satisfying (.1): If f : R + R d R, then, provided the derivatives and integrals exist, (.10) f(τ,u τ )= I α f α (ρ, U ρ ) ρ,τ + I α f α (,U ) ρ,τ, α A α B(A) where A is some hierarchical set. Now, for γ =0.5, 1.0, 1.5,..., we set { A γ = α l(α)+n(α) γ or l(α) =n(α) =γ + 1 } (.11). We call the stochastic Taylor expansion with A = A γ the stochastic Taylor expansion to order γ. We will make use of the following lemmas in the succeeding sections. In each of them, U is the solution to (.1), and t k = k t for k =0, 1,...,N is an equipartition of 0,T, so that t N = T ; we partly rely on 6 for the proofs. Lemma.1. Suppose Y n is a stochastic process adapted to the filtration F t at the equipartition (i.e., Y n is F tn -measurable), the function f satisfies f(t, x) f(t, y) K x y for all t 0,T and x, y R, and α is a multi-index with l(α) 1. Then m 1 I α f(t k,u tk ) f(t k,y k ) tk,t k+1 (.1) C t U tm Y m. 0 m k
5 (.13) (.14) STOCHASTIC SCHEMES IN GEOPHYSICS 61 Proof. Forα = (0), we have m 1 (f(t k,u tk ) f(t k,y k )) t For α =(j), t E f(t k,u tk ) f(t k,y k ) m 1 m n t E K U tk Y k K T t E E 0 m k U tm Y m. m 1 (f(t k,u tk ) f(t k,y k )) W j k 4E (f(t k,u tk ) f(t k,y k )) W j k 4 E f(t k,u tk ) f(t k,y k ) t 4K t E 0 m k U tk Y k. For longer α s, we just repeat the above two arguments as necessary. Lemma.. Suppose the function f satisfies f(t, x) K (1 + x ) for all t 0,T and x R, and that α is a multi-index with l(α) 1. Then (.15) E m 1 { C t I α f(,u ) tk,t k+1 (l(α) 1) (1 + E U 0 ) if l(α) =n(α), C t l(α)+n(α) 1 (1 + E U 0 ) if l(α) n(α). Proof. The ideas of this proof are the same as those in the proof of Lemma.1, along with the following bound on the solution U t (see equation (4.5.16) of 6): (.16) E t 0 s T U s C(1 + E U t0 ). If we apply Lemma of 6 with g(s) =f(s, U s ), we have m 1 I α f(,u ) tk,t k+1 tn (.17) C t (l(α) 1) f(s, U s ) dt t 0 t 0 s t tn C t l(α)+n(α) 1 f(s, U s ) dt t 0 s t t 0 if l(α) =n(α), if l(α) n(α).
6 6 BRIAN D. EWALD AND ROGER TÉMAM Here, the constant C depends only on the length of the time interval T t 0 and on α. We then apply (.16) and the growth condition on f, and we have the desired result. Lemma.3. Suppose that the sequence of positive numbers Z n for n =0, 1,...,N satisfies the inequality ( ) n (.18) Z n C t Z k + t γ for some positive constant C and some γ>0. Then Z N = O( t γ ) as t 0. Proof. Set ξ n = t n Z k,soz n = 1 t (ξ n ξ ), and we have 1 (.19) t (ξ n ξ ) Cξ n + C t γ. That is, (.0) (1 C t)ξ n ξ + C t γ+1. Therefore, (1 C t) n ξ n (1 C t) ξ +(1 C t) C t γ+1, (1 C t) ξ (1 C t) n ξ n +(1 C t) n C t γ+1,. (1 C t)ξ 1 C t γ+1, and, summing, (1 C t) n ξ n (1+(1 C t)+ +(1 C t) ))C t γ+1 1 (1 C t)n (.1) 1 (1 C t) C tγ+1 (1 (1 C t) n ) t γ. Since (1 C t) N e CT as N = T/ t, we see that ξ n C t γ for some (different) C. Thus, by (.18), (.) Z n C(ξ n + t γ ) C t γ. 3. A stochastic Adams Bashforth scheme. The deterministic Adams Bashforth scheme for the ordinary differential equation φ = F (φ) takes the form (3.1) φ n+1 = φ n + t 3F (φ n) F (φ ). This scheme is order t. We will derive a stochastic version of this scheme which maintains the same order. We begin with the stochastic Taylor expansion to order γ =.0: (3.) U t+ = U t + b j W j + a +L j1 b j I (j1,j ) + L 0 b j I (0,j) + L j ai (j,0) + L j1 L j b j3 I (j1,j,j 3) + 1 L0 a + L 0 L j1 b j I (0,j1,j ) + L j1 L 0 b j I (j1,0,j ) + L j1 L j ai (j1,j,0) + L j1 L j L j3 b j4 I (j1,j,j 3,j 4) + R.0(t) = U t + a + 1 L0 a + M (t),
7 STOCHASTIC SCHEMES IN GEOPHYSICS 63 where each coefficient is at the point (t, U t ), and each stochastic integral is from t to t +, = t. We have also used the Einstein summation convention. Similarly, for γ =1.5, we have (3.3) a(t +,U t+ )=a + L 0 a +N (t), where N (t) = 1 L0 L 0 a + L j a W j + L 0 L j ai (0,j) + L j L 0 ai (j,0) + L j1 L j ai (j1,j ) + L j1 L j L j3 ai (j1,j,j 3) + R 1.5(t), and, for γ =1.0, (3.4) L 0 a(t +,U t+ )=L 0 a + P (t), where P (t) =L 0 L 0 a +L j L 0 a W j + L j1 L j L 0 ai (j1,j ) + R 1.0(t). Combining these results, we get (3.5) U t+ = U t +αa(t +,U t+ )+(1 α)a ( ) 1 + α βl 0 a(t +,U t+ )+(1 β)l 0 a ( ) 1 α N (t) α β P (t)+m (t). In particular, if t = t n, = t, α =0,β =0, and writing U n = U tn, (3.6) U n+ = U n +a(t n,u n ) t +L 0 a(t n,u n ) t + M t (t n ), and if t = t n, = t, α = 3,β =0, (3.7) U n+1 = U n 3 a(t n+1,u n+1 ) t + 5 a(t n,u n ) t +L 0 a(t n,u n ) t + 3 N t (t n ) t + M t (t n ). Therefore, (3.8) U n+ = U n+1 +(U n+ U n ) (U n+1 U n ) 3 = U n+1 + a(t n+1,u n+1 ) 1 a(t n,u n ) t 3 tn t (t n )+(M t (t n ) M t (t n )). This leads us to consider the following stochastic Adams Bashforth (SAB) scheme: 3 Y n+ = Y n+1 + a(t n+1,y n+1 ) 1 a(t n,y n ) t (3.9) 3 ta n(t n,y n )+B n (t n,y n ), in which (3.10) A n (t, x) =L j a(t, x) W j + L j1 L j a(t, x)i (j1,j ),
8 64 BRIAN D. EWALD AND ROGER TÉMAM where the random intervals are from time t n to t n+1, and B n (t, x) =b j (t, x) W j + L 0 b j (t, x)i (0,j) + L j a(t, x)i (j,0) (3.11) + L j1 b j (t, x)i (j1,j ) + L 0 L j1 b j (t, x)i (0,j1,j ) + L j1 L 0 b j (t, x)i (j1,0,j ) + L j1 L j a(t, x)i (j1,j,0) + L j1 L j b j3 (t, x)i (j1,j,j 3) + L j1 L j L j3 b j4 (t, x)i (j1,j,j 3,j 4), where the random intervals are those from time t n to t n+ minus those from time t n to t n+1. We then have the following theorem. Theorem 3.1. Suppose that the coefficient functions f α satisfy (3.1) f α (t, x) f α (t, y) K x y for all α A.0, t 0,T, and x, y R d ; (3.13) f α C 1, and f α H α for all α A.0 B(A.0 ); and (3.14) f α (t, x) K(1 + x ) for all α A.0 B(A.0 ), t 0,T, and x R d. Choose t 1 and set N = T/ t, and define t n = n t for n =1,...,N. Suppose that Y 0 is some (nonrandom) initial condition and that some appropriate numerical scheme is used to generate Y 1 such that E U t1 Y 1 F 0 1 C t. Then (3.15) 1 U tn Y n F 0 C t. 0 n N Proof. First, we note that (3.16) 3 U n+ = U n+1 + a(t n+1u n+1 ) 1 a(t n,u n ) t 3 ta n(t n,u n )+B n (t n,u n )+R n, where R n = 3 t 1 L0 L 0 a(t n,u n ) t + L 0 L j a(t n,u n )I (0,j) + L j L 0 a(t n,u n )I (j,0) + L j1 L j L j3 a(t n,u n )I (j1,j,j 3) + R t 1.5(t n ) + R t.0 (t n ) R t.0(t n ).
9 (3.17) STOCHASTIC SCHEMES IN GEOPHYSICS 65 If we iterate (3.9) and (3.16), we arrive at n 3 U n = U 1 + t a(t k+1,u k+1 ) 1 a(t k,u k ) 3 n t n A k (t k,u k )+ n B k (t k,u k )+ n 3 Y n = Y 1 + t a(t k+1,y k+1 ) 1 a(t k,y k ) Set ζ n = U n Y n. Then 3 n t n A k (t k,y k )+ B k (t k,y k ). n 3 ζ n = ζ 1 + t (a(t k+1,u k+1 ) a(t k+1,y k+1 )) 1 (a(t k,u k ) a(t k,y k )) (3.18) 3 n t (A k (t k,u k ) A k (t k,y k )) n n + (B k (t k,u k ) B k (t k,y k )) + R k. Set Z n = E ζ m F 0. Then we have ( m Z n C E ζ 1 F 0 + t a(t k+1,u k+1 ) a(t k+1,y k+1 ) F 0 m + t a(t k,u k ) a(t k,y k ) F 0 + t m (L j a(t k,u k ) L j a(t k,y k )) W j k F 0 Terms m + t from (L j1 L j A k a(t k,u k ) L j1 L j a(t k,y k ))I (j1,j ) F 0 + m I αɛa α f α (t k,u k ) f α (t k,y k ) tk,t k+ F 0 Terms from + E αɛa + R k, m I α f α (t k,u k ) f α (t k,y k ) tk,t k+1 F 0 m R k F 0 ). B k
10 66 BRIAN D. EWALD AND ROGER TÉMAM (3.19) We go term-by-term: m t a(t k+1,u k+1 ) a(t k+1,y k+1 ) F 0 m t n E a(t k+1,u k+1 ) a(t k+1,y k+1 ) F 0 T t E Similarly, (3.0) t E n a(t k+1,u k+1 ) a(t k+1,y k+1 ) F 0 n C t E U k+1 Y k+1 F 0 C t Z k. m a(t k,u k ) a(t k,y k ) F 0 C t Next, we consider the terms from A k. From Lemma.1 with α =(j), m t (L j a(t k,u k ) L j a(t k,y k )) W j t k,t k+1 F 0 (3.1) n C t Z k, and from Lemma.1 with α =(j 1,j ), m t (L j1 L j a(t k,u k ) L j1 L j a(t k,y k ))I (j1,j )t k,t k+1 F 0 (3.) C t Z k. Now, we consider the terms from B k.forα A (i.e., α A,l(α) n(α)), m (3.3) I α f α (t k,u k ) f α (t k,y k ) tk,t k+1 F 0 C t Z k. The other terms from B k are similar. This leaves only the terms from R k. From Lemma. with α = (0), m t L 0 L 0 a(t n,u n ) t F 0 n t 4 (3.4) L 0 L 0 a(t n,u n ) F 0 C t 4 (1 + U 0 ). Z k.
11 STOCHASTIC SCHEMES IN GEOPHYSICS 67 For α =(0,j), (j, 0), (j 1,j,j 3 ), by Lemma. m t I α a α (t k,u k ) tk,t k+1 F 0 (3.5) C t 4 (1 + U 0 ), since l(α)+n(α) =3. If α B(Aγ) (here γ =1.5 or.0), we have by Lemma., m (3.6) I α α(,u ) tk,t k+1 F 0 C(1 + U 0 ) t γ. Therefore, we have (3.7) m R k F 0 C(1 + U 0 ) t 4. So, overall, we see that (3.8) Z n C Z 1 +(1+ U 0 ) t 4 + t Z k. (3.9) and The result then follows from Lemma.3. Remark 3.1. If we truncate A n and B n to A n (t, x) =L j a(t, x) W j (3.30) B n (t, x) =b j (t, x) W j + L 0 b j (t, x)i (0,j) + L j a(t, x)i (j,0) + L j1 b j (t, x)i (j1,j ) + L j1 L j b j3 (t, x)i (j1,j,j 3), the same proof will show that the convergence is now to order t 3. We note that although the order t SAB scheme seems to have no obvious advantages over the standard t strong one-step explicit scheme (as in 6), the order t 3 SAB scheme does have an advantage over the order t 3 strong one-step explicit scheme in that the former lacks the terms involving the second derivative of a which are present in the latter. Remark 3.. It can be shown that the scheme 3 Y n+ = Y n+1 + a(t n+1,y n+1 ) 1 a(t n,y n ) t (3.31) 3 ta n(t n,y n )+B n (t n,y n ), in which (3.3) A n (t, x) =L j a(t, x) W j, where the random intervals are from time t n to t n+1, and (3.33) B n (t, x) =b j (t, x) W j + L 0 b j (t, x)i (0,j) + L j a(t, x)i (j,0) + L j1 b j (t, x)i (j1,j ),
12 68 BRIAN D. EWALD AND ROGER TÉMAM where the random intervals are those from time t n to t n+ minus those from time t n to t n+1, converges to the Itô solution in the weak sense to order. As can be seen, this scheme is considerably simpler than the strong scheme, and it avoids the difficulties with generating the higher-order moments that the strong scheme has. 4. A stochastic implicit leapfrog scheme. The barotropic vorticity model plied to us by Cecile Penland and Prashant Sardeshmukh of the National Oceanic and Atmospheric Administration in Boulder, Colorado (see 11), takes the form ζ (4.1) t = (vζ)+s rξ κ 4 ξ, where ζ = ψ + f = ξ + f and v = ˆk ψ. Here, ζ is the total vorticity, v is the velocity vector, f is the Coriolis term, S is a (deterministic) forcing, r and κ are constants, and ξ is the local vorticity. The numerical scheme they provided for this uses spherical harmonics, and, writing F for (vζ), the equation becomes (4.) d dt ζm n = F m n + S m n rξ m n κ n(n +1) a ζ m n. Then the scheme has two steps. First, a leapfrog step, (4.3) ζ m n (t + t) =ζ m n (t t)+ tf m n (t)+s m n (t), followed by an implicit step, (4.4) ζ ζn m n m (t + t) (t + t) =. 1+ t r + κ n(n+1) a If we simplify notation and write a 1 for F + S and a for rξ κ 4 ξ, we see that this is just an implicit leapfrog scheme {Ỹ (t + t) =Y (t t)+ ta1 (t, Y (t)), (4.5) Y (t + t) =Ỹ (t + t)+ ta (t + t, Y (t + t)) for the equation (4.6) du(t) =a 1 (t, U(t)) + a (t, U(t)) dt. Therefore, we consider a stochastic differential equation of the form (4.7) du t =(a 1 (t, U t )+a (t, U t )) dt + b(t, U t ) dw t. Note that we have simply added a general diffusion term to the deterministic differential equation (4.6). We will consider the scheme {Ỹn+ = Y n +a 1 (t n+1,y n+1 ) t + M n (Y n )+M n+1 (Y n+1 ), (4.8) Y n+ = Ỹn+ +a (t n+,y n+ ) t, where (4.9) M n (y) =b(t n,y) W n + bb (t n,y)i (1,1),n.
13 STOCHASTIC SCHEMES IN GEOPHYSICS 69 Theorem 4.1. Suppose that the coefficient functions f α satisfy (4.10) f α (t, x) f α (t, y) K x y for all α A 1.0, t 0,T, and x, y R d ; (4.11) f α C 1, and f α H α for all α A 1.0 B(A 1.0 ); and (4.1) f α (t, x) K(1 + x ) for all α A 1.0 B(A 1.0 ), t 0,T, and x R d. Choose t 1 and set N = T/ t, and define t n = n t for n =1,...,N. Suppose that some appropriate numerical scheme is used to generate Y 1 such that E U t1 Y 1 F 0 1 C t. Then 1 U tn Y n (4.13) F 0 C t. 0 n N Proof. We note first that, by Itô s formula (i.e., the Taylor expansion with γ = 0.0), the solution U to (4.7) satisfies the following equations (where for notational simplicity we have written U n for U tn ): a 1 (t n+1,u n+1 )=a 1 (t n,u n )+R t,a1 0.0 (t n ); (4.14) a (t n+,u n+ )=a (t n,u n )+R t,a 0.0 (t n ); a (t n+,u n+ )=a (t n+1,u n+1 )+R t,a 0.0 (t n+1 ). Therefore we have U n+ = U n +(U n+ U n+1 )+(U n+1 U n ) = U n +b(t n+1,u n+1 ) W n+1 + a 1 (t n+1,u n+1 ) t (4.15) + a (t n+1,u n+1 ) t + bb (t n+1,u n+1 )I (1,1),n+1 + R1.0(t t n+1 )+b(t n,u n ) W n + a 1 (t n,u n ) t + a (t n,u n ) t + bb (t n,u n )I (1,1),n + R1.0(t t n ). After substituting (4.14) into this, we see that U n+ =U n +a 1 (t n+1,u n+1 ) t + M n (U n )+M n+1 (U n+1 ) (4.16) +a (t n+,u n+ ) t + R n ; here R n = R1.0(t t n )+R1.0(t t n+1 ) tr t,a1 0.0 (t n ) (4.17) + R t,a 0.0 (t n )+R t,a 0.0 (t n+1 ). If we iterate (4.16), we arrive at (4.18) U n = U n + t + n/ k=1 k=n M n (U n )+ n/ a 1 (t k 1+n,U k 1+n )+ a (t k+n,u k+n ) n/ here n is0ifn is even and 1 if n is odd. R k +n ; k=1
14 70 BRIAN D. EWALD AND ROGER TÉMAM (4.19) Similarly, we have for Y n/ n/ Y n = Y n + t a 1 (t k 1+n,Y k 1+n )+ a (t k+n,y k+n ) + k=1 k=n M n (Y n ). Let us set Z n = E U m Y m F 0. Then, by subtracting (4.19) from (4.18) and then squaring and taking expectations, we find (4.0) In (4.0), Z n CE + Z m + t A 1,m + t A,m ( m 1 k=m M k (U k ) M k (Y k ) ) ( m/ + k=1 R k +m ) F 0. (4.1) A 1,n = n/ k=1 a 1 (t k 1+n,U k 1+n ) a 1 (t k 1+n,Y k 1+n ) and (4.) A,n = n/ k=1 a (t k+n,u k+n ) a (t k+n,y k+n ). We then have the following estimates (where we omit the dependence on t when it is clear): t A 1,m (4.3) m m/ t a 1 (U k 1+n ) a 1 (Y k 1+n ) k=1 n n K t ne (U k Y k ) KT t Z k. k=1 k=1 (4.4) Similarly, t A,n KT t n Z k. k=1 There are two terms in M k. For the first one, from Lemma.1, with α = (1), we obtain m 1 n (4.5) b(u k ) b(y k ) W k C t Z k. k=m The second term is similar, with α =(1, 1).
15 STOCHASTIC SCHEMES IN GEOPHYSICS 71 Finally, we show a representative term from the remainder R. From Lemma., with α =(1, 0), m 1 tk+1 s L 0 (4.6) b(u r ) dr dw s C t (1 + U0 ). t k t k The remaining terms from R give similar bounds. Taking all of these estimates into account, we have the inequality n (4.7) Z n K Z 1 + t Z k + t. Therefore, from Lemma.3, we see that 1 (4.8) U n Y n F 0 = O( t). o n N Remark 4.1. It is possible to show that the scheme Ỹ n+ = Y n +a 1 (t n+1,y n+1 ) t + b(t n,y n ) W n (4.9) + b(t n+1,y n+1 ) W n+1, Y n+ = Ỹn+ +a (t n+,y n+ ) t converges to the Itô solution in the weak sense to order 1. Again, it can be seen that the weak scheme is simpler than the strong scheme to the same order. However, we also note that the weak order 1 Stratonovich scheme Ỹ n+ = Y n +a 1 (t n+1,y n+1 ) t + bb (t n+1,y n+1 )( W n+1 ) (4.30) + b(t n,y n ) W n + b(t n+1,y n+1 ) W n+1, Y n+ = Ỹn+ +a (t n+,y n+ ) t is not appreciably simpler than the strong order 1 Stratonovich scheme {Ỹn+ = Y n +a 1 (t n+1,y n+1 ) t + M n (Y n )+M n+1 (Y n+1 ), (4.31) where Y n+ = Ỹn+ +a (t n+,y n+ ) t, (4.3) M n (y) =b(t n,y) W n + 1 bb (t n,y)( W n ). 5. Discretization of spatial derivatives by finite differences. It sometimes happens in applications that the functions a and b may only be known empirically (i.e., in tables) rather than analytically. In such cases, analytic derivatives of these functions can be difficult to obtain. It is therefore useful to replace these derivatives by discrete approximations. As a first example, consider this modification of Mil shtein s scheme: d Ŷn+1 k = Ŷ n k + b k,j (Ŷn) Wn j + a k (Ŷn) t (5.1) + j=1 d j 1,j,l=1 1 x bl,j1 (Ŷn)(b k,j (Ŷn + xe l ) b k,j (Ŷn))I (j1,j ),n,
16 7 BRIAN D. EWALD AND ROGER TÉMAM where e l is the vector (0,...,0, 1, 0,...,0), with 1 in the lth position, and we have chosen x >0. We have also pressed the dependence of a and b on time to simplify notation. We then have the following theorem. Theorem 5.1. Suppose that a and b have the regularity required for Mil shtein s scheme to converge to the solution U to order t. Then 1 U n Ŷn 1 (5.) = O(max{ t, x t }). 0 n N Note that if we want to maintain the order of convergence of Mil shtein s scheme, we need that x = O( t 1 ). Proof. We denote Mil shtein s scheme by Y and recall that it satisfies 1 Y n U n (5.3) = O( t). 0 n N First, we see that (using Einstein s summation convention on repeated indices) Ŷn+1 k Yn+1 k = Ŷ n k Yn k +(b k,j (Ŷn) b k,j (Y n )) Wn j +(a k (Ŷn) a k (Y n )) t (5.4) + b l,j1 (Ŷn) 1 x (bk,j (Ŷn + xe l ) b k,j (Ŷn)) b l,j1 (Y n ) bk,j x l (Y n) I (j1,j ),n. Iterating this, we have (5.5) Ŷn k Yn k = Ŷ 0 k Y0 k + (b k,j (Ŷi) b k,j (Y i )) W j i + (a k (Ŷi) a k (Y i )) t + b l,j1 (Ŷi) 1 x (bk,j (Ŷi + xe l ) b k,j (Ŷi)) b l,j1 (Y i ) bk,j x l (Y i) I (j1,j ),i. Set (5.6) Z n = Ŷn Y n. We then have the estimates m 1 E (b k,j (Ŷi) b k,j (Y i )) W j i 4E (b k,j (Ŷi) b k,j (Y i )) W j i 4E b k,j (Ŷi) b k,j (Y i ) t C t Z i,
17 STOCHASTIC SCHEMES IN GEOPHYSICS 73 and E m 1 (a k (Ŷi) a k (Y i )) t t ne a k (Ŷi) a k (Y i ) C t Z i, and, finally, E m 1 b l,j1 (Ŷi) 1 x bk,j (Ŷi + xe l ) b k,j (Ŷi) b l,j1 (Y i ) bk,j x l (Y i) I (j1,j ),i C t E b l,j1 (Ŷi) 1 x bk,j (Ŷi + xe l ) b k,j (Ŷi) b l,j1 (Y i ) bk,j x l (Y i) C t E b l,j1 (Ŷi) bk,j x l (Ŷi) b l,j1 (Y i ) bk,j x l (Y i) 1 + C t b l,j1 (Ŷi) x (bk,j (Ŷi + xe l ) b k,j (Ŷi)) bk,j x l C t Z i + C t x. Therefore, altogether we have (Ŷi) (5.7) Z n C t Z k + C t x, and an application of Lemma.3 implies that (5.8) 1 Y n Ŷn 1 = O( x t ). 0 n N We can apply a similar idea to the SAB scheme. That is, if we replace A n and B n in the order t 3 scheme (see Remark 3.1) by, for instance (with, again, Einstein s summation convention in effect), (5.9) A n (t, x) =b kj (t, x) 1 x a(t, x + xek ) a(t, x) W j, where e k =(0,...,0, 1, 0,...,0), with the 1 in the kth position, and
18 74 BRIAN D. EWALD AND ROGER TÉMAM B n (t, x) =b j (t, x) W j + 1 t bj (t + t, x) b j (t, x)i (0,j),n + a k (t, x) 1 x bj (t, x + xe k ) b j (t, x)i (0,j),n + b ki (t, x)b li 1 (t, x) 8 x bj (t, x + x(e i + e l )) b j (t, x + x(e i e l )) b j (t, x + x(e l e i )) + b j (t, x x(e i + e l ))I (0,j),n 1 + b kj1 (t, x) x bj (t, x + xe k ) b j (t, x xe k )I (j1,j ),n + b kj (t, x) 1 x a(t, x + xek ) a(t, x)i (j,0),n 1 + b k1j1 (t, x) 4 x bkj (t, x + xe k1 ) b kj (t x + xe k1 ) b j3 (t, x + xe k ) b j3 (t, x xe k1 )I (j1,j,j 3),n 1 + b k1j1 (t, x)b k,j (t, x) 4 x b j3 (t, x + x(e k1 + e k )) b j3 (t, x + x(e k1 e k )) b j3 (t, x + x(e k e k1 )) + b j3 (t, x x(e k1 + e k ))I (j1,j,j 3),n, we could then prove that this scheme converges to order max{ t 3, t 1 x} in a similar fashion. 6. Numerical simulation. The object of this section is to test numerically the accuracy of the scheme of section 3 and compare it to the theoretical result above (i.e., O( t ) accuracy) and to the accuracy of the Euler and Mil shtein schemes (respectively, O( t 1 ) and O( t)). All the numerical results below are consistent with the theoretical ones. We consider the following equation: (6.1) dx t = β sinh X t cosh X t dt + β cosh X t dw t, with β = 1 10 and X 0 = 1. This has the exact solution (6.) X t = arctanh(βw t + tanh X 0 ), respectively. This can be easily verified using Itô s formula and is just one of many possible examples listed in 6. We computed approximate solutions Y n using the Euler and Mil shtein schemes and the SAB scheme from section 3. Then we computed the following error: ( (6.3) e = X n Y n ). 0 n N To estimate the mean value needed, we used 500 sample trajectories. In Figure 6.1, the order of each scheme is given by the slope of the corresponding line. So we can see that the orders are 1 for Euler, 1 for Mil shtein, and for the SAB of section 3.
19 STOCHASTIC SCHEMES IN GEOPHYSICS slope = 1/ 10 slope = 1 log(error) 1 14 slope = 16 SAB Euler Milshtein log(time step) Fig Results obtained with the stochastic equation (6.1). Note that for the SAB scheme, the stochastic integral I (0,1,1) (which is difficult to generate) was approximated by a normal law. The results tend to show that this does not affect the accuracy (at least in these two cases). We will try to improve this point, which seems to raise interesting probabilistic questions, as already mentioned in the introduction. Acknowledgments. The authors are very grateful to Cecile Penland for bringing these issues to their attention, and they acknowledge very useful discussions with her and with Prashant Sardeshmukh. They are also indebted to Arnaud Debussche for several improvements on an earlier draft, and to Sylvain Faure, who provided the numerical simulations of section 6. REFERENCES 1 V. Bally and D. Talay, The law of the Euler scheme for stochastic differential equations, I. Convergence rate of the distribution function, Probab. Theory Related Fields, 104 (1996), pp V. Bally and D. Talay, The law of the Euler scheme for stochastic differential equations, II. Convergence rate of the density, Monte Carlo Methods Appl., (1996), pp B. D. Ewald, Numerical Methods for Stochastic Differential Equations in the Geosciences, Ph.D. dissertation, Indiana University, Bloomington, B. D. Ewald, C. Penland, and R. Temam, Accurate integration of stochastic climate models with application to El Niño, Monthly Weather Rev., 13 (004), pp J. G. Gaines and T. J. Lyons, Random generation of stochastic area integrals, SIAM J. Appl. Math., 54 (1994), pp P. E. Kloeden and E. Platen, Numerical Solution of Stochastic Differential Equations, Appl. Math. 3, Springer-Verlag, Berlin, 199.
20 76 BRIAN D. EWALD AND ROGER TÉMAM 7 P. E. Kloeden, E. Platen, and H. Schurz, Numerical Solution of SDE through Computer Experiments, Universitext, Springer-Verlag, Berlin, G. N. Mil shtein, Approximate integration of stochastic differential equations, Theory Probab. Appl., 19 (1974), pp C. Penland and P. D. Sardeshmukh, The optimal growth of tropical sea surface temperatures, J. Climate, 8 (1995), pp W. Rümelin, Numerical treatment of stochastic differential equations, SIAM J. Numer. Anal., 19 (198), pp P. D. Sardeshmukh and B. J. Hoskins, The generation of global rotational flow by steady idealized tropical divergence, J. Atmospheric Sci., 45 (1988), pp D. Talay, Simulation of stochastic differential equations, in Probabilistic Methods in Applied Physics, P. Krée and W. Wedig, eds., Springer-Verlag, Berlin, 1995, pp
WEAK 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 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 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 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 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 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 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 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 informationOn the Strong Approximation of Jump-Diffusion Processes
QUANTITATIV FINANC RSARCH CNTR QUANTITATIV FINANC RSARCH CNTR Research Paper 157 April 25 On the Strong Approximation of Jump-Diffusion Processes Nicola Bruti-Liberati and ckhard Platen ISSN 1441-81 www.qfrc.uts.edu.au
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 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 informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More 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 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 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 informationc 2012 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 50, No. 4, pp. 849 860 c 0 Society for Industrial and Applied Mathematics TWO RESULTS CONCERNING THE STABILITY OF STAGGERED MULTISTEP METHODS MICHELLE GHRIST AND BENGT FORNBERG
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
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 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 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 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 informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More informationRegularity of the density for the stochastic heat equation
Regularity of the density for the stochastic heat equation Carl Mueller 1 Department of Mathematics University of Rochester Rochester, NY 15627 USA email: cmlr@math.rochester.edu David Nualart 2 Department
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationNumerical Solutions for Stochastic Differential Equations and Some Examples
Brigham Young University BYU ScholarsArchive All Theses and Dissertations 2009-07-06 Numerical Solutions for Stochastic Differential Equations and Some Examples Yi Luo Brigham Young University - Provo
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 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 informationReview Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester
HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE
More informationA scaling limit from Euler to Navier-Stokes equations with random perturbation
A scaling limit from Euler to Navier-Stokes equations with random perturbation Franco Flandoli, Scuola Normale Superiore of Pisa Newton Institute, October 208 Newton Institute, October 208 / Subject of
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 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 informationApproximation of BSDEs using least-squares regression and Malliavin weights
Approximation of BSDEs using least-squares regression and Malliavin weights Plamen Turkedjiev (turkedji@math.hu-berlin.de) 3rd July, 2012 Joint work with Prof. Emmanuel Gobet (E cole Polytechnique) Plamen
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 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 informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationYURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL
Journal of Comput. & Applied Mathematics 139(2001), 197 213 DIRECT APPROACH TO CALCULUS OF VARIATIONS VIA NEWTON-RAPHSON METHOD YURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL Abstract. Consider m functions
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationOn pathwise stochastic integration
On pathwise stochastic integration Rafa l Marcin Lochowski Afican Institute for Mathematical Sciences, Warsaw School of Economics UWC seminar Rafa l Marcin Lochowski (AIMS, WSE) On pathwise stochastic
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 informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
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 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 informationMulti-Factor Finite Differences
February 17, 2017 Aims and outline Finite differences for more than one direction The θ-method, explicit, implicit, Crank-Nicolson Iterative solution of discretised equations Alternating directions implicit
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 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 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 methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
More informationItô s formula. Samy Tindel. Purdue University. Probability Theory 2 - MA 539
Itô s formula Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Itô s formula Probability Theory
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 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 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 information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 20
2.29 Numerical Fluid Mechanics Fall 2011 Lecture 20 REVIEW Lecture 19: Finite Volume Methods Review: Basic elements of a FV scheme and steps to step-up a FV scheme One Dimensional examples d x j x j 1/2
More informationA Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit
A Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit David Vener Department of Mathematics, MIT May 5, 3 Introduction In 8.366, we discussed the relationship
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
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 informationMath 660 Lecture 4: FDM for evolutionary equations: ODE solvers
Math 660 Lecture 4: FDM for evolutionary equations: ODE solvers Consider the ODE u (t) = f(t, u(t)), u(0) = u 0, where u could be a vector valued function. Any ODE can be reduced to a first order system,
More informationEuler s Method, Taylor Series Method, Runge Kutta Methods, Multi-Step Methods and Stability.
Euler s Method, Taylor Series Method, Runge Kutta Methods, Multi-Step Methods and Stability. REVIEW: We start with the differential equation dy(t) dt = f (t, y(t)) (1.1) y(0) = y 0 This equation can be
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 informationIndependence of some multiple Poisson stochastic integrals with variable-sign kernels
Independence of some multiple Poisson stochastic integrals with variable-sign kernels Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological
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 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 informationAn Uncertain Control Model with Application to. Production-Inventory System
An Uncertain Control Model with Application to Production-Inventory System Kai Yao 1, Zhongfeng Qin 2 1 Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China 2 School of Economics
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 informationApproximation of random dynamical systems with discrete time by stochastic differential equations: I. Theory
Random Operators / Stochastic Eqs. 15 7, 5 c de Gruyter 7 DOI 1.1515 / ROSE.7.13 Approximation of random dynamical systems with discrete time by stochastic differential equations: I. Theory Yuri A. Godin
More informationNumerische Mathematik
Numer. Math. 2007 07:669 705 DOI 0.007/s002-007-004-z Numerische Mathematik A fourth-order numerical method for the planetary geostrophic equations with inviscid geostrophic balance Roger Samelson Roger
More informationAccurate Integration of Stochastic Climate Models with Application to El Niño
154 MONTHLY WEATHER REVIEW Accurate Integration of Stochastic Climate Models with Application to El Niño BRIAN EWALD * Institute for Scientific Computing and Applied Mathematics, Department of Mathematics,
More informationQuasi-invariant measures on the path space of a diffusion
Quasi-invariant measures on the path space of a diffusion Denis Bell 1 Department of Mathematics, University of North Florida, 4567 St. Johns Bluff Road South, Jacksonville, FL 32224, U. S. A. email: dbell@unf.edu,
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 informationWeighted Sums of Orthogonal Polynomials Related to Birth-Death Processes with Killing
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 2, pp. 401 412 (2013) http://campus.mst.edu/adsa Weighted Sums of Orthogonal Polynomials Related to Birth-Death Processes
More informationRounding Error in Numerical Solution of Stochastic Differential Equations
STOCHASTIC ANALYSIS AND APPLICATIONS Vol. 21, No. 2, pp. 281 300, 2003 Rounding Error in Numerical Solution of Stochastic Differential Equations Armando Arciniega* and Edward Allen Department of Mathematics
More informationOrdinary Differential Equations (ODEs)
Ordinary Differential Equations (ODEs) 1 Computer Simulations Why is computation becoming so important in physics? One reason is that most of our analytical tools such as differential calculus are best
More informationmultistep methods Last modified: November 28, 2017 Recall that we are interested in the numerical solution of the initial value problem (IVP):
MATH 351 Fall 217 multistep methods http://www.phys.uconn.edu/ rozman/courses/m351_17f/ Last modified: November 28, 217 Recall that we are interested in the numerical solution of the initial value problem
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 informationCOMPUTER-AIDED MODELING AND SIMULATION OF ELECTRICAL CIRCUITS WITH α-stable NOISE
APPLICATIONES MATHEMATICAE 23,1(1995), pp. 83 93 A. WERON(Wroc law) COMPUTER-AIDED MODELING AND SIMULATION OF ELECTRICAL CIRCUITS WITH α-stable NOISE Abstract. The aim of this paper is to demonstrate how
More informationChapter 2 Event-Triggered Sampling
Chapter Event-Triggered Sampling In this chapter, some general ideas and basic results on event-triggered sampling are introduced. The process considered is described by a first-order stochastic differential
More informationApproximation of Top Lyapunov Exponent of Stochastic Delayed Turning Model Using Fokker-Planck Approach
Approximation of Top Lyapunov Exponent of Stochastic Delayed Turning Model Using Fokker-Planck Approach Henrik T. Sykora, Walter V. Wedig, Daniel Bachrathy and Gabor Stepan Department of Applied Mechanics,
More informationOptimal Prediction for Radiative Transfer: A New Perspective on Moment Closure
Optimal Prediction for Radiative Transfer: A New Perspective on Moment Closure Benjamin Seibold MIT Applied Mathematics Mar 02 nd, 2009 Collaborator Martin Frank (TU Kaiserslautern) Partial Support NSF
More informationDissipative quasi-geostrophic equations with L p data
Electronic Journal of Differential Equations, Vol. (), No. 56, pp. 3. ISSN: 7-669. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) Dissipative quasi-geostrophic
More informationMultistep Methods for IVPs. t 0 < t < T
Multistep Methods for IVPs We are still considering the IVP dy dt = f(t,y) t 0 < t < T y(t 0 ) = y 0 So far we have looked at Euler s method, which was a first order method and Runge Kutta (RK) methods
More informationCBMS Lecture Series. Recent Advances in. the Numerical Approximation of Stochastic Partial Differential Equations
CBMS Lecture Series Recent Advances in the Numerical Approximation of Stochastic Partial Differential Equations or more accurately Taylor Approximations of Stochastic Partial Differential Equations A.
More informationt f(u)g (u) g(u)f (u) du,
Chapter 2 Notation. Recall that F(R 3 ) denotes the set of all differentiable real-valued functions f : R 3 R and V(R 3 ) denotes the set of all differentiable vector fields on R 3. 2.1 7. The problem
More informationRough paths methods 4: Application to fbm
Rough paths methods 4: Application to fbm Samy Tindel Purdue University University of Aarhus 2016 Samy T. (Purdue) Rough Paths 4 Aarhus 2016 1 / 67 Outline 1 Main result 2 Construction of the Levy area:
More informationSeries Solutions. 8.1 Taylor Polynomials
8 Series Solutions 8.1 Taylor Polynomials Polynomial functions, as we have seen, are well behaved. They are continuous everywhere, and have continuous derivatives of all orders everywhere. It also turns
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 informationON THE UNIQUENESS OF SOLUTIONS TO THE GROSS-PITAEVSKII HIERARCHY. 1. Introduction
ON THE UNIQUENESS OF SOLUTIONS TO THE GROSS-PITAEVSKII HIERARCHY SERGIU KLAINERMAN AND MATEI MACHEDON Abstract. The purpose of this note is to give a new proof of uniqueness of the Gross- Pitaevskii hierarchy,
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 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 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 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 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 informationThis ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0
Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x
More informationJim Lambers MAT 772 Fall Semester Lecture 21 Notes
Jim Lambers MAT 772 Fall Semester 21-11 Lecture 21 Notes These notes correspond to Sections 12.6, 12.7 and 12.8 in the text. Multistep Methods All of the numerical methods that we have developed for solving
More informationRunge Kutta methods for numerical solution of stochastic dierential equations
Journal of Computational and Applied Mathematics 18 00 19 41 www.elsevier.com/locate/cam Runge Kutta methods for numerical solution of stochastic dierential equations A.Tocino a; ; 1, R.Ardanuy b a Departamento
More informationChapter 5 Exercises. (a) Determine the best possible Lipschitz constant for this function over 2 u <. u (t) = log(u(t)), u(0) = 2.
Chapter 5 Exercises From: Finite Difference Methods for Ordinary and Partial Differential Equations by R. J. LeVeque, SIAM, 2007. http://www.amath.washington.edu/ rjl/fdmbook Exercise 5. (Uniqueness for
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 information