Abstract

We consider numerical schemes for computing the linear response of steady-state averages with respect to a perturbation of the drift part of the stochastic differential equation. The schemes are based on the Girsanov change-of-measure theory in order to reweight trajectories with factors derived from a linearization of the Girsanov weights. The resulting estimator is the product of a time average and a martingale correlated to this time average. We investigate both its discretization and finite-time approximation errors. The designed numerical schemes are shown to be of a bounded variance with respect to the integration time which is desirable feature for long time simulations. We also show how the discretization error can be improved to second-order accuracy in the time step by modifying the weight process in an appropriate way.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.