<< Chapter < Page Chapter >> Page >

We begin with the simplest of discrete-time systems, where X = C N and Y = C M . In this case a linear operator is just an M × N matrix. We can generalize this concept by letting M and N go to , in which case we can think of a linear operator L : 2 ( Z ) 2 ( Z ) as an infinite matrix.

Consider the shift operator Δ k : 2 ( Z ) 2 ( Z ) that takes a sequence and shifts it by k . As an example, Δ 1 can be viewed as the infinite matrix given by

y - 1 y 0 y 1 = 0 0 1 0 0 1 0 0 1 0 0 x - 1 x 0 x 1

Note that Δ k 2 = 1 (for any k and p ) since the delay doesn't change the norm of x . The delay operator is also an example of a linear shift-invariant (LSI) system.

Definition 1

An operator L : 2 ( Z ) 2 ( Z ) is called shift-invariant if L ( Δ k ( x ) ) = Δ k ( L ( x ) ) for all x 2 ( Z ) and for any k Z .

Observe that Δ k 1 ( Δ k 2 ( x ) ) = Δ k 1 + k 2 ( x ) so that Δ k itself is an LSI operator.

Lets take a closer look at the structure of an LSI system by viewing it as an infinite matrix. In this case we write y = H x to denote

y - 1 y 0 y 1 = | | | h - 1 h 0 h 1 | | | x - 1 x 0 x 1

Suppose we want to figure out the column of H corresponding to h 0 . What input x could help us determine h 0 ? Consider the vector

x = 0 1 0 ,

i.e., x = δ [ n ] . For this input y = H x = h 0 . What about h 1 ? Δ 1 ( x ) = δ [ n - 1 ] would yield h 1 . In general Δ k ( x ) = δ [ n - k ] tell us the column h k . But , if H is LSI, then

h k = H ( Δ k ( δ [ n ] ) ) = Δ k ( H ( δ [ n ] ) ) = Δ k ( h 0 )

This means that each column is just a shifted version of h 0 , which is usually called the impulse response.

Now just to keep notation clean, let h = h 0 denote the impulse response. Can we get a simple formula for the output y in terms of h and x ? Observe that we can write

y - 1 y 0 y 1 = h 0 h - 1 h - 2 h 1 h 0 h - 1 h 2 h 1 h 0 x - 1 x 0 x 1

Each column is just shifted down one. (Each successive row is also shifted right one.) Looking at y - 1 , y 0 and y 1 , we can rewrite this formula as

y [ - 1 ] y [ 0 ] y [ 1 ] = + x [ - 1 ] h [ 0 ] h [ 1 ] h [ 2 ] + x [ 0 ] h [ - 1 ] h [ 0 ] h [ 1 ] + x [ 1 ] h [ - 2 ] h [ - 1 ] h [ 0 ] +

From this we can observe the general pattern

y [ n ] = + x [ - 1 ] h [ n + 1 ] + x [ 0 ] h [ n + 0 ] + x [ 1 ] h [ n - 1 ] +

or more concisely

y [ n ] = k = - x [ k ] h [ n - k ] .

Does this look familiar? It is simply the formula for the discrete-time convolution of x and h , i.e.,

y = x * h .

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Digital signal processing. OpenStax CNX. Dec 16, 2011 Download for free at http://cnx.org/content/col11172/1.4
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Digital signal processing' conversation and receive update notifications?

Ask