Mon Apr 9 EP 7.6 Convolutions and Laplace transforms. Announcements: Warm-up Exercise:

Similar documents
Hint: There's a table of particular solutions at the end of today's notes.

e a s a f t dt f t dt = p = p. t = a

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring 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...

Chapter 2. First Order Scalar Equations

Announcements: Warm-up Exercise:

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

Some Basic Information about M-S-D Systems

Notes 04 largely plagiarized by %khc

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

KEY. Math 334 Midterm III Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

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

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow

Two Coupled Oscillators / Normal Modes

Mon Apr 2: Laplace transform and initial value problems like we studied in Chapter 5

Math 334 Fall 2011 Homework 11 Solutions

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.

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

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

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

ENGI 9420 Engineering Analysis Assignment 2 Solutions

Chapter 6. Systems of First Order Linear Differential Equations

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

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

6.003 Homework #9 Solutions

EECE 301 Signals & Systems Prof. Mark Fowler

Chapter Three Systems of Linear Differential Equations

8. Basic RL and RC Circuits

Math 333 Problem Set #2 Solution 14 February 2003

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

e 2t u(t) e 2t u(t) =?

t )? How would you have tried to solve this problem in Chapter 3?

ME 452 Fourier Series and Fourier Transform

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff

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.

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

Solutions to Assignment 1

6.2 Transforms of Derivatives and Integrals.

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

6.003 Homework #9 Solutions

EECE 301 Signals & Systems Prof. Mark Fowler

Section 7.4 Modeling Changing Amplitude and Midline

6.003 Homework #8 Solutions

Math 10B: Mock Mid II. April 13, 2016

t, x 2 t above, using input-output modeling. Assume solute concentration is uniform in each tank. If x 1 0 = b 1

ψ ( t) = c n ( t) t n ( )ψ( ) t ku t,t 0 ψ I V kn

Application 5.4 Defective Eigenvalues and Generalized Eigenvectors

M x t = K x F t x t = A x M 1 F t. M x t = K x cos t G 0. x t = A x cos t F 0

Echocardiography Project and Finite Fourier Series

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis

MA 366 Review - Test # 1

SOLUTIONS TO ECE 3084

Laplace Transform and its Relation to Fourier Transform

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

System of Linear Differential Equations

ψ ( t) = c n ( t ) n

δ (τ )dτ denotes the unit step function, and

Chapter 7: Solving Trig Equations

10. State Space Methods

Differential Equations

Math 116 Practice for Exam 2

Hamilton Jacobi equations

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

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

Solutions Problem Set 3 Macro II (14.452)

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

EECS 2602 Winter Laboratory 3 Fourier series, Fourier transform and Bode Plots in MATLAB

18.03SC Unit 3 Practice Exam and Solutions

ME 391 Mechanical Engineering Analysis

( ) = b n ( t) n " (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2.

1 Review of Zero-Sum Games

Math Final Exam Solutions

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

6.003: Signals and Systems. Relations among Fourier Representations

Math 315: Linear Algebra Solutions to Assignment 6

From Complex Fourier Series to Fourier Transforms

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

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

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

LAPLACE TRANSFORM AND TRANSFER FUNCTION

ES.1803 Topic 22 Notes Jeremy Orloff

2. Nonlinear Conservation Law Equations

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

MA Study Guide #1

6.003: Signal Processing

FITTING EQUATIONS TO DATA

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

The Laplace Transform

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

k B 2 Radiofrequency pulses and hardware

Solutions from Chapter 9.1 and 9.2

Matlab and Python programming: how to get started

Predator - Prey Model Trajectories and the nonlinear conservation law

Math 106: Review for Final Exam, Part II. (x x 0 ) 2 = !

EXERCISES FOR SECTION 1.5

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

ln 2 1 ln y x c y C x

Math Spring 2015 PRACTICE FINAL EXAM (modified from Math 2280 final exam, April 29, 2011)

Transcription:

Mah 225-4 Week 3 April 9-3 EP 7.6 - convoluions; 6.-6.2 - eigenvalues, eigenvecors and diagonalizabiliy; 7. - sysems of differenial equaions. Mon Apr 9 EP 7.6 Convoluions and Laplace ransforms. Announcemens: Warm-up Exercise:

The convoluion Laplace ransform able enry says i's possible o find he inverse Laplace ransform of a produc of Laplace ransforms. The answer is NOT he produc of he inverse Laplace ransforms, bu a more complicaed expression known as a "convoluion". If you've had a mulivariable Calculus class in which you sudied ineraed inegrals and changing he order of inegraion, you can verify ha his able enry is rue - I've included he proof as he las page of oday's noes. f g f g d F s G s convoluion inegrals o inver Laplace ransform producs Overview: This somewha amazing able enry allows us o wrie down soluion formulas for any consan coefficien nonhomogeneous differenial equaion, no maer how complicaed he righ hand side is. Le's focus our discussion on he sor of differenial equaion ha arises in Mah 225, namely he second order case x a x b x = f x = x x = v. Then in Laplace land, his equaion is equivalen o s 2 X s sx v a s X s x b X s = F s X s s 2 a s b = F s x s v a x. x s v a x X s = F s s 2 a s b s 2. a s b The inverse Laplace ransform of he second fracion conains he iniial value informaion, and is inverse Laplace ransform will be a homogeneous soluion for he differenial equaion, and will be zero if x = v =. (Noe ha he Chaper 5 characerisic polynomial is exacly p r = r 2 a r b, which coincides wih he denominaor p s = s 2 a s b.) The firs fracion is a produc of wo Laplace ransforms F s s 2 a s b = F s W s for W s s 2 a s b. and so we can use he convoluion able enry o wrie down an (inegral) formula for he inverse Laplace ransform. No maer wha forcing funcion f appears on he righ side of he differenial equaion, he corresponding soluion (o he IVP wih x = v = ) is always given by he inegral where x = f w = w f = w = W s. w f d,

Le's look more closely a ha soluion formula: The soluion o x a x b x = f x = x = is given by x = f w = w f = w f d = f w d. Exercise ) The funcion w = W s is called he "weigh funcion" for he differenial equaion. Verify ha i is a soluion o he homogeneous DE IVP w a w b w = w = w =. Exercise 2 Inerpre he convoluion formula x = f w d in erms of x being a resul of he forces f for, and how responsive he sysem hrough he weigh funcion, and he corresponding imes. This is relaed o a general principle known as "Duhamal's Principal", ha applies in a variey of seings. (See wikipedia.)

Exercise 3. Le's play he resonance game and pracice convoluion inegrals, firs wih an old friend, bu hen wih non-sinusoidal forcing funcions. We'll sick wih our earlier swing, bu consider various forcing periodic funcions f. x x = f x = x = a) Find he weigh funcion w. b) Wrie down he soluion formula for x as a convoluion inegral. c) Work ou he special case of X s when f = cos, and verify ha he convoluion formula reproduces he answer we would've goen from he able enry 2 k sin k s s 2 k 2 2

Exercise 4 The soluion x o x a x b x = f x = x = is given by x = f w = w f = w f d = f w d. We worked ou ha he soluion o our DE IVP will be Example ) A square wave forcing funcion wih ampliude and period 2. Le's alk abou how we came up wih he formula (which works unil = ). > wih plos : > x = sin f d Since he unforced sysem has a naural angular frequency =, we expec resonance when he forcing funcion has he corresponding period of T = 2 = 2. We will discover ha here is he possibiliy w for resonance if he period of f is a muliple of T. (Also, forcing a he naural period doesn' guaranee resonance...i depends wha funcion you force wih.) f 2 n = n Heaviside n Pi : ploa plo f, =..3, color = green : display ploa, ile = `square wave forcing a naural period` ; square wave forcing a naural period 2 3 ) Wha's your voe? Is his square wave going o induce resonance, i.e. a response wih linearly growing ampliude?

> x sin f d : plob plo x, =..3, color = black : display ploa, plob, ile = `resonance response?` ; resonance response? 2 3

Example 2) A riangle wave forcing funcion, same period > f2 f s ds.5 : # his aniderivaive of square wave should be riangle wave plo2a plo f2, =..3, color = green : display plo2a, ile = `riangle wave forcing a naural period` ; 2) Resonance? riangle wave forcing a naural period 2 3

> x2 sin f2 d : plo2b plo x2, =..3, color = black : display plo2a, plo2b, ile = `resonance response?` ; 5 5 5 5 resonance response? 2 3

Example 3) Forcing no a he naural period, e.g. wih a square wave having period T = 2. > f3 2 2 n = n Heaviside n : plo3a plo f3, =..2, color = green : display plo3a, ile = `periodic forcing, no a he naural period` ; 3) Resonance? periodic forcing, no a he naural period 5 5 2

> x3 sin f3 d : plo3b plo x3, =..2, color = black : display plo3a, plo3b, ile = `resonance response?` ;.5.5 resonance response? 5 5 2

Example 4) Forcing no a he naural period, e.g. wih a paricular wave having period T = 6. > f4 n = Heaviside 6 n Heaviside 6 n : plo4a plo f4, =..5, color = green : display plo4a, ile = sporadic square wave wih period 6 ; sporadic square wave wih period 6 5 5 4) Resonance?

> > x4 sin f4 d : plo4b plo x4, =..5, color = black : display plo4a, plo4b, ile = `resonance response?` ; 5 5 5 5 resonance response? 5 5

Hey, wha happened???? How do we need o modify our hinking if we force a sysem wih somehing which is no sinusoidal, in erms of worrying abou resonance? In he case ha his was modeling a swing (pendulum), how is i geing pushed? Precise Answer: I urns ou ha any periodic funcion wih period P is a (possibly infinie) superposiion P of a consan funcion wih cosine and sine funcions of periods P, 2, P 3, P,.... Equivalenly, hese 4 funcions in he superposiion are, cos, sin, cos 2, sin 2, cos 3, sin 3,... wih ω = 2. This is he P heory of Fourier series, which you will sudy in oher courses, e.g. Mah 35, Parial Differenial Equaions. If he given periodic forcing funcion f has non-zero erms in his superposiion for which n = (he naural angular frequency) (equivalenly P n = 2 = T ), here will be resonance; oherwise, no resonance. We could already have undersood some of his in Chaper 5, for example Exercise 5) The naural period of he following DE is (sill) T = 2 forcing funcion below is T = 6. Noice ha he period of he firs and ha he period of he second one is T = T = 2. Ye, i is he firs DE whose soluions will exhibi resonance, no he second one. Explain, using Chaper 5 superposiion ideas: a) x x = cos sin. 3 b) x x = cos 2 3 sin 3.