Tropical differential equations

Tropical differential equations are introduced and an algorithm is designed which tests solvability of a system of tropical linear differential equations within the complexity polynomial in the size of the system and in the absolute values of its coefficients. Moreover, we show that there exists a minimal solution, and the algorithm constructs it (in case of solvability). This extends a similar complexity bound established for tropical linear systems. In case of tropical linear differential systems in one variable a polynomial complexity algorithm for testing its solvability is designed.We prove also that the problem of solvability of a system of tropical non-linear differential equations in one variable is NP-hard, and this problem for arbitrary number of variables belongs to NP. Similar to tropical algebraic equations, a tropical differential equation expresses the (necessary) condition on the dominant term in the issue of solvability of a differential equation in power series.     

Galerkin approximations with embedded boundary conditions for retarded delay differential equations

Finite-dimensional approximations are developed for retarded delay differential equations (DDEs). The DDE system is equivalently posed as an initial-boundary value problem consisting of hyperbolic partial differential equations (PDEs). By exploiting the equivalence of partial derivatives in space and time, we develop a new PDE representation for the DDEs that is devoid of boundary conditions. The resulting boundary condition-free PDEs are discretized using the Galerkin method with Legendre polynomials as the basis functions, whereupon we obtain a system of ordinary differential equations (ODEs) that is a finite-dimensional approximation of the original DDE system. We present several numerical examples comparing the solution obtained using the approximate ODEs to the direct numerical simulation of the original non-linear DDEs. Stability charts developed using our method are compared to existing results for linear DDEs. The presented results clearly demonstrate that the equivalent boundary condition-free PDE formulation accurately captures the dynamic behaviour of the original DDE system and facilitates the application of control theory developed for systems governed by ODEs.     

A collocation approach for solving linear complex differential equations in rectangular domains

In this paper, a collocation method is presented to find the approximate solution of high‐order linear complex differential equations in rectangular domain. By using collocation points defined in a rectangular domain and the Bessel polynomials, this method transforms the linear complex differential equations into a matrix equation. The matrix equation corresponds to a system of linear equations with the unknown Bessel coefficients. The proposed method gives the analytic solution when the exact solutions are polynomials. Numerical examples are included to demonstrate the validity and applicability of the technique and the comparisons are made with existing results. The results show the efficiency and accuracy of the present work. All of the numerical computations have been performed on a computer using a program written in MATLAB v7.6.0 (R2008a). Copyright © 2012 John Wiley & Sons, Ltd.     

Isomonodromic differential equations and differential categories

We study isomonodromicity of systems of parameterized linear differential equations and related conjugacy properties of linear differential algebraic groups by means of differential categories. We prove that isomonodromicity is equivalent to isomonodromicity with respect to each parameter separately under a filtered-linearly closed assumption on the field of functions of parameters. Our result implies that one does not need to solve any non-linear differential equations to test isomonodromicity anymore. This result cannot be further strengthened by weakening the requirement on the parameters as we show by giving a counterexample. Also, we show that isomonodromicity is equivalent to conjugacy to constants of the associated parameterized differential Galois group, extending a result of P. Cassidy and M. Singer, which we also prove categorically. We illustrate our main results by a series of examples, using, in particular, a relation between the Gauss–Manin connection and parameterized differential Galois groups.     

Spectral methods for some singularly perturbed third order ordinary differential equations

Spectral methods with interface point are presented to deal with some singularly perturbed third order boundary value problems of reaction-diffusion and convection-diffusion types. First, linear equations are considered and then non-linear equations. To solve non-linear equations, Newton’s method of quasi-linearization is applied. The problem is reduced to two systems of ordinary differential equations. And, then, each system is solved using spectral collocation methods. Our numerical experiments show that the proposed methods are produce highly accurate solutions in little computer time when compared with the other methods available in the literature.      

Stein implicit Runge-Kutta methods with high stage order for large-scale ordinary differential equations

The Runge-Kutta method is one of the most popular implicit methods for the solution of stiff ordinary differential equations. For large problems, the main drawback of such methods is the cost required at each integration step for computing the solution of a nonlinear system of equations. In this paper, we propose to reduce the cost of the computation by transforming the linear systems arising in the application of Newton's method to Stein matrix equations. We propose an iterative projection method onto block Krylov subspaces for solving numerically such Stein matrix equations. Numerical examples are given to illustrate the performance of our proposed method.     

Solving systems of fractional differential equations using differential transform method

This paper presents approximate analytical solutions for systems of fractional differential equations using the differential transform method. The fractional derivatives are described in the Caputo sense. The application of differential transform method, developed for differential equations of integer order, is extended to derive approximate analytical solutions of systems of fractional differential equations. The solutions of our model equations are calculated in the form of convergent series with easily computable components. Some examples are solved as illustrations, using symbolic computation. The numerical results show that the approach is easy to implement and accurate when applied to systems of fractional differential equations. The method introduces a promising tool for solving many linear and nonlinear fractional differential equations.     

Fractional-order Legendre functions for solving fractional-order differential equations

In this article, a general formulation for the fractional-order Legendre functions (FLFs) is constructed to obtain the solution of the fractional-order differential equations. Fractional calculus has been used to model physical and engineering processes that are found to be best described by fractional differential equations. Therefore, an efficient and reliable technique for the solution of them is too important. For the concept of fractional derivative we will adopt Caputo’s definition by using Riemann–Liouville fractional integral operator. Our main aim is to generalize the new orthogonal functions based on Legendre polynomials to the fractional calculus. Also a general formulation for FLFs fractional derivatives and product operational matrices is driven. These matrices together with the Tau method are then utilized to reduce the solution of this problem to the solution of a system of algebraic equations. The method is applied to solve linear and nonlinear fractional differential equations. Illustrative examples are included to demonstrate the validity and applicability of the presented technique.     

Infinitely many solutions for perturbed impulsive fractional differential systems

In this paper, the existence of infinitely many solutions for perturbed systems of impulsive non-linear fractional differential equations including Lipschitz continuous non-linear terms is discussed. The approach is based on variational methods. In addition, examples are presented to illustrate the feasibility and effectiveness of the main results.     

Solving differential equations using modified Picard iteration

Many classes of differential equation are shown to be open to solution through a method involving a combination of a direct integration approach with suitably modified Picard iterative procedures. The classes of differential equations considered include typical initial value, boundary value and eigenvalue problems arising in physics and engineering and include non-linear as well as linear differential equations. Examples involving partial as well as ordinary differential equations are presented. The method is easy to implement on a computer and the solutions so obtained are essentially power series. With its conceptual clarity (differential equations are integrated directly), its uniform methodology (the overall approach is the same in all cases) and its straightforward computer implementation (the integration and iteration procedures require only standard commercial software), the modified Picard methods offer obvious benefits for the teaching of differential equations as well as presenting a basic but flexible tool-kit for the solution process itself.     

A class of linear differential dynamical systems with fuzzy matrices

This paper investigates the first order linear fuzzy differential dynamical systems with fuzzy matrices. We use a complex number representation of the α-level sets of the fuzzy system, and obtain the solution by employing such representation. It is applicable to practical computations and has also some implications for the theory of fuzzy differential equations. We then present some properties of the 2-dimensional dynamical systems and their phase portraits. Some examples are considered to show the richness of the theory and we can clearly see that new behaviors appear. We finally present some conclusions and new directions for further research in the area of fuzzy dynamical systems.     

Systems of matrix rational differential equations arising in connection with linear stochastic systems with Markovian jumping

In this paper, a class of systems of matrix nonlinear differential equations containing as particular cases the systems of coupled Riccati differential equations arising in connection with control of some linear stochastic systems is considered.The system of differential equations considered in this paper are converted in a suitable nonlinear differential equation on a finite-dimensional Hilbert space adequately choosen.This allows us to use the positivity properties of the linear evolution operator defined by the linear differential equations of Lyapunov type.Our aim is to investigate properties of stabilizing and bounded solutions of the considered differential equations and to obtain some conditions ensuring the existence of such solutions.Conditions providing the existence of a maximal solution (minimal solution respectively) with respect to some classes of global solutions are presented. It is shown that if the coefficients of the equations are periodic functions all these special solutions (stabilizing, maximal, minimal) are periodic functions, too.Whenever possible the probabilistic arguments were avoided and so the results proved in the paper appear as results in the field of differential equations with interest in themselves.     

A numerical scheme based on Bernoulli wavelets and collocation method for solving fractional partial differential equations with Dirichlet boundary conditions

In this paper, an efficient and accurate numerical method is presented for solving two types of fractional partial differential equations. The fractional derivative is described in the Caputo sense. Our approach is based on Bernoulli wavelets collocation techniques together with the fractional integral operator, described in the Riemann‐Liouville sense. The main characteristic behind this approach is to reduce such problems to those of solving systems of algebraic equations, which greatly simplifies the problem. By using Newton's iterative method, this system is solved and the solution of fractional partial differential equations is achieved. Some results concerning the error analysis are obtained. The validity and applicability of the method are demonstrated by solving four numerical examples. Numerical examples are presented in the form of tables and graphs to make comparisons with the results obtained by other methods and with the exact solutions much easier.     

Parameter estimation in stochastic differential equations with Markov chain Monte Carlo and non-linear Kalman filtering

This paper is concerned with parameter estimation in linear and non-linear Itô type stochastic differential equations using Markov chain Monte Carlo (MCMC) methods. The MCMC methods studied in this paper are the Metropolis–Hastings and Hamiltonian Monte Carlo (HMC) algorithms. In these kind of models, the computation of the energy function gradient needed by HMC and gradient based optimization methods is non-trivial, and here we show how the gradient can be computed with a linear or non-linear Kalman filter-like recursion. We shall also show how in the linear case the differential equations in the gradient recursion equations can be solved using the matrix fraction decomposition. Numerical results for simulated examples are presented and discussed in detail.     

常系数非齐次线性微分方程组的初等解法

利用初等变换将常系数非齐次线性微分方程组化为由若干个相互独立的高阶常系数非齐次线性微分方程组成的方程组,再利用高阶常系数齐次线性微分方程的特征根法和非齐次方程的待定系数法求该方程组的基本解组及特解,最后通过初等变换求原方程组的基本解组及特解,从而可求出其通解.     

A new Chebyshev polynomial approximation for solving delay differential equations

The purpose of this study is to give a Chebyshev polynomial approximation for the solution of mth-order linear delay differential equations with variable coefficients under the mixed conditions. For this purpose, a new Chebyshev collocation method is introduced. This method is based on taking the truncated Chebyshev expansion of the function in the delay differential equations. Hence, the resulting matrix equation can be solved, and the unknown Chebyshev coefficients can be found approximately. In addition, examples that illustrate the pertinent features of the method are presented, and the results of this investigation are discussed.     

An integral operational matrix based on Jacobi polynomials for solving fractional-order differential equations

This paper aims to construct a general formulation for the Jacobi operational matrix of fractional integral operator. Fractional calculus has been used to model physical and engineering processes that are found to be best described by fractional differential equations. Therefore, a reliable and efficient technique for the solution of them is too important. For the concept of fractional derivative we will adopt Caputo’s definition by using Riemann–Liouville fractional integral operator. Our main aim is to generalize the Jacobi integral operational matrix to the fractional calculus. These matrices together with the Tau method are then utilized to reduce the solution of this problem to the solution of a system of algebraic equations. The method is applied to solve linear and nonlinear fractional differential equations. Illustrative examples are included to demonstrate the validity and applicability of the presented technique.     

Comparison between variational iteration method and homotopy perturbation method for linear and nonlinear partial differential equations with the nonhomogeneous initial conditions

In this study, linear and nonlinear partial differential equations with the nonhomogeneous initial conditions are considered. We used Variational iteration method (VIM) and Homotopy perturbation method (HPM) for solving these equations. Both methods are used to obtain analytic solutions for different types of differential equations. Four examples are presented to show the application of the present techniques. In these schemes, the solution takes the form of a convergent series with easily computable components. The present methods perform extremely well in terms of efficiency and simplicity. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

An exponential matrix method for solving systems of linear differential equations

This paper presents an exponential matrix method for the solutions of systems of high‐order linear differential equations with variable coefficients. The problem is considered with the mixed conditions. On the basis of the method, the matrix forms of exponential functions and their derivatives are constructed, and then by substituting the collocation points into the matrix forms, the fundamental matrix equation is formed. This matrix equation corresponds to a system of linear algebraic equations. By solving this system, the unknown coefficients are determined and thus the approximate solutions are obtained. Also, an error estimation based on the residual functions is presented for the method. The approximate solutions are improved by using this error estimation. To demonstrate the efficiency of the method, some numerical examples are given and the comparisons are made with the results of other methods. Copyright © 2012 John Wiley & Sons, Ltd.     

Magnus integrators for solving linear-quadratic differential games

We consider Magnus integrators to solve linear-quadratic N$N$-player differential games. These problems require to solve, backward in time, non-autonomous matrix Riccati differential equations which are coupled with the linear differential equations for the dynamic state of the game, to be integrated forward in time. We analyze different Magnus integrators which can provide either analytical or numerical approximations to the equations. They can be considered as time-averaging methods and frequently are used as exponential integrators. We show that they preserve some of the most relevant qualitative properties of the solution for the matrix Riccati differential equations as well as for the remaining equations. The analytical approximations allow us to study the problem in terms of the parameters involved. Some numerical examples are also considered which show that exponential methods are, in general, superior to standard methods.