Analysis of estimation accuracy of the first moments of a Monte Carlo solution to an SDE with Wiener and Poisson components

In this paper, the estimation accuracy of the first moments of a numerical solution to an SDE with Wiener and Poisson components is investigated by a generalized explicit Euler method. Exact expressions for the mathematical expectation and variance of a test SDE solution are obtained. These expressions allow us to investigate the estimation accuracy obtained by a Monte Carlo method versus the SDE parameters, the integration step, and the size of the ensemble of simulated trajectories of the solution. The results of test numerical experiments are presented.     

Exact Solutions and Integrability of the Duffing – Van der Pol Equation

The force-free Duffing–Van der Pol oscillator is considered. The truncated expansions for finding the solutions are used to look for exact solutions of this nonlinear ordinary differential equation. Conditions on parameter values of the equation are found to have the linearization of the Duffing–Van der Pol equation. The Painlevé test for this equation is used to study the integrability of the model. Exact solutions of this differential equation are found. In the special case the approach is simplified to demonstrate that some well-known methods can be used for finding exact solutions of nonlinear differential equations. The first integral of the Duffing–Van der Pol equation is found and the general solution of the equation is given in the special case for parameters of the equation. We also demonstrate the efficiency of the method for finding the first integral and the general solution for one of nonlinear second-order ordinary differential equations.     

Numerical analysis of stochastic oscillators on supercomputers

The problem of numerical analysis of stochastic differential equations (SDEs) with oscillating solutions is investigated. The expectation and variance of SDE numerical solutions are shown as functions of the mesh size of integrating the generalized Euler method. Results of some numerical experiments on the simulation of linear and nonlinear stochastic oscillators on the supercomputer of the Siberian Supercomputer Center are presented.     

Haar wavelet solutions of nonlinear oscillator equations

In this paper, we present a numerical scheme using uniform Haar wavelet approximation and quasilinearization process for solving some nonlinear oscillator equations. In our proposed work, quasilinearization technique is first applied through Haar wavelets to convert a nonlinear differential equation into a set of linear algebraic equations. Finally, to demonstrate the validity of the proposed method, it has been applied on three type of nonlinear oscillators namely Duffing, Van der Pol, and Duffing–van der Pol. The obtained responses are presented graphically and compared with available numerical and analytical solutions found in the literature. The main advantage of uniform Haar wavelet series with quasilinearization process is that it captures the behavior of the nonlinear oscillators without any iteration. The numerical problems are considered with force and without force to check the efficiency and simple applicability of method on nonlinear oscillator problems.     

On the limit cycles,period-doubling,and quasi-periodic solutions of the forced Van der Pol-Duffing oscillator

In this paper, the limit cycles, period-doubling, and quasi-periodic solutions of the forced Van der Pol oscillator and the forced Van der Pol-Duffing oscillator are studied by combining the homotopy analysis method (HAM) with the multi-scale method analytically. Comparisons of the obtained solutions and numerical results show that this method is effective and convenient even when t is large enough, and the convergence of the approximate solutions is discussed by the so-called discrete square residual error. This method is a capable tool for solving this kind of nonlinear problems.     

Energy methods for stochastic differential equations

This article is devoted to the existence of strong solutions to stochastic differential equations (SDEs). Compared with Ito's theory, we relax the assumptions on the volatility term and replace the global Lipschitz continuity condition with a local Lipschitz continuity condition and a Hoelder continuity condition. In particular, our general SDE covers the Cox–Ingersoll–Ross SDE as a special case. We note that the general weak existence theory presumably extends to our general SDE (although the explicit time dependence of the drift term and the volatility term might require some extra considerations). However, avoiding weak existence theory we prove the existence of a strong solution directly using a priori estimates (the so-called energy estimates) derived from the SDE. The benefit of this approach is that the argument only requires some basic knowledge about stochastic and functional analysis. Moreover, the underlying principle has developed to become one of the cornerstones of the modern theory of partial differential equations (PDEs). In this sense, the general goal of this article is not just to establish the existence of a strong solution to the SDE under consideration but rather to introduce a new principle in the context of SDEs that has already proven to be successful in the context of PDEs.     

Weak antithetic MLMC estimation of SDEs with the Milstein scheme for low-dimensional Wiener processes

In this paper, we implement a weak Milstein Scheme to simulate low-dimensional stochastic differential equations (SDEs). We prove that combining the antithetic multilevel Monte-Carlo (MLMC) estimator introduced by Giles and Szpruch with the MLMC approach for weak SDE approximation methods by Belomestny and Nagapetyan, we can achieve a quadratic computational complexity in the inverse of the Root Mean Square Error (RMSE) when estimating expected values of smooth functionals of SDE solutions, without simulating Lévy areas and without requiring any strong convergence of the underlying SDE approximation method. By using appropriate discrete variables this approach allows us to calculate the expectation on the coarsest level of resolution by enumeration, which, for low-dimensional problems, results in a reduced computational effort compared to standard MLMC sampling. These theoretical results are also confirmed by a numerical experiment.     

A posteriori error analysis of round-off errors in the numerical solution of ordinary differential equations

We prove sharp, computable error estimates for the propagation of errors in the numerical solution of ordinary differential equations. The new estimates extend previous estimates of the influence of data errors and discretization errors with a new term accounting for the propagation of numerical round-off errors, showing that the accumulated round-off error is inversely proportional to the square root of the step size. As a consequence, the numeric precision eventually sets the limit for the pointwise computability of accurate solutions of any ODE. The theoretical results are supported by numerically computed solutions and error estimates for the Lorenz system and the van der Pol oscillator.     

Oscillatory region and asymptotic solution of fractional van der Pol oscillator via residue harmonic balance technique

In this paper, a powerfully analytical technique is proposed for predicting and generating the steady state solution of the fractional differential system based on the method of harmonic balance. The zeroth-order approximation using just one Fourier term is applied to predict the parametric function for the boundary between oscillatory and non-oscillatory regions of the fractional van der Pol oscillator. The unbalanced residues due to Fourier truncation are considered iteratively by solving linear algebraic equations to improve the accuracy of the solutions successively. The highly accurate solutions to the angular frequency and limit cycle of fractional van der Pol oscillator are obtained and compared. The results reveal that the technique described in this paper is very effective and simple for obtaining asymptotic solution of nonlinear system having fractional order derivative.     

Parameter estimation for a price series model by solving a linear SDE with a Poisson component

A model of a series of price increments with jumps is constructed based on a linear stochastic differential equation with a Poisson component. Some estimates of unknown parameters of the model and SDE are obtained by using the method of moments. A statistical simulation algorithm for solving an SDE with a Poisson component in general form is proposed. Results of numerical experiments are given.     

A non-standard-Euler–Maruyama scheme

We construct a non-standard finite difference numerical scheme to approximate stochastic differential equations (SDEs) using the idea of weighed step introduced by R.E. Mickens. We prove the strong convergence of our scheme under locally Lipschitz conditions of a SDE and linear growth condition. We prove the preservation of domain invariance by our scheme under a minimal condition depending on a discretization parameter and unconditionally for the expectation of the approximate solution. The results are illustrated through the geometric Brownian motion. The new scheme shows a greater behaviour compared with the Euler–Maruyama scheme and balanced implicit methods which are widely used in the literature and applications.     

A novel approach to construct numerical methods for stochastic differential equations

In this paper we propose a new numerical method for solving stochastic differential equations (SDEs). As an application of this method we propose an explicit numerical scheme for a super linear SDE for which the usual Euler scheme diverges.     

Numerical solution of stochastic differential equations on Supercomputers

This paper deals with some problems of accuracy of algorithms for the numerical solutions of stochastic differential equations (SDEs) versus the size of the ensemble of trajectories simulated. The problems of accuracy arise in estimating functionals of SDE solutions with increasing variance, highly asymmetric distributions, and an indefinite time of arrival of trajectories of solutions at the boundaries of given domains. Some ways of parallelization of statistical algorithms on a multiprocessor cluster are described. Results of numerical experiments performed on a supercomputer available at the Siberian Supercomputer Center are presented.     

Numerical study of the controlled Van der Pol oscillator in Chebyshev series

Summary The optimal control of the Van der Pol oscillator is investigated by a numerical method based on the expansion of the state function and the control strategy in Chebyshev series. The optimal control problem is reduced to a parameter optimization problem.

---

Resumé On étudie le contrôle optimal de l'oscillateur de Van der Pol à l'aide d'une méthode numérique qui utilise le développement des fonctions de position et de contrôle en série de Chebyshev. Le problème de contrôle optimal est réduit à un problème d'optimisation de paramètres.

     

Asymptotic expansion for forward-backward SDEs with jumps

This work provides a semi-analytic approximation method for decoupled forward-backward SDEs (FBSDEs) with jumps. In particular, we construct an asymptotic expansion method for FBSDEs driven by the random Poisson measures with σ-finite compensators as well as the standard Brownian motions around the small-variance limit of the forward SDE. We provide a semi-analytic solution technique as well as its error estimate for which we only need to solve essentially a system of linear ODEs. In the case of a finite jump measure with a bounded intensity, the method can also handle state-dependent and hence non-Poissonian jumps, which are quite relevant for many practical applications.     

Global error estimates for a uniparametric family of stiffly accurate Runge-Kutta collocation methods on singularly perturbed problems

An error analysis for a newly defined uniparametric family of stiffly accurate Runge-Kutta collocation methods when applied to initial value problems for singularly perturbed differential equations is carried out. The so-called SAFERK methods possess a first internal stage of explicit type and are based on collocation nodes. Sharp convergence results are obtained for these methods through the analysis of a sequence of higher index Differential Algebraic Equations. A numerical test with the Van der Pol oscillator reveals that the proposed error estimates are realistic whenever the stepsize h is large enough compared to the stiffness parameter ε.     

A Small Noise Asymptotic Expansion for Young SDE Driven by Fractional Brownian Motion: A Sharp Error Estimate With Malliavin Calculus

This article shows an analytically tractable small noise asymptotic expansion with a sharp error estimate for the expectation of the solution to Young’s pathwise stochastic differential equations (SDEs) driven by fractional Brownian motions with the Hurst index H > 1/2. In particular, our asymptotic expansion can be regarded as small noise and small time asymptotics by the error estimate with Malliavin culculus. As an application, we give an expansion formula in one-dimensional general Young SDE driven by fractional Brownian motion. We show the validity of the expansion through numerical experiments.     

Minimizing the variance of a mathematical expectation estimate for a diffusion process functional based on a parametric transformation of a parabolic boundary value problem

This paper deals with finding ways of reducing the variance of a mathematical expectation estimate for the functional of a diffusion process moving in a domain with an absorbing boundary. The estimate of mathematical expectation of the functional is obtained based on a numerical solution of stochastic differential equations (SDEs) by using the Euler method. A formula of the limiting variance is derived with decreasing integration step in the Euler method. A method of reducing the variance value of the estimate based on transformation of the parabolic boundary value problem corresponding to the diffusion process is proposed. Some numerical results are presented.     

A homotopy analysis method for limit cycle of the van der Pol oscillator with delayed amplitude limiting

In this work, a powerful analytical method, called Liao’s homotopy analysis method is used to study the limit cycle of a two-dimensional nonlinear dynamical system, namely the van der Pol oscillator with delayed amplitude limiting. It is shown that the solutions are valid for a wide range of variation of the system parameters. Comparison of the obtained solutions with those achieved by numerical solutions and by other perturbation techniques shows that the utilized method is effective and convenient to solve this type of problems with the desired order of approximation.