Search for: All records

This article investigates the ‐optimal estimation problem of a class of linear system with delays in states, disturbance input, and outputs. The estimator uses an extended Luenberger estimator format which estimates both the present and history states. The estimator is designed using an equivalent Partial Integral Equation (PIE) representation of the coupled nominal system. The advantage of the resulting PIE representation is compact and delay free—obviating the need for commonly used bounding technique such as integral inequalities which typically introduces conservatism into the resulting optimization problem. The ‐optimal estimator synthesis problem is then reformulated as a Linear Partial Inequality (LPI)—a form of convex optimization using operator variables and inequlities. Such LPI‐based optimization problems can be solved using semidefinite programming via the PIETOOLS toolbox in Matlab. Compared with previous work, the proposed method simplifies the analysis and computation process and resulting in observers which are non‐conservtism to 4 decimal places when compared with Pad‐based ODE observer design methodologies. Numerical examples and simulation results are given to illustrate the effectiveness and scalability of the proposed approach.

 

It has recently been shown that the evolution of a linear Partial Differential Equation (PDE) can be more conveniently represented in terms of the evolution of a higher spatial derivative of the state. This higher spatial derivative (termed the `fundamental state') lies in $L_2$ - requiring no auxiliary boundary conditions or continuity constraints. Such a representation (termed a Partial Integral Equation or PIE) is then defined in terms of an algebra of bounded integral operators (termed Partial Integral (PI) operators) and is constructed by identifying a unitary map from the fundamental state to the state of the original PDE. Unfortunately, when the PDE is nonlinear, the dynamics of the associated fundamental state are no longer parameterized in terms of PI operators. However, in this paper we show that such dynamics can be compactly represented using a new tensor algebra of partial integral operators acting on the tensor product of the fundamental state. We further show that this tensor product of the fundamental state forms a natural distributed equivalent of the monomial basis used in representation of polynomials on a finite-dimensional space. This new representation is then used to provide a simple SDP-based Lyapunov test of stability of quadratic PDEs. The test is applied to three illustrative examples of quadratic PDEs. 

Modern control theory provides us with a spectrum of methods for studying the interconnection of dynamic systems using input-output properties of the interconnected subsystems. Perhaps the most advanced framework for such inputoutput analysis is the use of Integral Quadratic Constraints (IQCs), which considers the interconnection of a nominal linear system with an unmodelled nonlinear or uncertain subsystem with known input-output properties. Although these methods are widely used for Ordinary Differential Equations (ODEs), there have been fewer attempts to extend IQCs to infinitedimensional systems. In this paper, we present an IQC-based framework for Partial Differential Equations (PDEs) and Delay Differential Equations (DDEs). First, we introduce infinitedimensional signal spaces, operators, and feedback interconnections. Next, in the main result, we propose a formulation of hard IQC-based input-output stability conditions, allowing for infinite-dimensional multipliers. We then show how to test hard IQC conditions with infinite-dimensional multipliers on a nominal linear PDE or DDE system via the Partial Integral Equation (PIE) state-space representation using a sufficient version of the Kalman-Yakubovich-Popov lemma (KYP). The results are then illustrated using four example problems with uncertainty and nonlinearity. 

In this paper, we present the Partial Integral Equation (PIE) representation of linear Partial Differential Equations (PDEs) in one spatial dimension, where the PDE has spatial integral terms appearing in the dynamics and the boundary conditions. The PIE representation is obtained by performing a change of variable where every PDE state is replaced by its highest, well-defined derivative using the Fundamental Theorem of Calculus to obtain a new equation (a PIE). We show that this conversion from PDE representation to PIE representation can be written in terms of explicit maps from the PDE parameters to PIE parameters. Lastly, we present numerical examples to demonstrate the application of the PIE representation by performing stability analysis of PDEs via convex optimization methods. 

Recently, a broad class of linear delayed and ODE-PDEs systems was shown to have an equivalent representation using Partial Integral Equations (PIEs). In this paper, we use this PIE representation, combined with algorithms for convex optimization of Partial Integral (PI) operators to bound the H2-norm for input-output systems of this class. Specifically, the methods proposed here apply to delayed and ODE-PDE systems (including delayed PDE systems) in one or two spatial variables where the disturbance does not enter through the boundary. For such systems, we define a notion of H2-norm using an initial state-to-output framework and show that this notion reduces to more traditional concepts under the assumption of existence of a strongly continuous semigroup. Next, we consider input-output systems for which there exists a PIE representation and for such systems show that computing a minimal upper bound on the H2-norm of delayed and PDE systems can be equivalently formulated as a convex optimization problem subject to linear PI operator inequalities (LPIs). We convert, then, these optimization problems to Semi-Definite Programming (SDP) problems using the PIETOOLS toolbox. Finally, we apply the results to several numerical examples – focusing on time-delay systems (TDS) for which comparable H2 approximation results are available in the literature. The numerical results demonstrate the accuracy of the computed upper bound on the H2-norm.