21 5 SMPLD Fns
21 5 SMPLD Fns
Introduction
A sequence can be obtained by sampling a continuous function or signal and in this Section we
show first of all how to extend our knowledge of z-transforms so as to be able to deal with sampled
signals. We then show how the z-transform of a sampled signal is related to the Laplace transform
of the unsampled version of the signal.
HELM (2008): 85
Section 21.5: Sampled Functions
1. Sampling theory
If a continuous-time signal f (t) is sampled at terms t = 0, T, 2T, . . . nT, . . . then a sequence of
values
{f (0), f (T ), f (2T ), . . . f (nT ), . . .}
is obtained. The quantity T is called the sample interval or sample period.
f (t)
------
t
T 2T nT
Figure 18
In the previous Sections of this Workbook we have used the simpler notation {fn } to denote a
sequence. If the sequence has actually arisen by sampling then fn is just a convenient notation for
the sample value f (nT ).
Most of our previous results for z-transforms of sequences hold with only minor changes for sampled
signals.
So consider a continuous signal f (t); its z-transform is the z-transform of the sequence of sample
values i.e.
X∞
Z{f (t)} = Z{f (nT )} = f (nT )z −n
n=0
We shall briefly obtain z-transforms of common sampled signals utilizing results obtained earlier. You
may assume that all signals are sampled at 0, T, 2T, . . . nT, . . .
Unit step function
1 t≥0
u(t) =
0 t<0
Since the sampled values here are a sequence of 1’s,
1
Z{u(t)} = Z{un } =
1 − z −1
z
= |z| > 1
z−1
86 HELM (2008):
Workbook 21: z-Transforms
®
Ramp function
t t≥0
r(t) =
0 t<0
The sample values here are
{r(nT )} = {0, T, 2T, . . .}
z
The ramp sequence {rn } = {0, 1, 2, . . .} has z-transform .
(z − 1)2
Tz
Hence Z{r(nT )} = since {r(nT )} = T {rn }.
(z − 1)2
Task
Obtain the z-transform of the exponential signal
−αt
e t≥0
f (t) =
0 t < 0.
[Hint: use the z-transform of the geometric sequence {an }.]
Your solution
Answer
The sample values of the exponential are
{1, e−αT , e−α2T , . . . , e−αnT , . . .}
i.e. f (nT ) = e−αnT = (e−αT )n .
z
But Z{an } =
z−a
z 1
∴ Z{(e−αT )n } = −αT
= −αT
z−e 1−e z −1
HELM (2008): 87
Section 21.5: Sampled Functions
Sampled sinusoids
Earlier in this Workbook we obtained the z-transform of the sequence {cos ωn} i.e.
z 2 − z cos ω
Z{cos ωn} =
z 2 − 2z cos ω + 1
Hence, since sampling the continuous sinusoid
f (t) = cos ωt
yields the sequence {cos nωT } we have, simply replacing ω by ωT in the z-transform:
Task
Obtain the z-transform of the sampled version of the sine wave f (t) = sin ωt.
Your solution
Answer
z sin ω
Z{sin ωn} =
z 2 − 2z cos ω + 1
∴ Z{sin ωt} = Z{sin nωT }
z sin ωT
= 2
z − 2z cos ωT + 1
88 HELM (2008):
Workbook 21: z-Transforms
®
Shift theorems
These are similar to those discussed earlier in this Workbook but for sampled signals the shifts are
by integer multiples of the sample period T . For example a simple right shift, or delay, of a sampled
signal by one sample period is shown in the following figure:
f (nT )
t
T 2T 3T
f (nT − T )
t
T 2T 3T 4T
Figure 19
The right shift properties of z-transforms can be written down immediately. (Look back at the shift
properties in Section 21.2 subsection 5, if necessary:)
If y(t) has z-transform Y (z) which, as we have seen, really means that its sample values {y(nT )}
give Y (z), then for y(t) shifted to the right by one sample interval the z-transform becomes
Z{y(t − T )} = y(−T ) + z −1 Y (z)
The proof is very similar to that used for sequences earlier which gave the result:
Z{yn−1 } = y−1 + z −1 Y (z)
Task
Using the result
Z{yn−2 } = y−2 + y−1 z −1 + z −2 Y (z)
write down the result for Z{y(t − 2T )}
Your solution
Answer
Z{y(t − 2T )} = y(−2T ) + y(−T )z −1 + z −2 Y (z)
These results can of course be generalised to obtain Z{y(t − mT )} where m is any positive integer.
In particular, for causal or one-sided signals y(t) (i.e. signals which are zero for t < 0):
Z{y(t − mT )} = z −m Y (z)
Note carefully here that the power of z is still z −m not z −mT .
HELM (2008): 89
Section 21.5: Sampled Functions
Examples:
For the unit step function we saw that:
z 1
Z{u(t)} = =
z−1 1 − z −1
Hence from the shift properties above we have immediately, since u(t) is certainly causal,
zz −1 z −1
Z{u(t − T )} = =
z−1 1 − z −1
zz −3 z −3
Z{u(t − 3T )} = =
z−1 1 − z −1
and so on.
u(t − T )
t
T 2T 3T
u(t − 3T )
t
T 2T 3T 4T 5T
Figure 20
while the Laplace transform written symbolically as L{f (t)} is an integral, defined for a continuous
time function f (t), t ≥ 0 as
Z ∞
F (s) = f (t)e−st dt.
0
90 HELM (2008):
Workbook 21: z-Transforms
®
1
L{f (t)} = F (s) =
s+α
which has a (simple) pole at s = −α = s1 say.
As we have seen, sampling f (t) gives the sequence {f (nT )} = {e−αnT } with z-transform
1 z
F (z) = = .
1 − e−αT z −1 z − e−αT
The z-transform has a pole when z = z1 where
z1 = e−αT = es1 T
[Note the abuse of notations in writing both F (s) and F (z) here since in fact these are different
functions.]
Task
The continuous time function f (t) = te−αt has Laplace transform
1
F (s) =
(s + α)2
Firstly write down the pole of this function and its order:
Your solution
Answer
1
F (s) = has its pole at s = s1 = −α. The pole is second order.
(s + α)2
Now obtain the z-transform F (z) of the sampled version of f (t), locate the pole(s) of F (z) and
state the order:
Your solution
HELM (2008): 91
Section 21.5: Sampled Functions
Answer
Consider f (nT ) = nT e−αnT = (nT )(e−αT )n
Tz
The ramp sequence {nT } has z-transform
(z − 1)2
∴ f (nT ) has z-transform
T zeαT T ze−αT
F (z) = = (see Key Point 8)
(zeαT − 1)2 (z − e−αT )2
This has a (second order) pole when z = z1 = e−αT = es1 T .
We have seen in both the above examples a close link between the pole s1 of the Laplace transform
of f (t) and the pole z1 of the z-transform of the sampled version of f (t) i.e.
z1 = es1 T (1)
where T is the sample interval.
Multiple poles lead to similar results i.e. if F (s) has poles s1 , s2 , . . . then F (z) has poles z1 , z2 , . . .
where zi = esi T .
The relation (1) between the poles is, in fact, an example of a more general relation between the
values of s and z as we shall now investigate.
Key Point 19
The unit impulse function δ(t) can be defined informally as follows:
P! (t)
1
!
t
!
Figure 21
1
The rectangular pulse P (t) of width ε and height shown in Figure 21 encloses unit area and has
ε
Laplace transform
Z ε
1 −st 1
Pε (s) = e = (1 − e−εs ) (2)
0 ε εs
As ε becomes smaller Pε (t) becomes taller and narrower but still encloses unit area. The unit impulse
function δ(t) (sometimes called the Dirac delta function) can be defined as
92 HELM (2008):
Workbook 21: z-Transforms
®
The Laplace transform, say ∆(s), of δ(t) can be obtained correspondingly by letting → 0 in (2),
i.e.
1
∆(s) = lim (1 − e−εs )
ε→0 εs
(εs)2
1 − (1 − εs + − . . .)
= lim 2! (Using the Maclaurin seies expansion of e−εs )
ε→0 εs
(εs)2 (εs)3
εs − + + ...
= lim 2! 3!
ε→0 εs
= 1
Task
A shifted unit impulse δ(t − nT ) is defined as lim Pε (t − nT ) as illustrated below.
ε→0
P! (t − nT )
1
!
t
nT nT + !
Your solution
HELM (2008): 93
Section 21.5: Sampled Functions
Answer
Z nT +ε nT +ε
1 −st 1 −st
L{Pε (t − nT )} = e dt = −e
nT ε εs nT
1 −snT
− e−s(nT +ε)
= e
εs
1 −snT
= e (1 − e−sε ) → e−snT as ε → 0
εs
Hence L{δ(t − nT )} = e−snT (4)
which reduces to the result (3)
L{δ(t)} = 1 when n = 0
These results (3) and (4) can be compared with the results
Z{δn } = 1
Z{δn−m } = z −m
for discrete impulses of height 1.
Now consider a continuous function f (t). Suppose, as usual, that this function is sampled at t = nT
for n = 0, 1, 2, . . .
f (t)
------ t
T 2T 3T 4T
Figure 22
This sampled equivalent of f (t), say f∗ (t) can be defined as a sequence of equidistant impulses, the
‘strength’ of each impulse being the sample value f (nT )i.e.
∞
X
f∗ (t) = f (nT )δ(t − nT )
n=0
This function is a continuous-time signal i.e. is defined for all t. Using (4) it has a Laplace transform
∞
X
F∗ (s) = f (nT )e−snT (5)
n=0
If, in this sum (5) we replace esT by z we obtain the z-transform of the sequence {f (nT )} of samples:
∞
X
f (nT )z −n
n=0
94 HELM (2008):
Workbook 21: z-Transforms
®
Key Point 20
The Laplace transform
∞
X
F (s) = f (nT )e−snT
n=0
of a sampled function is equivalent to the z-transform F (z) of the sequence {f (nT )} of sample
values with z = esT .
z
1 1 1
z−1
z
t nT 1
(z − 1)2
2 2 T 2 z(z + 1)
t (nT ) 1
(z − 1)3
z
e−αt e−αnT |e−αT |
z − e−αT
z sin ωT
sin ωt sin nωT 1
z2 − 2z cos ωT + 1
z(z − cos ωT )
cos ωt cos nωT 1
z 2 − 2z cos ωT + 1
T ze−αT
te−αt nT e−αnT |e−αT |
(z − e−αT )2
e−αT z −1 sin ωT
e−αt sin ωt e−αnT sin ωnT |e−αT |
1 − 2e−αT z −1 cos ωT + e−2aT z −2
1 − e−αT z −1 cos ωT
e−αT cos ωt e−αnT cos ωnT |e−αT |
1 − 2e−αT z −1 cos ωT + e−2aT z −2
Note: R is such that the closed forms of F (z) (those listed in the above table) are valid for |z| > R.
HELM (2008): 95
Section 21.5: Sampled Functions
Table of z-transforms
fn F (z) Name
δn 1 unit impulse
δn−m z −m
z
un unit step sequence
z−1
z
an geometric sequence
z−a
z
eαn
z − eα
z sinh α
sinh αn
z2 − 2z cosh α + 1
z 2 − z cosh α
cosh αn
z 2 − 2z cosh α + 1
z sin ω
sin ωn
z2 − 2z cos ω + 1
z 2 − z cos ω
cos ωn
z 2 − 2z cos ω + 1
ze−α sin ω
e−αn sin ωn
z 2 − 2ze−α cos ω + e−2α
z 2 − ze−α cos ω
e−αn cos ωn
z 2 − 2ze−α cos ω + e−2α
z
n ramp sequence
(z − 1)2
z(z + 1)
n2
(z − 1)3
z(z 2 + 4z + 1)
n3
(z − 1)4
z
n
a fn F
a
dF
n fn −z
dz
96 HELM (2008):
Workbook 21: z-Transforms