Univ.-Prof. Dr.in Evelyn Buckwar
Tel.: +43 732 2468 4161
Evelyn Buckwar und Thorsten Sickenberger
Software (MATLAB® files) für
Asymptotic mean-square stability analysis of linear systems of stochastic differential equations (SDEs) and their numerical approximation
basierend auf dem Manuskript:
E. Buckwar, T. Sickenberger:
A structural analysis of asymptotic mean-sqaure stability for multi-dimensional linear stochastic differential systems,
Heriot-Watt Mathematics Report HWM10-21, 2010.
Abstract: We are concerned with a linear mean-square stability analysis of numerical methods applied to systems of stochastic differential equations (SDEs) and, in particular, consider the Theta-Maruyama and the Theta-Milstein method in this context. We propose a technique, based on the vectorisation of matrices and the Kronecker product, that allows us to deal with the matrix expressions arising in this analysis and which provides the explicit structure of the stability matrices in the general case of linear systems of SDEs. For a set of simple test SDE systems, incorporating different noise structures but only a few parameters, we apply the general results and provide visual and numerical comparisons of the stability properties of the two methods.
Keywords: Asymptotic mean-square stability • Theta-Maruyama method • Theta-Milstein method • Systems of stochastic differential equations • Linear stability analysis
Software download: You can test the asymptotic mean-square stability of your own d-dimensional linear system of SDEs having m multiplicative noise sources. Download the following files and insert your drift and diffusion matrices (see the "else" part of the first loop in the main file). Alternatively you can run the test SDE systems of our paper.
• main_stability.m - main program file (MATLAB®),
defines test-equations and writes outputs
• stabanalysis.m - function file (MATLAB®),
analyses the system and computes conditions for the asymptotic mean-square stability of the zero solution
• stabcalculation.m - function file (MATLAB®),
computes conditions for the asymptotic mean-square stability of Theta-Maruyama and Theta-Milstein approximations
for different values of the step-size and the method parameter Theta.
• Evelyn Buckwar, e-mail: Evelyn.Buckwar[at]jku.at
• Thorsten Sickenberger [Homepage]
MATLAB® is a registered trademark of The MathWorks Inc. Click here to go to The Mathworks home page.