Numerical solution of Volterra integro-differential equations modelling thalamo-cortical systems

Our study concerns thalamo-cortical systems which are modelled by nonlinear systems of Volterra integro-differential equations of convolution type. The thalamo-cortical systems describe a new architecture for a neurocomputer. Such a computer employs principles of human brain. It consists of oscillators which have different frequencies and are weakly connected via a common medium forced by an external input. Since a neurocomputer consists of many interconnected oscillators (referred also as neurons), the thalamo-cortical systems include large numbers of Volterra integro-differential equations. Solving such systems numerically is expensive not only because of their large dimensions but also because of many kernel evaluations which are needed over the whole interval from the initial point, where the initial condition is imposed, up to the present point, where the computations are currently executed. Moreover, the whole computed history of the solution has to be stored in the memory of the computing machine. Therefore, robust and efficient numerical algorithms are needed for computer simulations for the solutions to the thalamocortical systems. In this paper, we illustrate an iteration technique to solve the thalamo-cortical systems. The proposed successive iterates are vector functions of time, which change the original problems into systems of easier and separated equations. Such separated equations can then be solved in parallel computing environments. Results of numerical experiments are presented for large numbers of oscillators. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the structure of additive systems of integers

A sum-and-distance system is a collection of finite sets of integers such that the sums and differences formed by taking one element from each set generate a prescribed arithmetic progression. Such systems, with two component sets, arise naturally in the study of matrices with symmetry properties and consecutive integer entries. Sum systems are an analogous concept where only sums of elements are considered. We establish a bijection between sum systems and sum-and-distance systems of corresponding size, and show that sum systems are equivalent to principal reversible cuboids, which are tensors with integer entries and a symmetry of ‘reversible square’ type. We prove a structure theorem for principal reversible cuboids, which gives rise to an explicit construction formula for all sum systems in terms of joint ordered factorisations of their component set cardinalities.

     

On the optimal control of stochastic systems with an exponential-of-integral performance index

In this paper a theory of optimal control is developed for stochastic systems whose performance is measured by the exponential of an integral form. Such a formulation of the cost function is shown to be not only general and useful but also analytically tractable. Starting with very general classes of stochastic systems, optimality conditions are obtained which exploit the multiplicative decomposability of the exponential-of-integral form. Specializing to partially observed systems of stochastic differential equations with Brownian Motion disturbances, optimality conditions are obtained which parallel those for systems with integral costs. Also treated are the special cases of linear systems with exponential of quadratic costs for which explicit optimal controls are obtainable. In addition, several general results of independent interest are obtained, which concern optimality of stochastic systems.     

Stability of Linear Parameter Varying and Linear Switching Systems

We consider stability of families of linear time‐varying systems, that are determined by a set of time‐varying parameters which adhere to certain rules. The conditions are general enough to encompass on the one hand stability questions for systems that are frequently called linear parameter varying systems in the literature and on the other hand also linear switching systems, in which parameter variations are allowed to have discontinuities. Combinations of these two sets of assumptions are also possible within the framework studied here. Under the assumption of irreducibility of the sets of system matrices, we show how to construct parameter dependent Lyapunov functions for the systems under consideration that exactly characterize the exponential growth rate. It is clear that such Lyapunov functions do not exist in general. But every system of our class can be reduced to a finite number of subsystems for which irreducibility holds.     

A theory for synchronization of dynamical systems

In this paper, a theory for synchronization of multiple dynamical systems under specific constraints is developed from a theory of discontinuous dynamical systems. The concepts on synchronization of two or more dynamical systems to specific constraints are presented. The synchronization, desynchronization and penetration of multiple dynamical systems to multiple specified constraints are discussed, and the necessary and sufficient conditions for such synchronicity are developed. The synchronicity of two dynamical systems to a single specific constraint and to multiple specific constraints is investigated. Finally, the synchronization and the corresponding complexity for multiple slave systems with multiple master systems are discussed briefly. The meaning of synchronization for dynamical systems with constraints is extended as a generalized, universal concept. The theory presented in this paper may be as a universal theory for dynamical systems. The paper provides a theoretic frame work in order to control slave systems which can be synchronized with master systems through specific constraints in a general sense.     

Approximate analytical solutions of Schnakenberg systems by homotopy analysis method

In this paper, the homotopy analysis method (HAM) has been employed to obtain analytical solution of a two reaction–diffusion systems of fractional order (fractional Schnakenberg systems) which has been modeling morphogen systems in developmental biology. Different from all other analytic methods, HAM provides us with a simple way to adjust and control the convergence region of solution series by choosing proper values for auxiliary parameter h. The fractional derivative is described in the Caputo sense. The reason of using fractional order differential equations (FOD) is that FOD are naturally related to systems with memory which exists in most biological systems. Also they are closely related to fractals which are abundant in biological systems. The results derived of the fractional system are of a more general nature. Respectively, solutions of FOD spread at a faster rate than the classical differential equations, and may exhibit asymmetry. However, the fundamental solutions of these equations still exhibit useful scaling properties that make them attractive for applications.     

On the basis property of a two-part trigonometric series

Conditions are presented under which two-part trigonometric systems arising in mixed type equations form a Riesz basis in the space of Lebesgue square integrable functions. For such systems, biorthogonal systems can be obtained in explicit form. As a result, an integral representation of the solution to the Frankl problem in a special domain can be found. The results are extended to two-part systems of broader functions.     

Output feedback stabilization for a class of nonlinear time-evolution systems

A variety of problems in nonlinear time-evolution systems such as communication networks, computer networks, manufacturing, traffic management, etc., can be modelled as min–max-plus systems in which operations of min, max and addition appear simultaneously. Systems with only maximum (or minimum) constraints can be modelled as max-plus system and handled by max-plus algebra which changes the original nonlinear system in the traditional sense into linear system in this framework. Min-max-plus systems are extensions of max-plus systems and nonlinear even in the max-plus algebra view. Output feedback stabilization for min–max-plus systems with min–max-plus inputs and max-plus outputs is considered in this paper. Max-plus projection representation for the closed-loop system with min–max-plus output feedback is introduced and the formula to calculate the cycle time is presented. Stabilization of reachable systems with at least one observable state and a further result for reachable and observable systems are worked out, during which max-plus output feedbacks are used to stabilize the systems. The method based on the max-plus algebra is constructive in nature.     

Galerkin's method and stability

Approximation in least squares by Galerkin's method leads to a consideration of strongly minimal systems. Theorems are derived which permit the recognition of systems which are not strongly minimal from the characteristics of the elements themselves. Normalised systems cannot be strongly minimal without their eigenvalues being bounded above. Speared systems, which have desirable properties, are introduced and their main features determined. Convergence earmarks and error bounds are exposed. A new definition of stability, which is self-checking in a computation, is suggested and its attributes delineated. The extension of the theory to equations involving positive-definite operators is mentioned.     

Linking systems of difference sets

A linking system of difference sets is a collection of mutually related group difference sets, whose advantageous properties have been used to extend classical constructions of systems of linked symmetric designs. The central problems are to determine which groups contain a linking system of difference sets, and how large such a system can be. All previous constructive results for linking systems of difference sets are restricted to 2‐groups. We use an elementary projection argument to show that neither the McFarland/Dillon nor the Spence construction of difference sets can give rise to a linking system of difference sets in non‐2‐groups. We make a connection to Kerdock and bent sets, which provides large linking systems of difference sets in elementary abelian 2‐groups. We give a new construction for linking systems of difference sets in 2‐groups, taking advantage of a previously unrecognized connection with group difference matrices. This construction simplifies and extends prior results, producing larger linking systems than before in certain 2‐groups, new linking systems in other 2‐groups for which no system was previously known, and the first known examples in nonabelian groups.     

Weierstrass integrability of differential equations

The integrability problem consists of finding the class of functions a first integral of a given planar polynomial differential system must belong to. We recall the characterization of systems which admit an elementary or Liouvillian first integral. We define Weierstrass integrability and we determine some Weierstrass integrable systems which are Liouvillian integrable. Inside this new class of integrable systems there are non-Liouvillian integrable systems.     

Asymptotic Behavior of Solutions of Dynamic Equations

We consider linear dynamic systems on time scales, which contain as special cases linear differential systems, difference systems, or other dynamic systems. We give an asymptotic representation for a fundamental solution matrix that reduces the study of systems in the sense of asymptotic behavior to the study of scalar dynamic equations. In order to understand the asymptotic behavior of solutions of scalar linear dynamic equations on time scales, we also investigate the behavior of solutions of the simplest types of such scalar equations, which are natural generalizations of the usual exponential function.     

The stability of a specific class of mechanical systems

A special class of mechanical systems is considered, the linearized equations of which either belong to the class of time-varying systems, reducible to stationary systems using constructive Lyapunov transformations or to systems close to these. A method of decomposing of the matrices of a system, which differs from the traditional method, is proposed for investigating of the stability of motion. It is shown that the conclusions concerning the stability are more complete in the case of this decomposition of the system matrix. A number of problems on the stability of motion of various mechanical systems is considered as examples.     

Irreducible Lie-Yamaguti algebras of generic type

Lie-Yamaguti algebras (or generalized Lie triple systems) are binary-ternary algebras intimately related to reductive homogeneous spaces. The Lie-Yamaguti algebras which are irreducible as modules over their inner derivation algebras are the algebraic counterparts of the isotropy irreducible homogeneous spaces.These systems splits into three disjoint types: adjoint type, non-simple type and generic type. The systems of the first two types were classified in a previous paper through a generalized Tits Construction of Lie algebras. In this paper, the Lie-Yamaguti algebras of generic type are classified by relating them to several other nonassociative algebraic systems: Lie and Jordan algebras and triple systems, Jordan pairs or Freudenthal triple systems.     

Controllability and stabilization of programmed motions of reversible mechanical and electromechanical systems

Problems of controllability and methods of stabilizing programmed motions of a large class of mechanical and electromechanical systems which are reversible with respect to the control are considered. Criteria of the controllability and stabilizability of reversible systems are obtained. Programmed motions and algorithms of programmed control are designed in analytical form and algorithms of programmed motions for non-linear reversible systems are synthesized.     

On the equivalence of differential systems

In this article. First, we construct some nonlinear differential systems which are equivalent to some known systems. Second, we discuss, in a different method, the equivalence between some linear differential systems. And then we apply the obtained results to the study of the qualitative properties of these systems simultaneously.     

On the modeling of nonlinear interactions in large complex systems

This work deals with the modeling of large systems of interacting entities in the framework of the mathematical kinetic theory for active particles. The contents are specifically focused on the modeling of nonlinear interactions which is one of the most important issues in the mathematical approach to modeling and simulating complex systems, and which includes a learning–hiding dynamics. Applications are focused on the modeling of complex biological systems and on immune competition.     

Solution spaces of homogeneous linear difference systems with coefficient matrices from commutative groups

We analyse the solution spaces of limit periodic homogeneous linear difference systems, where the coefficient matrices of the considered systems are taken from a commutative group which does not need to be bounded. In particular, we study such systems whose fundamental matrices are not asymptotically almost periodic or which have solutions vanishing at infinity. We identify a simple condition on the matrix group which guarantees that the studied systems form a dense subset in the space of all considered systems. The obtained results improve previously known theorems about non-almost periodic and non-asymptotically almost periodic solutions. Note that the elements of the coefficient matrices are taken from an infinite field with an absolute value and that the corresponding almost periodic case is treated as well.