Input-Output (I/O) Stability

Total Page:16

File Type:pdf, Size:1020Kb

Input-Output (I/O) Stability I/O Stability ME689 Lecture Notes by B.Yao Input-Output (I/O) Stability -Stability of a System Outline: ¾ Introduction • White Boxes and Black Boxes • Input-Output Description ¾ Formalization of the Input-Output View • Signals and Signal Spaces • The Notions of Gain and Phase • The Notion of Passivity • The BIBO Stability Concept ¾ I/O Stability Theorems • The Small Gain Theorem • The Passivity Theorem • Positive Real Functions and Kalman-Yakubovich Lemma For LTI System 1 I/O Stability ME689 Lecture Notes by B.Yao 1. Introduction Systems can be modeled from two points of view: the internal (or state-space) approach and the external (or input-output) approach. The state-space approach takes the white box view of a system, and is based on a detailed description of the inner structure of the system. In the input-output approach, a system is considered to be a black box that transforms inputs to outputs. As a consequence, the system is modeled as an operator. The operator can be represented by either a verbal description of the input-output relationship, or a table look-up, or an abstract mathematical mapping that maps an input signal (or a function) in the input signal space (or a functional space) to an output signal (or a function) in the output signal space. 2. Formalization of I/O View 2.1 Signals and Signal Spaces 2.1.1 Single-Input Single-Output (SISO) Systems In SISO systems, the input and output signals can be described by real valued time functions, e.g., an input signal is denoted by a time-function u in R, i.e., u =Œ{}ut(), t R Note that to specify an input signal u, it is necessary to give its value at any time, not values at a point or an interval. The input and output signal spaces can be described by normed spaces; certain norms have to be defined to measure the size of a signal. For example, • The super norm or L• -norm u suput ( ) (I.1) • = t≥0 Application examples of using the super norm are the problem of finding maximum absolute tracking error and the problem of control saturation. • L2 -norm 1 • 2 u = Ê utdt2 () ˆ (I.2) 2 ˯Ú0 If we are interested in the energy of a signal, L2 -norm is a suitable norm to choose. • Lp -norm 1 • pp u = Ê ut() dtˆ , 1 £<•p (I.3) p ˯Ú0 Corresponding to different norms, different norm spaces are defined: • L• -space the set of all bounded functions, i.e., the set of all functions with a finite L• norm. LxxM=£: for some M >0 (I.4) • {}• 2 I/O Stability ME689 Lecture Notes by B.Yao • Lp -space the set of all functions with a finite Lp -norm., i.e., ϸ1 ÔÔ• p p LxxM=£: or Ê utdtM( )ˆ £ for some M >0 (I.5) p Ì˝p ËÚ0 ¯ Ó˛ÔÔ In general, these normed spaces may be too restrictive and may not include certain physically meaningful signals. For example, the L2 -space does not include the common sinusoid signals (why?). To be practically meaningful, the following extended space is defined to avoid this restrictiveness: Definition [ Extended Space] If X is a normed linear subspace of Y, then the extended space X e is the set XxYxXeT=Œ{}: Œ for any fixed T ≥0 , (I.6) where xT represents the truncation of x at T defined as Ïx(tt ) 0 ££T xtT ()= Ì (I.7) Ó0 tT> The extended L1 space is denoted as L1e , the extended L2 space as L2e , " . Example: u(t )=fiŒœ sinw t uL22e but uL y(tt )=fiŒœ yL••e but yL 2.1.2 Multi-Input Multi-Output (MIMO) Systems Multi inputs and multi outputs can be represented by vectors of time functions. For example, m inputs can be represented by a m-dimension of time functions u, i.e., u :0[ •Æ] Rm with È˘ut1() u()tRÍ˙# m (I.8) =ŒÍ˙ Î˚Í˙utm () Norms corresponding to (I.1)-(I.3) are defined as • The super norm u suput ( ) (I.9) • = t≥0 where u()t represents the norm of ut()in Rm space1. 1 ∞ m = The superior -norm of u(t) in R is u(t) ∞ max ui (t) . i 3 I/O Stability ME689 Lecture Notes by B.Yao • The L2 -norm ••2 u ==ut() dt uT () tutd () t (I.10) 22ÚÚ00 • The Lp -norm 1 • pp u = Ê ut() dtˆ , 1 £<•p (I.11) p ˯Ú0 where u()t represents the norm of ut()in Rm space2. m m m The corresponding norm spaces are represented by L• , L2 and Lp respectively. The extended m m m spaces are L•e , L2e and Lpe respectively. 2.1.3 Useful Facts The following lemmas state some useful facts about Lp -spaces. Lemma I.1 [Ref.1] (Holder’s Inequality) 11 If pq,1,Œ•, and +=1, then, fLŒ , g ŒL imply that fg ŒL , and [ ] pq p q 1 fg1 £ fpq g (I.12) i.e., 11 •••pqpq f() t g () t dt£ ÊÊ f () t dtˆˆ g () t dt (I.13) ÚÚÚ000Ë˯¯ Note: When pq==2 , the Holder’s inequality becomes the Cauchy-Schwartz inequality, i.e., 11 •••2222 f() t g () t dt£ ÊÊ f () t dtˆˆ g () t dt (I.14) ÚÚÚ000Ë˯¯ An extension of the Lemma I.1 to Lpe –space can be stated as Lemma I.2 11 If pq,1,Œ•, and +=1, then, fLŒ , g ŒL imply that fgLŒ , and [ ] pq pe qe 1e fg£ f g , "≥T 0 (I.15) ( )T 1 TTpq or 11 TTTpqpq ftgtdt() ()£ ÊÊ ft () dtˆˆ gt () dt "≥T 0 (I.16) ÚÚÚ000Ë˯¯ 1 m p p 2 The p-norm of u(t) in R m is u(t) = u (t) . p ∑ i i=1 4 I/O Stability ME689 Lecture Notes by B.Yao Lemma I.3 [Ref.1] (Minkowski Inequality) For p Œ•[1, ] , fg, Œ Lp imply that fgL+Œp , and fg+£ppp f + g (I.17) Note: (a) A Banach space is a normed space which is complete in the metric defined by its norm; this means that every Cauchy sequence is required to converge. It can be shown that Lp -space is a Banach space. (b) Power Signals: [Ref.5] The average power of a signal u over a time span T is 1 T u2 ()t dt (I.18) T Ú0 The signal u will be called a power signal if the limit of (I.18) exists as T Æ•, and then the square root of the average power will be denoted by pow(u ) : 1 2 ʈ1 T 2 pow(u )= Á˜ limutdt ( ) (I.19) ˯TÆ• T Ú0 (c) The sizes of common Lp spaces are graphically illustrated below [Ref.5] L pow 2 L∞ Note: L12«ÃLL• L1 Example: (verify by yourself) 1 ft()= fi f ()tLLŒ«2 , but f ()t œL1 , f ()tLŒ 1e 1+ t • Ï0 if t £ 0 Ô Ô 1 ft( )=<<Ì 0 t 1 fi f ()t ŒL1 , but f ()tLœ • , f ()tLœ 2 Ô t ÓÔ0 t > 1 5 I/O Stability ME689 Lecture Notes by B.Yao 2.2 Gain of A System 2.2.1 Definitions From an I/O point of view, a system is simply an operator S such that maps a input signal u into an output signal y. Thus, it can be represented by S : XXieÆ oe with y = Su (I.20) where the extended norm space X ie represents the input signal space (i.e., u Œ X ie ) and X oe is the output signal space. The gain of a linear system can thus be introduced through the induced norm of a linear operator defined below: Definition [Ref.2] [Gain of A Linear System] The gain g ()S of a linear system S is defined as Su g ()S = sup (I.21) uXŒ e u uπ0 where u is the input signal to the system. In other words, the gain g ()S is the smallest value such that SuS£ g ()u , "Œu X e (I.22) Nonlinear systems may have static offsets (i.e., Su π 0 when u = 0 ). To accommodate this particular phenomenon, a bias term is added in the above definition of gains: Definition [Gain of A Nonlinear System] The gain g ()S of a nonlinear system S is defined to be the smallest value such that Su£+g () S u b , "Œu X e (I.23) for some non-negative constant b . Definition [ Finite Gain Input-Output Stability] A nonlinear system S is called finite-gain Lp input-output stable if the gain g ()S defined in (I.23) is bounded (or finite), in which the Lp norm is used for input and output signals. Note: When p=∞, the above finite gain Lp stability, i.e., L∞ stability, results in bounded-input bounded-output (BIBO) stability. The converse is in general not true. For example, the static mapping y=u2 is BIBO stable but does not have a finite gain. 2.2.2 Linear Time Invariant Systems An LTI system with a TF G(s) can be described in input/output forms by Y ()sGsU= () ()s in s-domain t y()tgutgtud==- ( * )() (ttt ) ( ) in time-domain (I.24) Ú0 where g()tLGs= -1 {} () is the unit impulse response. 6 I/O Stability ME689 Lecture Notes by B.Yao Theorem I.1 [Ref.1] Let G(s) be a strictly proper rational function of s. (a) The system is stable (i.e., G(s) is analytic in Re[s ]≥ 0 or has all poles in LHP) if and only if gLŒ 1 .
Recommended publications
  • Control Theory
    Control theory S. Simrock DESY, Hamburg, Germany Abstract In engineering and mathematics, control theory deals with the behaviour of dynamical systems. The desired output of a system is called the reference. When one or more output variables of a system need to follow a certain ref- erence over time, a controller manipulates the inputs to a system to obtain the desired effect on the output of the system. Rapid advances in digital system technology have radically altered the control design options. It has become routinely practicable to design very complicated digital controllers and to carry out the extensive calculations required for their design. These advances in im- plementation and design capability can be obtained at low cost because of the widespread availability of inexpensive and powerful digital processing plat- forms and high-speed analog IO devices. 1 Introduction The emphasis of this tutorial on control theory is on the design of digital controls to achieve good dy- namic response and small errors while using signals that are sampled in time and quantized in amplitude. Both transform (classical control) and state-space (modern control) methods are described and applied to illustrative examples. The transform methods emphasized are the root-locus method of Evans and fre- quency response. The state-space methods developed are the technique of pole assignment augmented by an estimator (observer) and optimal quadratic-loss control. The optimal control problems use the steady-state constant gain solution. Other topics covered are system identification and non-linear control. System identification is a general term to describe mathematical tools and algorithms that build dynamical models from measured data.
    [Show full text]
  • Notes on Allpass, Minimum Phase, and Linear Phase Systems
    DSPDSPDSPDSPDSPDSPDSPDSPDSPDSPDSP SignalDSP Processing (18-491 ) Digital DSPDSPDSPDSPDSPDSPDSP/18-691 S pring DSPDSP DSPSemester,DSPDSPDSP 2021 Allpass, Minimum Phase, and Linear Phase Systems I. Introduction In previous lectures we have discussed how the pole and zero locations determine the magnitude and phase of the DTFT of an LSI system. In this brief note we discuss three special cases: allpass, minimum phase, and linear phase systems. II. Allpass systems Consider an LSI system with transfer function z−1 − a∗ H(z) = 1 − az−1 Note that the system above has a pole at z = a and a zero at z = (1=a)∗. This means that if the pole is located at z = rejθ, the zero would be located at z = (1=r)ejθ. These locations are illustrated in the figure below for r = 0:6 and θ = π=4. Im[z] Re[z] As usual, the DTFT is the z-transform evaluated along the unit circle: j! H(e ) = H(z)jz=0 18-491 ZT properties and inverses Page 2 Spring, 2021 It is easy to obtain the squared magnitude of the frequency response by multiplying H(ej!) by its complex conjugate: 2 (e−j! − re−j!)(ej! − rej!) H(ej!) = H(ej!)H∗(ej!) = (1 − rejθe−j!)(1 − re−jθej!) (1 − rej(θ−!) − re−j(θ−!) + r2) = = 1 (1 − rej(θ−!) − re−j(θ−!) + r2) Because the squared magnitude (and hence the magnitude) of the transfer function is a constant independent of frequency, this system is referred to as an allpass system. A sufficient condition for the system to be allpass is for the poles and zeros to appear in \mirror-image" locations as they do in the pole-zero plot on the previous page.
    [Show full text]
  • EECS 451 LECTURE NOTES BASIC SYSTEM PROPERTIES What: Input X[N] → |SYSTEM| → Y[N] Output
    EECS 451 LECTURE NOTES BASIC SYSTEM PROPERTIES What: Input x[n] → |SYSTEM| → y[n] output. Why: Design the system to filter input x[n]. DEF: A system is LINEAR if these two properties hold: 1. Scaling: If x[n] → |SYS| → y[n], then ax[n] → |SYS| → ay[n] for: any constant a. NOT true if a varies with time (i.e., a[n]). 2. Superposition: If x1[n] → |SYS| → y1[n] and x2[n] → |SYS| → y2[n], Then: (ax1[n]+ bx2[n]) → |SYS| → (ay1[n]+ by2[n]) for: any constants a, b. NOT true if a or b vary with time (i.e., a[n], b[n]). EX: y[n]=3x[n−2]; y[n]= x[n+1]−nx[n]+2x[n−1]; y[n] = sin(n)x[n]. NOT: y[n]= x2[n]; y[n] = sin(x[n]); y[n]= |x[n]|; y[n]= x[n]/x[n − 1]. NOT: y[n]= x[n] + 1 (try it). This is called an affine system. HOW: If any nonlinear function of x[n], not linear. Nonlinear of just n OK. DEF: A system is TIME-INVARIANT if this property holds: If x[n] → |SYS| → y[n], then x[n − N] → |SYS| → y[n − N] for: any integer time delay N. NOT true if N varies with time (e.g., N(n)). EX: y[n]=3x[n − 2]; y[n] = sin(x[n]); y[n]= x[n]/x[n − 1]. NOT: y[n]= nx[n]; y[n]= x[n2]; y[n]= x[2n]; y[n]= x[−n]. HOW: If n appears anywhere other than in x[n], not time-invariant.
    [Show full text]
  • Method for Undershoot-Less Control of Non- Minimum Phase Plants Based on Partial Cancellation of the Non-Minimum Phase Zero: Application to Flexible-Link Robots
    Method for Undershoot-Less Control of Non- Minimum Phase Plants Based on Partial Cancellation of the Non-Minimum Phase Zero: Application to Flexible-Link Robots F. Merrikh-Bayat and F. Bayat Department of Electrical and Computer Engineering University of Zanjan Zanjan, Iran Email: [email protected] , [email protected] Abstract—As a well understood classical fact, non- minimum root-locus method [11], asymptotic LQG theory [9], phase zeros of the process located in a feedback connection waterbed effect phenomena [12], and the LTR problem cannot be cancelled by the corresponding poles of controller [13]. In the field of linear time-invariant (LTI) systems, since such a cancellation leads to internal instability. This the source of all of the above-mentioned limitations is that impossibility of cancellation is the source of many the non-minimum phase zero of the given process cannot limitations in dealing with the feedback control of non- be cancelled by unstable pole of the controller since such a minimum phase processes. The aim of this paper is to study cancellation leads to internal instability [14]. the possibility and usefulness of partial (fractional-order) cancellation of such zeros for undershoot-less control of During the past decades various methods have been non-minimum phase processes. In this method first the non- developed by researchers for the control of processes with minimum phase zero of the process is cancelled to an non-minimum phase zeros (see, for example, [15]-[17] arbitrary degree by the proposed pre-compensator and then and the references therein for more information on this a classical controller is designed to control the series subject).
    [Show full text]
  • Digital Signal Processing Lecture Outline Discrete-Time Systems
    Lecture Outline ESE 531: Digital Signal Processing ! Discrete Time Systems ! LTI Systems ! LTI System Properties Lec 3: January 17, 2017 ! Difference Equations Discrete Time Signals and Systems Penn ESE 531 Spring 2017 - Khanna Penn ESE 531 Spring 2017 - Khanna 2 Discrete Time Systems Discrete-Time Systems Penn ESE 531 Spring 2017 - Khanna 4 System Properties Examples ! Causality ! Causal? Linear? Time-invariant? Memoryless? " y[n] only depends on x[m] for m<=n BIBO Stable? ! Linearity ! Time Shift: " Scaled sum of arbitrary inputs results in output that is a scaled sum of corresponding outputs " x[n]y[n =] x[n-m]= x[n − m] " Ax [n]+Bx [n] # Ay [n]+By [n] 1 2 1 2 ! Accumulator: ! Memoryless " y[n] n " y[n] depends only on x[n] y[n] = x[k] ! Time Invariance ∑ k=−∞ " Shifted input results in shifted output " x[n-q] # y[n-q] ! Compressor (M>1): ! BIBO Stability y[n] x[Mn] " A bounded input results in a bounded output (ie. max signal value = exists for output if max ) Penn ESE 531 Spring 2017 - Khanna 5 Penn ESE 531 Spring 2017 - Khanna 6 1 Non-Linear System Example Spectrum of Speech ! Median Filter " y[n]=MED{x[n-k], …x[n+k]} Speech " Let k=1 " y[n]=MED{x[n-1], x[n], x[n+1]} Corrupted Speech Penn ESE 531 Spring 2017 - Khanna 7 Penn ESE 531 Spring 2017 - Khanna 8 Low Pass Filtering Low Pass Filtering Corrupted Speech LP-Filtered Speech Penn ESE 531 Spring 2017 - Khanna 9 Penn ESE 531 Spring 2017 - Khanna 10 Median Filtering LTI Systems Corrupted Speech Med-Filter Speech Penn ESE 531 Spring 2017 - Khanna 11 Penn ESE 531 Spring 2017 - Khanna
    [Show full text]
  • Feed-Forward Compensation of Non-Minimum Phase Systems
    Wright State University CORE Scholar Browse all Theses and Dissertations Theses and Dissertations 2018 Feed-Forward Compensation of Non-Minimum Phase Systems Venkatesh Dudiki Wright State University Follow this and additional works at: https://corescholar.libraries.wright.edu/etd_all Part of the Electrical and Computer Engineering Commons Repository Citation Dudiki, Venkatesh, "Feed-Forward Compensation of Non-Minimum Phase Systems" (2018). Browse all Theses and Dissertations. 2200. https://corescholar.libraries.wright.edu/etd_all/2200 This Thesis is brought to you for free and open access by the Theses and Dissertations at CORE Scholar. It has been accepted for inclusion in Browse all Theses and Dissertations by an authorized administrator of CORE Scholar. For more information, please contact [email protected]. FEED-FORWARD COMPENSATION OF NON-MINIMUM PHASE SYSTEMS A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Electrical Engineering by VENKATESH DUDIKI BTECH, SRM University, 2016 2018 Wright State University Wright State University GRADUATE SCHOOL November 30, 2018 I HEREBY RECOMMEND THAT THE THESIS PREPARED UNDER MY SUPER- VISION BY Venkatesh Dudiki ENTITLED FEED-FORWARD COMPENSATION OF NON-MINIMUM PHASE SYSTEMS BE ACCEPTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF Master of Science in Electrical En- gineering. Pradeep Misra, Ph.D. THESIS Director Brian Rigling, Ph.D. Chair Department of Electrical Engineering College of Engineering and Computer Science Committee on Final Examination Pradeep Misra, Ph.D. Luther Palmer,III, Ph.D. Kazimierczuk Marian K, Ph.D Barry Milligan, Ph.D Interim Dean of the Graduate School ABSTRACT Dudiki, Venkatesh.
    [Show full text]
  • MTHE/MATH 332 Introduction to Control
    1 Queen’s University Mathematics and Engineering and Mathematics and Statistics MTHE/MATH 332 Introduction to Control (Preliminary) Supplemental Lecture Notes Serdar Yuksel¨ April 3, 2021 2 Contents 1 Introduction ....................................................................................1 1.1 Introduction................................................................................1 1.2 Linearization...............................................................................4 2 Systems ........................................................................................7 2.1 System Properties...........................................................................7 2.2 Linear Systems.............................................................................8 2.2.1 Representation of Discrete-Time Signals in terms of Unit Pulses..............................8 2.2.2 Linear Systems.......................................................................8 2.3 Linear and Time-Invariant (Convolution) Systems................................................9 2.4 Bounded-Input-Bounded-Output (BIBO) Stability of Convolution Systems........................... 11 2.5 The Frequency Response (or Transfer) Function of Linear Time-Invariant Systems..................... 11 2.6 Steady-State vs. Transient Solutions............................................................ 11 2.7 Bode Plots................................................................................. 12 2.8 Interconnections of Systems..................................................................
    [Show full text]
  • EE 225A Digital Signal Processing Supplementary Material 1. Allpass
    EE 225A — DIGITAL SIGNAL PROCESSING — SUPPLEMENTARY MATERIAL 1 EE 225a Digital Signal Processing Supplementary Material 1. Allpass Sequences A sequence ha()n is said to be allpass if it has a DTFT that satisfies jw Ha()e = 1 . (1) Note that this is true iff jw * jw Ha()e Ha()e = 1 (2) which in turn is true iff * ha()n *ha(–n)d= ()n (3) which in turn is true iff * 1 H ()z H æö---- = 1. (4) a aèø* z For rational Z transforms, the latter relation implies that poles (and zeros) of Ha()z are cancelled * 1 * 1 by zeros (and poles) of H æö---- . In other words, if H ()z has a pole at zc= , then H æö---- must aèø a aèø z* z* have have a zero at the same zc= , in order for (4) to be possible. Note further that a zero of * 1 1 H æö---- at zc= is a zero of H ()z at z = ----- . Consequently, (4) implies that poles (and zeros) aèø a z* c* PROF. EDWARD A. LEE, DEPARTMENT OF EECS, UNIVERSITY OF CALIFORNIA, BERKELEY EE 225A — DIGITAL SIGNAL PROCESSING — SUPPLEMENTARY MATERIAL 2 have zeros (and poles) in conjugate-reciprocal locations, as shown below: 1/c* c Intuitively, since the DTFT is the Z transform evaluated on the unit circle, the effect on the mag- nitude due to any pole will be canceled by the effect on the magnitude of the corresponding zero. Although allpass implies conjugate-reciprocal pole-zero pairs, the converse is not necessarily true unless the appropriate constant multiplier is selected to get unity magnitude.
    [Show full text]
  • En-Sci-175-01
    UNCLASSIFIED/UNLIMITED Performance Limits in Control with Application to Communication Constrained UAV Systems Prof. Dr. Richard H. Middleton ARC Centre for Complex Dynamic Systems and Control The University of Newcastle, Callaghan NSW, Australia, 2308 [email protected] ABSTRACT The study of performance limitations in feedback control systems has a long history, commencing with Bode’s gain-phase and sensitivity integral formulae. This area of research focuses on questions of what aspects of closed loop performance can we reasonably demand, and what are the trade-offs or costs of achieving this performance. In particular, it considers concepts related to the inherent costs of stabilisation, limitations due to ‘non-minimum phase’ behaviour, time delays, bandwidth limitations etc. There are a range of results developed in the 90’s for multivariable linear systems, together with extensions to sampled data, non-linear systems. More recently, this work has been extended to consider communication limitations, such as bit-rate or signal to noise ratio limitations and quantization effects. This paper presents: (i) A review of work on performance limitations, including recent results; (ii) Results on impacts of communication constraints on feedback control performance; and (iii) Implications of these results on control of distributed UAVs in formation. 1.0 INTRODUCTION In many control applications it is crucial that the factors limiting the performance of the system are identified. For example, closed loop control performance may be limited due to the fact that the actuators are too slow, or too inaccurate (perhaps due to complex nonlinear effects, such as backlash and hysteresis), to support the desired objectives.
    [Show full text]
  • A Note on BIBO Stability
    1 A Note on BIBO Stability Michael Unser Abstract—The statements on the BIBO stability of continuous- imposes that the result of the convolution be continuous, anda time convolution systems found in engineering textbooks are second (Theorem 3) that characterizes the BIBO-stable filters often either too vague (because of lack of hypotheses) or in full generality. The mathematical derivations are presented mathematically incorrect. What is more troubling is that they usually exclude the identity operator. The purpose of this note in Section IV, where we also make the connection with known is to clarify the issue while presenting some fixes. In particular, results in harmonic analysis. we show that a linear shift-invariant system is BIBO-stable in We like to mention a similar clarification effort by Hans Fe- L the ∞-sense if and only if its impulse response is included in ichtinger, who proposes to limit the framework to convolution the space of bounded Radon measures, which is a superset of operators that are operating on C0(R) (a well-behaved subclass L1(R) (Lebesgue’s space of absolutely integrable functions). As we restrict the scope of this characterization to the convolution of bounded functions) in order to avoid pathologies [3]. This operators whose impulse response is a measurable function, we is another interesting point of view that is complementary to recover the classical statement. ours, as discussed in Section IV. II. BIBO STABILITY:THE CLASSICAL FORMULATION I. INTRODUCTION The convolution of two functions h,f : R → R is the A statement that is made in most courses on the theory of function usually specified by linear systems as well as in the English version of Wikipedia1 △ is that a convolution operator is stable in the BIBO sense t 7→ (h ∗ f)(t) = h(τ)f(t − τ)dτ (1) (bounded input and bounded output) if and only if its impulse ZR response is absolutely summable/integrable.
    [Show full text]
  • Controllability, Observability, Realizability, and Stability of Dynamic Linear Systems
    Electronic Journal of Differential Equations, Vol. 2009(2009), No. 37, pp. 1–32. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu CONTROLLABILITY, OBSERVABILITY, REALIZABILITY, AND STABILITY OF DYNAMIC LINEAR SYSTEMS JOHN M. DAVIS, IAN A. GRAVAGNE, BILLY J. JACKSON, ROBERT J. MARKS II Abstract. We develop a linear systems theory that coincides with the ex- isting theories for continuous and discrete dynamical systems, but that also extends to linear systems defined on nonuniform time scales. The approach here is based on generalized Laplace transform methods (e.g. shifts and con- volution) from the recent work [13]. We study controllability in terms of the controllability Gramian and various rank conditions (including Kalman’s) in both the time invariant and time varying settings and compare the results. We explore observability in terms of both Gramian and rank conditions and establish related realizability results. We conclude by applying this systems theory to connect exponential and BIBO stability problems in this general setting. Numerous examples are included to show the utility of these results. 1. Introduction In this paper, our goal is to develop the foundation for a comprehensive linear systems theory which not only coincides with the existing canonical systems theories in the continuous and discrete cases, but also to extend those theories to dynamical systems with nonuniform domains (e.g. the quantum time scale used in quantum calculus [9]). We quickly see that the standard arguments on R and Z do not go through when the graininess of the underlying time scale is not uniform, but we overcome this obstacle by taking an approach rooted in recent generalized Laplace transform methods [13, 19].
    [Show full text]
  • Supplementary Notes for ELEN 4810 Lecture 13 Analysis of Discrete-Time LTI Systems
    Supplementary Notes for ELEN 4810 Lecture 13 Analysis of Discrete-Time LTI Systems John Wright Columbia University November 21, 2016 Disclaimer: These notes are intended to be an accessible introduction to the subject, with no pretense at com- pleteness. In general, you can find more thorough discussions in Oppenheim’s book. Please let me know if you find any typos. Reading suggestions: Oppenheim and Schafer Sections 5.1-5.7. In this lecture, we discuss how the poles and zeros of a rational Z-transform H(z) combine to shape the frequency response H(ej!). 1 Frequency Response, Phase and Group Delay Recall that if x is the input to an LTI system with impulse response h[n], we have 1 X y[n] = x[k]h[n − k]; (1.1) k=−∞ a relationship which can be studied through the Z transform identity Y (z) = H(z)X(z) (1.2) (on ROC fhg \ ROC fxg), or, if h is stable, through the DTFT identity Y (ej!) = H(ej!)X(ej!): (1.3) This relationship implies that magnitudes multiply and phases add: jY (ej!)j = jH(ej!)jjX(ej!)j; (1.4) j! j! j! \Y (e ) = \H(e ) + \X(e ): (1.5) 1 The phase \· is only defined up to addition by an integer multiple of 2π. In the next two lectures, we will occasionally need to be more precise about the phase. The text uses the notation ARG[H(ej!)] j! for the unique value \H(e ) + 2πk lying in (π; π]: −π < ARG[H] ≤ π: (1.6) This gives a well-defined way of removing the 2π ambiguity in the phase.
    [Show full text]