System Control. Lesson #19a. BME 333 Biomedical Signals and Systems - J.Schesser

Similar documents
Introduction to Control Systems

ELEC 372 LECTURE NOTES, WEEK 4 Dr. Amir G. Aghdam Concordia University

State space systems analysis

Automatic Control Systems

Time Response. First Order Systems. Time Constant, T c We call 1/a the time constant of the response. Chapter 4 Time Response

CONTROL SYSTEMS. Chapter 7 : Bode Plot. 40dB/dec 1.0. db/dec so resultant slope will be 20 db/dec and this is due to the factor s

Lecture 30: Frequency Response of Second-Order Systems

ECEN620: Network Theory Broadband Circuit Design Fall 2014

Last time: Ground rules for filtering and control system design

Chapter (a) ζ. ω. 5 2 (a) Type 0 (b) Type 0 (c) Type 1 (d) Type 2 (e) Type 3 (f) Type 3. (g) type 2 (h) type (a) K G s.

The Performance of Feedback Control Systems

Lecture 11. Course Review. (The Big Picture) G. Hovland Input-Output Limitations (Skogestad Ch. 3) Discrete. Time Domain

EECE 301 Signals & Systems Prof. Mark Fowler

Chapter #5 EEE Control Systems

CONTROL ENGINEERING LABORATORY

Curve Sketching Handout #5 Topic Interpretation Rational Functions

ECM Control Engineering Dr Mustafa M Aziz (2013) SYSTEM RESPONSE

High-Speed Serial Interface Circuits and Systems. Lect. 4 Phase-Locked Loop (PLL) Type 1 (Chap. 8 in Razavi)

Brief Review of Linear System Theory

MAHALAKSHMI ENGINEERING COLLEGE TIRUCHIRAPALLI

STABILITY OF THE ACTIVE VIBRATION CONTROL OF CANTILEVER BEAMS

ELEC 372 LECTURE NOTES, WEEK 1 Dr. Amir G. Aghdam Concordia University

ECE 422 Power System Operations & Planning 6 Small Signal Stability. Spring 2015 Instructor: Kai Sun

Adaptive control design for a Mimo chemical reactor

Fig. 1: Streamline coordinates

Isolated Word Recogniser

Module 4: Time Response of discrete time systems Lecture Note 1

ELEC 372 LECTURE NOTES, WEEK 6 Dr. Amir G. Aghdam Concordia University

COMPARISONS INVOLVING TWO SAMPLE MEANS. Two-tail tests have these types of hypotheses: H A : 1 2

Reliable Decentralized PID Stabilization of MIMO Systems

EE Control Systems

M227 Chapter 9 Section 1 Testing Two Parameters: Means, Variances, Proportions

COMM 602: Digital Signal Processing

1the 1it is said to be overdamped. When 1, the roots of

Chapter #3 EEE Subsea Control and Communication Systems

1it is said to be overdamped. When 1, the roots of

x z Increasing the size of the sample increases the power (reduces the probability of a Type II error) when the significance level remains fixed.

Chapter 2 Feedback Control Theory Continued

Lecture 4. Chapter 11 Nise. Controller Design via Frequency Response. G. Hovland 2004

Chapter #2 EEE Subsea Control and Communication Systems

Last time: Completed solution to the optimum linear filter in real-time operation

u t u 0 ( 7) Intuitively, the maximum principles can be explained by the following observation. Recall

VIII. Interval Estimation A. A Few Important Definitions (Including Some Reminders)

SOLUTION: The 95% confidence interval for the population mean µ is x ± t 0.025; 49

Q-BINOMIALS AND THE GREATEST COMMON DIVISOR. Keith R. Slavin 8474 SW Chevy Place, Beaverton, Oregon 97008, USA.

6.003 Homework #3 Solutions

376 CHAPTER 6. THE FREQUENCY-RESPONSE DESIGN METHOD. D(s) = we get the compensated system with :

Statistical Inference Procedures

Answer: 1(A); 2(C); 3(A); 4(D); 5(B); 6(A); 7(C); 8(C); 9(A); 10(A); 11(A); 12(C); 13(C)

Analog Filter Design. Part. 3: Time Continuous Filter Implementation. P. Bruschi - Analog Filter Design 1

Performance-Based Plastic Design (PBPD) Procedure

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

Professor: Mihnea UDREA DIGITAL SIGNAL PROCESSING. Grading: Web: MOODLE. 1. Introduction. General information

Chapter 5. Root Locus Techniques

1. (a) If u (I : R J), there exists c 0 in R such that for all q 0, cu q (I : R J) [q] I [q] : R J [q]. Hence, if j J, for all q 0, j q (cu q ) =

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

Chap8 - Freq 1. Frequency Response

DISCRETE MELLIN CONVOLUTION AND ITS EXTENSIONS, PERRON FORMULA AND EXPLICIT FORMULAE

Chapter 4 : Laplace Transform

ECE 3510 Root Locus Design Examples. PI To eliminate steady-state error (for constant inputs) & perfect rejection of constant disturbances

Lab(8) controller design using root locus

MEM 355 Performance Enhancement of Dynamical Systems Root Locus Analysis

Assignment 1 - Solutions. ECSE 420 Parallel Computing Fall November 2, 2014

REVIEW OF SIMPLE LINEAR REGRESSION SIMPLE LINEAR REGRESSION

University of California at Berkeley College of Engineering Department of Electrical Engineering and Computer Sciences

Introduction to Signals and Systems, Part V: Lecture Summary

MAT 271 Project: Partial Fractions for certain rational functions

Feedback Control Systems (FCS)

Course Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B.

Series: Infinite Sums

THE CONCEPT OF THE ROOT LOCUS. H(s) THE CONCEPT OF THE ROOT LOCUS

EE 508 Lecture 6. Dead Networks Scaling, Normalization and Transformations

CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems

The Root Locus Method

IntroEcono. Discrete RV. Continuous RV s

Gain-scheduling of Acceleration Estimator for Low-velocity Measurement with Encoders

Erick L. Oberstar Fall 2001 Project: Sidelobe Canceller & GSC 1. Advanced Digital Signal Processing Sidelobe Canceller (Beam Former)

In algebra one spends much time finding common denominators and thus simplifying rational expressions. For example:

Z - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.

ECE 522 Power Systems Analysis II 3.2 Small Signal Stability

On the number of sums of three unit fractions

Digital Signal Processing, Fall 2010

MAS160: Signals, Systems & Information for Media Technology. Problem Set 5. DUE: November 3, (a) Plot of u[n] (b) Plot of x[n]=(0.

Math 213b (Spring 2005) Yum-Tong Siu 1. Explicit Formula for Logarithmic Derivative of Riemann Zeta Function

EE 508 Lecture 6. Scaling, Normalization and Transformation

Application of Kronecker Summation Method in Computation of Robustly Stabilizing PI Controllers for Interval Plants

Frequency Domain Filtering

Heat Equation: Maximum Principles

MA Lesson 26 Notes Graphs of Rational Functions (Asymptotes) Limits at infinity

Dorf, R.C., Wan, Z., Johnson, D.E. Laplace Transform The Electrical Engineering Handbook Ed. Richard C. Dorf Boca Raton: CRC Press LLC, 2000

Numerical Integration Formulas

Homework 12 Solution - AME30315, Spring 2013

100(1 α)% confidence interval: ( x z ( sample size needed to construct a 100(1 α)% confidence interval with a margin of error of w:

Special Notes: Filter Design Methods

2.004 Dynamics and Control II Spring 2008

Root Locus Contents. Root locus, sketching algorithm. Root locus, examples. Root locus, proofs. Root locus, control examples

a 1 = 1 a a a a n n s f() s = Σ log a 1 + a a n log n sup log a n+1 + a n+2 + a n+3 log n sup () s = an /n s s = + t i

MM1: Basic Concept (I): System and its Variables

TESTS OF SIGNIFICANCE

Active reduction of identified machine drive system vibrations in the form of multi-stage gear units

Transcription:

Sytem Cotrol Leo #9a 76

Sytem Cotrol Baic roblem Say you have a ytem which you ca ot alter but it repoe i ot optimal Example Motor cotrol for exokeleto Robotic cotrol roblem that ca occur Utable Traiet Repoe Ocillatio Overhoot Rie Time Output Error 77

Ucotrolled Sytem Uit Step Repoe G ope () 0 0 Output() 0 0 ( 7.)(.8) 7..8 3 3 0 0 0 0 0.05 (.8) 7. 7.(7..8) 0.03 ( 7.).8.8(.8 7.) 0.08 Output(t) (0.03e 0.08e.8t 0.05)u(t) Output value i i error: -0.05=95% errror Rie Time: ecod G ope () 0 00 Output() 0 00 ( 5 j3.)( 5 j3.) 5 j3. 5 j3. 3 3 0 00 0 00 0.005 ( 5 j3.) 5 j3.5 5 j3.(5 j3. 5 j3.) 0.007.8 Output(t) (0.0054co(3.t.8)e 5t 0.005)u(t) roblem Output value i i error:-0.005=99.5% error Overhoot & Ocillatio: 0.006-0.005=0.00 0% Rie time: ecod 78

Ucotrolled Sytem G ope () Output() a b c a b c r r Output( t) ( e e ) u( t) where root are real ad uequal rt rt 3 [( t ) e ] ut ( ) where root are real ad equal rr r 3 rt 3 [e co( t t 3 ut r j r r ) ] ( ) where root are imagiary ad uequal ; * Fial Output Value A t, Output( t) 3 0 a b c c Root b b 4ac r, a Let' chooe a, b, ad c, r For the uderdamped cae,, r j., 79

ID Cotrol dif()= Iput()-Output() Iput() + - Σ G ID () G ope () Output() Take our ytem, G ope (), add a ID cotroller, G ID (), i erie ad put them a uity egative feedback arragemet 80

ID Cotroller Elemet D D dif () Derivative dif () I Itegral I dif () + + + Σ ID output () ( I D)dif () D I dif () dif () roportioal 8

ID G ope () a b c G ID () I D D I G G CLOSED () ope ()G ID () H()G ope ()G ID () H() G ope ()G ID () a b c D I D I a 3 b c G CLOSED () D I a 3 b c D I a 3 b c D I a 3 b c D I 8

ID Cotrol oly G G CLOSED CLOSED CONTROLLER () a b c () Output() I D I 3 ( D) ( ) I D 0 ( ) a b c a b c r r ( ) ( ) ( rt r t ) ( ) where root are real ad uequ Output t e e u t 3 [( t ) e ] ut ( ) where root are real ad equal rr r rt 3 3 [ e co( t) ] u( t) where root are imagiary ad uequal r j; r r * t 3 al 83

ID Cotrol oly CONTROLLER G CLOSED () Output() ( ) a b c a b c r r ( ) 3 Fial Output Value A t, Output( t) a a b c c Root b b 4 a( c ) r, a Let' chooe a, b, ad c r 3 0 ( ) ( ), For the uderdamped cae,, r j ( ),, icreae the frequecy of damped ocillatio. For the critically damped or over-damped cae,, r, ( ), icreae oe of the pole ad may decreae rie time. 84

ID Cotrol oly G CLOSED () CONTROLLER a b (c ) Output() a b (c ) Output(t) [ e t co(t ) 3 ]u(t) where root are imagiary ad uequal r j;r r *. 85

G G CLOSED CLOSED Output() ID D CONTROL D CONTROLLER () a b c I 0 D I 3 ( D) ( ) I () a ( b ) ( c ) a ( b ) ( c ) D D 3 D D a b c r r D ( D) ( ) () ( rt rt Output t e e 3,, ) u( t) where root are real ad uequal [( t ) e ] ut ( ) where root are real ad equal rr r rt 3 [ e co( t) ] u( t) where root are imagiary ad uequal r j; r r * t 3 Fial Output Value A t, Output( t) a Root D 3 3 0 a ( b D) ( c ) c ( b D) ( b D) 4 a( c ) r, a Let' chooe a, b, ad c r r ( D) ( D) 4( ) D ( ) D D ( ) ( ) ha the potetial to make the traiet become critically dampled or overdamped. 3 86

ID D CONTROL G CLOSED () D CONTROLLER D a 3 (b D ) (c ) D a (b D ) (c ) Output() D a (b D ) (c ) 3 r r Output(t) ( e r t e r t 3 )u(t) where root are real ad uequal 87

ID I Cotrol Oly I CONTROLLER G CLOSED () D I a 3 (b D ) (c ) I D 0 I G CLOSED () a 3 b c I Output() I a 3 b c I 3 4 r r r 3 Output(t) ( e rt e rt 3 e r3t 4 )u(t) where root are real ad uequal [( t )e rt 3 e r 3 t 4 ]u(t) where root are real ad equal r r r [( t t 3 )e rt 4 ]u(t) where root are real ad equal r r r r 3 [ e t co(t ) 3 e r 3 t 4 ]u(t) where root are imagiary ad uequal r j;r r * Fial Output Value A t, Output(t) 4 I I 88

ID G CLOSED () a b c D I 3 D I Output() a b c r r r D I 3 D I rt rt rt 3 3 4 3 4 3 Output( t) ( e e e ) u( t) where root are real ad uequal rt rt 3 [( t ) e e ] ut ( ) where root are real ad equal rrr 3 4 [( t t ) e ] ut ( ) where root are real ad equal rr r r rt 3 4 3 l r j; r r * t rt 3 [e co( t ) e 3 4] ut ( ) where root are imagiary ad uequa Fial Output Value A t, Output( t) D I I 4 3 0 a b c D I I 89

ID 90

ID 9

Geeral Tip for Deigig a ID Cotroller. Obtai a ope-loop repoe ad determie what eed to be improved. Add a proportioal cotrol to improve the rie time 3. Add a derivative cotrol to improve the overhoot 4. Add a itegral cotrol to elimiate the teady-tate error Adjut each of, I, ad D util you obtai a deired overall repoe. It i ot eceary to implemet all three cotroller (proportioal, derivative, ad itegral) ito a igle ytem. eep the cotroller a imple a poible. Itroductio: ID Cotroller Deig, Uiverity of Michiga, http://ctm.egi.umich.edu/ctms/idex.php?example=itroductio&ectio=cotrolid 9

Matlab Cotrol Sytem Toolbox Trafer Fuctio tf(um,de) Cotruct a trafer fuctio Iput are array which repreet the coefficiet of the umerator ad deomiator of the deired trafer fuctio 3 5 3 4 um=[ 3]; deom=[5 0 4]; tf(um,deom) Trafer fuctio: + 3 --------------- 5 ^3 + + 4 93

Matlab Cotrol Sytem Toolbox Step Fuctio Repoe tep(tf,t) Calculate the repoe due to a tep fuctio Iput are deired trafer fuctio ad timeframe 5 4 um=[]; deom=[5 4]; tettf=tf(um,deom) time=0:0.:0; tep(tettf,time) Trafer fuctio: --------------- 5 ^ + + 4 94

Matlab Cotrol Sytem Toolbox Serie erie(tf,tf) Cotruct a trafer fuctio from trafer fuctio um=[]; deom=[5 4]; Iput are trafer fuctio TF 5 4 3 TF 5 3 4 tettf=tf(um,deom); um=[3 ]; deom=[5 0 4]; tettf=tf(um,deom); erie(tettf,tettf) Trafer fuctio: 3 + --------------------------------------------- 5 ^5 + 0 ^4 + 30 ^3 + 4 ^ + 6 + 6 95

Matlab Cotrol Sytem Toolbox Feedback feedback(opelooptf,feedbacktf) Cotruct a trafer fuctio of the feedback arragemet Iput are trafer fuctio of the opeloop ytem ad the feedback ytem; for uity feedback ue a the ecod iput. Opetf 5 4 Cloedloop 5 4 5 4 5 5 um=[]; deom=[5 4]; tettf=tf(um,deom) feedback(tettf,) Trafer fuctio: --------------- 5 ^ + + 5 96

Matlab Cotrol Sytem Toolbox Root root(array) Calculate the root of a polyomial Iput are array which repreet the coefficiet of the polyomial 5 3 4 root([5 0 4]) 0.3930 + 0.99i 0.3930-0.99i -0.7860 97

Matlab Cotrol Sytem Toolbox SISO ID GUI Tool 98

Homework Uig the Matlab cotrol ytem toolbox fuctio, chooe D,, ad I to yield a tep repoe that ha mall output error, o overhoot, ad a rie time of le tha ecod for the followig trafer fuctio: 000 Choe each ID eparately to how what it doe to the ytem. Show the root for each cae. Mot likely you ll have to iterate your awer to yield a good repoe. 99