Transfer Function Analysis

Similar documents
School of Mechanical Engineering Purdue University. ME375 Transfer Functions - 1

Dynamic Response of Linear Systems

Chapter 4 : Laplace Transform

School of Mechanical Engineering Purdue University. ME375 Frequency Response - 1

(s)h(s) = K( s + 8 ) = 5 and one finite zero is located at z 1

University of California at Berkeley College of Engineering Department of Electrical Engineering and Computer Sciences

Class 07 Time domain analysis Part II 2 nd order systems

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

FIR Filters. Lecture #7 Chapter 5. BME 310 Biomedical Computing - J.Schesser

The state space model needs 5 parameters, so it is not as convenient to use in this control study.

Chapter 7: The z-transform. Chih-Wei Liu

Solution of Linear Constant-Coefficient Difference Equations

Definition of z-transform.

Generalizing the DTFT. The z Transform. Complex Exponential Excitation. The Transfer Function. Systems Described by Difference Equations

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

Bernoulli Polynomials Talks given at LSBU, October and November 2015 Tony Forbes

Bode Diagrams School of Mechanical Engineering ME375 Frequency Response - 29 Purdue University Example Ex:

Z - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.

Chapter 1. Complex Numbers. Dr. Pulak Sahoo

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

M2.The Z-Transform and its Properties

Answer: 1(A); 2(C); 3(A); 4(D); 5(B); 6(A); 7(C); 8(C); 9(A); 10(A); 11(A); 12(C); 13(C)

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

DIGITAL SIGNAL PROCESSING LECTURE 3

Digital Signal Processing

Lecture 20 - Wave Propagation Response

Chapter 2. Asymptotic Notation

6.003 Homework #3 Solutions

Lab(8) controller design using root locus

5.6 Binomial Multi-section Matching Transformer

Vasyl Moisyshyn*, Bogdan Borysevych*, Oleg Vytyaz*, Yuriy Gavryliv*

Chapter 3. z-transform

Time-Domain Representations of LTI Systems

School of Mechanical Engineering Purdue University. ME375 Dynamic Response - 1

Discrete-Time Systems, LTI Systems, and Discrete-Time Convolution

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science. BACKGROUND EXAM September 30, 2004.

: Transforms and Partial Differential Equations

Solutions - Homework # 1

19.1 The dictionary problem

2D DSP Basics: Systems Stability, 2D Sampling

CEMTool Tutorial. The z-transform

Exponential Moving Average Pieter P

Principle Of Superposition

Digital signal processing: Lecture 5. z-transformation - I. Produced by Qiangfu Zhao (Since 1995), All rights reserved

Define a Markov chain on {1,..., 6} with transition probability matrix P =

Chapter 7 z-transform

AH Checklist (Unit 3) AH Checklist (Unit 3) Matrices

ECE 308 Discrete-Time Signals and Systems

Lebesgue Constant Minimizing Bivariate Barycentric Rational Interpolation

Solutions of Chapter 5 Part 1/2

EECE 301 Signals & Systems

School of Mechanical Engineering Purdue University

f(w) w z =R z a 0 a n a nz n Liouville s theorem, we see that Q is constant, which implies that P is constant, which is a contradiction.

The Non-homogeneous Diffusion Equation

The Performance of Feedback Control Systems

MA Lesson 26 Notes Graphs of Rational Functions (Asymptotes) Limits at infinity

Acoustic Field inside a Rigid Cylinder with a Point Source

State space systems analysis

PARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES

Engineering Mechanics Dynamics & Vibrations. Engineering Mechanics Dynamics & Vibrations Plane Motion of a Rigid Body: Equations of Motion

ELEG 4603/5173L Digital Signal Processing Ch. 1 Discrete-Time Signals and Systems

ME 375 FINAL EXAM Friday, May 6, 2005

EECE 301 Signals & Systems Prof. Mark Fowler

Sinusoidal Steady-state Analysis

ELEG3503 Introduction to Digital Signal Processing

1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y

Discrete population models

ECE Spring Prof. David R. Jackson ECE Dept. Notes 20

Most text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t

Chapter 10 z Transform

The Binomial Multi-Section Transformer

Computing the output response of LTI Systems.

6.003: Signals and Systems. Feedback, Poles, and Fundamental Modes

A Generalization of Ince s Equation

The z-transform can be used to obtain compact transform-domain representations of signals and systems. It

Summer MA Lesson 13 Section 1.6, Section 1.7 (part 1)

(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is

Course Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B.

SOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.

Difference Equation Construction (1) ENGG 1203 Tutorial. Difference Equation Construction (2) Grow, baby, grow (1)

Signals and Systems. Problem Set: From Continuous-Time to Discrete-Time

The Z-Transform. (t-t 0 ) Figure 1: Simplified graph of an impulse function. For an impulse, it can be shown that (1)

Solution of EECS 315 Final Examination F09

Exam. Notes: A single A4 sheet of paper (double sided; hand-written or computer typed)

6. Uniform distribution mod 1

Multinomial. Notations. Primary definition. Specific values. Traditional name. Traditional notation. Mathematica StandardForm notation

Unit 4: Polynomial and Rational Functions

GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL OF ELECTRICAL AND COMPUTER ENGINEERING

5.6 Binomial Multi-section Matching Transformer

COMM 602: Digital Signal Processing

Discrete-Time Signals and Systems. Discrete-Time Signals and Systems. Signal Symmetry. Elementary Discrete-Time Signals.

Application of Homotopy Analysis Method for Solving various types of Problems of Ordinary Differential Equations

Olli Simula T / Chapter 1 3. Olli Simula T / Chapter 1 5

1the 1it is said to be overdamped. When 1, the roots of

Signal Processing. Lecture 02: Discrete Time Signals and Systems. Ahmet Taha Koru, Ph. D. Yildiz Technical University.

DETERMINATION OF NATURAL FREQUENCY AND DAMPING RATIO

Chapter 4 Postulates & General Principles of Quantum Mechanics

j=1 dz Res(f, z j ) = 1 d k 1 dz k 1 (z c)k f(z) Res(f, c) = lim z c (k 1)! Res g, c = f(c) g (c)

Generating Functions. 1 Operations on generating functions

e to approximate (using 4

Transcription:

Trasfer Fuctio Aalysis Free & Forced Resposes Trasfer Fuctio Syste Stability ME375 Trasfer Fuctios - Free & Forced Resposes Ex: Let s s look at a stable first order syste: τ y + y = Ku Take LT of the I/O odel ad reeber to keep tracks of the ICs: L[ τ y + y] = L[ Ku] τ ( ) + = K Rearrage ters s.t. the output Y(s) ) ters are o oe side ad the iput U(s) ) ad IC ters are o the other: Ys () = Us () + y() ( ) ( ) ( ) Factor out the output side of the equatio: Y ( s) = U ( s) + y() ME375 Trasfer Fuctios - 2

Free & Forced Resposes Free Respose (u(t)( ) = & ozero ICs) The respose of a syste to zero iput ad ozero iitial coditios. Ca be obtaied by Let u(t) ) = ad use LT ad ILT to solve for the free respose. Forced Respose (zero( ICs & ozero u(t)) The respose of a syste to ozero iput ad zero iitial coditios. Ca be obtaied by Assue zero ICs ad use LT ad ILT to solve for the forced respose (replace differetiatio with s i the I/O ODE). ME375 Trasfer Fuctios - 3 I Class Exercise Fid the free ad forced resposes of the followig I/O odel: y + 4 y y + 5y = 2u u + u ME375 Trasfer Fuctios - 4 2

Trasfer Fuctio Give a geeral th order syste odel: ( ) ( ) ( ) ( ) a y + a y + + ay + a y = bu + b u + + bu + bu The forced respose (zero ICs) of the syste due to iput u(t) ) is: Takig the LT of the ODE: ( ) L y = s Y( s) ( WHY? ) asys () + a s Ys () + + asys () + ays () = bsus ( ) + b s Us ( ) + + bsus ( ) + bus ( ) () as + a s + + as+ a Ys = bs + b s + + bs+ b Us () Ds ( ) Ns ( ) bs + b s + + bs + b Ns () Ys () = Us () = Us () = Gs () Us () as + a s + + as + a Ds () Gs () Trasfer Fuctio ME375 Trasfer Fuctios - 5 Exaples () Recall the first order syste: τ y + y = K u Fid the trasfer fuctio of the syste. Takig LT of the ODE: (2) For the followig 2d order syste: y + 2 ζ ω y + ω y = Kω u Fid the trasfer fuctio of the syste. Takig LT of the ODE: 2 2 ME375 Trasfer Fuctios - 6 3

Trasfer Fuctio Give a geeral th order syste: ay + a y + + ay + ay= bu + b u + + bu + bu ( ) ( ) ( ) ( ) The trasfer fuctio of the syste is: bs + b s + + bs+ b Gs () as + a s + + as+ a The trasfer fuctio ca be iterpreted as: u(t) Iput Differetial Equatio y(t) Output U(s) Iput G(s) Y(s) Output Tie Doai s -Doai ME375 Trasfer Fuctios - 7 Poles ad Zeros Give a trasfer fuctio (TF) of a syste: bs + b s + + bs + b N() s Gs () = as + a s + + as + a Ds () Poles Zeros The roots of the deoiator of the The roots of the uerator of the TF. TF, i.e. the roots of the characteristic equatio. Ns () = bs + b s + + bs + b Ds () = as + a s + + as + a = b( s z)( s z2) ( s z) = = a( s p)( s p2) ( s p) = z, z2,, z : zeros of the TF p, p,, p : poles of the TF 2 bs + b s + + bs+ b Ns ( ) Gs () = = as + a s + + as+ a Ds () ME375 Trasfer Fuctios - 8 4

Static Gai Static Gai ( G() ) The value of the trasfer fuctio whe s =. If bs + b s + + bs + b N() s Gs () = as + a s + + as + a Ds () N() b KS = G() = = D() a The static gai K S ca be iterpreted as the steady state value of the uit step respose. Ex: For a secod order syste: 2 2 y+ 2ζ ωy + ω y = Ksω u Fid the trasfer fuctio ad the static gai. Ex: Fid the steady state value of the syste y + 3 y + 5y + 7y = u + 2u + u to a step iput of agitude 2. ME375 Trasfer Fuctios - 9 I Class Exercise Give a I/O odel of a 2d order syste: 5 y + 2 y + 4 y = 6 u () Fid the trasfer fuctio of the syste (2) Fid the poles ad zeros of the syste (3) What is the syste forced respose Y(s) ) to a uit step iput u(t) ) =? (4) As tie goes to ifiity, what is the steady state value of the uit step respose? ME375 Trasfer Fuctios - 5

A Closer Look at Free Respose Give a geeral th order syste odel: ( ) ( ) ( ) ( ) ay + a y + + ay + ay = bu + b u + + bu + bu The free respose (zero iput) of the syste due to ICs is: ( ) ( ) Takig the LT of the Hoogeeous ODE: ay + a y + + ay + ay= ( ) ( ) L y L y ( ) 2 ( 2) a s Y( s) s y() y () + a s Y( s) s y() y () [ ] + + a sy( s) y() + a Y( s) = L y L y as + a s + + as+ a YFree ( s) 2 ( ) = ( as + a s + + a) y() + + y () Y ( s) = F( s) F( s) A Polyoial of s That depeds o ICs Free as + a s + + as+ a ME375 Trasfer Fuctios - A Closer Look at Free Respose Ex: : Fid the free respose of the followig syste: 5 y + 3y + 5y = 2 u + u Ex: : Perfor partial fractio expasio (PFE) of the above free respose whe: y ( ) = ad y( ) = ME375 Trasfer Fuctios - 2 6

Free Respose ad Pole Positio The free respose of a syste ca be represeted by: Fs () Fs () YFree() s = = as + a s + + as + a a( s p)( s p2) ( s p) A A2 A = + + + s p s p2 s p Assue p p p i.e. distict poles R S T 2 p t p2 t y ( t) = Y ( s) = A e + A e + + A e p p i i is real = σ + jβ Free R S T i R S T pi < pi = p > σ < σ = σ > p t L Free 2 Ig. Real ME375 Trasfer Fuctios - 3 Coplete Respose U(s) N() s Y(s) Gs () Iput D() s Output ( ) ICs: y(), y (),, y () Coplete Respose Y ( s) = Y ( s) + Y ( s) = G ( s) U ( s) + Y ( s) Forced Free Free = U( s) + Q: What part of the syste affects both the free ad forced respose? Q: If U(s) U ) = ad there are o-zero ICs, what will guaratee that y(t) y? ME375 Trasfer Fuctios - 4 7

Stability Stability Cocept Describes the ability of a syste to stay at its equilibriu positio (for liear systes: all state variables = or y(t) ) = ) i the absece of ay iputs. A liear tie ivariat (LTI) syste is stable if ad oly if (iff) its free respose coverges to zero. Ex: Pedulu Ball o curved surface ME375 Trasfer Fuctios - 5 Stability of LTI Systes Stability Criterio for LTI Systes ( ) ( ) ( ) ( ) ay + a y + + ay + ay = bu + b u + + bu + bu Stable all poles of Ds ( ) = as + a s + + as + a lie i the left-half coplex plae Coets o LTI Stability Stability of a LTI syste does ot deped o the iput. (why?) For st ad 2d order systes, stability is guarateed if all the e coefficiets of the characteristic polyoial are positive. D( s) = as + a : Stable ai > i or ai < i 2 Ds ( ) = as + as+ a : Stable a> i or a< i 2 Characteristic Polyoial Effect of Poles ad Zeros o Stability Stability of a syste depeds oly o its poles. Zeros do ot affect syste stability. Zeros affect the trasiet respose of the syste. i i LHP ME375 Trasfer Fuctios - 6 8

I Class Exercises () Fid the trasfer fuctio of the followig I/O equatio: y 2 y 5y = 2u + u (2) Deterie the syste s s stability. (3) Plot the poles ad zeros of the syste o the coplex plae. () Fid the trasfer fuctio of the followig I/O equatio: y + y + 6y = u 3u + 4u (2) Deterie the syste s s stability. (3) Plot the poles ad zeros of the syste o the coplex plae. ME375 Trasfer Fuctios - 7 Exaple (Iverted) Pedulu () Derive a atheatical odel for a pedulu. (2) Fid the equilibriu positios. (3) Discuss the stability of the equilibriu positios. ME375 Trasfer Fuctios - 8 9

Exaple (Iverted Pedulu) ME375 Trasfer Fuctios - 9