MCE693/793: Analysis and Control of Nonlinear Systems
|
|
- Emmeline Wade
- 5 years ago
- Views:
Transcription
1 MCE693/793: Analysis and Control of Nonlinear Systems Introduction to Describing Functions Hanz Richter Mechanical Engineering Department Cleveland State University
2 Introduction Frequency domain methods are essential for analysis and design of linear systems. Stability, robustness and performance can be studied using classical and modern tools such as 1. Bode, Nyquist and Nichols plots 2. Singular values and sensitivity plots (multivariable systems) 3. H control for performance 4. µ-synthesis analysis and control for robust performance 5. Quantitative feedback theory (frequency domain design on the Nichols chart) Describing functions provide an approximate method for frequency-domain analysis of nonlinear systems. We will focus on limit cycle detection/creation or avoidance. 2 / 31
3 Basic Idea In linear systems, transfer functions become frequency response functions when the Laplace variable s is restricted to the imaginary axis. That is, the Fourier transform is the Laplace transform restricted to the imaginary axis. Frequency response functions in linear systems are functions of frequency only: G(s) = s 1 1 w2, G(jw) = s+1 1+w 2 + 2w 1+w 2j In control systems containing nonlinearities, we seek to find a frequency response function for the nonlinear block. Due to nonlinearity, this function will be amplitude-dependent. 0 φ(e) e G(s) { N(A,w) { G(jw) 3 / 31
4 Memoryless vs. Memory, Soft vs. Hard Nonlinearities A memoryless nonlinearity produces an output which depends only on the current value of the input. The nonlinearity must then be a function of its argument: φ = φ(e(t)). Example: saturation. The output of a nonlineariry with memory depends on the current and past values of the inputs. The nonlinearity is not a function. Sometimes the terminology multiple-valued function is used. Example: hysteresis. A soft nonlinearity is continuous. The terminology is vague, sometimes soft nonlinearities are also required to be differentiable. Example: φ(e) = e 3. A hard nonlinearity has discontinuities. In some cases, a continuous function which is not differentiable everywhere is called a hard nonlinearity. Example: saturation. 4 / 31
5 Describing Function Concept When a sinuosoidal function is applied to a nonlinearity, harmonic distortion occurs and the output is no longer a sinusoid. Assuming the output remains periodic (always true for static nonlinearities, still true for some others), there is a Fourier (trigonometric) decomposition. The decomposition contains sinuosoidal components with frequencies which are integer multiples of the input frequency, and possibly a bias term (zero for odd functions). In particular, the fundamental component will have an amplitude which depends nonlinearly on the input amplitude. The principle behind describing functions is to retain only this fundamental output component as an approximation. A transfer function is formed by taking the ratio of the output component (expressed as a complex exponential) to the input sinusoid (also expressed as a complex exponential). This defines the sinusoidal input describing function (SIDF). 5 / 31
6 Calculating a simple SIDF As a first example, consider the cubic nonlinearity: φ(u) = u 3. Let the input be u(t) = Asin(wt) and express this as the complex exponential U(A,w) = Ae jwt = Acos(wt)+Ajsin(wt). The output of the nonlinearity is φ(u) = A 3 sin 3 (wt). We must find the fundamental component of the Fourier series expansion of this function. Recall that the expansion has the form: φ(t) = a n=1 a n cos( 2πnt T )+b nsin( 2πnt T ), where a n = 2 T b n = 2 T x0 +T x 0 x0 +T x 0 φ(t) cos( 2πnt T ) dt φ(t) sin( 2πnt T ) dt and a 0 /2 is the mean of φ(t) in one period. 6 / 31
7 Calculating a simple SIDF... Because φ is odd, a n = 0 for all n. We only need to calculate b 1. With T = 2π/w, the result is b 1 = 3 4 A3 The output sinusoid is then Φ = 3A3 4 ejwt Therefore the describing function for the cubic nonlinearity is N(A,w) = N(A) = 3A2 4 e0 All odd static nonlinearities will have a purely real and frequency-independent SIDF. 7 / 31
8 Example: Frequency-dependent SIDF Following the introductory example in S&L on application of the DF, consider the nonlinearity φ(x,ẋ) = x 2 ẋ associated with the van der Pol equation. Assume X = Ae jwt. Then the output of the nonlinearity is φ(x,ẋ) = A 3 wsin 2 (wt)cos(wt). Using trigonometric identities this is the same as φ(x,ẋ) = A3 w 4 (cos(wt) cos(3wt)) By the definition of SIDF, we retain only the fundamental component. Therefore Φ = A3 w 4 ej(wt π/2) The SIDF is then N(A,w) = Φ/U = A2 w 4 e jπ/2 = A2 w 4 j. In S&L the input is taken as ( x), so their N(A,w) has the opposite sign 8 / 31
9 Example: Memory Relay with Deadzone In class, we obtain the SIDF for a relay having the characteristic shown below: out D δ δ in D Show that N(A,w) = 4D πa 1 ( δ A )2 j 4Dδ. Because the nonlinearity is not πa 2 dynamic, frequency is absent. 9 / 31
10 Limit Cycle Analysis with Describing Functions If N(A,w) were a truly linear transfer function, we could investigate the presence of oscillations by using the idea of marginal stability and the Nyquist criterion. For example, linear analysis can provide the value of K at which the following system oscillates: 1 Nyquist Diagram 15 Root Locus ) 10 0 k 16 (s+2)(s 2 +2s+8) Imaginary Axis Imaginary Axis (seconds Real Axis Real Axis (seconds -1 ) 10 / 31
11 Review of the Nyquist Stability Criterion Complex Mapping Theorem: Let F(s) : C C be a mapping defined on a region of the complex plane and let C be a closed curve contained in that region. Assume: 1. F(s) has no more zeroes than poles and C does not pass through any of them. 2. F( ) is finite. 3. C encloses Z zeroes and P poles of F(s) Suppose C is parameterized with s = s(t) for t = [0,T] so that the curve is traversed clockwise. Then the image curve F(s) has N = Z P clockwise encirclements of the origin. (Animation) 11 / 31
12 Review of the Nyquist Stability Criterion... Nyquist s idea was to take C as an infinite half-circle encompassing the entire Im axis and the right half of the complex plane (where unstable poles lie). Since F( ) is a constant, only the Im axis portion determines how/if F(s) encircles the origin. Also, the Complex Mapping Theorem can be shifted : if 1+F(s) has Z P clockwise encirclements of the origin, then F(s) has the same number of encirclements of ( 1,0). 12 / 31
13 Review of the Nyquist Stability Criterion... In a unity negative feedback control system, the open-loop transfer function is L(s) = G(s)K(s) and the closed-loop transfer function is T(s) = L(s) 1+L(s) The closed-loop system is stable whenever L(s) doesn t have poles inside the Nyquist contour. Key idea: when L(s) is a rational function (L(s) = n(s) d(s) with n(s) and d(s) polynomials in s), the characteristic equation for T(s) has the form: L(s) = n(s) d(s) = 1 *Note*: The closed-loop poles are the zeroes of 1+L(s) and the open-loop poles (poles of L(s)) are also poles of 1+L(s). 13 / 31
14 Review of the Nyquist Stability Criterion... Apply the shifted Complex Mapping Theorem to 1+L(s): N is the number of clockwise origin encirclements of 1+L(s), or equivalently (we actually use this form) the number of clockwise encirclements of (-1,0). Note we evaluate L(s) on the Im axis only (this is called a Nyquist plot) Z is the number of zeroes of 1 + L(s) inside the Nyquist contour (unstable closed-loop poles) P is the number of poles of 1+L(s) (equivalently of L(s) and we use this) inside the Nyquist contour (unstable open-loop poles) Then Z = N +P 14 / 31
15 Example We use the Nyquist criterion to predict the stability of a unity feedback control 8 system with K(s) = k and G(s) =. Can the system be (s 2)(s 2 +2s+8) stabilized with proportional control? 0.25 Nyquist Diagram Imaginary Axis Real Axis Confirm analysis using the Routh-Hurwitz criterion. 15 / 31
16 Practical Significance of the Nyquist Stability Criterion Many systems are open-loop stable, so P = 0 is known even without a math model. We can use the Nyquist criterion to predict the critical proportional gain for closed-loop stability without having to model the system. Since P = 0, the system can be operate safely in open-loop. We can use a dynamic signal analyzer or offline signal processing to obtain Nyquist plot data. Then the stability of the closed-loop system can be studied using data from the open-loop system only. Specifically, we can find the critical proportional gain that scales the Nyquist plot to touch the (-1,0) point. Also, since pure phase lags (such as those introduced by delays) rotate the Nyquist plot, we can determine the maximum amount of delay that the system can tolerate before going unstable. 16 / 31
17 Nyquist-like analysis with Describing Functions Several conditions have to be satisfied for a DF-based frequency analysis to produce correct results (see S&L Sect ): 1. There s a single nonlinearity in the loop (unless 2 or more can be blended into one) 2. The nonlinearity is time-invariant (a relay with a time-dependent deadzone wouldn t) 3. Higher harmonics produced by the nonlinearity are attenuated by the linear dynamics following the nonlinearity (the linear TF must be low-pass). 4. The nonlinearity is odd, so that the Fourier expansion of its output has zero bias. Notes: Even when these conditions hold, the DF technique may produce erroneous results, particularly when predicting limit cycles. There are extended DF techniques relaxing assumptions 1,2 and 4, but they are cumbersome. 17 / 31
18 Example: van der Pol oscillator See S&L Sect The van der Pol equation is re-arranged as a negative unity feedback loop, where the nonlinear block is x 2 ẋ, followed by the unstable transfer function µ G(s) = s 2 µs+1 We had determined that the SIDF of the nonlinear block is N(A,w) = j A2 w 4. We can solve the characteristic equation for A and w: 1+N(A,w)G(jw) = 0 The solution is A = 2 and w = 1. Remarkably, the solution is independent of µ. The characteristic equation solution corresponds to a Nyquist plot of N(2,w)G(jw) crossing through ( 1,0) at a frequency of w = 1, as we can verify by obtaining a computer plot. 18 / 31
19 Observations We have solved the characteristic equation directly, but in the general case it won t be easy to find solutions for both A and w. A more exciting graphical method will be used. The existence of the limit cycle was correctly predicted (by the existence of solutions for A and w), but the accuracy in predicting w is reduced as µ is increased. Likewise, the limit cycle predicted by the DF is sinusoidal (circular phase orbit), while the shape of the actual cycle is not circular unless µ = / 31
20 Observations... The stability of the limit cycle can also be studied using the DF method. In S&L, the poles of the closed-loop characteristic equation are shown to be 1 8 µ(a2 4)± 1 64 µ2 (A 2 4) 2 1 If A > 2, the poles are stable, indicating L.C. stability (approaching the L.C.). If A < 2, the poles are unstable, indicating that trajectories grow towards the L.C. This indicates the L.C. is stable. Alternatively, from the Nyquist diagram, A < 2 causes the plot to enclose ( 1,0) (trajectories grow towards A = 2), while A > 2 stays clear of ( 1,0) (trajectories decrease towards A = 2). 20 / 31
21 Extended Nyquist Criterion To allow non-unity loops (for instance with sensor dynamics H(s)), the Nyquist criterion is simply applied to L(s) = K(s)G(s)H(s) Note that S&L uses L(s) = G(s)H(s), takes K(s) = K (constant) and uses the critical point 1/K, but K is allowed to be a complex number. But Figure 5.22 does not include K or H anymore. Our version: The characteristic equation is 1 + N(A, w)l(jw) = 0. Include everything linear in L(s) (L = KGH) Limit cycle detection is based on L(jw) = 1/N(A,w) Follow S&L Sect by thinking of their G as our L. 21 / 31
22 Limit Cycle Detection with DFs Since N(A,w) is both amplitude and frequency dependent in general, there will be a family of Nyquist plots of 1/N(A,w), where each curve has fixed w and is parameterized by A. Graph L(jw) and the family of 1/N(A,w) and look for intersections where frequencies coincide. When N(A,w) = N(A) (frequency-independent), there is only one curve parameterized by A, so the process is simpler. 22 / 31
23 Example Consider the nonlinear 3rd-order diff.eq.: d 3 x dt 3 +ẍ+2(1+kx2 )ẋ+3(1+x 2 )x = u with u = u 0 (constant). There is an additional condition for a periodic solution to be possible with u 0 0 (harmonic balance). Represent the equation as a unity negative feedback loop with a cascade of linear dynamics and nonlinearity: u 0 1 s 3 +s 2 +2s+3 x φ(x) 2kẋx 2 +3x 3 23 / 31
24 Example... We note that the nonlinearity is odd and the linear transfer function has low-pass properties. Let x = x 0 +Asin(wt). Work with trigonometric identities and leave out second and higher harmonics to find: ẋx 2 Awcos(wt)(x A2 4 ) x 3 3(x A2 4 )Asin(wt)+x 0(x A2 ) The output of the non-bias portion of nonlinearity is then φ(x,ẋ) = 9(x A2 4 )A[(2k/9)wcos(wt)+sin(wt)] As an operator on y = sin(wt) this is: [ 9(x A2 4 )A (2k/9) dy ] dt +y 24 / 31
25 Example... That is Φ(s) = 9(x A2 4 )A[(2k/9)s+1] In the frequency domain: N(A,w) = 9(x A2 4 )A[(2k/9)jw+1] Suppose that u 0 = 0. A DC harmonic balance will show that x 0 must be zero. To find potential limit cycles we solve (graphically or numerically): G(jw) = 1/N(A,w) Taking x 0 = 0 and k = 6, a solution is found: w = and A = / 31
26 Example... 1 Nyquist Diagram 0.8 w= Imaginary Axis w= Real Axis 26 / 31
27 Example Nyquist Diagram Imaginary Axis System: sysg Real: Imag: Frequency (rad/s): Real Axis 27 / 31
28 Example... DC Harmonic Balance Once the nonlinearity has been approximated with the SIDF, substitution of the candidate periodic solution into the diff. eq. gives: Aw 3 cos(wt) Aw 2 sin(wt)+2awcos(wt)+3(x 0 +Asin(wt)) = u 0 2kAw(x A2 4 ) 9(x2 0 + A2 4 )Asin(wt) 3x 0(x A2 ) In particular, the DC components must be matched: u 0 = 3x 0 (1+x A2 ) This provides information about the possible combinations of u 0 and x 0 associated with a periodic solution of amplitude A. Coefficient matching for sin(wt) and cos(wt) on both sides of the equation is essentially the same as the frequency domain analysis with the Nyquist criterion. We obtain a system of nonlinear equations. 28 / 31
29 Example: Limit Cycle Stability and Routh-Hurwitz Analysis The nonlinearity was replaced by the amplitude-dependent linear operator N(A, s). A characteristic equation can be written for the closed-loop system: [ s 3 +s k(x ] [ 4 A2 ) s+3 1+3(x ] 4 A2 ) = 0 The classical Routh-Hurwitz criterion can provide an equation for the onset of oscillations. Using the critical stability condition and the DC harmonic balance leads to the system of nonlinear equations: (2k 9)(x A2 ) = 1 3x 0 (1+x A2 ) = u 0 From the first equation, it s clear that k > 4.5 is necessary for a limit cycle to exist. The R-H criterion can also indicate the stability of any limit cycles. 29 / 31
30 DF Analysis of Forced Sinusoidal Input Response - Jump Resonance Consider the nonlinear system with sinusoidal input: ẋ = f(x)+u(t) u(t) = Ue jwt Use SIDF as necessary to represent f as f = N(X,w)x with x = Xe jwt and X i = A i e jφ i, also U = [u 1 u 2...u n ] T. Then, for each frequency w, the following must hold: X = [jwi N(X,w)] 1 U This is a system of 2n nonlinear equations with 2n unknowns that can be solved for the A i and φ i. When more than one solution exists for each frequency, the system exhibits jump resonance. 30 / 31
31 Example Consider the mass-spring-damper system with nonlinear stiffness: ẍ+2ζẋ+x(1+µx 2 ) = u Let the input be u(t) = U sin(wt). In class, we show that the output-input amplitude ratio satisfies X U = 1 (1+ 3µ 4 X2 w 2 )+(2ζw) 2 Up to 3 solutions are found for X in a frequency range. Amplitude ratio A B D C frequency 31 / 31
Control Systems. Frequency Method Nyquist Analysis.
Frequency Method Nyquist Analysis chibum@seoultech.ac.kr Outline Polar plots Nyquist plots Factors of polar plots PolarNyquist Plots Polar plot: he locus of the magnitude of ω vs. the phase of ω on polar
More informationNyquist Stability Criteria
Nyquist Stability Criteria Dr. Bishakh Bhattacharya h Professor, Department of Mechanical Engineering IIT Kanpur Joint Initiative of IITs and IISc - Funded by MHRD This Lecture Contains Introduction to
More informationAnalysis of SISO Control Loops
Chapter 5 Analysis of SISO Control Loops Topics to be covered For a given controller and plant connected in feedback we ask and answer the following questions: Is the loop stable? What are the sensitivities
More informationTable of Laplacetransform
Appendix Table of Laplacetransform pairs 1(t) f(s) oct), unit impulse at t = 0 a, a constant or step of magnitude a at t = 0 a s t, a ramp function e- at, an exponential function s + a sin wt, a sine fun
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 23: Drawing The Nyquist Plot Overview In this Lecture, you will learn: Review of Nyquist Drawing the Nyquist Plot Using
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 informationLinear Control Systems Lecture #3 - Frequency Domain Analysis. Guillaume Drion Academic year
Linear Control Systems Lecture #3 - Frequency Domain Analysis Guillaume Drion Academic year 2018-2019 1 Goal and Outline Goal: To be able to analyze the stability and robustness of a closed-loop system
More informationRobust Control 3 The Closed Loop
Robust Control 3 The Closed Loop Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /2/2002 Outline Closed Loop Transfer Functions Traditional Performance Measures Time
More informationControl Systems I. Lecture 9: The Nyquist condition
Control Systems I Lecture 9: The Nyquist condition Readings: Åstrom and Murray, Chapter 9.1 4 www.cds.caltech.edu/~murray/amwiki/index.php/first_edition Jacopo Tani Institute for Dynamic Systems and Control
More informationEC6405 - CONTROL SYSTEM ENGINEERING Questions and Answers Unit - I Control System Modeling Two marks 1. What is control system? A system consists of a number of components connected together to perform
More informationFrequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability
Lecture 6. Loop analysis of feedback systems 1. Motivation 2. Graphical representation of frequency response: Bode and Nyquist curves 3. Nyquist stability theorem 4. Stability margins Frequency methods
More informationControls Problems for Qualifying Exam - Spring 2014
Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function
More informationChemical Process Dynamics and Control. Aisha Osman Mohamed Ahmed Department of Chemical Engineering Faculty of Engineering, Red Sea University
Chemical Process Dynamics and Control Aisha Osman Mohamed Ahmed Department of Chemical Engineering Faculty of Engineering, Red Sea University 1 Chapter 4 System Stability 2 Chapter Objectives End of this
More informationIC6501 CONTROL SYSTEMS
DHANALAKSHMI COLLEGE OF ENGINEERING CHENNAI DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING YEAR/SEMESTER: II/IV IC6501 CONTROL SYSTEMS UNIT I SYSTEMS AND THEIR REPRESENTATION 1. What is the mathematical
More informationKINGS 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 informationVALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur
VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur 603 203. DEPARTMENT OF ELECTRONICS & COMMUNICATION ENGINEERING SUBJECT QUESTION BANK : EC6405 CONTROL SYSTEM ENGINEERING SEM / YEAR: IV / II year
More information2006 Fall. G(s) y = Cx + Du
1 Class Handout: Chapter 7 Frequency Domain Analysis of Feedback Systems 2006 Fall Frequency domain analysis of a dynamic system is very useful because it provides much physical insight, has graphical
More informationFrequency Response Techniques
4th Edition T E N Frequency Response Techniques SOLUTION TO CASE STUDY CHALLENGE Antenna Control: Stability Design and Transient Performance First find the forward transfer function, G(s). Pot: K 1 = 10
More informationD(s) G(s) A control system design definition
R E Compensation D(s) U Plant G(s) Y Figure 7. A control system design definition x x x 2 x 2 U 2 s s 7 2 Y Figure 7.2 A block diagram representing Eq. (7.) in control form z U 2 s z Y 4 z 2 s z 2 3 Figure
More informationEngineering Tripos Part IIB Nonlinear Systems and Control. Handout 3: Describing Functions
Engineering Tripos Part IIB Module 4F2 Nonlinear Systems and Control Handout 3: Describing Functions 1 Harmonic balance The describing function method (also called the method of harmonic balance) uses
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 23: Drawing The Nyquist Plot Overview In this Lecture, you will learn: Review of Nyquist Drawing the Nyquist Plot Using the
More informationDigital Control Systems
Digital Control Systems Lecture Summary #4 This summary discussed some graphical methods their use to determine the stability the stability margins of closed loop systems. A. Nyquist criterion Nyquist
More informationControl Systems I. Lecture 9: The Nyquist condition
Control Systems I Lecture 9: The Nyquist condition adings: Guzzella, Chapter 9.4 6 Åstrom and Murray, Chapter 9.1 4 www.cds.caltech.edu/~murray/amwiki/index.php/first_edition Emilio Frazzoli Institute
More informationRichiami 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 informationROOT 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 informationAndrea 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 informationTime Response Analysis (Part II)
Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary
More informationFATIMA MICHAEL COLLEGE OF ENGINEERING & TECHNOLOGY
FATIMA MICHAEL COLLEGE OF ENGINEERING & TECHNOLOGY Senkottai Village, Madurai Sivagangai Main Road, Madurai - 625 020. An ISO 9001:2008 Certified Institution DEPARTMENT OF ELECTRONICS AND COMMUNICATION
More informationThe Nyquist criterion relates the stability of a closed system to the open-loop frequency response and open loop pole location.
Introduction to the Nyquist criterion The Nyquist criterion relates the stability of a closed system to the open-loop frequency response and open loop pole location. Mapping. If we take a complex number
More informationRecitation 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 informationChapter 7 : Root Locus Technique
Chapter 7 : Root Locus Technique By Electrical Engineering Department College of Engineering King Saud University 1431-143 7.1. Introduction 7.. Basics on the Root Loci 7.3. Characteristics of the Loci
More informationControl Systems. Root Locus & Pole Assignment. L. Lanari
Control Systems Root Locus & Pole Assignment L. Lanari Outline root-locus definition main rules for hand plotting root locus as a design tool other use of the root locus pole assignment Lanari: CS - Root
More informationStability of Feedback Control Systems: Absolute and Relative
Stability of Feedback Control Systems: Absolute and Relative Dr. Kevin Craig Greenheck Chair in Engineering Design & Professor of Mechanical Engineering Marquette University Stability: Absolute and Relative
More informationMEM 355 Performance Enhancement of Dynamical Systems
MEM 355 Performance Enhancement of Dynamical Systems Frequency Domain Design Intro Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /5/27 Outline Closed Loop Transfer
More informationUnit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace
Unit : Modeling in the Frequency Domain Part : Engineering 81: Control Systems I Faculty of Engineering & Applied Science Memorial University of Newfoundland January 1, 010 1 Pair Table Unit, Part : Unit,
More informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems 1 Analysis of Discrete-Time Systems Overview Stability Sensitivity and Robustness Controllability, Reachability, Observability, and Detectabiliy TU Berlin Discrete-Time
More informationControl Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 7: Feedback and the Root Locus method Readings: Jacopo Tani Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 2, 2018 J. Tani, E. Frazzoli (ETH) Lecture 7:
More informationLaplace 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 informationECEN 605 LINEAR SYSTEMS. Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability 1/27
1/27 ECEN 605 LINEAR SYSTEMS Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability Feedback System Consider the feedback system u + G ol (s) y Figure 1: A unity feedback system
More informationControl Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 6: Poles and Zeros Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 27, 2017 E. Frazzoli (ETH) Lecture 6: Control Systems I 27/10/2017
More informationCONTROL SYSTEM STABILITY. CHARACTERISTIC EQUATION: The overall transfer function for a. where A B X Y are polynomials. Substitution into the TF gives:
CONTROL SYSTEM STABILITY CHARACTERISTIC EQUATION: The overall transfer function for a feedback control system is: TF = G / [1+GH]. The G and H functions can be put into the form: G(S) = A(S) / B(S) H(S)
More information(Continued on next page)
(Continued on next page) 18.2 Roots of Stability Nyquist Criterion 87 e(s) 1 S(s) = =, r(s) 1 + P (s)c(s) where P (s) represents the plant transfer function, and C(s) the compensator. The closedloop characteristic
More informationControl Systems. University Questions
University Questions UNIT-1 1. Distinguish between open loop and closed loop control system. Describe two examples for each. (10 Marks), Jan 2009, June 12, Dec 11,July 08, July 2009, Dec 2010 2. Write
More informationGATE EE Topic wise Questions SIGNALS & SYSTEMS
www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)
More informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems TU Berlin Discrete-Time Control Systems 2 Stability Definitions We define stability first with respect to changes in the initial conditions Analysis of Discrete-Time
More informationECE 345 / ME 380 Introduction to Control Systems Lecture Notes 8
Learning Objectives ECE 345 / ME 380 Introduction to Control Systems Lecture Notes 8 Dr. Oishi oishi@unm.edu November 2, 203 State the phase and gain properties of a root locus Sketch a root locus, by
More informationSoftware Engineering/Mechatronics 3DX4. Slides 6: Stability
Software Engineering/Mechatronics 3DX4 Slides 6: Stability Dr. Ryan Leduc Department of Computing and Software McMaster University Material based on lecture notes by P. Taylor and M. Lawford, and Control
More informationLecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.
ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition
More informationMAK 391 System Dynamics & Control. Presentation Topic. The Root Locus Method. Student Number: Group: I-B. Name & Surname: Göksel CANSEVEN
MAK 391 System Dynamics & Control Presentation Topic The Root Locus Method Student Number: 9901.06047 Group: I-B Name & Surname: Göksel CANSEVEN Date: December 2001 The Root-Locus Method Göksel CANSEVEN
More informationRoot Locus Methods. The root locus procedure
Root Locus Methods Design of a position control system using the root locus method Design of a phase lag compensator using the root locus method The root locus procedure To determine the value of the gain
More informationCourse roadmap. Step response for 2nd-order system. Step response for 2nd-order system
ME45: Control Systems Lecture Time response of nd-order systems Prof. Clar Radcliffe and Prof. Jongeun Choi Department of Mechanical Engineering Michigan State University Modeling Laplace transform Transfer
More informationLecture 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 informationCONTROL * ~ SYSTEMS ENGINEERING
CONTROL * ~ SYSTEMS ENGINEERING H Fourth Edition NormanS. Nise California State Polytechnic University, Pomona JOHN WILEY& SONS, INC. Contents 1. Introduction 1 1.1 Introduction, 2 1.2 A History of Control
More informationECE317 : Feedback and Control
ECE317 : Feedback and Control Lecture : Routh-Hurwitz stability criterion Examples Dr. Richard Tymerski Dept. of Electrical and Computer Engineering Portland State University 1 Course roadmap Modeling
More informationCONTROL SYSTEMS ENGINEERING Sixth Edition International Student Version
CONTROL SYSTEMS ENGINEERING Sixth Edition International Student Version Norman S. Nise California State Polytechnic University, Pomona John Wiley fir Sons, Inc. Contents PREFACE, vii 1. INTRODUCTION, 1
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 informationDESIGN USING TRANSFORMATION TECHNIQUE CLASSICAL METHOD
206 Spring Semester ELEC733 Digital Control System LECTURE 7: DESIGN USING TRANSFORMATION TECHNIQUE CLASSICAL METHOD For a unit ramp input Tz Ez ( ) 2 ( z ) D( z) G( z) Tz e( ) lim( z) z 2 ( z ) D( z)
More informationEEE 184: Introduction to feedback systems
EEE 84: Introduction to feedback systems Summary 6 8 8 x 7 7 6 Level() 6 5 4 4 5 5 time(s) 4 6 8 Time (seconds) Fig.. Illustration of BIBO stability: stable system (the input is a unit step) Fig.. step)
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 22: The Nyquist Criterion Overview In this Lecture, you will learn: Complex Analysis The Argument Principle The Contour
More informationAlireza Mousavi Brunel University
Alireza Mousavi Brunel University 1 » Control Process» Control Systems Design & Analysis 2 Open-Loop Control: Is normally a simple switch on and switch off process, for example a light in a room is switched
More informationSTABILITY ANALYSIS TECHNIQUES
ECE4540/5540: Digital Control Systems 4 1 STABILITY ANALYSIS TECHNIQUES 41: Bilinear transformation Three main aspects to control-system design: 1 Stability, 2 Steady-state response, 3 Transient response
More informationPrüfung Regelungstechnik I (Control Systems I) Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 29. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid
More informationSTABILITY OF CLOSED-LOOP CONTOL SYSTEMS
CHBE320 LECTURE X STABILITY OF CLOSED-LOOP CONTOL SYSTEMS Professor Dae Ryook Yang Spring 2018 Dept. of Chemical and Biological Engineering 10-1 Road Map of the Lecture X Stability of closed-loop control
More informationCourse Summary. The course cannot be summarized in one lecture.
Course Summary Unit 1: Introduction Unit 2: Modeling in the Frequency Domain Unit 3: Time Response Unit 4: Block Diagram Reduction Unit 5: Stability Unit 6: Steady-State Error Unit 7: Root Locus Techniques
More informationBasic Procedures for Common Problems
Basic Procedures for Common Problems ECHE 550, Fall 2002 Steady State Multivariable Modeling and Control 1 Determine what variables are available to manipulate (inputs, u) and what variables are available
More informationDesign and Tuning of Fractional-order PID Controllers for Time-delayed Processes
Design and Tuning of Fractional-order PID Controllers for Time-delayed Processes Emmanuel Edet Technology and Innovation Centre University of Strathclyde 99 George Street Glasgow, United Kingdom emmanuel.edet@strath.ac.uk
More informationH(s) = s. a 2. H eq (z) = z z. G(s) a 2. G(s) A B. s 2 s(s + a) 2 s(s a) G(s) 1 a 1 a. } = (z s 1)( z. e ) ) (z. (z 1)(z e at )(z e at )
.7 Quiz Solutions Problem : a H(s) = s a a) Calculate the zero order hold equivalent H eq (z). H eq (z) = z z G(s) Z{ } s G(s) a Z{ } = Z{ s s(s a ) } G(s) A B Z{ } = Z{ + } s s(s + a) s(s a) G(s) a a
More informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Frequency Response-Design Method
.. AERO 422: Active Controls for Aerospace Vehicles Frequency Response- Method Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. ... Response to
More informationSome special cases
Lecture Notes on Control Systems/D. Ghose/2012 87 11.3.1 Some special cases Routh table is easy to form in most cases, but there could be some cases when we need to do some extra work. Case 1: The first
More informationCHAPTER # 9 ROOT LOCUS ANALYSES
F K א CHAPTER # 9 ROOT LOCUS ANALYSES 1. Introduction The basic characteristic of the transient response of a closed-loop system is closely related to the location of the closed-loop poles. If the system
More informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries
. AERO 632: of Advance Flight Control System. Preliminaries Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Preliminaries Signals & Systems Laplace
More information20. The pole diagram and the Laplace transform
95 0. The pole diagram and the Laplace transform When working with the Laplace transform, it is best to think of the variable s in F (s) as ranging over the complex numbers. In the first section below
More informationECE 486 Control Systems
ECE 486 Control Systems Spring 208 Midterm #2 Information Issued: April 5, 208 Updated: April 8, 208 ˆ This document is an info sheet about the second exam of ECE 486, Spring 208. ˆ Please read the following
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Lyapunov Stability - I Hanz Richter Mechanical Engineering Department Cleveland State University Definition of Stability - Lyapunov Sense Lyapunov
More informationCourse Outline. Closed Loop Stability. Stability. Amme 3500 : System Dynamics & Control. Nyquist Stability. Dr. Dunant Halim
Amme 3 : System Dynamics & Control Nyquist Stability Dr. Dunant Halim Course Outline Week Date Content Assignment Notes 1 5 Mar Introduction 2 12 Mar Frequency Domain Modelling 3 19 Mar System Response
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 5. 2. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -
More informationOutline. Control systems. Lecture-4 Stability. V. Sankaranarayanan. V. Sankaranarayanan Control system
Outline Control systems Lecture-4 Stability V. Sankaranarayanan Outline Outline 1 Outline Outline 1 2 Concept of Stability Zero State Response: The zero-state response is due to the input only; all the
More information1-DOF Vibration Characteristics. MCE371: Vibrations. Prof. Richter. Department of Mechanical Engineering. Handout 7 Fall 2017
MCE371: Vibrations Prof. Richter Department of Mechanical Engineering Handout 7 Fall 2017 Free Undamped Vibration Follow Palm, Sect. 3.2, 3.3 (pp 120-138), 3.5 (pp 144-151), 3.8 (pp. 167-169) The equation
More informationECE317 : Feedback and Control
ECE317 : Feedback and Control Lecture : Stability Routh-Hurwitz stability criterion Dr. Richard Tymerski Dept. of Electrical and Computer Engineering Portland State University 1 Course roadmap Modeling
More informationEC 8391-CONTROL SYSTEMS ENGINEERING. Questions and Answers PART-A. Unit - I Systems Components And Their Representation
EC 8391-CONTROL SYSTEMS ENGINEERING Questions and Answers PART-A Unit - I Systems Components And Their Representation 1. What is control system? A system consists of a number of components connected together
More information1-DOF Forced Harmonic Vibration. MCE371: Vibrations. Prof. Richter. Department of Mechanical Engineering. Handout 8 Fall 2011
MCE371: Vibrations Prof. Richter Department of Mechanical Engineering Handout 8 Fall 2011 Harmonic Forcing Functions Transient vs. Steady Vibration Follow Palm, Sect. 4.1, 4.9 and 4.10 Harmonic forcing
More informationSome solutions of the written exam of January 27th, 2014
TEORIA DEI SISTEMI Systems Theory) Prof. C. Manes, Prof. A. Germani Some solutions of the written exam of January 7th, 0 Problem. Consider a feedback control system with unit feedback gain, with the following
More informationReturn Difference Function and Closed-Loop Roots Single-Input/Single-Output Control Systems
Spectral Properties of Linear- Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2018! Stability margins of single-input/singleoutput (SISO) systems! Characterizations
More informationFEEDBACK and CONTROL SYSTEMS
SCHA UM'S OUTLINE OF THEORY AND PROBLEMS OF FEEDBACK and CONTROL SYSTEMS Second Edition CONTINUOUS (ANALOG) AND DISCRETE (DIGITAL) JOSEPH J. DiSTEFANO, III, PhD. Departments of Computer Science and Mediane
More informationVALLIAMMAI 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 informationSuppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber
Suppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber J.C. Ji, N. Zhang Faculty of Engineering, University of Technology, Sydney PO Box, Broadway,
More informationMAE 143B - Homework 9
MAE 143B - Homework 9 7.1 a) We have stable first-order poles at p 1 = 1 and p 2 = 1. For small values of ω, we recover the DC gain K = lim ω G(jω) = 1 1 = 2dB. Having this finite limit, our straight-line
More informationFeedback Control of Dynamic Systems
THIRD EDITION Feedback Control of Dynamic Systems Gene F. Franklin Stanford University J. David Powell Stanford University Abbas Emami-Naeini Integrated Systems, Inc. TT Addison-Wesley Publishing Company
More informationTest 2 SOLUTIONS. ENGI 5821: Control Systems I. March 15, 2010
Test 2 SOLUTIONS ENGI 5821: Control Systems I March 15, 2010 Total marks: 20 Name: Student #: Answer each question in the space provided or on the back of a page with an indication of where to find the
More informationControl Systems Engineering ( Chapter 8. Root Locus Techniques ) Prof. Kwang-Chun Ho Tel: Fax:
Control Systems Engineering ( Chapter 8. Root Locus Techniques ) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 02-760-4253 Fax:02-760-4435 Introduction In this lesson, you will learn the following : The
More informationIndex. Index. More information. in this web service Cambridge University Press
A-type elements, 4 7, 18, 31, 168, 198, 202, 219, 220, 222, 225 A-type variables. See Across variable ac current, 172, 251 ac induction motor, 251 Acceleration rotational, 30 translational, 16 Accumulator,
More information10 Transfer Matrix Models
MIT EECS 6.241 (FALL 26) LECTURE NOTES BY A. MEGRETSKI 1 Transfer Matrix Models So far, transfer matrices were introduced for finite order state space LTI models, in which case they serve as an important
More informationCollocated versus non-collocated control [H04Q7]
Collocated versus non-collocated control [H04Q7] Jan Swevers September 2008 0-0 Contents Some concepts of structural dynamics Collocated versus non-collocated control Summary This lecture is based on parts
More informationK(s +2) s +20 K (s + 10)(s +1) 2. (c) KG(s) = K(s + 10)(s +1) (s + 100)(s +5) 3. Solution : (a) KG(s) = s +20 = K s s
321 16. Determine the range of K for which each of the following systems is stable by making a Bode plot for K = 1 and imagining the magnitude plot sliding up or down until instability results. Verify
More informationThe Nyquist Stability Test
Handout X: EE24 Fall 2002 The Nyquist Stability Test.0 Introduction With negative feedback, the closed-loop transfer function A(s) approaches the reciprocal of the feedback gain, f, as the magnitude of
More informationEC CONTROL SYSTEM UNIT I- CONTROL SYSTEM MODELING
EC 2255 - CONTROL SYSTEM UNIT I- CONTROL SYSTEM MODELING 1. What is meant by a system? It is an arrangement of physical components related in such a manner as to form an entire unit. 2. List the two types
More informationINSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad
INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad - 500 043 Electrical and Electronics Engineering TUTORIAL QUESTION BAN Course Name : CONTROL SYSTEMS Course Code : A502 Class : III
More informationIntroduction to Nonlinear Control Lecture # 4 Passivity
p. 1/6 Introduction to Nonlinear Control Lecture # 4 Passivity È p. 2/6 Memoryless Functions ¹ y È Ý Ù È È È È u (b) µ power inflow = uy Resistor is passive if uy 0 p. 3/6 y y y u u u (a) (b) (c) Passive
More informationCh 14: Feedback Control systems
Ch 4: Feedback Control systems Part IV A is concerned with sinle loop control The followin topics are covered in chapter 4: The concept of feedback control Block diaram development Classical feedback controllers
More information