Adaptive stepsize based on control theory for stochastic differential equations

The numerical solution of stochastic differential equations (SDEs) has been focussed recently on the development of numerical methods with good stability and order properties. These numerical implementations have been made with fixed stepsize, but there are many situations when a fixed stepsize is not appropriate. In the numerical solution of ordinary differential equations, much work has been carried out on developing robust implementation techniques using variable stepsize. It has been necessary, in the deterministic case, to consider the best choice for an initial stepsize, as well as developing effective strategies for stepsize control-the same, of course, must be carried out in the stochastic case. In this paper, proportional integral (PI) control is applied to a variable stepsize implementation of an embedded pair of stochastic Runge-Kutta methods used to obtain numerical solutions of nonstiff SDEs. For stiff SDEs, the embedded pair of the balanced Milstein and balanced implicit method is implemented in variable stepsize mode using a predictive controller for the stepsize change. The extension of these stepsize controllers from a digital filter theory point of view via PI with derivative (PID) control will also be implemented. The implementations show the improvement in efficiency that can be attained when using these control theory approaches compared with the regular stepsize change strategy. (C) 2004 Elsevier B.V. All rights reserved.

Student Difficulties with units in differential equations in modelling contexts

First-year undergraduate engineering students' understanding of the units of factors and terms in first-order ordinary differential equations used in modelling contexts was investigated using diagnostic quiz questions. Few students appeared to realize that the units of each term in such equations must be the same, or if they did, nevertheless failed to apply that knowledge when needed. In addition, few students were able to determine the units of a proportionality factor in a simple equation. These results indicate that lecturers of modelling courses cannot take this foundational knowledge for granted and should explicitly include it in instruction.

Practical Stability and Vector-Lyapunov Functions for Impulsive Differential Equations with "Supremum"

Stability of nonlinear impulsive differential equations with "supremum" is studied. A special type of stability, combining two different measures and a dot product on a cone, is defined. Perturbing cone-valued piecewise continuous Lyapunov functions have been applied. Method of Razumikhin as well as comparison method for scalar impulsive ordinary differential equations have been employed.

Stability in Terms of Two Measures for Differential Equations with “Maxima”

Атанаска Георгиева, Стела Глухчева, Снежана Христова - Изследвана е устойчивостта на нелинейни диференциални уравнения с “максимуми” по отношение на две мерки. Приложени са две различни мерки за началните условия и за решението. Използван е методът на Разумихин, а също така и методът на сравнението на обикновени скаларни диференциални уравнения. Приложението на получените резултати и достатъчни условия за устойчивост е илюстрирано с пример.

Modeling some real phenomena by fractional differential equations

This paper deals with fractional differential equations, with dependence on a Caputo fractional derivative of real order. The goal is to show, based on concrete examples and experimental data from several experiments, that fractional differential equations may model more efficiently certain problems than ordinary differential equations. A numerical optimization approach based on least squares approximation is used to determine the order of the fractional operator that better describes real data, as well as other related parameters.

EXISTENCE RESULTS FOR NEUTRAL INTEGRO-DIFFERENTIAL EQUATIONS WITH UNBOUNDED DELAY

In this paper we discuss the existence of mild, strict and classical solutions for a class of abstract integro-differential equations in Banach spaces. Some applications to ordinary and partial integro-differential equations are considered.

Global solutions for abstract neutral differential equations

We study the existence of global solutions for a class of abstract neutral differential equation defined on the whole real axis. Some concrete applications related to ordinary and partial differential equations are considered. (C) 2009 Elsevier Ltd. All rights reserved.

Global solutions for abstract impulsive neutral differential equations

We study the existence of mild solutions for a class of impulsive neutral functional differential equation defined on the whole real axis. Some concrete applications to ordinary and partial neutral differential equations with impulses are considered. (C) 2010 Elsevier Ltd. All rights reserved.

Upwind solution of singular differential equations arising from steady channel flows

We study ordinary nonlinear singular differential equations which arise from steady conservation laws with source terms. An example of steady conservation laws which leads to those scalar equations is the Saint–Venant equations. The numerical solution of these scalar equations is sought by using the ideas of upwinding and discretisation of source terms. Both the Engquist–Osher scheme and the Roe scheme are used with different strategies for discretising the source terms.

Discontinuous local semiflows for Kurzweil equations leading to LaSalle`s invariance principle for differential systems with impulses at variable times

In this paper, we consider an initial value problem for a class of generalized ODEs, also known as Kurzweil equations, and we prove the existence of a local semidynamical system there. Under certain perturbation conditions, we also show that this class of generalized ODEs admits a discontinuous semiflow which we shall refer to as an impulsive semidynamical system. As a consequence, we obtain LaSalle`s invariance principle for such a class of generalized ODEs. Due to the importance of LaSalle`s invariance principle in studying stability of differential systems, we include an application to autonomous ordinary differential systems with impulse action at variable times. (C) 2011 Elsevier Inc. All rights reserved.

A multirate ordinary differential equation integrator /

Originally presented as the author's thesis, University of Illinois at Urbana-Champaign.

Theory of differential equations.

pt. I. (Vol. I) Exact equations and Pfaff's problem. 1890.--pt. II. (Vol. II-III) Ordinary equations, not linear. 1900.--pt. III. (Vol. IV) Ordinary equations. 1902.--pt. IV (vol. V-VI) Partial differential equations. 1906.

Variational mean-field algorithm for efficient inference in large systems of stochastic differential equations

This work introduces a Gaussian variational mean-field approximation for inference in dynamical systems which can be modeled by ordinary stochastic differential equations. This new approach allows one to express the variational free energy as a functional of the marginal moments of the approximating Gaussian process. A restriction of the moment equations to piecewise polynomial functions, over time, dramatically reduces the complexity of approximate inference for stochastic differential equation models and makes it comparable to that of discrete time hidden Markov models. The algorithm is demonstrated on state and parameter estimation for nonlinear problems with up to 1000 dimensional state vectors and compares the results empirically with various well-known inference methodologies.

Solving Differential Equations by Parallel Laplace Method with Assured Accuracy

The paper has been presented at the 12th International Conference on Applications of Computer Algebra, Varna, Bulgaria, June, 2006

Oscillation Theorems for Perturbed Second Order Nonlinear Differential Equations with Damping

Some oscillation criteria for solutions of a general perturbed second order ordinary differential equation with damping (r(t)x′ (t))′ + h(t)f (x)x′ (t) + ψ(t, x) = H(t, x(t), x′ (t)) with alternating coefficients are given. The results obtained improve and extend some existing results in the literature.