Chapter 2 z-transform. X(z) = Z[x(t)] = Z[x(kT)] = Z[x(k)] x(kt)z k = x(kt)z k = x(k)z k. X(z)z k 1 dz 2πj c = 1

Similar documents
It is common to think and write in time domain. creating the mathematical description of the. Continuous systems- using Laplace or s-

EE480.3 Digital Control Systems. Part 2. z-transform

Automatique. A. Hably 1. Commande d un robot mobile. Automatique. A.Hably. Digital implementation

Lecture 04: Discrete Frequency Domain Analysis (z-transform)

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

Chapter 7: The z-transform

1. Z-transform: Initial value theorem for causal signal. = u(0) + u(1)z 1 + u(2)z 2 +

Discrete-Time Signals and Systems. The z-transform and Its Application. The Direct z-transform. Region of Convergence. Reference: Sections

Signals and Systems Lecture 8: Z Transform

SIGNALS AND SYSTEMS. Unit IV. Analysis of DT signals

Digital Signal Processing

5 Z-Transform and difference equations

Introduction to the z-transform

The Convolution Sum for Discrete-Time LTI Systems

X (z) = n= 1. Ã! X (z) x [n] X (z) = Z fx [n]g x [n] = Z 1 fx (z)g. r n x [n] ª e jnω

Signals & Systems Handout #4

Chapter 13 Z Transform

y[n] = = h[k]x[n k] h[k]z n k k= 0 h[k]z k ) = H(z)z n h[k]z h (7.1)

ESE 531: Digital Signal Processing

DIGITAL SIGNAL PROCESSING. Chapter 3 z-transform

Methods for analysis and control of. Lecture 6: Introduction to digital control

III. Time Domain Analysis of systems

ECE503: Digital Signal Processing Lecture 4

Z-Transform. x (n) Sampler

Sampling and Discrete Time. Discrete-Time Signal Description. Sinusoids. Sampling and Discrete Time. Sinusoids An Aperiodic Sinusoid.

DSP-I DSP-I DSP-I DSP-I

3.2 Complex Sinusoids and Frequency Response of LTI Systems

Z-Transform. The Z-transform is the Discrete-Time counterpart of the Laplace Transform. Laplace : G(s) = g(t)e st dt. Z : G(z) =

ESE 531: Digital Signal Processing

Use: Analysis of systems, simple convolution, shorthand for e jw, stability. Motivation easier to write. Or X(z) = Z {x(n)}

EL 625 Lecture 11. Frequency-domain analysis of discrete-time systems

Discrete Time Signals and Systems Time-frequency Analysis. Gloria Menegaz

Section 3 Discrete-Time Signals EO 2402 Summer /05/2013 EO2402.SuFY13/MPF Section 3 1

Identification Methods for Structural Systems

Signal Analysis, Systems, Transforms

EE 438 Essential Definitions and Relations

1 Lecture 8: Interpolating polynomials.

Interchange of Filtering and Downsampling/Upsampling

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014

Models for Sampled Data Systems

EE Homework 5 - Solutions

UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH

ELEG 305: Digital Signal Processing

EE 224 Signals and Systems I Review 1/10

6.003: Signals and Systems

CS-9645 Introduction to Computer Vision Techniques Winter 2018

Elementary Differential Equations, Section 2 Prof. Loftin: Practice Test Problems for Test Find the radius of convergence of the power series

Discrete-time linear systems

Conformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.

Chapter 12 - Solved Problems

NAME: 13 February 2013 EE301 Signals and Systems Exam 1 Cover Sheet

PUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES. Notes

Linear System Theory. Wonhee Kim Lecture 1. March 7, 2018

Z Transform (Part - II)

DESIGN OF CMOS ANALOG INTEGRATED CIRCUITS

The Z-Transform. Fall 2012, EE123 Digital Signal Processing. Eigen Functions of LTI System. Eigen Functions of LTI System

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

ADVANCED ENGINEERING MATHEMATICS MATLAB

ELEG 305: Digital Signal Processing

Digital Filters Ying Sun

LTI Systems (Continuous & Discrete) - Basics

ELEC2400 Signals & Systems

Like bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform.

Need for transformation?

GATE EE Topic wise Questions SIGNALS & SYSTEMS

A.1 THE SAMPLED TIME DOMAIN AND THE Z TRANSFORM. 0 δ(t)dt = 1, (A.1) δ(t)dt =

UNIT-II Z-TRANSFORM. This expression is also called a one sided z-transform. This non causal sequence produces positive powers of z in X (z).

2. CONVOLUTION. Convolution sum. Response of d.t. LTI systems at a certain input signal

Communication Signals (Haykin Sec. 2.4 and Ziemer Sec Sec. 2.4) KECE321 Communication Systems I

Review of Linear Time-Invariant Network Analysis

Lecture 2: The z-transform

SEL4223 Digital Signal Processing. Inverse Z-Transform. Musa Mohd Mokji

The z-transform Part 2

The Generalized Laplace Transform: Applications to Adaptive Control*

1.1 SPECIAL FUNCTIONS USED IN SIGNAL PROCESSING. δ(t) = for t = 0, = 0 for t 0. δ(t)dt = 1. (1.1)

Automatic Control Systems theory overview (discrete time systems)

ECE-S Introduction to Digital Signal Processing Lecture 4 Part A The Z-Transform and LTI Systems

Module 4 : Laplace and Z Transform Problem Set 4

Definition of the Laplace transform. 0 x(t)e st dt

QUESTION BANK SIGNALS AND SYSTEMS (4 th SEM ECE)

# FIR. [ ] = b k. # [ ]x[ n " k] [ ] = h k. x[ n] = Ae j" e j# ˆ n Complex exponential input. [ ]Ae j" e j ˆ. ˆ )Ae j# e j ˆ. y n. y n.

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

Introduction to Signals and Systems Lecture #4 - Input-output Representation of LTI Systems Guillaume Drion Academic year

Review of Discrete-Time System

8 Periodic Linear Di erential Equations - Floquet Theory

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

The Laplace Transform

Table 1: Properties of the Continuous-Time Fourier Series. Property Periodic Signal Fourier Series Coefficients

Digital Signal Processing. Lecture Notes and Exam Questions DRAFT

Fourier series for continuous and discrete time signals

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

Discrete Systems. Step response and pole locations. Mark Cannon. Hilary Term Lecture

EEM 409. Random Signals. Problem Set-2: (Power Spectral Density, LTI Systems with Random Inputs) Problem 1: Problem 2:

Chapter 7. Digital Control Systems

Linear Systems Theory

ECE 301 Fall 2011 Division 1. Homework 1 Solutions.

Generalizing the DTFT!

! z-transform. " Tie up loose ends. " Regions of convergence properties. ! Inverse z-transform. " Inspection. " Partial fraction

Review of Frequency Domain Fourier Series: Continuous periodic frequency components

EE 225D LECTURE ON DIGITAL FILTERS. University of California Berkeley

Transcription:

One-sided -transform Two-sided -transform Chapter 2 -TRANSFORM X() Z[x(t)] Z[x(kT)] Z[x(k)] x(kt) k x(k) k X() k0 k x(kt) k k0 k x(k) k Note that X() x(0) + x(t ) + x(2t ) 2 + + x(kt) k + Inverse -transform Z [X()] x(kt) x(k) X() k d 2πj c Where c is a circle with its center at the origin of the plane such that all poles of X() k are inside it Z Transform of elementary functions: Unit step function { (t) 0 t x(t) 0 t<0 X() Z[(t)] k k k0 k0 Region of convergence > Geometric series a + ar + ar 2 + ar 3 + a r r < Exponential Function x(t) { e at 0 t 0 t<0 x(kt) e akt, k 0,, 2, X() Z[e at ] x(kt) k e akt k k0 k0 e at e at

See table of -transforms on page 29 and 30 (new edition), or page 49 and 50 (old edition). The -transform X() and its inverse x(k) have a one-to-one correspondence, however, the -transform X() and its inverse -transform x(t) do not have a unique correspondence. Properties and theorems of the -transform Multiplication by a constant: Z[ax(t)] ax() Linearity: Z[αf(k) +βg(k)] αf () +βg() Multiplication by a k : Z[a k x(k)] X(a ) Real translation theorem (shifting theorem): If x(t) 0for t<0 and Z[x(t nt )] n X() Z[x(t + nt )] n [X() n k0 x(kt) k ] Initial value theorem: x(0) lim X() Final value theorem: lim x(k) lim [( k )X()] Real convolution Theorem: let x (t) 0for t<0 x 2 (t) 0for t<0 then k X () X 2 () Z[ x (ht ) x 2 (kt ht )] h0 2

INVERSE TRANSFORM Different Methods. Direct division method (Power Series Method) 2. Computational method 3. Partial-fraction-expansion method 4. Inversion integral method Direct division method Express X() inpowersof Example Find Z of X() +2 +3 2 +4 3 Solution: x(0) ; x() 2; x(2) 3; x(3) 4 Example 2 Solution: Find Z of X() 0 +5 ( )( 0.2) X() 0 +5 2.2 +0.2 2 X() 0 +7 2 +8.4 3 +8.68 4 + x(0) 0 x() 0 x(2) 7 x(3) 8.4 x(4) 8.68 Computational method X() 0 +5 ( )( 0.2) Solution: 0 +5 Let X() 2.2 +0.2 U() where U() now, U() u(0) + u() + u(2) 2 + + u(k) k + 3

for U() u(0) u(k) 0, for k, 2, 3 Converting to difference equation now, x(k +2).2 x(k +)+0.2 x(k) 0u(k +)+5u(k) let k 2 ( ) x(0).2 x( ) + 0.2 x( 2) 0 u( ) + 5 u( 2) now, x( ) x( 2) 0 and u( ) u( 2) 0 x(0) 0 Similarly, we find x() 0 We may continue the process to find x(k), k 2, 3, using ( ) Partial Fraction Expansion To find the Z X(), we may expand X() or X() into partial fractions. X() is expanded since each of the expanded terms is generally available in -transform tables. Alternatively, X() may be expanded and use of the shifting theorem may be made. Example let Y () X() a X() a now, thus, Z {Y ()} y(k) a k X() Y () Z {X()} x(k) y(k ) a k { a k k, 2, 3, x(k) 0 k 0 4

General procedure for partial fraction expansion: let Given X(), find X() X() If M>N, no adjustment need be made to X(), If N>M, we divide through X() a 0 + a + + a N N b 0 + b + + b M M () c N M N M + c N M N M + + c + c 0 + d 0 + d + + d M M } b 0 + b + + b M {{ M } ψ() Factoring ψ() where we have one repeated pole of order k, callit r,and the rest unique, k+, k+2,, m ψ() A k ( r ) k + A k ( r ) k + + A r + M jk+ A j j (2) Where A j (k j)! [ dk j d k j ( j ) k ψ()] j, j, 2,,k A j ( j ) ψ() j, j k +,k+2,,m Substituting (3) into (2) and multiplying by and taking inverse transform gives us: Z [X()] x(n) Z [ c N M N M+ + c N M N M + + c 2 + c 0 ] k + Z A j M + Z A j ( r ) j r j jk+ x(n) N nm C N M δ(n +(N M +)) + [A N r + A 2 n n r + M jk+ A j n j ] u(n) + + A k n(n ) (n (k 2)) r n k+ (k )! Where the following has been used { } Z n (n ) (n (k 2)) an k+ u(n) ( a) k (k )! where u(n) is the unit step function. 5

Example Find Z {X()} where, X() 4 + 2 ( 2 )( 4 ) Solution now, Where Thus X() 3 + 2 3 + 4 8 23 3 6 32 2 3 + 4 8 A A 2 X() 23 6 3 32 4 23 3 6 32 2 + 3 4 + 23 6 3 32 2 3 4 + 8 A 2 4 2 + 3 4 + 5 2 2 + A 2 4 5 8 4 7 64 4 5 2 7 6 7 6 4 x(n) Z [ 2 + 3 4 ] + Z [ 5 2 2 7 ] 6 4 δ(n +2)+ 3 4 δ(n +)+ [ 5 2 ( 2 )n 7 6 ( 4 )n ] u(n) 6

Inversion integral method Background material: Suppose 0 is an isolated singular point (pole) of F(). Expand F() in a Laurent series about 0 F () a n ( 0 ) n b n + n0 n ( 0 ) n where a n 2πj b n 2πj Γ Γ 2 F () d n 0,, 2, ( 0 ) n+ F () d n, 2, 3, ( 0 ) n+ where Γ and Γ 2 are closed paths around 0 and b F () d() 2πj Γ where Γ is any closed path within and on which F() is analytic except at 0,andb is called the residue of F() at the pole 0. Now F () d Γ F () d + Γ F () d + + Γ 2 F () d Γ m 2πj(b + b 2 + + b m ) Residue theorem 7

Inversion integral X() x(kt) k x(0) + x(t ) + x(2t ) 2 + + x(kt) k + k0 X() k x(0) k + x(t ) k 2 + x(2t ) k 3 + + x(kt) + Note, this is the Laurent series expression of X() k around point 0, and x(kt) is the residue x(kt) 2πj c X()k d the inverse integral for the -transform Inverse transform using inversion integral M x(k) x(kt) [residue of X() k at pole i of X() k ] i assuming M poles. The residue K, for simple pole is given by K lim i [( i )X() k ] The residue K, for multiple pole j of order q is given by K (q )! lim j d q d q [( j) q X() k ] 8

Example Find Z [X()], where X() 2 ( ) 2 ( e at ) Solution: X() k k+ ( ) 2 ( e at ) Simple pole at e at Double pole at [ 2 k+ ] x(k) residue of i ( ) 2 ( e at ) at pole i K + K 2 where [ K lim ( e at k+ ] ) e a(k+)t e at ( ) 2 ( e at ) ( e at ) 2 K 2 x(kt) (2 )! lim d d [ ( ) 2 k+ ( ) 2 ( e at ) k e e at see steps below at ( e at ) 2 kt T ( e at ) ( e akt ) e at k 0,, 2, ( e at ) 2 ] Steps d v u udv vdu u 2 lim d ( k+ ) d e at (k +) k ( e at ) k+ lim ( e at ) 2 [ (k +) k lim e k+ ] at ( e at ) 2 k + e at ( e at ) 2 k e at + e at ( e at ) 2 ( e at ) 2 k e e at at ( e at ) 2 9

Pulse-Transfer Function Difference equation: x(k)+a x(k ) + + a n x(k n) b 0 u(k)+b u(k ) + + b n u(k n) Taking transform X()+a X()+ + a n n X() b 0 U()+b U()+ + b n n U() G() X() U() b 0 + b + + b n n +a + + a n n Pulse Transfer Function Now, Kronecker delta function δ 0 (kt) δ 0 (kt) { for k 0 0 for k 0 Z [δ 0 (kt)] G() is the transform of the response to δ 0 (kt). It is called the pulse transfer function g(k) Z {G()} is called the weighting sequence. 0

transform method of solving difference equations Example Solve: x(k +2)+3x(k +)+2x(k) 0; x(0) 0, x() Solution taking the transform 2 X() 2 x(0) x() + 3 X() 3 x(0) + 2X() 0 Substituting initial data X() 2 +3 + ( +)( +2) + + +2 +2 [ ] [ ] Z ( ) k, Z ( 2) k + +2 x(k) ( ) k ( 2) k, k 0,, 2,