A comparative study of some computer algebra packages which determine the lie point symmetries of differential equations

The study seeks to determine which of five computer algebra packages is best at finding the Lie point symmetries of systems of partial differential equations with minimal user intervention. The chosen packages are LIEPDE and DIMSYM for REDUCE, LIE and BIGLIE for MUMATH, DESOLV for MAPLE, and MATHLIE for MATHEMATICA. A series of systems of partial differential equations are used in the study. The paper concludes that while all of the computer packages are useful, DESOLV appears to be the most successful system at determining the complete set of Lie point symmetries of systems of partial differential equations.

Similarity solutions of partial differential equations using DESOLV

We present and describe new reduction routines included in DESOLV which, in many cases, may allow the complete automation of the determination of similarity solutions of partial differential equations.

Lie symmetries of partial differential equations using symbolic computing

This study presents a theoretical basis for and outlines the method of finding the Lie point symmetries of systems of partial differential equations. It seeks to determine which of five computer algebra packages is best at finding these symmetries. The chosen packages are LIEPDE and DIMSYM for REDUCE, LIE and BIGLIE for MUMATH, DESOLV for MAPLE, and MATHLIE for MATHEMATICA. This work concludes that while all of the computer packages are useful, DESOLV appears to be the most successful system at determining the complete set of Lie symmetries. Also, the study describes REDUCEVAR, a new package for MAPLE, that reduces the number of independent variables in systems of partial differential equations, using particular Lie point symmetries. It outlines the results of some testing carried out on this package. It concludes that REDUCEVAR is a very useful tool in performing the reduction of independent variables according to Lie's theory and is highly accurate in identifying cases where the symmetries are not suitable for finding S/G equations.

Finding higher symmetries of differential equations using the MAPLE package DESOLVII

We present and describe, with illustrative examples, the MAPLE computer algebra package DESOLVII, which is a major upgrade of DESOLV. DESOLVII now includes new routines allowing the determination of higher symmetries (contact and Lie-Backlund) for systems of both ordinary and partial differential equations.

Applications of symbolic computing for symmetry analysis of differential equations

 This thesis presents a number of applications of symbolic computing to the study of differential equations. In particular, three packages have been produced for the computer algebra system MAPLE and used to find a variety of symmetries (and corresponding invariant solutions) for a range of differential systems.

FracSym: automated symbolic computation of Lie symmetries of fractional differential equations

In this paper, we present an algorithm for the systematic calculation of Lie point symmetries for fractional order differential equations (FDEs) using the method as described by Buckwar & Luchko (1998) and Gazizov, Kasatkin & Lukashchuk (2007, 2009, 2011). The method has been generalised here to allow for the determination of symmetries for FDEs with n independent variables and for systems of partial FDEs. The algorithm has been implemented in the new MAPLE package FracSym (Jefferson and Carminati 2013) which uses routines from the MAPLE symmetry packages DESOLVII (Vu, Jefferson and Carminati, 2012) and ASP (Jefferson and Carminati, 2013). We introduce FracSym by investigating the symmetries of a number of FDEs; specific forms of any arbitrary functions, which may extend the symmetry algebras, are also determined. For each of the FDEs discussed, selected invariant solutions are then presented. © 2013 Elsevier B.V. All rights reserved.

Identification of nonlinear systems using hybrid functions

Most real systems have nonlinear behavior and thus model linearization may not produce an accurate representation of them. This paper presents a method based on hybrid functions to identify the parameters of nonlinear real systems. A hybrid function is a combination of two groups of orthogonal functions: piecewise orthogonal functions (e.g. Block-Pulse) and continuous orthogonal functions (e.g. Legendre polynomials). These functions are completed with an operational matrix of integration and a product matrix. Therefore, it is possible to convert nonlinear differential and integration equations into algebraic equations. After mathematical manipulation, the unknown linear and nonlinear parameters are identified. As an example, a mechanical system with single degree of freedom is simulated using the proposed method and the results are compared against those of an existing approach.

New generalized Halanay inequalities with applications to stability of nonlinear non-autonomous time-delay systems

In this paper, a general class of Halanay-type non-autonomous functional differential inequalities is considered. A new concept of stability, namely global generalized exponential stability, is proposed. We first prove some new generalizations of the Halanay inequality. We then derive explicit criteria for global generalized exponential stability of nonlinear non-autonomous time-delay systems based on our new generalized Halanay inequalities. Numerical examples and simulations are provided to illustrate the effectiveness of the obtained results.

An observer design procedure for a class of nonlinear time-delay systems

This paper considers a class of uncertain, nonlinear differential state delayed control systems and presents a reduced-order observer design procedure to asymptotically estimate any vector state functionals. The method proposed involves decomposition of the delayed portion of the system into two parts: a matched and mismatched part. Provided that the rank of the mismatched part is less than the number of the outputs, a reduced-order linear functional observer, with any prescribed stability margin, can be constructed by using a simple procedure. A numerical example is given to illustrate the new design procedure and its features.

Semi-hyperbolic mappings in Banach spaces

The definition of semi-hyperbolic dynamical systems generated by Lipschitz continuous and not necessarily invertible mappings in Banach spaces is presented in this thesis. Like hyperbolic mappings, they involve a splitting into stable and unstable spaces, but a slight leakage from the strict invariance of the spaces is possible and the unstable subspaces are assumed to be finite dimensional. Bi-shadowing is a combination of the concepts of shadowing and inverse shadowing and is usually used to compare pseudo-trajectories calculated by a computer with the true trajectories. In this thesis, the concept of bi-shadowing in a Banach space is defined and proved for semi-hyperbolic dynamical systems generated by Lipschitz mappings. As an application to the concept of bishadowing, linear delay differential equations are shown to be bi-shadowing with respect to pseudo-trajectories generated by nonlinear small perturbations of the linear delay equation. This shows robustness of solutions of the linear delay equation with respect to small nonlinear perturbations. Complicated dynamical behaviour is often a consequence of the expansivity of a dynamical system. Semi-hyperbolic dynamical systems generated by Lipschitz mappings on a Banach space are shown to be exponentially expansive, and explicit rates of expansion are determined. The result is applied to a nonsmooth noninvertible system generated by delay differential equation. It is shown that semi-hyperbolic mappings are locally φ-contracting, where -0 is the Hausdorff measure of noncompactness, and that a linear operator is semi-hyperbolic if and only if it is φ-contracting and has no spectral values on the unit circle. The definition of φ-bi-shadowing is given and it is shown that semi-hyperbolic mappings in Banach spaces are φ-bi-shadowing with respect to locally condensing continuous comparison mappings. The result is applied to linear delay differential equations of neutral type with nonsmooth perturbations. Finally, it is shown that a small delay perturbation of an ordinary differential equation with a homoclinic trajectory is ‘chaotic’.

Application of real and functional analysis to solve boundary value problems

This thesis is about using appropriate tools in functional analysis arid classical analysis to tackle the problem of existence and uniqueness of nonlinear partial differential equations. There being no unified strategy to deal with these equations, one approaches each equation with an appropriate method, depending on the characteristics of the equation. The correct setting of the problem in appropriate function spaces is the first important part on the road to the solution. Here, we choose the setting of Sobolev spaces. The second essential part is to choose the correct tool for each equation. In the first part of this thesis (Chapters 3 and 4) we consider a variety of nonlinear hyperbolic partial differential equations with mixed boundary and initial conditions. The methods of compactness and monotonicity are used to prove existence and uniqueness of the solution (Chapter 3). Finding a priori estimates is the main task in this analysis. For some types of nonlinearity, these estimates cannot be easily obtained, arid so these two methods cannot be applied directly. In this case, we first linearise the equation, using linear recurrence (Chapter 4). In the second part of the thesis (Chapter 5), by using an appropriate tool in functional analysis (the Sobolev Imbedding Theorem), we are able to improve previous results on a posteriori error estimates for the finite element method of lines applied to nonlinear parabolic equations. These estimates are crucial in the design of adaptive algorithms for the method, and previous analysis relies on, what we show to be, unnecessary assumptions which limit the application of the algorithms. Our analysis does not require these assumptions. In the last part of the thesis (Chapter 6), staying with the theme of choosing the most suitable tools, we show that using classical analysis in a proper way is in some cases sufficient to obtain considerable results. We study in this chapter nonexistence of positive solutions to Laplace's equation with nonlinear Neumann boundary condition. This problem arises when one wants to study the blow-up at finite time of the solution of the corresponding parabolic problem, which models the heating of a substance by radiation. We generalise known results which were obtained by using more abstract methods.

Exponential stabilization of non-autonomous delayed neural networks via Riccati equations

This paper concerns with the problem of exponential stabilization for a class of non-autonomous neural networks with mixed discrete and distributed time-varying delays. Two cases of discrete time-varying delay, namely (i) slowly time-varying; and (ii) fast time-varying, are considered. By constructing an appropriate Lyapunov-Krasovskii functional in case (i) and utilizing the Razumikhin technique in case (ii), we establish some new delay-dependent conditions for designing a memoryless state feedback controller which stabilizes the system with an exponential convergence of the resulting closed-loop system. The proposed conditions are derived through solutions of some types of Riccati differential equations. Applications to control a class of autonomous neural networks with mixed time-varying delays are also discussed in this paper. Some numerical examples are provided to illustrate the effectiveness of the obtained results.