# Design and Tuning of Fractional-order PID Controllers for Time-delayed Processes

Size: px
Start display at page:

Transcription

1 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 Reza Katebi Industrial Control Centre University of Strathclyde Glasgow, United Kingdom Abstract Frequency domain based design methods are investigated for the design and tuning of fractional-order PID for scalar applications. Since Ziegler-Nichol s tuning rule and other algorithms cannot be applied directly to tuning of fractional-order controllers, a new algorithm is developed to handle the tuning of these fractional-order PID controllers based on a single frequency point just like Ziegler-Nichol s rule for inter order PID. Critical parameters of the system are obtained at the ultimate point and the controller parameters are calculated from these critical measurements to meet design specifications. Thereafter, fractional order is obtained to meet a specified robustness criteria which is the phase-invariability against gain variations around the phase cross-over frequency. Results are simulated on second order plus dead time plant to demonstrate both performance and robustness. Keywords robustness; performance; tuning; PID I. INTRODUCTION In recent years, interesting control engineering applications have been developed using differential equations of generalized order α where the derivative order α can be any non-integer. Control systems utilizing fractional dynamics termed fractional order control systems have been found to yield better performance than integer-order controllers under fair comparison []. A. Some Definitions of Fractional-order Derivative Consider the generalised-order differential equation () shown below: y(t) = D α u(t) + u(t) () where t > ; n < α < n; n R Equation is a fractional order differential equation as far as the order "α" is non-integer. There are two important definitions that allow for computation of fractional order derivative namely: Caputo s Fractional order derivative Riemann-Liouville Definition. Fractional order derivative with order (α) is defined by Riemann-Liouville as: D α f(t) D n l n α f(t) = dn t dt n [ Γ(n α) f(τ) (t τ) α n+ ] dτ. () where t > ; n < α < n; n R. An alternative definition was also given by Caputo [] as shown in equation 3: D α f(t) l n α D n f(t) = [ Γ(n α) t > ; n < α < n; n R t f n (τ) (t τ) α n+ dτ] (3) However, the inclusion of the nth order derivative of f(t) in Caputo s definition imposes restriction compared to Riemann- Liouville form [3]. B. Fractional-order PID The theory of fractional dynamics has been utilized to design PIDs with non-integer order - PI D. Fractional order PID Controllers (PI D ) have been extensively tested in demanding applications especially in mechatronic and automatic control applications and are found to yield very good results []. Basically, PI D is of the form: C(s) = K p + K i s λ + K ds μ (4) and are the fractional orders of integral and derivative parts of the controller respectively. Therefore there are five parameters that are to be determined: K P, K I, K D,, and. It can be observed that classical integer PID controller is a subset of this generalized-order PI D where and are set to one. In terms of frequency domain analysis, given the nonrestriction of the orders of integral and derivative terms of the controller, more parameters are available for the manipulation of the slope of the magnitude curve and phase at both high and low frequencies. In terms of robustness to plant gain variations and uncertainties in the plant model, the fractional order PID controller is known to give better performance than the classical PID controller due to this inherent frequency domain property as well as availability of more degrees of freedom for the design engineer []. Also the high frequency noises, and load disturbances associated with plant models can be effectively handled using fractional controller [4].

2 II. SPECIFICATION OF DESIGN OBJECTIVES The first step is to define design specifications expected to be satisfied by the fractional order controller. These design objectives are defined in frequency domain in order to take care of important control objectives: stability, performance and robustness. Each controller parameter can be tuned to satisfy each design specification. A. Gain Margin and Phase Cross-over Frequency Specification. Gain margin is a primary index of relative stability in classical control design. A pre-defined margin can be used to formulate a robustness constraint on the system gain. Equation I defines the relationship between gain margin and phase cross over frequency while equation 5 defines the gain margin constraint: C(jω cg )G(jw cg ) db = db (5) C(jω cp )G(jw cp ) db = ω cg - The gain crossover frequency w cp Phase cross over frequency A m The Gain margin. A m (6) B. Phase Margin and Gain Crossover Frequency Specification Phase margin and gain margin are extended to fractional order control as useful measures of robust stability. It is chosen to satisfy equation 7 arg (C(jω cg )G(jw cg )) = π + m (7) C. Robustness against Plant s Gain Variations Bode s ideal loop defines the criteria for absolutely stable SISO loop. As long as the open loop gain is defined by a constant phase around the crossover frequency (useful band), robust stability against gain variations is guaranteed within that frequency range of constant phase. d{arg(c(jω)g(jω))} = around ω dω cg. (8) The phase of the forward loop function will be flat around the cross over frequency. This improves robustness against gainlike variations in plant and the overshoot is nearly constant within that given frequency range. D. High frequency noise rejection specification In order to ensure satisfactory measurement noise rejection, appropriate bound for the complementary sensitivity function has to be defined: T(jω) db A db C(jω)G(jω) A db for ω ω +C(jω)G(jω) t rad/s. (9) db T(jω) is the complementary sensitivity function. A db = the specified attenuation level in db for the band ω ω t rad s. E. Output disturbance rejection specification. Sensitivity bound can be defined to ensure satisfactory output disturbance rejection by the controller: +C(jω)G(jω) db B db () B db = the specified magnitude of the sensitivity function in db for the band ω ω s rad s. In addition, the nth-fractional integral term ensures steady state error rejection as effectively as an equivalent integer-order integral controller due to the implementation as s s n. Five parameters are available for tuning in order to meet five specifications. The proportional gain, integral gain, derivative gain, derivative (fractional) order and the fractional order of the integral term are available parameters that can be optimised to meet the defined specifications [3]. This is the primary reason behind the superior control action of fractional order controller type over the integer order types more parameters are available for optimization to meet more control objectives. Monje [3, 5] used MATLAB s optimization tool (fmincon) to formulate an optimization problem that solves for these five parameters and it was a constraint non-linear optimization problem with five unknown variables. Another way of obtaining robust values for PID controller gains is to use sustained oscillation method as described next while the fractional order for the integral controller term can be chosen to satisfy robust performance criteria. This method reduces computational cost. III. THE ULTIMATE POINT PARAMETER METHOD In industrial circles, Ziegler-Nichols PID tuning rule and similar PID design methods are still dominant for simple process models. In this section, a similar technique is developed for SISO processes. The basic principle involves using a proportional controller in cascade with the plant under a closed loop configuration. Thereafter, the proportional gain is systematically increased from very small values until sustained oscillation or continuous cycling is observed at the process output. The value of the proportional gain that yields this sustained oscillation is recorded as the critical gain (K u ) while the period of oscillation is noted as the ultimate period T u. Several other important information about the process are obtainable from these critical point measurements. For instance, the phase cross over frequency can be obtained as: ω pco = π T u () The fractional PI controller is then obtained by relocating this ultimate frequency point on the complex plane to a more desirable point. This is illustrated on the Nyquist diagram of figure. The critical point or ultimate point is the stability limit signified by the point where the Nyquist curve intersects the unit circle (-,).

3 Imaginary Axis Critical Point (-,) FOPDT Plant at Critical Point Recall the controller structure in frequency domain as given in equation : C(jω) = K p + (jω) α C(jω) = K( + (τ i )(jω) α) C(jω) = K( + (cos απ απ j sin ) ) (8) ω α τ i K( + ω α (cos απ απ j sin τ i ) ) = r B K u (cos B + j sin B ). Comparing magnitudes: r B K u cos B = K ( + ω α (cos απ )) (9) τ i K i Fig.. Nyquist Diagram indicating Plant s Critical Point A. Design of the Fractional-order PI controller Consider a fractional PI controller of the form: C(s) = K p + K i () sα K i C(jω) = K p + (jω) α K i = K p T i A desirable point on the complex plane is chosen as design point - B given as (r B e j(π+ B) ). The controller is expected to move the ultimate point ( K u, ) to this desirable point. On the Nyquist plot, the ultimate point is - point A (r A e jπ) ). r A = K u (3) Finally, Let the frequency characteristic of the controller C(s) be r C e j(π+ C). For this controller to move the ultimate point to the desirable point on the complex plane, equation 4 must be satisfied: r A r C e jπ e j(π+ C) = r B e j(π+ B) (4). r C e j(π+ C) = r B r A (cos B + j sin B ) (5) r C e j(π+ C) = r B K u (cos B + j sin B ) Therefore, the fractional-order PI controller is chosen to satisfy the magnitude and phase conditions below in 6 and 7: r C = r B r A (6) C = B (7) In the same vein, arguments are compared and the phase is expected to fulfil equation below: tan B = If we set y = τ i ; y = ω α τ i sin απ/ ( + ω α τ i cos απ/) () tan B ω α (tan B cos απ + sin απ/) () B. Justification of Design Point The desired point B is chosen such that r B =.9 and B = 46 deg. This is a very desirable point because it translates to -.-j. point. The distance from this point to the stability limit (-,j) is d = (. + ) + (.) which is approximately.9. Closed loop sensitivity can be deduced from this information on the Nyquist curve of the forward loop because the maximum sensitivity M s is the reciprocal of the shortest distance from the Nyquist curve to the critical point. By definition: M s = max. ω +C(s)G(s) However, on the Nyquist curve, M s =. d M s =.9 =. This is very desirable for robustness. For the closed loop system to be robust against variations in process dynamics, maximum sensitivity of the closed loop system is specified: M s <. Reasonable values of M s ranges from.5 to.95. This indicates good robustness for many practical control applications and it is therefore recommended as the design point as far as single frequency point method is concerned throughout this paper. This design point is also justified as far as closed loop stability is concerned. It meets Nyquist stability criterion as the critical stability (-,) point will not be encircled. Therefore, the fractional order PI controller can be obtained from and 3. K =.K u ( τ u π )α cos απ ()

4 τ i =.966 ( τ u π )α (sin απ απ +.36cos ) (3) Fractional order is selected to meet the required resonant peak M r condition as given by Monje for SISO-FOPDT systems. This is sufficient for SISO control problems with FOPDT process model. G(s) = Ye Ls. Relative dead-time is: T = L sτ+ L+τ (4). TABLE I. RELATIVE DEAD TIME AND FRACTIONAL-ORDER Cases - order Relative Dead-time T..7 T <...9. T < T < T.6 For Second order processes such as SOPDT (Second Order Plus Dead-Time model), derivative component can be included for a more damped response. One possible solution is to design the fractional controller in the form of PI μ D-structure. It is however tuned in a similar fashion as the fractional PI case using ultimate parameters. C. The PI μ D - Structure The structure of the controller is given as: K i C(jω) = K p + (jω) α + K D(jω) (5) C(jω) = K p ( + τ i (cos απ απ + j sin ) + τ D(jω) (6) ωα Let the controller C(jω) be characterised by r C e j(π+ C) in frequency domain and the desirable point on the Nyquist plane chosen as point B that is: (r B e j(π+ B) ) as earlier given. The controller is expected to move the ultimate point (, ) K u to this desirable point as earlier explained in previous section. On the complex plane, the ultimate point is - point A (r A e jπ) ). Therefore, equations 4, 5 and 6 must be fulfilled as earlier shown. It follows that: (K p ( + τ i (cos απ απ + j sin ) + τ D(jω) ) ωα = r B cos( r B A ) (7) A Consider the argument: ωτ D ω tan( B A ) = α sin( απ/) τ i + ω α cos(απ/) τ i Given that the ultimate point is to be moved; A =. tan B ( + ω α cos ( απ τ i )) = ωτ D ω α sin ( απ/) τ i Making τ i subject of formular and substituting ω = π yields: τ i = τ u α (tan B cos απ απ + sin ) (π) α ( πτ (8) D tan τ B ) u K p = r B K u cos B τ u α cos απ ((π) α (9) τ i ) Given that the order α is obtainable as before from the M r table, there are now three parameters left to be obtained from two equations. This is clearly insufficient. One way to solve it is to treat the ratio of τ D to τ i as a constant [6]. For instance: τ D =.5τ i. All four parameters can therefore be calculated. The second structure of fractional controller with derivative component is of PID α form. D. The PID α Controller form. C(jω) = K p + K i jω + K D(jω) α (3) C(jω) = K p ( j τ i ω + τ D(jω) α ) The controller structure is given in equation 3. It is expected to move the ultimate point to our new design point B implying that equations 4, 5 and 6 applies again as stated in previous sections. Therefore: (K p ( + ω α τ D cos(απ/)) = r B cos( r B A ) A r B cos( ra B A ) K p = + ω α τ D cos απ K p = r BK u cos( B ) + ω α τ D cos απ (3) Considering the argument as before: ω α τ D sin( απ ) ωτ tan( B ) = i + ω α τ D cos απ (3) The ratio of integral time to derivative time is set to constant as before: τ D =.5τ i and equations 3 and 3 gives the two controller parameters. Fractional order are obtained in similar fashion to previous cases using the peak resonant condition. IV. RESULTS OF SISO PLANT SIMULATIONS Consider a FOPDT plant given as G(s) =.8e s 6.7s + Three controllers C, C and C 3 are designed in the following forms: PI μ, PI μ D and PID α respectively. This is done using the sustained oscillation method as described in this work: T u

5 Imaginary Axis Imaginary Axis Imaginary Axis FOPDT with C Gain Margin: 3. db Phase Margin: 44.7 deg Delay Margin: Fig.. Nyquist Diagram of FOPDT Plant with controller C C =.4( +.s.7) s.7 is implemented as s.3 s using the well-known ORA approximation (Oustaloup Recursive Approximation) inorder to obtain integer order form of the controller with the same frequency domain properties [7, 8]. ORA has been reported by many authors as a good approximation to the fractional order. Therefore, it is used throughout this work to realise a rational and implementable function. C =.3( s) s.7 C 3 =.8( + 3.6s +.8s.7 ) It is observed that the controller moves the design point of the plant in figure to a desired (stable) region as shown in Nyquist diagrams shown next: Consider a second example - a SOPDT plant below: FOPDT with C Gain Margin: 7. db Phase Margin: 69.7 deg Delay Margin: 4.5 e.s G (s) = s(.5s + ) The controller realized using PI μ D structure as earlier discussed is: C (s) = s. The fractional controller of the form PID α is also realized as C given below: C = s.9. s Performance is compared using Bode diagram of figure 5and step response shown in figure 6. s Fig. 3. Nyquist Diagram of FOPDT Plant with controller C FOPDT with C3 Magnitude (db) C GAIN MARGIN:. C PHASE MARGIN: 7.3 C DELAY MARGIN:. Control of SOPDT plant c c Gain Margin: 3.3 db Phase Margin: 76.7 deg Delay Margin: 3.3 Phase (deg) C GAIN MARGIN: 7.4 C PHASE MARGIN: 53.4 C DELAY MARGIN: Frequency (rad/s) Fig. 4. Nyquist Diagram of FOPDT Plant with controller C3 Fig. 5. Bode Diagram of SOPDT Plant with controller C and C

6 V. DISCUSSION OF RESULTS AND CONCLUSION These SISO based methods yield excellent response for the plants examined as reflected in the gain margin, phase margin and delay margin. The delay margin is an important robustness measure which gives information about how much time delay can be introduced to the system before it goes unstable. In addition, very little computation is required to yield the desired fractional controller.in summary, the main contribution of this paper is development of a new design method for fractional-order PID (FPID) controller. Extended simulations has been carried out on delayed processes with FOPDT and SOPDT dynamics. It has been found to yield control actions with good compromise between robustness and performance..5.5 Step Response -Fractional Control of SOPDT Plant P-Plus-Fractional-I-plus-D P-Plus-I-plus-Fractional-D Fig. 6. Diagram of step response SOPDT Plant with controller C and C VI. REFERENCES [] Z. Li and Y. Chen, "Ideal, Simplified and Inverted Decoupling of Fractional order TITO Processes," in IFAC, Cape Town, 4. [] I. Podlubny, "Fractional-Order Sysmtems and PID Controllers," IEEE Transactions on Automatic Control, vol. 44, no., pp. 8-4, 999. [3] C. Monje, Y. Chen, B. D. Vinagre and V. Feliu, Fractional-order Systems and Controls: Fundamentals and Applications, London: Springer-Verlag,. [4] G. Vinagre, D. Valerio and J. Costa, "Rule-Tuned PIDs and Fractional PIDs for a Three-Tank Liquid System," in Symposium on Fractional Signals and Systems, Lisbon, Portugal, 9. [5] C. Monje, Y. Chen, D. Vinagre and V. Feliu, "Tuning and Auto-tuning of Fractional order Controllers for Industry Applications," Control Engineering Practice, vol. 6, pp , 8. [6] K. Astrom and T. Hagglund, Advanced PID Control, Durham, USA.: ISA, 6. [7] D. Valerio and S. Costa, "Time domain implementation of Fractional order controllers," IEEE Transaction on Control Theory Application, vol. 5, no. 5, pp , 5. [8] F. Merrikh-Bayat, "Rules for Selecting the Parameters of Oustaloup Recurssive Approximation for the Simulation of Linear Feedback Systems containing Fractional PID controller," Nonlinear Science and Numerical Simulation, vol. 7, pp ,.

### TUNING-RULES FOR FRACTIONAL PID CONTROLLERS

TUNING-RULES FOR FRACTIONAL PID CONTROLLERS Duarte Valério, José Sá da Costa Technical Univ. of Lisbon Instituto Superior Técnico Department of Mechanical Engineering GCAR Av. Rovisco Pais, 49- Lisboa,

### Robust PID and Fractional PI Controllers Tuning for General Plant Model

2 مجلة البصرة للعلوم الهندسية-المجلد 5 العدد 25 Robust PID and Fractional PI Controllers Tuning for General Plant Model Dr. Basil H. Jasim. Department of electrical Engineering University of Basrah College

### Frequency 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

### Tuning of fractional PI controllers for fractional order system models with and without time delays

2 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 3-July 2, 2 FrC2. Tuning of fractional PI controllers for fractional order system models with and without time delays Anuj Narang,

### Part II. Advanced PID Design Methods

Part II Advanced PID Design Methods 54 Controller transfer function C(s) = k p (1 + 1 T i s + T d s) (71) Many extensions known to the basic design methods introduced in RT I. Four advanced approaches

### ECE 388 Automatic Control

Lead Compensator and PID Control Associate Prof. Dr. of Mechatronics Engineeering Çankaya University Compulsory Course in Electronic and Communication Engineering Credits (2/2/3) Course Webpage: http://ece388.cankaya.edu.tr

### Control System Design

ELEC ENG 4CL4: Control System Design Notes for Lecture #11 Wednesday, January 28, 2004 Dr. Ian C. Bruce Room: CRL-229 Phone ext.: 26984 Email: ibruce@mail.ece.mcmaster.ca Relative Stability: Stability

### 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

### The loop shaping paradigm. Lecture 7. Loop analysis of feedback systems (2) Essential specifications (2)

Lecture 7. Loop analysis of feedback systems (2). Loop shaping 2. Performance limitations The loop shaping paradigm. Estimate performance and robustness of the feedback system from the loop transfer L(jω)

### Analysis 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

### Iterative 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

### CHAPTER 7 FRACTIONAL ORDER SYSTEMS WITH FRACTIONAL ORDER CONTROLLERS

9 CHAPTER 7 FRACTIONAL ORDER SYSTEMS WITH FRACTIONAL ORDER CONTROLLERS 7. FRACTIONAL ORDER SYSTEMS Fractional derivatives provide an excellent instrument for the description of memory and hereditary properties

### MULTILOOP PI CONTROLLER FOR ACHIEVING SIMULTANEOUS TIME AND FREQUENCY DOMAIN SPECIFICATIONS

Journal of Engineering Science and Technology Vol. 1, No. 8 (215) 113-1115 School of Engineering, Taylor s University MULTILOOP PI CONTROLLER FOR ACHIEVING SIMULTANEOUS TIME AND FREQUENCY DOMAIN SPECIFICATIONS

### IMC 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,

### Chapter 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

### Lecture 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

### 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

### Cascade Control of a Continuous Stirred Tank Reactor (CSTR)

Journal of Applied and Industrial Sciences, 213, 1 (4): 16-23, ISSN: 2328-4595 (PRINT), ISSN: 2328-469 (ONLINE) Research Article Cascade Control of a Continuous Stirred Tank Reactor (CSTR) 16 A. O. Ahmed

### Robust Loop Shaping Controller Design for Spectral Models by Quadratic Programming

Robust Loop Shaping Controller Design for Spectral Models by Quadratic Programming Gorka Galdos, Alireza Karimi and Roland Longchamp Abstract A quadratic programming approach is proposed to tune fixed-order

### Chapter 2. Classical Control System Design. Dutch Institute of Systems and Control

Chapter 2 Classical Control System Design Overview Ch. 2. 2. Classical control system design Introduction Introduction Steady-state Steady-state errors errors Type Type k k systems systems Integral Integral

### 1 Loop Control. 1.1 Open-loop. ISS0065 Control Instrumentation

Lecture 4 ISS0065 Control Instrumentation 1 Loop Control System has a continuous signal (analog) basic notions: open-loop control, close-loop control. 1.1 Open-loop Open-loop / avatud süsteem / открытая

### Prü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

### Control Systems I Lecture 10: System Specifications

Control Systems I Lecture 10: System Specifications Readings: Guzzella, Chapter 10 Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 24, 2017 E. Frazzoli (ETH) Lecture

### CHBE320 LECTURE XI CONTROLLER DESIGN AND PID CONTOLLER TUNING. Professor Dae Ryook Yang

CHBE320 LECTURE XI CONTROLLER DESIGN AND PID CONTOLLER TUNING Professor Dae Ryook Yang Spring 2018 Dept. of Chemical and Biological Engineering 11-1 Road Map of the Lecture XI Controller Design and PID

### Classify a transfer function to see which order or ramp it can follow and with which expected error.

Dr. J. Tani, Prof. Dr. E. Frazzoli 5-059-00 Control Systems I (Autumn 208) Exercise Set 0 Topic: Specifications for Feedback Systems Discussion: 30.. 208 Learning objectives: The student can grizzi@ethz.ch,

### Design Methods for Control Systems

Design Methods for Control Systems Maarten Steinbuch TU/e Gjerrit Meinsma UT Dutch Institute of Systems and Control Winter term 2002-2003 Schedule November 25 MSt December 2 MSt Homework # 1 December 9

### CDS 101/110a: Lecture 8-1 Frequency Domain Design

CDS 11/11a: Lecture 8-1 Frequency Domain Design Richard M. Murray 17 November 28 Goals: Describe canonical control design problem and standard performance measures Show how to use loop shaping to achieve

### Closed Loop Identification Of A First Order Plus Dead Time Process Model Under PI Control

Dublin Institute of Technology RROW@DIT Conference papers School of Electrical and Electronic Engineering -6- Closed Loop Identification Of First Order Plus Dead Time Process Model Under PI Control Tony

### Lecture 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

### 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

### Controls 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

### Raktim 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

### 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

### Robust 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

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

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

### THE ball and beam system is one of the most popular

Vol:8, No:7, 24 Fractional Order Feedback Control of a Ball and Beam System Santosh Kr. Choudhary International Science Index, Electrical and Computer Engineering Vol:8, No:7, 24 waset.org/publication/99999

### Linear 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

### Outline. Classical Control. Lecture 1

Outline Outline Outline 1 Introduction 2 Prerequisites Block diagram for system modeling Modeling Mechanical Electrical Outline Introduction Background Basic Systems Models/Transfers functions 1 Introduction

### 7.2 Controller tuning from specified characteristic polynomial

192 Finn Haugen: PID Control 7.2 Controller tuning from specified characteristic polynomial 7.2.1 Introduction The subsequent sections explain controller tuning based on specifications of the characteristic

### DESIGN 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)

### Robust fixed-order H Controller Design for Spectral Models by Convex Optimization

Robust fixed-order H Controller Design for Spectral Models by Convex Optimization Alireza Karimi, Gorka Galdos and Roland Longchamp Abstract A new approach for robust fixed-order H controller design by

### Ü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 -

### MEM 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

### Digital 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

### Frequency 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

### K c < K u K c = K u K c > K u step 4 Calculate and implement PID parameters using the the Ziegler-Nichols tuning tables: 30

1.5 QUANTITIVE PID TUNING METHODS Tuning PID parameters is not a trivial task in general. Various tuning methods have been proposed for dierent model descriptions and performance criteria. 1.5.1 CONTINUOUS

### MAS107 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

### Robust QFT-based PI controller for a feedforward control scheme

Integral-Derivative Control, Ghent, Belgium, May 9-11, 218 ThAT4.4 Robust QFT-based PI controller for a feedforward control scheme Ángeles Hoyo José Carlos Moreno José Luis Guzmán Tore Hägglund Dep. of

### A NEW APPROACH TO MIXED H 2 /H OPTIMAL PI/PID CONTROLLER DESIGN

Copyright 2002 IFAC 15th Triennial World Congress, Barcelona, Spain A NEW APPROACH TO MIXED H 2 /H OPTIMAL PI/PID CONTROLLER DESIGN Chyi Hwang,1 Chun-Yen Hsiao Department of Chemical Engineering National

### Model-based PID tuning for high-order processes: when to approximate

Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 25 Seville, Spain, December 2-5, 25 ThB5. Model-based PID tuning for high-order processes: when to approximate

### Ü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. 8. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid

### Robust and Optimal Control, Spring A: SISO Feedback Control A.1 Internal Stability and Youla Parameterization

Robust and Optimal Control, Spring 2015 Instructor: Prof. Masayuki Fujita (S5-303B) A: SISO Feedback Control A.1 Internal Stability and Youla Parameterization A.2 Sensitivity and Feedback Performance A.3

### Control System Design

ELEC ENG 4CL4: Control System Design Notes for Lecture #15 Friday, February 6, 2004 Dr. Ian C. Bruce Room: CRL-229 Phone ext.: 26984 Email: ibruce@mail.ece.mcmaster.ca (3) Cohen-Coon Reaction Curve Method

### Control 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

### 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.

### MEM 355 Performance Enhancement of Dynamical Systems

MEM 355 Performance Enhancement of Dynamical Systems Frequency Domain Design Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 5/8/25 Outline Closed Loop Transfer Functions

### An 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.

### ECE 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

### Improved Identification and Control of 2-by-2 MIMO System using Relay Feedback

CEAI, Vol.17, No.4 pp. 23-32, 2015 Printed in Romania Improved Identification and Control of 2-by-2 MIMO System using Relay Feedback D.Kalpana, T.Thyagarajan, R.Thenral Department of Instrumentation Engineering,

### Intro to Frequency Domain Design

Intro to Frequency Domain Design MEM 355 Performance Enhancement of Dynamical Systems Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Closed Loop Transfer Functions

### PID control of FOPDT plants with dominant dead time based on the modulus optimum criterion

Archives of Control Sciences Volume 6LXII, 016 No. 1, pages 5 17 PID control of FOPDT plants with dominant dead time based on the modulus optimum criterion JAN CVEJN The modulus optimum MO criterion can

### Lecture 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

### EE3CL4: 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

### Ü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 -

### CDS 101/110a: Lecture 10-1 Robust Performance

CDS 11/11a: Lecture 1-1 Robust Performance Richard M. Murray 1 December 28 Goals: Describe how to represent uncertainty in process dynamics Describe how to analyze a system in the presence of uncertainty

### Improved Autotuning Using the Shape Factor from Relay Feedback

Article Subscriber access provided by NATIONAL TAIWAN UNIV Improved Autotuning Using the Shape Factor from Relay Feedback T. Thyagarajan, and Cheng-Ching Yu Ind. Eng. Chem. Res., 2003, 42 (20), 4425-4440

### VALLIAMMAI 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

### Process Identification for an SOPDT Model Using Rectangular Pulse Input

Korean J. Chem. Eng., 18(5), 586-592 (2001) SHORT COMMUNICATION Process Identification for an SOPDT Model Using Rectangular Pulse Input Don Jang, Young Han Kim* and Kyu Suk Hwang Dept. of Chem. Eng., Pusan

### Experimental Study of Fractional Order Proportional Integral (FOPI) Controller for Water Level Control

Proceedings of the 47th IEEE Conference on Decision and Control Cancun, Mexico, Dec. 9-11, 2008 Experimental Study of Fractional Order Proportional Integral (FOPI) Controller for Water Level Control Varsha

### A unified approach for proportional-integral-derivative controller design for time delay processes

Korean J. Chem. Eng., 32(4), 583-596 (2015) DOI: 10.1007/s11814-014-0237-6 INVITED REVIEW PAPER INVITED REVIEW PAPER pissn: 0256-1115 eissn: 1975-7220 A unified approach for proportional-integral-derivative

### 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

### 100 (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

### (a) Find the transfer function of the amplifier. Ans.: G(s) =

126 INTRDUCTIN T CNTR ENGINEERING 10( s 1) (a) Find the transfer function of the amplifier. Ans.: (. 02s 1)(. 001s 1) (b) Find the expected percent overshoot for a step input for the closed-loop system

### Analysis and Synthesis of Single-Input Single-Output Control Systems

Lino Guzzella Analysis and Synthesis of Single-Input Single-Output Control Systems l+kja» \Uja>)W2(ja»\ um Contents 1 Definitions and Problem Formulations 1 1.1 Introduction 1 1.2 Definitions 1 1.2.1 Systems

### 9. Two-Degrees-of-Freedom Design

9. Two-Degrees-of-Freedom Design In some feedback schemes we have additional degrees-offreedom outside the feedback path. For example, feed forwarding known disturbance signals or reference signals. In

### 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

### Frequency 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

### Iterative Feedback Tuning

Iterative Feedback Tuning Michel Gevers CESAME - UCL Louvain-la-Neuve Belgium Collaboration : H. Hjalmarsson, S. Gunnarsson, O. Lequin, E. Bosmans, L. Triest, M. Mossberg Outline Problem formulation Iterative

### Tuning Rules for Proportional Resonant Controllers

Tuning Rules for Proportional Resonant Controllers Luís Fernando Alves Pereira and Alexandre Sanfelice Bazanella, Senior Member, IEEE Abstract In this paper we propose a particular structure for resonant

### Ch 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

### Control of integral processes with dead time Part IV: various issues about PI controllers

Control of integral processes with dead time Part IV: various issues about PI controllers B. Wang, D. Rees and Q.-C. Zhong Abstract: Various issues about integral processes with dead time controlled by

### 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

### 6.1 Sketch the z-domain root locus and find the critical gain for the following systems K., the closed-loop characteristic equation is K + z 0.

6. Sketch the z-domain root locus and find the critical gain for the following systems K (i) Gz () z 4. (ii) Gz K () ( z+ 9. )( z 9. ) (iii) Gz () Kz ( z. )( z ) (iv) Gz () Kz ( + 9. ) ( z. )( z 8. ) (i)

### Auto-tuning Fractional Order Control of a Laboratory Scale Equipment

Auto-tuning Fractional Order Control of a Laboratory Scale Equipment Rusu-Both Roxana * and Dulf Eva-Henrietta Automation Department, Technical University of Cluj-Napoca, Memorandumului Street, No.28 Cluj-Napoca,

### EE3CL4: Introduction to Linear Control Systems

1 / 30 EE3CL4: Introduction to Linear Control Systems Section 9: of and using Techniques McMaster University Winter 2017 2 / 30 Outline 1 2 3 4 / 30 domain analysis Analyze closed loop using open loop

### Introduction. Performance and Robustness (Chapter 1) Advanced Control Systems Spring / 31

Introduction Classical Control Robust Control u(t) y(t) G u(t) G + y(t) G : nominal model G = G + : plant uncertainty Uncertainty sources : Structured : parametric uncertainty, multimodel uncertainty Unstructured

### University of Science and Technology, Sudan Department of Chemical Engineering.

ISO 91:28 Certified Volume 3, Issue 6, November 214 Design and Decoupling of Control System for a Continuous Stirred Tank Reactor (CSTR) Georgeous, N.B *1 and Gasmalseed, G.A, Abdalla, B.K (1-2) University

### FREQUENCY-RESPONSE DESIGN

ECE45/55: Feedback Control Systems. 9 FREQUENCY-RESPONSE DESIGN 9.: PD and lead compensation networks The frequency-response methods we have seen so far largely tell us about stability and stability margins

### 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 7 Root Locus 2 Assign

### PID CONTROLLER SYNTHESIS FREE OF ANALYTICAL MODELS

PID CONTROLLER SYNTHESIS FREE OF ANALYTICAL MODELS L.H. Keel,1 S.P. Bhattacharyya, Tennessee State University, Nashville, Tennessee, USA Texas A&M University, College Station, Texas, USA Abstract: In this

### Exercises for lectures 13 Design using frequency methods

Exercises for lectures 13 Design using frequency methods Michael Šebek Automatic control 2016 31-3-17 Setting of the closed loop bandwidth At the transition frequency in the open loop is (from definition)

### CHAPTER 6 CLOSED LOOP STUDIES

180 CHAPTER 6 CLOSED LOOP STUDIES Improvement of closed-loop performance needs proper tuning of controller parameters that requires process model structure and the estimation of respective parameters which

### Generalized Fractional Order Reset Element (GFrORE)

ENOC 217, June 25-3, 217, Budapest, Hungary Generalized Fractional Order Reset Element (GFrORE) Niranjan Saikumar and Hassan HosseinNia Precision and Microsystems Engineering, 3ME, TU Delft, The Netherlands

### MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall 2007

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering.4 Dynamics and Control II Fall 7 Problem Set #9 Solution Posted: Sunday, Dec., 7. The.4 Tower system. The system parameters are

### MAE 143B - Homework 9

MAE 43B - Homework 9 7.2 2 2 3.8.6.4.2.2 9 8 2 2 3 a) G(s) = (s+)(s+).4.6.8.2.2.4.6.8. Polar plot; red for negative ; no encirclements of, a.s. under unit feedback... 2 2 3. 4 9 2 2 3 h) G(s) = s+ s(s+)..2.4.6.8.2.4

### 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

### Robust and Gain-Scheduled PID Controller Design for Condensing Boilers by Linear Programming

Robust and Gain-Scheduled PID Controller Design for Condensing Boilers by Linear Programming Vinicius de Oliveira and Alireza Karimi Laboratoire d Automatque École Polytechnique Fédérale de Lausanne (EPFL)

### CHAPTER 10: STABILITY &TUNING

When I complete this chapter, I want to be able to do the following. Determine the stability of a process without control Determine the stability of a closed-loop feedback control system Use these approaches