Convolution. Lecture #6 2CT.3 8. BME 333 Biomedical Signals and Systems - J.Schesser

Similar documents
Signals and Systems Linear Time-Invariant (LTI) Systems

Continuous Time Linear Time Invariant (LTI) Systems. Dr. Ali Hussein Muqaibel. Introduction

Signals and Systems Profs. Byron Yu and Pulkit Grover Fall Midterm 1 Solutions

h[n] is the impulse response of the discrete-time system:

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

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

EE 313 Linear Signals & Systems (Fall 2018) Solution Set for Homework #8 on Continuous-Time Signals & Systems

Theory of! Partial Differential Equations!

Chapter 8 The Complete Response of RL and RC Circuits

Theory of! Partial Differential Equations-I!

SOLUTIONS TO ECE 3084

The Laplace Transform

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

ME 452 Fourier Series and Fourier Transform

Homework 4 SOLUTION EE235, Summer 2012

LAPLACE TRANSFORM AND TRANSFER FUNCTION

EE 301 Lab 2 Convolution

5. Response of Linear Time-Invariant Systems to Random Inputs

Notes 04 largely plagiarized by %khc

Lecture #6: Continuous-Time Signals

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler

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

Representing a Signal. Continuous-Time Fourier Methods. Linearity and Superposition. Real and Complex Sinusoids. Jean Baptiste Joseph Fourier

Solutions - Midterm Exam

EE 224 Signals and Systems I Complex numbers sinusodal signals Complex exponentials e jωt phasor addition

Ordinary dierential equations

Homework 6 AERE331 Spring 2019 Due 4/24(W) Name Sec. 1 / 2

DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2008

Linear Time-invariant systems, Convolution, and Cross-correlation

MAE143A Signals & Systems - Homework 2, Winter 2014 due by the end of class Thursday January 23, 2014.

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

Laplace Transforms. Examples. Is this equation differential? y 2 2y + 1 = 0, y 2 2y + 1 = 0, (y ) 2 2y + 1 = cos x,

9/9/99 (T.F. Weiss) Signals and systems This subject deals with mathematical methods used to describe signals and to analyze and synthesize systems.

in Engineering Prof. Dr. Michael Havbro Faber ETH Zurich, Switzerland Swiss Federal Institute of Technology

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

EECS20n, Solution to Midterm 2, 11/17/00

Chapter 1 Fundamental Concepts

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

System Processes input signal (excitation) and produces output signal (response)

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL

Laplace Transform and its Relation to Fourier Transform

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

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

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

THE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI

4/9/2012. Signals and Systems KX5BQY EE235. Today s menu. System properties

8. Basic RL and RC Circuits

Signal and System (Chapter 3. Continuous-Time Systems)

F This leads to an unstable mode which is not observable at the output thus cannot be controlled by feeding back.

Matlab and Python programming: how to get started

10. State Space Methods

CE 395 Special Topics in Machine Learning

EECE.3620 Signal and System I

Characteristics of Linear System

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Final Exam Solutions

6.2 Transforms of Derivatives and Integrals.

Calculus Tricks #1. So if you understand derivatives, you ll understand the course material much better. a few preliminaries exponents

Lecture 2: Optics / C2: Quantum Information and Laser Science

6.003 Homework #13 Solutions

Topic Astable Circuits. Recall that an astable circuit has two unstable states;

EECE 301 Signals & Systems Prof. Mark Fowler

KEEE313(03) Signals and Systems. Chang-Su Kim

Linear Circuit Elements

Communication Systems, 5e

10/10/2011. Signals and Systems EE235. Today s menu. Chicken

2 int T. is the Fourier transform of f(t) which is the inverse Fourier transform of f. i t e

Solutions to Assignment 1

Delhi Noida Bhopal Hyderabad Jaipur Lucknow Indore Pune Bhubaneswar Kolkata Patna Web: Ph:

EECE 301 Signals & Systems Prof. Mark Fowler

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...

6.003: Signals and Systems. Relations among Fourier Representations

6.003 Homework #8 Solutions

6.003 Homework #9 Solutions

Module 4: Time Response of discrete time systems Lecture Note 2

Q1) [20 points] answer for the following questions (ON THIS SHEET):

11/16/99 (T.F. Weiss) Lecture #18: Continuous time periodic signals. Outline: Fourier series of periodic functions. Motivation:

on the interval (x + 1) 0! x < ", where x represents feet from the first fence post. How many square feet of fence had to be painted?

6.003 Homework #9 Solutions

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal?

Tutorial Sheet #2 discrete vs. continuous functions, periodicity, sampling

Time Domain Transfer Function of the Induction Motor

6/27/2012. Signals and Systems EE235. Chicken. Today s menu. Why did the chicken cross the Möbius Strip? To get to the other er um

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

+ γ3 A = I + A 1 (1 + γ + γ2. = I + A 1 ( t H O M E W O R K # 4 Sebastian A. Nugroho October 5, 2017

6.003: Signals and Systems

Lectures 29 and 30 BIQUADRATICS AND STATE SPACE OP AMP REALIZATIONS. I. Introduction

Sub Module 2.6. Measurement of transient temperature

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

Unsteady Flow Problems

Math 333 Problem Set #2 Solution 14 February 2003

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

Chapter 4 The Fourier Series and Fourier Transform

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

Basic definitions and relations

ADDITIONAL PROBLEMS (a) Find the Fourier transform of the half-cosine pulse shown in Fig. 2.40(a). Additional Problems 91

6.003: Signal Processing

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Transcription:

Convoluion Lecure #6 C.3 8

Deiniion When we compue he ollowing inegral or τ and τ we say ha he we are convoluing wih g d his says: ae τ, lip i convolve in ime -τ, hen displace i in ime by seconds -τ, and muliply i by τ. Finally, inegrae he produc over all τ.

Properies o Convoluion Firs some shorhand: Commuaive: Associaive: Disribuive: d ] [ ] [ 3 3 ] [ ] [ ] [ 3 3

A Graphical Example o How o Perorm a Convoluion =A- τ/ A B hτ C h d h-τ B h-τ B τ τ - τ -- τ We need o loo a 5 cases: <, < <, 3 < <, 4 < < 3, & 5 > 3 For case and case 5, C= since here is no overlap. B h-τ, < A h-τ, >3 A B -- -- τ 3 τ 3

Graphical Convoluion Example Coninued Case : < < Case 3: < < Case 4: < < 3 AB/ 3 h-τ, < < B A -- h-τ, < < 3 A B -- 3 h-τ, < < B A -- C τ C τ τ AB d AB AB / AB d C AB[ AB / AB AB d ] 4 AB

Malab Code clear all; endpulse=; s=.; endpoin=; n=-endpoin:s:endpoin; nn=-endpoin*:s:endpoin*; pulser=n>=&n<=endpulse; pulse=*pulser; ripulse=n>=&n<=; ri=-n.*ripulse; subplo,, plon,ri,'r',n,pulse,'b'; ile'signals'; xlabel'ime Seconds'; axis[- min[minri minpulse].*max[maxri maxpulse]]; c=convpulse,ri*s; subplo,, plonn,c; ile'convoluion'; xlabel'shi Seconds'; axis[- minc.*maxc];.5-3 4 5 6 7 8 9 ime Seconds Convoluion.4. Signals - 3 4 5 6 7 8 9 Shi Seconds 5

Inegraion o Convoluion Inegral Case : C AB d AB AB { } AB { } AB { } AB { } AB{ } Case 3: From he inegraion rom Case C AB d AB AB{ } AB{ } AB Case 4: 3 From he inegraion rom Case C AB d AB AB { } AB { } AB { } AB { } AB{ } 6

7 Convoluion and Sysems For an LI sysem, le s deine h as he sysem response o a uni impulse source, δ. hen he ollowing mus be rue: superposiion scalar ime invariance h x x h x x h h y x

Convoluion and Sysems Coninued Consruc x as he sum o uni impulse slices. he irs expression represens he slices o he source oaled as approaches ininiy. Le Side Equals lim x approaches x d x x x h x x h d y Consruc y as he sum o slices o he response due o an uni impulse uncion: xδh-δ. he inegral on he righ is he convoluion o x and h. Righ Side Equals lim x h approaches x h d his resul is very imporan since i says ha i one nows he impulse response o a sysem hen he oupu response or any given inpu source can be ound by convolving he inpu wih he impulse response. 8

Impulse Response and Causaliy y o does no depend on x i ha occur a imes aer o, i > o. y h x d h x d h x d o o o o For his inegral h x d, is posiive,, and or x x we have i o i o o o o For his inegral h x d, is negaive,, and or x x we have i o i o hereore or y o be causal his h x d mus be zero. For his o happen h,. o o o o 9

Calculaing he Uni Impulse Response, h R Vs i L Le s irs loo a mehods:. Narrow Pulse approximaion. Diereniaing u Beer Mehods o Come 3

Narrow Pulse Approximaion ε = widh o pulse ε heigh o pulse As ε, he narrow pulse = /ε [u-u-ε] δ So hen le s irs loo a he response o u: di Vs u i R L d i = Vs/R -e -R/L u Now we consruc he narrow pulse response: i = Vs/R /ε {-e -R/L u - e -R/L-ε u-ε} 3

Narrow Pulse Approximaion Coninued i = Vs/R/ε [-e -R/L u -e -R/L-ε u-ε] OR, or < i= Vs/R/ε [-e -R/L ] or < ε Vs/R/ε [-e -R/L ] -e -R/L-ε ], or > ε OR, or < i= Vs/R/ε [-e -R/L ], or < ε Vs/R/ε [e R/Lε e -R/L ], or > ε 3

aylor Series Approximaion or e x ax ax e! or x, e ax ax! 3 ax ax ax...!! 3! we can drop he higher order erms : ax! ax 33

Narrow Pulse Approximaion Coninued Applying he approximaion or e x, Vs= he uni impulse uncion and α = R/L, or < i= Vs/R/ε [-e -R/L ]= Vs/R/εR/L=Vs/Rα/ε, or < ε Vs/R/ε [e +R/Lε e -R/L ]= Vs/R/εαεe -α, or > ε, or < i= /Rα/ε=α/R, or <ε /R αe -α, or > ε h lim i e R u 34

Diereniaing he Uni Sep uncion he response due o a Uni Sep uncion is i = Vs/R -e -α u and since du, hen d h di d h [ e u e ] R [ e e u ] R e R u 35

Convoluion or Discree Sysems For an LI sysem, le s deine h[n] as he sysem response o a uni impulse source, δ[n]. δ[n]=, n= and or n We have: x[n]= Σ x[m]δ[n-m] y[n]=σ x[m] h[n-m] In addiion he same convoluion properies hold: Commuaive Associaive Disribuive [ n] [ n] [ n] [ n] [ n] { [ n] 3[ n]} { [ n] [ n]} 3[ n] [ n] { [ n] 3[ n]} { [ n] [ n]} { [ n] 3[ ]} 36

Sabiliy o Sysems Coninuous Sysems I a sysem is sable, i.e., Bounded Inpu, Bounded Oupu BIBO, he ollowing mus be rue: x, hen y y h x d h x d x max max OR However, his is no always he case. Posiive Feedbac causes insabiliy h d h d 37

Wha is needed or BIBO For a coninuous ime sysem, he poles o Hp mus lie wihin he le hand complex plane such ha Re s i < where s i are he poles o Hp. his will assure ha he ree response will be damped and no grow exponenially. HIS IS WHY WE SUDIED SOLUIONS OF LINEAR ODE IN ERMS OF SOURCE-FREE AND SOURCE COMPONENS. HIS IMPLIES HA Hp AND HE IMPULSE RESPONSE, h, MAY BE RELAED. 38

Sabiliy o Sysems Discree Sysems I a sysem is sable, i.e., Bounded Inpu, Bounded Oupu BIBO, he ollowing mus be rue: x[n], hen y[n] y[n] h[m]x[n m] x max h[m] OR h[m] However, his is no always he case. Posiive Feedbac causes insabiliy 39

Wha is needed or BIBO For a discree ime sysem, he eigenvalues o hn [ ] mus lie wihin he uni circle n such ha z where z are he eigenvalues o h[ n] Az. i i i i i his will assure ha ree response will no diverge and hn [ ] Using he ormula or he parial sums o a geomeric series, where N is he number o roos o he Characerisic Equaion L L N N L N L N n n i lim hn [ ] lim zi lim zi L L L L n n i i n i zi i z i z, lim lim i L N N L n i zi L L n i i zi z z angle z i i i L L z z Langle z i i i L n n i A z i n i <. z lim ; provided zi z L lim z lim{ z Langle z } lim{ z } { Langle z } i i i i L L L L lim{ z } ; z andlim{ z } ; z L i i i i L L L i i 4

Homewor Convoluion Veriy your all your resuls o hese convoluion problems using Malab and is conv uncion. Problem Assume ha a sysem response is given by he ollowing: h Sech he response o a u, b u-u-a or a=.5, a=, and a=5, and c evaluae e - u a = and = Problem Assume ha a sysem response is given by he ollowing: h /3 3 Evaluae he response o e - u a = and =3

Homewor Sabiliy Deermine he sabiliy o he ollowing sysems wih poles in he complex plane, describe he orm o he ransien response: Imaginary axis Imaginary axis Imaginary axis Real axis Real axis Real axis C.3., C.3. 4