22 The z-Transform

Recommended
Problems
P22.1
An LTI system has an impulse response h[n] for which the z-transform is

H(z) = h[n]z- z-1 , |

(a) Plot the pole-zero pattern for H(z).
(b) Using the fact that signals of the form z" are eigenfunctions of LTI systems,
determine the system output for all n if the input x[n] is
x[n] = (?)" + 3(2)n

P22.2
Consider the sequence x[n] = 2"u[n].
(a) Is x[n] absolutely summable?
(b) Does the Fourier transform of x[n] converge?
(c) For what range of values of r does the Fourier transform of the sequence
r-"x[n] converge?
(d) Determine the z-transform X(z) of x[n], including a specification of the ROC.
(e) X(z) for z = 3e'a can be thought of as the Fourier transform of a sequence x,[n],
i.e.,
2"u[n] X(z),
xi[n] 5 X(3ei") = XI(ei")
Determine x1 [n].

P22.3
Shown in Figure P22.3 is the pole-zero plot for the z-transform X(z) of a sequence
x[nI.

Im
z plane

X X XRe
1 2 2
3 3

Figure P22.3
P22-1
Signals and Systems
P22-2

Determine what can be inferred about the associated region of convergence from
each of the following statements.
(a) x[n] is right-sided.
(b) The Fourier transform of x[n] converges.
(c) The Fourier transform of x[n] does not converge.
(d) x[n] is left-sided.

P22.4
(a) Determine the z-transforms of the following two signals. Note that the z-trans-
forms for both have the same algebraic expression and differ only in the ROC.
(i) xi[n] = (i)"u[n]
(ii) X2[n] = -(I)"n[-U - 1]
(b) Sketch the pole-zero plot and ROC for each signal in part (a).
(c) Repeat parts (a) and (b) for the following two signals:
(i) x 3[n] = 2u[n]
(ii) x 4 [n] = -(2)"u[-n - 1]
(d) For which of the four signals xi[n], x 2[n], xs[n], and x 4[n] in parts (a) and (c)
does the Fourier transform converge?

P22.5

Consider the pole-zero plot of H(z) given in Figure P22.5, where H(a/2) = 1.

z plane

Figure P22.5
(a) Sketch |H(e'")I as the number of zeros at z = 0 increases from 1 to 5.
(b) Does the number of zeros affect <'H(ej")? If so, specifically in what way?
(c) Find the region of the z plane where IH(z) I = 1.
 The z-Transform / Problems
P22-3

P22.6
Determine the z-transform (including the ROC) of the following sequences. Also
sketch the pole-zero plots and indicate the ROC on your sketch.
(a) (I)"U [n]
(b) 6[n + 1]

P22.7
For each of the following z-transforms determine the inverse z-transform.
1 1
(a) X(z) = 1± zIZ Iz| >
1 -z 1z 1
(b) X(z) = 1 z-2 , Iz| >

1 -az- 1
(c) X(z)= z _-' IzI > iaa

Optional
Problems
P22.8
In this problem we study the relation between the z-transform, the Fourier trans-
form, and the ROC.
(a) Consider the signal x[n] = u[n]. For which values of r does r-"x[n] have a con-
verging Fourier transform?
(b) In the lecture, we discussed the relation between X(z) and I{r-"x[nj}. For each
of the following values of r, sketch where in the z plane X(z) equals the Fourier
transform of r-x[n].
(i) r=1
(ii) r=
(iii) r= 3
(c) From your observations in parts (a) and (b), sketch the ROC of the z-transform
of u[n].

P22.9
(a) Suppose X(z) on the circle z = 2e'j is given by
1
X(2ej') = 1 _

Using the relation X(reju) = 5{r~"x[n]}, find 2-x[n] and then x[n], the inverse
z-transform of X(z).
(b) Find x[n] from X(z) below using partial fraction expansion, where x[n] is
known to be causal, i.e., x[n] = 0 for n < 0.
3 + 2z~1
X(z) -2 + 3z-' + z- 2
Signals and Systems
P22-4

P22.10
A discrete-time system with the pole-zero pattern shown in Figure P22.10-1 is
referred to as a first-order all-pass system because the magnitude of the frequency
response is a constant, independent of frequency.

Urn
Unit circle

I ROC: I z I > a

Figure P22.10-1

(a) Demonstrate algebraically that IH(eu) I is constant.

(b) To demonstrate the same property geometrically, consider the vector diagram
in Figure P22.10-2. Show that the length of v 2 is proportional to the length of vi
independent of 0 by following these two steps:
(i) Express the length of v, using the law of cosines and the fact that it is
one leg of a triangle for which the other two legs are the unit vector and
a vector of length a.
(ii) In a manner similar to that in step (i), determine the length of v 2 and show
that it is proportional in length to vi independent of Q.

9JM

Unit circle

Figure P22.10-2

P22.11

Parts (a)-(e) (Figures P22.11-1 to P22.11-5) give pole-zero plots, and parts (i)-(iv)
(Figures P22.11-6 to P22.11-9) give sketches of possible Fourier transform magni-
tudes. Assume that for all the pole-zero plots, the ROC includes the unit circle. For
each pole-zero plot (a)-(e), specify which one if any of the sketches (i)-(iv) could
represent the associated Fourier transform magnitude. More than one pole-zero plot
may be associated with the same sketch.
 The z-Transform / Problems
P22-5

(a) (b)

Im z plane

F PRer

Figure P22.11-1

(c) (d)

(e)

(ii)

7T

Figure P22.11-7
Signals and Systems
P22-6

(iii) (iv)

iTV

ir
Figure P22.11-8 Figure P22.11-9

P22.12
Determine the z-transform for the following sequences. Express all sums in closed
form. Sketch the pole-zero plot and indicate the ROC. Indicate whether the Fourier
transform of the sequence exists.
(a) (i)"{u[n] - u[n - 10]}
(b) (i)II

(c) 7 cos [ + r u[n]

0, n < 0
(d) x[n]= 1, 0 s n - 9
0, 9< n

P22.13

Using the power-series expansion

log(1 - W) = - Wi , w l< 1,

determine the inverse of the following z-transforms.

(a) X(z) = log(1 - 2z), IzI < i
(b) X(z) = log(1 - z- 1), |z > i
MIT OpenCourseWare
http://ocw.mit.edu

Resource: Signals and Systems

Professor Alan V. Oppenheim

The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
23 Mapping Continuous-Time Filters
to Discrete-Time Filters
Recommended
Problems
P23.1
For each of the following sequences, determine the associated z-transform, includ-
ing the ROC. Make use of properties of the z-transform wherever possible.
(a) x1[n] = (I)"u[n]
(b) X2[n] = (-3)"u[-n]
(c) xsn] = xi[n] + x2[n]
(d) X4[n] = xl[n - 5]
(e) x5[n] = xi[n + 51
(f) eX[n] = (I)"u[n]
(g) x 7[n] = xi[n] * x 6 [n]

P23.2
A causal LTI system is described by the difference equation
y[n] - 3y[n - 1] + 2y[n - 2] = x[n]
(a) Find H(z) = Y(z)/X(z). Plot the poles and zeros and indicate the ROC.
(b) Find the unit sample response. Is the system stable? Justify your answer.
(c) Find y[n] if x[n] = 3"u[n].
(d) Determine the system function, associated ROC, and impulse response for all
LTI systems that satisfy the preceding difference equation but are not causal.
In each case, specify whether the corresponding system is stable.

P23.3
Carry out the proof for the following properties from Table 10.1 of the text (page
654).
(a) 10.5.2
(b) 10.5.3
(c) 10.5.6 [Hint: Consider
dX(z) x n
dz dzn=-_ ]

P23.4

Consider the second-order system with the pole-zero plot given in Figure P23.4. The
poles are located at z = re", z = re- ", and H(1) = 1.

P23-1
Signals and Systems
P23-2

Im
z plane

/ r \
Re

Figure P23.4

(a) Sketch |H(e'")I as 0 is kept constant at w/ 4 and with r = 0.5, 0.75, and 0.9.
(b) Sketch H(ein) as r is kept constant at r = 0.75 and with 0 = r/4, 27r/4, and
37r/ 4.

P23.5

Consider the system function

H(z) =z
(Z - = )(z - 2)

(a) Sketch the pole-zero locations.

(b) Sketch the ROC assuming the system is causal. Is the system stable?
(c) Sketch the ROC assuming the system is stable. Is the system causal?
(d) Sketch the remaining possible ROC. Is the corresponding system either stable or
causal?

P23.6
Consider the continuous-time LTI system described by the following equation:
d 2 y(t) dy(t) dx(t)
+5 +6ytt)=x(t)+2 d
dt2 dt dt
(a) Determine the system function He(s) and the impulse response he(t).
(b) Determine the system function Hd(z) of a discrete-time LTI system obtained
from He(s) through impulse invariance.
(c) For T = 0.01 determine the impulse response associated with Hd(z), hd[n].
(d) Verify that hd[n] = hc(nT) for all n.
 Mapping Continuous-Time Filters to Discrete-Time Filters / Problems
P23-3

Optional
Problems
P23.7

Consider a continuous-time LTI filter described by the following differential

equation:
dy(t) + 0.5y(t) = x(t)
dt
A discrete-time filter is obtained by replacing the derivative by a first forward dif-
ference to obtain the difference equation
y[n + 1] - y[n] + 0.5y[n] = x[n]
T
Assume that the resulting system is causal.
(a) Determine and sketch the magnitude of the frequency response of the continu-
ous-time filter.
(b) Determine and sketch the magnitude of the frequency response of the discrete-
time filter for T = 2.
(c) Determine the range of values of T (if any) for which the discrete-time filter is
unstable.

P23.8
Consider an even sequence x[n] (i.e., x[n] = x[-n]) with rational z-transform X(z).
(a) From the definition of the z-transform show that

Xz) = X (1)
(b) From your result in part (a), show that if a pole (or a zero) of X(z) occurs at
z = zo, then a pole (or a zero) must also occur at z = 1/zo.
(c) Verify the result in part (b) for each of the following sequences:
(i b[n + 1] + b[n - 11
(ii) 6[n + 1] - !Rkn] + b[n - 1]
(d) Consider a real-valued sequence y[n] with rational z-transform Y(z).
(i) Show that Y(z) = Y*(z*).
(ii) From part (i) show that if a pole (or a zero) of Y(z) occurs at z = zo, then
a pole (or a zero) must also occur at z = z*.
(e) By combining your result in part (b) with that in part (d), show that for a real,
even sequence, if there is a pole (or a zero) of H(z) at z = pej", then there is also
a pole (or a zero) of H(z) at z = pe-j", at z = (1/p)e'", and at z = (1/p)e -".
Signals and Systems
P23-4

P23.9

In Section 10.5.5 of the text we stated the convolution property for the z-transform.
To prove this property, we begin with the convolution sum expressed as

x 3 [n] = x 1 [n] * x 2[n] = E x 1 [k]x2[n - k] (P23.9-1)

k= -o

(a) By taking the z-transform of eq. (P23.9-1) and using eq. (10.3) of the text (page
630), show that

XA(z) = x1[k]X 2(z),

k= -o0

where X 2(z) = Z{x2[n - k]}.

(b) Using the result in part (a) and the time-shifting property of z-transforms, show
that

XA(z) = X 2(z) x1[k]z-

k=-o

(c) From part (b), show that

XA(z) = X 1(z)X 2(z).

P23.10

Consider a signal x[n] that is absolutely summable and its associated z-transform
X(z). Show that the z-transform of y[n] = x[n]u[n] can have poles only at the poles
of X(z) that are inside the unit circle.

P23.11

Let hc(t), sc(t), and He(s) denote the impulse response, step response, and system
function, respectively, of a continuous-time, linear, time-invariant filter.
Let hd[n], sd[n], and Hd(z) denote the unit sample response, step response, and
system function, respectively, of a discrete-time, linear, time-invariant filter.

(a) If hd[n] = hc(nT), does Sd[n] = 2 hc(kT)?

k=

(b) If Sd[n] = sc(nT), does hd[n] = hc(nT)?

P23.12

Consider a continuous-time filter with input xc(t) and output yc(t) that is described
by a linear constant-coefficient differential equation of the form
N dk yc(0 M dkxc(t) (P23.12-i)
=0
a k (k=O bk dtk
 Mapping Continuous-Time Filters to Discrete-Time Filters / Problems
P23-5

The filter is to be mapped to a discrete-time filter with input x[n] and output y[n]
by replacing derivatives with central differences. Specifically, let V("k{x[n]} denote
the kth central difference of x[n], defined as follows:
V(43{x[n]} = x[n]
V('[x[n]} x[n + 1] - x[n - 1]

V(k){xln]} V(l){V(k-'){x[nJ}}

The difference equation for the digital filter obtained from the differential equation
(P23.12-1) is then
N M

E
k=O
akV k){y[n]) = Z bkV(k){x[nI}
k=0

(a) If the transfer function of the continuous-time filter is Hc(s) and if the transfer
function of the corresponding discrete-time filter is Hd(z), determine how Hd(z)
is related to He(s).
(b) For the continuous-time frequency response He(jw), as indicated in Figure
P23.12-1, sketch the discrete-time frequency response Hd(ej0 ) that would result
from the mapping determined in part (a).

- 2 4

Figure P23.12

(c) Assume that Hc(s) corresponds to a causal stable filter. If the region of conver-
gence of Hd(z) is specified to include the unit circle, will Hd(z) necessarily cor-
respond to a causal filter?

P23.13
In discussing impulse invariance in Section 10.8.1 of the text, we considered Hc(s)
to be of the form of eq. (10.84) of the text with only first-order poles. In this problem
we consider how the presence of a second-order pole in eq. (10.84) would be
reflected in eq. (10.87) of the text. Toward this end, consider Hc(s) to be
A
He(s) (S - so)2
(a) By referring to Table 9.2 of the text (page 604), determine hc(t). (Assume
causality.)
Signals and Systems
P23-6

(b) Determine hd[n] defined as hd[n] = hc(nT).

(c) By referring to Table 10.2 of the text (page 655), determine Hd(z), the z-trans-
form of hd[n].
(d) Determine the system function and pole-zero pattern for the discrete-time sys-
tem obtained by applying impulse invariance to the following continuous-time
system:
1
He(s) =
(s + 1)(s + 2)2
MIT OpenCourseWare
http://ocw.mit.edu

Resource: Signals and Systems

Professor Alan V. Oppenheim

The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
22 The z-Transform
Solutions to
Recommended Problems
S22.1

(a) The z-transform H(z) can be written as

H(z) = z
z -2
Setting the numerator equal to zero to obtain the zeros, we find a zero at z = 0.
Setting the denominator equal to zero to get the poles, we find a pole at z = 1.
The pole-zero pattern is shown in Figure S22.1.

z plane

ROC: IzI > R

2

Figure S22.1

(b) Since H(z) is the eigenvalue of the input z' and the system is linear, the output
is given by

y[n] = 1 3) +3 1 (2)"
= 3()" + 4(2)"

S22.2
(a) To see if x[n] is absolutely summable, we form the sum
N-I N-i

SN = x[nhl 2' = 1 2
n=o n=O

Since limN-oSN diverges, x[n] is not absolutely summable.

(b) Since x[n] is not absolutely summable, the Fourier transform of x[n] does not
converge.
N

N(N- "
(C) SN
n=o(2

S22-1
Signals and Systems
S22-2

limN-oN is finite for Ir > 2. Therefore, the Fourier transform of r-"x[n] con-
verges for Ir| > 2.

(d) X(z) = 2"z~" = E (2z- 1)

n=O n.O

__ 1
= 2z- for 12z~'I < 1
Therefore, the ROC is IzI > 2.
1
(e) X 1(e") =
1 e -s
Therefore, x 1 [n] = (2)"u[n].

S22.3
(a) Since x[n] is right-sided, the ROC is given by Iz I > a. Since the ROC cannot
include poles, for this case the ROC is given by Iz I > 2.
(b) The statement implies that the ROC includes the unit circle Iz = 1. Since the
ROC is a connected region and bounded by poles, the ROC must be
5 < IzI < 2
(c) For this situation there are three possibilities:
(i) Izi < 3
(ii) -1 < I z| < 2

(iii) IzI > 2

(d) This statement implies that the ROC is given by |z < 1.

S22.4

(a) (i) X, (z) = x I~ -n T (-2 )f

n=-oon=O

1
1 - '

with an ROC of < 1, or IzI > 2.

-1 2z
(ii) X2 (z) = E (i)z~

Letting n = -m, we have

X 2(z) = - l()--m'
m= 1
GO 2z
= - L (2z) m
= - 1-2z
m1

1 - z '
with an ROC of 12z| < 1, or Iz < i.
 The z-Transform / Solutions
S22-3

(b) (i)
Im (z

ROC

Re lz

Figure S22.4-1

(ii)
ImZl

ROC
Re (z)

lzl= 1

Figure S22.4-2

2z
(c) (i) X 3(z) 2 z- = 2 (1 > 1, as shown
=
n=O
-1z) -Z-1 .The ROC is Izi
in Figure S22.4-3.

-1i

(ii) X 4(z) = - 2"z~ = - Z 2-"z

n=

z/2 z
n=Al
1 -(z/2) z -2'
with an ROC of Iz/21 < 1, or IzI < 2, shown in Figure S22.4-4.
Signals and Systems
S22-4

Im

/j

A. Re Iz)

Izl=2

Figure S22.4-4
(d) For the Fourier transform to converge, the ROC of the z-transform must include
the unit circle. Therefore, for x 1 [n] and x 4[n], the corresponding Fourier trans-
forms converge.

S22.5
Consider the pole-zero plot of H(z) given in Figure S22.5-1, where H(a/2) = 1.

z plane

K ' 1 to 5zeros
Figure S22.5-1

(a) When H(z) = z/(z - a), i.e., the number of zeros is 1, we have

H(ej) = cos 9 + j sin 0_

(cos 0 - a) + j sin 12
Therefore,
2
|H(e )| =
1 + a - 2a cos 0'

and we can plot IH(ejQ)I as in Figure S22.5-2.

 The z-Transform / Solutions
S22-5

IH(en2)I

Figure S22.5-2

When H(z) =z/(z - a), i.e., the number of zeros is 2, we have

- cos 20 + j sin 22
(cos 0 - a) + j sin 2
Therefore,

|H(es")| = 1 + a2 - 2a cos 0
Hence, we see that the magnitude of H(eia) does not change as the number of
zeros increases.
(b) For one zero at z = 0, we have

H(z)= z-a'

H(e") = e
e- a

We can calculate the phase of H(e ") by [Q - 4 (denominator)]. For two zeros
at 0, the phase of H(ein) is [29 - 4 (denominator)]. Hence, the phase changes
by a linear factor with the number of zeros.
(c) The region of the z plane where IH(z) I = 1 is indicated in Figure S22.5-3.

A z plane

+- Re z =

Figure S22.5-3
Signals and Systems
S22-6

S22.6

(a) (3)"[n]
n=O

= (3z)- 1 z
n=O i e z Z
Therefore, there is a zero at z = 0 and a pole at z = 1, and the ROC is

-<1 or IzI> ,
3,
as shown in Figure S22.6-1.

z plane

WZ
/
A, i - - - - -

7
37

I Figure S22.6-1

(b) b[n +1] L bS[n + 1]z-" = z,

with the ROC comprising the entire z-plane, as shown in Figure S22.6-2.

Im z plane

Re

Figure S22.6-2

S22.7

(a) Using long division, we have

X(z) = 1 -!Z- + !z - z- 3 +
We recognize that
x[n] = (-i)"u[n]
 The z-Transform / Solutions
S22-7

(b) Proceeding as we did in part (a), we have

x[n] = (-I"u[n]

(c) X(z) = -a
-az

S_ a2)

(1
a 1-az

1 a2)z-
1 a2
a 1-a-IZ-1
Therefore,

a-[n] u[n
x[n) = - - aa I -

Solutions to
Optional Problems
S22.8
(a) 1{r-"x[n]} = 5{r-u[n]}

= r~"e -jan
n=O

= Z (re+yu-n
n=O

For the sum to converge, we must have

< 1
Thus, Ir| > 1.

(b) (i)
Im

r =1z plane

Re

Figure S22.8-1
Signals and Systems
S22-8

(ii)
Im

r= z plane

Re

Figure S22.8-2

(iii)

r= 3 z plane

Figure S22.8-3

(c)

Figure S22.8-4

S22.9
(a) The inverse transform of
1
1 -i -
 The z-Transform / Solutions
S22-9

is (i)"u[n]. But from the given relation, we have

2-"x[n] = (L)"u[n],

x[n] = (I)"u[n]
3 + 2z- 1 3 + 2z-'
(b) 2 + 3z-' + z-2
(2 + z-1)(1 + z-1)
1 1
2+z- 1 + z-'
2 1
1+ iz-' 1 + z- 1
= -1)u[n] + (-1)"u[n]

S22.10
Az - a)
(a) H(z) = ( - with A a constant
(1 - azTe)
Therefore,
H A(e -j' - a)
(e' ) 1 - ae j
I A 2(e ' a)(e" - a)
=
|H~es")2 a) = A2
(1 - ae-3)(1 - aej")

and thus,
|H(e'")| = Al
(b) (i) |v 112 = 1 + a2 - 2a cos 9
1 2
(ii) |v 2 12 = 1 + -2 - - cos 9
a a
1
= - (a 2 + 1 - 2a cos Q)
1
= -| vil2

S22.11

In all the parts of this problem, draw the vectors from the poles or zeros to the unit
circle. Then estimate the frequency response from the magnitudes of these vectors,
as was done in the lecture. The following rough association can be made:
(a) (i)
(b) (ii)
(c) (iv)
(d) (iii)
(e) (iv)
Signals and Systems
S22-10

S22.12
(a) x[n] = (I)"[u[n] - u[n - 101]. Therefore,
9
X(z) = E(i)"Z--"
n=o
1 0
-1(zY
= (2z)" = 1 (2z)-
n=O 1 (2z)-1
10 11

= |z| > 0,
zl
z9(z- (1)10I">0
- 2)

shown in Figure S22.12-1. The Fourier transform exists.

(b) x[n] = (.)In' = (I)"u[n] + (.) -u[-n - 1]. But
(l\n z z 1
u[n] , IzI> ,

and
(1)-n
u[-n - ]- z < 2
2 ~z - 2'
Summing the two z-transforms, we have
-2z
X(z) = 2 1 < |zI < 2
(z 2)D(z - 2) 2
(See Figure S22.12-2.) The Fourier transform exists.

(c) x[n] = 7()" cos + u[n]

Therefore,

X(z) = 76 cos + 4z

7 00 n j(2r/6)n+(,r/4)] j(2w/6)n+(-/4)]
Z
- l + e n
2n=0 3

[ej/4 ei(2/6) -1 -jr/

4
1 -j(2r/6)-1)]

7 e _ e
6__ jw/4

le(2r/6)Z- (2v/6) 1
7z e e -j4

,r
2z cosir
(2 / 6)

2 cos (2,r
z
Z ire-
,r)
(e j(21/6)
7

7z 4 3 6 4 1
z- I (2r/6>)(z - j(2/6) where |z| >

The pole-zero plot and ROC are shown in Figure S22.12-3. Clearly, the Fourier
transform exists.
 The z-Transform / Solutions
S22-11

9th order pole

pole -zero
cancel

Figure S22.12-1

z plane

Figure S22.12-2

Figure S22.12-3
Signals and Systems
S22-12

1 - z' 0 z10 -1
(d) X(z) = z- - - 1 = z -
.=, 1 - Z' z"(z - 1

The ROC is all z except z = 0, shown in Figure S22.12-4. The Fourier transform
exists.

9th-order pole
.1.0 xPC 1 -1 1-

pole -zero
cancel

Figure S22.12-4

S22.13
(a) From

log(1 - w)= IwI < 1,

i=1
we find
(2z)'
log(1 - 2z) = 12z| < 1

z", |z| < ,

2-n
x[n] = s
0, n >- 0
(b) We solve this similarly to the way we solved part (a).
(Oz-')'f
log 1 12 - 1-1 Z11<1
1
z-< 1
n 2 2

01
n > 0,
x[n]
=t
0, n s 0
MIT OpenCourseWare
http://ocw.mit.edu

Resource: Signals and Systems

Professor Alan V. Oppenheim

The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
23 Mapping Continuous-Time Filters
to Discrete-Time Filters
Solutions to
Recommended Problems
S23.1

(a) X,(z) = - 1 -2z 2Z <1,

so the ROC is IzI > i.

0
(b) X 2(z) = (-3)"z-

= (-3)-"z" 1 z | 3-'z| < 1,

n=O +3

so the ROC is IzI < 3.

We can also show this by using the property that

x[-n] < X(z-1)

Letting x[-n] = x 2[n], we have

x[-n] = (-3)"u[-n],
x[+n] = ()[n
1
X(z) = + iz'

Therefore,
1
X 2(z) = 1+ 1z3

and the ROC is IzI < 3.

(c) Using linearity we see that

X 3(z) = X 1(z) + X 2(z) 1 - iz- 1 + z

The ROC is the common ROC for XI(z) and X 2(z), which is - < z < 3, as shown
in Figure S23.1.

S23-1
Signals and Systems
S23-2

z plane

-3

Figure S23.1

(d) Using the time-shifting property

x[n - n 0] z "OX(z),
we have

X 4(z) = z- 5X 1(z) = -Z

Delaying the sequence does not affect the ROC of the corresponding z-trans-
form, so the ROC is Izi > i.
(e) Using the time-shifting property, we have

X 5(z) = z

and the ROC is Iz > .

(f) X 6(z) = =
n=0 3 1
and the ROC is I z-'| < 1, or Izi >i.
(g) Using the convolution property, we have

X7(z) =X1(z)X
6(z)= (1-1i) 1-z
and the ROC is Iz I > i, corresponding to ROC1 n ROC.

S23.2

(a) We have
y[n] - 3y[n - 1] + 2y[n - 21 = x[n]
Taking the z-transform of both sides, we obtain
Y(z)[1 - 3z-' + 2z -2 =X(z),

Y(z) _ 1 -2 _Z2 _z_2 z_2__

X(z) 1 - 3z- + 2z - 3z + 2 (z - 2)(z - 1)

 Mapping Continuous-Time Filters to Discrete-Time Filters / Solutions
S23-3

and the ROC is outside the outermost pole for the causal (and therefore right-
sided) system, as shown in Figure S23.2.

z plane

Figure S23.2

(b) Using partial fractions, we have

1 2 -1
H(z) = 1 = +
(1 - 2z- )(1 - z-1) 1 - 2z -1 1- z
By inspection we recognize that the corresponding causal h[n] is the sum of two
terms:
h[n] = (2)2"u[n] + (-1)1"u[n]
= 2"*'u[n] - u[n]
= (2"*' - 1)u[n].

The system is not stable because the ROC does not include the unit circle. We
can also conclude this from the fact that

E
n=O
12n+ 1 = o

(c) Since x[n] = 3"u[n],

1
X(z) = 3 1'
-
1 - 3z IzI > 3,

1 1
Y(z) = H(z)X(z) = _
(1 - 3z 1) (1 - 2z-1)(1 - z 1)

Using partial fractions, we have

-4 i
Y(z) = + + IzI > 3
1 - 3z-' 1 - 2z-' 1 - Z-1'
since the output is also causal. Therefore,
y[n] = (2)3"u[n] - (4)2"u[n] + in[n]
(d) There are two other possible impulse responses for the same
1
1
H(z) =
1 - 3z-' + 2z2
Signals and Systems
S23-4

corresponding to different ROCs. For the ROC Iz < 1 the system impulse
response is left-sided. Therefore, since
2 -1
H(z) = 1- 2z 1 z-'

then
h[n] = - (2)2"u[-n - 11 + (1)u[-n - 1]
= -2" 1u[-n - 1] + u[-n - 1]
For the ROC 1 < Iz I < 2, which yields a two-sided impulse response, we have
h[n] = -2"*lu[-n - 1] - u[n]
since the second term corresponding to -1/(1 - z-) has the ROC 1 < Iz|.
Neither system is stable since the ROCs do not include the unit circle.

S23.3
(a) Consider

X 1(z) = ( x[n - noiz-"

n= 0

Letting m = n - no, we have

X 1(z) = E x[m]z-m+ no)

m= -o
00

= z-"i' ( x [m]z-m
M = -_0

= z-oX(z)
It is clear that the ROC of Xi(z) is identical to that of X(z) since both require
that IU _0 x[n]z-" converge in the ROC.
(b) Property 10.5.3 corresponds to multiplication of x[n] by a real or complex expo-
nential. There are three cases listed in the text, which we consider separately
here.

(i) Xi(z) = E e x[n]Z--

- : x[n] (ze -. 7 )- n

= X(ze --io),

with the same ROC as for X(z).

(ii) Now suppose that

X 2(z)= zox[n]z

x[n] -

=X
 Mapping Continuous-Time Filters to Discrete-Time Filters / Solutions
S23-5

Letting z' = z/zo, we see that the ROC for X 2(z) are those values of z such
that z' is in the ROC of X(z'). If the ROC of X(z) is RO < Iz| < R 1, then
the ROC of X 2(z) isRozol < |zi < RjIzo|.
(iii) This proof is the same as that for part (ii), with a = zo.
(c) We want to show that
z dX(z)
nx[n] dz
Consider

X(z) = ( x[njz-"
n = -oo

Then
dX(z)
dz)= 700n--nx[niz-*-
dz n=_-o

= -z- n 1= -o ( nx[n]z-",

so
dX(z) = n
-Z dz nx[niz-',
dz n=_-o

which is what we wanted to show. The ROC is the same as for X(z) except for
possible trouble due to the presence of the z' term.

S23.4

(a)

IH(ein)|

r = 0.9
-r= 0.75
r = 0.5

0 iT 7 21r
4 2
Figure S23.4-1
Signals and Systems
S23-6

(b)

S23.5
(a)
Im z plane

()F X )( Re
3

Figure S23.5-1
 Mapping Continuous-Time Filters to Discrete-Time Filters / Solutions
S23-7

(b)

The ROC is IzI > 2. The system is not stable because the ROC does not include
the unit circle.
(c)

z plane

Figure S23.5-3

The ROC is 2 > Iz I > i, which for this case includes the unit circle. The corre-
sponding impulse response is two-sided because the ROC is annular. Therefore,
the system is not causal.
(d)
Im z plane

Re

Figure S23.5-4
Signals and Systems
S23-8

The remaining ROC does not include the unit circle and is not outside the out-
ermost pole. Therefore, the system is not stable and not causal.

S23.6
daytt) dytt)
(a) + 5 + 6y(t) = x(t) + 2 dxt)
dx(t)
dt2 dt dt
is the system differential equation. Taking Laplace transforms of both sides, we
have
Y(s)(s 2 + 5s + 6) = X(s)(1 + 2s),
so

Y(s) 1 + 2s
2
X(s) s + 5s + 6
1+2s 5 -3
(s + 3)(s + 2) s + 3 s + 2
Assuming the system is causal, we obtain by inspection
h,(t) = 5e - 3 u(t) - 3e -2'u(t)
(b) Using the fact that the continuous-time system function Ak/(s - s) maps to the
discrete-time system function Ak/(l - e*kz- 1) (see page 662 of the text), we
have
5 3
Hd(z) 1 e "z- - e -2TZI

(c) Suppose T = 0.01. Then

5 3
Hd(Z)
1 - e -0 0 3
z-'1- e ~o 02z- 1
0
Letting a = e - .
03
, b = e- 002, we have
5 _ -
33
Hd(z) =
1 - az' 1- bz-

So by inspection, assuming causality,

hd[n] = 5a'u[n] - 3b u[n]
(d) From part (a), we have
hc(t) = 5e - 3 tu(t) - 3e 2
u(t)
Replacing t by nT = 0.01n, we have
hc(nT) = 5e 0~ 0 3 u(0.01n) - 3e 0 02
"u(0.02n)
Letting a = e -0.03 and b = e -0.02 yields
he(nT) = 5a'u[n] - 3b u[n],
which agrees with the result in part (c).
 Mapping Continuous-Time Filters to Discrete-Time Filters / Solutions
S23-9

Solutions to
Optional Problems
S23.7
(a) The differential equation is
dy(t) + 0.5y(t) x(t)
dt
Taking the Laplace transform yields
Y(s)[s + 0.5] = X(s),

H(s) = Y(S) = 1
X(s) s+ 0.5 '
H(w) = 0
jF + 0.5'
which is sketched in Figure S23.7-1.

20 logH(w)
20lgI H(0)

slope
0- -20 dB/decade

-20-. i W
.5 5 50 ...

Figure S23.7-1

(b) y[n + 1] - y[n] + 0.5y[n] = x[n]

T
Taking the z-transform of both sides yields
1
1 (z - 1)Y(z) + 0.5Y(z) = X(z),
T

Y(z) ( 0.5 + zT X(z)

Letting T = 2 yields
Y(z) 2
= Hd(z) = -
Z(z) z
Now since

|Hd(ej0 )I = |Hd(z)I I
Signals and Systems
S23-10

we have
|Hd(e j)| = 2, for all Q,
which is an all-pass filter and is sketched in Figure S23.7-2.

(c) HA(z) =
1 1
0.5- + - z
T T

(0.5T - 1) + z
The pole is located at zo = -(0.5T - 1) and, since we assume causality, we
require that the ROC be outside this pole. When the pole moves onto or outside
the unit circle, stability does not exist. The filter is unstable for
Izo l 1 or I-(0.5T - 1)1 1,
|0.5T - 11 1,
T 4
Therefore, for T > 4, the system is not stable.

S23.8

(a) X(z) = x[n]z-

X(z-1) = Z x[nZ"

Letting m = -n, we have

X(z 1) =
M= -0
x[-m]z- = ( x[mjz-"' = X(z)

f (z - ak)
(b) X(z) = A k
J (z - bk)

from the definition of a rational z-transform. Now

J1 (z-1 - ak)
X(z 1) = A k
1 (z- - bk)
k
Each pole (or zero) at zo in X(z) goes to a pole (or zero) z-1 inX(z-1). This implies
that zo = 1 or that X(z) must have another pole (or zero) at z--1.
 Mapping Continuous-Time Filters to Discrete-Time Filters / Solutions
S23-11

(c) (i) x[n] = b[n + 1] + b[n - 1],

_ _ z2 + _ (z + j)(z - j)
_ +
z z
The zeros are at z = j, 1/j, and the poles are at z = 0, z = oo.
(ii) x[n] = b[n + 1] - U[n] + 6[n - 1],
X(z) = z - 2 + z-
2 z + 1 (z - 1)(z - 2)

The zeros are at z = i, 2, and the poles are at 0, oo.

(d) (i) Y(z) = E y[n]z-",

Y*(z) = y[n]z- = L y[n](z*)" = Y(z*),

Y*(z*) = [ y[n]z-" = Y(z)

n= -00
(ii) Since Y(z) is rational,
fl(z - ak)
Y(z) = A k

f(z - bk)

Now if a term such as (z - ak) appears in Y(z), a term such as (z* - ak)
must also appear in Y(z). For example,
Y(z) = (z - ak)(z - ak),
*(z*) = [(z* - ak)(z - ak)*]
= (z - a*)(z - ak) = Yz)

So if a pole (or zero) appears at z = ak, a pole (or zero) must also appear
at z = a* because
(z - ak) = 0 = z = ak

(e) Both conditions discussed in parts (b) and (d) hold, i.e., a real, even sequence is
considered. A pole at z = z, implies a pole at 1/z, from part (b). The poles at
z = z, and z = 1/z, imply poles at z = z* and z = (1/zp)* from part (d). There-
fore, if z, = pej", poles exist at

pei"
1 =1
p
(pej")* = pe -',
( *1pei p

S23.9

(a) X 3 (z) = X3[n]z~

00 00

-
k=-
Y x1[k] [
( x2[n
n=--oo
- k]z-"
oO

= x1[kli2(z),
k= -oo
Signals and Systems
S23-12

where

X 2(z) = E X2[n - kz--

= Z{x2[n - k]}
(b) Z{x 2 [n - k) = z-kX 2(z) from the time-shifting property of the z-transform, so

XA(z) = ( x 1 [k]z-kX2(z)
k=-oo

(c) XA(z) = X 2(z) T x 1[k]z k

k=-o
= X 1 (z)X 2(z)

S23.10
Consider x[n] to be composed of a causal and an anticausal part:
x[n] = x[n]u[-n - 1] + x[n]u[n]
Let
x 1 [n] = x[n]u[-n - 1],
X2[n] = x[n]u[n],
so that
x[n] = x 1 [n] + x 2[n]
and
X(z) = X 1(z) + X 2(z)
It is clear that every pole of X 2(z) is also a pole of X(z). The only way for this not
to be true is by pole cancellation from X1 (z). But pole cancellation cannot happen
because a pole ak that appears in X 2(z) yields a contribution (ak)"u[n], which cannot
be canceled by terms of x 1[n] that are of the form (bk )U[-n - 11.
From the linearity property of z-transforms, if
y[n] = y 1 [n] + y 2[n]
then
Y(z) = Y1(z) + Y 2(z),
with the ROC of Y(z) being at least the intersection of the ROC of Yi(z) and the ROC
of Y 2(z). The "at least" specification is required because of possible pole cancella-
tion. In our case, pole cancellation cannot occur, so the ROC of X(z) is exactly the
intersection of the ROC of X 1(z) and the ROC of X 2(z).
Now suppose X 2(z) has a pole outside the unit circle. Since x2[n] is causal, the
ROC of X 2(z) must be outside the unit circle, which implies that the ROC of X(z)
must be outside the unit circle. This is a contradiction, however, because x[n] is
assumed to be absolutely summable, which implies that X(z) has an ROC that
includes the unit circle.
Therefore, all poles of the z-transform of x[n]u[n] must be within the unit
circle.
 Mapping Continuous-Time Filters to Discrete-Time Filters / Solutions
S23-13

S23.11
(a) If h[n] = hc(nT), then

Sd[n] = E hc(kT)
k= -o

The proof follows.

Sd[n] = ( u[n - k]hd[k]

k= -o

= E hd[k],
k= -o

but hd[k] = he(kT), so

n

Sd[n] = E h,(kT)
k= -00

(b) If Sd[n] = se(nT), then hd[n] does not necessarily equal he(nT). For example,
hc(t) = e ~4'U(t),
sc(t) = f e -au(r)u(t -r) d-

= t - "'dr = - (1 - e -a'), t : 0
o 1a

Sd[n] = sc(nT) = - (1 - e -a"T), n >- 0

a
However,
n

Sd[n] = Z
k= -oo
hd[k],

so
Sd[n] - sd[n - 1] = hd[n]

and, in our case, for n 0,

1 1
hd[n] = - (1 - e -a"T) - e(1
e ~-"a"-)T)
a a
= 1 e -an(ea" - 1)
a
But, for n > 0,
he(nT) = e -a"Tu(nT)
1
#- e -anT(eae 1)
a

S23.12
(a) From the differential equation

( aks) Y(s) = ( bksk)X(S),

k-0 k=0
Signals and Systems
S23-14

we have
M
Y) bksk
Ys) = He(s) = k=O
X(s) N k
( aks"
k=0

Now consider

yldn] V"{x[n]} -x[n + 1] - x[n - 1]

2
Yi(z) = Z{y 1 [n]} = Z{Vo> {x[n]}) z -z X(z),

Y2[nJ = V(2>{X[n]= - y 1 [n + 1 ~ y1[n - 1]

2
Y 2 (z) = Z Z Yi(z) = (z 1)2X(z)
By induction,

Z{Vlk){X[n} =(z z1 X(Z)

Therefore,
N _(- k M b(k
a Y(z) =Yb ( X) )
k=0 k=

z - z-
b'
Hd(z) = Yz
X(z) N _(zk
(ak
k=0

=He,(s)

(b) He(s) = Hd(Z)

s =(z -z-/
from part (a). Consider s = jw, z = ej'. So

.' e j" - e -j"

= 2
and, thus, w = sin 0 is the mapping between discrete-time and continuous-time
frequencies. Since H(w) = w for IwI < 1, Hd(ej") is as indicated in Figure
S23.12.
 Mapping Continuous-Time Filters to Discrete-Time Filters / Solutions
S23-15

Hd(ejR)

7T T7T iT
2 -2

-1-

Figure S23.12

(c) From part (a) we see that

Hd(z) = Hd(z )

and that Hd(z) is a rational z-transform.

P
Afl (z - zo.)Mi

Hd(z) = Q'
] N,(z - zpi)Ni

Therefore, if a term such as (z - zo,) appears, (z-' - zo,) must also appear. If
Hd(z) has a pole within the unit circle, it must also have a pole outside the unit
circle. If the ROC includes the unit circle, it is therefore not outside the outer-
most pole (which lies outside the unit circle) and, therefore, Hd(z) does not cor-
respond to a causal filter.
Consider
1
He(s) =

corresponding to a stable, causal h,(t).

1
Hd(z) = He(s) =s=(z-z1)2
Z - Z 1
2 2
2z 2z
2
z +Z 1 [)(-1 + V5) (-1 - V5)
2 1[2

so poles of z are at 0.618, -1.618. Therefore, Hd(z) is not causal if it is assumed

stable because stability and causality require that all poles be inside the unit
circle.
Signals and Systems
S23-16

S23.13
(a) We are given that
A
He(s) = (S 2

From Table 9.2 of the text (page 604), we see that hc(t) = Ate'Ou(t).
To verify, consider

S O= f .e'O'utt)e -"dt,
d
1 d * sotu(t)e Stdt]
ds (s s) ds u

=SS) teso'ult)e~"dt

Therefore,
.Cr
tesotu(t) ( 2
(s -so

(b) hd[n] = hc(nT) = AnTeson Tu[n]

(c) Hd(z) = L hd[nlz~" = AT : nesonTz-n

From Table 10.2 of the text (page 655),
z az-
na~u[n]~--
(1 - az-1)2
-

This can be verified:

1 _

1 _ = a u[n]z-"
d 1) d *0
-z = C-z E a u[n]z-"
dz -- 1 z dz =
-az-2

= (-n)anu[n]z-"-
(1 - az )

= natu[n]z"
(1 -z) ,= _

In our case, a = e'oT, so

ATe soTz-
Hd(Z) = ( ez)2

(d) Hc(s) = = + +
(s+ 1)(s + 2)2 s + 1 s + 2 (s + 2)2
Using the first-order pole result for 1/(s + 1) and - 1/(s + 2) and the second-
order pole result for - 1/(s + 2)2, we have
1 1 Te 2Tz-
Hd(z) = 1 - e -Tz 1 - e 2 Tz- - e-2rz 1)2
1 e-T- (1
After some algebra, we obtain
z[z( -e -2T + e- - Te-2) + e -4T - e -3 T+ Te-3T]
Hd(z) = (z - e-T)(z - e -2T)2
 Mapping Continuous-Time Filters to Discrete-Time Filters / Solutions
S23-17

The corresponding pole-zero pattern is shown in Figure S23.13.

Im z plane

double pole

cl 1 C41 vRe
e-2T e-T

e-4T _ e-3 T (1 - T)
zero at z = -
2
e-T e T (1+T)

Figure S23.13
MIT OpenCourseWare
http://ocw.mit.edu

Resource: Signals and Systems

Professor Alan V. Oppenheim

The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.