# Control Systems! Copyright 2017 by Robert Stengel. All rights reserved. For educational use only.

Save this PDF as:

Size: px
Start display at page:

Download "Control Systems! Copyright 2017 by Robert Stengel. All rights reserved. For educational use only."

## Transcription

1 Control Systems Robert Stengel Robotics and Intelligent Systems MAE 345, Princeton University, 2017 Analog vs. digital systems Continuous- and Discretetime Dynamic Models Frequency Response Transfer Functions Bode Plots Root Locus Proportional-Integral- Derivative (PID) Control Copyright 2017 by Robert Stengel. All rights reserved. For educational use only. 1 Analog vs. Digital Signals Signal: A physical indicator that conveys information, e.g., position, voltage, temperature, or pressure Analog signal: A signal that is continuous in time, with infinite precision and infinitesimal time spacing between indication points Discrete-time signal: A signal with infinite precision and discrete spacing between indication points Digital signal: A signal with finite precision and discrete spacing between indication points 2

2 Analog vs. Digital Systems System: Assemblage of parts with structure, connectivity, and behavior that responds to input signals and produces output signals Analog system: A system that operates continuously, with infinite precision and infinitesimal time spacing between signaling points Discrete-time system: A system that operates continuously, with infinite precision and discrete spacing between signaling points Digital system: A system that operates continuously, with finite precision and discrete spacing between signaling points 3 Frequency Response of a Continuous- Time Dynamic System Sinusoidal command input (i.e., Desired rotational angle) Long-term output of the dynamic system: Sinusoid with same frequency as input Output/input amplitude ratio dependent on input frequency Output/input phase shift dependent on input frequency Bandwidth: Input frequency below which amplitude and phase angle variations are negligible < n = n > n Output Angle, deg Output Rate, deg/s Input Angle, deg : Input frequency, rad/s n : Natural frequency of system, rad/s Hz Hz : Hertz, frequency, cycles per sec = ( rad / s) 2 4

3 Amplitude Ratio and Phase Angle Amplitude Ratio (AR) = y Output Peak y Input Peak Phase Angle = 360 t Input Peak t Output Peak Period of Input, deg 5 Amplitude-Ratio and Phase-Angle Frequency Response of a 2 nd -Order System Output Angle AR Straight-line asymptotes at high and low frequency 6 db Bandwidth, Hz Output Rate AR 20log AR and vs. log Straight-line asymptotes at high and low frequency Logarithm of transfer function Amplitude ratio and phase angle are orthogonal Decibels (db): 20 log AR Decade: power of 10 6

4 Asymptotes of 2 nd -Order Frequency Response from Bode Plot Natural Frequency 20 log 10 ( n ) has slope of 20n db/decade on Bode plot 20log 10 H CLi (j ) = 20log 10 AR i ( ) e j i ( ) = 20 log 10 AR i + log 10 e j i ( ) = 20log 10 AR i + j i ( ) 20log 10 e 7 F x(t) Fourier Transform of a Scalar Variable [ ] = x( j ) = x(t)e jt dt = frequency, rad / s j i 1 x(t) x( j) = a() + jb() x(t) :real variable x( j ) :complex variable = a( ) + jb( ) j ( ) = A( )e A :amplitude :phase angle 8

5 Laplace Transforms of Scalar Variables Laplace transform of a scalar variable is a complex number s is the Laplace operator, a complex variable [ ] = x = x(t)e st dt L x(t), s = + j 0 x(t) :real variable x : complex variable = a( ) + jb( ) j ( ) = A( )e Multiplication by a constant L a x(t) Sum of Laplace transforms [ ] = x 1 + x 2 L x 1 (t) + x 2 (t) [ ] = a x 9 Laplace Transforms of Vectors and Matrices Laplace transform of a vector variable Laplace transform of a matrix variable L[ x(t) ] = x = x 1 x 2... L[ A(t) ] = A = a 11 a a 21 a Laplace transform of a derivative w.r.t. time L dx(t) dt = sx x(0) Laplace transform of an integral over time L x( )d ( = x s

6 Laplace Transforms of the System Equations Time-Domain System Equations x(t) = Fx(t) + Gu(t) y(t) = H x x(t) + H u u(t) Dynamic Equation Output Equation Laplace Transforms of System Equations sx x(0) = Fx + Gu y = H x x + H u u Dynamic Equation Output Equation 11 Laplace Transform of State Response to Initial Condition and Control Adj si F Rearrange sx Fx = x(0) + Gu [ si F]x = x(0) + Gu [ si F] 1 = [ ] 1 [ x(0) + Gu ] x = si F Matrix inverse si F Adj si F (n n) : Adjoint matrix (n n) : Transpose of matrix of cofactors si F = det( si F) : Determinant ( 11) 12

7 Laplace Transform of System Output Transfer function relates control input to system output y = H x { x } + H u u { 1 [ Gu + x(0) ]} + H u u = H x si F Let: x(0) = 0 [initial condition, not control] H u = 0 [control input not contained in the output measurement] y = H x ( si F) 1 G u H u 13 Transfer Function Matrix, H [from u to y] Dimension of y = r 1 H H x ( si F) 1 G G 1 Adj si F = H x si F Dimension of the Matrix ( r n) n n 11 Dimension of x = n 1 Dimension of u = m 1 ( n m) = ( r m) 14

8 Characteristic Polynomial of a Dynamic System Characteristic matrix ( si F) = ( s f 11 ) f f 1n f 21 s f f 2n f n1 f n2... s f nn (n ( n) Characteristic polynomial, si F = det( si F) = s n + a n1 s n a 1 s + a 0 15 Eigenvalues (or Roots) of the Dynamic System Characteristic equation = s n + a n1 s n a 1 s + a 0 = 0 ( s 2 )... = s 1 = 0 s n i are eigenvalues of F or roots of the characteristic polynomial, s is the denominator of the transfer function matrix s 16

9 Numerator of the Matrix Inverse [ si F] 1 = Matrix Inverse si F Adj si F (n n) Adjoint matrix is the transpose of the matrix of cofactors * Adj( si F) = C T (n n) ( si F) = 2 x 2 example ( s f 11 ) f 12 f 21 s f 22 Adj( si F) = = ( s f 22 ) f 21 f 12 s f 11 ( s f 22 ) f 12 f 21 s f 11 T * Cofactors = signed minor determinants of the matrix 17 Example of System Transformation DC Motor Equations x 1 (t) x 2 (t) = 0 1 x 1 (t) 0 0 x 2 (t) + 0 u(t) 1/ J y 1 (t) y 2 (t) = 1 0 x 1 (t) 0 1 x 2 (t) + 0 x 1 (t) u(t) = 0 x 2 (t) Laplace Transforms of Motor Equations Dynamic Equation Output Equation sx 1 x 1 (0) sx 2 x 2 (0) y 1 y 2 = = x 1 x 2 x 1 x / J u u = x 1 x 2 Dynamic Equation Output Equation 18

10 2 x 2 Eigenvalue Example Characteristic matrix ( si F) = ( s f 11 ) f 12 f 21 s f 22 Determinant of characteristic matrix si F = ( s f 11 ) f 12 f 21 s f 22 ( s f 22 ) f 12 f 21 = s f 11 = s 2 ( f 11 + f 22 )s + f 11 f 22 f 12 f Factors of the 2 nd - Degree Characteristic Equation x = s 2 f 12 + f 21 s s + f 11 f 22 f 12 f 21 = 0 x x = cos 1 = ( s 1 )( s 2 ) = 0 Solutions of the equation are eigenvalues of the system, either 2 real roots, or Complex-conjugate pair x s Plane 1 = 1 + j 1, 2 = 1 j 1 1 = 1, 2 = 2 ( s) = s n s + n 20

11 Eigenvalues Determine the Stability of the LTI System Positive real part represents instability Envelope of time response converges or diverges x x x 1 (t) = e nt cos n 1 2 t + 1 ( ) x 1(0) s Plane Same criterion for real roots x 1 (t) = e nt x 1 (0) 21 Open-Loop Transfer Function Matrix for DC Motor Transfer Function Matrix H OL = H x [ si F OL ] 1 G (r m) where F OL = H x = ; [ si F OL ] 1 = Adj si F OL = si F OL [ si F OL ] = = I 2 ; G = Matrix Inverse s 1 0 s J = s 1 0 s = s 2 s 1 0 s 1/ s 1/ s 2 0 1/ s 22

12 Simplify H OL = H x [ si F OL ] 1 Adj si F G = H OL x si F OL = s 1 0 s 0 s 2 1/ J = G Dimension = 2 x 1 1/ Js 2 1/ Js y 1 u y 2 u H 1 H 2 ( s) ( s) = 1 Js 2 1 Js Angle = Double integral of input torque, u Angular rate = Integral of input torque, u 23 Closed-Loop System x 1 (t) x 2 (t) Closed-loop dynamic equation = 0 1 x 1 (t) c 1 / J c 2 / J x 2 (t) 0 1 2)( n ( n 2 x 1 (t) x 2 (t) + 0 c 1 / J + 0 ( n 2 y C y C Natural Frequency: n = c 1 J Damping Ratio: = ( c 2 J ) 2 n = ( c 2 J ) 2 c 1 J 24

13 Closed-Loop Transfer Function Matrix H CL = H x [ si F CL ] 1 Adj si F G = H CL x si F CL = s + c 2 J c 1 J s s 2 + c 2 J s + c ( 1 + ) * J, - 0 c 1 / J Scalar transfer functions differ only in numerators G H 1 H 2 ( s) ( s) = y 1 y C y 2 y C c 1 J s( c 1 J ) = s 2 + c 2 J s + c 1 * ( ) J +, n 1 n 2 ( s) ( s) - s n 1 n 2 ( s) ( s) d s 25 H 1 ( j ) = Angle Frequency Response = Substitute s = j in transfer functions Output Angle Commanded Angle = y ( j ) ( c 2 J )( j ) + c 1 J y C ( j ) = c 1 J j 2 n n ( j ) + a 2 1( ) + jb 1 ( ) AR 1 ( ) e j 1( ) n 26

14 Angular Rate Frequency Response H 2 ( j ) = Output Rate Commanded Angle = y 2( j ) 2 ( j ) = n n j c 1 J 2 + ( c 2 J )( j ) + c 1 J y C ( j ) = j j + n 2 AR 2( ) e j 2 ( ) 27 y 1 y C y 2 y C Root (Eigenvalue) Locus c 1 J ( c 1 J )s = s 2 + c 2 J s + c = 1 * ( ) J +, 2 - n - 2 n s s n s + - n = ( s n s +1) Variation of roots as scalar gain, k, goes from 0 to With nominal gains, c 1 and c 2, 1 s n = c 1 J = 1 rad / s, Natural frequency = c 2 J 2 n = 0.707, Damping ratio Root Locus of DC Motor Angle Control F = [0 1; ]; G = [0;1]; Hx1 = [1 0]; Angle Output Hx2 = [0 1]; Angular Rate Output Sys1 Sys2 = ss(f,g,hx1,0); = ss(f,g,hx2,0); rlocus(sys1), grid figure rlocus(sys2), grid 28

15 Root (Eigenvalue) Locus Increase Angle Feedback Gain, Sys1 Increase Rate Feedback Gain, Sys2 Initial Root Initial Root Zero Initial Root Initial Root See for rules of root locus construction 29 Classical Control System Design Criteria Step Response Frequency Response Eigenvalue Location 30

16 Single-Axis Angular Control of Non-Spinning Spacecraft Pitching motion (about the y axis) is to be controlled (t) q(t) = (t) q(t) / I yy M (t) y (t) q(t) = Pitch Angle Pitch Rate Identical to the DC Motor control problem 31 Proportional-Integral- Derivative (PID) Controller Control Error Control Command Control Law Transfer Function (w/common denominator) e = C e u = c P e + c I + c D se s u e = c Ps + c I + c D s 2 s Proportional term weights control error directly Integrator compensates for persistent (bias) disturbance Differentiator produces rate term for damping Ziegler-Nichols PID Tuning Method Nichols_method 32

17 Open-Loop Angle Transfer Function g A = Actuator gain H OL ( s) = e = c + c s + c I P D s2 s g A I yy s 2 33 Closed-Loop Spacecraft Control Transfer Function w/pid Control = H OL ( s)e = H OL s ( s) +1 = H OL ( s) c H OL H CL ( s) = = c = H OL s 1+H OL s [ ] c e 1+ e 34

18 Closed-Loop Frequency Response w/pid Control c I + c P s + c D s 2 g A * ) s I yy s 2, H CL ( s) = () +, c = D s 2 + c P s + c I 1+ c + c s + c I P D s2 g A * I yy s 3 g A + c D s 2 + c P s + c I ) s I yy s 2, () +, Let s = j. As 0 As H CL ( j ) = H CL ( j ) = ( j ) c ( j ) c I c I = 1 Steady-state output = desired steady-state input ( j ) c ( j ) c D 2 ji yy 3 g A = c D ji yy g A = jc D I yy g A AR c D I yy g A; 90 deg High-frequency response rolls off and lags input 35 Digital Control Stengel, Optimal Control and Estimation, 1994, pp , E-Reserve 36

19 From Continuous- to Discrete- Time Systems Continuous-time systems are described by differential equations, e.g., x(t) = Fx(t) + Gu(t) Discrete-time systems are described by difference equations, e.g., x(t k +1 ) = x(t k ) + u(t k ) Discrete-time systems that are meant to describe continuous-time systems are called sampled-data systems = fcn( F); = fcn( F,G) 37 Integration of Linear, Time-Invariant (LTI) Systems x(t) = Fx(t) + Gu(t) t x(t) = x(0) + xd = x(0) + [ Fx + Gu]d 0 Homogeneous solution (no input), u = 0 x(t) = x(0) + [ Fx]d = e F(t0) x(0) = ( t,0)x(0) t 0 ( t,0) = e F(t0) = e Ft = State transition matrix from 0 to t t 0 38

20 Equivalence of 1 st -Order Continuousand Discrete-Time Systems Continuous- and discrete-time values are identical at time instants, kt x(t) = ax(t), x(0) given t x(t) = x(0) + x(t)dt = e at x(0) (t)x(0) 0 Choose small time interval, t Propagate state to time, t + nt 39 Equivalence of 1 st -Order Continuous- and Discrete-Time Systems x 0 = x(0) x 1 = x(t) = e at x 0 = (t)x 0 x 2 = x(2t) = e 2at x 0 = e at x 1 = (t)x 1... x k = x(kt) = e kat x 0 = e at x k1 = (t)x k1... x n = x(nt) = x(t) = e nat x 0 = (t)x n1 40

21 Input Response of LTI System Inhomogeneous (forced) solution t ( x(t) = e Ft x(0) + e F(t ) Gu d = ( t 0)x(0) + t 0 0 t t ) 0 Gu( ) ( d 41 Sampled-Data System (Stepwise application of prior results) t k+1 t k Assume (t k t k+1 ) = t = constant = e F(t k+1t k ) = e F(t ) = t Assume u = constant from t k to t k+1 x(t k+1 ) = (t)x(t k ) + t t t Evaluate the integral ) G ( d u k x(t k+1 ) = (t)x(t k ) + (t) I 1 (t) F 1 Gu t k = (t)x(t k ) + ((t)u( t k ) 0 42

22 Digital Feedback Control System x(t k+1 ) = (t)x(t k ) + (t)u( t k ) Linear feedback control law with command input u( t k ) = u c ( t k ) C k x( t k ) = t x t k+1 Closed-loop system x(t k ) + t u c ( t k ) C k x( t k ) ( = (t) (t)c k ( x(t k ) + (t)u c t k 43 Digital Flight Control Historical notes: NASA Fly-By-Wire F-8C (Apollo GNC system, 1972) 1 st conventional aircraft with digital flight control NASA Armstrong (Dryden) Research Center Princetons Variable-Response Research Aircraft (Z-80 microprocessor, Spring, 1978) 2 nd conventional aircraft with digital flight control Princeton University s Flight Research Laboratory 44

23 Next Time: Sensors and Actuators 45 SupplementaL Material 46

24 Princetons Variable- Response Research Aircraft 47 Eigenvalues, Damping Ratio, and Natural Frequency Eigenvalues, Damping Ratio, and Natural Frequency of Angle Control F1 = [0 1;-1 0]; G1 = [0;1]; F1a = [0 1; ]; F1b = [0 1; ]; Hx = [1 0;0 1]; Eigenvalues 0 + 1i 0-1i Damping Ratio, Natural Frequency 1, 2 n 0 1 Sys1 Sys2 Sys3 eig(f1) eig(f1a) eig(f1b) damp(sys1) damp(sys2) = ss(f1,g1,hx,0); = ss(f1a,g1,hx,0); = ss(f1b,g1,hx,0); i i Overdamped 48

25 Bode Plots of Proportional, Integral, and Derivative Compensation H( j) =1 H( j) =10 H( j) =100 H( j) = 1 j H( j) = 10 j H( j) = j H( j) =10 j Amplitude ratio, AR(dB) = Slope of AR(dB) = Slope of AR(dB) = constant 20 db/decade +20 db/decade Phase Angle, ( deg) = 0 Phase Angle, = 90 deg Phase Angle, = +90 deg 49 Equilibrium Response of Linear Discrete-Time Model Dynamic model x k +1 = x k + u k + w k At equilibrium x k +1 = x k = x* = constant Equilibrium response ( I )x* = u*+w * 1 ( u*+w *) x* = I (.)* = constant 50

26 Discrete-Time and Sampled-Data Models Linear, time-invariant, discrete-time model x k +1 = x k + u k + w k Discrete-time model is a sampled-data model if it represents a continuous-time system 51 Sampled-Data (Digital) Feedback Control Digital-to-Analog Conversion Analog-to-Digital Conversion 52

27 Sampling Effect on Continuous Signal Two waves, different frequencies, indistinguishable in periodically sampled data Frequencies above the sampling frequency aliased to appear at lower frequency (frequency folding) Solutions: Either Sampling at higher rate or Analog low-pass filtering Folding Frequency Aliased Spectrum Actual Spectrum 53 Quantization and Delay Effects Continuous signal sampled with finite precision (quantized) Solution: More bits in A/D and D/A conversion Effective delay of sampled signal Described as phase shift or lag Solution: Higher sampling rate or Lead compensation Analog-to-Digital Converter r h illustrates zero-order hold effect 54

28 Time Delay Effect Control command delayed by sec Laplace transform of pure time delay L u( t ) = e s u( s) Delay introduces frequency-dependent phase lag with no change in amplitude AR e j = 1 = e j Phase lag reduces closed-loop stability 55 Effect of Closed-Loop Time Delay on Step Response With no delay n = 1 rad / s, = With varying delays 56

29 Discrete-Time Frequency Response can be Evaluated Using the z Transform x k+1 = x k + u k z transform is the Laplace transform of a periodic sequence System z transform is zx(z) x(0) = x(z) + u(z) which leads to [ ] x(z) = ( zi ) 1 x(0) + u(z) Further details are beyond the present scope 57 Factors That Complicate Precise Control 58

30 Error, Uncertainty, and Incompleteness May Have Significant Effects Scale factors and biases: y = kx + b Disturbances: Constant/variable forces Measurement error: Noise, bias Modeling error Parameters: e.g., Inertia Structure: e.g., Higher-order modes 59 Saturation (limiting) Displacement limit, maximum force or rate Threshold Dead zone, slop Preload Breakout force Friction Sliding surface resistance Rectifier Hard constraint, absolute value (double) Hysteresis Backlash, slop Higher-degree terms Cubic spring, separation Nonlinearity Complicates Response 60

31 Second-Order Example: Airplane Roll Motion p ( ( = Roll rate, rad/s ( Roll angle, rad ( ) A = Aileron deflection, rad L p : Rolling motion of an airplane, continuous-time p ( ( = L p 0 p ( ( 1 0 ( ( + L ) A 0 roll damping coefficient L A : aileron control effect ( ) A ( Rolling motion of an airplane, discrete-time p k k ( ( = e L pt 0 ( e LpT )1) 1 L p ( ( ( ( p k)1 k)1 ( ( + ~ L T * A 0 ( * A k)1 ( T = sampling interval, s 61

### Translational and Rotational Dynamics!

Translational and Rotational Dynamics Robert Stengel Robotics and Intelligent Systems MAE 345, Princeton University, 217 Copyright 217 by Robert Stengel. All rights reserved. For educational use only.

More information

### First-Order Low-Pass Filter!

Filters, Cost Functions, and Controller Structures! Robert Stengel! Optimal Control and Estimation MAE 546! Princeton University, 217!! Dynamic systems as low-pass filters!! Frequency response of dynamic

More information

### Stochastic Optimal Control!

Stochastic Control! Robert Stengel! Robotics and Intelligent Systems, MAE 345, Princeton University, 2015 Learning Objectives Overview of the Linear-Quadratic-Gaussian (LQG) Regulator Introduction to Stochastic

More information

### Linearized Equations of Motion!

Linearized Equations of Motion Robert Stengel, Aircraft Flight Dynamics MAE 331, 216 Learning Objectives Develop linear equations to describe small perturbational motions Apply to aircraft dynamic equations

More information

### SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015

FACULTY OF ENGINEERING AND SCIENCE SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015 Lecturer: Michael Ruderman Problem 1: Frequency-domain analysis and control design (15 pt) Given is a

More information

### Dr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review

Week Date Content Notes 1 6 Mar Introduction 2 13 Mar Frequency Domain Modelling 3 20 Mar Transient Performance and the s-plane 4 27 Mar Block Diagrams Assign 1 Due 5 3 Apr Feedback System Characteristics

More information

### Feedback Control of Linear SISO systems. Process Dynamics and Control

Feedback Control of Linear SISO systems Process Dynamics and Control 1 Open-Loop Process The study of dynamics was limited to open-loop systems Observe process behavior as a result of specific input signals

More information

### Learn2Control Laboratory

Learn2Control Laboratory Version 3.2 Summer Term 2014 1 This Script is for use in the scope of the Process Control lab. It is in no way claimed to be in any scientific way complete or unique. Errors should

More information

### Exam. 135 minutes + 15 minutes reading time

Exam January 23, 27 Control Systems I (5-59-L) Prof. Emilio Frazzoli Exam Exam Duration: 35 minutes + 5 minutes reading time Number of Problems: 45 Number of Points: 53 Permitted aids: Important: 4 pages

More information

### Optimal Control and Estimation MAE 546, Princeton University Robert Stengel, Preliminaries!

Optimal Control and Estimation MAE 546, Princeton University Robert Stengel, 2017 Copyright 2017 by Robert Stengel. All rights reserved. For educational use only. http://www.princeton.edu/~stengel/mae546.html

More information

### Dynamic Effects of Feedback Control!

Dynamic Effecs of Feedback Conrol! Rober Sengel! Roboics and Inelligen Sysems MAE 345, Princeon Universiy, 2017 Inner, Middle, and Ouer Feedback Conrol Loops Sep Response of Linear, Time- Invarian (LTI)

More information

### INTRODUCTION TO DIGITAL CONTROL

ECE4540/5540: Digital Control Systems INTRODUCTION TO DIGITAL CONTROL.: Introduction In ECE450/ECE550 Feedback Control Systems, welearnedhow to make an analog controller D(s) to control a linear-time-invariant

More information

### Recitation 11: Time delays

Recitation : Time delays Emilio Frazzoli Laboratory for Information and Decision Systems Massachusetts Institute of Technology November, 00. Introduction and motivation. Delays are incurred when the controller

More information

### APPLICATIONS FOR ROBOTICS

Version: 1 CONTROL APPLICATIONS FOR ROBOTICS TEX d: Feb. 17, 214 PREVIEW We show that the transfer function and conditions of stability for linear systems can be studied using Laplace transforms. Table

More information

### Autonomous Mobile Robot Design

Autonomous Mobile Robot Design Topic: Guidance and Control Introduction and PID Loops Dr. Kostas Alexis (CSE) Autonomous Robot Challenges How do I control where to go? Autonomous Mobile Robot Design Topic:

More information

### Topic # Feedback Control. State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback

Topic #17 16.31 Feedback Control State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback Back to reality Copyright 21 by Jonathan How. All Rights reserved 1 Fall

More information

### Today (10/23/01) Today. Reading Assignment: 6.3. Gain/phase margin lead/lag compensator Ref. 6.4, 6.7, 6.10

Today Today (10/23/01) Gain/phase margin lead/lag compensator Ref. 6.4, 6.7, 6.10 Reading Assignment: 6.3 Last Time In the last lecture, we discussed control design through shaping of the loop gain GK:

More information

### Automatique. A. Hably 1. Commande d un robot mobile. Automatique. A.Hably. Digital implementation

A. Hably 1 1 Gipsa-lab, Grenoble-INP ahmad.hably@grenoble-inp.fr Commande d un robot mobile (Gipsa-lab (DA)) ASI 1 / 25 Outline 1 2 (Gipsa-lab (DA)) ASI 2 / 25 of controllers Signals must be sampled and

More information

### Lecture 6 Classical Control Overview IV. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Lecture 6 Classical Control Overview IV Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lead Lag Compensator Design Dr. Radhakant Padhi Asst.

More information

### Nonlinear State Estimation! Extended Kalman Filters!

Nonlinear State Estimation! Extended Kalman Filters! Robert Stengel! Optimal Control and Estimation, MAE 546! Princeton University, 2017!! Deformation of the probability distribution!! Neighboring-optimal

More information

### Aircraft Flight Dynamics!

Aircraft Flight Dynamics Robert Stengel MAE 331, Princeton University, 2016 Course Overview Introduction to Flight Dynamics Math Preliminaries Copyright 2016 by Robert Stengel. All rights reserved. For

More information

### Richiami di Controlli Automatici

Richiami di Controlli Automatici Gianmaria De Tommasi 1 1 Università degli Studi di Napoli Federico II detommas@unina.it Ottobre 2012 Corsi AnsaldoBreda G. De Tommasi (UNINA) Richiami di Controlli Automatici

More information

### Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 3.. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -

More information

### Example: Modeling DC Motor Position Physical Setup System Equations Design Requirements MATLAB Representation and Open-Loop Response

Page 1 of 5 Example: Modeling DC Motor Position Physical Setup System Equations Design Requirements MATLAB Representation and Open-Loop Response Physical Setup A common actuator in control systems is the

More information

### DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2010

[E2.5] IMPERIAL COLLEGE LONDON DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2010 EEE/ISE PART II MEng. BEng and ACGI SIGNALS AND LINEAR SYSTEMS Time allowed: 2:00 hours There are FOUR

More information

### Lecture 1: Feedback Control Loop

Lecture : Feedback Control Loop Loop Transfer function The standard feedback control system structure is depicted in Figure. This represend(t) n(t) r(t) e(t) u(t) v(t) η(t) y(t) F (s) C(s) P (s) Figure

More information

### Time Varying Circuit Analysis

MAS.836 Sensor Systems for Interactive Environments th Distributed: Tuesday February 16, 2010 Due: Tuesday February 23, 2010 Problem Set # 2 Time Varying Circuit Analysis The purpose of this problem set

More information

### BIBO STABILITY AND ASYMPTOTIC STABILITY

BIBO STABILITY AND ASYMPTOTIC STABILITY FRANCESCO NORI Abstract. In this report with discuss the concepts of bounded-input boundedoutput stability (BIBO) and of Lyapunov stability. Examples are given to

More information

### Review of Linear Time-Invariant Network Analysis

D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x

More information

### The Relation Between the 3-D Bode Diagram and the Root Locus. Insights into the connection between these classical methods. By Panagiotis Tsiotras

F E A T U R E The Relation Between the -D Bode Diagram and the Root Locus Insights into the connection between these classical methods Bode diagrams and root locus plots have been the cornerstone of control

More information

### Minimization of Static! Cost Functions!

Minimization of Static Cost Functions Robert Stengel Optimal Control and Estimation, MAE 546, Princeton University, 2017 J = Static cost function with constant control parameter vector, u Conditions for

More information

### Control System Design

ELEC4410 Control System Design Lecture 19: Feedback from Estimated States and Discrete-Time Control Design Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science

More information

### KINGS COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING

KINGS COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING QUESTION BANK SUB.NAME : CONTROL SYSTEMS BRANCH : ECE YEAR : II SEMESTER: IV 1. What is control system? 2. Define open

More information

### Rotary Motion Servo Plant: SRV02. Rotary Experiment #01: Modeling. SRV02 Modeling using QuaRC. Student Manual

Rotary Motion Servo Plant: SRV02 Rotary Experiment #01: Modeling SRV02 Modeling using QuaRC Student Manual SRV02 Modeling Laboratory Student Manual Table of Contents 1. INTRODUCTION...1 2. PREREQUISITES...1

More information

### Lab Experiment 2: Performance of First order and second order systems

Lab Experiment 2: Performance of First order and second order systems Objective: The objective of this exercise will be to study the performance characteristics of first and second order systems using

More information

### Laplace Transform Analysis of Signals and Systems

Laplace Transform Analysis of Signals and Systems Transfer Functions Transfer functions of CT systems can be found from analysis of Differential Equations Block Diagrams Circuit Diagrams 5/10/04 M. J.

More information

### Steady State Frequency Response Using Bode Plots

School of Engineering Department of Electrical and Computer Engineering 332:224 Principles of Electrical Engineering II Laboratory Experiment 3 Steady State Frequency Response Using Bode Plots 1 Introduction

More information

### KINGS COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING

KINGS COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING QUESTION BANK SUB.NAME : CONTROL SYSTEMS STAFF NAME: Mr. P.NARASIMMAN BRANCH : ECE Mr.K.R.VENKATESAN YEAR : II SEMESTER

More information

### CONTROL OF DIGITAL SYSTEMS

AUTOMATIC CONTROL AND SYSTEM THEORY CONTROL OF DIGITAL SYSTEMS Gianluca Palli Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI) Università di Bologna Email: gianluca.palli@unibo.it

More information

### CHAPTER 1 Basic Concepts of Control System. CHAPTER 6 Hydraulic Control System

CHAPTER 1 Basic Concepts of Control System 1. What is open loop control systems and closed loop control systems? Compare open loop control system with closed loop control system. Write down major advantages

More information

### FEEDBACK CONTROL SYSTEMS

FEEDBAC CONTROL SYSTEMS. Control System Design. Open and Closed-Loop Control Systems 3. Why Closed-Loop Control? 4. Case Study --- Speed Control of a DC Motor 5. Steady-State Errors in Unity Feedback Control

More information

### Frequency Response Analysis

Frequency Response Analysis Consider let the input be in the form Assume that the system is stable and the steady state response of the system to a sinusoidal inputdoes not depend on the initial conditions

More information

### EECS C128/ ME C134 Final Thu. May 14, pm. Closed book. One page, 2 sides of formula sheets. No calculators.

Name: SID: EECS C28/ ME C34 Final Thu. May 4, 25 5-8 pm Closed book. One page, 2 sides of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 4 2 4 3 6 4 8 5 3

More information

### INSTRUMENTAL ENGINEERING

INSTRUMENTAL ENGINEERING Subject Code: IN Course Structure Sections/Units Section A Unit 1 Unit 2 Unit 3 Unit 4 Unit 5 Unit 6 Section B Section C Section D Section E Section F Section G Section H Section

More information

### Dynamic System Response. Dynamic System Response K. Craig 1

Dynamic System Response Dynamic System Response K. Craig 1 Dynamic System Response LTI Behavior vs. Non-LTI Behavior Solution of Linear, Constant-Coefficient, Ordinary Differential Equations Classical

More information

### Dr. Ian R. Manchester

Dr Ian R. Manchester Week Content Notes 1 Introduction 2 Frequency Domain Modelling 3 Transient Performance and the s-plane 4 Block Diagrams 5 Feedback System Characteristics Assign 1 Due 6 Root Locus

More information

### Dynamic Optimal Control!

Dynamic Optimal Control! Robert Stengel! Robotics and Intelligent Systems MAE 345, Princeton University, 2017 Learning Objectives Examples of cost functions Necessary conditions for optimality Calculation

More information

### Control Systems I. Lecture 2: Modeling. Suggested Readings: Åström & Murray Ch. 2-3, Guzzella Ch Emilio Frazzoli

Control Systems I Lecture 2: Modeling Suggested Readings: Åström & Murray Ch. 2-3, Guzzella Ch. 2-3 Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich September 29, 2017 E. Frazzoli

More information

### EECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators.

Name: SID: EECS C28/ ME C34 Final Wed. Dec. 5, 2 8- am Closed book. Two pages of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 2 2 6 3 4 4 5 6 6 7 8 2 Total

More information

### Physics 401. Fall 2017 Eugene V. Colla. 10/9/2017 Physics 401 1

Physics 41. Fall 17 Eugene V. Colla 1/9/17 Physics 41 1 Electrical RLC circuits Torsional Oscillator Damping Data Analysis 1/9/17 Physics 41 V R +V L +V C =V(t) If V(t)= R d d q(t) q(t) dt dt C C L q(t)

More information

### DC Motor Position: System Modeling

1 of 7 01/03/2014 22:07 Tips Effects TIPS ABOUT BASICS INDEX NEXT INTRODUCTION CRUISE CONTROL MOTOR SPEED MOTOR POSITION SUSPENSION INVERTED PENDULUM SYSTEM MODELING ANALYSIS DC Motor Position: System

More information

### Power System Control

Power System Control Basic Control Engineering Prof. Wonhee Kim School of Energy Systems Engineering, Chung-Ang University 2 Contents Why feedback? System Modeling in Frequency Domain System Modeling in

More information

### Chapter 6: The Laplace Transform. Chih-Wei Liu

Chapter 6: The Laplace Transform Chih-Wei Liu Outline Introduction The Laplace Transform The Unilateral Laplace Transform Properties of the Unilateral Laplace Transform Inversion of the Unilateral Laplace

More information

### Modeling and Control Overview

Modeling and Control Overview D R. T A R E K A. T U T U N J I A D V A N C E D C O N T R O L S Y S T E M S M E C H A T R O N I C S E N G I N E E R I N G D E P A R T M E N T P H I L A D E L P H I A U N I

More information

### Methods for analysis and control of. Lecture 6: Introduction to digital control

Methods for analysis and of Lecture 6: to digital O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.inpg.fr www.lag.ensieg.inpg.fr/sename 6th May 2009 Outline Some interesting books:

More information

### Pitch Rate CAS Design Project

Pitch Rate CAS Design Project Washington University in St. Louis MAE 433 Control Systems Bob Rowe 4.4.7 Design Project Part 2 This is the second part of an ongoing project to design a control and stability

More information

### Additional Closed-Loop Frequency Response Material (Second edition, Chapter 14)

Appendix J Additional Closed-Loop Frequency Response Material (Second edition, Chapter 4) APPENDIX CONTENTS J. Closed-Loop Behavior J.2 Bode Stability Criterion J.3 Nyquist Stability Criterion J.4 Gain

More information

### Comparison of four state observer design algorithms for MIMO system

Archives of Control Sciences Volume 23(LIX), 2013 No. 2, pages 131 144 Comparison of four state observer design algorithms for MIMO system VINODH KUMAR. E, JOVITHA JEROME and S. AYYAPPAN A state observer

More information

### Dynamic circuits: Frequency domain analysis

Electronic Circuits 1 Dynamic circuits: Contents Free oscillation and natural frequency Transfer functions Frequency response Bode plots 1 System behaviour: overview 2 System behaviour : review solution

More information

### Transient response via gain adjustment. Consider a unity feedback system, where G(s) = 2. The closed loop transfer function is. s 2 + 2ζωs + ω 2 n

Design via frequency response Transient response via gain adjustment Consider a unity feedback system, where G(s) = ωn 2. The closed loop transfer function is s(s+2ζω n ) T(s) = ω 2 n s 2 + 2ζωs + ω 2

More information

### Solutions to Skill-Assessment Exercises

Solutions to Skill-Assessment Exercises To Accompany Control Systems Engineering 4 th Edition By Norman S. Nise John Wiley & Sons Copyright 2004 by John Wiley & Sons, Inc. All rights reserved. No part

More information

### State Space Representation

ME Homework #6 State Space Representation Last Updated September 6 6. From the homework problems on the following pages 5. 5. 5.6 5.7. 5.6 Chapter 5 Homework Problems 5.6. Simulation of Linear and Nonlinear

More information

### Homework 7 - Solutions

Homework 7 - Solutions Note: This homework is worth a total of 48 points. 1. Compensators (9 points) For a unity feedback system given below, with G(s) = K s(s + 5)(s + 11) do the following: (c) Find the

More information

### (b) A unity feedback system is characterized by the transfer function. Design a suitable compensator to meet the following specifications:

1. (a) The open loop transfer function of a unity feedback control system is given by G(S) = K/S(1+0.1S)(1+S) (i) Determine the value of K so that the resonance peak M r of the system is equal to 1.4.

More information

### Feedback Basics. David M. Auslander Mechanical Engineering University of California at Berkeley. copyright 1998, D.M. Auslander

Feedback Basics David M. Auslander Mechanical Engineering University of California at Berkeley copyright 1998, D.M. Auslander 1 I. Feedback Control Context 2 What is Feedback Control? Measure desired behavior

More information

### PID controllers, part I

Faculty of Mechanical and Power Engineering Dr inŝ. JANUSZ LICHOTA CONTROL SYSTEMS PID controllers, part I Wrocław 2007 CONTENTS Controller s classification PID controller what is it? Typical controller

More information

### ROOT LOCUS. Consider the system. Root locus presents the poles of the closed-loop system when the gain K changes from 0 to. H(s) H ( s) = ( s)

C1 ROOT LOCUS Consider the system R(s) E(s) C(s) + K G(s) - H(s) C(s) R(s) = K G(s) 1 + K G(s) H(s) Root locus presents the poles of the closed-loop system when the gain K changes from 0 to 1+ K G ( s)

More information

### Automatic Control 2. Loop shaping. Prof. Alberto Bemporad. University of Trento. Academic year

Automatic Control 2 Loop shaping Prof. Alberto Bemporad University of Trento Academic year 21-211 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 21-211 1 / 39 Feedback

More information

### Singular Value Analysis of Linear- Quadratic Systems!

Singular Value Analyi of Linear- Quadratic Sytem! Robert Stengel! Optimal Control and Etimation MAE 546! Princeton Univerity, 2017!! Multivariable Nyquit Stability Criterion!! Matrix Norm and Singular

More information

### Distributed Real-Time Control Systems

Distributed Real-Time Control Systems Chapter 9 Discrete PID Control 1 Computer Control 2 Approximation of Continuous Time Controllers Design Strategy: Design a continuous time controller C c (s) and then

More information

### Dynamic Response. Assoc. Prof. Enver Tatlicioglu. Department of Electrical & Electronics Engineering Izmir Institute of Technology.

Dynamic Response Assoc. Prof. Enver Tatlicioglu Department of Electrical & Electronics Engineering Izmir Institute of Technology Chapter 3 Assoc. Prof. Enver Tatlicioglu (EEE@IYTE) EE362 Feedback Control

More information

### VALLIAMMAI ENGINEERING COLLEGE

VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur 603 203 DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING QUESTION BANK V SEMESTER IC650 CONTROL SYSTEMS Regulation 203 Academic Year 207 8 Prepared

More information

### Analysis and Design of Control Systems in the Time Domain

Chapter 6 Analysis and Design of Control Systems in the Time Domain 6. Concepts of feedback control Given a system, we can classify it as an open loop or a closed loop depends on the usage of the feedback.

More information

### Professor Fearing EE C128 / ME C134 Problem Set 7 Solution Fall 2010 Jansen Sheng and Wenjie Chen, UC Berkeley

Professor Fearing EE C8 / ME C34 Problem Set 7 Solution Fall Jansen Sheng and Wenjie Chen, UC Berkeley. 35 pts Lag compensation. For open loop plant Gs ss+5s+8 a Find compensator gain Ds k such that the

More information

### Advanced Digital Controls

Advanced Digital Controls University of California, Los Angeles Department of Mechanical and Aerospace Engineering Report Author David Luong Winter 8 Table of Contents Abstract 4 Introduction..5 Physical

More information

### Mechatronics Engineering. Li Wen

Mechatronics Engineering Li Wen Bio-inspired robot-dc motor drive Unstable system Mirko Kovac,EPFL Modeling and simulation of the control system Problems 1. Why we establish mathematical model of the control

More information

### Time Response of Systems

Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =

More information

### Step Response for the Transfer Function of a Sensor

Step Response f the Transfer Function of a Sens G(s)=Y(s)/X(s) of a sens with X(s) input and Y(s) output A) First Order Instruments a) First der transfer function G(s)=k/(1+Ts), k=gain, T = time constant

More information

### Suboptimal adaptive control system for flight quality improvement

Suboptimal adaptive control system for flight uality improvement Andrzej Tomczyk Department of Avionics and Control, Faculty of Mechanical Engineering and Aeronautics Rzeszów University of Technology,

More information

### EE 224 Signals and Systems I Review 1/10

EE 224 Signals and Systems I Review 1/10 Class Contents Signals and Systems Continuous-Time and Discrete-Time Time-Domain and Frequency Domain (all these dimensions are tightly coupled) SIGNALS SYSTEMS

More information

### Dynamic Systems. Simulation of. with MATLAB and Simulink. Harold Klee. Randal Allen SECOND EDITION. CRC Press. Taylor & Francis Group

SECOND EDITION Simulation of Dynamic Systems with MATLAB and Simulink Harold Klee Randal Allen CRC Press Taylor & Francis Group Boca Raton London New York CRC Press is an imprint of the Taylor & Francis

More information

### Control System. Contents

Contents Chapter Topic Page Chapter- Chapter- Chapter-3 Chapter-4 Introduction Transfer Function, Block Diagrams and Signal Flow Graphs Mathematical Modeling Control System 35 Time Response Analysis of

More information

### Application Note #3413

Application Note #3413 Manual Tuning Methods Tuning the controller seems to be a difficult task to some users; however, after getting familiar with the theories and tricks behind it, one might find the

More information

### Performance of Feedback Control Systems

Performance of Feedback Control Systems Design of a PID Controller Transient Response of a Closed Loop System Damping Coefficient, Natural frequency, Settling time and Steady-state Error and Type 0, Type

More information

### Control Systems I. Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback. Readings: Emilio Frazzoli

Control Systems I Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 13, 2017 E. Frazzoli (ETH)

More information

### Iterative Learning Control (ILC)

Department of Automatic Control LTH, Lund University ILC ILC - the main idea Time Domain ILC approaches Stability Analysis Example: The Milk Race Frequency Domain ILC Example: Marine Vibrator Material:

More information

### Open Loop Tuning Rules

Open Loop Tuning Rules Based on approximate process models Process Reaction Curve: The process reaction curve is an approximate model of the process, assuming the process behaves as a first order plus

More information

### State feedback, Observer, and State feedback using an observer

Control and system theory 9/5/28 State feedback, Observer, and State feedback using an observer KOSEKI, Takafumi The University of Tokyo Fundamentals i What is a state observer? ii State feedback and classical

More information

### Andrea Zanchettin Automatic Control AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Linear systems (frequency domain)

1 AUTOMATIC CONTROL Andrea M. Zanchettin, PhD Spring Semester, 2018 Linear systems (frequency domain) 2 Motivations Consider an LTI system Thanks to the Lagrange s formula we can compute the motion of

More information

### ELECTRICAL ENGINEERING

ELECTRICAL ENGINEERING Subject Code: EE Course Structure Sections/Units Section A Unit 1 Unit 2 Unit 3 Unit 4 Unit 5 Unit 6 Unit 7 Section B Section C Section D Section E Section F Section G Section H

More information

### ScicosLab/Scicos for dummies like us Simone Mannori ScicosLab/Scicos developer

ScicosLab/Scicos for dummies like us 1 Digital control systems design and simulation: the Bermuda Triangle for engineers Computer science Control systems Physics 2 What ScicosLab is? An interpreted language

More information

### Singular Value Decomposition Analysis

Singular Value Decomposition Analysis Singular Value Decomposition Analysis Introduction Introduce a linear algebra tool: singular values of a matrix Motivation Why do we need singular values in MIMO control

More information

### Quanser NI-ELVIS Trainer (QNET) Series: QNET Experiment #02: DC Motor Position Control. DC Motor Control Trainer (DCMCT) Student Manual

Quanser NI-ELVIS Trainer (QNET) Series: QNET Experiment #02: DC Motor Position Control DC Motor Control Trainer (DCMCT) Student Manual Table of Contents 1 Laboratory Objectives1 2 References1 3 DCMCT Plant

More information

### Control Systems II. ETH, MAVT, IDSC, Lecture 4 17/03/2017. G. Ducard

Control Systems II ETH, MAVT, IDSC, Lecture 4 17/03/2017 Lecture plan: Control Systems II, IDSC, 2017 SISO Control Design 24.02 Lecture 1 Recalls, Introductory case study 03.03 Lecture 2 Cascaded Control

More information

### Laboratory III: Operational Amplifiers

Physics 33, Fall 2008 Lab III - Handout Laboratory III: Operational Amplifiers Introduction Operational amplifiers are one of the most useful building blocks of analog electronics. Ideally, an op amp would

More information

### FEEDBACK AND STABILITY

FEEDBCK ND STBILITY THE NEGTIVE-FEEDBCK LOOP x IN X OUT x S + x IN x OUT Σ Signal source _ β Open loop Closed loop x F Feedback network Output x S input signal x OUT x IN x F feedback signal x IN x S x

More information

### GATE 2017 ECE Session 3 Answer Key

1. Three pair cubical dice are thrown simultaneously. What is the probability that all three dice have same number is... Ans. ( 1 36 ). The one of the eigen value is real, rest are imaginary the real value

More information

### The different frequency components in a sound signal are processed first by the sensory system in the ear then some specialized cortical areas in the

Signals Represented in Terms of Magnitude Phase The physical interpretation of the Fourier transform of a time signal is its decomposition into a linear combination of frequency components with different

More information

### Vehicle longitudinal speed control

Vehicle longitudinal speed control Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin February 10, 2015 1 Introduction 2 Control concepts Open vs. Closed Loop Control

More information

### Digital Control Semester Project

Digital Control Semester Project Part I: Transform-Based Design 1 Introduction For this project you will be designing a digital controller for a system which consists of a DC motor driving a shaft with

More information