LECTURE 18: THE JOY OF CONVOLUTION

Similar documents
ECE-202 FINAL April 30, 2018 CIRCLE YOUR DIVISION

Ordinary differential equations

Lecture 9 Time-domain properties of convolution systems

Lecture 13: H(s) Poles-Zeros & BIBO Stability. (a) What are poles and zeros? Answer: HS math and calculus review.

Convolution and Linear Systems

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

MAE143 A - Signals and Systems - Winter 11 Midterm, February 2nd

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

06/12/ rws/jMc- modif SuFY10 (MPF) - Textbook Section IX 1

GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM. COURSE: ECE 3084A (Prof. Michaels)

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

x(t) = t[u(t 1) u(t 2)] + 1[u(t 2) u(t 3)]

Chapter 2: Linear systems & sinusoids OVE EDFORS DEPT. OF EIT, LUND UNIVERSITY

Figure 1 A linear, time-invariant circuit. It s important to us that the circuit is both linear and time-invariant. To see why, let s us the notation

Volterra Series. Nick Gamroth. January 2005

Ordinary Differential Equations

Lecture 4 : Transform Properties and Interpretations. Continued to the Next (Higher) Level. 1. Example 1. Demo of the mult-by-t property.

ENGIN 211, Engineering Math. Laplace Transforms

2 Classification of Continuous-Time Systems

EE Control Systems LECTURE 9

LTI Systems (Continuous & Discrete) - Basics

Problem Weight Score Total 100

EE-202 Exam II March 3, 2008

Digital Signal Processing Prof. T. K. Basu Department of Electrical Engineering Indian Institute of Technology, Kharagpur

Solving a RLC Circuit using Convolution with DERIVE for Windows

(f g)(t) = Example 4.5.1: Find f g the convolution of the functions f(t) = e t and g(t) = sin(t). Solution: The definition of convolution is,

EECE 3620: Linear Time-Invariant Systems: Chapter 2

Dr. Ian R. Manchester

PROFESSOR: WELCOME BACK TO THE LAST LECTURE OF THE SEMESTER. PLANNING TO DO TODAY WAS FINISH THE BOOK. FINISH SECTION 6.5

20.6. Transfer Functions. Introduction. Prerequisites. Learning Outcomes

1.17 : Consider a continuous-time system with input x(t) and output y(t) related by y(t) = x( sin(t)).

Signal Processing and Linear Systems1 Lecture 4: Characterizing Systems

PHYSICS 107. Lecture 5 Newton s Laws of Motion

Computing inverse Laplace Transforms.

2.161 Signal Processing: Continuous and Discrete Fall 2008

1. SINGULARITY FUNCTIONS

EE102 Homework 2, 3, and 4 Solutions

MITOCW watch?v=rf5sefhttwo

Chapter 3 Convolution Representation

EE-202 Exam III April 6, 2017

Unit 2: Electrostatics

Homework Solutions:

p(t)dt a p(τ)dτ , c R.

Time Response of Systems

Midterm 1 Review. Distance = (x 1 x 0 ) 2 + (y 1 y 0 ) 2.

HW 9, Math 319, Fall 2016

Instructor (Brad Osgood)

Continuous-Time Frequency Response (II) Lecture 28: EECS 20 N April 2, Laurent El Ghaoui

TIME: 45 minutes. LESSON: Curious About Clouds GRADE: 1 st SUMMARY:

CH.3 Continuous-Time Linear Time-Invariant System

Lecture 29. Convolution Integrals and Their Applications

1 otherwise. Note that the area of the pulse is one. The Dirac delta function (a.k.a. the impulse) can be defined using the pulse as follows:

Discussion Section #2, 31 Jan 2014

Communication Engineering Prof. Surendra Prasad Department of Electrical Engineering Indian Institute of Technology, Delhi

Math Shifting theorems

GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM. COURSE: ECE 3084A (Prof. Michaels)

ASTRO 114 Lecture Okay. We re now gonna continue discussing and conclude discussing the entire

Notes 17 largely plagiarized by %khc

The Convolution Sum for Discrete-Time LTI Systems

Predator - Prey Model Trajectories and the nonlinear conservation law

Using MATLAB with the Convolution Method

First Derivative Test

Systems Analysis and Control

MITOCW big_picture_derivatives_512kb-mp4

Laplace Transforms and use in Automatic Control

Linear dynamical systems with inputs & outputs

LESSON 21: DIFFERENTIALS OF MULTIVARIABLE FUNCTIONS OCTOBER 18, 2017

Several Things to Talk about

Frequency Response (III) Lecture 26:

Introduction & Laplace Transforms Lectures 1 & 2

LESSON 21: DIFFERENTIALS OF MULTIVARIABLE FUNCTIONS MATH FALL 2018

a_n x^(n) a_0 x = 0 with initial conditions x(0) =... = x^(n-2)(0) = 0, x^(n-1)(0) = 1/a_n.

Linear Filters and Convolution. Ahmed Ashraf

3. Solutions of first order linear ODEs

1 Work, Power, and Machines

L2 gains and system approximation quality 1

Therefore the new Fourier coefficients are. Module 2 : Signals in Frequency Domain Problem Set 2. Problem 1

MITOCW watch?v=hdyabia-dny

Convolution. Define a mathematical operation on discrete-time signals called convolution, represented by *. Given two discrete-time signals x 1, x 2,

Amy Ant. Formica Mica Grant. Dr. Sheila Grant (Mica s mom)

Nonlinear Control Lecture # 14 Input-Output Stability. Nonlinear Control

MITOCW R9. Generalized Forces

Generalized sources (Sect. 6.5). The Dirac delta generalized function. Definition Consider the sequence of functions for n 1, Remarks:

4 th Grade PSI. Slide 1 / 107 Slide 2 / 107. Slide 3 / 107. Slide 4 / 107. Slide 5 / 107. Slide 6 / 107. The History of Planet Earth

Lecture 8 ELE 301: Signals and Systems

Today s goals So far Today 2.004

Lecture 6.1 Work and Energy During previous lectures we have considered many examples, which can be solved using Newtonian approach, in particular,

Warm-up: Are accuracy and precision the same thing? (If so do you want to bet the house on it?)

{ }. The dots mean they continue in that pattern.

2.161 Signal Processing: Continuous and Discrete Fall 2008

Lecture 18: Reflection and Impedance

Math 3313: Differential Equations Laplace transforms

LECTURE 31: Basics of Magnetically Coupled Circuits Part 1

Signals, Instruments, and Systems W5. Introduction to Signal Processing Sampling, Reconstruction, and Filters

MITOCW ocw-5_60-s08-lec19_300k

Chapter 6: The Laplace Transform. Chih-Wei Liu

Talk Science Professional Development

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 2

Algebra. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.

Methods of Mathematics

Transcription:

LECTURE 18: THE JOY OF CONVOLUTION Part 1: INTRO 1. Just because no one understands you, doesn t mean you re an artist. Rich Hebda. I am reminded of this quote every time I talk about convolution. 2. WORDS: CONVOLUTED: folded in curved or tortuous windings. CONVOLVE: to role together to writh (in agony over HW problems). PRESCRIPTION: read the notes first. 3. RECALL LECTURE 1. Now you will understand why we avoided staying in the t-world and did intensive counseling in the s-world to solve our circuits problems. Input Signal (voltage/current source) Circuit /Model h(t) Impulse Response Output Signal (voltage/current)!!! L Input Signal L Circuit /Model H(s) Transfer Function L Output Signal

Lecture 17 Sp 18 2 R. A. DeCarlo COMMENT: Convolution allows us to traverse the top part of the diagram for a much larger class of signals than the one-sided Laplace transform what would one expect since the transform is SOOOO one-sided. DEFINITION: Let f (t) be an input signal possibly nonzero for t < and y(t) the circuit output signal. Suppose h(t) is the impulse response of the circuit. THEN y(t) is the convolution of h(t) with f (t), denoted as y(t) = h(t)* f (t) = f (t)* h(t), and defined by the integral y(t) = h(t τ ) f (τ )dτ = h(τ ) f (t τ )dτ 4. What is the meaning of the convolution integral? y(t) = h(t τ ) f (τ )dτ (i) t is the time at which the output y(t) is observed/measured. (ii) τ is the time when a piece of the input signal is applied. (iii) t τ is the time that has elapsed between the application of the input signal and the observation/measurement of the output signal.

Lecture 17 Sp 18 3 R. A. DeCarlo (iv) h(t τ ) measures the effect that a nano piece of the input signal at time τ HAS on the output at time t, i.e., the effect at time t due to the input f (τ ) is precisely h(t τ ) f (τ ). (v) By superposition, the integral sums up each and every little crumb of a contribution. Comment: Imagine ants climbing out of their nest at the dining hall and crawling along the broken road of glass, cake crumbs, spills, footprints, and dirt. Some ants get lost along the way. But every student is but a northern star, a sign, pointing straight to the loving arms of the outside sunshine and warm temperatures (somewhere in Florida). So only a handful, some fraction, say h(t τ ), of the original nest escapees wipe their brow (and antennae) and push on through as they march on home into. (Apologies to Rascal Flats.) PART 2: GENERAL EXAMPLES Example 1. y(t) = h(t)*[k δ (t T )]. From the sifting property of the delta function, y(t) = h(t)*[k δ (t T )] = K h(t T ).

Lecture 17 Sp 18 4 R. A. DeCarlo Example 2. Show that r(t) = u(t)*u(t). (VERY IMPPORTANT) Step 1. Write down the convolution integral and make a simple observation. u(t)*u(t) = since u(τ ) = 1 and is nonzero only for τ. u(t τ )u(τ )dτ = u(t τ )dτ Step 2. My dear Watson, consider the case of the negative t, i.e., t <. From step 1, τ. And so, t τ < implies u(t τ ) =, and so for t <, y(t) =. Step 3. My dear Watson, now consider the case of the non-negative t, i.e., t. From step 1, τ. And so, t τ only when τ t and so u(t τ ) = 1 only when τ t in which case: for t, u(t)*u(t) = t u(t τ )dτ = tu(t) = r(t) EOE Student: Hey Professor Ray, a question. Ray: Yes.

Lecture 17 Sp 18 5 R. A. DeCarlo Student: Ya know, it is a heck of a lot easier to use the bottom half of your diagram since L{ u(t)*u(t) } = L{ u(t) } L{ u(t) } = 1 s 2 and 1 L 1 Ray: Well, I was just trying to see if you were paying attention. s 2 = r(t). Student: Yeah Right Professor Ray. Now if I were you and you were me, I would have said that I was trying to teach you how to solve a convolution integral even though it is easier to solve the problem in the s-world. You were always telling us how easy it is to solve things in the s-world, right? Ray: Ha ha. Some Properties of Convolution: 1. Commutativety: h* f = f * h 2. Associativity: h*( f * g) = (h* f )* g 3. Distributive: h*( f + g) = h* f + h* g 4. TIME SHIFT PROPERTY (VERY VERY IMPPORTANT): if y(t) = h(t)* f (t) then h( t T 1 )* f ( t T 2 ) = y( t T 1 T 2 ).

Lecture 17 Sp 18 6 R. A. DeCarlo Example 3. Compute y(t) = h(t)* f (t) where h(t) and f (t) are shown in the figures below. 2 1.8 1.6 1.4 f(t) blue and h(t) red 1.2 1.8.6.4.2 1.5 1.5.5 1 1.5 2 2.5 time Strategy: combine the results of Example 2 with the properties of convolution, especially the time shift property. y(t) = h(t)* f (t) = f (t)* h(t) = u(t +1)* 2u(t) u(t 1) u(t 2) u(t 1)* 2u(t) u(t 1) u(t 2) = 2r(t +1) r(t) r(t 1) 2r(t 1) r(t 2) r(t 3) = 2r(t +1) r(t) 3r(t 1) + r(t 2) + r(t 3)

Lecture 17 Sp 18 7 R. A. DeCarlo 3 2.5 2 y(t)=h(t)*f(t) 1.5 1.5 2 1 1 2 3 4 Time Example 4. Compute y(t) = h(t)* f (t) where h(t) = u(t) u(t T ), T > and f (t) = u( t), an input from the past. Solution. Step 1. Again y(t) = h(t)* f (t) = f (t)* h(t): Step 2. Small details. y(t) = u(τ t) u(τ ) u(τ T ) dτ u(τ ) u(τ T ) = 1 τ < T otherwise Thus

Lecture 17 Sp 18 8 R. A. DeCarlo T y(t) = u(τ t) dτ Step 3. Case 1: t <, τ < T, implies u(τ t) = 1 implies that for t <, y(t) = T dτ = T Step 4. Case 2: t < T, τ < T implies u(τ t) = 1 τ t otherwise y(t) = T t dτ = T t Step 5. Case 3: t > T, τ < T, implies τ < t implies τ t < implies u(τ t) = implies y(t) =. Step 6. Summary: y(t) = T t < T t t < T T t The plot below is obtained by setting T = 1.

Lecture 17 Sp 18 9 R. A. DeCarlo 1.4 1.2 1 u( t)*[u(t) u(t T), T = 1.8.6.4.2 2 1 1 2 3 4 Time Example 5. Compute y(t) = h(t)* f (t) where h(t) = e t u( t), an impulse response from the first age of middle earth, and f (t) = e t u(t), Step 1. Write down the specific convolution as per the definition of convolution: y(t) = e t τ u(τ t)e τ u(τ ) dτ = e t e 2τ u(τ t) dτ since u(τ ) = for τ <.

Lecture 17 Sp 18 1 R. A. DeCarlo Step 2. Break the integral up into cases. Case 1: t <, τ. (i) u(τ t) = 1 and (ii) y(t) = e t e 2τ dτ = e t e 2τ 2 =.5e t, t < Case 2: t, τ : (i) u(τ t) = 1 τ t t > τ and (ii) y(t) = e t e 2τ u(τ t) dτ = e t e 2τ dτ = e t e 2τ 2 t t =.5e t, t Step 3. Summary: y(t) =.5e t t <.5e t t =.5e t u( t) +.5e t u(t)

Lecture 17 Sp 18 11 R. A. DeCarlo.5.45.4.35 The Convolution.3.25.2.15.1.5 3 2 1 1 2 3 time Student: Hey Professor Ray, this convolution stuff is like a 5% off sale at Sofa King Cool.