Invasive weed optimization for model order reduction of linear MIMO systems

In this work, a model order reduction (MOR) technique for a linear multivariable system is proposed using invasive weed optimization (IWO). This technique is applied with the combined advantages of retaining the dominant poles and the error minimization. The state space matrices of the reduced order system are chosen such that the dominant eigenvalues of the full order system are unchanged. The other system parameters are chosen using the invasive weed optimization with objective function to minimize the mean squared errors between the outputs of the full order system and the outputs of the reduced order model when the inputs are unit step. The proposed algorithm has been applied successfully, a 10th order Multiple-Input–Multiple-Output (MIMO) linear model for a practical power system was reduced to a 3rd order and compared with recently published work.     

On Padé-type model order reduction of J-Hermitian linear dynamical systems

A simple, yet powerful approach to model order reduction of large-scale linear dynamical systems is to employ projection onto block Krylov subspaces. The transfer functions of the resulting reduced-order models of such projection methods can be characterized as Padé-type approximants of the transfer function of the original large-scale system. If the original system exhibits certain symmetries, then the reduced-order models are considerably more accurate than the theory for general systems predicts. In this paper, the framework of J-Hermitian linear dynamical systems is used to establish a general result about this higher accuracy. In particular, it is shown that in the case of J-Hermitian linear dynamical systems, the reduced-order transfer functions match twice as many Taylor coefficients of the original transfer function as in the general case. An application to the SPRIM algorithm for order reduction of general RCL electrical networks is discussed.     

Reduced order modeling of nonlinear time periodic systems subjected to external periodic excitations

In this work, new methodologies for order reduction of nonlinear systems with periodic coefficients subjected to external periodic excitations are presented. The periodicity of the linear terms is assumed to be non-commensurate with the periodicity of forcing vector. The dynamical equations of motion are transformed using the Lyapunov–Floquet (L–F) transformation such that the linear parts of the resulting equations become time-invariant while the forcing and nonlinearity takes the form of quasiperiodic functions. The techniques proposed here construct a reduced order equivalent system by expressing the non-dominant states as time-varying functions of the dominant (master) states. This reduced order model preserves stability properties and is easier to analyze, simulate and control since it consists of relatively small number of states in comparison with the large scale system.Specifically, two methods are discussed to obtain the reduced order model. First approach is a straightforward application of linear method similar to the ‘Guyan reduction’. The second novel technique proposed here extends the concept of ‘invariant manifolds’ for the forced problem to construct the fundamental solution. Order reduction approach based on this extended invariant manifold technique yields unique ‘reducibility conditions’. If these ‘reducibility conditions’ are satisfied only then an accurate order reduction via extended invariant manifold approach is possible. This approach not only yields accurate reduced order models using the fundamental solution but also explains the consequences of various ‘primary’ and ‘secondary resonances’ present in the system. One can also recover ‘resonance conditions’ associated with the fundamental solution which could be obtained via perturbation techniques by assuming weak parametric excitation. This technique is capable of handling systems with strong parametric excitations subjected to periodic and quasi-periodic forcing. It is anticipated that these order reduction techniques will provide a useful tool in the analysis and control system design of large-scale parametrically excited nonlinear systems subjected to external periodic excitations.     

Low‐gain adaptive stabilization of semilinear second‐order hyperbolic systems

In this paper low‐gain adaptive stabilization of undamped semilinear second‐order hyperbolic systems is considered in the case where the input and output operators are collocated. The linearized systems have an infinite number of poles and zeros on the imaginary axis. The adaptive stabilizer is constructed by a low‐gain adaptive velocity feedback. The closed‐loop system is governed by a non‐linear evolution equation. First, the well‐posedness of the closed‐loop system is shown. Next, an energy‐like function and a multiplier function are introduced and the exponential stability of the closed‐loop system is analysed. Some examples are given to illustrate the theory. Copyright © 2004 John Wiley & Sons, Ltd.     

Stabilization of infinite‐dimensional second‐order system by adaptive PI‐controllers

In this paper adaptive stabilization of infinite‐dimensional undamped second‐order systems is considered in the case where the input and output operators are collocated. The systems have an infinite number of poles and zeros on the imaginary axis. The adaptive stabilizer is constructed by an adaptive PI‐controller (proportional plus integral controller). An energy‐like function and a multiplier function are introduced and adaptive stabilization of the linear second‐order systems is analyzed. Copyright © 2001 John Wiley & Sons, Ltd.     

Uniform exponential attractors for first order non-autonomous lattice dynamical systems

In Abdallah (2008, 2009) [2] and [3], we have investigated the existence of exponential attractors for first and second order autonomous lattice dynamical systems. Within this work, in l², we carefully study the existence of a uniform exponential attractor for the family of processes associated with an abstract family of first order non-autonomous lattice dynamical systems with quasiperiodic symbols acting on a closed bounded set.     

Using the phase plane and the computer for the study of systems dynamics

Because of their graphics capabilities, it is possible to use personal computers to analyse the behaviour of dynamical systems. We therefore require a way to introduce pedagogically both this technique and the important aspects of the study of non‐linear systems. This paper considers first the formal points of the phase plane method, and then presents three relevant examples: (a) second order linear systems; (b) the Van der Pol oscillator; and (c) the forced pendulum. In each case programs in Turbo‐Pascal are given for plotting the trajectories.     

On macromodeling of nonlinear systems with time periodic coefficients

Some new techniques for reduced order (macro) modeling of nonlinear systems with time periodic coefficients are discussed in this paper. The dynamical evolution equations are transformed using the Lyapunov–Floquet (L–F) transformation such that the linear parts of the new set of equations become time-invariant. The techniques presented here reduce the order of this transformed system and all original states are obtained via the appropriate transformations. This macromodel preserves the desired stability and bifurcation characteristics of the original large-scale system and due to relatively few states; it is suitable for simulation and controller design.In this work, methodologies based on linear and nonlinear projections as well as ‘time periodic invariant manifold’ idea are presented. The invariant manifold technique yields a ‘reducibility condition’ that determines when an accurate nonlinear order reduction is possible. A comparative study of these order reduction methods is also included. These techniques are compared by means of time traces and Poincaré maps. A numerical error analysis is also included and advantages and limitations are discussed by means of a practical example.     

A theoretical analysis of a new second order finite volume approximation based on a low-order scheme using general admissible spatial meshes for the one dimensional wave equation

This paper details our note [6] and it is an extension of our previous works  and  which dealt with first order (both in time and space) and second order time accurate (second order in time and first order in space) implicit finite volume schemes for second order hyperbolic equations with Dirichlet boundary conditions on general nonconforming multidimensional spatial meshes introduced recently in [14]. We aim in this work (and some forthcoming studies) to get higher order (both in time and space) finite volume approximations for the exact solution of hyperbolic equations using the class of spatial generic meshes introduced recently in [14] on low order schemes from which the matrices used to compute the discrete solutions are sparse. We focus in the present contribution on the one dimensional wave equation and on one of its implicit finite volume schemes described in [4]. The implicit finite volume scheme approximating the one dimensional wave equation we consider (hereafter referred to as the basic finite volume scheme) yields linear systems to be solved successively. The matrices involved in these linear systems are tridiagonal, symmetric and definite positive. The finite volume approximate solution of the basic finite volume scheme is of order h+k$h+k$, where h (resp. k  ) is the mesh size of the spatial (resp. time) discretization. We construct a new finite volume approximation of order (h+k)²${(h+k)}^2$ in several discrete norms which allow us to get approximations of order two for the exact solution and its first derivatives. This new high-order approximation can be computed using linear systems whose matrices are the same ones used to compute the discrete solution of the basic finite volume scheme while the right hand sides are corrected. The construction of these right hand sides includes the approximation of some high order spatial derivatives of the exact solution. The computation of the approximation of these high order spatial derivatives can be performed using the same matrices stated above with another two tridiagonal matrices. The manner by which this new high-order approximation is constructed can be repeated to compute successively finite volume approximations of arbitrary order using the same matrices stated above. These high-order approximations can be obtained on any one dimensional admissible finite volume mesh in the sense of [13] without any condition. To reach the above results, a theoretical framework is developed and some numerical examples supporting the theory are presented. Some of the tools of this framework are new and interesting and they are stated in the one space dimension but they can be extended to several space dimensions. In particular a new and useful a prior estimate for a suitable discrete problem is developed and proved. The proof of this a prior estimate result is based essentially on the decomposition of the solution of the discrete problem into the solutions of two suitable discrete problems. A new technique is used in order to get a convenient finite volume approximation whose discrete time derivatives of order up to order two are also converging towards the solution of the wave equation and their corresponding time derivatives.     

Numerical controllability of fractional dynamical systems

In this paper, we provide a computational procedure for controlled state and steering control for linear and nonlinear fractional dynamical systems of order 1 < α ≤ 2 in finite-dimensional spaces using the Mittag-Leffler matrix function and the iterative technique. Some numerical examples are provided to illustrate the results.     

A kind of nonnegative matrices and its application on the stability of discrete dynamical systems

In this paper we introduce a new kind of nonnegative matrices which is called (sp) matrices. We show that the zero solutions of a class of linear discrete dynamical systems are asymptotically stable if and only if the coefficient matrices are (sp) matrices. To determine that a matrix is (sp) matrix or not is very simple, we need only to verify that some elements of the coefficient matrices are zero or not. According to the result above, we obtain the conditions for the stability of several classes of discrete dynamical systems.     

Second order parabolic equations in Banach spaces with dynamic boundary conditions

In this paper, we exhibit a unified treatment of the mixed initial boundary value problem for second order (in time) parabolic linear differential equations in Banach spaces, whose boundary conditions are of a dynamical nature. Results regarding existence, uniqueness, continuous dependence (on initial data) and regularity of classical and strict solutions are established. Moreover, several examples are given as samples for possible applications.

     

Mesh‐independent convergence of the modified inexact Newton method for a second order non‐linear problem

In this paper, we consider an inexact Newton method applied to a second order non‐linear problem with higher order non‐linearities. We provide conditions under which the method has a mesh‐independent rate of convergence. To do this, we are required, first, to set up the problem on a scale of Hilbert spaces and second, to devise a special iterative technique which converges in a higher than first order Sobolev norm. We show that the linear (Jacobian) system solved in Newton's method can be replaced with one iterative step provided that the initial non‐linear iterate is accurate enough. The closeness criteria can be taken independent of the mesh size. Finally, the results of numerical experiments are given to support the theory. Published in 2005 by John Wiley & Sons, Ltd.     

Convergence and stability of extended block boundary value methods for Volterra delay integro-differential equations

In this paper, we construct a class of extended block boundary value methods (B₂VMs) for Volterra delay integro-differential equations and analyze the convergence and stability of the methods. It is proven under the classical Lipschitz condition that an extended B₂VM is convergent of order p if the underlying boundary value methods (BVM) has consistent order p. The analysis shows that a B₂VM extended by an A-stable BVM can preserve the delay-independent stability of the underlying linear systems. Moreover, under some suitable conditions, the extended B₂VMs can also keep the delay-dependent stability of the underlying linear systems. In the end, we test the computational effectiveness by applying the introduced methods to the Volterra delay dynamical model of two interacting species, where the theoretical precision of the methods is further verified.     

Monotonicity and positive invariance of linear systems over dioids

In this paper, we make connections between two apparently different concepts. The first concept is the (linear) monotonicity of a given matrix which is usually used in order to compare Markov chains. This concept is involved in the simplification of complex stochastic systems in order to control the approximation error made. The second concept is the positive invariance of sets by a (linear) map. The properties of positively invariant sets are involved in many different problems in classical control theory, such as constrained control, robustness analysis, optimisation, and also in aggregation of Markov chains (namely strong lumpability and coherency).

In the context of linear dynamical systems over semirings which play an important role in the study of discrete event systems, we establish links between monotone (or isotone) linear maps and linear maps which admit some special families of positively invariant sets.     

On Fourier series for higher order (partial) derivatives of functions

This paper is focused on higher order differentiation of Fourier series of functions. By means of Stokes's transformation, the recursion relations between the Fourier coefficients in Fourier series of different order (partial) derivatives of the functions as well as the general formulas for Fourier series of higher order (partial) derivatives of the functions are acquired. And then, the sufficient conditions for term‐by‐term differentiation of Fourier series of the functions are presented. These findings are subsequently used to reinvestigate the Fourier series methods for linear elasto‐dynamical systems. The results given in this paper on the constituent elements, together with their combinatorial modes and numbering, of the sets of coefficients concerning 2rth order linear differential equation with constant coefficients are found to be different from the results deduced by Chaudhuri back in 2002. And it is also shown that the displacement solution proposed by Li in 2009 is valid only when the second order mixed partial derivative of the displacement vanishes at all of the four corners of the rectangular plate. Copyright © 2016 John Wiley & Sons, Ltd.     

Robust tube-based MPC of constrained piecewise affine systems with bounded additive disturbances

In this article, we propose a robust tube-based MPC formulation for a class of hybrid systems, namely autonomously switched PWA systems, with bounded additive disturbances. The term tube-based refers to those control techniques whose objective is to maintain all possible trajectories of the uncertain system inside a tube which is a set around the nominal (or reference) system trajectory, that is free from disturbances. Common methods in tube-based control systems consider an error dynamical system as the difference between the state of the nominal system and the state of the perturbed system. However, this definition of the error dynamical system leads to a complicated switched affine system for PWA systems. Therefore, we use a new notion of the reference system similar to the nominal system except that the switching between the various modes of the PWA system is driven by the state of the real system. Using this reference system instead of the nominal system leads us to an error dynamical system that can be modeled as a switched linear system. We employ a switched linear controller to stabilize this error system under arbitrary switching. This auxiliary controller forces the states of the uncertain system to remain in a tube confined to the invariant set around the state of the reference system. We add new constraints and tighten some other constraints of the nominal hybrid MPC for the reference system, in order to ensure convergence of the uncertain system and to guarantee robust exponential stability of the closed-loop system.     

Some Local Properties of Spectrum of Linear Dynamical Systems in Hilbert Space

We prove some local properties of the spectrum of a linear dynamical system in Hilbert space. The semigroup generator, the control operator and the observation operator may be unbounded. We consider (i) the PBH test, (ii) the correspondence between the poles of the resolvent of the semigroup generator and the poles of the transfer function, and (iii) pole-zero cancellation between two transfer functions of the cascade connection of two dynamical systems. For our investigation we take well-posed linear systems and a subclass of them called weakly regular systems as the most general setting.     

Generalized model matching and (F,G)-invariant submodules for linear systems over rings

A generalized version of the exact model matching problem (GEMMP) is considered for linear multivariable systems over an arbitrary commutative ring K with identity. Reduced forms of this problem are introduced, and a characterization of all solutions and minimal order solutions is given, both with and without the properness constraint on the solutions, in terms of linear equations over K and K-modules. An approach to the characterization of all stable solutions is presented which, under a certain Bezout condition and a freeness condition, provides a parametrization of all stable solutions. The results provide an explicit parametrization of all solutions and all stable solutions in case K is a field, without the Bezout condition. This is achieved through a very simple characterization and a generalization to an arbitrary field K of the “fixed poles” of the model matching problem in terms of invariant factors of a certain polynomial matrix. The results also show that whenever the GEMMP has a solution, there exist solutions whose poles can be chosen arbitrarily as far as they contain the “fixed poles” with the right multiplicities (in the algebraic closure of K). Implications of these results in regard to inverse systems are shown. Equivalent simpler forms (in state space form) of the problem are shown to be obtainable. A theory of finitely generated (F,G)-invariant submodules for linear systems over rings is developed, and the geometric equivalent of the model matching problem—the dynamic cover problem—is formulated, to which the results of the previous sections provide a solution in the reduced case.     

Synthesis of minimum-order robust inverters

In the present paper, we consider the inversion problem for dynamical systems, that is, the problem of reconstruction of the unknown input signal ξ(t) of a given system on the basis of known information (about either the complete phase vector or a measurable output of the system). An auxiliary dynamical system forming the desired estimate of the signal ξ(t) is called an inverter.In earlier papers of the authors, attention was mainly paid to the possibility of inversion of a dynamical system in different cases in principle. In this relation, a model of dynamical systems with some stabilizing control was used as an inverter for the solution of the problem; moreover, this control was often designed with the use of an additional dynamical system, an observer of the phase vector of the original system or the system in deviations. Thus, a dynamical system whose dimension either coincides with the dimension of the original system or exceeds it was considered as an inverter.In the solution of practical problems, it is often required to synthesize inverters of minimal order. (This requirement is related to constraints on the complexity, cost, and operation speed of automated control systems.) In the present paper, we consider the problem on the possible reduction of the order of the inverter in various cases and the problem on the construction of inverters of minimal order.