CHEE 319 Tutorial 3 Solutions. 1. Using partial fraction expansions, find the causal function f whose Laplace transform. F (s) F (s) = C 1 s + C 2

Similar documents
Math 353 Lecture Notes Week 6 Laplace Transform: Fundamentals

Laplace Transforms Chapter 3

CH.6 Laplace Transform

ENGIN 211, Engineering Math. Laplace Transforms

Lecture 7: Laplace Transform and Its Applications Dr.-Ing. Sudchai Boonto

Name: Solutions Final Exam

Math 266 Midterm Exam 2

A sufficient condition for the existence of the Fourier transform of f : R C is. f(t) dt <. f(t) = 0 otherwise. dt =

Introduction & Laplace Transforms Lectures 1 & 2

Homework 3 Solutions Math 309, Fall 2015

Honors Differential Equations

9.5 The Transfer Function

Control Systems. Frequency domain analysis. L. Lanari

2.161 Signal Processing: Continuous and Discrete Fall 2008

Laplace Transforms. Chapter 3. Pierre Simon Laplace Born: 23 March 1749 in Beaumont-en-Auge, Normandy, France Died: 5 March 1827 in Paris, France

Time Response of Systems

Chemical Engineering 436 Laplace Transforms (1)

Basic Procedures for Common Problems

Formulation of Linear Constant Coefficient ODEs

The Laplace Transform

HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS II: Nonhomogeneous Equations. David Levermore Department of Mathematics University of Maryland

Math 307 Lecture 19. Laplace Transforms of Discontinuous Functions. W.R. Casper. Department of Mathematics University of Washington.

Unit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace

MathQuest: Differential Equations

Math 216 Second Midterm 16 November, 2017

Math 216 Second Midterm 19 March, 2018

Differential Equations Practice: 2nd Order Linear: Nonhomogeneous Equations: Undetermined Coefficients Page 1

EE102 Homework 2, 3, and 4 Solutions

Lecture Discrete dynamic systems

Control Systems. Laplace domain analysis

APPM 2360: Midterm 3 July 12, 2013.

Laplace Transform Part 1: Introduction (I&N Chap 13)

Partial Fractions and the Coverup Method Haynes Miller and Jeremy Orloff

INC 341 Feedback Control Systems: Lecture 2 Transfer Function of Dynamic Systems I Asst. Prof. Dr.-Ing. Sudchai Boonto

Ordinary Differential Equations. Session 7

Name: Solutions Exam 3

Solutions to Math 53 Math 53 Practice Final

EE/ME/AE324: Dynamical Systems. Chapter 7: Transform Solutions of Linear Models

Solution: Homework 3 Biomedical Signal, Systems and Control (BME )

The Laplace Transform and the IVP (Sect. 6.2).

Transform Solutions to LTI Systems Part 3

Mathematics 3 Differential Calculus

JUST THE MATHS UNIT NUMBER LAPLACE TRANSFORMS 3 (Differential equations) A.J.Hobson

Queen s University at Kingston. CHEE Winter Process Dynamics and Control. M. Guay. Quiz 1

The Laplace Transform

Module 4. Related web links and videos. 1. FT and ZT

MA 201, Mathematics III, July-November 2016, Laplace Transform

1 Notes on Laplace circuit analysis

Control Systems I. Lecture 5: Transfer Functions. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

One-Sided Laplace Transform and Differential Equations

+ + LAPLACE TRANSFORM. Differentiation & Integration of Transforms; Convolution; Partial Fraction Formulas; Systems of DEs; Periodic Functions.

Solution of ODEs using Laplace Transforms. Process Dynamics and Control

CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang

Chapter 6: The Laplace Transform. Chih-Wei Liu

2.3 Oscillation. The harmonic oscillator equation is the differential equation. d 2 y dt 2 r y (r > 0). Its solutions have the form

MAT389 Fall 2016, Problem Set 11

Dynamic Response. Assoc. Prof. Enver Tatlicioglu. Department of Electrical & Electronics Engineering Izmir Institute of Technology.

Fourier Integral. Dr Mansoor Alshehri. King Saud University. MATH204-Differential Equations Center of Excellence in Learning and Teaching 1 / 22

Control Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS II: Nonhomogeneous Case and Vibrations

ELG 3150 Introduction to Control Systems. TA: Fouad Khalil, P.Eng., Ph.D. Student

ECEN 420 LINEAR CONTROL SYSTEMS. Lecture 2 Laplace Transform I 1/52

Section 6.4 DEs with Discontinuous Forcing Functions

Lecture 23 Using Laplace Transform to Solve Differential Equations

3.5 Undetermined Coefficients

Contents Partial Fraction Theory Real Quadratic Partial Fractions Simple Roots Multiple Roots The Sampling Method The Method of Atoms Heaviside s

37. f(t) sin 2t cos 2t 38. f(t) cos 2 t. 39. f(t) sin(4t 5) 40.

EXAM 2 MARCH 17, 2004

Solution via Laplace transform and matrix exponential

Ch 6.2: Solution of Initial Value Problems

Differential Equations, Math 315 Midterm 2 Solutions

Chapter 31. The Laplace Transform The Laplace Transform. The Laplace transform of the function f(t) is defined. e st f(t) dt, L[f(t)] =

Examination paper for TMA4130 Matematikk 4N: SOLUTION

Review Sol. of More Long Answer Questions

SOLUTIONS TO PRACTICE EXAM 3, SPRING 2004

Course roadmap. ME451: Control Systems. Example of Laplace transform. Lecture 2 Laplace transform. Laplace transform

spring mass equilibrium position +v max

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

Finding Transfer Functions of Circuits Using State-Space

f(t)e st dt. (4.1) Note that the integral defining the Laplace transform converges for s s 0 provided f(t) Ke s 0t for some constant K.

APPM 2360: Midterm exam 3 April 19, 2017

Periodic functions: simple harmonic oscillator

Lecture 7 - Separable Equations

GATE EE Topic wise Questions SIGNALS & SYSTEMS

The Method of Laplace Transforms.

Systems Engineering/Process Control L4

Chapter 10: Sinusoids and Phasors

Solutions for homework 5

7.2 Relationship between Z Transforms and Laplace Transforms

Work sheet / Things to know. Chapter 3

COMPLEX NUMBERS AND DIFFERENTIAL EQUATIONS

Response to a pure sinusoid

Step Response Analysis. Frequency Response, Relation Between Model Descriptions

{ sin(t), t [0, sketch the graph of this f(t) = L 1 {F(p)}.

ELEC2400 Signals & Systems

Sample Quiz 8, Problem 1. Solving Higher Order Constant-Coefficient Equations

e st f (t) dt = e st tf(t) dt = L {t f(t)} s

Answers to Problem Set Number MIT (Fall 2005).

Short Solutions to Review Material for Test #2 MATH 3200

George Mason University Signals and Systems I Spring 2016

Transcription:

CHEE 39 Tutorial 3 Solutions. Using partial fraction expansions, find the causal function f whose Laplace transform is given by: F (s) 0 f(t)e st dt (.) F (s) = s(s+) ; Solution: Note that the polynomial s(s + ) in the denominator of F (s) is of the form n i= (s p i), where the p i are all distinct. We can therefore use the so-called cover-up method. We want to write s + C s +. Let s find the constants C and C. We start with C. We have s + sc s +. It follows that C = (sf (s)), so that s C = s(s + ) = (s + ) = =. Similarly, we have (s + )F (s) = (s + )C s + C. It follows that C = ((s + )F (s)) s=, so that (s + ) C = s(s + ) = s= s =. s= Hence, F (s) = s s +. Using the Laplace table, as well as linearity of the inverse Laplace transform, we see that f is the causal function { e t, t 0 (.) F (s) = 0 s(s+)(s+0) ;

Solution: Since the polynomial s(s + )(s + 0) in the denominator of F (s) is of the form n i= (s p i), where the p i are all distinct, we can use the cover-up method as we did in (.). We want to write s + C s + + C 3 s + 0. Let s find the constants C, C, and C 3, starting with C. We have s + sc s + + sc 3 s + 0. It follows that C = (sf (s)), so that 0s C = s(s + )(s + 0) = 0 (s + )(s + 0) =. Similarly, for C we have (s + )F (s) = (s + )C s + C + (s + )C 3 s + 0, It follows that C = ((s + )F (s)) s=, so that 0(s + ) C = 0 s(s + )(s + 0) = s= s(s + 0) = 0 s= 9. Finally, for C 3 we have (s + 0)F (s) = (s + 0)C s + (s + 0)C s + + C 3. It follows that C 3 = ((s + 0)F (s)) s= 0, so that 0(s + 0) C 3 = 0 s(s + )(s + 0) = s= 0 s(s + ) = 0 s= 0 90 = 9. Hence, F (s) = s 0 9 s + + 9 s + 0. As in (.), using the Laplace table, as well as linearity of the inverse Laplace transform, we see that f is the causal function { 0 9 e t + 9 e 0t, t 0 (.3) F (s) = 3s+ s +4s+0 ; Solution: Denoting + 4i by z, we see that the denominator of F (s) factors as s + 4s + 0 = (s + z)(s + z). (Here the overline denotes complex conjugation. That is, if

z = a + bi, then z a bi.) Since z, z are distinct, we can proceed as in (.) and (.). We want to write s + z + C s + z. Let s find C. We have C = ((s + z)f (s)), so that (s + z)(3s + ) C = (s + z)(s + z) = 3s + s + z = 3z + z + z or C = 3/ i/. Similarly, C = ((s + z)f (s)), so that (s + z)(3s + ) C = (s + z)(s + z) = 3s + s + z = 3z + z + z or C = 3/ + i/. Hence = 3( + 4i) + 4i + 4i = 3( 4i) + + 4i + + 4i 4 i =, 8i 4 + i =, 8i F (s) = 3/ i/ s + z + 3/ + i/. s + z As in (.) and (.), it follows that f is the function given by (3/ i/)e zt + (3/ + i/)e zt = (3/ i/)e ( 4i)t + (3/ + i/)e ( +4i)t = e t 3/ (e i4t + e i4t ) i/ (e i4t e i4t )] = e t 3 ei4t + e i4t i e i4t e i4t ] = e t 3 ei4t + e i4t + i ei4t e i4t ] = e t 3 ei4t + e i4t = 3e t cos(4t) e t sin(4t) ei4t e i4t i when t 0. (Here we have used Euler s identity for the sine and cosine functions in the last line.) Hence f is the causal function OR Solution : Noting the Laplace transform pairs { 3e t cos(4t) e t sin(4t), t 0 e at sin ωt ω (s + a) + ω e at cos ωt s + a (s + a) + ω ] 3

and completing the square on the denominator of F (s): we see that F (s) = 3s + F (s) = (s + ) + 4 = 3 ] s + = 3 (s + ) + 4 and the solution f(t) follows. 3s + (s + ) + 4 s + /3 (s + ) + 4 ] 4 (s + ) + 4 ] s + 4/3 = 3 (s + ) + 4 (.4) F (s) = 3s +9s+ (s+)(s +5s+) ; Solution: The denominator of F (s) factors as (s + )(s + 5s + ) = (s + )(s + z)(s + z), where z 5 + i 9. Since, z, z are distinct, we can proceed as in (.) (.3). We want to write s + + C s + z + C 3 s + z. Let s find C. We have C = ((s + )F (s)) s=, so that C = (s + )(3s + 9s + ) (s + )(s + 5s + ) = 3s + 9s + s= s + 5s + = 6 s= 5. For C, we have C = ((s + z)f (s)), so that C = (s + z)(3s + 9s + ) (s + )(s + z)(s + z) = 3s + 9s + (s + )(s + z) = 9 0 5 90 i 9. Similarly, C 3 = ((s + z)f (s)), so that C 3 = (s + z)(3s + 9s + ) (s + )(s + z)(s + z) = 3s + 9s + (s + )(s + z) = 9 0 + 5 90 i 9. It follows that f is the function given by 6 ( 9 5 e t + 0 5 ) ( 9 90 i 9 e zt + 0 + 5 = 6 ( 9 5 e t + 0 5 ) 90 i 9 e ( 5/ i 9/)t + ) 90 i 9 e zt ( 9 0 + 5 90 i 9 ) e ( 5/+i 9/)t = 6 9 5 e t + e 5t/ 0 (ei 9t/ + e i 9t/ ) + 5 90 i 9(e i 9t/ e )] i 9t/ = 6 5 e t + e 5t/ 9 5 ei 9t/ + e i 9t/ 5 ] 9 ei 9t/ e i 9t/ 95 i = 6 5 e t + 9 5 e 5t/ cos( 9t/) 5 9 e 5t/ sin( 9t/) 95 4

when t 0. Hence f is the causal function { 6 5 e t + 9 5 e 5t/ cos( 9t/) 5 9 95 e 5t/ sin( 9t/), t 0 (.5) F (s) =. s 6 Solution: We can simply look this one up in the table, which tells us that f is the causal function { t 5 5!, t 0. Find the causal function f whose Laplace transform is given by: (.) F (s) = s(s+) ; Solution: Note that the denominator s(s + ) of F (s) cannot be written in the form n i= (s p i), where the p i are all distinct. As we have seen in the lectures, we must make use of Heaviside s theorem to find the partial fraction expansion of F (s). We want to write s + C s + + C 3 (s + ). The constants C, C and C 3 are given by C = (sf (s)) = (s + ) = 4 ; C = d ds (s + ) F (s) = d s= ds s= s = 4 ; C 3 = ((s + ) F (s)) s= = s = s=, so that F (s) = /4 s /4 s + / (s + ) Consulting the table, we see that f is the causal function { 4 4 e t te t, t 0 (.) F (s) = (s+)(s+5) (s+)(s +4) ; Solution: As in (.), it is clear from looking at the denominator of F (s) that we must make use of Heaviside s theorem. The denominator of F (s) can be written as (s+)(s + 4) = (s + )(s + 4)(s + 4) = (s + )(s + i) (s i). We want to write s + + C s + i + C 3 (s + i) + C 4 s i + C 5 (s i). 5

The constants C, C, C 3, C 4 and C 5 are given by (s + )(s + 5) C = ((s + )F (s)) s= = (s + 4) = 3 s= 5 ; C = d ds (s + i) F (s) = 6 s= i 5 + 579 00 i; C 3 = ((s + i) (s + )(s + 5) F (s)) s= i = (s + )(s i) = 83 s= i 0 + 39 0 i; C 4 = d ds (s i) F (s) = 6 s=i 5 579 00 i; C 5 = ((s i) (s + )(s + 5) F (s)) s=i = (s + )(s + i) = 83 s=i 0 39 0 i, so that F (s) = 3/5 6/5 + 579i/00 + + s + s + i 83/0 39i/0 + (s i). 83/0 + 39i/0 6/5 579i/00 (s + i) + s i After performing some straightforward manipulations as in (.4), and consulting the table, we see that f is the function given by 3 5 e t 3 ( 579 cos(t) + 5 00 sin(t) + 83 0 + 39 ) ( 0 i te it + 83 0 39 ) 0 i te it when t 0. fashion as Using Euler s identity again, we can write this in a more recognizable 3 5 e t 3 579 cos(t) + 5 00 Hence f is the causal function { 3 5 e t 3 5 83 39 sin(t) t cos(t) 0 0 t sin(t). 579 83 39 cos(t) + 00 sin(t) 0t cos(t) 0 t sin(t), t 0 (.3) F (s) = s (s +). Solution: Proceeding as in (.), we write the denominator as (s + ) = (s + i) (s i). We want to write s + i + C (s + i) + C 3 s i + C 4 (s i). 6

The constants C, C, C 3 and C 4 are given by C = d ds (s + i) F (s) = 0; s= i C = ((s + i) F (s)) s= i = s (s i) = s= i ; C 3 = d ds (s i) F (s) = 0; s=i C 4 = ((s i) F (s)) s=i = s (s + i) = s=i, so that F (s) = / (s + i) + / (s i). Consulting the table, we see that f is the function given by te it + ( e it + e it ) teit = t = t cos(t) when t 0. Hence f is the causal function { t cos(t), t 0 3. Solve the following scalar initial-value problems using the Laplace transform: (3.) ÿ(t) + ẏ(t) + 4y(t) = 0, with initial conditions y(0) = and ẏ(0) = ; Solution: Taking the Laplace transform of both sides, and denoting the Laplace transform of y by Y, we have s Y (s) sy(0) ẏ(0) + sy (s) y(0) + 4Y (s) = 0. Plugging in the initial conditions, we obtain Isolating Y (s), we obtain where s Y (s) s + sy (s) + 4Y (s) = 0. Y (s) = s + 3 s + s + 4 = s + 3 (s + z)(s + z), z + 5 i. 7

Let s use a partial fraction expansion to find y, as we did in Question. Obviously we can use the cover-up method. We want to write Y (s) = C s + z + C s + z. As in Question, we can perform straightforward calculations to show that the constants C and C are given by and that C = ((s + z)y (s)) = 5 6 i; C = ((s + z)y (s)) = + 5 6 i, y(t) = e t/ cos( 5t/) 5 3 e t/ sin( 5t/). (3.) ÿ(t) + ẏ(t) = t, with initial conditions y(0) = and ẏ(0) =. Solution: Taking the Laplace transform of both sides as in (3.), we have s Y (s) sy(0) ẏ(0) + sy (s) y(0) = s. Plugging in the initial conditions, we obtain s Y (s) s + sy (s) = s. Isolating Y (s) and rearranging terms, we obtain Y (s) = s + s s + s = s3 + s 3 (s + ) = s + + s 3 (s + ). For the partial fraction expansion of the second term of Y (s), we want to write s 3 (s + ) = C s + C s + C 3 s 3 + C 4 s +. The constants C, C, C 3 and C 4 are given by C = d ds s + = (s + ) 3 = ; C = d ds s + = (s + ) = ; C 3 = s + = ; C 4 = ((s + ) s 3 (s + ) ) = s= s s= 3 =, 8

so that Consulting the table, we see that Y (s) = s + + s s + s 3 s + = s s + s 3. y(t) = t + t. 9