Lecture 6 Describing function analysis
|
|
- Godfrey Warren
- 6 years ago
- Views:
Transcription
1 Lecture 6 Describing function analysis Today s Goal: To be able to Derive describing functions for static nonlinearities Predict stability and existence of periodic solutions through describing function analysis Material: Chapter 14 in Glad & Ljung Chapter 7.2 (pp.28 29) in Khalil Lecture notes
2 Course Outline Lecture 1-3 Lecture 2-6 Lecture 7-8 Lecture 9-13 Lecture 14 Modelling and basic phenomena (linearization, phase plane, limit cycles) Analysis methods (Lyapunov, circle criterion, describing functions) Common nonlinearities (Saturation, friction, backlash, quantization) Design methods (Lyapunov methods, Backstepping, Optimal control) Summary
3 Example: saturated sinusoidals input signal output signal time [s] time [s] Sine Wave Saturation Scope The effective gain (the ratio sat(asinωt) ) varies with the Asinωt input signal amplitude A.
4 Motivating Example r e u y G(s) 2 1 u y G(s)= 4 andu=sat(e) gives stable oscillation forr=. s(s+1) 2 How can the oscillation be predicted? Q: What is the amplitude/topvalue ofuandy? What is the frequency?
5 Recall the Nyquist Theorem Assume G(s) stable, and k is positive gain. The closed-loop system is unstable if the point 1/k is encircled by G(iω) The closed-loop system is stable if the point 1/k is not encircled by G(iω) e u y k G(s) 1/k G(iω)
6 Heuristic reasoning: Motivating Example (cont d) For what frequency and what amplitude is "the loop gain"f G= 1? Introduce N(A) as an amplitude dependent approximation of the nonlinearity f( ). G(iω) e u y f( ) G(s) 1/N(A) A y=g(iω)u G(iω)N(A)y G(iω)= 1 N(A)
7 Motivating Example (cont d) 2 1/N(A) 4 6 G(iω) Introduce N(A) as an amplitude dependent gain-approximation of the nonlinearity f( ). Heuristic reasoning: For what frequency and what amplitude is "the loop gain"n(a) G(iw)= 1? The intersection of the 1/N(A) and the Nyquist curveg(iω) predicts amplitude and frequency.
8 How do we derive the describing function N(A)? Does the intersection predict a stable oscillation? Are the estimated amplitude and frequency accurate?
9 Fourier Series Every periodic functionu(t)=u(t+t) has a Fourier series expansion u(t)= a 2 + (a n cosnωt+b n sinnωt) n=1 = a 2 + a 2 n+b 2 nsin[nωt+arctan(a n /b n )] n=1 where ω=2π/t and a n = 2 T T u(t)cosnωtdt b n = 2 T T u(t)sinnωtdt Note: Sometimes we make the change of variablet φ/ω
10 The Fourier Coefficients are Optimal The finite expansion solves û k (t)= a k 2 + (a n cosnωt+b n sinnωt) 2 min {â n,ˆb n } 1 n k T T n=1 if{a n,b n } are the Fourier coefficients. [ u(t) ûk (t) ] 2 dt T= 2π ω
11 The Key Idea e u y N.L. G(s) Assumee(t)=Asinωt andu(t) periodic. Then u(t)= a 2 + a 2 n +b2 n sin[nωt+arctan(a n/b n )] n=1 If G(inω) G(iω) forn=2,3,... anda =, then y(t) G(iω) a 2 1 +b2 1 sin[ωt+arctan(a 1/b 1 )+argg(iω)] Find periodic solution by matching coefficients in y = e.
12 Definition of Describing Function The describing function of e(t) N.L. u(t) is a 1 (ω):= ω π 2π/ω N(A,ω):= b 1(ω)+ia 1 (ω) A u(t)cos(ωt)dt b 1 (ω):= ω π 2π/ω u(t) sin(ω t)dt where u(t) is the output corresponding to e(t):= A sin(ω t) IfGis low pass anda =, then û 1 (t)= N(A,ω) Asin[ωt+argN(A,ω)] can be used instead ofu(t) to analyze the system. Amplitude dependent gain and phase shift!
13 Idea: Use the describing function to approximate the part of the signal coming out from the nonlinearity which will survive through the low-pass filtering linear system. e(t)=asinωt=im(ae iωt ) e(t) u(t) N.L. u(t)= a 2 + (a n cosnωt+b n sinnωt) n=1 e(t) N(A,ω) u 1 (t) u 1 (t)=a 1 cos(ωt)+b 1 sin(ωt) =Im(N(A,ω)Ae iωt ) where the describing function is defined as N(A,ω)= b 1(ω)+ia 1 (ω) A = U(iω) N(A,ω)E(iω)
14 Existence of Limit Cycles G(iω) e u y f( ) G(s) 1/N(A) A y=g(iω)u G(iω)N(A)y G(iω)= 1 N(A) The intersections ofg(iω) and 1/N(A) give ω anda for possible limit cycles.
15 Describing Function for a Relay H u 1.5 e u e H φ a 1 = 1 π 2π 2π u(φ)cos φdφ= π b 1 = 1 π u(φ)sin φdφ= 2 π Hsinφdφ= 4H π The describing function for a relay is thusn(a)= 4H πa.
16 Describing Function for Odd Static Nonlinearities Assumef( ) and ( ) are odd static nonlinearities (i.e., f( e)= f(e)) with describing functionsn f andn. Then, ImN f (A,ω)=, coeff.(a 1 ) N f (A,ω)=N f (A) N αf (A)= αn f (A) N f+ (A)=N f (A)+N (A)
17 Limit Cycle in Relay Feedback System r e u y G(s) / N(A) G(iω) G(s)= 3 (s+1) 3 with feedback u= sgny 3/8= 1/N(A)= πa/4 A=12/8π.48 G(iω)= 3/8 ω 1.7, T=2π/ω=3.7
18 Limit Cycle in Relay Feedback System (cont d) The prediction via the describing function agrees very well with the true oscillations: 1.5 y u G filters out almost all higher-order harmonics.
19 Describing Function for a Saturation u H e D D e.4.2 u.2.4 H.6.8 φ Lete(t)=Asinωt=Asinφ. First seth=d. IfA Dthen N(A)=1, ifa>dthen for φ (,π) { Asinφ, φ (,φ ) (π φ u(φ)=,π) D, φ (φ,π φ ) where φ =arcsind/a.
20 Describing Function for a Saturation (cont d) a 1 = 1 π b 1 = 1 π 2π 2π φ u(φ)cos φdφ= u(φ)sin φdφ= 4 π = 4A sin 2 φdφ+ 4D π π = A ) (2φ +sin2φ π π/2 π/2 u(φ)sin φdφ φ sinφdφ
21 Describing Function for a Saturation (cont d) IfH=D N(A)= 1 ) (2φ +sin2φ, A D π ForH =Dthe rulen αf (A)= αn f (A) gives N(A)= H ) (2φ +sin2φ, A D Dπ N(A) forh=d=1 NOTE: dependance of A shows up in φ =arcsind/a
22 3 minute exercise: What oscillation amplitude and frequency do the describing function analysis predict for the Motivating Example?
23 Solution: Find ω anda such thatg(iω) N(A)= 1; AsN(A) is positive and real valued, find ω s.t. argg(iω)= π= π 2 2arctan ω 1 = π= ω=1. which gives a period time of 6.28 seconds. G(1.i) N(A)=2 N(A)=1= N(A)=.5. To find the amnplitude A, either ) Alt. 1 Solve (numerically)n(a)= 1 1 π (2φ +sin2φ =.5, where φ =arcsin(1/a) Alt. 2 From the diagram ofn(a) one can finda N(A) forh=d=
24 The Nyquist Theorem Assume G(s) stable, and k is positive gain. The closed-loop system is unstable if the point 1/k is encircled by G(iω) The closed-loop system is stable if the point 1/k is not encircled by G(iω) e u y k G(s) 1/k G(iω)
25 How to Predict Stability of Limit Cycles AssumeG(s) stable. For a givena=a : A increases if the point 1/N f (A ) is encircled byg(iω) A decreases otherwise e u y f G(s) G(iω) 1/N(A) A stable limit cycle is predicted
26 How to Predict Stability of Limit Cycles G(Ω) 1/ N(A) An unstable limit cycle is predicted An intersection with amplitudea is unstable ifa<a gives decreasing amplitude anda>a gives increasing.
27 Stable Periodic Solution in Relay System r e u y G(s) G(iω) 1/N(A) G(s)= (s+1)2 (s+1) 3 with feedback u= sgny gives one stable and one unstable limit cycle. The left most intersection corresponds to the stable one.
28 Periodic Solutions in Relay System The relay gainn(a) is higher for smalla: Big amplitudes Small relay gain No encirclement Shrinking amplitudes Growing relay gain One encirclement Growing amplitudes Shrinking relay gain Small amplitudes High relay gain No encirclement Shrinking amplitudes Growing relay gain Stable periodic orbit Unstable periodic orbit
29 Σ Automatic Tuning of PID Controller Period and amplitude of relay feedback limit cycle can be used for autotuning. PID A T u Process y Relay 1
30 Describing Function for a dead-zone relay D u e e.4.2 u D Lete(t)=Asinωt=Asinφ. Then for φ (,π) {, φ (,φ ) u(φ)= D, φ (φ,π φ ) where φ =arcsind/a (ifa D)
31 Describing Function for a dead-zone relay cont d. a 1 = 1 π b 1 = 1 π = 4D π 2π 2π u(φ)cos φdφ= u(φ)sinφdφ= 4 π cosφ = 4D π π/2 1 D 2 /A 2 φ Dsinφdφ {, A<D N(A)= 4 1 D 2 /A 2, A D πa
32 Plot of Describing Function for dead-zone relay N(A) ford= Notice that N(A) 1.3/A for large amplitudes
33 Pitfalls Describing function analysis can give erroneous results. DF analysis may predict a limit cycle, even if it does not exist. A limit cycle may exist, even if DF analysis does not predict it. The predicted amplitude and frequency are only approximations and can be far from the true values.
34 Example The control of output powerx(t) from a mobile telephone is critical for good performance. One does not want to use too large power since other channels are affected and the battery length is decreased. Information about received power is sent back to the transmitter and is used for power control. A very simple scheme is given by ẋ(t)= αu(t) u(t)= signy(t L), α, β> y(t)= βx(t). Use describing function analysis to predict possible limit cycle amplitude and period of y. (The signals have been transformed so x=corresponds to nominal output power) y(t)=βx(t) x(t) y(t) y(t L) α x(t) s β e sl y(t L) u(t)
35 The system can be written as a negative feedback loop with G(s)= e sl α β s and a relay with amplitude 1. The describing function of a relay satisfies 1/N(A)= πa/4 hence we are interesting ing(iω) on the negative real axis. A stable intersection is given by π= argg(iω)= π/2 ωl which gives ω= π/(2l).this gives πa 4 = G(iω) = α β ω =2Lα β π and hencea=8lα β/π 2. The period is given by T=2π/ω=4L. (More exact analysis gives the true values A= α βl andt=4l, so the prediction is quite good.)
36 Accuracy of Describing Function Analysis Control loop with frictionf= sgny: F Friction y ref _ C u _ G y Corresponds to G 1+GC = s(s b) s 3 +2s 2 +2s+1 with feedback u= sgny The oscillation depends on the zero ats=b.
37 Accuracy of Describing Function Analysis 1.2 b=1/3 1.8 b=4/ DF predicts period times and ampl.(t,a) b=4/3 =(11.4,1.) and(t,a) b=1/3 =(17.3,.23) y b=4/3 Simulation: y b=1/ (T,A) b=4/3 =(12,1.1) (T,A) b=1/3 =(22,.28) Accurate results only if y is sinusoidal!
38 Analysis of Oscillations A summary There exist both time-domain and frequency-domain methods to analyze oscillations. Time-domain: Poincaré maps and Lyapunov functions Rigorous results but hard to use for large problems Frequency-domain: Describing function analysis Approximate results Powerful graphical methods
39 Today s Goal To be able to Derive describing functions for static nonlinearities Predict stability and existence of periodic solutions through describing function analysis
40 Next Lecture Saturation and antiwindup compensation Lyapunov analysis of phase locked loops Friction compensation
Lecture 6 Describing function analysis
Lecture 6 Describing function analysis Today s Goal: To be able to Derive describing functions for static nonlinearities Predict stability and existence of periodic solutions through describing function
More informationLecture 7: Anti-windup and friction compensation
Lecture 7: Anti-windup and friction compensation Compensation for saturations (anti-windup) Friction models Friction compensation Material Lecture slides Course Outline Lecture 1-3 Lecture 2-6 Lecture
More informationLecture 5 Input output stability
Lecture 5 Input output stability or How to make a circle out of the point 1+0i, and different ways to stay away from it... k 2 yf(y) r G(s) y k 1 y y 1 k 1 1 k 2 f( ) G(iω) Course Outline Lecture 1-3 Lecture
More informationAutomatic Control III (Reglerteknik III) fall Nonlinear systems, Part 3
Automatic Control III (Reglerteknik III) fall 20 4. Nonlinear systems, Part 3 (Chapter 4) Hans Norlander Systems and Control Department of Information Technology Uppsala University OSCILLATIONS AND DESCRIBING
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 informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 21: Stability Margins and Closing the Loop Overview In this Lecture, you will learn: Closing the Loop Effect on Bode Plot Effect
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 informationLecture 9 Nonlinear Control Design
Lecture 9 Nonlinear Control Design Exact-linearization Lyapunov-based design Lab 2 Adaptive control Sliding modes control Literature: [Khalil, ch.s 13, 14.1,14.2] and [Glad-Ljung,ch.17] Course Outline
More informationOverview Lecture 1. Nonlinear Control and Servo systems Lecture 1. Course Material. Course Goal. Course Material, cont. Lectures and labs.
Overview Lecture Nonlinear Control and Servo sstems Lecture Anders Rantzer Practical information Course contents Nonlinear control sstems phenomena Nonlinear differential equations Automatic Control LTH
More informationStep Response Analysis. Frequency Response, Relation Between Model Descriptions
Step Response Analysis. Frequency Response, Relation Between Model Descriptions Automatic Control, Basic Course, Lecture 3 November 9, 27 Lund University, Department of Automatic Control Content. Step
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Introduction to Describing Functions Hanz Richter Mechanical Engineering Department Cleveland State University Introduction Frequency domain methods
More informationLecture 9 Nonlinear Control Design. Course Outline. Exact linearization: example [one-link robot] Exact Feedback Linearization
Lecture 9 Nonlinear Control Design Course Outline Eact-linearization Lyapunov-based design Lab Adaptive control Sliding modes control Literature: [Khalil, ch.s 13, 14.1,14.] and [Glad-Ljung,ch.17] Lecture
More informationEL1820 Modeling of Dynamical Systems
EL1820 Modeling of Dynamical Systems Lecture 10 - System identification as a model building tool Experiment design Examination and prefiltering of data Model structure selection Model validation Lecture
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 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 informationIterative Controller Tuning Using Bode s Integrals
Iterative Controller Tuning Using Bode s Integrals A. Karimi, D. Garcia and R. Longchamp Laboratoire d automatique, École Polytechnique Fédérale de Lausanne (EPFL), 05 Lausanne, Switzerland. email: alireza.karimi@epfl.ch
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 informationComputer Aided Control Design
Computer Aided Control Design Project-Lab 3 Automatic Control Basic Course, EL1000/EL1100/EL1120 Revised August 18, 2008 Modified version of laboration developed by Håkan Fortell and Svante Gunnarsson
More information2.3 Oscillation. The harmonic oscillator equation is the differential equation. d 2 y dt 2 r y (r > 0). Its solutions have the form
2. Oscillation So far, we have used differential equations to describe functions that grow or decay over time. The next most common behavior for a function is to oscillate, meaning that it increases and
More informationCHAPTER 4 FOURIER SERIES S A B A R I N A I S M A I L
CHAPTER 4 FOURIER SERIES 1 S A B A R I N A I S M A I L Outline Introduction of the Fourier series. The properties of the Fourier series. Symmetry consideration Application of the Fourier series to circuit
More informationAn Improved Relay Auto Tuning of PID Controllers for SOPTD Systems
Proceedings of the World Congress on Engineering and Computer Science 7 WCECS 7, October 4-6, 7, San Francisco, USA An Improved Relay Auto Tuning of PID Controllers for SOPTD Systems Sathe Vivek and M.
More informationContents lecture 6 2(17) Automatic Control III. Summary of lecture 5 (I/III) 3(17) Summary of lecture 5 (II/III) 4(17) H 2, H synthesis pros and cons:
Contents lecture 6 (7) Automatic Control III Lecture 6 Linearization and phase portraits. Summary of lecture 5 Thomas Schön Division of Systems and Control Department of Information Technology Uppsala
More informationExercises Automatic Control III 2015
Exercises Automatic Control III 205 Foreword This exercise manual is designed for the course "Automatic Control III", given by the Division of Systems and Control. The numbering of the chapters follows
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 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 informationFrequency domain analysis
Automatic Control 2 Frequency domain analysis Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
More informationFourier Series and Fourier Transforms
Fourier Series and Fourier Transforms EECS2 (6.082), MIT Fall 2006 Lectures 2 and 3 Fourier Series From your differential equations course, 18.03, you know Fourier s expression representing a T -periodic
More informationEL2520 Control Theory and Practice
So far EL2520 Control Theory and Practice r Fr wu u G w z n Lecture 5: Multivariable systems -Fy Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden SISO control revisited: Signal
More informationControl 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 informationExam in Systems Engineering/Process Control
Department of AUTOMATIC CONTROL Exam in Systems Engineering/Process Control 7-6- Points and grading All answers must include a clear motivation. Answers may be given in English or Swedish. The total number
More informationHow many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation?
How many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation? (A) 0 (B) 1 (C) 2 (D) more than 2 (E) it depends or don t know How many of
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 informationSinusoidal Forcing of a First-Order Process. / τ
Frequency Response Analysis Chapter 3 Sinusoidal Forcing of a First-Order Process For a first-order transfer function with gain K and time constant τ, the response to a general sinusoidal input, xt = A
More information( ) Frequency Response Analysis. Sinusoidal Forcing of a First-Order Process. Chapter 13. ( ) sin ω () (
1 Frequency Response Analysis Sinusoidal Forcing of a First-Order Process For a first-order transfer function with gain K and time constant τ, the response to a general sinusoidal input, xt = A tis: sin
More information221B Lecture Notes on Resonances in Classical Mechanics
1B Lecture Notes on Resonances in Classical Mechanics 1 Harmonic Oscillators Harmonic oscillators appear in many different contexts in classical mechanics. Examples include: spring, pendulum (with a small
More informationSolutions to Problems in Chapter 4
Solutions to Problems in Chapter 4 Problems with Solutions Problem 4. Fourier Series of the Output Voltage of an Ideal Full-Wave Diode Bridge Rectifier he nonlinear circuit in Figure 4. is a full-wave
More informationNotes on the Periodically Forced Harmonic Oscillator
Notes on the Periodically orced Harmonic Oscillator Warren Weckesser Math 38 - Differential Equations 1 The Periodically orced Harmonic Oscillator. By periodically forced harmonic oscillator, we mean the
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 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 informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 8: Response Characteristics Overview In this Lecture, you will learn: Characteristics of the Response Stability Real
More informationStability of Parameter Adaptation Algorithms. Big picture
ME5895, UConn, Fall 215 Prof. Xu Chen Big picture For ˆθ (k + 1) = ˆθ (k) + [correction term] we haven t talked about whether ˆθ(k) will converge to the true value θ if k. We haven t even talked about
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 informationEE 435. Lecture 28. Data Converters Linearity INL/DNL Spectral Performance
EE 435 Lecture 8 Data Converters Linearity INL/DNL Spectral Performance Performance Characterization of Data Converters Static characteristics Resolution Least Significant Bit (LSB) Offset and Gain Errors
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 2: Drawing Bode Plots, Part 2 Overview In this Lecture, you will learn: Simple Plots Real Zeros Real Poles Complex
More informationCMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation
CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University September 26, 2005 1 Sinusoids Sinusoids
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 informationUnstable Oscillations!
Unstable Oscillations X( t ) = [ A 0 + A( t ) ] sin( ω t + Φ 0 + Φ( t ) ) Amplitude modulation: A( t ) Phase modulation: Φ( t ) S(ω) S(ω) Special case: C(ω) Unstable oscillation has a broader periodogram
More informationModern Control Systems
Modern Control Systems Matthew M. Peet Illinois Institute of Technology Lecture 18: Linear Causal Time-Invariant Operators Operators L 2 and ˆL 2 space Because L 2 (, ) and ˆL 2 are isomorphic, so are
More informationControl Systems I. Lecture 5: Transfer Functions. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 5: Transfer Functions Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 20, 2017 E. Frazzoli (ETH) Lecture 5: Control Systems I 20/10/2017
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 informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationTHE FOURIER TRANSFORM (Fourier series for a function whose period is very, very long) Reading: Main 11.3
THE FOURIER TRANSFORM (Fourier series for a function whose period is very, very long) Reading: Main 11.3 Any periodic function f(t) can be written as a Fourier Series a 0 2 + a n cos( nωt) + b n sin n
More informationBIBO 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 informationStability and Robustness 1
Lecture 2 Stability and Robustness This lecture discusses the role of stability in feedback design. The emphasis is notonyes/notestsforstability,butratheronhowtomeasurethedistanceto instability. The small
More informationThis is the number of cycles per unit time, and its units are, for example,
16 4. Sinusoidal solutions Many things in nature are periodic, even sinusoidal. We will begin by reviewing terms surrounding periodic functions. If an LTI system is fed a periodic input signal, we have
More informationTheoretical physics. Deterministic chaos in classical physics. Martin Scholtz
Theoretical physics Deterministic chaos in classical physics Martin Scholtz scholtzzz@gmail.com Fundamental physical theories and role of classical mechanics. Intuitive characteristics of chaos. Newton
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 information1) SIMPLE HARMONIC MOTION/OSCILLATIONS
1) SIMPLE HARMONIC MOTION/OSCILLATIONS 1.1) OSCILLATIONS Introduction: - An event or motion that repeats itself at regular intervals is said to be periodic. Periodicity in Space is the regular appearance
More informationSSY155 Applied Mechatronics Examination date
Chalmers University of Technology Department of Signals and Systems SSY155 Applied Mechatronics Examination date 937 Time: 8:3-12:3 Teacher: Jonas Sjöberg, tel 31-772 1855. Allowed material during the
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 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 informationSinusoids. Amplitude and Magnitude. Phase and Period. CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation
Sinusoids CMPT 889: Lecture Sinusoids, Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University September 6, 005 Sinusoids are
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 8: Response Characteristics Overview In this Lecture, you will learn: Characteristics of the Response Stability Real Poles
More informationPeriodogram of a sinusoid + spike Single high value is sum of cosine curves all in phase at time t 0 :
Periodogram of a sinusoid + spike Single high value is sum of cosine curves all in phase at time t 0 : X(t) = µ + Asin(ω 0 t)+ Δ δ ( t t 0 ) ±σ N =100 Δ =100 χ ( ω ) Raises the amplitude uniformly at all
More informationMAS107 Control Theory Exam Solutions 2008
MAS07 CONTROL THEORY. HOVLAND: EXAM SOLUTION 2008 MAS07 Control Theory Exam Solutions 2008 Geir Hovland, Mechatronics Group, Grimstad, Norway June 30, 2008 C. Repeat question B, but plot the phase curve
More informationFourier Transform. Find the Fourier series for a periodic waveform Determine the output of a filter when the input is a periodic function
Objectives: Be able to Fourier Transform Find the Fourier series for a periodic waveform Determine the output of a filter when the input is a periodic function Filters with a Single Sinusoidal Input: Suppose
More informationNonlinear Control. Nonlinear Control Lecture # 18 Stability of Feedback Systems
Nonlinear Control Lecture # 18 Stability of Feedback Systems Absolute Stability + r = 0 u y G(s) ψ( ) Definition 7.1 The system is absolutely stable if the origin is globally uniformly asymptotically stable
More informationEE3CL4: Introduction to Linear Control Systems
1 / 17 EE3CL4: Introduction to Linear Control Systems Section 7: McMaster University Winter 2018 2 / 17 Outline 1 4 / 17 Cascade compensation Throughout this lecture we consider the case of H(s) = 1. We
More informationLecture XXVI. Morris Swartz Dept. of Physics and Astronomy Johns Hopkins University November 5, 2003
Lecture XXVI Morris Swartz Dept. of Physics and Astronomy Johns Hopins University morris@jhu.edu November 5, 2003 Lecture XXVI: Oscillations Oscillations are periodic motions. There are many examples of
More informationAn anti-windup control of manufacturing lines: performance analysis
An anti-windup control of manufacturing lines: performance analysis A. Pogromsky Systems Engineering Group Eindhoven University of Technology The Netherlands In cooperation with R. A. van den Berg, W.
More informationContents lecture 5. Automatic Control III. Summary of lecture 4 (II/II) Summary of lecture 4 (I/II) u y F r. Lecture 5 H 2 and H loop shaping
Contents lecture 5 Automatic Control III Lecture 5 H 2 and H loop shaping Thomas Schön Division of Systems and Control Department of Information Technology Uppsala University. Email: thomas.schon@it.uu.se,
More informationLecture Note #6 (Chap.10)
System Modeling and Identification Lecture Note #6 (Chap.) CBE 7 Korea University Prof. Dae Ryook Yang Chap. Model Approximation Model approximation Simplification, approximation and order reduction of
More informationProcess Control & Instrumentation (CH 3040)
First-order systems Process Control & Instrumentation (CH 3040) Arun K. Tangirala Department of Chemical Engineering, IIT Madras January - April 010 Lectures: Mon, Tue, Wed, Fri Extra class: Thu A first-order
More informationAutomatic Control (MSc in Mechanical Engineering) Lecturer: Andrea Zanchettin Date: Student ID number... Signature...
Automatic Control (MSc in Mechanical Engineering) Lecturer: Andrea Zanchettin Date: 29..23 Given and family names......................solutions...................... Student ID number..........................
More informationIMC based automatic tuning method for PID controllers in a Smith predictor configuration
Computers and Chemical Engineering 28 (2004) 281 290 IMC based automatic tuning method for PID controllers in a Smith predictor configuration Ibrahim Kaya Department of Electrical and Electronics Engineering,
More informationCircuit Analysis-III. Circuit Analysis-II Lecture # 3 Friday 06 th April, 18
Circuit Analysis-III Sinusoids Example #1 ü Find the amplitude, phase, period and frequency of the sinusoid: v (t ) =12cos(50t +10 ) Signal Conversion ü From sine to cosine and vice versa. ü sin (A ± B)
More informationSignals and systems Lecture (S3) Square Wave Example, Signal Power and Properties of Fourier Series March 18, 2008
Signals and systems Lecture (S3) Square Wave Example, Signal Power and Properties of Fourier Series March 18, 2008 Today s Topics 1. Derivation of a Fourier series representation of a square wave signal
More informationExam in Systems Engineering/Process Control
Department of AUTOMATIC CONTROL Exam in Systems Engineering/Process Control 27-6-2 Points and grading All answers must include a clear motivation. Answers may be given in English or Swedish. The total
More informationExercises in Nonlinear Control Systems
Exercises in Nonlinear Control Systems Mikael Johansson Bo Bernhardsson Karl Henrik Johansson Created: 998 Latest update: November 9, 26 Introduction The exercises are divided into problem areas that roughly
More informationMATH4406 (Control Theory) Unit 1: Introduction Prepared by Yoni Nazarathy, July 21, 2012
MATH4406 (Control Theory) Unit 1: Introduction Prepared by Yoni Nazarathy, July 21, 2012 Unit Outline Introduction to the course: Course goals, assessment, etc... What is Control Theory A bit of jargon,
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 informationFinal Exam December 20, 2011
Final Exam December 20, 2011 Math 420 - Ordinary Differential Equations No credit will be given for answers without mathematical or logical justification. Simplify answers as much as possible. Leave solutions
More informationA conjecture on sustained oscillations for a closed-loop heat equation
A conjecture on sustained oscillations for a closed-loop heat equation C.I. Byrnes, D.S. Gilliam Abstract The conjecture in this paper represents an initial step aimed toward understanding and shaping
More informationLecture 5 Classical Control Overview III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore
Lecture 5 Classical Control Overview III Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore A Fundamental Problem in Control Systems Poles of open
More informationExam in Automatic Control II Reglerteknik II 5hp (1RT495)
Exam in Automatic Control II Reglerteknik II 5hp (1RT495) Date: August 4, 018 Venue: Bergsbrunnagatan 15 sal Responsible teacher: Hans Rosth. Aiding material: Calculator, mathematical handbooks, textbooks
More informationNON-LINEAR SYSTEMS. SUBJECT CODE: EIPE204 M.Tech, Second Semester. Prepared By. Dr. Kanhu Charan Bhuyan
NON-LINEAR SYSTEMS SUBJECT CODE: EIPE24 M.Tech, Second Semester Prepared By Dr. Kanhu Charan Bhuyan Asst. Professor Instrumentation and Electronics Engineering COLLEGE OF ENGINEERING AND TECHNOLOGY BHUBANESWAR
More informationExam. 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 informationAdvanced Aerospace Control. Marco Lovera Dipartimento di Scienze e Tecnologie Aerospaziali, Politecnico di Milano
Advanced Aerospace Control Dipartimento di Scienze e Tecnologie Aerospaziali, Politecnico di Milano ICT for control systems engineering School of Industrial and Information Engineering Aeronautical Engineering
More informationBioControl - Week 2, Lecture 1
BioControl - Week 2, Lecture Goals of this lecture: CDS tools for bio-molecular systems Study of equilibria Dynamic performance Periodic behaviors Main eample: Negative autoregulation Suggested readings
More informationLecture V : Oscillatory motion and spectral analysis
Lecture V : Oscillatory motion and spectral analysis I. IDEAL PENDULUM AND STABILITY ANALYSIS Let us remind ourselves of the equation of motion for the pendulum. Remembering that the external torque applied
More informationMAE143a: Signals & Systems (& Control) Final Exam (2011) solutions
MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions Question 1. SIGNALS: Design of a noise-cancelling headphone system. 1a. Based on the low-pass filter given, design a high-pass filter,
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 informationNyquist Criterion For Stability of Closed Loop System
Nyquist Criterion For Stability of Closed Loop System Prof. N. Puri ECE Department, Rutgers University Nyquist Theorem Given a closed loop system: r(t) + KG(s) = K N(s) c(t) H(s) = KG(s) +KG(s) = KN(s)
More information8.2 Harmonic Regression and the Periodogram
Chapter 8 Spectral Methods 8.1 Introduction Spectral methods are based on thining of a time series as a superposition of sinusoidal fluctuations of various frequencies the analogue for a random process
More informationr + - FINAL June 12, 2012 MAE 143B Linear Control Prof. M. Krstic
MAE 43B Linear Control Prof. M. Krstic FINAL June, One sheet of hand-written notes (two pages). Present your reasoning and calculations clearly. Inconsistent etchings will not be graded. Write answers
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 informationEngineering Tripos Part IIB Nonlinear Systems and Control. Handout 4: Circle and Popov Criteria
Engineering Tripos Part IIB Module 4F2 Nonlinear Systems and Control Handout 4: Circle and Popov Criteria 1 Introduction The stability criteria discussed in these notes are reminiscent of the Nyquist criterion
More informationLecture 5: Frequency domain analysis: Nyquist, Bode Diagrams, second order systems, system types
Lecture 5: Frequency domain analysis: Nyquist, Bode Diagrams, second order systems, system types Venkata Sonti Department of Mechanical Engineering Indian Institute of Science Bangalore, India, 562 This
More informationChapter 6 - Solved Problems
Chapter 6 - Solved Problems Solved Problem 6.. Contributed by - James Welsh, University of Newcastle, Australia. Find suitable values for the PID parameters using the Z-N tuning strategy for the nominal
More information100 (s + 10) (s + 100) e 0.5s. s 100 (s + 10) (s + 100). G(s) =
1 AME 3315; Spring 215; Midterm 2 Review (not graded) Problems: 9.3 9.8 9.9 9.12 except parts 5 and 6. 9.13 except parts 4 and 5 9.28 9.34 You are given the transfer function: G(s) = 1) Plot the bode plot
More information