Liapunov equations for time-varying linear systems

Liapunov equations for time-varying linear systems in Hilbert space are given. This generalizes Datko's result for semigroups. Time-varying linear stochastic systems are also considered and analogous results are obtained.     

The variational iteration method for solving linear and nonlinear Volterra integral and integro-differential equations

The variational iteration method is used for solving the linear and nonlinear Volterra integral and integro-differential equations. The method is reliable in handling Volterra equations of the first kind and second kind in a direct manner without any need for restrictive assumptions. The method significantly reduces the size of calculations.     

On the numerical solution of nonlinear systems of Volterra integro-differential equations with delay arguments

Particular cases of nonlinear systems of delay Volterra integro-differential equations (denoted by DVIDEs) with constant delay τ > 0, arise in mathematical modelling of ‘predator–prey’ dynamics in Ecology. In this paper, we give an analysis of the global convergence and local superconvergence properties of piecewise polynomial collocation for systems of this type. Then, from the perspective of applied mathematics, we consider the Volterra’s integro-differential system of ‘predator–prey’ dynamics arising in Ecology. We analyze the numerical issues of the introduced collocation method applied to the ‘predator–prey’ system and confirm that we can achieve the expected theoretical orders of convergence.      

On the error analysis of the numerical solution of linear Volterra integral equations of the second kind

In the numerical solution of linear Volterra integral equations, two kinds of errors occur. If we use the collocation method, these errors are the collocation and numerical quadrature errors. Each error has its own effect in the accuracy of the obtained numerical solution. In this study we obtain an error bound that is sum of these two errors and using this error bound the relation between the smoothness of the kernel in the equation and also the length of the integration interval and each of these two errors are considered. Concluded results also are observed during the solution of some numerical examples.     

A new analytical method for solving systems of Volterra integral equations

In this paper, a new modified homotopy perturbation method (NHPM) is introduced for solving systems of Volterra integral equations of the second kind. Theorems of existence and uniqueness of the solutions to these equations are presented. Comparison of the results of applying the NHPM with those of the homotopy perturbation method and Adomian's decomposition method leads to significant consequences. Several examples, including the system of linear and nonlinear Volterra integral equations, are given to demonstrate the efficiency of the new method.     

Liapunov and lagrange stability: Inverse theorems for discontinuous systems

The main result of this paper is a converse Liapunov theorem which applies to systems of ordinary differential equations with a discontinuous righthand side. We treat both the problem of local stability of an equilibrium position and the problem of boundedness of solutions. In particular, we show that in order to achieve a necessary and sufficient condition in terms of continuous Liapunov functions, the classical definitions need to be strengthened in a convenient way. This work was motivated by the recently renewed interest in stabilization by discontinuous feedback and analysis of the state evolution with respect to bounded inputs. To achieve a more general treatment, the exposition is developed in the framework of differential inclusions theory.     

Bessel polynomial solutions of high-order linear Volterra integro-differential equations

In this study, a practical matrix method, which is based on collocation points, is presented to find approximate solutions of high-order linear Volterra integro-differential equations (VIDEs) under the mixed conditions in terms of Bessel polynomials. Numerical examples are included to demonstrate the validity and applicability of the technique and comparisons are made with the existing results. The results show the efficiency and accuracy of the present work. All of the numerical computations have been performed on the computer using a program written in MATLAB v7.6.0 (R2008a).     

Approximate solution of systems of Volterra integral equations with error analysis

This paper describes a procedure for solving the system of linear Volterra integral equations by means of the Sinc collocation method. A convergence and an error analysis are given; it is shown that the Sinc solution produces an error of order O(exp(?c N ^1/2)), where c>0 is a constant. This approximation reduces the system of integral equations to an explicit system of algebraic equations. The analytical results are illustrated with numerical examples that exhibit the exponential convergence rate.     

Stability of partial difference equations governing control gains in infinite-dimensional backstepping

We examine the stability properties of a class of LTV difference equations on an infinite-dimensional state space that arise in backstepping designs for parabolic PDEs. The nominal system matrix of the difference equation has a special structure: all of its powers have entries that are −1, 0, or 1, and all of the eigenvalues of the matrix are on the unit circle. The difference equation is driven by initial conditions, additive forcing, and a system matrix perturbation, all of which depend on problem data (for example, viscosity and reactivity in the case of a reaction–diffusion equation), and all of which go to zero as the discretization step in the backstepping design goes to zero. All of these observations, combined with the fact that the equation evolves only in a number of steps equal to the dimension of its state space, combined with the discrete Gronwall inequality, establish that the difference equation has bounded solutions. This, in turn, guarantees the existence of a state-feedback gain kernel in the backstepping control law. With this approach we greatly expand, relative to our previous results, the class of parabolic PDEs to which backstepping is applicable.     

A numerical method for solving systems of linear and nonlinear integral equations of the second kind by hat basis functions

In this paper, we use hat basis functions to solve the system of Fredholm integral equations (SFIEs) and the system of Volterra integral equations (SVIEs) of the second kind. This method converts the system of integral equations into a linear or nonlinear system of algebraic equations. Also, we consider the order of convergence of the method and show that it is O(h²). Application of the method on some examples show its accuracy and efficiency.     

A two step method for the numerical integration of stiff differential equations

We examine a single-step implicit-integration algorithm which is obtained by a modification of the well-known Simpson rule. The accuracy and stability properties of these methods are investigated. The obtained new method is a fourth-order numerical process and preserves the property of A-stability of the Simpson rule. Numerical results for the solution of certain stiff equations.     

Necessary first- and second-order optimality conditions in problems of control described by a system of Volterra difference equations

The terminal problem of optimal control described by a system of Volterra difference equations is considered. Constructively verifiable necessary first- and second-order optimality conditions are obtained under the assumption that the domain of control is open.     

A dissipative exponentially-fitted method for the numerical solution of the Schrödinger equation and related problems

In this paper a dissipative exponentially-fitted method for the numerical integration of the Schrödinger equation and related problems is developed. The method is called dissipative since is a nonsymmetric multistep method. An application to the the resonance problem of the radial Schrödinger equation and to other well known related problems indicates that the new method is more efficient than the corresponding classical dissipative method and other well known methods. Based on the new method and the method of Raptis and Cash a new variable-step method is obtained. The application of the new variable-step method to the coupled differential equations arising from the Schrödinger equation indicates the power of the new approach.     

Stability and control of nonlinear systems described by retarded functional equations: a review of recent results

This paper reports on recent results in a series of the work of the authors on the stability and nonlinear control for general dynamical systems described by retarded functional differential and difference equations. Both internal and external stability properties are studied. The corresponding Lyapunov and Razuminkhin characterizations for input-to-state and input-to-output stabilities are proposed. Necessary and sufficient Lyapunov-like conditions are derived for robust nonlinear stabilization. In particular, an explicit controller design procedure is developed for a new class of nonlinear time-delay systems. Lastly, sufficient assumptions, including a small-gain condition, are presented for guaranteeing the input-to-output stability of coupled systems comprised of retarded functional differential and difference equations.     

Approximate solution of a class of linear integro-differential equations by Taylor expansion method

In this paper, a simple and effective Taylor expansion method is presented for solving a class of linear integro-differential equations including those of Fredholm and of Volterra types. By means of the nth-order Taylor expansion of an unknown function at an arbitrary point, a linear integro-differential equation can be converted approximately to a system of linear equations for the unknown function itself and its first n derivatives under initial conditions. The nth-order approximate solution is exact for a polynomial of degree equal to or less than n. Some examples are given to illustrate the accuracy of this method.     

STWS approach for Hammerstein system of nonlinear Volterra integral equations of the second kind

In this paper, single-term Walsh series (STWS) method is applied to obtain the numerical solutions of Hammerstein systems of nonlinear Volterra integral equations of second kind (HSNVIES). Using the properties of the STWS method, HSNVIES can be easily converted into solvable recursive system of algebraic equations. Solutions obtained from the recursive system of algebraic equations are the solutions of the HSNVIES. Illustrative examples are provided with numerical solutions and the efficiency of this STWS method is also compared with the existing methods.     

A new approach for numerical solution of integro-differential equations via Haar wavelets

A new method is proposed for numerical solution of Fredholm and Volterra integro-differential equations of second kind. The proposed method is based on Haar wavelets approximation. Special characteristics of Haar wavelets approximation has been used in the derivation of this method. The new method is the extension of the recent work [Aziz and Siraj-ul-Islam, New algorithms for numerical solution of nonlinear Fredholm and Volterra integral equations using Haar wavelets, J. Comput. Appl. Math. 239 (2013), pp. 333–345] from integral equations to integro-differential equations. The method is specifically derived for nonlinear problems. Two new algorithms are also proposed based on this new method, one each for numerical solution of Fredholm and Volterra integro-differential equations. The proposed algorithms are generic and are applicable to all types of both nonlinear Fredholm and Volterra integro-differential equations of second kind. The cost of the new algorithms is considerably reduced by using the Broyden's method instead of Newton's method for solution of system of nonlinear equations. Most of the numerical methods designed for solution of integro-differential equations rely on some other technique for numerical integration. The advantage of our method is that it does not use numerical integration. The integrand is approximated using Haar wavelets approximation and then exact integration is performed. The method is tested on number of problems and numerical results are compared with existing methods in the literature. The numerical results indicate that accuracy of the obtained solutions is reasonably high even when the number of collocation points is small.