, the. L and the L. x x. max. i n. It is easy to show that these two norms satisfy the following relation: x x n x = (17.3) max

Similar documents
u(t) Figure 1. Open loop control system

To become more mathematically correct, Circuit equations are Algebraic Differential equations. from KVL, KCL from the constitutive relationship

Let. x y. denote a bivariate time series with zero mean.

CONTROL SYSTEMS. Chapter 10 : State Space Response

6.302 Feedback Systems Recitation : Phase-locked Loops Prof. Joel L. Dawson

EECE 301 Signals & Systems Prof. Mark Fowler

18.03SC Unit 3 Practice Exam and Solutions

Problem Set If all directed edges in a network have distinct capacities, then there is a unique maximum flow.

CHAPTER 7: SECOND-ORDER CIRCUITS

Introduction to SLE Lecture Notes

The Residual Graph. 12 Augmenting Path Algorithms. Augmenting Path Algorithm. Augmenting Path Algorithm

Chapter 7: Inverse-Response Systems

Laplace Transform. Inverse Laplace Transform. e st f(t)dt. (2)

The Residual Graph. 11 Augmenting Path Algorithms. Augmenting Path Algorithm. Augmenting Path Algorithm

1 Motivation and Basic Definitions

s-domain Circuit Analysis

Flow networks. Flow Networks. A flow on a network. Flow networks. The maximum-flow problem. Introduction to Algorithms, Lecture 22 December 5, 2001

Network Flows: Introduction & Maximum Flow

Sample Final Exam (finals03) Covering Chapters 1-9 of Fundamentals of Signals & Systems

Additional Methods for Solving DSGE Models

Design of Controller for Robot Position Control

6.8 Laplace Transform: General Formulas

Fractional Ornstein-Uhlenbeck Bridge

Control Systems -- Final Exam (Spring 2006)

Lecture 20: Riccati Equations and Least Squares Feedback Control

EE Control Systems LECTURE 2

CSC 364S Notes University of Toronto, Spring, The networks we will consider are directed graphs, where each edge has associated with it

CS4445/9544 Analysis of Algorithms II Solution for Assignment 1

Stability in Distribution for Backward Uncertain Differential Equation

18 Extensions of Maximum Flow

Generalized Orlicz Spaces and Wasserstein Distances for Convex-Concave Scale Functions

Research Article Existence and Uniqueness of Solutions for a Class of Nonlinear Stochastic Differential Equations

Chapter 6. Laplace Transforms

CHAPTER. Forced Equations and Systems { } ( ) ( ) 8.1 The Laplace Transform and Its Inverse. Transforms from the Definition.

FIXED POINTS AND STABILITY IN NEUTRAL DIFFERENTIAL EQUATIONS WITH VARIABLE DELAYS

6 December 2013 H. T. Hoang - www4.hcmut.edu.vn/~hthoang/ 1

NECESSARY AND SUFFICIENT CONDITIONS FOR LATENT SEPARABILITY

ANSWERS TO ODD NUMBERED EXERCISES IN CHAPTER

2. VECTORS. R Vectors are denoted by bold-face characters such as R, V, etc. The magnitude of a vector, such as R, is denoted as R, R, V

The generalized deniion allow u for inance o deal wih parial averaging of yem wih diurbance (for ome claical reul on parial averaging ee [3, pp

Control Systems. Lecture 9 Frequency Response. Frequency Response

Introduction to Congestion Games

Instrumentation & Process Control

Exponential Sawtooth

Convergence of the gradient algorithm for linear regression models in the continuous and discrete time cases

Linear Algebra Primer

Chapter 9 - The Laplace Transform

Main Reference: Sections in CLRS.

CHAPTER 7: UNCERTAINTY

Explicit form of global solution to stochastic logistic differential equation and related topics

Frequency Response. We now know how to analyze and design ccts via s- domain methods which yield dynamical information

Systems of nonlinear ODEs with a time singularity in the right-hand side

Chapter 6. Laplace Transforms

On the Exponential Operator Functions on Time Scales

Mathematische Annalen

CHAPTER 7. Definition and Properties. of Laplace Transforms

They were originally developed for network problem [Dantzig, Ford, Fulkerson 1956]

ESTIMATES FOR THE DERIVATIVE OF DIFFUSION SEMIGROUPS

Randomized Perfect Bipartite Matching

Optimal State-Feedback Control Under Sparsity and Delay Constraints

CHAPTER 3 SIGNALS & SYSTEMS. z -transform in the z -plane will be (A) 1 (B) 1 (D) (C) . The unilateral Laplace transform of tf() (A) s (B) + + (D) (C)

4. Advanced Stability Theory

Algorithms and Data Structures 2011/12 Week 9 Solutions (Tues 15th - Fri 18th Nov)

Group C*-algebras A.K. Recall 1 that if G is a nonempty set, the linear space

Please Complete Course Survey. CMPSCI 311: Introduction to Algorithms. Approximation Algorithms. Coping With NP-Completeness. Greedy Vertex Cover

FUZZY n-inner PRODUCT SPACE

Rough Paths and its Applications in Machine Learning

Physics 240: Worksheet 16 Name

Graphs III - Network Flow

5.2 GRAPHICAL VELOCITY ANALYSIS Polygon Method

Price of Stability and Introduction to Mechanism Design

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

NODIA AND COMPANY. GATE SOLVED PAPER Electrical Engineering SIGNALS & SYSTEMS. Copyright By NODIA & COMPANY

t )? How would you have tried to solve this problem in Chapter 3?

8.1. a) For step response, M input is u ( t) Taking inverse Laplace transform. as α 0. Ideal response, K c. = Kc Mτ D + For ramp response, 8-1

GLOBAL ANALYTIC REGULARITY FOR NON-LINEAR SECOND ORDER OPERATORS ON THE TORUS

CONTROL SYSTEMS. Chapter 3 Mathematical Modelling of Physical Systems-Laplace Transforms. Prof.Dr. Fatih Mehmet Botsalı

CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS

EECE 301 Signals & Systems Prof. Mark Fowler

Buckling of a structure means failure due to excessive displacements (loss of structural stiffness), and/or

CSE 521: Design & Analysis of Algorithms I

Discussion Session 2 Constant Acceleration/Relative Motion Week 03

Algorithmic Discrete Mathematics 6. Exercise Sheet

4/12/12. Applications of the Maxflow Problem 7.5 Bipartite Matching. Bipartite Matching. Bipartite Matching. Bipartite matching: the flow network

FLAT CYCLOTOMIC POLYNOMIALS OF ORDER FOUR AND HIGHER

Network Flow. Data Structures and Algorithms Andrei Bulatov

Approximate Controllability of Fractional Stochastic Perturbed Control Systems Driven by Mixed Fractional Brownian Motion

EE202 Circuit Theory II

A Theoretical Model of a Voltage Controlled Oscillator

Existence and Global Attractivity of Positive Periodic Solutions in Shifts δ ± for a Nonlinear Dynamic Equation with Feedback Control on Time Scales

Chapter 8 Objectives

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Identification of the Solution of the Burgers. Equation on a Finite Interval via the Solution of an. Appropriate Stochastic Control Problem

1 Adjusted Parameters

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and

Study of simple inductive-capacitive series circuits using MATLAB software package

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

5. Response of Linear Time-Invariant Systems to Random Inputs

Reminder: Flow Networks

Piecewise-Defined Functions and Periodic Functions

Transcription:

ecure 8 7. Sabiliy Analyi For an n dimenional vecor R n, he and he vecor norm are defined a: = T = i n i (7.) I i eay o how ha hee wo norm aify he following relaion: n (7.) If a vecor i ime-dependen, hen i runcaed funcional norm i defined a: ( τ ) = (7.3) τ Noe ha in (7.3), he vecor norm i ued inead of he vecor norm. For a mari A, vecor induced and mari norm will be denoed by A and A, correpondingly. Suppoe ha Λ = Λ K = Kˆ ˆ () (7.4) are he parameer bound for he Projecion Operaor ued in he adapive law (6.9). Fir, we ae and prove he following f abou inpu-o-ae abiliy of a linear imeinvarian (TI) yem. emma 7. There ei ricly poiive conan ( A, A) n ( ) =, he ae vecor () R k λ, uch ha aring a any iniial condiion of he TI yem wih a Hurwiz mari A and an eernal inpu m = A+ Bu (7.5) u R, aifie 73

B k A u λa () + () (7.6) Proof: Since A i Hurwiz hen here mu ei ricly poiive conan (, ) k λ uch A A ha e ka e λ. In hi cae, runcaed norm upper bound for he oluion of (7.5) can be found. () A ( τ ) A A A ( τ ) ( τ) = e + e Bu dτ λa λ ( τ) kae + k B e u( τ) dτ k + B e u d k + B e d u λa ( τ) λa ( τ) A ( τ ) τ A τ ( τ) τ B ka + u λa Thi complee he proof. A A (7.7) Ne, we find ufficien condiion o guaranee uniform ulimae boundedne (UUB) of all he ignal in he yem (5.), which i conrolled by a filered MRAC conroller, in he form of (6.5), (6.) and (6.9). Toward ha end, inroduce he o-called ideal erence model c = A + B r (7.8) and le e = (7.9) repreen he error beween he ual erence model (6.) and i ideal arge (7.8). Then he erence model error dynamic can be compued a: Uing (6.) and (7.9), give ( ) ˆ T e = A e + BΛˆ I K (7.) m m uc ( ) ˆ T ( ) e = A e + BΛˆ I K e + + e (7.) m m where ˆΛ, K ˆ, and e are UUB. We will ue Small ain Theorem ype argumen o find ufficien condiion for inpu-o-ae abiliy of (7.). Thi yem can be 74

repreened by he wo feedback-conneced yem H and H, a hown in he blockdiagram: H ( ) I m m Λˆ ( ) ( ) n n I A B e uc Kˆ T ( ) + e H Figure : Reference Model Error Dynamic According o he figure, hee wo yem are defined a and ( ) e = H u e = A e + BΛˆ I u (7.) c m m c ˆ T c c u = H u = K e + + e (7.3) Noe ha H i driven by he commanded conrol inpu ˆ T c while H i driven by he yem ae vecor u = K e + + e (7.4) ( ) = e + + e Bounded (7.5) e H denoe he gain of he yem H. One can eaily how ha he gain value i finie. In f, applying emma 7. o (7.) give 75

B Λ e () k e ( ) + ( ( ) I ) u A m m c λa B Λ k e ( ) + ( ) u + u A c c λ A uc k k e + B Λ + u A ( ) ( ) A c λa (7.6) where i he. The inequaliy in (7.6) give finie upper bound for he gain H of he yem H : gain of he proper and able filer H k ( ) A B Λ + (7.7) λa Furhermore, he relaion ˆ T uc = K ( e + + e) K e () () e() + + e (7.8) K e + K + e () ( ) imply ha he gain H of he yem H i alo finie: K (7.9) H Finally, uing Small ain Theorem allow o claim inpu-o-ae abiliy of he yem (7.), if < (7.) H H Becaue of (7.9), i i ufficien o chooe he erence model mari proper able filer ( ) uch ha he mall gain condiion A and he < (7.) H K i aified. 76

Remark 7. The upper bound in (7.) provide guideline for elecing he erence model mari. Baically, hee quaniie need o be choen o minimize he A and he filer gain H of he yem H in (7.). Baed on he above argumen and coninuing from he inequaliy (6.), one can ae and prove he following lemma. emma 7. e he proper able filer bound (, ), he Hurwiz mari A, and he Projecion Operaor K Λ be choen uch ha (7.) hold. Then all rajecorie of he cloed-loop yem (6.6), he erence model (6.), and he racking error dynamic are uniformly ulimaely bounded. Moreover, he racking error dynamic are globally aympoically able. Proof: Since he eernal command r i bounded hen i UUB. Alo, becaue of (6.), he racking error e i UUB. Furhermore, ince he aumed inequaliy (7.) implie (7.), he erence model error dynamic (7.) are inpu-o-ae able and, conequenly e i UUB. Hence, a i follow from (7.9), i UUB. Thu, he definiion of he racking error in (6.), implie ha i UUB. All hee f allow for he applicaion of he Barbala emma o (6.), proving UUB of all he ignal and aympoic convergence of he racking error e o he origin. The proof i complee. We now ummarize he filered MRAC deign and i abiliy properie. Theorem 7. The filered MRAC conrol archiecure coni of he filered feedback ignal (6.5), he erence model (6.), ˆ T u = K = u uc c = A + BΛˆ u u + B r c c and he Projecion Operaor baed adapive law (6.9). 77

T Proj (, ), ( ) T T ( ( c) ) Kˆ = Γ Kˆ e PB Kˆ = K ˆ Λ = Proj Λˆ, BPeu u Γ, Λ ˆ = I B If he conrol filer ( ), he erence mari (, ) A Λ m m K Λ are choen o aify he mall gain condiion in (7.) < H K, and he Projecion Operaor bound where H i he gain of he yem H in (7.) hen: ( ) e = H u e = A e + BΛˆ I u c m m c The yem ae aympoically rack he ae of he erence model (6.), when he laer i driven by any bounded erence command r( ). The yem ae vecor and he conrol inpu u are UUB. The correponding cloed-loop dynamic remain able wih all of i inernal ignal bounded, uniformly in ime. Eample 7. Conider open-loop unable ( A a ) b =Λ=. Alo uppoe ha = > calar TI yem dynamic and uppoe ha = a+ u A = a = a+ k < and ha a low-pa filer ( ) i choen in he form: u a = ( ) = u + a c 78

wih a >. Then he yem H dynamic are e = H u e = a e u e = u Conequenly, c c c + a + a + a H = ( + a )( + a ) The yem impule repone funcion h ( ) can be direcly compued via he invere aplace ranform: () h = ( + a )( + a ) = = a e a e a a + a a a a a a a ( ) + a In hi cae, he gain of he yem H i Thu, for a given a (7.) i aified: a a H = () h d = a e ae d a a a a, i i ufficien o chooe a uch ha he mall gain condiion a < a k Conequenly, he filer conan a mu be choen ufficienly large, o ha: a a > k or, equivalenly a > k + a = a+ k 3 79

Suppoe ha he plan conain a order uaor model, whoe ranfer funcion i: = a + a ( ) where a > i he invere of he uaor ime conan. In hi cae, he H dynamic become: a e = H u e = u ( + a )( + a )( + a ) c c The impule repone funcion of H i: a e a e a e h = a + + Conider he cae when a a a ( a a )( a a ) ( a a)( a a) ( a a)( a a) a < a < a Then he impule repone funcion become: a a a a e a e a e h = a + a a a a a a a a a a a a Now, he yem gain can be upper bounded: H = () h d a + + a a a a a a a a a a a a a = ( a a )( a a) and, hu H a ( a a )( a a) 8

Since our goal i o minimize he gain of H, we chooe a uch ha he upper bound above i minimized. The laer i equivalen o imizing he following epreion: I i eay o ee ha he opimal oluion i: ( a a )( a a) min a a a = + a In hi cae, he gain upper bound become: 8a H ( a a ) A he ame ime, he mall gain condiion require ha 8a < ( a ) a k When given he uaor dynamic, he above inequaliy conrain he elecion of he erence model dynamic which i achievable under he filered MRAC deign. 8