Solution of Linear Constant-Coefficient Difference Equations

Similar documents
ECE 308 Discrete-Time Signals and Systems

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

Discrete-Time Signals and Systems. Discrete-Time Signals and Systems. Signal Symmetry. Elementary Discrete-Time Signals.

DIGITAL SIGNAL PROCESSING LECTURE 3

ECE4270 Fundamentals of DSP. Lecture 2 Discrete-Time Signals and Systems & Difference Equations. Overview of Lecture 2. More Discrete-Time Systems

Z - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.

Signals and Systems. Problem Set: From Continuous-Time to Discrete-Time

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science. BACKGROUND EXAM September 30, 2004.

GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL OF ELECTRICAL AND COMPUTER ENGINEERING

ADVANCED DIGITAL SIGNAL PROCESSING

Classification of DT signals

Time-Domain Representations of LTI Systems

ELEG 4603/5173L Digital Signal Processing Ch. 1 Discrete-Time Signals and Systems

COMM 602: Digital Signal Processing

FIR Filters. Lecture #7 Chapter 5. BME 310 Biomedical Computing - J.Schesser

Discrete-Time Signals and Systems. Signals and Systems. Digital Signals. Discrete-Time Signals. Operations on Sequences: Basic Operations

Recurrence Relations

Recursive Algorithms. Recurrences. Recursive Algorithms Analysis

MAS160: Signals, Systems & Information for Media Technology. Problem Set 5. DUE: November 3, (a) Plot of u[n] (b) Plot of x[n]=(0.

Ch3 Discrete Time Fourier Transform

EECE 301 Signals & Systems

Computing the output response of LTI Systems.

LINEAR RECURSION RELATIONS - LESSON FOUR SECOND-ORDER LINEAR RECURSION RELATIONS

You may work in pairs or purely individually for this assignment.

Question1 Multiple choices (circle the most appropriate one):

2.4 - Sequences and Series

Discrete-Time Systems, LTI Systems, and Discrete-Time Convolution

Chapter 2 Systems and Signals

Lecture 3. Digital Signal Processing. Chapter 3. z-transforms. Mikael Swartling Nedelko Grbic Bengt Mandersson. rev. 2016

Analog and Digital Signals. Introduction to Digital Signal Processing. Discrete-time Sinusoids. Analog and Digital Signals

SOLUTION SET VI FOR FALL [(n + 2)(n + 1)a n+2 a n 1 ]x n = 0,

CSI 2101 Discrete Structures Winter Homework Assignment #4 (100 points, weight 5%) Due: Thursday, April 5, at 1:00pm (in lecture)

Signal Processing. Lecture 02: Discrete Time Signals and Systems. Ahmet Taha Koru, Ph. D. Yildiz Technical University.

Digital signal processing: Lecture 5. z-transformation - I. Produced by Qiangfu Zhao (Since 1995), All rights reserved

Engineering Analysis ( & ) Lec(7) CH 2 Higher Order Linear ODEs

Filter banks. Separately, the lowpass and highpass filters are not invertible. removes the highest frequency 1/ 2and

Appendix: The Laplace Transform

Olli Simula T / Chapter 1 3. Olli Simula T / Chapter 1 5

Chapter 4 : Laplace Transform

CHAPTER 1 SEQUENCES AND INFINITE SERIES

Signals & Systems Chapter3

Math 475, Problem Set #12: Answers

6.003 Homework #3 Solutions

Dynamic Response of Linear Systems

SOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.

Problem. Consider the sequence a j for j N defined by the recurrence a j+1 = 2a j + j for j > 0

Frequency Response of FIR Filters

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

Introduction to Signals and Systems, Part V: Lecture Summary

Review Problems 1. ICME and MS&E Refresher Course September 19, 2011 B = C = AB = A = A 2 = A 3... C 2 = C 3 = =

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

Unit 4: Polynomial and Rational Functions

EE Midterm Test 1 - Solutions

Find quadratic function which pass through the following points (0,1),(1,1),(2, 3)... 11

Chapter 7: The z-transform. Chih-Wei Liu

Unit 6: Sequences and Series

Solutions to Homework 1

MAXIMALLY FLAT FIR FILTERS

Solutions. Number of Problems: 4. None. Use only the prepared sheets for your solutions. Additional paper is available from the supervisors.

Section 11.8: Power Series

Zeros of Polynomials

Curve Sketching Handout #5 Topic Interpretation Rational Functions

Lecture 2 Linear and Time Invariant Systems

PAPER : IIT-JAM 2010

Generating Functions. 1 Operations on generating functions

: Transforms and Partial Differential Equations

CSE 1400 Applied Discrete Mathematics Number Theory and Proofs

Linear time invariant systems

Section A assesses the Units Numerical Analysis 1 and 2 Section B assesses the Unit Mathematics for Applied Mathematics

ELEG3503 Introduction to Digital Signal Processing

Generalizing the DTFT. The z Transform. Complex Exponential Excitation. The Transfer Function. Systems Described by Difference Equations

Chapter 3. z-transform

A 2nTH ORDER LINEAR DIFFERENCE EQUATION

2D DSP Basics: 2D Systems

Find a formula for the exponential function whose graph is given , 1 2,16 1, 6

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

Math 312 Lecture Notes One Dimensional Maps

[ 11 ] z of degree 2 as both degree 2 each. The degree of a polynomial in n variables is the maximum of the degrees of its terms.

PROBLEMS AND SOLUTIONS 2

Created by T. Madas SERIES. Created by T. Madas

Linear regression. Daniel Hsu (COMS 4771) (y i x T i β)2 2πσ. 2 2σ 2. 1 n. (x T i β y i ) 2. 1 ˆβ arg min. β R n d

It is often useful to approximate complicated functions using simpler ones. We consider the task of approximating a function by a polynomial.

Linear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy

k-generalized FIBONACCI NUMBERS CLOSE TO THE FORM 2 a + 3 b + 5 c 1. Introduction

Solutions of Chapter 5 Part 1/2

10.1 Sequences. n term. We will deal a. a n or a n n. ( 1) n ( 1) n 1 2 ( 1) a =, 0 0,,,,, ln n. n an 2. n term.

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Chapter 4. Fourier Series

Mathematics review for CSCI 303 Spring Department of Computer Science College of William & Mary Robert Michael Lewis

Solution of EECS 315 Final Examination F09

Transfer Function Analysis

The z-transform can be used to obtain compact transform-domain representations of signals and systems. It

Mathematics 116 HWK 21 Solutions 8.2 p580

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Polynomials with Rational Roots that Differ by a Non-zero Constant. Generalities

Module 18 Discrete Time Signals and Z-Transforms Objective: Introduction : Description: Discrete Time Signal representation

Solutions - Homework # 1

e to approximate (using 4

Transcription:

ECE 38-9 Solutio of Liear Costat-Coefficiet Differece Equatios Z. Aliyazicioglu Electrical ad Computer Egieerig Departmet Cal Poly Pomoa Solutio of Liear Costat-Coefficiet Differece Equatios Example: Determie the respose y ( ), of the system described by the secod-order differece equatio y ( ).7 y ( ). y ( ) x ( ) x ( ) to the iput x( ) 4 u( ) The homogeous solutio is y ( ).7 y ( ). y ( ) λ λ λ.7. ( ) λ λ λ.7. λ ad.5 λ. y ( ) c.5 c. h ECE 38-9

Solutio of Liear Costat-Coefficiet Differece Equatios Particular Solutio: y ( ) K4 u( ) p K u K u K u u u 4 ( ).7 4 ( ). 4 ( ) ()4 ( ) 4 ( ) K4.7K4.K4 (4) 4 6K.8K.K 3 y ( ).33(4) u( ) p The total solutio 3 K.33 3.3 y ( ) c.5 c..33(4) u ( ) ECE 38-9 3 Solutio of Liear Costat-Coefficiet Differece Equatios To fid c ad c For : From differece equatio, y().7 y( ). y( ) x() x( ) y() From the total solutio, y() c c.33, For : From differece equatio, y().7 y( ). y( ) x() x( ) y().4 8 9.4 From differece equatio, y().5c.c 9.3 Therefore, c c.33 9.4.5c.c 9.3 c.87 c.466 Total Solutio y ( ).466(.5).87(.).33(4) u ( ) ECE 38-9 4

The Impulse Respose of a LTI recursive system I geeral case If the iput y ( ) h( ) x( ), zs x ( ) δ ( ), the we obtai y ( ) h( ) zs The impulse respose ca be obtaied from the liear costatcoefficiet differece equatio. That is the solutio of homogeeous equatio ad particular solutio to the excitatio fuctio. I the case where the excitatio fuctio is a impulse fuctio. The particular solutio is zero y ( p ), sice x( ) for >. ECE 38-9 5 The Impulse Respose of a LTI recursive system Example: Fid the impulse respose form the followig equatio, y ( ) 3 y ( ) 4 y ( ) x ( ) x ( ) The homogeous solutio is y ( ) C ( ) C 4 h The particular solutio is zero whe To fid C ad C, we evaluate differece equatio ad homogeous solutio for ad.( y( ), y( ), sice the system must be relaxed) y() y() 3 y() 5 The impulse respose is y() C C y() C 4C x( ) δ ( ) C 6 5 6 h ( ) ( ) 4 u ( ) 5 5 C 5 ECE 38-9 6

Implemetatio of Discrete-Time Systems A system ca be described by a liear costat-coefficiet differece equatio. Let s cosider the first order system y ( ) ay ( ) bx ( ) bx ( ) The system ca be described by two systems i cascade. The first is a orecursive system described by the equatio v ( ) bx ( ) bx ( ) The secod part is recursive system x() y ( ) ay ( ) v ( ) b z - z -a - b v() y() This is called a direct form I structure. ECE 38-9 7 Implemetatio of Discrete-Time Systems Usig covolutio properties, we ca iterchage the order of the recursive ad orecursive system x() w() b y() z - z - -a b w(-) w(-) w ( ) aw ( ) x ( ) y ( ) bw ( ) bw ( ) x() w() z - b -a b w(-) y() Usig oly oe delay. It is more efficiet i terms of memory requiremet. It is called the direct form II structure ECE 38-9 8

Implemetatio of Discrete-Time Systems I geeral form, The differece equatio is give by y( ) a y( ) b x( ) The orecursive system is The recursive system is v ( ) bx ( ) y ( ) ay ( ) v ( ) Direct form II structure, the recursive system is w ( ) aw ( ) x ( ) The orecursive system is y ( ) bw ( ) ECE 38-9 9 Implemetatio of Discrete-Time Systems x() b v() z - z - b -a z - z - -a b y() x() w() z - b -a b z - -a b y() z - z - -a b M z - -a - b M z - -a w(-) Step- Step-3 ECE 38-9

Implemetatio of Discrete-Time Systems A special case of geeral case a,,,..., y ( ) bx ( ) which is a orecursive LTI system. Such a system coefficiets b b M h ( ) otherwise The secod part, set M to obtai the geeral case differece equatio y ( ) ay ( ) bx ( ) ECE 38-9 Crosscorrelatio ad Autocorrelatio Sequece The crosscorelatio of x() ad y() is s sequece defied as Or r ( l) x( ) y( l), l, ±, ±,... xy r () l x( l) y( ), l, ±, ±,... xy The reverse crosscorrelatio is Or r ( l) y( ) x( l), l, ±, ±,... yx r ( l) y( l) x( ), l, ±, ±,... yx rxy () l is r () l r ( l) xy y x ECE 38-9

Crosscorrelatio ad Autocorrelatio Sequece I special case, we have autocorrelatio, which is defied as So, that r ( ) ( ) ( ) xx l xx l rxx () l x ( ) x ( ) Or Example: x ( ) {...,,,,3,7,,, 3,,,...} y ( ) {...,,,,,,4,,,5,,,...} xy xy r () x( ) y( ) 6 4 4 6 7 r () x( ) y( ) 3 4 8 3 3 rxy ( l ) {, 9,9,36, 4,33,,7,3, 8,6, 7,5, 3} ECE 38-9 3 Problem 43.a Problem Solutios Determie the direct form II realizatio for the followig LTI system y ( ) y ( ) 4 y ( 3) x ( ) 3 x ( 5) 3 y ( ) y ( ) y ( 3) x ( ) x ( 5) x() / z - y() -/ z - z - z - z - 3/ ECE 38-9 4

Problem 44.a Problem Solutios Compute the first samples of its impulse respose x() y() / z - x( ) δ ( ) x ( ) {,,,... } y ( ) y ( ) x ( ) x ( ) y() x() 3 y() y() x() x() 3 y() y() x() x() 4 3 y(3) y() x(3) x() 8 3 3 3 3 3 3 y ( ),,,,,,,... 4 8 `6 3 64 ECE 38-9 5 Problem Solutios Problem 44.b Fid the iput output relatio y ( ) y ( ) x ( ) x ( ) Problem 44.c x ( ) {,,,... } The iput y ( ) y ( ) x ( ) x ( ) y() x() 5 y() y() x() x() 3 y() y() x() x() 4 9 y(3) y() x(3) x() 8 ECE 38-9 6

Problem 44.d Problem Solutios Use covolutio y ( ) u ( * h ( ) uh ( ) ( ) h ( ) y() h() 5 y() h() h() 3 y() h() h() h() 4 9 y(3) h() h() h() h(3) 8 ECE 38-9 7 Problem 54 Problem Solutios Fid the y() for the followig equatio y ( ) 4 y ( ) 4 y ( ) x ( ) x ( ) x( ) ( ) u( ) The characteristic equatio λ 4λ 4 λ, λ The homogeous solutio is y ( ) c c The particular solutio is h y ( ) ( ) u( ) p u u u u u ( ) ( ) 4 ( ) ( ) 4 ( ) ( ) ( ) ( ) ( ) ( ) For 4 4 9 ECE 38-9 8

Problem Solutios The total solutio is y ( ) yh( ) yp( ) c c ( ) ( ) 9 u Usig the iitial coditio, y( ) y( ) we ca obtai from differece equatio at y() 4 y( ) 4 y( ) x() x( ) y() y() 4 y() 4 y( ) x() x() y() From the total solutio y() c y() c c 9 9 7 9 c c 3 The total solutio y() c c 9 7 y ( ) ( ) u ( ) 9 3 9 ECE 38-9 9 Problem Solutios Problem 55 Fid the impulse respose h() for the followig equatio y ( ) 4 y ( ) 4 y ( ) x ( ) x ( ) The homogeous solutio is So system respose To fid the costat y ( ) c c h y() 4 y( ) 4 y( ) δ () δ ( ) y() h ( ) c c u ( ) y() 4 y() 4 y( ) δ () δ () y() 3 From the system respose h() y() c y() c c 3 c h ( ) u ( ) ECE 38-9