A MATLAB Toolbox for Teaching Model Order Reduction Techniques

Size: px
Start display at page:

Download "A MATLAB Toolbox for Teaching Model Order Reduction Techniques"

Transcription

1 A MATLAB Toolbox for Teaching Model Order Reduction Techniques Authors: Ali Eydgahi, Department of Engineering and Aviation Sciences, University of Maryland Eastern Shore, Princess Anne, MD 853, Jalal Habibi, Department of Electrical and Computer Engineering, University of Tehran, Tehran, I.R. IRAN, Behzad Moshiri, Control & Intelligent Processing, Center of Excellency, ECE Department, Faculty of Engineering, University of Tehran, Tehran, I.R. IRAN, Abstract This paper presents a MATLAB-based toolbox with a Graphical User Interface (GUI), which can be used to compute reduced models of a large system by using one of the twenty order reduction techniques available in the toolbox. The methods that have been implemented in the toolbox include the Padé, Routh, Cauer, Continued-Fraction Expansion, and algorithms that provide mixtures of these techniques. Upon execution of the toolbox, a GUI will appear with four frames named Methods, High Order System, Output Options, and Results. In Methods frame, one or more of the reduction techniques can be selected. The high order system is defined by its transfer function numerator and denominator in High Order System frame. In this frame, a number of high order benchmark systems, in a popup menu, are provided that can be selected by the user for testing purposes. In Output Options frame, the order of the reduced model and the types of the desired output plots are selected. The impulse response and step response are available as an output. In the output plot, the software shows and compares the reduced models responses with nominal system response. The graphs of step and/or impulse responses will be appeared in two different windows. The numerator and denominator of the computed reduced order models will be shown in Results frame. Index Terms Control systems, graphical user interface, MATLAB toolbox, model order reduction, simulation. INTRODUCTION The modeling of complex dynamic systems is one of the most important subjects in engineering. A model is often too complicated to be used in real problems, so approximation procedures based on physical considerations or mathematical approaches are used to achieve simpler models than the original one. The subject of model reduction is very important to engineers and scientists working in many fields of engineering, especially, for those who work in the process control area. In control engineering field, model reduction techniques are fundamental for the design of controllers where particular numerically heavy procedures are involved. This would provide the designer with low order controllers that have less hardware requirements. Efforts towards obtaining low-order models from high-degree systems are related to the aims of deriving stable reduced-order models from stable original ones and assuring that the reduced-order model matches some quantities of the original one. In past decades, many papers in the literature have efficiently addressed these two objectives[]-[7]. Polynomial reduction methods are one of the important groups of the reduction techniques that are applied in the frequency domain. They are used to approximate low-order transfer functions where the model coefficients are chosen according to various criteria. The frequency domain model-reduction techniques that are mainly based on polynomial manipulations are Padé-type approximations, continued fraction expansion methods, model reduction using the Routh stability criterion, model reduction based on the Routh table criterion, model order reduction algorithms using stability equations, mixed methods of approximation, and energy-based methods. Popular software packages such as MATLAB, MAPLE, MATHEMATICA, and PSPICE are essential tools for analysis and design of an engineering system. The benefits of using these packages has been confirmed by the number of new and revised textbooks that have incorporated new exercises and problems based on them. These packages have been found to greatly enhance readers understanding of the materials and can prepare them for design and development of more complex system, much like the ones that are seen in industry. In this paper, a MATLAB-based toolbox for model order reduction with a Graphical User Interface is presented. The toolbox has options for twenty order reduction methods such as Padé approximation, Cauer continued-fraction expansion, July 5, 003, Valencia, Spain.

2 Routh stability criterion, differntiation method, and a number of different mixed methods of frequency techniques. Some MATLAB codes of order reduction methods have been adopted from []. MODEL ORDER REDUCTION METHODS In this section, brief description of the order reduction techniques that have been implemented in the toolbox along with their advantages and disadvantages are presented. Padé technique[3]: This approach is based on matching a few coefficients of the power series expansion, about s 0, of the reduced-order model with the corresponding coefficients of the original model. The power series expansion of a transfer function about s 0 is an expansion of the form: H ( c0 cs cs L () This method is computationally simple and provides perfect agreement between the steady state of the output of the system and the model for polynomial inputs belonging to the class of form α t i. Serious drawback of Padé-based approximants i is that unstable reduced-order models can be generated from original asymptotically stable high order system. Cauer first, second, and third form[-6]: Continued-fraction expansion method is one of the most attractive methods for the order reduction of transfer functions. It has many useful properties such as computational simplicity, the fitting of the time moments, and the preservation of the steady-state responses for polynomial inputs of the form α it i. There are three basic forms for continued fraction representation of a transfer function. Each form concentrates on a different frequency range for the approximation of the original transfer function. They are referred to as the three Cauer forms. The Cauer first form is represented as: H ( () fs f f 3s f O Reduced order model can be obtained by truncating the last terms of expansion. It has been shown that the reduced order model obtained from the Cauer first form cannot preserve the steady state value of the original model. The restriction to the original transfer function is that it requires the order of the numerator to be one less than that of the denominator. Otherwise, the Cauer first form does not exist. The reduced order model from the Cauer first form is the matching of the original transfer function by the Markov parameters, which is the power series expansion over s. The Cauer second form is given as: H ( (3) h h s h3 h s O It has been shown that this kind of continued-fraction is equivalent to the Maclaurin expansion about s 0. This form mathes the steady state reponses, but does not provide a good match for the initial transient response. The Cauer third form also referred to as the mixed Cauer expansion combines the Cauer first and second forms. It is given by: July 5, 003, Valencia, Spain.

3 H ( () h fs h f s h3 f 3s h f s O The Cauer third form can only be applied when the order of the numerator of the original transfer function is one less than that of the denominator. The mixed form is equivalent to the Maclaurin series about s 0 and s. The reduced order model matches the first m time moments and m Markov parameters of the original model, where m is the order of reduced order model. This method has one serious disadvantage. It may produce unstable (non-minimum phase) reduced models even though when the original high-order system is stable (minimum phase). Order reduction based on stability equation[7]: This method is based on forming stability equations of the numerator and denominator polynomials of high-order system, discarding the non-dominant poles and zeros, and obtaining the reducedorder model. This approach preserves the stability of reduced-order model. The disadvantages are related to the fact that this method is founded on heuristic considerations and besides the structural stability that is preserved in the reduced-order model, further correlations between the original and the reduced model are lost. Differentiation method[8]: In this approach, the reciprocals of numerator and denominator polynomials of the highorder transfer function are differentiated suitably many times to yield the coefficients of the reduced order transfer function. This method is computationally simple and is equally applicable to unstable and nonminimum phase systems. Aggregation method[9]: The aggregation method is based on the concept of combining a defined set of state variables of the original system, involving chosen weighting factors. The set of state variables (z), often computed from the original system state vector (x), using an aggregation matrix (L), by z Lx. The adequacy of finding the suitable aggregation state vector, and therefore the matrix L, depends upon the objective that the designer wants to reach such as open-loop or closed-loop performance, optimisation of certain indexes, etc. Order reduction based on the Routh stability criterion[0]: The system with transfer function M ( can be expanded into the so-called α and β parameters as: M ( β f( β f( f( L β n f( f( L fn( (5) where n is the degree of high-order system, β are constants, and the functions f i ( are defined as: i fi( (6) αis αi O αn s αns the first term is ( α s ) instead of α s. Order reduction based on the Routh stability criterion is to compute αi and βi parameters and, then, reconstructing the reduced-order models with certain algorithms after dropping out a few αi and βi parameters. This method cannot be applied to non-asymptotically stable systems. SRAM (Simplified Routh Approximation Method)[5-7]: This method obtains the reduced-order models using only α parameters, where α parameters are defined above. The simulated version of SRAM in this toolbox also retains some time moments of the full system in the reduced model. Order reduction based on real Schur form decomposition[]: Model order reduction, based on measures of controllability and observability, has been proposed by Moore[] and extended and generalized by other researchers. In this approach, the system is transformed using balancing transformation such that the input-to-state coupling and the state-tooutput coupling are weighted equally so that state components which are weakly coupled to both input and output are discarded and ones which are strongly coupled to both input and output are retained and used as a low-order approximation for the high order system. 3 July 5, 003, Valencia, Spain.

4 Calculating the balancing transformation in Moore s method, tends to be ill-conditioned, especially when the original system is non-minimal or when it has very nearly uncontrollable or unobservable modes. So, an alternative algorithm has been proposed based on finding the eigenspaces associated with large eigenvalues of the cross-gramian matrix using the real Schur-form decomposition. This method avoids calculating the balancing transformation. Stable Padé[3]: This approach combines the Routh-Hurwitz array method and the Padé approximation technique to obtain stable reduced-order models. It is computationally efficient and does not require computation of the poles of the original system for dominant pole retention.this method is not the only method that solves the stability problem of Padé approximation. Condensed continued-fraction method[]: Much attention has been diverted to overcoming stability problem of order reduction techniques based on continued-fraction expansion. Condensed CFE proposed to improve the chances of stability in the reduced model. One of the main causes of instability in the reduced-order models obtained by Cauer continued fractons is thought to be the dependence of the reduced denominator on the full numerator. This is especially true if the numerator has zeros in the right half-plane. Condensed CFE method proposed to use Cauer CFE only on D (, (where D( is the denominator) instead of G ( N( D(. The effect of N( on the reduced denuminator is thus removed, and the k th-order approximants k of D( can then be multiplied by N (, truncating the numerator after the term in s. This latter operation ensures that the first k time moments are preserved in the reduced model. However it should be remembered that unstable poles might still occur, and experiences have shown that lightly damped, highly oscillatory systems tend to cause most of the problems. Moment matching technique[5]: The moment matching technique is in fact the matching between the time moments of the impulse response of the original model with those of the reduced order model. This reduction method is not restricted to rational transfer functions. Since other functions, once written in the Maclaurin form, can also be reduced using the moment matching method. This approach, though simple in theory and computations, seems to have serious drawbacks. For example stable reduced-order models may be obtained by high-order unstable system. Truncation approach[6, 7]: From the viewpoint of condensed CFE order reduction approach, if [ p, q] is taken to be the Padé approximant about s 0 whose numerator and denominator degrees are p and q respectively, then the following result holds for a k -th order approximants of D (. If D s a a s a s (7) n ( ) 0 L n then the [ m, k] approximant of D (, where 0 m k, is equivalent to the two-point Padé/continued-fraction approximant of D( matching ( k m) terms about s 0 and ( k m) terms about s. This follows from the fact that the first ( n ) terms of D( about s are all zero. In the special case of m 0, let the approximant be given by e e s L e s (8) k k k The series expansion of this must be identical to that of D( for all terms in s. By equating the two series it is easily seen that e i ai for i 0,,, K, k. Forming the final reduced transfer function by multiplying N( and truncating after k the term in s shows that the model is identical to that obtained by simple truncation of terms in the full denominator and numerator, respecively.yeung in [7] has shown that for the truncation method if the poles of the system are well damped (i.e. real part > imaginary part ), then stability is guaranteed in the reduced model. Factor Divivision (Padé with dominant mode retention)[8]: This approach obtains the numerator of the reducedmodel in a direct way. It avoids to find system time moments and does not solve the Padé equaions, but obtains the reduced models that still retain the initial time moments of the full system. This method retains the dominant modes of the system and has been claimed to be often quicker to use than the similar method of [9] at least when the poles of the high-order system are known. Mixed method [0]: Generally speaking, mixed methods are developed to ensure the stability of reduced models. They often use criteria which assure the stability of the reduced-order model if the original system has this property. The numerator is obtained by using Padé method or the continued-fraction expansion. July 5, 003, Valencia, Spain.

5 In first mixed method, the denominator of the reduced-order transfer function is formed by the dominant poles of highorder system, while the numerator is obtained by fitting the first r time moments of reduced-order system and the high-order one. Mixed method []: In this approach, Routh stability criterion is used to compute the denominator of the reduced model. The numerator is obtained as in the first mixed method. This approach avoids the drawback of the first method in its need to the computation of system poles. Mixed method 3[]: This method is based on determining the denominator using the Routh array method. The numerator is obtaind by using Padé approximant technique. Mixed method [3]: This method takes into account both the adavantages of the stability-equations method for obtaining a reduced-order stable system, and the Padé method for finding the numerator. So this way it fits a number of the Markov parameters and time moments of the original transfer function. Mixed method 5[]: In this method, the differentiation procedure[8] is used to obtain the reduced-order denominator, thus ensuring stability preservation, while the numerator is obtained by matching the first MacLaurin expansion coefficients of the original system and model transfer functions using Padé method. DESCRIPTION OF THE TOOLBOX Upon execution of the toolbox, a GUI will appear with four frames named 'Methods', 'High Order System', 'Output Options', and 'Results'. Figure shows the workspace of GUI. In 'Methods' frame, one or more of the reduction techniques can be selected. The high order system is defined by its transfer function numerator and denominator in 'High Order System' frame. In this frame, a number of high order benchmark systems, in a popup menu, are provided that can be selected by the user for testing purposes. In 'Output Options' frame, the order of the reduced model and the types of the desired output plots are selected. The impulse response and step response are available as an output. In the output plot, the software shows and compares the responses of the reduced models with nominal system response. By clicking 'Run' push button, the software will start the computation of the reduced order model. Step and/or impulse responses graphs will be appeared in two different windows. The numerator and denominator of the computed reduced order models will be shown in 'Results' frame. Computing another model or changing some coefficients of the original model can easily be done in GUI. A 'Reset' push button is also provided to clear all frames and start computations for a new system. An example is provided to demonstrate the performance of the order reduction algorithms of the toolbox. EXAMPLE Here we present an illustrative example to show the operation of the order reduction toolbox. Some methods inherently put some restrictions on high order system. For example Cauer first and third forms restrict the relative degree of high order system to be equal to one. So, as notified previously, some standard high order systems have been considered in the toolbox as benchmark. Assume the high-order system from benchmark of the form: 3 ( s 35s 9s 093s 700 H (9) s 9s 66s 9s 09s 5s 68s 5856s 69s 700 and let s use Cauer second form, Routh stability criterion, Schur decomposition, and condensed CFE methods to find the reduced-order models of order 5. Figure shows the impulse responses for the nominal system and for reduced-order models obtained using selected methods. It can be seen that the Schur decomposition method produces most accurate model. REFERENCES [] Fortuna, L., Nunnari, G. and Gallo, A., "Model order reduction techniques with application in electrical engineering", Springer-Verlag, 99. [] Eydgahi, A.M., A mixed method for the simplification of high order systems, Proceedings of American Control Conference, Atlanta, GA, June 5-7, 988, pp July 5, 003, Valencia, Spain.

6 [3] Bultheel, A. and Van Barel, M., " Padé techniques for model reduction in linear system theory: a survey", Journal of Computational and Applied Mathematics, Vol., 986, pp [] Chen, C. F. and Shieh, L. S., "A novel approach to linear model simplification", Int. J. Control, Vol. 8, 968, pp [5] Shieh, L. S. and Goldman, M.J., "Continued Fraction Expansion and Inversion of the Cauer Third Form", IEEE Trans. Circuits and Systems, May 97, pp [6] Rathore, T. S., Singhi, B. M. and Kibe, A. V., "Continued fraction inversion and expansion", IEEE Trans. Automat. Contr., Vol., No., April 979, pp [7] Chen, T. C., Chang, C. Y. and Han, K. W., "Reduction of transfer functions by the stability-equation method", J. Franklin Inst., 979, pp [8] Gutman, P., Mannerfelt, C. F., Molander, P., "Contributions to the model reduction problem", IEEE Trans. Automat. Contr., Vol. 7, April 98, pp [9] Aoki, M., "Some approximation methods for estimation and control of large scale systems", IEEE Trans. Automat. Contr., Vol. 3, February 978, pp [0] Hutton, M. F. and Friedland, B., "Routh approximations for reducing order of linear time invariant systems", IEEE Trans. Automat. Contr., Vol. 0, 975, pp [] Aldhaheri, R.W., "Model order reduction via real Schur-form decomposition", Int. J. Control, Vol. 53, No. 3, 99, pp [] Moore, B. C., "Principal component analysis in linear systems: controllability, observability and model reduction", IEEE Trans. Automat. Contr., Vol. 6, January 98, pp [3] Pal, J., Stable reduced-order Padé approximants using the Routh-Hurwitz array, Electronic Letters, Vol. 5, No. 8, April 979, pp [] Lucas, T. N., "Model reduction by condensed continued-fraction method", Electronic Letters, Vol., No.6, August 985, pp [5] Lal, M. and Mitra, R., "Simplification of large system dynamics using a moment evaluation algorithm ", [6] Shamash, Y., "Truncation method of reduction: a viable alternative", Electronic Letters, Vol. 7, 98, pp [7] Yeung, K. S., "Stability of reduced model obtained by truncation", ibid., Vol. 7, 98, pp [8] Lucas, T. N., "Factor Division: a useful algorithm in model reduction", IEE Proceedings, Vol. 30, Pt. D, No. 6, November 983, pp [9] Shamash, Y., "Linear system reduction by Padé approximation to allow retention of dominant modes", Int. J. Control, Vol., 975, pp [0] Davison, E. J., Chidambara, M. R., "On a method for simplifying linear dynamic systems", IEEE Trans. Automat. Contr., 967, pp. 9-. [] Shamash, Y., "Model reduction using the Routh stability criterion and the Padé approximation technique", Int. J. Contr., Vol., No. 3, 975, pp [] Pal, J. and Ray, L.N., "Stable Padé approximant in multivariable systems using a mixed method", Proc. IEEE, Vol. 68, No., 980, pp [3] Chen, T. C., Chang, C. Y. and Han, K. W., "Model reduction using the stability-equation and the Padé approximation method", J. Franklin Inst., 980. [] Lepschy, A., Viaro, U., "A note on model reduction problem", IEEE Tans. Automat. Contr., Vol. 8, April 983, pp [5] Sastry, G. V. K. R. and Krishnamurthy, V., "Relative stability using simplified Routh approximation method (SRAM)", J. Inst. Electron. & Telecommun. Eng., Vol. 33, 987. [6] Sastry, G. V. K. R. and Krishnamurthy, V., "Biased model reduction by simplified Routh approximation method", Electronic Letters, Vol. 3, No. 0, September 987, pp [7] Sastry, G. V. K. R. and Krishnamurthy, V., "State-space models using simplified Routh approximation method", Electronic Letters, Vol. 3, No., November 987, pp July 5, 003, Valencia, Spain.

7 FIGURES AND TABLES FIGURE. WORKSPACE OF THE GRAPHICAL USER INTERFACE OF THE TOOLBOX. FIGURE. IMPULSE RESPONSES FOR HIGH-ORDER SYSTEM OF EXAMPLE AND REDUCED-ORDER MODELS USING FOUR METHODS OF TOOLBOX. 7 July 5, 003, Valencia, Spain.

Model Order Reduction of Interval Systems using Mihailov Criterion and Factor Division Method

Model Order Reduction of Interval Systems using Mihailov Criterion and Factor Division Method Model Order Reduction of Interval Systems using Mihailov Criterion Factor Division Method D. Kranthi Kumar Research Scholar S. K. Nagar Professor J. P. Tiwari Professor ABSTRACT This paper presents a mixed

More information

TWO- PHASE APPROACH TO DESIGN ROBUST CONTROLLER FOR UNCERTAIN INTERVAL SYSTEM USING GENETIC ALGORITHM

TWO- PHASE APPROACH TO DESIGN ROBUST CONTROLLER FOR UNCERTAIN INTERVAL SYSTEM USING GENETIC ALGORITHM International Journal of Electrical and Electronics Engineering Research (IJEEER) ISSN:2250-155X Vol.2, Issue 2 June 2012 27-38 TJPRC Pvt. Ltd., TWO- PHASE APPROACH TO DESIGN ROBUST CONTROLLER FOR UNCERTAIN

More information

7.4 STEP BY STEP PROCEDURE TO DRAW THE ROOT LOCUS DIAGRAM

7.4 STEP BY STEP PROCEDURE TO DRAW THE ROOT LOCUS DIAGRAM ROOT LOCUS TECHNIQUE. Values of on the root loci The value of at any point s on the root loci is determined from the following equation G( s) H( s) Product of lengths of vectors from poles of G( s)h( s)

More information

Poles, Zeros and System Response

Poles, Zeros and System Response Time Response After the engineer obtains a mathematical representation of a subsystem, the subsystem is analyzed for its transient and steady state responses to see if these characteristics yield the desired

More information

Model Order Reduction of Interval Systems Using Mihailov Criterion and Cauer Second Form

Model Order Reduction of Interval Systems Using Mihailov Criterion and Cauer Second Form Model Order Reduction of Interval Systems Using Mihailov Criterion and Cauer Second Form D. Kranthi Kumar Research Scholar S. K. Nagar Professor J. P. Tiwari Professor ABSTRACT This paper presents a new

More information

DSP First. Laboratory Exercise #10. The z, n, and ˆω Domains

DSP First. Laboratory Exercise #10. The z, n, and ˆω Domains DSP First Laboratory Exercise #10 The z, n, and ˆω Domains 1 Objective The objective for this lab is to build an intuitive understanding of the relationship between the location of poles and zeros in the

More information

MIMO Toolbox for Matlab

MIMO Toolbox for Matlab MIMO Toolbox for Matlab Oskar Vivero School of Electric and Electronic Engineering University of Manchester, Manchester, UK M60 1QD Email: oskar.vivero@postgrad.manchester.ac.uk Jesús Liceaga-Castro Departamento

More information

LMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System

LMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System LMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System Gholamreza Khademi, Haniyeh Mohammadi, and Maryam Dehghani School of Electrical and Computer Engineering Shiraz

More information

ECE317 : Feedback and Control

ECE317 : 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 information

Computing Consecutive-Type Reliabilities Non-Recursively

Computing Consecutive-Type Reliabilities Non-Recursively IEEE TRANSACTIONS ON RELIABILITY, VOL. 52, NO. 3, SEPTEMBER 2003 367 Computing Consecutive-Type Reliabilities Non-Recursively Galit Shmueli Abstract The reliability of consecutive-type systems has been

More information

II. GENERALISED PROPOSED REDUCTION TECHNIQUE Considering the transfer function of an nth order discrete time interval system be denoted as

II. GENERALISED PROPOSED REDUCTION TECHNIQUE Considering the transfer function of an nth order discrete time interval system be denoted as Analysis and Design of High Order Discrete Time Interval Systems using New Order Reduction Technique via Least Squares Method and Time Moments Technique Kiran Kumar Kalyana 1, Sangeeta Kurman 2 1,2 Department

More information

EEE 184 Project: Option 1

EEE 184 Project: Option 1 EEE 184 Project: Option 1 Date: November 16th 2012 Due: December 3rd 2012 Work Alone, show your work, and comment your results. Comments, clarity, and organization are important. Same wrong result or same

More information

Automatic Control Systems theory overview (discrete time systems)

Automatic Control Systems theory overview (discrete time systems) Automatic Control Systems theory overview (discrete time systems) Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Motivations

More information

Graphical User Interface for Design Stabilizing Controllers

Graphical User Interface for Design Stabilizing Controllers Graphical User Interface for Design Stabilizing Controllers Petr Urban 1,2, Michael Šebek1 1 Department of Control Engineering, Faculty of Electrical Engineering, Czech Technical University, Prague 2 Institute

More information

EE 422G - Signals and Systems Laboratory

EE 422G - Signals and Systems Laboratory EE 4G - Signals and Systems Laboratory Lab 9 PID Control Kevin D. Donohue Department of Electrical and Computer Engineering University of Kentucky Lexington, KY 40506 April, 04 Objectives: Identify the

More information

Signal Processing First Lab 11: PeZ - The z, n, and ˆω Domains

Signal Processing First Lab 11: PeZ - The z, n, and ˆω Domains Signal Processing First Lab : PeZ - The z, n, and ˆω Domains The lab report/verification will be done by filling in the last page of this handout which addresses a list of observations to be made when

More information

DECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION

DECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION Journal of ELECTRICAL ENGINEERING, VOL. 58, NO. 6, 2007, 307 312 DECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION Szabolcs Dorák Danica Rosinová Decentralized control design approach based on partial

More information

1 Introduction & Objective. 2 Warm-up. Lab P-16: PeZ - The z, n, and O! Domains

1 Introduction & Objective. 2 Warm-up. Lab P-16: PeZ - The z, n, and O! Domains DSP First, 2e Signal Processing First Lab P-6: PeZ - The z, n, and O! Domains The lab report/verification will be done by filling in the last page of this handout which addresses a list of observations

More information

ECEN 605 LINEAR SYSTEMS. Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability 1/27

ECEN 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 information

Chapter 5 THE APPLICATION OF THE Z TRANSFORM. 5.3 Stability

Chapter 5 THE APPLICATION OF THE Z TRANSFORM. 5.3 Stability Chapter 5 THE APPLICATION OF THE Z TRANSFORM 5.3 Stability Copyright c 2005- Andreas Antoniou Victoria, BC, Canada Email: aantoniou@ieee.org February 13, 2008 Frame # 1 Slide # 1 A. Antoniou Digital Signal

More information

POSITIVE REALNESS OF A TRANSFER FUNCTION NEITHER IMPLIES NOR IS IMPLIED BY THE EXTERNAL POSITIVITY OF THEIR ASSOCIATE REALIZATIONS

POSITIVE REALNESS OF A TRANSFER FUNCTION NEITHER IMPLIES NOR IS IMPLIED BY THE EXTERNAL POSITIVITY OF THEIR ASSOCIATE REALIZATIONS POSITIVE REALNESS OF A TRANSFER FUNCTION NEITHER IMPLIES NOR IS IMPLIED BY THE EXTERNAL POSITIVITY OF THEIR ASSOCIATE REALIZATIONS Abstract This letter discusses the differences in-between positive realness

More information

IC6501 CONTROL SYSTEMS

IC6501 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 information

Module 3F2: Systems and Control EXAMPLES PAPER 2 ROOT-LOCUS. Solutions

Module 3F2: Systems and Control EXAMPLES PAPER 2 ROOT-LOCUS. Solutions Cambridge University Engineering Dept. Third Year Module 3F: Systems and Control EXAMPLES PAPER ROOT-LOCUS Solutions. (a) For the system L(s) = (s + a)(s + b) (a, b both real) show that the root-locus

More information

DSP First Lab 11: PeZ - The z, n, and ωdomains

DSP First Lab 11: PeZ - The z, n, and ωdomains DSP First Lab : PeZ - The, n, and ωdomains The lab report/verification will be done by filling in the last page of this handout which addresses a list of observations to be made when using the PeZ GUI.

More information

Automatic Control Systems

Automatic Control Systems Automatic Control Systems Edited by Dr. Yuri Sokolov Contributing Authors: Dr. Victor Iliushko, Dr. Emaid A. Abdul-Retha, Mr. Sönke Dierks, Dr. Pascual Marqués. Published by Marques Aviation Ltd Southport,

More information

FUZZY LOGIC CONTROL Vs. CONVENTIONAL PID CONTROL OF AN INVERTED PENDULUM ROBOT

FUZZY LOGIC CONTROL Vs. CONVENTIONAL PID CONTROL OF AN INVERTED PENDULUM ROBOT http:// FUZZY LOGIC CONTROL Vs. CONVENTIONAL PID CONTROL OF AN INVERTED PENDULUM ROBOT 1 Ms.Mukesh Beniwal, 2 Mr. Davender Kumar 1 M.Tech Student, 2 Asst.Prof, Department of Electronics and Communication

More information

Feedback Control of Linear SISO systems. Process Dynamics and Control

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

More information

Software Engineering/Mechatronics 3DX4. Slides 6: Stability

Software 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 information

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

(b) A unity feedback system is characterized by the transfer function. Design a suitable compensator to meet the following specifications: 1. (a) The open loop transfer function of a unity feedback control system is given by G(S) = K/S(1+0.1S)(1+S) (i) Determine the value of K so that the resonance peak M r of the system is equal to 1.4.

More information

Control Systems Engineering ( Chapter 6. Stability ) Prof. Kwang-Chun Ho Tel: Fax:

Control Systems Engineering ( Chapter 6. Stability ) Prof. Kwang-Chun Ho Tel: Fax: Control Systems Engineering ( Chapter 6. Stability ) 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 : How to determine

More information

CHAPTER 5 ROBUSTNESS ANALYSIS OF THE CONTROLLER

CHAPTER 5 ROBUSTNESS ANALYSIS OF THE CONTROLLER 114 CHAPTER 5 ROBUSTNESS ANALYSIS OF THE CONTROLLER 5.1 INTRODUCTION Robust control is a branch of control theory that explicitly deals with uncertainty in its approach to controller design. It also refers

More information

CONTROL SYSTEMS ENGINEERING Sixth Edition International Student Version

CONTROL 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

CYBER EXPLORATION LABORATORY EXPERIMENTS

CYBER EXPLORATION LABORATORY EXPERIMENTS CYBER EXPLORATION LABORATORY EXPERIMENTS 1 2 Cyber Exploration oratory Experiments Chapter 2 Experiment 1 Objectives To learn to use MATLAB to: (1) generate polynomial, (2) manipulate polynomials, (3)

More information

REVISED GCSE Scheme of Work Mathematics Higher Unit 4. For First Teaching September 2010 For First Examination Summer 2011

REVISED GCSE Scheme of Work Mathematics Higher Unit 4. For First Teaching September 2010 For First Examination Summer 2011 REVISED GCSE Scheme of Work Mathematics Higher Unit 4 For First Teaching September 2010 For First Examination Summer 2011 Version 1: 28 April 10 Version 1: 28 April 10 Unit T4 Version 1: 28 April 10 Unit

More information

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

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

More information

FREQUENCY DOMAIN IDENTIFICATION TOOLBOX

FREQUENCY DOMAIN IDENTIFICATION TOOLBOX NASA Technical Memorandum 109039 FREQUENCY DOMAIN IDENTIFICATION TOOLBOX Lucas G. Horta and Jer-Nan Juang NASA Langley Research Center Hampton, Virginia 23681-0001 Chung-Wen Chen North Carolina State University

More information

Some special cases

Some 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 information

SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015

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

More information

PM diagram of the Transfer Function and its use in the Design of Controllers

PM diagram of the Transfer Function and its use in the Design of Controllers PM diagram of the Transfer Function and its use in the Design of Controllers Santiago Garrido, Luis Moreno Abstract This paper presents the graphical chromatic representation of the phase and the magnitude

More information

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

Stability Analysis of Nonlinear Systems with Transportation Lag

Stability Analysis of Nonlinear Systems with Transportation Lag ABSTRACT Stability Analysis of Nonlinear Systems with Transportation Lag Authors Dr.K.Sreekala 1, Dr.S.N.Sivanandam 2 1 Professor, Sreenarayana Gurukulam College of Engineering, Kadayiruppu,Kolenchery

More information

HANDOUT E.22 - EXAMPLES ON STABILITY ANALYSIS

HANDOUT E.22 - EXAMPLES ON STABILITY ANALYSIS Example 1 HANDOUT E. - EXAMPLES ON STABILITY ANALYSIS Determine the stability of the system whose characteristics equation given by 6 3 = s + s + 3s + s + s + s +. The above polynomial satisfies the necessary

More information

EEE 480 LAB EXPERIMENTS. K. Tsakalis. November 25, 2002

EEE 480 LAB EXPERIMENTS. K. Tsakalis. November 25, 2002 EEE 480 LAB EXPERIMENTS K. Tsakalis November 25, 2002 1. Introduction The following set of experiments aims to supplement the EEE 480 classroom instruction by providing a more detailed and hands-on experience

More information

AEAN ITEqATIVE TECHNIQUE TO STABILIZE A LINEAR TMEv VISWANATHAN SANKARAN NRC POSTDOCTORAL RESEARCH ASSOCIATE FLIGHT DYNAMICS AND CONTROL DIVISION

AEAN ITEqATIVE TECHNIQUE TO STABILIZE A LINEAR TMEv VISWANATHAN SANKARAN NRC POSTDOCTORAL RESEARCH ASSOCIATE FLIGHT DYNAMICS AND CONTROL DIVISION AD-A280 164 IIImmI UlI NASA TECHNICAL NASA TM X-71986 MEMORANDUM COPY NO. 00 0" (7% '-4 AEAN ITEqATIVE TECHNIQUE TO STABILIZE A LINEAR TMEv INVARIANT MULTIVARIABLE SYSTEM WITH UJTPUJT FEEDBACK VISWANATHAN

More information

80DB A 40DB 0DB DB

80DB A 40DB 0DB DB Stability Analysis On the Nichols Chart and Its Application in QFT Wenhua Chen and Donald J. Ballance Centre for Systems & Control Department of Mechanical Engineering University of Glasgow Glasgow G12

More information

CONTROL * ~ SYSTEMS ENGINEERING

CONTROL * ~ 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 information

Stability of Feedback Control Systems: Absolute and Relative

Stability 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 information

FIT100 Spring 01. Project 2. Astrological Toys

FIT100 Spring 01. Project 2. Astrological Toys FIT100 Spring 01 Project 2 Astrological Toys In this project you will write a series of Windows applications that look up and display astrological signs and dates. The applications that will make up the

More information

Software Engineering 3DX3. Slides 8: Root Locus Techniques

Software Engineering 3DX3. Slides 8: Root Locus Techniques Software Engineering 3DX3 Slides 8: Root Locus Techniques Dr. Ryan Leduc Department of Computing and Software McMaster University Material based on Control Systems Engineering by N. Nise. c 2006, 2007

More information

1. Introduction. Hang Qian 1 Iowa State University

1. Introduction. Hang Qian 1 Iowa State University Users Guide to the VARDAS Package Hang Qian 1 Iowa State University 1. Introduction The Vector Autoregression (VAR) model is widely used in macroeconomics. However, macroeconomic data are not always observed

More information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 12: Overview In this Lecture, you will learn: Review of Feedback Closing the Loop Pole Locations Changing the Gain

More information

Learn2Control Laboratory

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

More information

An Optimum Fitting Algorithm for Generation of Reduced-Order Models

An Optimum Fitting Algorithm for Generation of Reduced-Order Models An Optimum Fitting Algorithm for Generation of Reduced-Order Models M.M. Gourary 1, S.G. Rusakov 1, S.L. Ulyanov 1, M.M. Zharov 1, B.J. Mulvaney 2 1 IPPM, Russian Academy of Sciences, Moscow 1523, e-mail:

More information

CERTAIN ALGEBRAIC TESTS FOR ANALYZING APERIODIC AND RELATIVE STABILITY IN LINEAR DISCRETE SYSTEMS

CERTAIN ALGEBRAIC TESTS FOR ANALYZING APERIODIC AND RELATIVE STABILITY IN LINEAR DISCRETE SYSTEMS CERTAIN ALGEBRAIC TESTS FOR ANALYZING APERIODIC AND RELATIVE STABILITY IN LINEAR DISCRETE SYSTEMS Ramesh P. and Manikandan V. 2 Department of Electrical and Electronics Engineering, University College

More information

EE Experiment 11 The Laplace Transform and Control System Characteristics

EE Experiment 11 The Laplace Transform and Control System Characteristics EE216:11 1 EE 216 - Experiment 11 The Laplace Transform and Control System Characteristics Objectives: To illustrate computer usage in determining inverse Laplace transforms. Also to determine useful signal

More information

ECE317 : Feedback and Control

ECE317 : 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 information

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

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

Remember that : Definition :

Remember that : Definition : Stability This lecture we will concentrate on How to determine the stability of a system represented as a transfer function How to determine the stability of a system represented in state-space How to

More information

State will have dimension 5. One possible choice is given by y and its derivatives up to y (4)

State will have dimension 5. One possible choice is given by y and its derivatives up to y (4) A Exercise State will have dimension 5. One possible choice is given by y and its derivatives up to y (4 x T (t [ y(t y ( (t y (2 (t y (3 (t y (4 (t ] T With this choice we obtain A B C [ ] D 2 3 4 To

More information

INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad

INSTITUTE 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 information

Appendix 3B MATLAB Functions for Modeling and Time-domain analysis

Appendix 3B MATLAB Functions for Modeling and Time-domain analysis Appendix 3B MATLAB Functions for Modeling and Time-domain analysis MATLAB control system Toolbox contain the following functions for the time-domain response step impulse initial lsim gensig damp ltiview

More information

Results on stability of linear systems with time varying delay

Results on stability of linear systems with time varying delay IET Control Theory & Applications Brief Paper Results on stability of linear systems with time varying delay ISSN 75-8644 Received on 8th June 206 Revised st September 206 Accepted on 20th September 206

More information

Nonlinear System Analysis

Nonlinear System Analysis Nonlinear System Analysis Lyapunov Based Approach Lecture 4 Module 1 Dr. Laxmidhar Behera Department of Electrical Engineering, Indian Institute of Technology, Kanpur. January 4, 2003 Intelligent Control

More information

Solving Polynomial and Rational Inequalities Algebraically. Approximating Solutions to Inequalities Graphically

Solving Polynomial and Rational Inequalities Algebraically. Approximating Solutions to Inequalities Graphically 10 Inequalities Concepts: Equivalent Inequalities Solving Polynomial and Rational Inequalities Algebraically Approximating Solutions to Inequalities Graphically (Section 4.6) 10.1 Equivalent Inequalities

More information

CANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM

CANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM CANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM Ralf L.M. Peeters Bernard Hanzon Martine Olivi Dept. Mathematics, Universiteit Maastricht, P.O. Box 616, 6200 MD Maastricht,

More information

Table of Laplacetransform

Table 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 information

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014 Prof. Dr. Eleni Chatzi System Stability - 26 March, 24 Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can

More information

Control of Electromechanical Systems

Control of Electromechanical Systems Control of Electromechanical Systems November 3, 27 Exercise Consider the feedback control scheme of the motor speed ω in Fig., where the torque actuation includes a time constant τ A =. s and a disturbance

More information

Chemical 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 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 information

Quantitative Feedback Theory based Controller Design of an Unstable System

Quantitative Feedback Theory based Controller Design of an Unstable System Quantitative Feedback Theory based Controller Design of an Unstable System Chandrima Roy Department of E.C.E, Assistant Professor Heritage Institute of Technology, Kolkata, WB Kalyankumar Datta Department

More information

Dynamic circuits: Frequency domain analysis

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

More information

Stability and composition of transfer functions

Stability and composition of transfer functions Problem 1.1 Stability and composition of transfer functions G. Fernández-Anaya Departamento de Ciencias Básicas Universidad Iberoaméricana Lomas de Santa Fe 01210 México D.F. México guillermo.fernandez@uia.mx

More information

Model reduction in the physical domain

Model reduction in the physical domain 48 Model reduction in the physical domain A Y Orbak, O S Turkay2*,**, E Eskinat and K Youcef-Toumi3 Department of Industrial Engineering, Uludag University, Bursa, Turkey 2Department of Mechanical Engineering,

More information

Automatic Control Systems, 9th Edition

Automatic Control Systems, 9th Edition Chapter 7: Root Locus Analysis Appendix E: Properties and Construction of the Root Loci Automatic Control Systems, 9th Edition Farid Golnaraghi, Simon Fraser University Benjamin C. Kuo, University of Illinois

More information

SRV02-Series Rotary Experiment # 1. Position Control. Student Handout

SRV02-Series Rotary Experiment # 1. Position Control. Student Handout SRV02-Series Rotary Experiment # 1 Position Control Student Handout SRV02-Series Rotary Experiment # 1 Position Control Student Handout 1. Objectives The objective in this experiment is to introduce the

More information

Identification of Dominant and Non-Dominant States of Large Scale Discrete-Time Linear Systems and its Applications to Reduced Order Models

Identification of Dominant and Non-Dominant States of Large Scale Discrete-Time Linear Systems and its Applications to Reduced Order Models Applied Mathematical Sciences, Vol. 6, 2012, no. 116, 5755-5762 Identification of Dominant and Non-Dominant States of Large Scale Discrete-Time Linear Systems and its Applications to Reduced Order Models

More information

Research Article On Complementary Root Locus of Biproper Transfer Functions

Research Article On Complementary Root Locus of Biproper Transfer Functions Mathematical Problems in Engineering Volume 9, Article ID 779, pages doi:.55/9/779 Research Article On Complementary Root Locus of Biproper Transfer Functions Marcelo C. M. Teixeira, Edvaldo Assunção,

More information

Structural Properties of LTI Singular Systems under Output Feedback

Structural Properties of LTI Singular Systems under Output Feedback Structural Properties of LTI Singular Systems under Output Feedback Runyi Yu Dept. of Electrical and Electronic Eng. Eastern Mediterranean University Gazimagusa, Mersin 10, Turkey runyi.yu@emu.edu.tr Dianhui

More information

Fast Step-response Evaluation of Linear Continuous-time Systems With Time Delay in the Feedback Loop

Fast Step-response Evaluation of Linear Continuous-time Systems With Time Delay in the Feedback Loop Proceedings of the 7th World Congress The International Federation of Automatic Control Seoul, Korea, July 6-, 28 Fast Step-response Evaluation of Linear Continuous-time Systems With Time Delay in the

More information

Parallel Properties of Poles of. Positive Functions and those of. Discrete Reactance Functions

Parallel Properties of Poles of. Positive Functions and those of. Discrete Reactance Functions International Journal of Mathematical Analysis Vol. 11, 2017, no. 24, 1141-1150 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ima.2017.77106 Parallel Properties of Poles of Positive Functions and

More information

A new multidimensional Schur-Cohn type stability criterion

A new multidimensional Schur-Cohn type stability criterion A new multidimensional Schur-Cohn type stability criterion Ioana Serban, Mohamed Najim To cite this version: Ioana Serban, Mohamed Najim. A new multidimensional Schur-Cohn type stability criterion. American

More information

Research Article Convex Polyhedron Method to Stability of Continuous Systems with Two Additive Time-Varying Delay Components

Research Article Convex Polyhedron Method to Stability of Continuous Systems with Two Additive Time-Varying Delay Components Applied Mathematics Volume 202, Article ID 689820, 3 pages doi:0.55/202/689820 Research Article Convex Polyhedron Method to Stability of Continuous Systems with Two Additive Time-Varying Delay Components

More information

SECTION 5: ROOT LOCUS ANALYSIS

SECTION 5: ROOT LOCUS ANALYSIS SECTION 5: ROOT LOCUS ANALYSIS MAE 4421 Control of Aerospace & Mechanical Systems 2 Introduction Introduction 3 Consider a general feedback system: Closed loop transfer function is 1 is the forward path

More information

Further Results on Bose s 2D Stability Test

Further Results on Bose s 2D Stability Test Further Results on Bose s D Stability Test Minoru Yamada, Li Xu, Osami Saito Gifu National College of Technology, Shinsei-cho, Gifu, 50-0495 Japan, Faculty of System Science and Technology, Akita Prefecture

More information

56 CHAPTER 3. POLYNOMIAL FUNCTIONS

56 CHAPTER 3. POLYNOMIAL FUNCTIONS 56 CHAPTER 3. POLYNOMIAL FUNCTIONS Chapter 4 Rational functions and inequalities 4.1 Rational functions Textbook section 4.7 4.1.1 Basic rational functions and asymptotes As a first step towards understanding

More information

Analysis and Design of a Controller for an SMIB Power System via Time Domain Approach

Analysis and Design of a Controller for an SMIB Power System via Time Domain Approach Indian Journal of Science and Technology, Vol 9(44), 10.17485/ijst/2016/v9i44/102985, November 2016 ISSN (Print) : 0974-6846 ISSN (Online) : 0974-5645 Analysis and Design of a Controller for an SMIB Power

More information

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

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

More information

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

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

More information

Uncertainty and Robustness for SISO Systems

Uncertainty and Robustness for SISO Systems Uncertainty and Robustness for SISO Systems ELEC 571L Robust Multivariable Control prepared by: Greg Stewart Outline Nature of uncertainty (models and signals). Physical sources of model uncertainty. Mathematical

More information

Lyapunov Stability of Linear Predictor Feedback for Distributed Input Delays

Lyapunov Stability of Linear Predictor Feedback for Distributed Input Delays IEEE TRANSACTIONS ON AUTOMATIC CONTROL VOL. 56 NO. 3 MARCH 2011 655 Lyapunov Stability of Linear Predictor Feedback for Distributed Input Delays Nikolaos Bekiaris-Liberis Miroslav Krstic In this case system

More information

CHAPTER # 9 ROOT LOCUS ANALYSES

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

KINGS COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING

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

More information

LMI based Stability criteria for 2-D PSV system described by FM-2 Model

LMI based Stability criteria for 2-D PSV system described by FM-2 Model Vol-4 Issue-018 LMI based Stability criteria for -D PSV system described by FM- Model Prashant K Shah Department of Electronics Engineering SVNIT, pks@eced.svnit.ac.in Abstract Stability analysis is the

More information

Discretization of Fractional-Order Differentiators and Integrators

Discretization of Fractional-Order Differentiators and Integrators Preprints of the 19th World Congress The International Federation of Automatic Control Cape Town, South Africa. August 24-29, 14 Discretization of Fractional-Order Differentiators and Integrators Reyad

More information

FATIMA MICHAEL COLLEGE OF ENGINEERING & TECHNOLOGY

FATIMA 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 information

MOMENT functions are used in several computer vision

MOMENT functions are used in several computer vision IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 13, NO. 8, AUGUST 2004 1055 Some Computational Aspects of Discrete Orthonormal Moments R. Mukundan, Senior Member, IEEE Abstract Discrete orthogonal moments

More information

Multivariable Calculus with MATLAB

Multivariable Calculus with MATLAB Multivariable Calculus with MATLAB Ronald L. Lipsman Jonathan M. Rosenberg Multivariable Calculus with MATLAB With Applications to Geometry and Physics Ronald L. Lipsman Department of Mathematics University

More information

Appendix A Prototypes Models

Appendix A Prototypes Models Appendix A Prototypes Models This appendix describes the model of the prototypes used in Chap. 3. These mathematical models can also be found in the Student Handout by Quanser. A.1 The QUANSER SRV-02 Setup

More information

Classroom Tips and Techniques: Series Expansions Robert J. Lopez Emeritus Professor of Mathematics and Maple Fellow Maplesoft

Classroom Tips and Techniques: Series Expansions Robert J. Lopez Emeritus Professor of Mathematics and Maple Fellow Maplesoft Introduction Classroom Tips and Techniques: Series Expansions Robert J. Lopez Emeritus Professor of Mathematics and Maple Fellow Maplesoft Maple has the ability to provide various series expansions and

More information

Stability Analysis Techniques

Stability Analysis Techniques Stability Analysis Techniques In this section the stability analysis techniques for the Linear Time-Invarient (LTI) discrete system are emphasized. In general the stability techniques applicable to LTI

More information