Lecture 13: H(s) Poles-Zeros & BIBO Stability. (a) What are poles and zeros? Answer: HS math and calculus review.

Similar documents
STABILITY. Have looked at modeling dynamic systems using differential equations. and used the Laplace transform to help find step and impulse

Time Response of Systems

Physics 6303 Lecture 22 November 7, There are numerous methods of calculating these residues, and I list them below. lim

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014

Unit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace

Basic Procedures for Common Problems

ECE 3620: Laplace Transforms: Chapter 3:

EE 380. Linear Control Systems. Lecture 10

Time Response Analysis (Part II)

Math Lecture 3 Notes

9.5 The Transfer Function

GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM. COURSE: ECE 3084A (Prof. Michaels)

LECTURE 18: THE JOY OF CONVOLUTION

EEE 184: Introduction to feedback systems

MITOCW watch?v=poho4pztw78

Control Systems Design

Single-Time-Constant (STC) Circuits This lecture is given as a background that will be needed to determine the frequency response of the amplifiers.

Introduction to the z-transform

Identification Methods for Structural Systems

Raktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response

Calculus: What is a Limit? (understanding epislon-delta proofs)

Design and Optimization of Energy Systems Prof. C. Balaji Department of Mechanical Engineering Indian Institute of Technology, Madras

Control Systems Engineering (Chapter 2. Modeling in the Frequency Domain) Prof. Kwang-Chun Ho Tel: Fax:

ECEN 605 LINEAR SYSTEMS. Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability 1/27

20. The pole diagram and the Laplace transform

ECE-202 FINAL April 30, 2018 CIRCLE YOUR DIVISION

27. The pole diagram and the Laplace transform

SECTION 9.2: ARITHMETIC SEQUENCES and PARTIAL SUMS

Math 31 Lesson Plan. Day 2: Sets; Binary Operations. Elizabeth Gillaspy. September 23, 2011

ECE317 : Feedback and Control

EE-202 Exam II March 3, 2008

Problem set 5 solutions 1

Module 4. Related web links and videos. 1. FT and ZT

Notes for ECE-320. Winter by R. Throne

CS 124 Math Review Section January 29, 2018

6.241 Dynamic Systems and Control

Advanced Analog Building Blocks. Prof. Dr. Peter Fischer, Dr. Wei Shen, Dr. Albert Comerma, Dr. Johannes Schemmel, etc

Lecture 9 Time-domain properties of convolution systems

Module 4 : Laplace and Z Transform Lecture 36 : Analysis of LTI Systems with Rational System Functions

MATH 12 CLASS 5 NOTES, SEP

Chapter 6: The Laplace Transform. Chih-Wei Liu

Math Lecture 4 Limit Laws

An Introduction to Control Systems

Math Calculus f. Business and Management - Worksheet 12. Solutions for Worksheet 12 - Limits as x approaches infinity

CH.6 Laplace Transform

Solutions to Math 53 Math 53 Practice Final

EE Experiment 11 The Laplace Transform and Control System Characteristics

Biostatistics and Design of Experiments Prof. Mukesh Doble Department of Biotechnology Indian Institute of Technology, Madras

12/20/2017. Lectures on Signals & systems Engineering. Designed and Presented by Dr. Ayman Elshenawy Elsefy

Ordinary Differential Equations Prof. A. K. Nandakumaran Department of Mathematics Indian Institute of Science Bangalore

MATH 1231 MATHEMATICS 1B CALCULUS. Section 4: - Convergence of Series.

Math 3361-Modern Algebra Lecture 08 9/26/ Cardinality

Lecture 10: Proportional, Integral and Derivative Actions

Control Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

Problem Weight Score Total 100

MITOCW MITRES18_005S10_DiffEqnsMotion_300k_512kb-mp4

Control Systems. System response. L. Lanari

ECEN 420 LINEAR CONTROL SYSTEMS. Lecture 2 Laplace Transform I 1/52

MAT01A1: Precise Definition of a Limit and Continuity

EE/ME/AE324: Dynamical Systems. Chapter 7: Transform Solutions of Linear Models

Review of Linear Time-Invariant Network Analysis

Some of the different forms of a signal, obtained by transformations, are shown in the figure. jwt e z. jwt z e

GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM. COURSE: ECE 3084A (Prof. Michaels)

ECE 388 Automatic Control

L2 gains and system approximation quality 1

Just-in-Time-Teaching and other gadgets. Richard Cramer-Benjamin Niagara University

SIGNALS AND SYSTEMS LABORATORY 4: Polynomials, Laplace Transforms and Analog Filters in MATLAB

Mathematics Algebra. It is used to describe the world around us.

CHAPTER 6 VECTOR CALCULUS. We ve spent a lot of time so far just looking at all the different ways you can graph

Solution of ECE 315 Test #1 Su03

Acceleration 1-D Motion for Calculus Students (90 Minutes)

MITOCW watch?v=4x0sggrxdii

PROFESSOR: WELCOME BACK TO THE LAST LECTURE OF THE SEMESTER. PLANNING TO DO TODAY WAS FINISH THE BOOK. FINISH SECTION 6.5

Math 1B Final Exam, Solution. Prof. Mina Aganagic Lecture 2, Spring (6 points) Use substitution and integration by parts to find:

Data Mining Prof. Pabitra Mitra Department of Computer Science & Engineering Indian Institute of Technology, Kharagpur

MATH 115, SUMMER 2012 LECTURE 12

so mathematically we can say that x d [n] is a discrete-time signal. The output of the DT system is also discrete, denoted by y d [n].

Copyright c 2008 Kevin Long

MITOCW watch?v=pqkyqu11eta

MITOCW ocw f99-lec23_300k

Laplace Transform Part 1: Introduction (I&N Chap 13)

EA2.3 - Electronics 2 1

7.1 Introduction. Apago PDF Enhancer. Definition and Test Inputs. 340 Chapter 7 Steady-State Errors

Laplace Transforms and use in Automatic Control

Solve Wave Equation from Scratch [2013 HSSP]

ECE : Linear Circuit Analysis II

YOU CAN BACK SUBSTITUTE TO ANY OF THE PREVIOUS EQUATIONS

EEE105 Teori Litar I Chapter 7 Lecture #3. Dr. Shahrel Azmin Suandi Emel:

MITOCW watch?v=fkfsmwatddy

MITOCW watch?v=b1ekhyc9tto

Discrete and continuous dynamic systems

2. FUNCTIONS AND ALGEBRA

MITOCW watch?v=hdyabia-dny

Definition of the Laplace transform. 0 x(t)e st dt

Position and Displacement

Lecture 11 - Basic Number Theory.

In other words, we are interested in what is happening to the y values as we get really large x values and as we get really small x values.

MA 3280 Lecture 05 - Generalized Echelon Form and Free Variables. Friday, January 31, 2014.

Circuit Theory Prof. S.C. Dutta Roy Department of Electrical Engineering Indian Institute of Technology, Delhi

Q You mentioned that in complex analysis we study analytic functions, or, in another name, holomorphic functions. Pray tell me, what are they?

Transcription:

1. Introduction Lecture 13: H(s) Poles-Zeros & BIBO Stability (a) What are poles and zeros? Answer: HS math and calculus review. (b) What are nice inputs? Answer: nice inputs are bounded inputs; if you excite a circuit with a nice bounded input, like a step function, then one hopes to get a nice (bounded) output signal. Bounded input Bounded output (BIBO) stability. (c) Hey Prof Ray, I read this article about meta stable which is in some books. Answer: Yeah right! No such thing. A flawed concept made up by those who never studied stability. (d) Prove it Prof Ray you are always talking about proving things. Huh! Answer: the Tacoma Narrows bridge the one that did the twist and shout (girders screaming as they twisted) to the symphony of the Tacoma narrows wind that bridge like London bridge it tumbled down. So much for metastable; what do ya think? (e) OK so what are transfer functions good for? 2. General Structure of a Rational H(s). (a) Assumption: Generally, in 202, H(s) = n(s) d(s) meaning that deg[n(s)] deg[d(s)]. is assumed to be PROPER (b) Structure of the transfer function:

H(s) = K sm + b 1 s m 1 +...+ b m s n + a 1 b n 1 +...+ a 1 = K (s z 1 )(s z 2 )...(s z m ) (s p 1 )(s p 2 )...(s p n ) (a) z i is a finite zero of H(s) ; (b) p i is a finite pole of H(s) (c) if lim s H(s) = 0 then is an infinite zero of H(s) ; (d) n m is the number of infinite zeros of H(s) ; (e) dc gain is H(0) = K( 1) n m m z i i=1 n p i i=1. Example 1. Find H(s) from the pole-zero plot below assuming H(0) = 5.

Step 1. H(s) = K (s + jb)(s jb) (s + a + jb)(s + a jb) = K s 2 + b 2 (s + a) 2 + b 2 Step 2. H(0) = 5 = K b2 a 2 + b 2 K = 5 a 2 + b 2 b 2 3. BIBO stability: nice input leads ALWAYS to a nice output. (a) What is nice? Think about what you consider a nice person and what others may consider a nice person. This is difficult to pin down. In mathematics we often define nice by what is called a norm. Probably the most familiar norm to most students in the Euclidean norm of a vector: x 1 x = x 2 then x = x 2 2 1 + x 2 2, the length of the vector. For us, nice refers to a property of the input signal and output signal. For us, nice means bounded. A function f (t) is bounded if f (t) K < for all t. So if nice in always leads to nice out means bounded inputs lead to bounded outputs. so called BIBO stability. Yeah right Ray. So how much time do we have to check all those bounded inputs to see if they ALL lead to bounded outputs??? No way. I ain t got time for that. I ll spend the rest of my life catching my breath, letting it go. Good point. That is exactly the problem. The idea is NOT to check all those input-output pairs but to characterize the property of BIBO stability in terms of the transfer function H(s). BIBO stability is true if and only if ALL poles of H(s) lie in the open left half complex plane. It s all so simple now.

Example 2. H(s) = V out = 12 V in s 3. Suppose v in (t) = u(t) V. Then 12 V out = s(s 3) = 4 s 4 s 3 to infinity as t. implies v out (t) = 4u(t) 4e3t u(t) which blows up Conclusion: poles in the right half complex plane are very bad. Example 3. H(s) = V out V in = 12 s. Suppose v in (t) = u(t) V. Then V out = 12 t. s 2 implies v out (t) = 4r(t) = tu(t) which blows up to infinity as Conclusion: poles on the imaginary axis are bad. Example 4. H(s) = V out = 12 V in s + 3. Suppose v in (t) = u(t) V. Then 12 V out = s(s + 3) = 4 s 4 s + 3 implies v out (t) = 4u(t) + 4e 3t u(t) which exponentially goes to zero as t. Conclusion: poles in the left half complex plane are very good. Big Conclusion: A circuit is BIBO stable if and only if all poles of H(s) lie in the open left half complex plane.

Example 5. The bridge over the narrows problem: Example 6. Find the range of a for instability. Strategy: nodal analysis.

1. Output node: 1 V out + 1 ( s V out V 1) = 0. Solving for V 1 yields: (s +1)V out = V 1 2. The V 1 node. ( V 1 V in ) av out + 1 ( s V 1 V out ) = 0 From 1, upon substitution (s +1)V out av out +V out = V in 3. Conclusion: H(s) = V out V in = 1 s + 2 a which requires that a 2 < 0 a < 2 for BIBO stability.

Worksheet 1. (s 1) 2 (s + 2) 2 H(s) = 36 (s +1) 2 (s + 4)(s 3) 2 1. How many finite zeros? 2. Location of finite zeros? 3. Multiplicity of finite zeros? 4. How many infinite zeros? 5. How many finite poles? 6. Location of finite poles? 7. Multiplicity of finite poles? 8. True-False. The impulse response contains a term of the form Kδ (t). 9. True-False. The circuit/system described by the above transfer function is BIBO stable.

Worksheet 2. Given H(1) = 15 and the pole-zero plot below, find H(s). 1. H(s) = K 2. H(1) = 15 = K Therefore K = 3. True-False: The circuit/system described by the above transfer function is BIBO stable.

Appendix to Lecture 13. Example 1. Two Component Model of the Ingestion and Metabolism of a Drug. Compartment 1 Model dm 1 dt = K 1 m 1 + r Compartment 2 Model dm 2 dt = K 1 m 1 K 2 m 2 Assumptions 1. Output = m 2 (t) = mass of drug in bloodstream at t. 2. All IC s are zero at t = 0.

Objective: Construct H(s) and discuss. Step 1. sm 1 (s) = K 1 M 1 (s) + R(s) implies M 1 (s) = Step 2. sm 2 (s) = K 1 M 1 (s) K 2 M 2 (s) M 2 (s) = K 1 M 1 (s) s + K 2. 1 s + K 1 R(s) Step 3. H(s) = M 2 (s) R(s) = K 1 (s + K 1 )(s + K 2 ) Step 4. Impulse Response: h(t) = K 1 K 2 K 1 e K 1 t u(t) K 1 K 2 K 1 e K 2 t u(t) Step 5. Pole-zero plot of H(s). Suppose K 2 > K 1.