d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

Similar documents
Some Basic Information about M-S-D Systems

10. State Space Methods

SOLUTIONS TO ECE 3084

CHAPTER 2 Signals And Spectra

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Signal and System (Chapter 3. Continuous-Time Systems)

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Chapter 2. First Order Scalar Equations

Module 4: Time Response of discrete time systems Lecture Note 2

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Final Spring 2007

Wall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.

Math 333 Problem Set #2 Solution 14 February 2003

INDEX. Transient analysis 1 Initial Conditions 1

LAPLACE TRANSFORM AND TRANSFER FUNCTION

For example, the comb filter generated from. ( ) has a transfer function. e ) has L notches at ω = (2k+1)π/L and L peaks at ω = 2π k/l,

15. Vector Valued Functions

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

Differential Equations

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

Chapter 5 Digital PID control algorithm. Hesheng Wang Department of Automation,SJTU 2016,03

= ( ) ) or a system of differential equations with continuous parametrization (T = R

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

The equation to any straight line can be expressed in the form:

Numerical Dispersion

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

STATE-SPACE MODELLING. A mass balance across the tank gives:

6.2 Transforms of Derivatives and Integrals.

4.6 One Dimensional Kinematics and Integration

Vehicle Arrival Models : Headway

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0.

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

Circuit Variables. AP 1.1 Use a product of ratios to convert two-thirds the speed of light from meters per second to miles per second: 1 ft 12 in

Y 0.4Y 0.45Y Y to a proper ARMA specification.

Chapter 8 The Complete Response of RL and RC Circuits

Linear Response Theory: The connection between QFT and experiments

2.9 Modeling: Electric Circuits

This document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC

Math Week 15: Section 7.4, mass-spring systems. These are notes for Monday. There will also be course review notes for Tuesday, posted later.

Two Coupled Oscillators / Normal Modes

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves

Stability and Bifurcation in a Neural Network Model with Two Delays

MATH 31B: MIDTERM 2 REVIEW. x 2 e x2 2x dx = 1. ue u du 2. x 2 e x2 e x2] + C 2. dx = x ln(x) 2 2. ln x dx = x ln x x + C. 2, or dx = 2u du.

2. Nonlinear Conservation Law Equations

Continuous Time Linear Time Invariant (LTI) Systems. Dr. Ali Hussein Muqaibel. Introduction

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

Sliding Mode Controller for Unstable Systems

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

Morning Time: 1 hour 30 minutes Additional materials (enclosed):

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

THE BERNOULLI NUMBERS. t k. = lim. = lim = 1, d t B 1 = lim. 1+e t te t = lim t 0 (e t 1) 2. = lim = 1 2.

EXERCISES FOR SECTION 1.5

Solutions to Assignment 1

The motions of the celt on a horizontal plane with viscous friction

3.6 Derivatives as Rates of Change

Problem Set #1. i z. the complex propagation constant. For the characteristic impedance:

EECE 301 Signals & Systems Prof. Mark Fowler

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term

Lectures 29 and 30 BIQUADRATICS AND STATE SPACE OP AMP REALIZATIONS. I. Introduction

4. Advanced Stability Theory

EE100 Lab 3 Experiment Guide: RC Circuits

Testing for a Single Factor Model in the Multivariate State Space Framework

Math Final Exam Solutions

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

Laplace Transform and its Relation to Fourier Transform

Time Domain Transfer Function of the Induction Motor

Homework 6 AERE331 Spring 2019 Due 4/24(W) Name Sec. 1 / 2

Chapter 4. Truncation Errors

Traveling Waves. Chapter Introduction

Sub Module 2.6. Measurement of transient temperature

CHAPTER 12 DIRECT CURRENT CIRCUITS

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Solutions Problem Set 3 Macro II (14.452)

- If one knows that a magnetic field has a symmetry, one may calculate the magnitude of B by use of Ampere s law: The integral of scalar product

5.1 - Logarithms and Their Properties

Economics 8105 Macroeconomic Theory Recitation 6

Basic Circuit Elements Professor J R Lucas November 2001

Ordinary dierential equations

RC, RL and RLC circuits

Physics 180A Fall 2008 Test points. Provide the best answer to the following questions and problems. Watch your sig figs.

2. The basics of differential equations 1 AGEC 642 Fall 2018

ADDITIONAL PROBLEMS (a) Find the Fourier transform of the half-cosine pulse shown in Fig. 2.40(a). Additional Problems 91

8. Basic RL and RC Circuits

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson

EE 435. Lecture 31. Absolute and Relative Accuracy DAC Design. The String DAC

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

EECE 301 Signals & Systems Prof. Mark Fowler

The complex Fourier series has an important limiting form when the period approaches infinity, i.e., T 0. 0 since it is proportional to 1/L, but

INVERSE RESPONSE COMPENSATION BY ESTIMATING PARAMETERS OF A PROCESS COMPRISING OF TWO FIRST ORDER SYSTEMS

Echocardiography Project and Finite Fourier Series

Unsteady Flow Problems

Chapter 3 Boundary Value Problem

The expectation value of the field operator.

HOMEWORK # 2: MATH 211, SPRING Note: This is the last solution set where I will describe the MATLAB I used to make my pictures.

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

Transcription:

and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or muliplied by a posiive number in order o simplify he subsequen numerical calculaion wihou alering he sabiliy conclusion. Rouh s sabiliy crierion saes ha he number of roos of Equaion ( 6) wih posiive real pars is equal o he number of changes in sign of he coefficiens of he firs column of he array. I should be noed ha he exac values of he erms in he firs column need no be known; insead, only he signs are needed. The necessary and sufficien condiion ha all roos of Equaion ( 6) lie in he lef-half s plane is ha all he coefficiens of Equaion ( 6) be posiive and all erms in he firs column of he array have posiive signs. EXAMPLE Le us apply Rouh s sabiliy crierion o he following hird-order polynomial: a + a + a s + a = where all he coefficiens are posiive numbers. The array of coefficiens becomes a a a a - a a a a a a The condiion ha all roos have negaive real pars is given by a a 7 a a EXAMPLE Consider he following polynomial: + + + 4s + = Le us follow he procedure jus presened and consruc he array of coefficiens. (The firs wo rows can be obained direcly from he given polynomial. The remaining erms are 4 Chaper / Transien and Seady-Sae Response Analyses

obained from hese. If any coefficiens are missing, hey may be replaced by zeros in he array.) - 6 4 6 6-4 The second row is divided by. In his example, he number of changes in sign of he coefficiens in he firs column i. This means ha here are wo roos wih posiive real pars. Noe ha he resul is unchanged when he coefficiens of any row are muliplied or divided by a posiive number in order o simplify he compuaion. Special Cases. If a firs-column erm in any row is zero, bu he remaining erms are no zero or here is no remaining erm, hen he zero erm is replaced by a very small posiive number and he res of he array is evaluaed. For example, consider he following equaion: + + s + = ( 6) The array of coefficiens is L If he sign of he coefficien above he zero ( ) is he same as ha below i, i indicaes ha here are a pair of imaginary roos. Acually, Equaion ( 6) has wo roos a s=; j. If, however, he sign of he coefficien above he zero ( ) is opposie ha below i, i indicaes ha here is one sign change. For example, for he equaion he array of coefficiens is - s + = (s - ) (s + ) = One sign change: L - One sign change: - - There are wo sign changes of he coefficiens in he firs column. So here are wo roos in he righ-half s plane. This agrees wih he correc resul indicaed by he facored form of he polynomial equaion. Secion 6 / Rouh s Sabiliy Crierion

If all he coefficiens in any derived row are zero, i indicaes ha here are roos of equal magniude lying radially opposie in he s plane ha is, wo real roos wih equal magniudes and opposie signs and/or wo conjugae imaginary roos. In such a case, he evaluaion of he res of he array can be coninued by forming an auxiliary polynomial wih he coefficiens of he las row and by using he coefficiens of he derivaive of his polynomial in he nex row. Such roos wih equal magniudes and lying radially opposie in he s plane can be found by solving he auxiliary polynomial, which is always even. For a n-degree auxiliary polynomial, here are n pairs of equal and opposie roos. For example, consider he following equaion: The array of coefficiens is The erms in he row are all zero. (Noe ha such a case occurs only in an oddnumbered row.) The auxiliary polynomial is hen formed from he coefficiens of he row. The auxiliary polynomial P(s) is which indicaes ha here are wo pairs of roos of equal magniude and opposie sign (ha is, wo real roos wih he same magniude bu opposie signs or wo complexconjugae roos on he imaginary axis). These pairs are obained by solving he auxiliary polynomial equaion P(s)=. The derivaive of P(s) wih respec o s is The erms in he row are replaced by he coefficiens of he las equaion ha is, 8 and 96. The array of coefficiens hen becomes s s s + + 4 + 48 - s - = 4 48 8 4.7 - - - d Auxiliary polynomial P(s) P(s) = + 48-4 48 96 - dp (s) ds - - = 8 + 96s d Coefficiens of dp (s) ds We see ha here is one change in sign in he firs column of he new array.thus, he original equaion has one roo wih a posiive real par. By solving for roos of he auxiliary polynomial equaion, we obain or + 48 - = =, =- s =;, s =;j 6 Chaper / Transien and Seady-Sae Response Analyses

These wo pairs of roos of P(s) are a par of he roos of he original equaion. As a maer of fac, he original equaion can be wrien in facored form as follows: Clearly, he original equaion has one roo wih a posiive real par. Relaive Sabiliy Analysis. Rouh s sabiliy crierion provides he answer o he quesion of absolue sabiliy. This, in many pracical cases, is no sufficien. We usually require informaion abou he relaive sabiliy of he sysem. A useful approach for examining relaive sabiliy is o shif he s-plane axis and apply Rouh s sabiliy crierion. Tha is, we subsiue ino he characerisic equaion of he sysem, wrie he polynomial in erms of ŝ; and apply Rouh s sabiliy crierion o he new polynomial in ŝ. The number of changes of sign in he firs column of he array developed for he polynomial in ŝ is equal o he number of roos ha are locaed o he righ of he verical line s= s.thus, his es reveals he number of roos ha lie o he righ of he verical line s= s. Applicaion of Rouh s Sabiliy Crierion o Conrol-Sysem Analysis. Rouh s sabiliy crierion is of limied usefulness in linear conrol-sysem analysis, mainly because i does no sugges how o improve relaive sabiliy or how o sabilize an unsable sysem. I is possible, however, o deermine he effecs of changing one or wo parameers of a sysem by examining he values ha cause insabiliy. In he following, we shall consider he problem of deermining he sabiliy range of a parameer value. Consider he sysem shown in Figure. Le us deermine he range of for sabiliy. The closed-loop ransfer funcion is The characerisic equaion is The array of coefficiens becomes (s + )(s - )(s + j)(s - j)(s + ) = s = ŝ - s (s = consan) C(s) = + + + s + = 7-9 7 sa + s + B(s + ) + + s( + s + ) (s + ) C(s) Figure Conrol sysem. Secion 6 / Rouh s Sabiliy Crierion 7

For sabiliy, mus be posiive, and all coefficiens in he firs column mus be posiive. Therefore, 4 9 7 7 When = 4 9, he sysem becomes oscillaory and, mahemaically, he oscillaion is susained a consan ampliude. Noe ha he ranges of design parameers ha lead o sabiliy may be deermined by use of Rouh s sabiliy crierion. 7 EFFECTS OF INTEGRAL AND DERIVATIVE CONTROL ACTIONS ON SYSTEM PERFORMANCE In his secion, we shall invesigae he effecs of inegral and derivaive conrol acions on he sysem performance. Here we shall consider only simple sysems, so ha he effecs of inegral and derivaive conrol acions on sysem performance can be clearly seen. Inegral Conrol Acion. In he proporional conrol of a plan whose ransfer funcion does no possess an inegraor s, here is a seady-sae error, or offse, in he response o a sep inpu. Such an offse can be eliminaed if he inegral conrol acion is included in he conroller. In he inegral conrol of a plan, he conrol signal he oupu signal from he conroller a any insan is he area under he acuaing-error-signal curve up o ha insan. The conrol signal u() can have a nonzero value when he acuaing error signal e() is zero, as shown in Figure 6(a).This is impossible in he case of he proporional conroller, since a nonzero conrol signal requires a nonzero acuaing error signal. (A nonzero acuaing error signal a seady sae means ha here is an offse.) Figure 6(b) shows he curve e() versus and he corresponding curve u() versus when he conroller is of he proporional ype. Noe ha inegral conrol acion, while removing offse or seady-sae error, may lead o oscillaory response of slowly decreasing ampliude or even increasing ampliude, boh of which are usually undesirable. Figure 6 (a) Plos of e() and u() curves showing nonzero conrol signal when he acuaing error signal is zero (inegral conrol); (b) plos of e() and u() curves showing zero conrol signal when he acuaing error signal is zero (proporional conrol). e() u() (a) e() u() (b) 8 Chaper / Transien and Seady-Sae Response Analyses

Figure 7 Proporional conrol sysem. E(s) C(s) + Ts + Proporional Plan conroller Proporional Conrol of Sysems. We shall show ha he proporional conrol of a sysem wihou an inegraor will resul in a seady-sae error wih a sep inpu.we shall hen show ha such an error can be eliminaed if inegral conrol acion is included in he conroller. Consider he sysem shown in Figure 7. Le us obain he seady-sae error in he uni-sep response of he sysem. Define Since E(s) he error E(s) is given by For he uni-sep inpu =/s, we have The seady-sae error is = - C(s) E(s) = + G(s) = + Ts + E(s) = G(s) = Ts + = - C(s) = + G(s) Ts + Ts + + s Ts + e ss = lim e() = lim se(s) = lim Sq s S s S Ts + + = + Such a sysem wihou an inegraor in he feedforward pah always has a seady-sae error in he sep response. Such a seady-sae error is called an offse. Figure 8 shows he uni-sep response and he offse. c() Offse Figure 8 Uni-sep response and offse. Secion 7 / Effecs of Inegral and Derivaive Conrol Acions on Sysem Performance 9