Intrinsic diculties in using the. control theory. 1. Abstract. We point out that the natural denitions of stability and

Similar documents
only nite eigenvalues. This is an extension of earlier results from [2]. Then we concentrate on the Riccati equation appearing in H 2 and linear quadr

A note on continuous behavior homomorphisms

PARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN. H.T. Banks and Yun Wang. Center for Research in Scientic Computation

A converse Lyapunov theorem for discrete-time systems with disturbances

1 Introduction We study innite-dimensional systems that can formally be written in the familiar form x 0 (t) = Ax(t) + Bu(t); y(t) = Cx(t) + Du(t); t

Discussion on: Measurable signal decoupling with dynamic feedforward compensation and unknown-input observation for systems with direct feedthrough

and the nite horizon cost index with the nite terminal weighting matrix F > : N?1 X J(z r ; u; w) = [z(n)? z r (N)] T F [z(n)? z r (N)] + t= [kz? z r

Chapter 7 Interconnected Systems and Feedback: Well-Posedness, Stability, and Performance 7. Introduction Feedback control is a powerful approach to o

u e G x = y linear convolution operator. In the time domain, the equation (2) becomes y(t) = (Ge)(t) = (G e)(t) = Z t G(t )e()d; and in either domains

The norms can also be characterized in terms of Riccati inequalities.

Vector Space Basics. 1 Abstract Vector Spaces. 1. (commutativity of vector addition) u + v = v + u. 2. (associativity of vector addition)

A NOTE ON MATRIX REFINEMENT EQUATIONS. Abstract. Renement equations involving matrix masks are receiving a lot of attention these days.

Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space

DISTANCE BETWEEN BEHAVIORS AND RATIONAL REPRESENTATIONS

arxiv: v3 [math.oc] 1 Sep 2018

Ole Christensen 3. October 20, Abstract. We point out some connections between the existing theories for

G : Statistical Mechanics Notes for Lecture 3 I. MICROCANONICAL ENSEMBLE: CONDITIONS FOR THERMAL EQUILIBRIUM Consider bringing two systems into

C.I.BYRNES,D.S.GILLIAM.I.G.LAUK O, V.I. SHUBOV We assume that the input u is given, in feedback form, as the output of a harmonic oscillator with freq

The model reduction algorithm proposed is based on an iterative two-step LMI scheme. The convergence of the algorithm is not analyzed but examples sho

H 1 optimisation. Is hoped that the practical advantages of receding horizon control might be combined with the robustness advantages of H 1 control.

ESC794: Special Topics: Model Predictive Control

CONTINUOUS TIME D=0 ZOH D 0 D=0 FOH D 0

BUMPLESS SWITCHING CONTROLLERS. William A. Wolovich and Alan B. Arehart 1. December 27, Abstract

L2 gains and system approximation quality 1

Control System. Contents

Digital Signal Processing, Homework 1, Spring 2013, Prof. C.D. Chung

Re-sampling and exchangeable arrays University Ave. November Revised January Summary

The Output Stabilization problem: a conjecture and counterexample

On at systems behaviors and observable image representations

Eigenvalue problems and optimization

Rearrangements and polar factorisation of countably degenerate functions G.R. Burton, School of Mathematical Sciences, University of Bath, Claverton D

1 The pendulum equation

Computability and Complexity

Functions: A Fourier Approach. Universitat Rostock. Germany. Dedicated to Prof. L. Berg on the occasion of his 65th birthday.

A Hilbert space approach to self-similar systems

University of California. Berkeley, CA fzhangjun johans lygeros Abstract

Experimental evidence showing that stochastic subspace identication methods may fail 1

Let T (N) be the algebra of all bounded linear operators of a Hilbert space L which leave invariant every subspace N in N, i.e., A T (N), AN N.

Control System Design

Null controllable region of LTI discrete-time systems with input saturation

A NEW APPROACH TO MIXED H 2 /H OPTIMAL PI/PID CONTROLLER DESIGN

QUANTITATIVE L P STABILITY ANALYSIS OF A CLASS OF LINEAR TIME-VARYING FEEDBACK SYSTEMS

EXTENDED ABSTRACT. 2;4 fax: , 3;5 fax: Abstract

INPUT-TO-STATE STABILITY OF EXPONENTIALLY. Lars Grune. Fachbereich Mathematik, AG 1.1. Johann Wolfgang Goethe-Universitat.

New concepts: Span of a vector set, matrix column space (range) Linearly dependent set of vectors Matrix null space

A Chaotic Phenomenon in the Power Swing Equation Umesh G. Vaidya R. N. Banavar y N. M. Singh March 22, 2000 Abstract Existence of chaotic dynamics in

INRIA Rocquencourt, Le Chesnay Cedex (France) y Dept. of Mathematics, North Carolina State University, Raleigh NC USA

Decentralized control with input saturation

Congurations of periodic orbits for equations with delayed positive feedback

Abstract An algorithm for Iterative Learning Control is developed based on an optimization principle which has been used previously to derive gradient

On the determinacy of repetitive structures

Information Structures Preserved Under Nonlinear Time-Varying Feedback

Richard DiSalvo. Dr. Elmer. Mathematical Foundations of Economics. Fall/Spring,

Abstract. Previous characterizations of iss-stability are shown to generalize without change to the

IMECE NEW APPROACH OF TRACKING CONTROL FOR A CLASS OF NON-MINIMUM PHASE LINEAR SYSTEMS

Input-Output Stability. A common task for an engineer is to design a system which reacts to stimuli

GATE EE Topic wise Questions SIGNALS & SYSTEMS

Discrete Mathematcs and Logic II.

To appear in IEEE Trans. on Automatic Control Revised 12/31/97. Output Feedback

Control Systems. EC / EE / IN. For

Eects of small delays on stability of singularly perturbed systems

Contents. 1 Introduction to Dynamics. 1.1 Examples of Dynamical Systems

An LQ R weight selection approach to the discrete generalized H 2 control problem

Optimal Polynomial Control for Discrete-Time Systems

The Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria

The Important State Coordinates of a Nonlinear System

Problem Set 5 Solutions 1

October 7, :8 WSPC/WS-IJWMIP paper. Polynomial functions are renable

Lifting to non-integral idempotents

Exam. 135 minutes, 15 minutes reading time

COPRIME FACTORIZATIONS AND WELL-POSED LINEAR SYSTEMS

Position in the xy plane y position x position

10 Transfer Matrix Models

Denition 2: o() is the height of the well-founded relation. Notice that we must have o() (2 ) +. Much is known about the possible behaviours of. For e

ẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)

2 C. A. Gunter ackground asic Domain Theory. A poset is a set D together with a binary relation v which is reexive, transitive and anti-symmetric. A s

Exhaustive Classication of Finite Classical Probability Spaces with Regard to the Notion of Causal Up-to-n-closedness

AN INTRODUCTION TO THE CONTROL THEORY

Time-Invariant Linear Quadratic Regulators!

Chapter 9 Observers, Model-based Controllers 9. Introduction In here we deal with the general case where only a subset of the states, or linear combin

growth rates of perturbed time-varying linear systems, [14]. For this setup it is also necessary to study discrete-time systems with a transition map

Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2015

Benchmark problems in stability and design of. switched systems. Daniel Liberzon and A. Stephen Morse. Department of Electrical Engineering

Math 1180, Notes, 14 1 C. v 1 v n v 2. C A ; w n. A and w = v i w i : v w = i=1

Math 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv

Disturbance Attenuation for a Class of Nonlinear Systems by Output Feedback

POLYNOMIAL EMBEDDING ALGORITHMS FOR CONTROLLERS IN A BEHAVIORAL FRAMEWORK

The Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria

Counting and Constructing Minimal Spanning Trees. Perrin Wright. Department of Mathematics. Florida State University. Tallahassee, FL

STABILITY OF INVARIANT SUBSPACES OF COMMUTING MATRICES We obtain some further results for pairs of commuting matrices. We show that a pair of commutin

Circuit depth relative to a random oracle. Peter Bro Miltersen. Aarhus University, Computer Science Department

Laplace Transforms and use in Automatic Control

Zeros. Carsten Scherer. Mathematisches Institut. Am Hubland. D-8700 Wurzburg

Module 4 : Laplace and Z Transform Problem Set 4

agree w/input bond => + sign disagree w/input bond => - sign

Problem Description The problem we consider is stabilization of a single-input multiple-state system with simultaneous magnitude and rate saturations,

New Lyapunov Krasovskii functionals for stability of linear retarded and neutral type systems

1 Introduction Travelling-wave solutions of parabolic equations on the real line arise in a variety of applications. An important issue is their stabi

QUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE. Denitizable operators in Krein spaces have spectral properties similar to those

Ridge analysis of mixture response surfaces

Transcription:

Intrinsic diculties in using the doubly-innite time axis for input-output control theory. Tryphon T. Georgiou 2 and Malcolm C. Smith 3 Abstract. We point out that the natural denitions of stability and causality in input-output control theory lead to certain inconsistencies when inputs and outputs are allowed to have support on the doublyinnite time-axis. In particular, linear time-invariant systems with right half plane poles cannot be considered to be both causal and stabilizable. In contrast, there is no such conict when the semi-innite time axis is used. Keywords. Linear systems, input-output stabilization, behaviour. Discussion Consider two systems P i (i = ; 2) dened in terms of convolution representations y(t) = Z? h i(t? )u()d = h i u where h (t) = e t for t 0 and zero otherwise, and h 2 (t) =?e t for t 0 and zero otherwise. Both systems have Laplace transfer functions =(s? ), but with diering regions of convergence. The rst system is unstable and Supported in part by the NSF and the AFOSR. 2 Department of Electrical Eng., University of Minnesota, Minneapolis, MN 55455, U.S.A. georgiou@ee.umn.edu, Tel. 62-625-3303, Fax: 62-625-4583 3 University of Cambridge, Department of Engineering, Cambridge CB2 PZ, U.K. mcs@eng.cam.ac.uk, Tel. +44-223-332475, Fax: +44-223-332662

causal and the second is stable and non-causal (in fact anticausal) according to the usual denitions. To view these systems abstractly in the input-output setting (or in the behavioural framework of [W]) we need to work out the system trajectories for signals in some function space, say L 2. This amounts to nding the graph, i.e. the set of input-output pairs in L 2. We mention that an explicit working out of such an approach for the semi-innite time-axis is given in [GS]. Below we consider the graphs of the two systems in the case of the doubly-innite time-axis. We consider rst P 2. From [R, p. 58, q. 4] f 2 L and g 2 L p implies f g 2 L p and kf gk p kfk kgk p : Thus, P 2 can be alternatively represented by the following graph on L 2 (?; ): ^GP 2 ;(?;) = L 2 (?j; j) () j? after transforming to the frequency domain. (As usual,^denotes the Fourier transform which maps L 2 (?; ) isometrically and isomorphically onto L 2 (?j; j) and L 2 [0; ) onto the Hardy space H 2 of the right half plane [R, Chapter 9]. Functions in H 2 can be extended analytically into the RHP by replacing j by s.) Now consider P. We will rst restrict the inputs to the space L 2 [0; ). Since y(t) = Z t 0 e t? u()d = e t Z t 0 e? u()d a necessary condition for y(t) 2 L 2 [0; ) is that Z t 0 e? u()d 0 as t, which is equivalent to he?t ; ui = 0. In fact, this is also sucient since then Z y(t) =?e t e? u()d = h 2 u: t Thus the domain of P in L 2 [0; ) is equal to the orthogonal complement of e?t C and, moreover, P coincides with P 2 on this domain. After taking 2

d + -? 6 u - P y y 2 C u 2? + d 2 Figure : Standard feedback conguration s? transforms, the orthogonal complement of C is H s+ s+ 2. Therefore the graph of P restricted to signals with support on [0; ) becomes (in the frequency domain): ^GP ;[0;) = s? s+ s+ H 2 =: GH 2 : Since P is shift invariant, the graph of P on L 2 (?; ) contains (again, expressed in the frequency domain) the innite union [ e st GH 2 : (2) T 0 The graph of P on L 2 (?; ) is actually slightly bigger than (2) (see below) but we will not need it for the next part of the discussion. We now turn our attention to the requirement of stabilizability in the feedback conguration of Fig.. Our denition of stability is the usual one: for all external L 2 disturbance inputs, there must be solutions of the feedback equations in L 2 so that the closed loop operators are norm bounded. If P is stabilized by some compensator C on L 2 (?; ), then it turns out that the graph of P must be closed. (This is a standard argument which proceeds as follows. Take a Cauchy sequence in the graph of P and set this equal to d d 2. Since by assumption the feedback equations have unique solutions 3

they must be given by: u = d, y = d 2, u 2 = 0, y 2 = 0. Now take the limit of the Cauchy sequence at the external inputs. Since the closed-loop operators are bounded, this pair of signals must also appear at u, y.) Thus, the graph of P must contain the closure of (2). But the closure of the innite union (2) equals the graph of P 2 Moreover, since P 2 is an anticausal operator, then there are L 2 (?; ) input-ouput pairs which satisfy the convolution representation for P 2 but not for P. Such a pair is: u(t) = e?t (t 0), 0 (t < 0) and y(t) =?e?jtj =2, which equals ^u ^y = s+ (s?)(s+) in the transform domain. Thus, it seems incorrect to close the graph of P and use the ordinary denition of input-output stability. On the other hand, stabilizability of P would require that any pair of possible disturbance signals d ; d 2 could act on the feedback system in Figure, and produce a bounded response. However, if a disturbance ^d ^d 2 = ^u^y ; with ^u; ^y as in (3), is allowed to act on the input ports in Figure, then there exists no solution which is consistent with the feedback equations and the integral representation of P. Of the following possible remedies none seems to be satisfactory:. To consider P to be non-stabilizable on L 2 (?; ), 2. To seek an alternative denition of closed-loop stability for the L 2 (?; ) case which would agree with the common belief that P is stabilized by proportional negative feedback of gain greater than one, 3. To identify the systems P ; P 2. It should be noted, that the use of extended spaces does not improve the situation, since the diculty is to determine the correct behaviour for signals in L 2 (?; ). Option () is a correct conclusion based on existing denitions. However, this would mean that the doubly-innite set-up is of limited interest for (3) 4

control purposes, since systems which in the usual sense are \open loop unstable" would have to be excluded. Option (2) does not seem to provide any satisfactory alternatives. It is, of course, possible to restrict attention to the subspace of L 2 consisting of functions which have support on some interval [T; ) for some arbitrary nite T. This would treat the graph of P precisely as (2) and would work with driving signals in S T 0 est H 2, which is not a closed subspace of L 2 2 (?j; j). However, this appears to be a rather cosmetic solution which more or less amounts to working on the half-line. A more natural avenue would be to consider the actual L 2 (?; ) graph of the convolution operator P. Similar reasoning to the above shows that this is the same as the graph of P 2, but with the restriction that the inputs satisfy: Z e? u()d = 0:? Again, this means that the graph P is not a closed subspace of L 2 2(?; ) and the problem is then to nd a suitable subspace for the external driving signals. The option of trying to work on some subspace of L 2 (?; ) with signals which \decay suciently fast" towards minus innity again does not seem to oer a satisfying resolution. Option (3), although unsatisfactory, is perhaps not as outrageous as would rst appear. If we consider the system represented by the dierential equation _y = y + u we can reproduce the trajectories of P by solving this equation forwards in time and those of P 2 by solving it backwards. This more or less amounts to abandonning any notion of causality (compare with the need to consider the anticausal trajectory (3) as a valid input-output pair for P ). However, this is not a natural option if the direction of time is well-dened (which is a basic assumption if we consider questions of control.) Conclusion The purpose of this note has been to point out certain features of inputoutput control theory on the doubly innite time axis which appear intrinsically unsatisfactory. The diculties are not of a mathematical nature a 5

consistent picture is obtained with a variety of denitions. The problem lies in trying to escape from conclusions which limit the engineering relevance of the theory, e.g. among the causal systems, only the stable ones are stabilizable. If the denitions lead to such a conclusion then it would not seem worthwhile to develop an elaborate theory of stabilization in that context. In contrast we remark that input-output systems theory on the doubly innite time axis does have many important uses e.g. in discussions of fundamental limitations in ltering imposed by causality conditions, or in system approximation using Hankel operators. References [R] W. Rudin, Real and Complex Analysis, 2nd Edition, McGraw-Hill, 982. [GS] T.T. Georgiou and M.C. Smith, Graphs, causality and stabilizability: linear, shift-invariant systems on L 2 [0; ), Mathematics of Control, Signals and Systems, 6(3): 95-223, 993. [W] J.C. Willems, Paradigms and puzzles in the theory of dynamical systems, IEEE Trans. on Automat. Contr., vol. 36, pp. 259{294, 99. 6