Stability. Outline Stability Sab Stability of Digital Systems. Stability for Continuous-time Systems. system is its stability:

Similar documents
State-Space Model. In general, the dynamic equations of a lumped-parameter continuous system may be represented by

TIME RESPONSE Introduction

Lecture 25 Outline: LTI Systems: Causality, Stability, Feedback

Analysis of Dynamic Systems

MODERN CONTROL SYSTEMS

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory

Chapter 2: Time-Domain Representations of Linear Time-Invariant Systems. Chih-Wei Liu

ME 3210 Mechatronics II Laboratory Lab 6: Second-Order Dynamic Response

ECE-314 Fall 2012 Review Questions

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002

Design of Controller for Robot Position Control

Economics 8723 Macroeconomic Theory Problem Set 3 Sketch of Solutions Professor Sanjay Chugh Spring 2017

Solutions to selected problems from the midterm exam Math 222 Winter 2015

Notes 03 largely plagiarized by %khc

Big O Notation for Time Complexity of Algorithms

The following Content. Objective. Introduction to Process Control. What have we talked in MM1-MM5? MM5? Software: Control Tutorials for Matlab

DETERMINATION OF PARTICULAR SOLUTIONS OF NONHOMOGENEOUS LINEAR DIFFERENTIAL EQUATIONS BY DISCRETE DECONVOLUTION

M. Rafeeyan. Keywords: MIMO, QFT, non-diagonal, control, uncertain

Hadamard matrices from the Multiplication Table of the Finite Fields

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4)

Math 2414 Homework Set 7 Solutions 10 Points

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws.

Review - Week 10. There are two types of errors one can make when performing significance tests:

EXISTENCE THEORY OF RANDOM DIFFERENTIAL EQUATIONS D. S. Palimkar

Automatic Control Systems

SHOCK AND VIBRATION RESPONSE SPECTRA COURSE Unit 21 Base Excitation Shock: Classical Pulse

u(t) Figure 1. Open loop control system

xp (X = x) = P (X = 1) = θ. Hence, the method of moments estimator of θ is

Electrical Engineering Department Network Lab.

An interesting result about subset sums. Nitu Kitchloo. Lior Pachter. November 27, Abstract

N! AND THE GAMMA FUNCTION

λiv Av = 0 or ( λi Av ) = 0. In order for a vector v to be an eigenvector, it must be in the kernel of λi

The Eigen Function of Linear Systems

Section 8 Convolution and Deconvolution

Theoretical Physics Prof. Ruiz, UNC Asheville, doctorphys on YouTube Chapter Q Notes. Laplace Transforms. Q1. The Laplace Transform.

CHAPTER 2 Quadratic diophantine equations with two unknowns

Fourier transform. Continuous-time Fourier transform (CTFT) ω ω

Lecture 15: Three-tank Mixing and Lead Poisoning

EECE 301 Signals & Systems Prof. Mark Fowler

S n. = n. Sum of first n terms of an A. P is

t = s D Overview of Tests Two-Sample t-test: Independent Samples Independent Samples t-test Difference between Means in a Two-sample Experiment

Extremal graph theory II: K t and K t,t

Ruled surfaces are one of the most important topics of differential geometry. The

LIMITS OF FUNCTIONS (I)

Suggested Solutions to Assignment 1 (REQUIRED)

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

10.3 Autocorrelation Function of Ergodic RP 10.4 Power Spectral Density of Ergodic RP 10.5 Normal RP (Gaussian RP)

e x x s 1 dx ( 1) n n!(n + s) + e s n n n=1 n!n s Γ(s) = lim

Manipulations involving the signal amplitude (dependent variable).

Vibration damping of the cantilever beam with the use of the parametric excitation

INVESTMENT PROJECT EFFICIENCY EVALUATION

arxiv:math/ v1 [math.fa] 1 Feb 1994

1 Notes on Little s Law (l = λw)

SUMMATION OF INFINITE SERIES REVISITED

Linear Time Invariant Systems

ECE 350 Matlab-Based Project #3

Comparison between Fourier and Corrected Fourier Series Methods

FIXED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE

Lecture 9: Polynomial Approximations

11. Adaptive Control in the Presence of Bounded Disturbances Consider MIMO systems in the form,

The universal vector. Open Access Journal of Mathematical and Theoretical Physics [ ] Introduction [ ] ( 1)

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

NEWTON METHOD FOR DETERMINING THE OPTIMAL REPLENISHMENT POLICY FOR EPQ MODEL WITH PRESENT VALUE

Calculus BC 2015 Scoring Guidelines

B. Maddah INDE 504 Simulation 09/02/17

2 f(x) dx = 1, 0. 2f(x 1) dx d) 1 4t t6 t. t 2 dt i)

The Signal, Variable System, and Transformation: A Personal Perspective

Time Dependent Queuing

12 Getting Started With Fourier Analysis

CS623: Introduction to Computing with Neural Nets (lecture-10) Pushpak Bhattacharyya Computer Science and Engineering Department IIT Bombay

CHATTERJEA CONTRACTION MAPPING THEOREM IN CONE HEPTAGONAL METRIC SPACE

Outline. simplest HMM (1) simple HMMs? simplest HMM (2) Parameter estimation for discrete hidden Markov models

, 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

Variational Iteration Method for Solving Differential Equations with Piecewise Constant Arguments

Communications II Lecture 4: Effects of Noise on AM. Professor Kin K. Leung EEE and Computing Departments Imperial College London Copyright reserved

PIECEWISE N TH ORDER ADOMIAN POLYNOMIAL STIFF DIFFERENTIAL EQUATION SOLVER 13

If boundary values are necessary, they are called mixed initial-boundary value problems. Again, the simplest prototypes of these IV problems are:

Homework 6 AERE331 Spring 2019 Due 4/24(W) Name Sec. 1 / 2

Math 6710, Fall 2016 Final Exam Solutions

EEC 483 Computer Organization

ECE 570 Session 7 IC 752-E Computer Aided Engineering for Integrated Circuits. Transient analysis. Discuss time marching methods used in SPICE

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

EECE 301 Signals & Systems Prof. Mark Fowler

Math-303 Chapter 7 Linear systems of ODE November 16, Chapter 7. Systems of 1 st Order Linear Differential Equations.

Chapter 9 - The Laplace Transform

Problems and Solutions for Section 3.2 (3.15 through 3.25)

6.003 Homework #5 Solutions

K3 p K2 p Kp 0 p 2 p 3 p

Pure Math 30: Explained!

Types Ideals on IS-Algebras

King Fahd University of Petroleum & Minerals Computer Engineering g Dept

ECE 145B / 218B, notes set 5: Two-port Noise Parameters

Digital Modulation Schemes

Modal Analysis of a Tight String

Sampling. AD Conversion (Additional Material) Sampling: Band limited signal. Sampling. Sampling function (sampling comb) III(x) Shah.

David Randall. ( )e ikx. k = u x,t. u( x,t)e ikx dx L. x L /2. Recall that the proof of (1) and (2) involves use of the orthogonality condition.

6.003: Signals and Systems

LINEAR APPROXIMATION OF THE BASELINE RBC MODEL JANUARY 29, 2013

RCT Worksheets/Quizzes 1.06 Radioactivity and Radioactive Decay

Solution of the Hyperbolic Partial Differential Equation on Graphs and Digital Spaces: a Klein Bottle a Projective Plane and a 4D Sphere

Transcription:

Oulie Sabiliy Sab Sabiliy of Digial Syem Ieral Sabiliy Exeral Sabiliy Example Roo Locu v ime Repoe Fir Orer Seco Orer Sabiliy e Jury e Rouh Crierio Example Sabiliy A very impora propery of a yamic yem i i abiliy: Ieral abiliy Exeral abiliy Ieral abiliy i cocere wih he repoe of all ieral variable a.k.a ae. Exeral abiliy uie he ipu-oupu p behaviour of a yamic yem. Chaper ME 54 Chaper ME 54 Exeral Sabiliy Mo commo efiiio of appropriae repoe for exeral abiliy i ha for every boue ipu BI, he yamic yem houl have a boue oupu BO. If hi coiio i aifie, he yem i ai o be BIBO able. Sabiliy for Coiuou-ime Syem A yem i aympoically i able if i reur o i iiial equilibrium ae afer he applicaio of a impule. Sabiliy i a propery of he yem a i i NO epee o a pecific ipu or iiial coiio. A coiuou-ime ime yem CS i able if a oly if he real par of i pole are all ricly egaive. Pole mu be locae a lef-ha ie of he -plae! Chaper ME 54 4 Chaper ME 54 5

Mappig bewee - a -plae Sabiliy of Corol Syem Coier he rafer fucio of he cloe-loop corol yem how: Le Sice σ ω where ˆ ˆ e σ ω σ ω e e e σ e{ co 444 ω 4 i 44 ω Raiu Depic ui circle Chaper ME 54 6 Y c p B R A where Hece, A um{ } N c p A a N... a a N Nr Nc A a ri pi ci pi ci 44 i 4 4 i 444444444 real pole where N N r N c i calle he yem orer. complex cougae pole Chaper ME 54 7 c p Sabiliy y Co For a BIBO able yem, all pole [i.e. roo of characeriic polyomial A] mu lie iie i he ui circle. Syem wih omia pole reiig iie lefha emi-circular regio exhibi force ocillaio a half he amplig frequecy. Very ueirable feaure! Example pe Deermie he abiliy of he followig yem: 7 Y.658 6.6 R.98.96 Soluio: 7.64.98 Rearrage he rafer fucio i power of o - : 7 7 Y.658 6.6.64 R.98.96.98 7 7 Chaper ME 54 8 Chaper ME 54 9

Soluio Roo Co Soluio Malab Scrip Co he characeriic polyomial A become A.98.96.98 Malab fucio roo come hay a hi poi: roo[ -.98.96 -.98] Hece, he roo of A are fou o be p.55; p,.967 ±.499499 Sice p i ouie he ui circle, he yem i UNSABLE! *** Defie he icree-ime ime F he la argume - ell Malab ha he amplig ime i upecifie. f[.658e-7 6.6e-7.64e-7],[ -.98.96 -.98],-; *** Calculae impule repoe for amplig ep... [y,k] impule,; *** Plo he reul emk,y,'.'; xlabel'ime Iex'; ylabel'impule Repoe'; Chaper ME 54 Chaper ME 54 Soluio Repoe Co epoe Impule Re.45.4.5..5..5..55 For a ui impule ipu, i oupu grow wihou a bou! 4 6 8 ime Iex Chaper ME 54 x eig Sabiliy Liear ime- ivaria Syem x Ieal Sampler y y* y y y o e he abiliy, a ui impule i applie o he yem a i ime repoe i oberve. Le u aume ha he impule repoe i alo ample for ake of argume. Chaper ME 54

ime Repoe of a Fir Orer CS ime Repoe v. Roo Locu of CS Le he rafer fucio of a fir orer yem be a y y y he yem ha oe real roo: p -a. Sice X L{δ}, Y X a 5 4 herefore, he ime repoe become y y y L { Y } L { } e a Chaper ME 54 4 Chaper ME 54 5 ime Repoe v. Roo Locu of DS Sice y e -a, he ample repoe y* become * y y k e Hece, i Z-raform lea o ak Y Z e e { y k } a a Noe ha X. he correpoig icree-ime yem ha oly oe pole a p e -a. Roo Locu of DS Co Im -plae y* y* y* y* Re y* Force ocillaio 5 4 wih f / are oberve i hi y* regio! - y* 5 4 Chaper ME 54 6 Chaper ME 54 7

-plae Aalyi of Seco Orer CS C. Co. ^ co he rafer fucio of a eco-orer yem i expree a Y X Kω ζω ω Kω ζω ω ζω ω where ζ i he ampig raio; ω i he aural frequecy; ω i calle ampe aural frequecy: ω ˆ ω ζ -plae Aalyi Co.56.8. If he complex cougae pole are locae alog he lie, he ampig raio oe o chage. If he pole are locae o a paricular emi-circle, he aural frequecy i kep coa. Chaper ME 54 8 Chaper ME 54 9 -plae Aalyi of a Seco Orer DS he rafer fucio of a eco-orer icree-ime yem become Y X e b b ζω ζω co ω e he roo of he characeriic polyomial are ζω p, e [co ω ± i ω ] -plae Aalyi Co ζ ω Mappig i lighly more complicae! Number alog he raial irecio how ecreaig ampig raio. Number ao alog he perimeer of he ui circle iicae coa aural frequecy. Chaper ME 54 Chaper ME 54

Malab Scrip o Prouce gri Sabiliy Crieria w lipace,pi,'; ea lipace,.9,'; p ero,; p p; qr-; Plo coa ampig coour for i : w*-eai/qr-eai^ ; w*-eai/qr-eai^ - ; p:,i exp; p:,i exp; e cloe all; plo[p p]; hol o xlabel'real'; ylabel'imagiary'; axi'quare' Plo coa w* coour ea lipace,,'; w lipacepi/,pi,'; for i : wi*-ea *qr-ea.^; wi*-ea - *qr-ea.^; p:,i exp; p:,i exp; e plo[p p]; gri Compario! Moer compuer program uch a Malab irecly give he pole of a paricular rafer fucio. Help u eermie he abiliy of ha yem righ away! Someime, he corol egieer wih o e for abiliy for a eire cla of yem: Coroller parameer gai are uually embee io he coefficie of he characeriic polyomial. he rage of corol parameer eig regio ha yiel a able yem are ough. I uch cae, abiliy e are uilie: Jury e w raform alog wih he Rouh-Hurwi Crierio Chaper ME 54 Chaper ME 54 Jury e Develope by Jury a Blachar i 96 for icree-ime yem. A Jury array i forme employig he coefficie of characeriic polyomial A. Mo of he array eleme are calculae via a umerical compuaio i proceure. Similar o he Rouh crierio, cerai eleme of he array mu be all o-egaive egaive o ge a able yem. Uforuaely, he calculaio of array eleme ca be quie complicae a eiou! Chaper ME 54 4 Rouh Crierio for Dicree-ime Syem o apply he well-kow Rouh Crierio, a coiuou-ime equivale for he icree-ime yem mu be obaie hrough ui raformaio: Some egieer replace by w a call he reulig operaio a w raform which ae he fac ha he equaio above o o irecly ake u back o he origial -plae! Ayhow, Y R F B' F' A' where F i he cloe-loop rafer fucio. he coefficie of A i ue o form Rouh array. Chaper ME 54 5

Example Example Co Deermie he abiliy of he followig yem wih ec uig Rouh crierio: Soluio: Y.967 F.484484 R.948 ui biliear raformaio lea o F F'.484.967.948 hu, F'.999 he characeriic polyomial become A'. Rouh Array: :. All fir-colum coefficie are : bigger ha ero yem i margially able. : ε Chaper ME 54 6 Chaper ME 54 7