期刊
SIAM JOURNAL ON SCIENTIFIC COMPUTING
卷 29, 期 1, 页码 245-264出版社
SIAM PUBLICATIONS
DOI: 10.1137/050646032
关键词
damped harmonic oscillators with noise; stationary distribution; stochastic Runge-Kutta methods; implicit midpoint rule; multiplicative noise
We seek numerical methods for second-order stochastic differential equations that reproduce the stationary density accurately for all values of damping. A complete analysis is possible for scalar linear second-order equations ( damped harmonic oscillators with additive noise), where the statistics are Gaussian and can be calculated exactly in the continuous-time and discrete-time cases. A matrix equation is given for the stationary variances and correlation for methods using one Gaussian random variable per timestep. The only Runge-Kutta method with a nonsingular tableau matrix that gives the exact steady state density for all values of damping is the implicit midpoint rule. Numerical experiments, comparing the implicit midpoint rule with Heun and leapfrog methods on nonlinear equations with additive or multiplicative noise, produce behavior similar to the linear case.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据