๐Ÿ•ท๏ธ Crawler Inspector

URL Lookup

Direct Parameter Lookup

Raw Queries and Responses

1. Shard Calculation

Query:
Response:
Calculated Shard: 0 (from laksa048)

2. Crawled Status Check

Query:
Response:

3. Robots.txt Check

Query:
Response:

4. Spam/Ban Check

Query:
Response:

5. Seen Status Check

โ„น๏ธ Skipped - page is already crawled

๐Ÿ“„
INDEXABLE
โœ…
CRAWLED
4 days ago
๐Ÿค–
ROBOTS ALLOWED

Page Info Filters

FilterStatusConditionDetails
HTTP statusPASSdownload_http_code = 200HTTP 200
Age cutoffPASSdownload_stamp > now() - 6 MONTH0.2 months ago
History dropPASSisNull(history_drop_reason)No drop reason
Spam/banPASSfh_dont_index != 1 AND ml_spam_score = 0ml_spam_score=0
CanonicalPASSmeta_canonical IS NULL OR = '' OR = src_unparsedNot set

Page Details

PropertyValue
URLhttps://pressbooks.library.torontomu.ca/controlsystems/chapter/1-4-laplace-transforms/
Last Crawled2026-04-04 02:15:48 (4 days ago)
First Indexed2023-03-02 19:22:30 (3 years ago)
HTTP Status Code200
Meta Title1.4 Laplace Transforms โ€“ Introduction to Control Systems
Meta Descriptionnull
Meta Canonicalnull
Boilerpipe Text
Self-Study: Review your ELE532 Notes and other resources. You can also refer to the review material on the course website. 1.4.1 Definitions F ( s ) = L [ f ( t ) ] = โˆซ 0 + โˆž f ( t ) e โˆ’ s t d t f ( t ) = L โˆ’ 1 [ F ( s ) ] = 1 2 ฯ€ j โˆซ ฯƒ โˆ’ j โˆž ฯƒ + j โˆž F ( s ) e s t d s 1.4.1.1 Final Value Theorem f s s = lim t โ†’ โˆž f ( t ) = lim s โ†’ 0 s F ( s ) 1.4.1.2 Initial Value Theorem f 0 = lim t โ†’ 0 + f ( t ) = lim s โ†’ โˆž s F ( s ) 1.4.1.3ย  Properties of Laplace transforms Table 1โ€‘1 Properties of Laplace Transform F ( s ) e โˆ’ T s f ( t โˆ’ T ) โ‹… 1 ( t ) F ( s + a ) f ( t ) e โˆ’ a t โ‹… 1 ( t ) s F ( s ) โˆ’ f ( 0 + ) d f ( t ) d t S 2 F ( s ) โˆ’ s f ( 0 + ) โˆ’ d f ( 0 + ) d t d 2 f ( t ) d t 2 1 s F ( s ) โˆซ 0 + + โˆž f ( t ) d t F 1 ( s ) โ‹… F 2 ( s ) f 1 ( t ) โˆ— f 2 ( t ) 1.4.2 Solving for System Response Parametric models cannot be developed without math. Laws of physics describe dynamic Linear Time-Invariant (LTI) systems using ordinary differential equations. To simplify their analysis, Laplace transform is used. Consider a certain LTI (Linear Time-Invariant), SISO (Single Input Single Output) system: Let the input โ€“ output relationship for the system be described by the following n th order differential equation: d n y d t n + a n โˆ’ 1 d n โˆ’ 1 y d t n โˆ’ 1 + . . . + a 1 d y d t + a 0 y = b m d m u d t m + b m โˆ’ 1 d m โˆ’ 1 u d t m โˆ’ 1 + . . . + b 1 d u d t + b 0 u The equation parameters relate to physical aspects of the system. The time domain description of systems is not convenient for quick paper-and-pencil speculations. To simplify math, Classical Control uses a Laplace Transform system description, which converts the differential equations into their algebraic equivalents in the s-domain. The solution for y(t) can then be found using inverse Laplace transformation to Y(s). 1.4.3 Two Transfer Functions Models: TF and ZPK In the transform domain, the input-output relationship of the system is defined by a transfer function G(s), defined as a ratio of the Laplace transform of the system output signal y(t), to the Laplace transform of the system input signal u(t), with any initial conditions in the system set to zero. The system transfer function G(s) can be thought of as a dynamic gain of the system: Block diagrams are used to graphically represent systems and their components, as shown above. In order to find G(s), a Laplace transform of the system differential equation in Equation 1โ€‘4 is taken: s n Y ( s ) + a n โˆ’ 1 s n โˆ’ 1 Y ( s ) + . . . + a 1 s Y ( s ) + a 0 Y ( s ) = b m s m U ( s ) + b m โˆ’ 1 s m โˆ’ 1 U ( s ) + . . . + b 1 s U ( s ) + b 0 U ( s ) G ( s ) = Y ( s ) U ( s ) = b m s m + b m โˆ’ 1 s m โˆ’ 1 + . . . + b 1 s + b 0 s n + a n โˆ’ 1 s n โˆ’ 1 + . . . + a 1 s + a 0 Equation 1โ€‘5 Transfer functions are ratios of polynomials written in terms of the s-operator. The resulting function in Equation 1โ€‘5 is a ratio of two polynomials, N(s) and D(s): G ( s ) = N ( s ) D ( s ) = b m s m + b m โˆ’ 1 s m โˆ’ 1 + . . . + b 1 s + b 0 s n + a n โˆ’ 1 s n โˆ’ 1 + . . . + a 1 s + a 0 Equation 1โ€‘6 Roots of the numerator polynomial of G(s) in Equation 1โ€‘6 are called system zeros, z i , and roots of the denominator polynomial are called system poles, p i . Note that in MATLAB a transfer function object in a polynomial form can be created by using the โ€œtfโ€ command. For example, consider the following transfer function in the polynomial ratio form: G ( s ) = Y ( s ) U ( s ) = 2 s + 20 s 2 + 4 s + 3 The same transfer function G ( s ) can be represented in the so-called ZPK form (factorized form): G ( s ) = K โˆ i m ( s โˆ’ z i ) โˆ j n ( s โˆ’ p j ) Equation 1โ€‘7 K is a multiplier. It is important to see the difference between K and K dc, which denotes the DC gain of the system (i.e. s=0): K d c = G ( 0 ) = N ( 0 ) D ( 0 ) = b 0 a 0 = K โˆ i m ( โˆ’ z i ) โˆ j n ( โˆ’ p j ) Equation 1โ€‘8 Our transfer function can be factorized and the multiplier gain K is equal to 2: G ( s ) = Y ( s ) U ( s ) = 2 ( s + 10 ) ( s + 1 ) ( s + 3 ) The DC gain of our transfer function is: K D C = G ( 0 ) = 2 ( 10 ) ( 1 ) ( 3 ) = 6.667 In MATLAB the transfer function can be shown in a factorized form by using the โ€œzpkโ€ command, and the DC gain can be found using โ€œdcgainโ€ command: In MATLAB the โ€œzpkโ€ command can be used to create a transfer function object in a ZPK form, and โ€œtfโ€ command to convert it to a polynomial form. Locations of the system poles and zeros can be presented graphically as the so-called Pole-Zero Map. 1.4.4 Partial Fractions Technique If a certain control system is described by a transfer function G(s), the system response can be found as Y ( s ) = U ( s ) โ‹… G ( s ) . Since Laplace Transform Tables do not provide exhaustive solutions, a technique of a Partial Fractions Expansion is used to find inverse Laplace Transforms for various time functions โ€“ see a table of basic Laplace โ€“ Time Domain Function pair shown in Table 1โ€‘2. 1.4.4.1 Residues โ€“ Distinct Roots Case Y ( s ) = N ( s ) โˆ i n ( s โˆ’ p i ) = N ( s ) D ( s ) Y ( s ) = K 1 s โˆ’ p 1 + K 2 s โˆ’ p 2 + . . . + K n s โˆ’ p n Equation 1โ€‘9 K 1 = N ( s ) ( s โˆ’ p 1 ) D ( s ) | s = p 1 K 2 = N ( s ) ( s โˆ’ p 2 ) D ( s ) | s = p 2 โ‹ฎ K n = N ( s ) ( s โˆ’ p n ) D ( s ) | s = p n ย  Equation 1โ€‘10 1.4.4.2 Residues โ€“ Multiple Roots Case This is the case for Laplace Transforms with multiple powers of some roots. Assume multiplicity of m for root r 1 : Y ( s ) = N ( s ) โˆ i n ( s โˆ’ p i ) = N ( s ) D ( s ) Y ( s ) = K 1 s โˆ’ r 1 + K 2 ( s โˆ’ r 1 ) 2 + . . . + K m ( s โˆ’ r 1 ) m + . . . + K n s โˆ’ p n Residues for distinct roots calculated as before: K m + 1 = N ( s ) ( s โˆ’ p 1 ) D ( s ) | s = p m + 1 โ‹ฎ K n = N ( s ) ( s โˆ’ p n ) D ( s ) | s = p n To calculate residues for a multiple root, multiply both sides of Equation 1โ€‘12 by ( s โˆ’ r 1 ) m and substitute the value of the root r 1 : N ( s ) D ( s ) = K 1 s โˆ’ r 1 + K 2 ( s โˆ’ r 1 ) 2 + . . . + K m ( s โˆ’ r 1 ) m + . . . + K n s โˆ’ p n N ( s ) ( s โˆ’ r 1 ) m D ( s ) = K 1 ( s โˆ’ r 1 ) m s โˆ’ r 1 + K 2 ( s โˆ’ r 1 ) m ( s โˆ’ r 1 ) 2 + . . . + K m + . . . + K n ( s โˆ’ r 1 ) m s โˆ’ p n N ( s ) ( s โˆ’ r 1 ) m D ( S ) | s = r 1 = K m To calculate the residues for the second multiplicity of the root r 1 : d d s ( N ( s ) ( s โˆ’ r 1 ) m D ( s ) ) = d d s ( K 1 ( s โˆ’ r 1 ) m s โˆ’ r 1 + K 2 ( s โˆ’ r 1 ) m ( s โˆ’ r 1 ) 2 + . . . + K m + . . . + K n ( s โˆ’ r 1 ) m s โˆ’ p n ) d d s ( N ( s ) ( s โˆ’ r 1 ) m D ( s ) ) | s = r 1 = 0 + . . . + K m โˆ’ 1 + 0 + . . . + 0 Laplace Transform Time Domain Function 1 ฯƒ ( t ) 1 s 1 ( t ) 1 s 2 t โ‹… 1 ( t ) 1 s k + 1 t k k ! 1 ( t ) 1 s + a e โˆ’ a t โ‹… 1 ( t ) 1 ( s + a ) 2 t e โˆ’ a t โ‹… 1 ( t ) a s ( s + a ) ( 1 โˆ’ e โˆ’ a t ) โ‹… ( t ) a s 2 + a 2 s i n ( a t ) โ‹… 1 ( t ) s s 2 + a 2 c o s ( a t ) โ‹… 1 ( t ) s + a ( s + a ) 2 + b 2 e โˆ’ a t โ‹… c o s ( b t ) โ‹… 1 ( t ) b ( s + a ) 2 + b 2 e โˆ’ a t โ‹… s i n ( b t ) โ‹… 1 ( t ) a 2 + b 2 s [ ( s + a ) 2 + b 2 ] ( 1 โˆ’ e โˆ’ a t โ‹… ( c o s ( b t ) + a b โ‹… s i n ( b t ) ) ) โ‹… 1 ( t ) ฯ‰ n 2 s 2 + 2 ฮถ ฯ‰ n s + ฯ‰ n 2 ฯ‰ n 1 โˆ’ ฮถ 2 e โˆ’ ฮถ ฯ‰ n t โ‹… s i n ( ฯ‰ n 1 โˆ’ ฮถ 2 t ) โ‹… 1 ( t ) ฯ‰ n 2 s ( s 2 + 2 ฮถ ฯ‰ n s + ฯ‰ n 2 ) ( 1 โˆ’ 1 1 โˆ’ ฮถ 2 e โˆ’ ฮถ ฯ‰ n t โ‹… s i n ( ฯ‰ n 1 โˆ’ ฮถ 2 t + c o s โˆ’ 1 ( ฮถ ) ) ) โ‹… 1 ( t ) To calculate the residues for the remaining multiplicities of the root r1, use this recursive formula: 1 2 ! d 2 d s 2 ( N ( s ) ( s โˆ’ r 1 ) m D ( s ) ) | s = r 1 = K m โˆ’ 2 1 ( m โˆ’ 1 ) ! d m โˆ’ 1 d s m โˆ’ 1 ( N ( s ) ( s โˆ’ r 1 ) m D ( s ) ) | s = r 1 = K m 1.4.5 Examples 1.4.5.1 Example Consider a system described by the following transfer function and find the pole-zero model of the transfer function and its DC gain. G ( s ) = 2 s + 3 S 2 + 3 s + 2 HINT: Use MATLAB software to check your results in this, and the remaining examples in this section. 1.4.5.2 Example Consider a system described by the following transfer function: 5 s 3 + 30 s 2 + 55 s + 30 s 5 + 14 s 4 + 62 s 3 + 110 s 2 + 153 s + 140 Find the pole-zero model of the transfer function and its DC gain. HINT: Use Matlab. 1.4.5.3 Example A certain LTI system is described as having one zero and four poles, as follows: z 1 = โˆ’ 2.2 , p 1 = โˆ’ 1 + j 1 , p 2 = โˆ’ 1 โˆ’ j 1 , p 3 = โˆ’ 10 , p 4 = โˆ’ 2 , It is also recorded that the system has a DC gain of 5. Write the complete transfer function of the system in a ZPK form. 1.4.5.4 Example Consider a system described by the following transfer function: G ( s ) = 10 s 2 + 30 s + 20 s 4 + 14 s 3 + 68 s 2 + 130 s + 75 Create an LTI object representing this system using both the transfer function model and the zero-pole-gain model. Extract zero-pole-gain data and numerator-denominator from the LTI object. Obtain the system dc gain and the pole-zero map of the transfer function. Obtain a minimum realization of this system. Use MATLAB to solve this problem. 1.4.5.5 Example Consider a system described by the following transfer function: G ( s ) = Y ( s ) U ( s ) = s 2 + 3 s + 3 s 3 + 6 s 2 + 11 s + 6 Find an analytical expression for an impulse response of the system. 1.4.5.6 Example A certain control system is described by the following transfer function: G ( s ) = Y ( s ) U ( s ) = 2 s + 8 s 3 + 5 s 2 + 8 s + 4 Find an analytical expression for a step response of the system. 1.4.5.7 Example A certain control system is described by the following closed loop transfer function: G c l ( s ) = 45 ( s + 6 ) ( s 2 + 65 s + 354 ) s Find an analytical expression for a step response of the system โ€“ note the integrator term in the denominator! 1.4.5.8 Example A certain control system is described by the following closed loop transfer function: G c l ( s ) = 45 ( s + 6 ) s 3 + 20 s 2 + 129 s + 270 One of the closed loop poles is at -5. Find an analytical expression for a step response of the system.
Markdown
[Skip to content](https://pressbooks.library.torontomu.ca/controlsystems/chapter/1-4-laplace-transforms/#content) [![Logo for Toronto Metropolitan University Pressbooks](https://pressbooks.library.torontomu.ca/controlsystems/wp-content/uploads/2022/04/TMUPressbooks.png)](https://pressbooks.library.torontomu.ca/) [Toggle Menu](https://pressbooks.library.torontomu.ca/controlsystems/chapter/1-4-laplace-transforms/#navigation) Primary Navigation - [Home](https://pressbooks.library.torontomu.ca/controlsystems) - [Read](https://pressbooks.library.torontomu.ca/controlsystems/front-matter/introduction/) - [Sign in](https://pressbooks.library.torontomu.ca/controlsystems/wp-login.php?redirect_to=https%3A%2F%2Fpressbooks.library.torontomu.ca%2Fcontrolsystems%2Fchapter%2F1-4-laplace-transforms%2F) Book Contents Navigation Contents 1. [Introduction](https://pressbooks.library.torontomu.ca/controlsystems/front-matter/introduction/) 2. [Acknowledgements](https://pressbooks.library.torontomu.ca/controlsystems/front-matter/acknowledgements/) 3. [Foreword](https://pressbooks.library.torontomu.ca/controlsystems/front-matter/foreword/) 4. Chapter 1 1. [1\.1 Definitions: System and Control](https://pressbooks.library.torontomu.ca/controlsystems/chapter/definitions-system-control/) 2. [1\.2 Types of Control Actions](https://pressbooks.library.torontomu.ca/controlsystems/chapter/1-2-types-of-control-actions/) 3. [1\.3 Control Objectives](https://pressbooks.library.torontomu.ca/controlsystems/chapter/1-3-control-objectives/) 4. [1\.4 Laplace Transforms](https://pressbooks.library.torontomu.ca/controlsystems/chapter/1-4-laplace-transforms/) 5. [1\.5 Transfer Function Representations of Simple Physical Systems](https://pressbooks.library.torontomu.ca/controlsystems/chapter/1-5-1-5transfer-function-representations-of-simple-physical-systems/) 6. [1\.6 Basic Block Diagrams](https://pressbooks.library.torontomu.ca/controlsystems/chapter/1-6-basic-block-diagrams/) 5. Chapter 2 1. [2\.1 General Definition of Stability](https://pressbooks.library.torontomu.ca/controlsystems/chapter/2-1-general-definition-of-stability/) 2. [2\.2 Locations in s-Plane vs. Time Response](https://pressbooks.library.torontomu.ca/controlsystems/chapter/2-2-locations-in-s-plane-vs-time-response/) 3. [2\.3 Stability in s-Domain: The Routh-Hurwitz Criterion of Stability](https://pressbooks.library.torontomu.ca/controlsystems/chapter/2-3stability-in-s-domain-the-routh-hurwitz-criterion-of-stability/) 4. [2\.4 Determining Stable Range for Proportional Controller Operations](https://pressbooks.library.torontomu.ca/controlsystems/chapter/2-4determining-stable-range-for-proportional-controller-operations/) 5. [2\.5 Relative Stability - Gain Margin](https://pressbooks.library.torontomu.ca/controlsystems/chapter/2-5-relative-stability-gain-margin/) 6. [2\.6 Examples](https://pressbooks.library.torontomu.ca/controlsystems/chapter/2-6-examples/) 6. Chapter 3 1. [3\.1 Basic Block Diagrams Continued](https://pressbooks.library.torontomu.ca/controlsystems/chapter/3-1basic-block-diagrams-continued/) 2. [3\.2 Signal Flow Graphs](https://pressbooks.library.torontomu.ca/controlsystems/chapter/signal-flow-graphs/) 3. [3\.3 Examples](https://pressbooks.library.torontomu.ca/controlsystems/chapter/3-3-examples/) 7. Chapter 4 1. [4\.1 Introduction](https://pressbooks.library.torontomu.ca/controlsystems/chapter/4-1-introduction/) 2. [4\.2 Standard Time Inputs](https://pressbooks.library.torontomu.ca/controlsystems/chapter/4-2-standard-time-inputs/) 3. [4\.3 Step response specifications - Definitions](https://pressbooks.library.torontomu.ca/controlsystems/chapter/4-3-step-response-specifications-definitions/) 4. [4\.4 Examples](https://pressbooks.library.torontomu.ca/controlsystems/chapter/4-4-examples/) 8. Chapter 5 1. [5\.1 Equivalent Unit Feedback Loop](https://pressbooks.library.torontomu.ca/controlsystems/chapter/5-1-equivalent-unit-feedback-loop/) 2. [5\.2 Steady State Error Analysis in an Equivalent Unit Feedback Loop](https://pressbooks.library.torontomu.ca/controlsystems/chapter/5-2steady-state-error-analysis-in-an-equivalent-unit-feedback-loop/) 3. [5\.3 Examples](https://pressbooks.library.torontomu.ca/controlsystems/chapter/5-3-examples/) 9. Chapter 6 1. [6\.1 First order systems](https://pressbooks.library.torontomu.ca/controlsystems/chapter/6-1-first-order-systems/) 2. [6\.2 Second Order Overdamped Systems](https://pressbooks.library.torontomu.ca/controlsystems/chapter/6-2second-order-overdamped-systems/) 10. Chapter 7 1. [7\.1 Second Order Underdamped Systems](https://pressbooks.library.torontomu.ca/controlsystems/chapter/7-1second-order-underdamped-systems/) 2. [7\.2 Response Specifications for the Second Order Underdamped System](https://pressbooks.library.torontomu.ca/controlsystems/chapter/7-2response-specifications-for-the-second-order-underdamped-system/) 3. [7\.3 Examples](https://pressbooks.library.torontomu.ca/controlsystems/chapter/7-3-examples/) 11. Chapter 8 1. [8\.1 Systems with Delay](https://pressbooks.library.torontomu.ca/controlsystems/chapter/8-1-systems-with-delay/) 2. [8\.2 Minimum Realizations and Reduced Order Models - Part 1](https://pressbooks.library.torontomu.ca/controlsystems/chapter/8-2minimum-realizations-and-reduced-order-models-part-1/) 3. [8\.3 Dominant System Dynamics and Reduced Order Models - Part 2](https://pressbooks.library.torontomu.ca/controlsystems/chapter/8-3dominant-system-dynamics-and-reduced-order-models-part-2/) 4. [8\.4 The Effect of an Additional Pole on the 2nd Order System Response](https://pressbooks.library.torontomu.ca/controlsystems/chapter/8-4the-effect-of-an-additional-pole-on-the-2nd-order-system-response/) 5. [8\.5 The Effect of an additional Zero on the 2nd Order System Response](https://pressbooks.library.torontomu.ca/controlsystems/chapter/8-5the-effect-of-an-additional-zero-on-the-2nd-order-system-response/) 6. [8\.6 The Effect of a Non-Minimum Phase Zero on the 2nd Order System Response](https://pressbooks.library.torontomu.ca/controlsystems/chapter/8-6the-effect-of-a-non-minimum-phase-zero-on-the-2nd-order-system-response/) 7. [8\.7 Examples](https://pressbooks.library.torontomu.ca/controlsystems/chapter/8-7examples/) 12. Chapter 9 1. [9\.1 Introduction](https://pressbooks.library.torontomu.ca/controlsystems/chapter/9-1introduction/) 2. [9\.2 Proportional Control](https://pressbooks.library.torontomu.ca/controlsystems/chapter/9-2proportional-control/) 3. [9\.3 Proportional + Derivative Control](https://pressbooks.library.torontomu.ca/controlsystems/chapter/9-3proportional-derivative-control/) 4. [9\.4 Proportional + Integral Control](https://pressbooks.library.torontomu.ca/controlsystems/chapter/9-4proportional-integral-control/) 5. [9\.5 PID Controller and Its Tuning](https://pressbooks.library.torontomu.ca/controlsystems/chapter/9-5pid-controller-and-its-tuning/) 6. [9\.6 Effect of PID Controller Modes on System Stability](https://pressbooks.library.torontomu.ca/controlsystems/chapter/9-6effect-of-pid-controller-modes-on-system-stability/) 7. [9\.7 Examples](https://pressbooks.library.torontomu.ca/controlsystems/chapter/9-7examples/) 13. Chapter 10 1. [10\.1 Introduction](https://pressbooks.library.torontomu.ca/controlsystems/chapter/10-1introduction/) 2. [10\.2 Evans' Root Locus Construction Rules - Introduction](https://pressbooks.library.torontomu.ca/controlsystems/chapter/10-2evans-root-locus-construction-rules-introduction/) 3. [10\.3 Evans Root Locus Construction Rule \# 1: Beginning, End and Symmetry](https://pressbooks.library.torontomu.ca/controlsystems/chapter/10-3evans-root-locus-construction-rule-1-beginning-end-and-symmetry/) 4. [10\.4 Evans Root Locus Construction Rule \# 2: Segments of Root Locus on Real Axis](https://pressbooks.library.torontomu.ca/controlsystems/chapter/10-4evans-root-locus-construction-rule-2-segments-of-root-locus-on-real-axis/) 5. [10\.5 Evans Root Locus Construction Rule \# 3: Asymptotic Angles and Centroid](https://pressbooks.library.torontomu.ca/controlsystems/chapter/10-5evans-root-locus-construction-rule-3-asymptotic-angles-and-centroid/) 6. [10\.6 Evans Root Locus Construction Rule \# 4: Break-Away and Break-In Points](https://pressbooks.library.torontomu.ca/controlsystems/chapter/10-6-evans-root-locus-construction-rule-4break-away-and-break-in-points/) 7. [10\.7 Evans Root Locus Construction Rule \# 5: Crossover with Imaginary Axis](https://pressbooks.library.torontomu.ca/controlsystems/chapter/10-7-evans-root-locus-construction-rule-5-crosspver-with-imaginary-axis/) 8. [10\.8 Evans Root Locus Construction Rule \#6: Angles of Departures/Arrivals at Complex Poles/Zeros](https://pressbooks.library.torontomu.ca/controlsystems/chapter/10-8-evans-root-locus-construction-rule-6-angles-of-departures-arrivals-at-complex-poles-zeros/) 9. [10\.9 Examples](https://pressbooks.library.torontomu.ca/controlsystems/chapter/10-9-examples/) 14. Chapter 11 1. [11\.1 Gain Margin from Bode Plot](https://pressbooks.library.torontomu.ca/controlsystems/chapter/11-1-gain-margin-from-bode-plot/) 2. [11\.2 Definition of Phase Margin](https://pressbooks.library.torontomu.ca/controlsystems/chapter/11-2definition-of-phase-margin/) 3. [11\.3 Examples](https://pressbooks.library.torontomu.ca/controlsystems/chapter/11-3-examples/) 15. Chapter 12 1. [12\.1 Model from Closed Loop Frequency Response](https://pressbooks.library.torontomu.ca/controlsystems/chapter/12-1-model-from-closed-loop-frequency-response/) 2. [12\.2 Model from Open Loop Frequency Response](https://pressbooks.library.torontomu.ca/controlsystems/chapter/12-2-model-from-open-loop-frequency-response/) 3. [12\.3 Summary](https://pressbooks.library.torontomu.ca/controlsystems/chapter/12-3summary/) 4. [12\.4 Examples](https://pressbooks.library.torontomu.ca/controlsystems/chapter/12-4examples/) 16. Chapter 13 1. [13\.1 Basic Rules - Summary](https://pressbooks.library.torontomu.ca/controlsystems/chapter/13-1-basic-rules-summary/) 2. [13\.2 Lead Controller](https://pressbooks.library.torontomu.ca/controlsystems/chapter/13-2lead-controller/) 3. [13\.3 Lead Controller Design โ€“ Solved Examples](https://pressbooks.library.torontomu.ca/controlsystems/chapter/13-3-lead-controller-design-solved-examples/) 4. [13\.4 Lag Controller](https://pressbooks.library.torontomu.ca/controlsystems/chapter/13-4-lag-controller/) 5. [13\.5 Lag Controller Design โ€“ Solved Example 1](https://pressbooks.library.torontomu.ca/controlsystems/chapter/13-5lag-controller-design-solved-example-1/) 6. [13\.6 Lead-Lag Controller](https://pressbooks.library.torontomu.ca/controlsystems/chapter/13-6lead-lag-controller/) 7. [13\.7 Examples](https://pressbooks.library.torontomu.ca/controlsystems/chapter/13-7examples/) 17. Chapter 14 1. [14\.1 Why We Need Another Criterion of Stability](https://pressbooks.library.torontomu.ca/controlsystems/chapter/14-1-why-we-need-another-criterion-of-stability/) 2. [14\.2 Polar Plots Revisited](https://pressbooks.library.torontomu.ca/controlsystems/chapter/14-2-polar-plots-revisited/) 3. [14\.3 Solved Examples for Polar Plots](https://pressbooks.library.torontomu.ca/controlsystems/chapter/14-3-solved-examples-for-polar-plots/) 4. [14\.4 Gain and Phase Margins vs. Polar Plots](https://pressbooks.library.torontomu.ca/controlsystems/chapter/14-4-1-gain-margin-from-polar-plot/) 5. [14\.5 Concept of Mapping](https://pressbooks.library.torontomu.ca/controlsystems/chapter/14-5-concept-of-mapping/) 6. [14\.6 Cauchy's Mapping Theorem](https://pressbooks.library.torontomu.ca/controlsystems/chapter/14-6-cauchys-mapping-theorem/) 7. [14\.7 Solved Examples of Nyquist Stability Criterion](https://pressbooks.library.torontomu.ca/controlsystems/chapter/14-7-solved-examples-of-nyquist-stability-criterion/) # [Introduction to Control Systems](https://pressbooks.library.torontomu.ca/controlsystems/) Chapter 1 # 1\.4 Laplace Transforms Self-Study: Review your ELE532 Notes and other resources. You can also refer to the review material on the course website. ### **1\.4.1 Definitions** | | |---| | F(s)\=L\[f(t)\]\=โˆซ\+โˆž0f(t)eโˆ’stdt F ( s ) \= L \[ f ( t ) \] \= โˆซ 0 \+ โˆž f ( t ) e โˆ’ s t d t f(t)\=Lโˆ’1\[F(s)\]\=12ฯ€jโˆซฯƒ\+jโˆžฯƒโˆ’jโˆžF(s)estds f ( t ) \= L โˆ’ 1 \[ F ( s ) \] \= 1 2 ฯ€ j โˆซ ฯƒ โˆ’ j โˆž ฯƒ \+ j โˆž F ( s ) e s t d s | #### **1\.4.1.1 Final Value Theorem** | | |---| | fss\=limtโ†’โˆžf(t)\=limsโ†’0sF(s) f s s \= lim t โ†’ โˆž f ( t ) \= lim s โ†’ 0 s F ( s ) | #### **1\.4.1.2 Initial Value Theorem** | | |---| | f0\=limtโ†’0\+f(t)\=limsโ†’โˆžsF(s) f 0 \= lim t โ†’ 0 \+ f ( t ) \= lim s โ†’ โˆž s F ( s ) | #### **1\.4.1.3 Properties of Laplace transforms** | | |---| | F(s)eโˆ’Ts F ( s ) e โˆ’ T s | ### **1\.4.2 Solving for System Response** Parametric models cannot be developed without math. Laws of physics describe dynamic Linear Time-Invariant (LTI) systems using ordinary differential equations. To simplify their analysis, Laplace transform is used. Consider a certain LTI (Linear Time-Invariant), SISO (Single Input Single Output) system: ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/12.png) Let the input โ€“ output relationship for the system be described by the following *n*th order differential equation: | | |---| | dnydtn\+anโˆ’1dnโˆ’1ydtnโˆ’1\+...\+a1dydt\+a0y\=bmdmudtm\+bmโˆ’1dmโˆ’1udtmโˆ’1\+...\+b1dudt\+b0u d n y d t n \+ a n โˆ’ 1 d n โˆ’ 1 y d t n โˆ’ 1 \+ . . . \+ a 1 d y d t \+ a 0 y \= b m d m u d t m \+ b m โˆ’ 1 d m โˆ’ 1 u d t m โˆ’ 1 \+ . . . \+ b 1 d u d t \+ b 0 u | The equation parameters relate to physical aspects of the system. The time domain description of systems is not convenient for quick paper-and-pencil speculations. To simplify math, Classical Control uses a Laplace Transform system description, which converts the differential equations into their algebraic equivalents in the s-domain. The solution for y(t) can then be found using inverse Laplace transformation to Y(s). ### **1\.4.3 Two Transfer Functions Models: TF and ZPK** In the transform domain, the input-output relationship of the system is defined by a transfer function G(s), defined as a ratio of the Laplace transform of the system output signal y(t), to the Laplace transform of the system input signal u(t), with any initial conditions in the system set to zero. The system transfer function G(s) can be thought of as a dynamic gain of the system: ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/111-300x104.png) Block diagrams are used to graphically represent systems and their components, as shown above. In order to find G(s), a Laplace transform of the system differential equation in Equation 1โ€‘4 is taken: | | |---| | snY(s)\+anโˆ’1snโˆ’1Y(s)\+...\+a1sY(s)\+a0Y(s)\=bmsmU(s)\+bmโˆ’1smโˆ’1U(s)\+...\+b1sU(s)\+b0U(s) s n Y ( s ) \+ a n โˆ’ 1 s n โˆ’ 1 Y ( s ) \+ . . . \+ a 1 s Y ( s ) \+ a 0 Y ( s ) \= b m s m U ( s ) \+ b m โˆ’ 1 s m โˆ’ 1 U ( s ) \+ . . . \+ b 1 s U ( s ) \+ b 0 U ( s ) G(s)\=Y(s)U(s)\=bmsm\+bmโˆ’1smโˆ’1\+...\+b1s\+b0sn\+anโˆ’1snโˆ’1\+...\+a1s\+a0 G ( s ) \= Y ( s ) U ( s ) \= b m s m \+ b m โˆ’ 1 s m โˆ’ 1 \+ . . . \+ b 1 s \+ b 0 s n \+ a n โˆ’ 1 s n โˆ’ 1 \+ . . . \+ a 1 s \+ a 0 | Transfer functions are ratios of polynomials written in terms of the s-operator. The resulting function in Equation 1โ€‘5 is a ratio of two polynomials, N(s) and D(s): | | |---| | G(s)\=N(s)D(s)\=bmsm\+bmโˆ’1smโˆ’1\+...\+b1s\+b0sn\+anโˆ’1snโˆ’1\+...\+a1s\+a0 G ( s ) \= N ( s ) D ( s ) \= b m s m \+ b m โˆ’ 1 s m โˆ’ 1 \+ . . . \+ b 1 s \+ b 0 s n \+ a n โˆ’ 1 s n โˆ’ 1 \+ . . . \+ a 1 s \+ a 0 | Roots of the numerator polynomial of G(s) in Equation 1โ€‘6 are called system zeros, zi z i, and roots of the denominator polynomial are called system poles, pi p i. ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/fig_MATLAB_LOGO.png)MATLAB Enabled Note that in MATLAB a transfer function object in a polynomial form can be created by using the โ€œtfโ€ command. For example, consider the following transfer function in the polynomial ratio form: ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/12-1.png) G(s)\=Y(s)U(s)\=2s\+20s2\+4s\+3 G ( s ) \= Y ( s ) U ( s ) \= 2 s \+ 20 s 2 \+ 4 s \+ 3 The same transfer function G(s) G ( s ) can be represented in the so-called ZPK form (factorized form): | | |---| | G(s)\=Kโˆmi(sโˆ’zi)โˆnj(sโˆ’pj) G ( s ) \= K โˆ i m ( s โˆ’ z i ) โˆ j n ( s โˆ’ p j ) | *K* is a multiplier. It is important to see the difference between *K* and *K*dc, which denotes the DC gain of the system (i.e. s=0): | | |---| | Kdc\=G(0)\=N(0)D(0)\=b0a0\=Kโˆmi(โˆ’zi)โˆnj(โˆ’pj) K d c \= G ( 0 ) \= N ( 0 ) D ( 0 ) \= b 0 a 0 \= K โˆ i m ( โˆ’ z i ) โˆ j n ( โˆ’ p j ) | Our transfer function can be factorized and the multiplier gain K is equal to 2: | | |---| | G(s)\=Y(s)U(s)\=2(s\+10)(s\+1)(s\+3) G ( s ) \= Y ( s ) U ( s ) \= 2 ( s \+ 10 ) ( s \+ 1 ) ( s \+ 3 ) | The DC gain of our transfer function is: | | |---| | KDC\=G(0)\=2(10)(1)(3)\=6\.667 K D C \= G ( 0 ) \= 2 ( 10 ) ( 1 ) ( 3 ) \= 6\.667 | ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/fig_MATLAB_LOGO.png) In MATLAB the transfer function can be shown in a factorized form by using the โ€œzpkโ€ command, and the DC gain can be found using โ€œdcgainโ€ command: ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/13.png) ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/fig_MATLAB_LOGO.png) In MATLAB the โ€œzpkโ€ command can be used to create a transfer function object in a ZPK form, and โ€œtfโ€ command to convert it to a polynomial form. ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/14.png) Locations of the system poles and zeros can be presented graphically as the so-called Pole-Zero Map. ### **1\.4.4 Partial Fractions Technique** If a certain control system is described by a transfer function G(s), the system response can be found as Y(s)\=U(s)โ‹…G(s) Y ( s ) \= U ( s ) โ‹… G ( s ). Since Laplace Transform Tables do not provide exhaustive solutions, a technique of a Partial Fractions Expansion is used to find inverse Laplace Transforms for various time functions โ€“ see a table of basic Laplace โ€“ Time Domain Function pair shown in Table 1โ€‘2. #### **1\.4.4.1 Residues โ€“ Distinct Roots Case** | | |---| | Y(s)\=N(s)โˆni(sโˆ’pi)\=N(s)D(s) Y ( s ) \= N ( s ) โˆ i n ( s โˆ’ p i ) \= N ( s ) D ( s ) Y(s)\=K1sโˆ’p1\+K2sโˆ’p2\+...\+Knsโˆ’pn Y ( s ) \= K 1 s โˆ’ p 1 \+ K 2 s โˆ’ p 2 \+ . . . \+ K n s โˆ’ p n | #### **1\.4.4.2 Residues โ€“ Multiple Roots Case** ``` ``` This is the case for Laplace Transforms with multiple powers of some roots. Assume multiplicity of *m* for root r1 r 1: | | |---| | Y(s)\=N(s)โˆni(sโˆ’pi)\=N(s)D(s) Y ( s ) \= N ( s ) โˆ i n ( s โˆ’ p i ) \= N ( s ) D ( s ) Y(s)\=K1sโˆ’r1\+K2(sโˆ’r1)2\+...\+Km(sโˆ’r1)m\+...\+Knsโˆ’pn Y ( s ) \= K 1 s โˆ’ r 1 \+ K 2 ( s โˆ’ r 1 ) 2 \+ . . . \+ K m ( s โˆ’ r 1 ) m \+ . . . \+ K n s โˆ’ p n | Residues for distinct roots calculated as before: | | |---| | Km\+1\=N(s)(sโˆ’p1)D(s)\|s\=pm\+1 K m \+ 1 \= N ( s ) ( s โˆ’ p 1 ) D ( s ) \| s \= p m \+ 1 โ‹ฎ โ‹ฎ Kn\=N(s)(sโˆ’pn)D(s)\|s\=pn K n \= N ( s ) ( s โˆ’ p n ) D ( s ) \| s \= p n | To calculate residues for a multiple root, multiply both sides of Equation 1โ€‘12 by (sโˆ’r1)m ( s โˆ’ r 1 ) mand substitute the value of the root r1 r 1: | | |---| | N(s)D(s)\=K1sโˆ’r1\+K2(sโˆ’r1)2\+...\+Km(sโˆ’r1)m\+...\+Knsโˆ’pn N ( s ) D ( s ) \= K 1 s โˆ’ r 1 \+ K 2 ( s โˆ’ r 1 ) 2 \+ . . . \+ K m ( s โˆ’ r 1 ) m \+ . . . \+ K n s โˆ’ p n N(s)(sโˆ’r1)mD(s)\=K1(sโˆ’r1)msโˆ’r1\+K2(sโˆ’r1)m(sโˆ’r1)2\+...\+Km\+...\+Kn(sโˆ’r1)msโˆ’pn N ( s ) ( s โˆ’ r 1 ) m D ( s ) \= K 1 ( s โˆ’ r 1 ) m s โˆ’ r 1 \+ K 2 ( s โˆ’ r 1 ) m ( s โˆ’ r 1 ) 2 \+ . . . \+ K m \+ . . . \+ K n ( s โˆ’ r 1 ) m s โˆ’ p n N(s)(sโˆ’r1)mD(S)\|s\=r1\=Km N ( s ) ( s โˆ’ r 1 ) m D ( S ) \| s \= r 1 \= K m | To calculate the residues for the second multiplicity of the root r1 r 1: | | |---| | dds(N(s)(sโˆ’r1)mD(s))\=dds(K1(sโˆ’r1)msโˆ’r1\+K2(sโˆ’r1)m(sโˆ’r1)2\+...\+Km\+...\+Kn(sโˆ’r1)msโˆ’pn) d d s ( N ( s ) ( s โˆ’ r 1 ) m D ( s ) ) \= d d s ( K 1 ( s โˆ’ r 1 ) m s โˆ’ r 1 \+ K 2 ( s โˆ’ r 1 ) m ( s โˆ’ r 1 ) 2 \+ . . . \+ K m \+ . . . \+ K n ( s โˆ’ r 1 ) m s โˆ’ p n ) dds(N(s)(sโˆ’r1)mD(s))\|s\=r1\=0\+...\+Kmโˆ’1\+0\+...\+0 d d s ( N ( s ) ( s โˆ’ r 1 ) m D ( s ) ) \| s \= r 1 \= 0 \+ . . . \+ K m โˆ’ 1 \+ 0 \+ . . . \+ 0 | | | | |---|---| | Laplace Transform | Time Domain Function | | 1 1 | | To calculate the residues for the remaining multiplicities of the root r1, use this recursive formula: | | |---| | 12\!d2ds2(N(s)(sโˆ’r1)mD(s))\|s\=r1\=Kmโˆ’2 1 2 \! d 2 d s 2 ( N ( s ) ( s โˆ’ r 1 ) m D ( s ) ) \| s \= r 1 \= K m โˆ’ 2 1(mโˆ’1)\!dmโˆ’1dsmโˆ’1(N(s)(sโˆ’r1)mD(s))\|s\=r1\=Km 1 ( m โˆ’ 1 ) \! d m โˆ’ 1 d s m โˆ’ 1 ( N ( s ) ( s โˆ’ r 1 ) m D ( s ) ) \| s \= r 1 \= K m | ### **1\.4.5 Examples** #### **1\.4.5.1 Example** Consider a system described by the following transfer function and find the pole-zero model of the transfer function and its DC gain. G(s)\=2s\+3S2\+3s\+2 G ( s ) \= 2 s \+ 3 S 2 \+ 3 s \+ 2 HINT: Use MATLAB software to check your results in this, and the remaining examples in this section. #### **1\.4.5.2 Example** Consider a system described by the following transfer function: ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/fig_MATLAB_LOGO.png) 5s3\+30s2\+55s\+30s5\+14s4\+62s3\+110s2\+153s\+140 5 s 3 \+ 30 s 2 \+ 55 s \+ 30 s 5 \+ 14 s 4 \+ 62 s 3 \+ 110 s 2 \+ 153 s \+ 140 Find the pole-zero model of the transfer function and its DC gain. **HINT: Use Matlab.** #### **1\.4.5.3 Example** A certain LTI system is described as having one zero and four poles, as follows: z1\=โˆ’2\.2 z 1 \= โˆ’ 2\.2, p1\=โˆ’1\+j1 p 1 \= โˆ’ 1 \+ j 1, p2\=โˆ’1โˆ’j1 p 2 \= โˆ’ 1 โˆ’ j 1, p3\=โˆ’10 p 3 \= โˆ’ 10, p4\=โˆ’2 p 4 \= โˆ’ 2, It is also recorded that the system has a DC gain of 5. Write the complete transfer function of the system in a ZPK form. #### **1\.4.5.4 Example** Consider a system described by the following transfer function: G(s)\=10s2\+30s\+20s4\+14s3\+68s2\+130s\+75 G ( s ) \= 10 s 2 \+ 30 s \+ 20 s 4 \+ 14 s 3 \+ 68 s 2 \+ 130 s \+ 75 Create an LTI object representing this system using both the transfer function model and the zero-pole-gain model. Extract zero-pole-gain data and numerator-denominator from the LTI object. Obtain the system dc gain and the pole-zero map of the transfer function. Obtain a minimum realization of this system. Use MATLAB to solve this problem. #### **1\.4.5.5 Example** Consider a system described by the following transfer function: G(s)\=Y(s)U(s)\=s2\+3s\+3s3\+6s2\+11s\+6 G ( s ) \= Y ( s ) U ( s ) \= s 2 \+ 3 s \+ 3 s 3 \+ 6 s 2 \+ 11 s \+ 6 Find an analytical expression for an impulse response of the system. #### **1\.4.5.6 Example** A certain control system is described by the following transfer function: G(s)\=Y(s)U(s)\=2s\+8s3\+5s2\+8s\+4 G ( s ) \= Y ( s ) U ( s ) \= 2 s \+ 8 s 3 \+ 5 s 2 \+ 8 s \+ 4 Find an analytical expression for a step response of the system. #### **1\.4.5.7 Example** A certain control system is described by the following closed loop transfer function: Gcl(s)\=45(s\+6)(s2\+65s\+354)s G c l ( s ) \= 45 ( s \+ 6 ) ( s 2 \+ 65 s \+ 354 ) s Find an analytical expression for a step response of the system โ€“ note the integrator term in the denominator\! #### **1\.4.5.8 Example** A certain control system is described by the following closed loop transfer function: Gcl(s)\=45(s\+6)s3\+20s2\+129s\+270 G c l ( s ) \= 45 ( s \+ 6 ) s 3 \+ 20 s 2 \+ 129 s \+ 270 One of the closed loop poles is at -5. Find an analytical expression for a step response of the system. ``` ``` Previous/next navigation [Previous: 1.3 Control Objectives](https://pressbooks.library.torontomu.ca/controlsystems/chapter/1-3-control-objectives/ "Previous: 1.3 Control Objectives") [Next: 1.5 Transfer Function Representations of Simple Physical Systems](https://pressbooks.library.torontomu.ca/controlsystems/chapter/1-5-1-5transfer-function-representations-of-simple-physical-systems/ "Next: 1.5 Transfer Function Representations of Simple Physical Systems") Back to top ## License ![Icon for the Creative Commons Attribution-NonCommercial 4.0 International License](https://pressbooks.library.torontomu.ca/controlsystems/wp-content/themes/pressbooks-book/packages/buckram/assets/images/cc-by-nc.svg) [Introduction to Control Systems](https://pressbooks.library.torontomu.ca/controlsystems) Copyright ยฉ by Malgorzata Zywno is licensed under a [Creative Commons Attribution-NonCommercial 4.0 International License](https://creativecommons.org/licenses/by-nc/4.0/), except where otherwise noted. [Powered by Pressbooks](https://pressbooks.com/) - [Guides and Tutorials](https://pressbooks.com/support/) - \|[Pressbooks Directory](https://pressbooks.directory/) [![YouTube](https://pressbooks.library.torontomu.ca/controlsystems/wp-content/themes/pressbooks-book/assets/images/yt_icon_mono_dark.png) Pressbooks on YouTube](https://www.youtube.com/user/pressbooks) [Pressbooks on Twitter](https://twitter.com/intent/follow?screen_name=pressbooks)
Readable Markdown
Self-Study: Review your ELE532 Notes and other resources. You can also refer to the review material on the course website. ### **1\.4.1 Definitions** F ( s ) \= L \[ f ( t ) \] \= โˆซ 0 \+ โˆž f ( t ) e โˆ’ s t d t f ( t ) \= L โˆ’ 1 \[ F ( s ) \] \= 1 2 ฯ€ j โˆซ ฯƒ โˆ’ j โˆž ฯƒ \+ j โˆž F ( s ) e s t d s #### **1\.4.1.1 Final Value Theorem** f s s \= lim t โ†’ โˆž f ( t ) \= lim s โ†’ 0 s F ( s ) #### **1\.4.1.2 Initial Value Theorem** f 0 \= lim t โ†’ 0 \+ f ( t ) \= lim s โ†’ โˆž s F ( s ) #### **1\.4.1.3 Properties of Laplace transforms** | | |---| | F ( s ) e โˆ’ T s | ### **1\.4.2 Solving for System Response** Parametric models cannot be developed without math. Laws of physics describe dynamic Linear Time-Invariant (LTI) systems using ordinary differential equations. To simplify their analysis, Laplace transform is used. Consider a certain LTI (Linear Time-Invariant), SISO (Single Input Single Output) system: ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/12.png) Let the input โ€“ output relationship for the system be described by the following *n*th order differential equation: d n y d t n \+ a n โˆ’ 1 d n โˆ’ 1 y d t n โˆ’ 1 \+ . . . \+ a 1 d y d t \+ a 0 y \= b m d m u d t m \+ b m โˆ’ 1 d m โˆ’ 1 u d t m โˆ’ 1 \+ . . . \+ b 1 d u d t \+ b 0 u The equation parameters relate to physical aspects of the system. The time domain description of systems is not convenient for quick paper-and-pencil speculations. To simplify math, Classical Control uses a Laplace Transform system description, which converts the differential equations into their algebraic equivalents in the s-domain. The solution for y(t) can then be found using inverse Laplace transformation to Y(s). ### **1\.4.3 Two Transfer Functions Models: TF and ZPK** In the transform domain, the input-output relationship of the system is defined by a transfer function G(s), defined as a ratio of the Laplace transform of the system output signal y(t), to the Laplace transform of the system input signal u(t), with any initial conditions in the system set to zero. The system transfer function G(s) can be thought of as a dynamic gain of the system: ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/111-300x104.png) Block diagrams are used to graphically represent systems and their components, as shown above. In order to find G(s), a Laplace transform of the system differential equation in Equation 1โ€‘4 is taken: | | |---| | s n Y ( s ) \+ a n โˆ’ 1 s n โˆ’ 1 Y ( s ) \+ . . . \+ a 1 s Y ( s ) \+ a 0 Y ( s ) \= b m s m U ( s ) \+ b m โˆ’ 1 s m โˆ’ 1 U ( s ) \+ . . . \+ b 1 s U ( s ) \+ b 0 U ( s ) G ( s ) \= Y ( s ) U ( s ) \= b m s m \+ b m โˆ’ 1 s m โˆ’ 1 \+ . . . \+ b 1 s \+ b 0 s n \+ a n โˆ’ 1 s n โˆ’ 1 \+ . . . \+ a 1 s \+ a 0 | Transfer functions are ratios of polynomials written in terms of the s-operator. The resulting function in Equation 1โ€‘5 is a ratio of two polynomials, N(s) and D(s): | | |---| | G ( s ) \= N ( s ) D ( s ) \= b m s m \+ b m โˆ’ 1 s m โˆ’ 1 \+ . . . \+ b 1 s \+ b 0 s n \+ a n โˆ’ 1 s n โˆ’ 1 \+ . . . \+ a 1 s \+ a 0 | Roots of the numerator polynomial of G(s) in Equation 1โ€‘6 are called system zeros, z i, and roots of the denominator polynomial are called system poles, p i. Note that in MATLAB a transfer function object in a polynomial form can be created by using the โ€œtfโ€ command. For example, consider the following transfer function in the polynomial ratio form: ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/12-1.png) G ( s ) \= Y ( s ) U ( s ) \= 2 s \+ 20 s 2 \+ 4 s \+ 3 The same transfer function G ( s ) can be represented in the so-called ZPK form (factorized form): | | |---| | G ( s ) \= K โˆ i m ( s โˆ’ z i ) โˆ j n ( s โˆ’ p j ) | *K* is a multiplier. It is important to see the difference between *K* and *K*dc, which denotes the DC gain of the system (i.e. s=0): | | |---| | K d c \= G ( 0 ) \= N ( 0 ) D ( 0 ) \= b 0 a 0 \= K โˆ i m ( โˆ’ z i ) โˆ j n ( โˆ’ p j ) | Our transfer function can be factorized and the multiplier gain K is equal to 2: G ( s ) \= Y ( s ) U ( s ) \= 2 ( s \+ 10 ) ( s \+ 1 ) ( s \+ 3 ) The DC gain of our transfer function is: K D C \= G ( 0 ) \= 2 ( 10 ) ( 1 ) ( 3 ) \= 6\.667 In MATLAB the transfer function can be shown in a factorized form by using the โ€œzpkโ€ command, and the DC gain can be found using โ€œdcgainโ€ command: ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/13.png) In MATLAB the โ€œzpkโ€ command can be used to create a transfer function object in a ZPK form, and โ€œtfโ€ command to convert it to a polynomial form. ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/14.png) Locations of the system poles and zeros can be presented graphically as the so-called Pole-Zero Map. ### **1\.4.4 Partial Fractions Technique** If a certain control system is described by a transfer function G(s), the system response can be found as Y ( s ) \= U ( s ) โ‹… G ( s ). Since Laplace Transform Tables do not provide exhaustive solutions, a technique of a Partial Fractions Expansion is used to find inverse Laplace Transforms for various time functions โ€“ see a table of basic Laplace โ€“ Time Domain Function pair shown in Table 1โ€‘2. #### **1\.4.4.1 Residues โ€“ Distinct Roots Case** | | |---| | Y ( s ) \= N ( s ) โˆ i n ( s โˆ’ p i ) \= N ( s ) D ( s ) Y ( s ) \= K 1 s โˆ’ p 1 \+ K 2 s โˆ’ p 2 \+ . . . \+ K n s โˆ’ p n | #### **1\.4.4.2 Residues โ€“ Multiple Roots Case** ``` ``` This is the case for Laplace Transforms with multiple powers of some roots. Assume multiplicity of *m* for root r 1: Y ( s ) \= N ( s ) โˆ i n ( s โˆ’ p i ) \= N ( s ) D ( s ) Y ( s ) \= K 1 s โˆ’ r 1 \+ K 2 ( s โˆ’ r 1 ) 2 \+ . . . \+ K m ( s โˆ’ r 1 ) m \+ . . . \+ K n s โˆ’ p n Residues for distinct roots calculated as before: K m \+ 1 \= N ( s ) ( s โˆ’ p 1 ) D ( s ) \| s \= p m \+ 1 โ‹ฎ K n \= N ( s ) ( s โˆ’ p n ) D ( s ) \| s \= p n To calculate residues for a multiple root, multiply both sides of Equation 1โ€‘12 by ( s โˆ’ r 1 ) mand substitute the value of the root r 1: N ( s ) D ( s ) \= K 1 s โˆ’ r 1 \+ K 2 ( s โˆ’ r 1 ) 2 \+ . . . \+ K m ( s โˆ’ r 1 ) m \+ . . . \+ K n s โˆ’ p n N ( s ) ( s โˆ’ r 1 ) m D ( s ) \= K 1 ( s โˆ’ r 1 ) m s โˆ’ r 1 \+ K 2 ( s โˆ’ r 1 ) m ( s โˆ’ r 1 ) 2 \+ . . . \+ K m \+ . . . \+ K n ( s โˆ’ r 1 ) m s โˆ’ p n N ( s ) ( s โˆ’ r 1 ) m D ( S ) \| s \= r 1 \= K m To calculate the residues for the second multiplicity of the root r 1: d d s ( N ( s ) ( s โˆ’ r 1 ) m D ( s ) ) \= d d s ( K 1 ( s โˆ’ r 1 ) m s โˆ’ r 1 \+ K 2 ( s โˆ’ r 1 ) m ( s โˆ’ r 1 ) 2 \+ . . . \+ K m \+ . . . \+ K n ( s โˆ’ r 1 ) m s โˆ’ p n ) d d s ( N ( s ) ( s โˆ’ r 1 ) m D ( s ) ) \| s \= r 1 \= 0 \+ . . . \+ K m โˆ’ 1 \+ 0 \+ . . . \+ 0 | | | |---|---| | Laplace Transform | Time Domain Function | | 1 | | To calculate the residues for the remaining multiplicities of the root r1, use this recursive formula: 1 2 \! d 2 d s 2 ( N ( s ) ( s โˆ’ r 1 ) m D ( s ) ) \| s \= r 1 \= K m โˆ’ 2 1 ( m โˆ’ 1 ) \! d m โˆ’ 1 d s m โˆ’ 1 ( N ( s ) ( s โˆ’ r 1 ) m D ( s ) ) \| s \= r 1 \= K m ### **1\.4.5 Examples** #### **1\.4.5.1 Example** Consider a system described by the following transfer function and find the pole-zero model of the transfer function and its DC gain. G ( s ) \= 2 s \+ 3 S 2 \+ 3 s \+ 2 HINT: Use MATLAB software to check your results in this, and the remaining examples in this section. #### **1\.4.5.2 Example** Consider a system described by the following transfer function: ![](http://pressbooks.library.ryerson.ca/controlsystems/wp-content/uploads/sites/75/2019/05/fig_MATLAB_LOGO.png) 5 s 3 \+ 30 s 2 \+ 55 s \+ 30 s 5 \+ 14 s 4 \+ 62 s 3 \+ 110 s 2 \+ 153 s \+ 140 Find the pole-zero model of the transfer function and its DC gain. **HINT: Use Matlab.** #### **1\.4.5.3 Example** A certain LTI system is described as having one zero and four poles, as follows: z 1 \= โˆ’ 2\.2, p 1 \= โˆ’ 1 \+ j 1, p 2 \= โˆ’ 1 โˆ’ j 1, p 3 \= โˆ’ 10, p 4 \= โˆ’ 2, It is also recorded that the system has a DC gain of 5. Write the complete transfer function of the system in a ZPK form. #### **1\.4.5.4 Example** Consider a system described by the following transfer function: G ( s ) \= 10 s 2 \+ 30 s \+ 20 s 4 \+ 14 s 3 \+ 68 s 2 \+ 130 s \+ 75 Create an LTI object representing this system using both the transfer function model and the zero-pole-gain model. Extract zero-pole-gain data and numerator-denominator from the LTI object. Obtain the system dc gain and the pole-zero map of the transfer function. Obtain a minimum realization of this system. Use MATLAB to solve this problem. #### **1\.4.5.5 Example** Consider a system described by the following transfer function: G ( s ) \= Y ( s ) U ( s ) \= s 2 \+ 3 s \+ 3 s 3 \+ 6 s 2 \+ 11 s \+ 6 Find an analytical expression for an impulse response of the system. #### **1\.4.5.6 Example** A certain control system is described by the following transfer function: G ( s ) \= Y ( s ) U ( s ) \= 2 s \+ 8 s 3 \+ 5 s 2 \+ 8 s \+ 4 Find an analytical expression for a step response of the system. #### **1\.4.5.7 Example** A certain control system is described by the following closed loop transfer function: G c l ( s ) \= 45 ( s \+ 6 ) ( s 2 \+ 65 s \+ 354 ) s Find an analytical expression for a step response of the system โ€“ note the integrator term in the denominator\! #### **1\.4.5.8 Example** A certain control system is described by the following closed loop transfer function: G c l ( s ) \= 45 ( s \+ 6 ) s 3 \+ 20 s 2 \+ 129 s \+ 270 One of the closed loop poles is at -5. Find an analytical expression for a step response of the system. ``` ```
Shard0 (laksa)
Root Hash7622135113927760600
Unparsed URLca,torontomu!library,pressbooks,/controlsystems/chapter/1-4-laplace-transforms/ s443