20/07/2020, 19:18z_transform
Page 1 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
Discrete-Time Systems and the Z-Transform
ColophonThe source code for this page is dt_systems/2/z_transform.ipynb (https://github.com/cpjobling/eg-247-textbook/dt_systems/2/z_transform.ipynb).You can view the notes for this presentation as a webpage (HTML (https://cpjobling.github.io/eg-247-textbook/dt_systems/2/z_transform.html)).This page is downloadable as a PDF (https://cpjobling.github.io/eg-247-textbook/dt_systems/2/z_transform.pdf) file.
Scope and Background ReadingThis session introduces the z-transform which is used in the analysis of discrete time systems. As for theFourier and Laplace transforms, we present the definition, define the properties and give someapplications of the use of the z-transform in the analysis of signals that are represented as sequencesand systems represented by difference equations.
The material in this presentation and notes is based on Chapter 9(https://ebookcentral.proquest.com/lib/swansea-ebooks/reader.action?docID=3384197&ppg=329) of{cite} karris from the Required Reading List. Additional coverage is to be found in Chapter 13 of{cite} benoit from the Recommended Reading List.
Agenda
Introduction
The Z-Transform
Properties of the Z-Transform
20/07/2020, 19:18z_transform
Page 2 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
Some Selected Z-Transforms
Relationship between Laplace and Z-Transform
Stability Regions
Introduction
In the remainder of this course we turn our attention to how we model and analyse the behaviour of thecentral block in this picture.
Nature of the signals
20/07/2020, 19:18z_transform
Page 3 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
The signals we process in discrete time systems are sequences of values where is an index.
A sequence can be obtained in real-time, e.g. as the output of a ADC, or can be stored in digital memory;processed and re-stored; or processed and output in real-time, for example in digital music.
!["] "
Nature of the systemsThe input to a discrete time system is a squence of values The output is also a sequence The block represents the operations that convert into .This processing takes the form of a difference equationThis is analogous to the representation of continuous-time operations by differential equations.
!["]#["]
!["] #["]
Transfer function model of a DT systemIn CT systems we use the Laplace transform to simplify the analysis of the differential equationsIn DT systems the z-Transform allows us to simplify the analysis of the difference equationsIn CT systems application of the Laplace transform allows us to represent systems as transferfunctions and solve convolution problems by multiplicationThe z-transform provides analogous (https://www.google.co.uk/search?q=define:analogous) toolsfor the analysis of DT systems.
The Z-Transform! {$["]} = % (&) = $["]!
"=0
!&""
{% (&)} = $["] = % (&) '&!"1 12() ! &*"1
Sampling and the Z-TransformIn the last lecture we showed that sampling could be represented as the multiplication of a CT signal by aperiodic train of impulses:
(+) = !(+) ,(+ " " )!- !"=0
!.-
20/07/2020, 19:18z_transform
Page 4 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
By the sampling property of ,(+)
(+) = !(" ),(+ " " )!- !"=0
!.- .-
Using the Laplace transform pairs and we obtain:,(+) # 1 ,(+ " . ) # /"-.
(-) = " { !(" ),(+ " " )} = !(" )0- !"=0
!.- .- !
"=0
!.- /""-.-
By substitution of and representing samples as sequence :& = /-.- !(" ).- !["]
0(&) = !["]!"=0
!&""
20/07/2020, 19:18z_transform
Page 5 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
Properties of the z-Transform
Property Discrete Time Domain Transform
1 Linearity
2 Shift of
3 Left shift
4 Right shift
5 Multiplication by
6Multiplication by
7 Multiplication by
8Multiplication by
9 Summation in time
10 Time convolution
11 Frequency convolution
12 Initial value theorem
13 Final value theorem
For proofs refer to Section 9.2 of {cite} karris .
!
1 ["] + 2 ["] + $$1 $2 1 (&) + 2 (&) + $%1 %2
!["] ["]30$[" " 4] [" " 4]30 % (&)&"4
$[" " 4] % (&) + $[" " 4]&"4!"=0
4"1&""
$[" + 4] % (&) + $[" + 4]&4!
"="4
"1&""
1" $["]1" % ( )&1
/""-.- $["]/""-.- % ( &)/-.-
" "$["] "& % (&)''&
"2 $["]"2 "& % (&) + % (&)''&
&2 '2
'&2
$[4]!4=0
"% (&)&
& " 1
["] % ["]$1 $2 (&) (&)%1 %2
["] ["]$1 $2 ! ! (5) ( ) '51
)2(%1 %2
&5
5"1
$[0] = % (&)lim&&!
$["] = (& " 1)% (&)lim"&!
lim&&1
Some Selected Common z-Transforms
The Geometric Sequence
$["] = { 01"
" = "1, "2, "3, …" = 0, 1, 2, 3, …
20/07/2020, 19:18z_transform
Page 6 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
% (&) = $["] = =!"=0
!&"" !
"=0
!1" &"" !
"=0
!
(1 )&"1 "
After some analysis , this can be shown to have a closed-form expression1 2
% (&) = =11 " 1&"1
&& " 1
Notes
1. See Karris pp 9-12—9-13{cite} karris for the details.2. This function converges only if
and the region of convergence is outside the circle centred at with radius
.
!&! < !1!
& = 0
6 = !1!
Region of convergence
20/07/2020, 19:18z_transform
Page 7 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
The Delta Sequence
,["] = { 10
" = 0otherwise
! {,["]} = '(&) = ,["] = 1 + 0 = 1!"=0
!&"" !
"=1
!&""
,["] # 1
The Unit Step["] = {30
01
" < 0" ( 0
Z-Transform of Unit Step
20/07/2020, 19:18z_transform
Page 8 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
! { ["]} (&) = ["] =30 70 !"=0
!30 &"" !
"=0
!&""
This is a special case of the geometric sequence with so1 = 1(&) = =70
11 " &"1
&& " 1
Region of convergence is!&! > 1
Exponontial Decay Sequence$["] = ["]/""1.- 30
% (&) = = 1 + + + + $!"=0
!/""1-.- &"" /"1.- &"1 /"21.- &"2 /"13.- &"3
This is a geometric sequence with , so1 = /"1.-
! { ["]} = =/"1.- 301
1 " /"1.- &"1&
& " /"1.-
Region of convergence is! ! < 1/"1.- &"1
Ramp Function$["] = " ["]30
! {" ["]} = " = 0 + + 2 + 3 + $30 !"=0
!&"" &"1 &"2 &"3
We recognize this as a signal multiplied by for which we have the property["]30 ""$["] # "& % (&)'
'&
20/07/2020, 19:18z_transform
Page 9 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
After applying the property and some manipulation, we arrive at:
" ["] #30&
(& " 1)2
20/07/2020, 19:18z_transform
Page 10 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
z-Transform TablesAs usual, we can rely on this and similar analysis to have been tabulated for us and in practice we canrely on tables of transform pairs, such as this one.
f[n] F(z)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
,["] 1
,[" " 4] &"4
["]1"30 !&! > 1&
& " 1
["]30 !&! > 1&& " 1
( ) ["]/"1".- 30 ! ! < 1&& " /"1.-
/"1.- &"1
(cos "1 ) ["].- 30 !&! > 1" & cos 1&2 .-
" 2& cos 1 + 1&2 .-
(sin "1 ) ["].- 30 !&! > 1& sin 1.-
" 2& cos 1 + 1&2 .-
( cos "1 ) ["]1" .- 30 !&! > 1" 1& cos 1&2 .-
" 21& cos 1 +&2 .- 12
( sin "1 ) ["]1" .- 30 !&! > 11& sin 1.-
" 21& cos 1 +&2 .- 12
["] " [" " 4]30 30" 1&4
(& " 1)&4"1
" ["]30&
(& " 1)2
["]"230&(& + 1)(& " 1)3
[" + 1] ["]30&2
(& " 1)2
" ["]1" 301&
(& " 1)2
["]1""2301&(& + 1)(& " 1)3
"[" + 1] ["]1" 3021&2
(& " 1)3
Relationship Between the Laplace and Z-Transform
20/07/2020, 19:18z_transform
Page 11 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
Given that we can represent a sampled signal in the complex frequency domain as the infinite sum ofeach sequence value delayed by an integer multiple of the sampling time:
% (-) = $["]!"=0
!/""-.-
And by definition, the z-transform of such a sequence is:
% (&) = $["]!"=0
!&""
It follows that& = /-.-
And
- = ln &1.-
Mapping of s to zSince and are both complex variables, is a mapping from the -domain to the -domainand is a mapping from the to -domain.
- & & = /-.- - && = ln &/.- & -
Now, since- = 8 + )9
20/07/2020, 19:18z_transform
Page 12 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
& = = = !&!/8 +)9.- .- /8.- /)9.- /):
where
and
!&! = /8.-
: = 9 ..-
Introduction of sampling frequencyNow, since then or and = 1/.- $- = 2(/9- $- = /(2()$- 9- = 2(/.- 9-
We let
: = 9 = 9 = 2(.-2(9-
99-
Hence by substitution:& = /8+/)2(9/9-
Interpretation of the mapping s to zThe quantity defines a unit-circle in the -plane centred at the origin.And of course the term represents the (stability) boundary between the left- and right-half planesof the -plane.What are the consequences of this?
/)2(9/9- &8
-
Case I:
When we see that from
that
The left-half plane of the -domain maps into the unit circle in the -plane.Different negative values of map onto concentric circles with radius less than unity.
8 < 08 < 0
!&! = /8.-
!&! < 1- &
8
20/07/2020, 19:18z_transform
Page 13 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
Case II: 8 > 0
When we see that from
that
The right-half plane of the -domain maps outside the unit circle in the -plane.Different positive values of map onto concentric circles with radius greater than unity.
8 > 0
!&! = /8.-
!&! > 1
- &8
Case III: 8 = 0
When ,
and
All values of lie on the circumference of the unit circle.
8 = 0
!&! = 1
: = 2(99-
9
20/07/2020, 19:18z_transform
Page 14 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
Stability Region - s-Plane
20/07/2020, 19:18z_transform
Page 15 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
Stability Region - z-Plane
Frequencies in the z-DomainAs a consequence of the result for Case III above, we can explore how frequencies (that is is valuesof ) map onto the -plane.We already know that these frequencies will map onto the unit circle and by the anglesare related to the sampling frequency.Let's see how
- = ±)9 &: = 2(9/9-
Mapping of multiples of sampling frequency
20/07/2020, 19:18z_transform
Page 16 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
[radians/sec] [radians]
0 1 0
1
1
1
1
1
1
1
1
9 !&! :
/89- (/4/49- (/2
3 /89- 3(/4/29- (
5 /89- 5(/43 /49- 3(/27 /89- 7(/4
9- 2(
Mapping of s-plane to z-plane
20/07/2020, 19:18z_transform
Page 17 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
Mapping z-plane to s-planeThere is no unique mapping of to since
but for a complex variable
This is in agreement with the theoretical idea that in the frequency domain, sampling creates an infinitenumber of spectra, each of which is centred around .
& -- = ln &1
.-
ln & = ln & ± )2"(
±"9-
Frequency aliasingIt's worth observing that any stable complex pole in the -plane will have complexconjugate pair .Providing these poles will be mapped to the upper and lower half-plane of the -planerespectively.If , an upper-half plane pole will be mapped to the lower-half plane and will have aneffective frequency of .Similarly, its conjugate pair will move into the upper-half plane.
This is another way of looking at aliasing.
- - = "8 + )9- = "8 " )9
9 < /29- &
9 > /29-/2 " 99-
Also, any poles with frequency will also be aliased back into into the unit circle.9 > 9-
SummaryIntroductionThe Z-TransformProperties of the Z-TransformSome Selected Z-TransformsRelationship Between Laplace and Z-TransformStability Regions
Next session
The Inverse Z-Transform – an examples class
20/07/2020, 19:18z_transform
Page 18 of 18http://127.0.0.1:8888/nbconvert/html/dt_systems/2/z_transform.ipynb?download=false
Homework
Problems 1 to 3 in Section 9.10 Exercises of {cite} Karris explore the z-Transform
ReferencesSee Bibliography (/zbib).