Cleveland State University MCE441: Intr. Linear Control Systems. Lecture 6: The Transfer Function Poles and Zeros Partial Fraction Expansions

Similar documents
Classroom Tips and Techniques: The Partial-Fraction Decomposition. Robert J. Lopez Emeritus Professor of Mathematics and Maple Fellow Maplesoft

Algebra I: A Fresh Approach. By Christy Walters

Algebra I: A Fresh Approach. By Christy Walters

1. A tendency to roll or heel when turning (a known and typically constant disturbance) 2. Motion induced by surface waves of certain frequencies.

Math 4. Unit 1: Conic Sections Lesson 1.1: What Is a Conic Section?

COMPLEX NUMBERS. Powers of j. Consider the quadratic equation; It has no solutions in the real number system since. j 1 ie. j 2 1

Process Dynamics, Operations, and Control Lecture Notes - 20

Lab 4: Root Locus Based Control Design

Compensator Design for Speed Control of DC Motor by Root Locus Approach using MATLAB

Operations on Radical Expressions; Rationalization of Denominators

EXPERIMENTAL STUDY ON THE HYDRODYNAMIC BEHAVIORS OF TWO CONCENTRIC CYLINDERS

HOW FAST/FAR DOES FLY LINE FALL? N. Perkins of the University of Michigan, March 2003

Novel empirical correlations for estimation of bubble point pressure, saturated viscosity and gas solubility of crude oils

The Simple Linear Regression Model ECONOMETRICS (ECON 360) BEN VAN KAMMEN, PHD

Advanced Hydraulics Prof. Dr. Suresh A. Kartha Department of Civil Engineering Indian Institute of Technology, Guwahati

OPTIMIZATION OF A WAVE CANCELLATION MULTIHULL SHIP USING CFD TOOLS

RATE OF CHANGE AND INSTANTANEOUS VELOCITY

Lecture 5. Optimisation. Regularisation

Section I: Multiple Choice Select the best answer for each problem.

Vibration of floors and footfall analysis

Control System Nur Istianah-THP-FTP-UB-2014

Chapter 12 Practice Test

Jasmin Smajic 1, Christian Hafner 2, Jürg Leuthold 2, March 16, 2015 Introduction to Finite Element Method (FEM) Part 1 (2-D FEM)

Applied Fluid Mechanics

Calculation of Trail Usage from Counter Data

SOFTWARE. Sesam user course. 12 May 2016 HydroD Hydrostatics & Stability. Ungraded SAFER, SMARTER, GREENER DNV GL 2016

DEVELOPMENT OF A THREE-DIMENSIONAL INVERSE SAIL DESIGN METHOD

APPROVED FACILITY SCHOOLS CURRICULUM DOCUMENT SUBJECT: Mathematics GRADE: 6. TIMELINE: Quarter 1. Student Friendly Learning Objective

ME 200 Thermodynamics I Spring 2010 (Last) (First) Thermo Number: CIRCLE YOUR LECTURE BELOW

MAE 322 Machine Design Lecture 5 Fatigue - 2. Dr. Hodge Jenkins Mercer University

INSTRUMENT INSTRUMENTAL ERROR (of full scale) INSTRUMENTAL RESOLUTION. Tutorial simulation. Tutorial simulation

Implementation of Height Measurement System Based on Pressure Sensor BMP085

Exercise (3): Open Channel Flow Rapidly Varied Flow

Wind Flow Model of Area Surrounding the Case Western Reserve University Wind Turbine

Math Spring Operational Geometry PBA Item #18 Classmates in the Pool VH003506

T H E R M O P T I M CALCULATION OF MOIST GAS FROM EXTERNAL CLASSES VERSION JAVA 1.5 R. GICQUEL MARCH 2007

Dynamic Component of Ship s Heeling Moment due to Sloshing vs. IMO IS-Code Recommendations

Engineering Note. Algorithms. Overview. Detailed Algorithm Description. NeoFox Calibration and Measurement. Products Affected: NeoFox

14 April 2016 Dr. John Hedengren 350 CB Brigham Young University Provo, UT 84606

The purpose of this experiment is to find this acceleration for a puck moving on an inclined air table.

The Lecture deals with:

LABORATORY EXERCISE 1 CONTROL VALVE CHARACTERISTICS

Using GPOPS-II to optimize sum of squared torques of a double pendulum as a prosthesis leg. Abstract

Chapter 7. SCIIB Pressure Sensor Performance Evaluations: Experiments, Results and Discussions

Gerald D. Anderson. Education Technical Specialist

International Journal of Technical Research and Applications e-issn: , Volume 4, Issue 3 (May-June, 2016), PP.

March Madness Basketball Tournament

Section 4.2 Objectives

A PHASE-AMPLITUDE ITERATION SCHEME FOR THE OPTIMIZATION OF DETERMINISTIC WAVE SEQUENCES

PSY201: Chapter 5: The Normal Curve and Standard Scores

DYNAMIC BEHAVIOR OF SPRING-LOADED PRESSURE RELIEF VALVE: NUMERICAL AND EXPERIMENTAL ANALYSIS

Simplifying Radical Expressions and the Distance Formula

Energy and mass transfer in gas-liquid reactors.

LEARNING OBJECTIVES. Overview of Lesson. guided practice Teacher: anticipates, monitors, selects, sequences, and connects student work

Aerodynamic Analyses of Horizontal Axis Wind Turbine By Different Blade Airfoil Using Computer Program

DECISION-MAKING ON IMPLEMENTATION OF IPO UNDER TOPOLOGICAL UNCERTAINTY

Development of Fluid-Structure Interaction Program for the Mercury Target

Organisation Internationale de Métrologie Légale

Navigate to the golf data folder and make it your working directory. Load the data by typing

Aerodynamic study of a cyclist s moving legs using an innovative approach

Transactions on Ecology and the Environment vol 12, 1996 WIT Press, ISSN

Midterm Exam 1, section 2. Thursday, September hour, 15 minutes

Preparation for Salinity Control ME 121

March Madness Basketball Tournament

2 are both ways of saying a ratio of 2 to 5

Incompressible Flow over Airfoils

Bridge Decomposition of Restriction Measures Tom Alberts joint work with Hugo Duminil (ENS) with lots of help from Wendelin Werner University of

Fluid-Structure Interaction Analysis of a Flow Control Device

J. Szantyr Lecture No. 21 Aerodynamics of the lifting foils Lifting foils are important parts of many products of contemporary technology.

Standing Waves in a String

Quality Planning for Software Development

Réf GB /d ALTERNATORS LSA Pole - Three phase. Electrical and mechanical data

Introductory Physics PHYS101

Figure Vapor-liquid equilibrium for a binary mixture. The dashed lines show the equilibrium compositions.

Name Date. 5. In each pair, which rational number is greater? Explain how you know.

Provably Secure Camouflaging Strategy for IC Protection

Quantification of the Effects of Turbulence in Wind on the Flutter Stability of Suspension Bridges

Kyoung-Su Im, Z.-C Zhang, and Grant Cook, Jr. LSTC, Livermore, CA USA. 06/14/2016, Dearborn, MI

CALCULATING THE SPEED OF SOUND IN NATURAL GAS USING AGA REPORT NO Walnut Lake Rd th Street Houston TX Garner, IA 50438

Human-Inspired Control of Bipedal Walking Robots

EDEXCEL NATIONALS UNIT 6 MECHANICAL PRINCIPLES and APPLICATIONS. ASSIGNMENT No. 4

CHAPTER 21 AN INCLINED-PLATE WAVE GENERATOR ABSTRACT

Announcements. Lecture 19: Inference for SLR & Transformations. Online quiz 7 - commonly missed questions

Chapter 11 Waves. Waves transport energy without transporting matter. The intensity is the average power per unit area. It is measured in W/m 2.

MODELLING OF WATER FLOW ON SMALL VESSEL S DECK

Fluid Flow. Link. Flow» P 1 P 2 Figure 1. Flow Model

Sizing Pulsation Dampeners Is Critical to Effectiveness

Domain Decomposition Method for 3-Dimensional Simulation of the Piston Cylinder Section of a Hermetic Reciprocating Compressor

AN ISOLATED SMALL WIND TURBINE EMULATOR

Waves. harmonic wave wave equation one dimensional wave equation principle of wave fronts plane waves law of reflection

CHAPTER IV FINITE ELEMENT ANALYSIS OF THE KNEE JOINT WITHOUT A MEDICAL IMPLANT

Quantitative Methods for Economics Tutorial 6. Katherine Eyal

Bioreactors Prof G. K. Suraishkumar Department of Biotechnology Indian Institute of Technology, Madras. Lecture - 15 Bioreactor env.par.

LONG WAVE RUN-UP OVER SUBMERGED REEF AND BREAKWATER

Development of Fish type Robot based on the Analysis of Swimming Motion of Bluefin Tuna Comparison between Tuna-type Fin and Rectangular Fin -

Addition and Subtraction of Rational Expressions

AP Physics 1 Summer Packet Review of Trigonometry used in Physics

67. Sectional normalization and recognization on the PV-Diagram of reciprocating compressor

09 - Choosing /sizing a cylinder and valve

3D CDF MODELING OF SHIP S HEELING MOMENT DUE TO LIQUID SLOSHING IN TANKS A CASE STUDY

Transcription:

Cleveland State University MCE441: Intr. Linear Control Systems Lecture 6: The and Zeros Partial Fraction Expansions Prof. Richter 1 / 12

System and Zeros and Zeros in Decomposition- Obtaining the - in In MCE441, we work with Input/Output systems (black boxes) We have learned how to represent the I/O dynamics of electromechanical systems with O.D.E. s. s carry the same information as the I/O diff. eq. in a more convenient format : The transfer function of a system is the ratio of Laplace transforms of output and input, with all initial conditions set to zero. 2 / 12

and Zeros and Zeros in Decomposition- Obtaining the - in U(s) G(s) = Y (s) U(s) Y (s) Does not provide information about internal system structure (nor does the I/O ODE) Extensively used for studying linear system properties and for design It equals the Laplace transform of the impulse response of a system May be experimentally obtained (impulse hammers) 3 / 12

Galvanometer Model and Zeros and Zeros in Decomposition- Obtaining the - in Let s decouple the galvanometer equations and find the T.F. from applied voltage (V i ) to angular displacement (θ) 4 / 12

and Zeros and Zeros in Decomposition- Obtaining the - in Check that the required transfer function is G(s) = Θ(s) V i (s) = α (Ls 2 + R)(Js 2 + bs + k) + α 2 s 5 / 12

and Zeros and Zeros and Zeros in Decomposition- Obtaining the - in Transfer functions of finite-dimensional (lumped-parameter systems) are always rational functions (ratio of polynomials) of s: G(s) = b 0s m +b 1 s m 1 +...+b m 1 s+b m a 0 s n +a 1 s n 1 +...+a n 1 s+a n : Roots of the denominator Zeros: Roots of the numerator and zeros are the fundamental indicators of dynamic response and stability Factored form: G(S) = K(s+z 1)(s+z 2 )...(s+z m ) (s+p 1 )(s+p 2 )...(s+p n ) 6 / 12

and Zeros in and Zeros and Zeros in Decomposition- Obtaining the - in We limit ourselves to causal systems: system order= n m. example: Find the poles and zeros of the TF: G(s) = s3 +5s 2 +9s+7 s 4 +3s 2 +2s+1 >>num=[1 5 9 7]; >>den=[1 0 3 2 1]; >>roots(num) >>roots(den) Alternatively, to find the gain as well: >>[z,p,k]=tf2zpk(num,den) 7 / 12

and Zeros and Zeros in Decomposition- Real Decomposition- Obtaining the - in Suppose G(s) = n(s) d(s), with G(s) strictly proper (n > m): To every simple real pole at s = p there corresponds a single term of the form a s p where a = lim s p (s p)g(p) To every multiple real pole of multiplicity k 2 there correspond k terms: c 1 s p + c 2 (s p) 2 +...+ c k (s p) k 8 / 12

Decomposition- and Zeros and Zeros in Decomposition- Complex Conjugate Obtaining the - in To every simple pair of complex conjugate poles at s = α±βi there correspond two terms of the form A(s α) (s α) 2 +β 2 + B (s α) 2 +β 2 The quantities a, c i, A and B are called residues. Why by hand? : For complex comjugate case, gives complex residues. We ll skip the case of multiple complex conjugate poles. Rarely found and may be solved using. 9 / 12

Obtaining the - and Zeros and Zeros in Decomposition- Obtaining the - in Find the poles and characterize them as real or complex, single or multiple. Determine the structure of the decomposition. If single real poles are found, find the residues now. Multiply through and equate the coefficients of the numerator to those of the original TF. Solve a linear system of equations. NOTE: If you found some residues for single real poles, you will have more equations than unknowns. Use the extra equations for verification as an alternative to recombining the. 10 / 12

and Zeros and Zeros in Decomposition- Obtaining the - in Decompose the following TFs into partial fractions: 1. G(s) = 2 (s+1)(s+2) 2. G(s) = 2 3. G(s) = (s+1) 3 (s+2) (s+1) (s 2 +2s+4)(s+2) 11 / 12

in and Zeros and Zeros in Decomposition- Obtaining the - in When the TF is proper (but not strictly proper), that is, n m, we need to perform polynomial division to obtain G(s) =constant + strictly proper TF. finds the residues, poles and multiplicities in one step >>[r,p,k]=residue(num,den) NOTE: There exist formulas for finding the residues directly. But overall, the method presented involves less computational work. 12 / 12