Lecture Slides For Signals and Systems-3.0
Lecture Slides For Signals and Systems-3.0
Edition 3.0
Michael D. Adams
To obtain the most recent version of these lecture slides (with functional hyperlinks) or for additional
. ......................
information and resources related to these slides (including video lectures and errata), please visit:
................
http://www.ece.uvic.ca/~mdadams/sigsysbook
If you like these lecture slides, please consider posting a review of them at:
https://play.google.com/store/search?q=ISBN:9781550586787 or
http://books.google.com/books?vid=ISBN9781550586787
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 iii
License II
License who has not previously violated the terms of this License with
respect to the Work, or who has received express permission from the
Licensor to exercise rights under this License despite a previous
violation.
h. "Publicly Perform" means to perform public recitations of the Work and
to communicate to the public those public recitations, by any means or
process, including by wire or wireless means or public digital
performances; to make available to the public Works in such a way that
members of the public may access these Works from a place and at a
place individually chosen by them; to perform the Work to the public
by any means or process and the communication to the public of the
performances of the Work, including by public digital performance; to
broadcast and rebroadcast the Work by any means including signs,
sounds or images.
i. "Reproduce" means to make copies of the Work by any means including
without limitation by sound or visual recordings and the right of
fixation and reproducing fixations of the Work, including storage of a
protected performance or phonogram in digital form or other electronic
medium.
2. Fair Dealing Rights. Nothing in this License is intended to reduce,
limit, or restrict any uses free from copyright or rights arising from
limitations or exceptions that are provided for in connection with the
copyright protection under copyright law or other applicable laws.
3. License Grant. Subject to the terms and conditions of this License,
Licensor hereby grants You a worldwide, royalty-free, non-exclusive,
perpetual (for the duration of the applicable copyright) license to
exercise the rights in the Work as stated below:
a. to Reproduce the Work, to incorporate the Work into one or more
Collections, and to Reproduce the Work as incorporated in the
Collections; and,
b. to Distribute and Publicly Perform the Work including as incorporated
in Collections.
The above rights may be exercised in all media and formats whether now
known or hereafter devised. The above rights include the right to make
such modifications as are technically necessary to exercise the rights in
other media and formats, but otherwise you have no rights to make
Adaptations. Subject to 8(f), all rights not expressly granted by Licensor
are hereby reserved, including but not limited to the rights set forth in
Section 4(d).
4. Restrictions. The license granted in Section 3 above is expressly made
subject to and limited by the following restrictions:
a. You may Distribute or Publicly Perform the Work only under the terms
of this License. You must include a copy of, or the Uniform Resource
Identifier (URI) for, this License with every copy of the Work You
Distribute or Publicly Perform. You may not offer or impose any terms
on the Work that restrict the terms of this License or the ability of
the recipient of the Work to exercise the rights granted to that
recipient under the terms of the License. You may not sublicense the
Work. You must keep intact all notices that refer to this License and
to the disclaimer of warranties with every copy of the Work You
Distribute or Publicly Perform. When You Distribute or Publicly
Perform the Work, You may not impose any effective technological
measures on the Work that restrict the ability of a recipient of the
Work from You to exercise the rights granted to that recipient under
the terms of the License. This Section 4(a) applies to the Work as
incorporated in a Collection, but this does not require the Collection
apart from the Work itself to be made subject to the terms of this
License. If You create a Collection, upon notice from any Licensor You
must, to the extent practicable, remove from the Collection any credit
as required by Section 4(c), as requested.
b. You may not exercise any of the rights granted to You in Section 3
above in any manner that is primarily intended for or directed toward
commercial advantage or private monetary compensation. The exchange of
the Work for other copyrighted works by means of digital file-sharing
or otherwise shall not be considered to be intended for or directed
toward commercial advantage or private monetary compensation, provided
there is no payment of any monetary compensation in connection with
the exchange of copyrighted works.
c. If You Distribute, or Publicly Perform the Work or Collections, You
must, unless a request has been made pursuant to Section 4(a), keep
intact all copyright notices for the Work and provide, reasonable to
the medium or means You are utilizing: (i) the name of the Original
Author (or pseudonym, if applicable) if supplied, and/or if the
Original Author and/or Licensor designate another party or parties
(e.g., a sponsor institute, publishing entity, journal) for
attribution ("Attribution Parties") in Licensor’s copyright notice,
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 vii
License VI
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 viii
License VII
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 xii
Part 0
Preface
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 xiii
About These Lecture Slides
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 xiv
Video Lectures
The author has prepared video lectures for some of the material covered
in this textbook.
All of the videos are hosted by YouTube and available through the author’s
YouTube channel:
2 https://www.youtube.com/iamcanadian1867
In a definition, the term being defined is often typeset in a font like this.
To emphasize particular words, the words are typeset in a font like this.
URLs are typeset like http://www.ece.uvic.ca/~mdadams.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 xvi
Part 1
Introduction
audio).
2 A signal with more than one independent variable is said to be
multi-dimensional (e.g., image).
Continuous or discrete independent variables:
2 A signal with continuous independent variables is said to be continuous
x(t) x(n)
3 3
2 2
1 1
t n
−30 −20 −10 0 10 20 30 −3 −2 −1 0 1 2 3
x0 y0
x1 y1
x2 System y2
.. .. .. ..
. . . .
xM yN
|{z} |{z}
Input Signals Output Signals
Number of inputs:
2 A system with one input is said to be single input (SI).
2 A system with more than one input is said to be multiple input (MI).
Number of outputs:
2 A system with one output is said to be single output (SO).
2 A system with more than one output is said to be multiple output (MO).
Discrete-Time Discrete-Time
Input Signal Signal Output
Continuous-Time Before After Continuous-Time
Signal Continuous-to- Processing Processing Discrete-to- Signal
Discrete-Time
Discrete-Time Continuous-Time
System
(C/D) Converter (D/C) Converter
Continuous-Time Continuous-Time
Input Signal Signal Output
Discrete-Time Before After Discrete-Time
Signal Discrete-to- Processing Processing Continuous-to- Signal
Continuous-Time
Continuous-Time Discrete-Time
System
(D/C) Converter (C/D) Converter
Estimate of
Message Transmitted Received Message
Signal Signal Signal Signal
Transmitter Channel Receiver
Reference
Input Error Output
+ Controller Plant
−
Sensor
Feedback
Signal
Preliminaries
For two integers a and b, we define the following notation for sets of
consecutive integers:
[a . . b] = {x ∈ Z : a ≤ x ≤ b},
[a . . b) = {x ∈ Z : a ≤ x < b},
(a . . b] = {x ∈ Z : a < x ≤ b}, and
(a . . b) = {x ∈ Z : a < x < b}.
In this notation, a and b indicate the endpoints of the range for the set,
and the type of brackets used (i.e., parenthesis versus square bracket)
indicates whether each endpoint is included in the set.
For example:
2 [0 . . 4] denotes the set of integers {0, 1, 2, 3, 4};
2 [0 . . 4) denotes the set of integers {0, 1, 2, 3}; and
[a, b] = {x ∈ R : a ≤ x ≤ b},
(a, b) = {x ∈ R : a < x < b},
[a, b) = {x ∈ R : a ≤ x < b}, and
(a, b] = {x ∈ R : a < x ≤ b}.
In this notation, a and b indicate the endpoints of the interval for the set,
and the type of brackets used (i.e., parenthesis versus square bracket)
indicate whether each endpoint is included in the set.
For example:
2 [0, 100] denotes the set of all real numbers from 0 to 100, including both 0
and 100;
2 (−π, π] denotes the set of all real numbers from −π to π, excluding −π but
including π; and
2 [−π, π) denotes the set of all real numbers from −π to π, including −π but
excluding π.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 17
Mappings
To express that two functions f and g are equal, we can write either:
1 f = g; or
2 f (t) = g(t) for all t .
Of the preceding two expressions, the first (i.e., f = g) is usually
preferable, as it is less verbose.
For the functions f and g and an operation ◦ that is defined pointwise for
functions (such as addition, subtraction, multiplication, and division), the
following relationship holds:
functions mapping R to R.
2 The system H maps a function to a function.
2 In particular, the domain and codomain are each F , which is a set of
functions.
2 The system H multiplies its input function x by a factor of 2 in order to
produce its output function Hx.
2 Note that Hx is a function, not a number.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 23
Remarks on Operator Notation for CT Systems
For a system operator H and a function x, Hx is the function produced as
the output of the system H when the input is the function x.
Brackets around the operand of an operator are often omitted when not
required for grouping.
For example, for an operator H, a function x, and a real number t , we
would normally prefer to write:
1 Hx instead of the equivalent expression H(x); and
Properties of Signals
2 2
1 1
t n
−3 −2 −1 1 2 3 −3 −2 −1 1 2 3
−1 −1
−2 −2
2 2
1 1
t n
−3 −2 −1 1 2 3 −3 −2 −1 1 2 3
−1 −1
−2 −2
4
... ... 3
t ··· 2 ···
−2T −T T 2T
1
n
−4 −3 −2 −1 0 1 2 3 4 5 6 7
The period of a periodic signal is not unique. That is, a signal that is
periodic with period T is also periodic with period kT , for every (strictly)
positive integer k.
x(t)
2T 2T
... ...
t
−2T −T T 2T
T T
Time shifting (also called translation) maps the input function x to the
output function y as given by
x(t)
t
−3 −2 −1 0 1 2 3
x(t − 1) x(t + 1)
3 3
2 2
1 1
t t
−3 −2 −1 0 1 2 3 −3 −2 −1 0 1 2 3
Time reversal (also known as reflection) maps the input function x to the
output function y as given by
y(t) = x(−t).
Geometrically, the output function y is a reflection of the input function x
about the (vertical) line t = 0.
x(t) x(−t)
3 3
2 2
1 1
t t
−3 −2 −1 0 1 2 3 −3 −2 −1 0 1 2 3
y(t) = x(at),
x(t)
t
−2 −1 0 1 2
1
x(2t) x 2t
1 1
t t
−2 −1 0 1 2 −2 −1 0 1 2
Time scaling maps the input function x to the output function y as given by
y(t) = x(at),
x(t) x(2t)
1 1
t t
−2 −1 0 1 2 −2 −1 0 1 2
1
x 2t x(−t)
1 1
t t
−2 −1 0 1 2 −2 −1 0 1 2
2 first, time scaling x by a, and then time shifting the result by b/a.
Note that the time shift is not by the same amount in both cases.
In particular, note that when time scaling is applied first followed by time
shifting, the time shift is by b/a, not b.
− 12 1
Given x as shown −1 1 2 3
t
−2 −1
2
1 3
2 2
t
below, find −1 −1
y(t) = x(2t − 1).
x(t)
1
1 1
− 12 1
t 2
t
−2 −1 1 2 −2 −1 1 3
2 2
−1 −1
1 1
t t
−3 −2 −1 1 2 3 −3 −2 −1 1 2 3
−1 −1
−2 −2
1 x(t) −2x(t)
2
2 2
1 1
t t
−3 −2 −1 1 2 3 −3 −2 −1 1 2 3
−1 −1
−2 −2
y(t) = x(t) + b,
1 1
t t
−3 −2 −1 1 2 3 −3 −2 −1 1 2 3
−1 −1
−2 −2
y(t) = ax(t) + b,
function by b; or
2 first amplitude shifting x by b/a, and then amplitude scaling the resulting
function by a.
Properties of Functions
Sums involving even and odd functions have the following properties:
2 The sum of two even functions is even.
2 The sum of two odd functions is odd.
2 The sum of an even function and odd function is neither even nor odd,
provided that neither of the functions is identically zero.
That is, the sum of functions with the same type of symmetry also has the
same type of symmetry.
Products involving even and odd functions have the following properties:
2 The product of two even functions is even.
2 The product of two odd functions is even.
2 The product of an even function and an odd function is odd.
That is, the product of functions with the same type of symmetry is even,
while the product of functions with opposite types of symmetry is odd.
···
t
t0
···
t
t0
x(t)
t
t0 t1
A function that is neither left sided nor right sided is said to be two sided.
An example of a two-sided function is shown below.
x(t)
···
···
Elementary Functions
A cos θ
x(t) = Aeλt ,
A A A
t t t
x(t) = Ae jωt ,
Thus, Re{x} and Im{x} are the same except for a time shift.
2π
Also, x is periodic with fundamental period T = |ω| and fundamental
frequency |ω|.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 63
Complex Sinusoidal Functions (Continued)
The graphs of Re{x} and Im{x} have the forms shown below.
|A| |A|
|A| cos θ
t t
|A| sin θ
− |A|
− |A|
−0.5 1
−1 2
−1
−1.5 3 t
−2 4
−2
2
Im
Re
−4 2
1 1.5
−3
−2 1
0.5
−1
ω = −2π
−0.5 1
−1 2
−1
−1.5 3 t
−2 4
−2
Thus, Re{x} and Im{x} are each the product of a real exponential and
real sinusoid.
One of three distinct modes of behavior is exhibited by x(t), depending on
the value of σ.
If σ = 0, Re{x} and Im{x} are real sinusoids.
If σ > 0, Re{x} and Im{x} are each the product of a real sinusoid and a
growing real exponential.
If σ < 0, Re{x} and Im{x} are each the product of a real sinusoid and a
decaying real exponential.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 66
General Complex Exponential Functions (Continued)
The three modes of behavior for Re{x} and Im{x} are illustrated below.
|A| eσt
|A| eσt
|A| eσt
t t t
Due to the manner in which u is used in practice, the actual value of u(0)
is unimportant. Sometimes values of 0 and 21 are also used for u(0).
A plot of this function is shown below.
u(t)
1 ···
··· −1 0 1
t
sgnt
1 ···
t
−3 −2 −1 1 2 3
··· −1
rectt
1
t
− 12 0 1
2
Functions and sequences that are one over some subset of their domain
and zero elsewhere appear very frequently in engineering (e.g., the
unit-step function and rectangular function).
Indicator function notation provides a concise way to denote such
functions and sequences.
The indicator function of a subset S of a set A, denoted χS , is defined as
(
1 if t ∈ S
χS (t) =
0 otherwise.
A rectangular pulse (defined on R) having an amplitude of 1, a leading
edge at a, and falling edge at b is χ[a,b] .
The unit-step function (defined on R) is χ[0,∞) .
The unit-rectangular pulse (defined on R) is χ[−1/2,1/2] .
The unit-step sequence (defined on Z) is χ[0..∞) .
trit
1
t
− 12 0 1
2
The floor function, denoted b·c, is a function that maps a real number x
to the largest integer not more than x.
In other words, the floor function rounds a real number to the nearest
integer in the direction of negative infinity.
For example,
1 1
− 2 = −1, 2 = 0, and b1c = 1.
The ceiling function, denoted d·e, is a function that maps a real number x
to the smallest integer not less than x.
In other words, the ceiling function rounds a real number to the nearest
integer in the direction of positive infinity.
For example,
1 1
− 2 = 0, 2 = 1, and d1e = 1.
t t
0 0 t0
t
− 2ε 0 ε
2
R∞
Clearly, for any choice of ε, −∞ gε (t)dt = 1.
The function δ can be obtained as the following limit:
δ(t) = lim gε (t).
ε→0
That is, δ can be viewed as a limiting case of a rectangular pulse where
the pulse width becomes infinitesimally small and the pulse height
becomes infinitely large in such a way that the integral of the resulting
function remains unity.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 78
Properties of the Unit-Impulse Function
δ(t − t0 )
x(t)
1
x(t0 )
t
t t0
t0
Time-Shifted Unit-Impulse
Function x
Function
x(t)δ(t − t0 )
x(t0 )
t
t0
Product
The idea from the previous slide can be extended to handle any function
that is defined in a piecewise manner (i.e., via an expression involving
multiple cases).
That is, by using unit-step functions, we can always collapse a formula
involving multiple cases into a single expression.
Often, simplifying a formula in this way can be quite beneficial.
y = Hx,
x → y.
Note that the symbols “→” and “=” have very different meanings.
The symbol “→” should be read as “produces” (not as “equals”).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 84
Block Diagram Representations
Often, a system defined by the operator H and having the input x and
output y is represented in the form of a block diagram as shown below.
Input Output
x System y
H
y = H2 H1 x.
A parallel connection ties the inputs of both systems together and sums
their outputs.
The overall parallel-connected system is described by the equation
y = H1 x + H2 x.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 86
Section 3.5
t
−∞ t0 ∞
A system H is said to be causal if, for every real constant t0 , Hx(t0 ) does
not depend on x(t) for some t > t0 .
In other words, a causal system is such that the value of its output at any
given point in time can depend on the value of its input at only the same or
earlier points in time (i.e., not later points in time).
If the independent variable t represents time, a system must be causal in
order to be physically realizable.
Noncausal systems can sometimes be useful in practice, however, since
the independent variable need not always represent time (e.g., the
independent variable might represent position).
A memoryless system is always causal, although the converse is not
necessarily true.
t ≤ t0
t
−∞ t0 ∞
A system H−1 being the inverse of H means that the following two
systems are equivalent (i.e., H−1 H is an identity):
x y x y
H H−1
A system H is said to be time invariant (TI) (or shift invariant (SI)) if,
for every function x and every real constant t0 , the following condition
holds:
Hx(t − t0 ) = Hx0 (t) for all t, where x0 (t) = x(t − t0 )
(i.e., H commutes with time shifts).
In other words, a system is time invariant if a time shift (i.e., advance or
delay) in the input always results only in an identical time shift in the
output.
A system that is not time invariant is said to be time varying.
In simple terms, a time invariant system is a system whose behavior does
not change with respect to time.
Practically speaking, compared to time-varying systems, time-invariant
systems are much easier to design and analyze, since their behavior
does not change with respect to time.
Let St0 denote an operator that applies a time shift of t0 to a function (i.e.,
St0 x(t) = x(t − t0 )).
A system H is time invariant if and only if the following two systems are
equivalent (i.e., H commutes with St0 ):
x y
St0 H x y
H St0
The system H is additive if and only if the following two systems are
equivalent (i.e., H commutes with addition):
x1 y x1 y
+ H H +
x2 x2
H
System 1: y = H(x1 + x2 )
System 2: y = Hx1 + Hx2
x y x y
a H H a
The system H is linear if and only if the following two systems are
equivalent (i.e., H commutes with linear combinations):
x1 y x1 y
a1 + H H a1 +
x2 x2
a2 H a2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 100
Part 4
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 101
Why Linear Time-Invariant (LTI) Systems?
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 102
Section 4.1
Convolution
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 103
CT Convolution
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 104
Practical Convolution Computation
we proceed as follows:
1 Plot x(τ) and h(t − τ) as a function of τ.
2 Initially, consider an arbitrarily large negative value for t . This will result in
h(t − τ) being shifted very far to the left on the time axis.
3 Write the mathematical expression for x ∗ h(t).
4 Increase t gradually until the expression for x ∗ h(t) changes form. Record
the interval over which the expression for x ∗ h(t) was valid.
5 Repeat steps 3 and 4 until t is an arbitrarily large positive value. This
corresponds to h(t − τ) being shifted very far to the right on the time axis.
6 The results for the various intervals can be combined in order to obtain an
expression for x ∗ h(t) for all t .
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 105
Properties of Convolution
The convolution operation is commutative. That is, for any two functions x
and h,
x ∗ h = h ∗ x.
The convolution operation is associative. That is, for any functions x, h1 ,
and h2 ,
(x ∗ h1 ) ∗ h2 = x ∗ (h1 ∗ h2 ).
The convolution operation is distributive with respect to addition. That is,
for any functions x, h1 , and h2 ,
x ∗ (h1 + h2 ) = x ∗ h1 + x ∗ h2 .
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 106
Representation of Functions Using Impulses
x ∗ δ = x.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 107
Periodic Convolution
(i.e., x0 (t) equals x(t) over a single period of x and is zero elsewhere).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 108
Section 4.2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 109
Impulse Response
The response h of a system H to the input δ is called the impulse
response of the system (i.e., h = Hδ).
For any LTI system with input x, output y, and impulse response h, the
following relationship holds:
y = x ∗ h.
In other words, a LTI system simply computes a convolution.
Furthermore, a LTI system is completely characterized by its impulse
response.
That is, if the impulse response of a LTI system is known, we can
determine the response of the system to any input.
Since the impulse response of a LTI system is an extremely useful
quantity, we often want to determine this quantity in a practical setting.
Unfortunately, in practice, the impulse response of a system cannot be
determined directly from the definition of the impulse response.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 110
Step Response
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 111
Block Diagram Representation of LTI Systems
x y
h
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 112
Interconnection of LTI Systems
x y x y
h1 h2 ≡ h1 ∗ h2
x y
h1 +
x y
≡ h1 + h2
h2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 113
Section 4.3
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 114
Memory
A LTI system with impulse response h is memoryless if and only if
h(t) = Kδ(t),
y = x ∗ (Kδ) = Kx
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 115
Causality
a causal function.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 116
Invertibility
h ∗ hinv = δ.
Consequently, a LTI system with impulse response h is invertible if and
only if there exists a function hinv such that
h ∗ hinv = δ.
Except in simple cases, the above condition is often quite difficult to test.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 117
BIBO Stability
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 118
Eigenfunctions of LTI Systems
As it turns out, every complex exponential is an eigenfunction of all LTI
systems.
For a LTI system H with impulse response h,
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 120
Part 5
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 121
Introduction
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 122
Section 5.1
Fourier Series
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 123
Harmonically-Related Complex Sinusoids
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 124
CT Fourier Series
A periodic (complex-valued) function x with fundamental period T and
fundamental frequency ω0 = 2π T can be represented as a linear
combination of harmonically-related complex sinusoids as
∞
x(t) = ∑ ck e jkω0t .
k=−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 125
CT Fourier Series (Continued)
1
Z
ck = x(t)e− jkω0t dt,
T T
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 126
Section 5.2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 127
Remarks on Equality of Functions
The equality of functions can be defined in more than one way.
Two functions x and y are said to be equal in the pointwise sense if
x(t) = y(t) for all t (i.e., x and y are equal at every point).
Two functions x and y are said to be equal in the mean-squared error
(MSE) sense if |x(t) − y(t)|2 dt = 0 (i.e., the energy in x − y is zero).
R
Again, let xN denote the Fourier series for the periodic function x truncated
after the N th harmonic components as given by
N
xN (t) = ∑ ck e jkω0t .
k=−N
If limN→∞ xN (t) = x(t) for all t (i.e., limN→∞ xN is equal to x in the
pointwise sense), the Fourier series is said to converge pointwise to x.
If convergence is pointwise and the rate of convergence is the same
everywhere, the convergence is said to be uniform.
If limN→∞ 1 2
T T |xN (t) − x(t)| dt = 0 (i.e., limN→∞ xN is equal to x in the
R
MSE sense), the Fourier series is said to converge to x in the MSE sense.
Pointwise convergence is a stronger condition than MSE convergence
(i.e., pointwise convergence implies MSE convergence, but the converse
is not true).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 130
Convergence of Fourier Series: Continuous Case
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 131
Convergence of Fourier Series: Finite-Energy Case
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 132
Dirichlet Conditions
The Dirichlet conditions for the periodic function x are as follows:
T |x(t)| dt < ∞);
1 over a single period, x is absolutely integrable (i.e.,
R
2 over a single period, x has a finite number of maxima and minima (i.e., x is
which is finite.
Examples of functions violating the Dirichlet conditions are shown below.
x(t)
x(t)
t −1
t
1 2
x(t)
··· ···
1
2
1
4
··· ···
t
−1 0 1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 133
Convergence of Fourier Series: Dirichlet Case
where x(ta− ) and x(ta+ ) denote the values of the function x on the left- and
right-hand sides of the discontinuity, respectively.
Since most functions tend to satisfy the Dirichlet conditions and the above
convergence result specifies the value of the Fourier series at every point,
this result is often very useful in practice.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 134
Gibbs Phenomenon
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 135
Gibbs Phenomenon: Periodic Square Wave Example
1.5 1.5
1 1
0.5 0.5
0 0
-1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1
-0.5 -0.5
-1 -1
-1.5 -1.5
Fourier series truncated after the Fourier series truncated after the
3rd harmonic components 7th harmonic components
1.5 1.5
1 1
0.5 0.5
0 0
-1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1
-0.5 -0.5
-1 -1
-1.5 -1.5
Fourier series truncated after the Fourier series truncated after the
11th harmonic components 101st harmonic components
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 136
Section 5.3
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 137
Properties of (CT) Fourier Series
CTFS CTFS
x(t) ←→ ak and y(t) ←→ bk
Property
Parseval’s Relation 1R
T T |x(t)|2 dt = ∑∞
k=−∞ |ak |
2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 138
Linearity
CTFS
Let x and y be two periodic functions with the same period. If x(t) ←→ ak
CTFS
and y(t) ←→ bk , then
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 139
Time Shifting (Translation)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 140
Frequency Shifting (Modulation)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 141
Time Reversal (Reflection)
CTFS
x(−t) ←→ c−k .
That is, time reversal of a function results in a time reversal of its Fourier
series coefficients.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 142
Conjugation
x∗ (t) ←→ c∗−k
CTFS
In other words, conjugating a function has the effect of time reversing and
conjugating the Fourier series coefficient sequence.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 143
Periodic Convolution
x ~ y(t) ←→ Tak bk .
CTFS
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 144
Multiplication
CTFS
Let x and y be two periodic functions with the same period. If x(t) ←→ ak
CTFS
and y(t) ←→ bk , then
∞
∑
CTFS
x(t)y(t) ←→ an bk−n
n=−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 145
Parseval’s Relation
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 146
Even and Odd Symmetry
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 147
Real Functions
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 148
Trigonometric Forms of a Fourier Series
where θk = arg ck .
The trigonometric form of a Fourier series has the appearance
∞
x(t) = c0 + ∑ [αk cos(kω0t) + βk sin(kω0t)] ,
k=1
where αk = 2 Re ck and βk = −2 Im ck .
Note that the trigonometric forms contain only real quantities.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 149
Other Properties of Fourier Series
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 150
Section 5.4
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 151
A New Perspective on Functions: The Frequency Domain
The Fourier series provides us with an entirely new way to view functions.
Instead of viewing a function as having information distributed with respect
to time (i.e., a function whose domain is time), we view a function as
having information distributed with respect to frequency (i.e., a function
whose domain is frequency).
This so called frequency-domain perspective is of fundamental
importance in engineering.
Many engineering problems can be solved much more easily using the
frequency domain than the time domain.
The Fourier series coefficients of a function x provide a means to quantify
how much information x has at different frequencies.
The distribution of information in a function over different frequencies is
referred to as the frequency spectrum of the function.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 152
Motivating Example
Consider the real 1-periodic function x having the Fourier series
representation
2 j − j10πt
x(t) = − 10j e− j14πt − 10 e 4 j − j6πt
− 10 e − 13 j − j2πt
10 e
+ 13 j j2πt
10 e
4 j j6πt
+ 10 e 2 j j10πt
+ 10 e + 10j e j14πt .
A plot of x is shown below.
x(t)
1
1
2
t
−1 − 12 1 1
2
− 12
−1
The terms that make the most dominant contribution to the overall sum
are the ones with the largest magnitude coefficients.
To illustrate this, we consider the problem of determining the best
approximation of x that keeps only 4 of the 8 terms in the Fourier series.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 153
Motivating Example (Continued)
approximation
1 x(t)
1
2
t
−1 − 12 1 1
2
− 12
−1
approximation
1 x(t)
1
2
t
−1 − 12 1 1
2
− 12
−1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 154
Fourier Series and Frequency Spectra
To gain further insight into the role played by the Fourier series coefficients
ck in the context of the frequency spectrum of the function x, it is helpful to
write the Fourier series with the ck expressed in polar form as follows:
∞ ∞
x(t) = ∑ ck e jkω0t = ∑ |ck | e j(kω0t+arg ck ) .
k=−∞ k=−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 155
Fourier Series and Frequency Spectra (Continued)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 156
Frequency Spectra of Real Functions
Recall that, for a real function x, the Fourier series coefficient sequence c
satisfies
ck = c∗−k
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 158
Frequency Response
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 159
Fourier Series and LTI Systems
Thus, if the input x to a LTI system is a Fourier series, the output y is also
CTFS CTFS
a Fourier series. More specifically, if x(t) ←→ ck then y(t) ←→ H(kω0 )ck .
The above formula can be used to determine the output of a LTI system
from its input in a way that does not require convolution.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 160
Filtering
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 161
Ideal Lowpass Filter
An ideal lowpass filter eliminates all frequency components with a
frequency whose magnitude is greater than some cutoff frequency, while
leaving the remaining frequency components unaffected.
Such a filter has a frequency response of the form
(
1 |ω| ≤ ωc
H(ω) =
0 otherwise,
where ωc is the cutoff frequency.
A plot of this frequency response is given below.
H(ω)
ω
−ωc ωc
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 162
Ideal Highpass Filter
An ideal highpass filter eliminates all frequency components with a
frequency whose magnitude is less than some cutoff frequency, while
leaving the remaining frequency components unaffected.
Such a filter has a frequency response of the form
(
1 |ω| ≥ ωc
H(ω) =
0 otherwise,
where ωc is the cutoff frequency.
A plot of this frequency response is given below.
H(ω)
1
··· ···
ω
−ωc ωc
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 163
Ideal Bandpass Filter
An ideal bandpass filter eliminates all frequency components with a
frequency whose magnitude does not lie in a particular range, while
leaving the remaining frequency components unaffected.
Such a filter has a frequency response of the form
(
1 ωc1 ≤ |ω| ≤ ωc2
H(ω) =
0 otherwise,
ω
−ωc2 −ωc1 ωc1 ωc2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 164
Part 6
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 165
Motivation for the Fourier Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 166
Section 6.1
Fourier Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 167
Development of the Fourier Transform [Aperiodic Case]
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 168
Development of the Fourier Transform [Aperiodic Case] (Continued)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 169
Generalized Fourier Transform
The
R ∞ classical Fourier transform for aperiodic functions does not exist (i.e.,
x(t)e − jωt dt fails to converge) for some functions of great practical
−∞
interest, such as:
2 a nonzero constant function;
2 a periodic function (e.g., a real or complex sinusoid);
2 the unit-step function (i.e., u); and
2 the signum function (i.e., sgn).
Fortunately, the Fourier transform can be extended to handle such
functions, resulting in what is known as the generalized Fourier
transform.
For our purposes, we can think of the classical and generalized Fourier
transforms as being defined by the same formulas.
Therefore, in what follows, we will not typically make a distinction between
the classical and generalized Fourier transforms.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 170
CT Fourier Transform (CTFT)
The (CT) Fourier transform of the function x, denoted Fx or X , is given
by
Z ∞
Fx(ω) = X(ω) = x(t)e− jωt dt.
−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 172
Remarks on Dot Notation
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 173
Remarks on Notational Conventions
Since dot notation is less frequently used by engineers, the author has
elected to minimize its use herein.
To avoid ambiguous notation, the following conventions are followed:
1 in the expression for the operand of a Fourier transform operator, the
independent variable is assumed to be the variable named “t” unless
otherwise indicated (i.e., in terms of dot notation, each “t ” is treated as if it
were a “·”)
2 in the expression for the operand of the inverse Fourier transform operator,
the independent variable is assumed to be the variable named “ω” unless
otherwise indicated (i.e., in terms of dot notation, each “ω” is treated as if it
were a “·”).
For example, with these conventions:
2 “F{cos(t − τ)}” denotes the function that is the Fourier transform of the
of the function V (ω) = δ(3ω − λ) (not the inverse Fourier transform of the
function V (λ) = δ(3ω − λ)).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 174
Section 6.2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 175
Convergence of the Fourier Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 176
Convergence of the Fourier Transform: Continuous Case
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 177
Convergence of the Fourier Transform: Finite-Energy Case
R∞ 2
If a function x is of finite energy (i.e., −∞ |x(t)| dt < ∞), then its Fourier
transform representation converges in the MSE sense.
In other words, if x is of finite energy, then the energy E in the difference
function x̃ − x is zero; that is,
Z ∞
E= |x̃(t) − x(t)|2 dt = 0.
−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 178
Dirichlet Conditions
The Dirichlet conditions for the function x are as follows:
R∞
−∞ |x(t)| dt < ∞);
1 the function x is absolutely integrable (i.e.,
2 on any finite interval, x has a finite number of maxima and minima (i.e., x is
10 ··· t
−1.25 −1 −0.75 −0.5 −0.25 0.25 0.5 0.75 1 1.25
···t −1
0.5 1
0.75
0.5
0.25
···
t
0.25 0.5 0.75 1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 179
Convergence of the Fourier Transform: Dirichlet Case
1
+
x̃(ta ) = 2 x(ta ) + x(ta− ) ,
where x(ta− ) and x(ta+ ) denote the values of the function x on the left- and
right-hand sides of the discontinuity, respectively.
Since most functions tend to satisfy the Dirichlet conditions and the above
convergence result specifies the value of the Fourier transform
representation at every point, this result is often very useful in practice.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 180
Section 6.3
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 181
Properties of the (CT) Fourier Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 182
Properties of the (CT) Fourier Transform (Continued)
Property
R∞ 2 1 R∞ 2
Parseval’s Relation −∞ |x(t)| dt = 2π −∞ |X(ω)| dω
Even Symmetry x is even ⇔ X is even
Odd Symmetry x is odd ⇔ X is odd
Real / Conjugate Symmetry x is real ⇔ X is conjugate symmetric
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 183
(CT) Fourier Transform Pairs
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 184
Linearity
CTFT CTFT
If x1 (t) ←→ X1 (ω) and x2 (t) ←→ X2 (ω), then
CTFT
a1 x1 (t) + a2 x2 (t) ←→ a1 X1 (ω) + a2 X2 (ω),
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 185
Time-Domain Shifting (Translation)
CTFT
If x(t) ←→ X(ω), then
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 186
Frequency-Domain Shifting (Modulation)
CTFT
If x(t) ←→ X(ω), then
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 187
Time- and Frequency-Domain Scaling (Dilation)
CTFT
If x(t) ←→ X(ω), then
CTFT 1 ω
x(at) ←→ X ,
|a| a
where a is an arbitrary nonzero real constant.
This is known as the dilation (or time/frequency-domain scaling)
property of the Fourier transform.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 188
Conjugation
CTFT
If x(t) ←→ X(ω), then
x∗ (t) ←→ X ∗ (−ω).
CTFT
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 189
Duality
CTFT
If x(t) ←→ X(ω), then
CTFT
X(t) ←→ 2πx(−ω)
This is known as the duality property of the Fourier transform.
This property follows from the high degree of symmetry in the forward and
inverse Fourier transform equations, which are respectively given by
Z ∞ Z ∞
X(λ) = x(θ)e− jθλ dθ and x(λ) = 2π 1
X(θ)e jθλ dθ.
−∞ −∞
That is, the forward and inverse Fourier transform equations are identical
except for a factor of 2π and different sign in the parameter for the
exponential function.
CTFT
Although the relationship x(t) ←→ X(ω) only directly provides us with the
Fourier transform of x(t), the duality property allows us to indirectly infer
the Fourier transform of X(t). Consequently, the duality property can be
used to effectively double the number of Fourier transform pairs that we
know.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 190
Time-Domain Convolution
CTFT CTFT
If x1 (t) ←→ X1 (ω) and x2 (t) ←→ X2 (ω), then
CTFT
x1 ∗ x2 (t) ←→ X1 (ω)X2 (ω).
This is known as the convolution (or time-domain convolution)
property of the Fourier transform.
In other words, a convolution in the time domain becomes a multiplication
in the frequency domain.
This suggests that the Fourier transform can be used to avoid having to
deal with convolution operations.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 191
Time-Domain Multiplication
CTFT CTFT
If x1 (t) ←→ X1 (ω) and x2 (t) ←→ X2 (ω), then
Z ∞
1 1
X1 (θ)X2 (ω − θ)dθ.
CTFT
x1 (t)x2 (t) ←→ 2π X1 ∗ X2 (ω) = 2π
−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 192
Time-Domain Differentiation
CTFT
If x(t) ←→ X(ω), then
dx(t) CTFT
←→ jωX(ω).
dt
This is known as the (time-domain) differentiation property of the
Fourier transform.
Differentiation in the time domain becomes multiplication by jω in the
frequency domain.
Of course,
by repeated application of the above property, we have that
d n
( jω)n X(ω).
CTFT
dt x(t) ←→
The above suggests that the Fourier transform might be a useful tool
when working with differential (or integro-differential) equations.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 193
Frequency-Domain Differentiation
CTFT
If x(t) ←→ X(ω), then
CTFT d
tx(t) ←→ j X(ω).
dω
This is known as the frequency-domain differentiation property of the
Fourier transform.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 194
Time-Domain Integration
CTFT
If x(t) ←→ X(ω), then
Z t
1
X(ω) + πX(0)δ(ω).
CTFT
x(τ)dτ ←→
−∞ jω
This is known as the (time-domain) integration property of the Fourier
transform.
Whereas differentiation in the time domain corresponds to multiplication
by jω in the frequency domain, integration in the time domain is
associated with division by jω in the frequency domain.
Since integration in the time domain becomes division by jω in the
frequency domain, integration can be easier to handle in the frequency
domain.
The above property suggests that the Fourier transform might be a useful
tool when working with integral (or integro-differential) equations.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 195
Parseval’s Relation
R∞ 2
Recall that the energy of a function x is given by −∞ |x(t)| dt .
CTFT
If x(t) ←→ X(ω), then
Z ∞ Z ∞
|x(t)|2 dt = 1
2π |X(ω)|2 dω
−∞ −∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 196
Even/Odd Symmetry
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 197
Real Functions
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 198
More Fourier Transforms
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 199
Section 6.4
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 200
Fourier Transform of Periodic Functions
The Fourier transform can be generalized to also handle periodic
functions.
Consider a periodic function x with period T and frequency ω0 = 2π
T .
Define the function xT as
(
x(t) − T2 ≤ t < T2
xT (t) =
0 otherwise.
(i.e., xT (t) is equal to x(t) over a single period and zero elsewhere).
Let a denote the Fourier series coefficient sequence of x.
Let X and XT denote the Fourier transforms of x and xT , respectively.
The following relationships can be shown to hold:
∞
X(ω) = ∑ ω0 XT (kω0 )δ(ω − kω0 ),
k=−∞
∞
ak = T1 XT (kω0 ), and X(ω) = ∑ 2πak δ(ω − kω0 ).
k=−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 201
Fourier Transform of Periodic Functions (Continued)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 202
Section 6.5
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 203
The Frequency-Domain Perspective on Functions
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 204
Fourier Transform and Frequency Spectra
To gain further insight into the role played by the Fourier transform X in
the context of the frequency spectrum of x, it is helpful to write the Fourier
transform representation of x with X(ω) expressed in polar form as
follows:
Z ∞ Z ∞
x(t) = 1
2π X(ω)e jωt dω = 1
2π |X(ω)| e j[ωt+arg X(ω)] dω.
−∞ −∞
In effect, the quantity |X(ω)| is a weight that determines how much the
complex sinusoid at frequency ω contributes to the integration result x.
The quantity arg X(ω) determines how the complex sinusoid at frequency
ω is shifted related to complex sinusoids at other frequencies.
Perhaps, this can be more easily seen if we express the above integral as
the limit of a sum, derived from an approximation of the integral using the
areas of rectangles, as shown on the next slide. [Recall that
R∞ ∞
−∞ f (x)dx = lim∆x→0 ∑k=−∞ ∆x f (k∆x).]
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 205
Fourier Transform and Frequency Spectra (Continued 1)
Expressing the integral (from the previous slide) as the limit of a sum, we
obtain
∞
x(t) = lim 1
∆ω→0 2π
∑ ∆ω |X(ω)| e j[ωt+arg X(ω)] ,
k=−∞
where ω = k∆ω.
In the above equation, the kth term in the summation corresponds to a
complex sinusoid with fundamental frequency ω = k∆ω that has had its
amplitude scaled by a factor of |X(ω)| and has been time shifted by an
amount that depends on arg X(ω).
For a given ω = k∆ω (which is associated with the kth term in the
summation), the larger |X(ω)| is, the larger the amplitude of its
corresponding complex sinusoid e jωt will be, and therefore the larger the
contribution the kth term will make to the overall summation.
In this way, we can use |X(ω)| as a measure of how much information a
function x has at the frequency ω.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 206
Fourier Transform and Frequency Spectra (Continued 2)
The Fourier transform X of the function x is referred to as the frequency
spectrum of x.
The magnitude |X(ω)| of the Fourier transform X is referred to as the
magnitude spectrum of x.
The argument arg X(ω) of the Fourier transform X is referred to as the
phase spectrum of x.
Since the Fourier transform is a function of a real variable, a function can
potentially have information at any real frequency.
Since the Fourier transform X of a periodic function x with fundamental
frequency ω0 and the Fourier series coefficient sequence a is given by
X(ω) = ∑∞ k=−∞ 2πak δ(ω − kω0 ), the Fourier transform and Fourier series
give consistent results for the frequency spectrum of a periodic function.
Since the frequency spectrum is complex (in the general case), it is
usually represented using two plots, one showing the magnitude
spectrum and one showing the phase spectrum.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 207
Frequency Spectra of Real Functions
Recall that, for a real function x, the Fourier transform X of x satisfies
X(ω) = X ∗ (−ω)
1 1
ω ω
−B B − B2 B
2
One can show that a function cannot be both time limited and
bandlimited.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 209
Energy-Density Spectra
By Parseval’s relation, the energy E in a function x with Fourier transform
X is given by
Z ∞
1
E= 2π Ex (ω)dω,
−∞
where
Ex (ω) = |X(ω)|2 .
We refer to Ex as the energy-density spectrum of the function x.
The function Ex indicates how the energy in x is distributed with respect to
frequency.
For example, the energy contributed by frequencies in the range [ω1 , ω2 ]
is given by
Z ω2
1
2π Ex (ω)dω.
ω1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 210
Section 6.6
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 211
Frequency Response of LTI Systems
Consider a LTI system with input x, output y, and impulse response h, and
let X , Y , and H denote the Fourier transforms of x, y, and h, respectively.
Since y(t) = x ∗ h(t), we have that
Y (ω) = X(ω)H(ω).
The function H is called the frequency response of the system.
A LTI system is completely characterized by its frequency response H .
The above equation provides an alternative way of viewing the behavior of
a LTI system. That is, we can view the system as operating in the
frequency domain on the Fourier transforms of the input and output
functions.
The frequency spectrum of the output is the product of the frequency
spectrum of the input and the frequency response of the system.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 212
Frequency Response of LTI Systems (Continued 1)
|Y (ω)| = |X(ω)| |H(ω)| and argY (ω) = arg X(ω) + arg H(ω).
The magnitude spectrum of the output equals the magnitude spectrum of
the input times the magnitude response of the system.
The phase spectrum of the output equals the phase spectrum of the input
plus the phase response of the system.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 213
Frequency Response of LTI Systems (Continued 2)
(i.e., the magnitude response |H(ω)| is even and the phase response
arg H(ω) is odd).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 214
Unwrapped Phase
For many types of analysis, restricting the range of a phase function to an
interval of length 2π (such as (−π, π]), often unnecessarily introduces
discontinuities into the function.
This motivates the notion of unwrapped phase.
The unwrapped phase is simply the phase defined in such a way so as
not to restrict the phase to an interval of length 2π and to keep the phase
function continuous to the greatest extent possible.
For example, the function H(ω) = e jπω has the unwrapped phase
Θ(ω) = πω.
Arg H(ω) Θ(ω) = πω
3π 3π
2π 2π
π π
ω ω
−4 −2 2 4 −4 −2 2 4
−π −π
−2π −2π
−3π −3π
2 translating by −
arg H(ω)
ω .
Therefore, the magnitude response determines how different complex
sinusoids are scaled (in amplitude) by the system.
Similarly, the phase response determines how different complex sinusoids
are translated (i.e., delayed/advanced) by the system.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 216
Magnitude Distortion
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 217
Phase Distortion
Recall that a LTI system H with frequency response H is such that
H{e jωt }(t) = |H(ω)| e jω[t−τp (ω)] where τp (ω) = − arg H(ω)
ω .
The function τp is known as the phase delay of the system.
If τp (ω) = td (where td is a constant), the system shifts all complex
sinusoids by the same amount td .
Since τp (ω) = td is equivalent to the (unwrapped) phase response being
of the form arg H(ω) = −td ω (which is a linear function with a zero
constant term), a system with a constant phase delay is said to have
linear phase.
In the case that τp (ω) = 0, the system is said to have zero phase.
If τp (ω) is not a constant, different complex sinusoids are shifted by
different amounts, resulting in what is known as phase distortion.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 218
Distortionless Transmission
Consider a LTI system H with input x and output y given by
y(t) = x(t − t0 ),
where t0 is a real constant.
That is, the output of the system is simply the input delayed by t0 .
This type of behavior is the ideal for which we strive in real-world
communication systems (i.e., the received signal y equals a delayed
version of the transmitted signal x).
Taking the Fourier transform of the preceding equation, we have
H(ω) = e− jωt0 .
Since the phase delay of the system is τp (ω) = − −ωt
ω
0
= t0 , the phase
delay is constant and the system has linear phase.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 219
Magnitude and Phase Distortion in Audio
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 220
Magnitude and Phase Distortion in Images
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 221
Example: Magnitude and Phase Distortion in Images (1)
Image A Image B
Magnitude Spectrum from Image B and Magnitude Spectrum from Image A and
Phase Spectrum from Image A Phase Spectrum from Image B
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 222
Example: Magnitude and Phase Distortion in Images (2)
Magnitude Spectrum from Image B and Magnitude Spectrum from Image A and
Phase Spectrum from Image A Phase Spectrum from Image B
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 223
Block Diagram Representations of LTI Systems
Consider a LTI system with input x, output y, and impulse response h, and
let X , Y , and H denote the Fourier transforms of x, y, and h, respectively.
Often, it is convenient to represent such a system in block diagram form in
the frequency domain as shown below.
X Y
H
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 224
Interconnection of LTI Systems
X Y
≡ H1 + H2
H2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 225
LTI Systems and Differential Equations
Many LTI systems of practical interest can be represented using an
Nth-order linear differential equation with constant coefficients.
Consider a system with input x and output y that is characterized by an
equation of the form
N M
d k d k
∑ bk dt y(t) = ∑ ak dt x(t),
k=0 k=0
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 226
Section 6.7
Application: Filtering
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 227
Filtering
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 228
Ideal Lowpass Filter
An ideal lowpass filter eliminates all frequency components with a
frequency whose magnitude is greater than some cutoff frequency, while
leaving the remaining frequency components unaffected.
Such a filter has a frequency response H of the form
(
1 |ω| ≤ ωc
H(ω) =
0 otherwise,
where ωc is the cutoff frequency.
A plot of this frequency response is given below.
H(ω)
ω
−ωc ωc
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 229
Ideal Highpass Filter
An ideal highpass filter eliminates all frequency components with a
frequency whose magnitude is less than some cutoff frequency, while
leaving the remaining frequency components unaffected.
Such a filter has a frequency response H of the form
(
1 |ω| ≥ ωc
H(ω) =
0 otherwise,
where ωc is the cutoff frequency.
A plot of this frequency response is given below.
H(ω)
1
··· ···
ω
−ωc ωc
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 230
Ideal Bandpass Filter
An ideal bandpass filter eliminates all frequency components with a
frequency whose magnitude does not lie in a particular range, while
leaving the remaining frequency components unaffected.
Such a filter has a frequency response H of the form
(
1 ωc1 ≤ |ω| ≤ ωc2
H(ω) =
0 otherwise,
where the limits of the passband are ωc1 and ωc2 .
A plot of this frequency response is given below.
H(ω)
ω
−ωc2 −ωc1 ωc1 ωc2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 231
Section 6.8
Application: Equalization
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 232
Equalization
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 233
Equalization (Continued)
Input Output x y
horig heq horig
Let Horig denote the frequency response of original system (i.e., without
equalization).
Let Hd denote the desired frequency response.
Let Heq denote the frequency response of the equalizer.
The new system with equalization has frequency response
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 235
Electronic Circuits
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 236
Resistors
R
i + −
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 237
Inductors
An inductor is a circuit element that converts an electric current into a
magnetic field and vice versa.
An inductor uses the energy stored in a magnetic field in order to oppose
changes in current (through the inductor).
An inductor is characterized by an equation of the form
Z t
d 1
v(t) = L dt i(t) (or equivalently, i(t) = L v(τ)dτ),
−∞
where L is a nonnegative real constant, and v and i respectively denote
the voltage across and current through the inductor as a function of time.
As a matter of terminology, the quantity L is known as the inductance of
the inductor.
Inductance is measured in units of henrys (H).
In circuit diagrams, an inductor is denoted by the symbol shown below.
L
i + −
v
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 238
Capacitors
A capacitor is a circuit element that stores electric charge.
A capacitor uses the energy stored in an electric field in order to oppose
changes in voltage (across the capacitor).
A capacitor is characterized by an equation of the form
Z t
1
v(t) = C i(τ)dτ (or equivalently, i(t) = C dtd v(t)),
−∞
where C is a nonnegative real constant, and v and i respectively denote
the voltage across and current through the capacitor as a function of time.
As a matter of terminology, the quantity C is known as the capacitance of
the capacitor.
Capacitance is measured in units of farads (F).
In circuit diagrams, a capacitor is denoted by the symbol shown below.
C
i + −
v
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 239
Circuit Analysis with the Fourier Transform
The Fourier transform is a very useful tool for circuit analysis.
The utility of the Fourier transform is partly due to the fact that the
differential/integral equations that describe inductors and capacitors are
much simpler to express in the Fourier domain than in the time domain.
Let v and i denote the voltage across and current through a circuit
element, and let V and I denote the Fourier transforms of v and i,
respectively.
In the frequency domain, the equations characterizing a resistor, an
inductor, and a capacitor respectively become:
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 240
Section 6.10
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 241
Motivation for Amplitude Modulation (AM)
In communication systems, we often need to transmit a signal using a
frequency range that is different from that of the original signal.
For example, voice/audio signals typically have information in the range of
0 to 22 kHz.
Often, it is not practical to transmit such a signal using its original
frequency range.
Two potential problems with such an approach are:
1 interference; and
2 constraints on antenna length.
Since many signals are broadcast over the airwaves, we need to ensure
that no two transmitters use the same frequency bands in order to avoid
interference.
Also, in the case of transmission via electromagnetic waves (e.g., radio
waves), the length of antenna required becomes impractically large for the
transmission of relatively low frequency signals.
For the preceding reasons, we often need to change the frequency range
associated with a signal before transmission.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 242
Trivial Amplitude Modulation (AM) System
x y y x̂
× ×
Transmitter Receiver
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 243
Trivial Amplitude Modulation (AM) System: Example
X(ω)
C1 (ω) C2 (ω)
1
2π 2π
ω
−ωb ωb
ω ω
ωc −ωc
Transmitter Input
Y (ω) X̂(ω)
1 1
ω ω
ωb ωc − ωb ωc ωc + ωb −ωb ωb
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 244
Double-Sideband Suppressed-Carrier (DSB-SC) AM
Transmitter Receiver
If ωb < ωc0 < 2ωc − ωb , we have X̂(ω) = X(ω) (implying x̂(t) = x(t)).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 245
DSB-SC AM: Transmitter
c(t) = cos(ωct)
x y
×
y(t) = cos(ωct)x(t)
X = Fx, Y = Fy
Y (ω) = F{cos(ωct)x(t)}(ω)
= F 21 e jωc t + e− jωc t x(t) (ω)
= 21 F{e jωc t x(t)}(ω) + F{e− jωc t x(t)}(ω)
= 12 [X(ω − ωc ) + X(ω + ωc )]
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 246
DSB-SC AM: Receiver
c(t) = cos(ωct)
2ωc0
h(t) = π sinc(ωc0t)
y v x̂
× h
2ωc0
v(t) = cos(ωct)y(t), h(t) = π sinc(ωc0t), x̂(t) = v ∗ h(t)
Y = Fy, V = Fv, H = Fh, X̂ = Fx̂
V (ω) = F{cos(ωct)y(t)}(ω)
= F 21 e jωc t + e− jωc t y(t) (ω)
= 12 F e jωc t y(t) (ω) + F e− jωc t y(t) (ω)
= 21 [Y (ω − ωc ) +Y (ω + ωc )]
n o
H(ω) = F 2ωπc0 sinc(ωc0t) (ω)
= 2 rect 2ωωc0
V (ω) = 12 [Y (ω − ωc ) +Y (ω + ωc )]
= 12 21 [X([ω − ωc ] − ωc ) + X([ω − ωc ] + ωc )] +
1
2 [X([ω + ωc ] − ωc ) + X([ω + ωc ] + ωc )]
= 12 X(ω) + 14 X(ω − 2ωc ) + 14 X(ω + 2ωc )
X̂(ω) = H(ω)V (ω)
= H(ω) 12 X(ω) + 14 X(ω − 2ωc ) + 14 X(ω + 2ωc )
= 12 H(ω)X(ω) + 14 H(ω)X(ω − 2ωc ) + 14 H(ω)X(ω + 2ωc )
= 21 [2X(ω)] + 14 (0) + 14 (0)
= X(ω)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 248
DSB-SC AM: Example
X(ω)
1 C(ω) H(ω)
π π 2
−ωb ωb ω
−ωc ωc ω −ωc0 ωc0 ω
Transmitter Input
Y (ω)
1
2
ω
−2ωc −ωc − ωb−ωc−ωc + ωb −ωb ωb ωc − ωb ωc ωc + ωb 2ωc
Transmitter Output
V (ω)
1
2
1
4
ω
−2ωc − ω−2ωc c + ωb −ωb ωb 2ωc − ωb 2ωc 2ωc + ωb
b −2ω
X̂(ω)
ω
−2ωc −ωc −ωb ωb ωc 2ωc
Receiver Output
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 249
Single-Sideband Suppressed-Carrier (SSB-SC) AM
Transmitter Receiver
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 250
SSB-SC AM: Example
X(ω) C(ω) G(ω) H(ω)
1
π π ··· 1 ··· 4
ω ωc ω ωc ω ωc0 ω
−ωb ωb −ωc −ωc −ωc0
Q(ω)
1
2
−ωc − ω−ω ω
−2ωc c c + ωb −ωb
b −ω ωb ωc − ωb ωc ωc + ωb 2ωc
Y (ω)
1
2
−ωc − ω−ω ω
−2ωc c c + ωb −ωb
b −ω ωb ωc − ωb ωc ωc + ωb 2ωc
V (ω)
1
2
1
4
ωb ω
ωb −2ω
−2ωc − −2ω c c + ωb −ωb 2ωc − ωb2ωc2ωc + ωb
X̂(ω)
1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 251
Section 6.11
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 252
Sampling and Interpolation
Note that, unless very special conditions are met, the sampling process
loses information (i.e., is not invertible).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 253
Periodic Sampling
Although sampling can be performed in many different ways, the most
commonly used scheme is periodic sampling.
With this scheme, a sequence y of samples is obtained from a function x
according to the relation
y(n) = x(T n) for all integer n,
where T is a (strictly) positive real constant.
As a matter of terminology, we refer to T as the sampling period, and
ωs = 2π
T as the (angular) sampling frequency.
An example of periodic sampling is shown below, where the function x has
been sampled with sampling period T = 10, yielding the sequence y.
x(t) y(n)
4 4
3 3
2 2
1 1
t n
0 10 20 30 40 50 60 70 0 1 2 3 4 5 6 7
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 254
Invertibility of Sampling
Unless constraints are placed on the functions being sampled, the
sampling process is not invertible.
In other words, in the absence of any constraints, a function cannot be
uniquely determined from a sequence of its equally-spaced samples.
Consider, for example, the functions x1 and x2 given by
x s convert from y
× impulse train
to sequence
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 256
Model of Sampling: Various Signals
x(t)
p(t)
4
3 1 1 1 1
2 ··· ···
1 t
0 T 2T 3T
t
0 T 2T 3T Periodic Impulse Train
Input Function
s(t)
y(n)
x(T )
4 x(T )
4
3 x(3T )
x(2T )
3 x(2T ) x(3T )
2 x(0)
2 x(0)
1
1
t
0 T 2T 3T n
0 1 2 3
Impulse-Sampled Function
Output Sequence (Discrete-Time)
(Continuous-Time)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 257
Model of Sampling: Invertibility of Sampling Revisited
ideal C/D converter
∞
p(t) = ∑ δ(t − kT )
k=−∞
x s convert from y
× impulse train
to sequence
Since sampling is not invertible and our model of sampling consists of only
two steps, at least one of these two steps must not be invertible.
Recall the two steps in our model of sampling are as follows (in order):
∞
1 x −→ s(t) = x(t)p(t) = ∑ x(nT )δ(t − nT ); and
n=−∞
∞
2 s(t) = ∑ x(nT )δ(t − nT ) −→ y(n) = x(nT ).
n=−∞
Step 1 cannot be undone (unless we somehow restrict which functions x
can be sampled).
Step 2 is always invertible.
Therefore, the fact that sampling is not invertible is entirely due to step 1.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 258
Model of Sampling: Characterization
ideal C/D converter
∞
p(t) = ∑ δ(t − kT )
k=−∞
x s convert from y
× impulse train
to sequence
∞
p(t) = ∑ δ(t − kT ), ωs = 2π
T
k=−∞
∞
p(t) = ∑ ck e jkωs t
k=−∞
Z T /2
ck = 1
T p(t)e− jkωs t dt
−T /2
Z T /2
= 1
T δ(t)e− jkωs t dt
−T /2
Z ∞
= 1
T δ(t)e− jkωs t dt
−∞
1
= T
ωs
= 2π
∞
p(t) = ωs
2π ∑ e jkωs t
k=−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 260
Sampling: Multiplication by a Periodic Impulse Train
ideal C/D converter
∞
p(t) = ∑ δ(t − kT )
k=−∞
x s convert from y
× impulse train
to sequence
∞
s(t) = p(t)x(t), p(t) = ∑ δ(t − kT ), ωs = 2π
T
k=−∞
∞
p(t) = ωs
2π ∑ e jkωs t
k=−∞
∞
s(t) = ωs
2π ∑ e jkωs t x(t)
k=−∞
X = Fx, S = Fs
∞
S(ω) = ωs
2π ∑ X(ω − kωs )
k=−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 261
Model of Sampling: Aliasing
Consider frequency spectrum S of the impulse-sampled function s given by
∞
S(ω) = ωs
2π ∑ X(ω − kωs ).
k=−∞
Spectrum of Input
1
Function
(Bandwidth ωm )
ω
−ωm 0 ωm
S(ω)
Spectrum of Impulse-
1
T Sampled Function:
··· ··· No Aliasing Case
(ωs > 2ωm )
ω
−ωs − ωm −ωs −ωs + ωm −ωm 0 ωm ωs − ωm ωs ωs + ωm
S(ω)
Spectrum of Impulse-
1
T
Sampled Function:
··· ··· Aliasing Case
ω
(ωs ≤ 2ωm )
0 ωm
ωs − ωm ωs
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 263
Model of Interpolation
y convert from s x̂
sequence to h
impulse train
π
h(t) = sinc Tt
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 264
Sampling Theorem
Sampling Theorem. Let x be a function with Fourier transform X , and
suppose that |X(ω)| = 0 for all ω satisfying |ω| > ωM (i.e., x is
bandlimited to frequencies [−ωM , ωM ]). Then, x is uniquely determined by
its samples y(n) = x(T n) for all integer n, if
ωs > 2ωM ,
where ωs = 2πT . The preceding inequality is known as the Nyquist
condition. If this condition is satisfied, we have that
∞ π
x(t) = ∑ y(n) sinc T (t − T n) ,
n=−∞
We call ω2s the Nyquist frequency and 2ωM the Nyquist rate.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 265
Part 7
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 266
Motivation Behind the Laplace Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 267
Motivation Behind the Laplace Transform (Continued)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 268
Section 7.1
Laplace Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 269
(Bilateral) Laplace Transform
The (bilateral) Laplace transform of the function x, denoted Lx or X , is
defined as
Z ∞
Lx(s) = X(s) = x(t)e−st dt.
−∞
The inverse Laplace transform of X , denoted L−1 X or x, is then given
by
Z σ+ j∞
1
L−1 X(t) = x(t) = X(s)est ds,
2π j σ− j∞
where Re(s) = σ is in the ROC of X . (Note that this is a contour
integration, since s is complex.)
We refer to x and X as a Laplace transform pair and denote this
relationship as
LT
x(t) ←→ X(s).
In practice, we do not usually compute the inverse Laplace transform by
directly using the formula from above. Instead, we resort to other means
(to be discussed later).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 270
Bilateral and Unilateral Laplace Transforms
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 271
Remarks on Operator Notation
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 272
Remarks on Dot Notation
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 273
Remarks on Notational Conventions
Since dot notation is less frequently used by engineers, the author has
elected to minimize its use herein.
To avoid ambiguous notation, the following conventions are followed:
1 in the expression for the operand of a Laplace transform operator, the
independent variable is assumed to be the variable named “t” unless
otherwise indicated (i.e., in terms of dot notation, each “t ” is treated as if it
were a “·”)
2 in the expression for the operand of the inverse Laplace transform operator,
the independent variable is assumed to be the variable named “s” unless
otherwise indicated (i.e., in terms of dot notation, each “s” is treated as if it
were a “·”).
For example, with these conventions:
2 “L{(t − τ)u(t − τ)}” denotes the function that is the Laplace transform of
the function v(t) = (t − τ)u(t − τ) (not the Laplace transform of the function
v(τ) = (t − τ)u(t − τ)).
2 “L−1 { 1 }” denotes the function that is the inverse Laplace transform of
s2 −λ
the function V (s) = { s2 1−λ } (not the inverse Laplace transform of the
function V (λ) = { s2 1−λ }).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 274
Relationship Between Laplace and Fourier Transforms
Let X and XF denote the Laplace and (CT) Fourier transforms of x,
respectively.
The function X evaluated at jω (where ω is real) yields XF (ω). That is,
X( jω) = XF (ω).
Due to the preceding relationship, the Fourier transform of x is sometimes
written as X( jω).
The function X evaluated at an arbitrary complex value s = σ + jω (where
σ = Re(s) and ω = Im(s)) can also be expressed in terms of a Fourier
transform involving x. In particular, we have
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 276
Section 7.2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 277
Left-Half Plane (LHP)
Re(s) < a
Re{s} Re{s}
a a
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 278
Right-Half Plane (RHP)
Re(s) > a
Re{s} Re{s}
a a
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 279
Intersection of Sets
1 1 1
Re Re Re
−3 −2 −1 1 2 3 −3 −2 −1 1 2 3 −3 −2 −1 1 2 3
−1 −1 −1
−2 −2 −2
R1 R2 R1 ∩ R2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 280
Adding a Scalar to a Set
S + a = {z + a : z ∈ S}
2 2
1 1
Re Re
−3 −2 −1 1 2 3 −3 −2 −1 1 2 3
−1 −1
−2 −2
R R+1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 281
Multiplying a Set by a Scalar
For a set S and a scalar constant a, aS denotes the set given by
aS = {az : z ∈ S}
(i..e, aS is the set formed by multiplying each element of S by a).
Multiplying z by a affects z by: scaling by |a| and rotating about the origin
by arg a.
So, effectively, multiplying a set by a scalar applies a scaling and/or
rotation to the region associated with the set.
An illustrative example is given below.
Im
2
Re
−5 −4 −3 −2 −1 1 2 3 4 5
−1
−2
Im
R Im
2 2
1 1
Re Re
−5 −4 −3 −2 −1 1 2 3 4 5 −5 −4 −3 −2 −1 1 2 3 4 5
−1 −1
−2 −2
2R −2R
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 282
Region of Convergence (ROC)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 283
ROC Property 1: General Form
Re Re Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 284
ROC Property 2: Rational Laplace Transforms
Re Re Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 285
ROC Property 3: Finite-Duration Functions
Re Re Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 286
ROC Property 4: Right-Sided Functions
Re Re Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 287
ROC Property 5: Left-Sided Functions
If a function x is left sided and the (vertical) line Re(s) = σ0 is in the ROC
of the Laplace transform X = Lx, then all values of s for which Re(s) < σ0
must also be in the ROC (i.e., the ROC includes a LHP containing
Re(s) = σ0 ).
Thus, if x is left sided but not right sided, the ROC of X is a LHP.
Some examples of sets that would be either valid or invalid as ROCs for
X , if x is left sided but not right sided, are shown below.
Im Im Im
Re Re Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 288
ROC Property 6: Two-Sided Functions
If a function x is two sided and the (vertical) line Re(s) = σ0 is in the ROC
of the Laplace transform X = Lx, then the ROC will consist of a strip in
the complex plane that includes the line Re(s) = σ0 .
Some examples of sets that would be either valid or invalid as ROCs for
X , if x is two sided, are shown below.
Im Im Im
Re Re Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 289
ROC Property 7: More on Rational Laplace Transforms
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 290
General Form of the ROC
To summarize the results of properties 3, 4, 5, and 6, if the Laplace
transform X of the function x exists, the ROC of X depends on the left-
and right-sidedness of x as follows:
x
left sided right sided ROC of X
no no strip
no yes RHP
yes no LHP
yes yes everywhere
Thus, we can infer that, if X exists, its ROC can only be of the form of a
LHP, a RHP, a vertical strip, or the entire complex plane.
For example, the sets shown below would not be valid as ROCs.
Im Im
Re Re
Invalid Invalid
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 291
Section 7.3
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 292
Properties of the Laplace Transform
Property
Initial Value Theorem x(0+ ) = lim sX(s)
s→∞
Final Value Theorem lim x(t) = lim sX(s)
t→∞ s→0
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 293
Laplace Transform Pairs
Pair x(t) X(s) ROC
1 δ(t) 1 All s
1
2 u(t) s Re(s) > 0
1
3 −u(−t) s Re(s) < 0
n!
4 t n u(t) sn+1
Re(s) > 0
n!
5 −t n u(−t) sn+1
Re(s) < 0
1
6 e−at u(t) s+a Re(s) > −a
1
7 −e−at u(−t) s+a Re(s) < −a
n!
8 t n e−at u(t) (s+a)n+1
Re(s) > −a
n!
9 −t n e−at u(−t) (s+a)n+1
Re(s) < −a
s
10 cos(ω0t)u(t) s2 +ω20
Re(s) > 0
ω0
11 sin(ω0t)u(t) s2 +ω20
Re(s) > 0
s+a
12 e−at cos(ω0t)u(t) (s+a)2 +ω20
Re(s) > −a
ω0
13 e−at sin(ω0t)u(t) (s+a)2 +ω20
Re(s) > −a
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 294
Linearity
LT LT
If x1 (t) ←→ X1 (s) with ROC R1 and x2 (t) ←→ X2 (s) with ROC R2 , then
LT
a1 x1 (t) + a2 x2 (t) ←→ a1 X1 (s) + a2 X2 (s) with ROC R containing R1 ∩ R2 ,
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 295
Time-Domain Shifting
LT
If x(t) ←→ X(s) with ROC R, then
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 296
Laplace-Domain Shifting
LT
If x(t) ←→ X(s) with ROC R, then
Re Re
σmin σmax σmin + Re(s0 ) σmax + Re(s0 )
R R0 = R + Re(s0 )
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 297
Time-Domain/Laplace-Domain Scaling
LT
If x(t) ←→ X(s) with ROC R, then
1 s
with ROC R0 = aR,
LT
x(at) ←→ X
|a| a
where a is a nonzero real constant.
This is known as the (time-domain/Laplace-domain) scaling property
of the Laplace transform.
As illustrated below, the ROC R is scaled and possibly flipped left to right.
Im Im Im
Re Re Re
σmin σmax aσmin aσmax aσmax aσmin
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 298
Conjugation
LT
If x(t) ←→ X(s) with ROC R, then
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 299
Time-Domain Convolution
LT LT
If x1 (t) ←→ X1 (s) with ROC R1 and x2 (t) ←→ X2 (s) with ROC R2 , then
LT
x1 ∗ x2 (t) ←→ X1 (s)X2 (s) with ROC R containing R1 ∩ R2 .
This is known as the time-domain convolution property of the Laplace
transform.
The ROC R always contains R1 ∩ R2 but can be larger than this
intersection (if pole-zero cancellation occurs).
Convolution in the time domain becomes multiplication in the Laplace
domain.
Consequently, it is often much easier to work with LTI systems in the
Laplace domain, rather than the time domain.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 300
Time-Domain Differentiation
LT
If x(t) ←→ X(s) with ROC R, then
dx(t) LT
←→ sX(s) with ROC R0 containing R.
dt
This is known as the time-domain differentiation property of the
Laplace transform.
The ROC R0 always contains R but can be larger than R (if pole-zero
cancellation occurs).
Differentiation in the time domain becomes multiplication by s in the
Laplace domain.
Consequently, it can often be much easier to work with differential
equations in the Laplace domain, rather than the time domain.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 301
Laplace-Domain Differentiation
LT
If x(t) ←→ X(s) with ROC R, then
LT dX(s)
−tx(t) ←→ with ROC R.
ds
This is known as the Laplace-domain differentiation property of the
Laplace transform.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 302
Time-Domain Integration
LT
If x(t) ←→ X(s) with ROC R, then
Z t
1
x(τ)dτ ←→ X(s) with ROC R0 containing R ∩ {Re(s) > 0}.
LT
−∞ s
This is known as the time-domain integration property of the Laplace
transform.
The ROC R0 always contains at least R ∩ {Re(s) > 0} but can be larger (if
pole-zero cancellation occurs).
Integration in the time domain becomes division by s in the Laplace
domain.
Consequently, it is often much easier to work with integral equations in the
Laplace domain, rather than the time domain.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 303
Initial Value Theorem
where x(0+ ) denotes the limit of x(t) as t approaches zero from positive
values of t .
This result is known as the initial value theorem.
In situations where X is known but x is not, the initial value theorem
eliminates the need to explicitly find x by an inverse Laplace transform
calculation in order to evaluate x(0+ ).
In practice, the values of functions at the origin are frequently of interest,
as such values often convey information about the initial state of systems.
The initial value theorem can sometimes also be helpful in checking for
errors in Laplace transform calculations.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 304
Final Value Theorem
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 305
More Laplace Transform Examples
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 306
Section 7.4
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 307
Finding Inverse Laplace Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 308
Section 7.5
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 309
System Function of LTI Systems
Consider a LTI system with input x, output y, and impulse response h. Let
X , Y , and H denote the Laplace transforms of x, y, and h, respectively.
Since y(t) = x ∗ h(t), the system is characterized in the Laplace domain by
Y (s) = X(s)H(s).
As a matter of terminology, we refer to H as the system function (or
transfer function) of the system (i.e., the system function is the Laplace
transform of the impulse response).
A LTI system is completely characterized by its system function H .
When viewed in the Laplace domain, a LTI system forms its output by
multiplying its input with its system function.
If the ROC of H includes the imaginary axis, then H( jω) is the frequency
response of the LTI system.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 310
Block Diagram Representations of LTI Systems
Consider a LTI system with input x, output y, and impulse response h, and
let X , Y , and H denote the Laplace transforms of x, y, and h, respectively.
Often, it is convenient to represent such a system in block diagram form in
the Laplace domain as shown below.
X Y
H
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 311
Interconnection of LTI Systems
X Y
≡ H1 + H2
H2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 312
Causality
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 313
BIBO Stability
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 314
Invertibility
H(s)Hinv (s) = 1,
1
Hinv (s) = .
H(s)
Since distinct systems can have identical system functions (but with
differing ROCs), the inverse of a LTI system is not necessarily unique.
In practice, however, we often desire a stable and/or causal system. So,
although multiple inverse systems may exist, we are frequently only
interested in one specific choice of inverse system (due to these
additional constraints of stability and/or causality).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 315
LTI Systems and Differential Equations
Many LTI systems of practical interest can be represented using an
Nth-order linear differential equation with constant coefficients.
Consider a system with input x and output y that is characterized by an
equation of the form
N M
d k d k
∑ bk dt y(t) = ∑ ak dt x(t),
k=0 k=0
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 316
Section 7.6
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 317
Electronic Circuits
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 318
Resistors
R
i + −
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 319
Inductors
An inductor is a circuit element that converts an electric current into a
magnetic field and vice versa.
An inductor uses the energy stored in a magnetic field in order to oppose
changes in current (through the inductor).
An inductor is characterized by an equation of the form
Z t
d 1
v(t) = L dt i(t) (or equivalently, i(t) = L v(τ)dτ),
−∞
where L is a nonnegative real constant, and v and i respectively denote
the voltage across and current through the inductor as a function of time.
As a matter of terminology, the quantity L is known as the inductance of
the inductor.
Inductance is measured in units of henrys (H).
In circuit diagrams, an inductor is denoted by the symbol shown below.
L
i + −
v
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 320
Capacitors
A capacitor is a circuit element that stores electric charge.
A capacitor uses the energy stored in an electric field in order to oppose
changes in voltage (across the capacitor).
A capacitor is characterized by an equation of the form
Z t
1
v(t) = C i(τ)dτ (or equivalently, i(t) = C dtd v(t)),
−∞
where C is a nonnegative real constant, and v and i respectively denote
the voltage across and current through the capacitor as a function of time.
As a matter of terminology, the quantity C is known as the capacitance of
the capacitor.
Capacitance is measured in units of farads (F).
In circuit diagrams, a capacitor is denoted by the symbol shown below.
C
i + −
v
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 321
Circuit Analysis with the Laplace Transform
The Laplace transform is a very useful tool for circuit analysis.
The utility of the Laplace transform is partly due to the fact that the
differential/integral equations that describe inductors and capacitors are
much simpler to express in the Laplace domain than in the time domain.
Let v and i denote the voltage across and current through a circuit
element, and let V and I denote the Laplace transforms of v and i,
respectively.
In the Laplace domain, the equations characterizing a resistor, an
inductor, and a capacitor respectively become:
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 322
Section 7.7
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 323
Control Systems
A control system manages the behavior of one or more other systems with
some specific goal.
Typically, the goal is to force one or more physical quantities to assume
particular desired values, where such quantities might include: positions,
velocities, accelerations, forces, torques, temperatures, or pressures.
The desired values of the quantities being controlled are collectively
viewed as the input of the control system.
The actual values of the quantities being controlled are collectively viewed
as the output of the control system.
A control system whose behavior is not influenced by the actual values of
the quantities being controlled is called an open loop (or non-feedback)
system.
A control system whose behavior is influenced by the actual values of the
quantities being controlled is called a closed loop (or feedback) system.
An example of a simple control system would be a thermostat system,
which controls the temperature in a room or building.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 324
Feedback Control Systems
Reference
Input Error Output
+ Controller Plant
−
Sensor
Feedback
Signal
input: desired value of the quantity to be controlled
output: actual value of the quantity to be controlled
error: difference between the desired and actual values
plant: system to be controlled
sensor: device used to measure the actual output
controller: device that monitors the error and changes the input of the
plant with the goal of forcing the error to zero
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 325
Stability Analysis of Feedback Systems
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 326
Stabilization Example: Unstable Plant
X Y
P
10
P(s) = s−1
ROC of P:
Im
1 Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 327
Stabilization Example: Using Pole-Zero Cancellation
X Y
W P
10 s−1
P(s) = s−1 , W (s) = 10(s+1)
ROC of H :
Im
−1 Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 328
Stabilization Example: Using Feedback (1)
feedback system (with causal LTI compensator and sensor):
X R Y
+ C P
−
10
P(s) = s−1 , C(s) = β, Q(s) = 1
system function H of feedback system:
C(s)P(s) 10β
H(s) = 1+C(s)P(s)Q(s) = s−(1−10β)
ROC of H :
Im
Re
1 − 10β
Y (s) = C(s)P(s)R(s)
= C(s)P(s)[X(s) − Q(s)Y (s)]
= C(s)P(s)X(s) −C(s)P(s)Q(s)Y (s)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 330
Stabilization Example: Using Feedback (3)
10
P(s) = s−1 , C(s) = β, Q(s) = 1
C(s)P(s)
H(s) =
1 +C(s)P(s)Q(s)
10
β( s−1 )
= 10
1 + β( s−1 )(1)
10β
=
s − 1 + 10β
10β
=
s − (1 − 10β)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 331
Remarks on Stabilization Via Pole-Zero Cancellation
Pole-zero cancellation is not achievable in practice, and therefore it cannot
be used to stabilize real-world systems.
The theoretical models used to represent real-world systems are only
approximations due to many factors, including the following:
2 Determining the system function of a system involves measurement, which
always has some error.
2 A system cannot be built with such precision that it will have exactly some
prescribed system function.
2 The system function of most systems will vary at least slightly with changes
in the physical environment.
2 Although a LTI model is used to represent a system, the likely reality is that
the system is not exactly LTI, which introduces error.
Due to approximation error, the effective poles and zeros of the system
function will only be approximately where they are expected to be.
Since pole-zero cancellation requires that a pole and zero be placed at
exactly the same location, any error will prevent this cancellation from
being achieved.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 332
Section 7.8
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 333
Unilateral Laplace Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 334
Inversion of the Unilateral Laplace Transform
L−1 −1
u {Lu x}(t) = Lu {L{xu}}(t)
= L−1 {L{xu}}(t)
= x(t)u(t)
(
x(t) t ≥ 0
=
0 t < 0.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 335
Unilateral Versus Bilateral Laplace Transform
Due to the close relationship between the unilateral and bilateral Laplace
transforms, these two transforms have some similarities in their properties.
Since these two transforms are not identical, however, their properties
differ in some cases, often in subtle ways.
In the unilateral case, we have that:
1 the time-domain convolution property has the additional requirement that
the functions being convolved must be causal;
2 the time/Laplace-domain scaling property has the additional constraint that
the scaling factor must be positive;
3 the time-domain differentiation property has an extra term in the expression
for Lu {Dx}(t), where D denotes the derivative operator (namely, −x(0− ));
4 the time-domain integration property has a different lower limit in the
time-domain integral (namely, 0− instead of −∞); and
5 the time-domain shifting property does not hold (except in special
circumstances).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 336
Properties of the Unilateral Laplace Transform
Property
Initial Value Theorem x(0+ ) = lim sX(s)
s→∞
Final Value Theorem lim x(t) = lim sX(s)
t→∞ s→0
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 337
Unilateral Laplace Transform Pairs
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 338
Solving Differential Equations [Using the Unilateral Laplace Transform]
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 339
Part 8
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 340
Section 8.1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 341
Time Shifting (Translation)
Time shifting (also called translation) maps the input sequence x to the
output sequence y as given by
where b is an integer.
Such a transformation shifts the sequence (to the left or right) along the
time axis.
If b > 0, y is shifted to the right by |b|, relative to x (i.e., delayed in time).
If b < 0, y is shifted to the left by |b|, relative to x (i.e., advanced in time).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 342
Time Shifting (Translation): Example
x(n)
n
−4 −3 −2 −1 0 1 2 3 4
x(n − 1) x(n + 1)
3 3
2 2
1 1
n n
−4 −3 −2 −1 0 1 2 3 4 −4 −3 −2 −1 0 1 2 3 4
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 343
Time Reversal (Reflection)
y(n) = x(−n).
Geometrically, the output sequence y is a reflection of the input sequence
x about the (vertical) line n = 0.
x(n) x(−n)
3 3
2 2
1 1
n n
−3 −2 −1 0 1 2 3 −3 −2 −1 0 1 2 3
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 344
Downsampling
3 3
2 2
1 1
n n
−4 −3 −2 −1 0 1 2 3 4 −4 −3 −2 −1 0 1 2 3 4
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 345
Upsampling
3 3
2 2
1 1
n n
−5 −4 −3 −2 −1 0 1 2 3 4 5 −5 −4 −3 −2 −1 0 1 2 3 4 5
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 346
Combined Independent-Variable Transformations
Consider a transformation that maps the input sequence x to the output
sequence y as given by
2 then, downsampling the result by |a| and, if a < 0, time reversing as well.
Properties of Sequences
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 348
Symmetry and Addition/Multiplication
Sums involving even and odd sequences have the following properties:
2 The sum of two even sequences is even.
2 The sum of two odd sequences is odd.
2 The sum of an even sequence and odd sequence is neither even nor odd,
provided that neither of the sequences is identically zero.
That is, the sum of sequences with the same type of symmetry also has
the same type of symmetry.
Products involving even and odd sequences have the following
properties:
2 The product of two even sequences is even.
2 The product of two odd sequences is even.
2 The product of an even sequence and an odd sequence is odd.
That is, the product of sequences with the same type of symmetry is even,
while the product of sequences with opposite types of symmetry is odd.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 349
Decomposition of a Sequence into Even and Odd Parts
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 350
Sum of Periodic Sequences
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 351
Right-Sided Sequences
A sequence x is said to be right sided if, for some (finite) integer constant
n0 , the following condition holds:
x(n) = 0 for all n < n0
(i.e., x is only potentially nonzero to the right of n0 ).
An example of a right-sided sequence is shown below.
x(n)
2 ···
n
−4 −3 −2 −1 0 1 2 3 4
x(n)
··· 2
n
−4 −3 −2 −1 0 1 2 3 4
n
−4 −3 −2 −1 0 1 2 3 4
A sequence that is neither left sided nor right sided is said to be two
sided.
An example of a two-sided sequence is shown below.
x(n)
··· 2 ···
n
−4 −3 −2 −1 0 1 2 3 4
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 354
Bounded Sequences
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 355
Energy of a Sequence
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 356
Section 8.3
Elementary Sequences
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 357
Real Sinusoidal Sequences
A real sinusoidal sequence is a sequence of the form
x(n) = A cos(Ωn + θ),
where A, Ω, and θ are real constants.
A real sinusoid is periodic if and only if Ω is a rational number, in which
2π
case the fundamental period is the smallest integer of the form 2πk
|Ω| where
k is a (strictly) positive integer.
For all integer k, xk (n) = A cos([Ω + 2πk]n + θ) is the same sequence.
An example of a periodic real sinusoid with fundamental period 12 is
shown plotted below.
π
x(n) = cos 6n
··· ···
n
−12 12
−1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 358
Oscillation Rate of Real Sinusoidal Sequences
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 359
Effect of Increasing Frequency on Oscillation Rate
0π
π
cos(0n) = cos 8 n cos 8n
1 1
n n
−16 16 −16 16
−1 −1
π
4π
8π
cos 2n = cos 8 n cos(πn) = cos 8 n
1 1
n n
−16 16 −16 16
−1 −1
15π
16π
cos 8 n cos(2πn) = cos 8 n
1 1
n n
−16 16 −16 16
−1 −1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 360
Complex Exponential Sequences
x(n) = can ,
x(n) = cebn ,
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 361
Real Exponential Sequences
x(n) = can ,
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 362
Real Exponential Sequences (Continued 1)
x(n) x(n)
3 3
2 2
1 ··· ··· 1
··· ···
n n
−4 −3 −2 −1 0 1 2 3 4 −4 −3 −2 −1 0 1 2 3 4
x(n)
1
··· ···
n
−4 −3 −2 −1 0 1 2 3 4
|a| = 1, a > 0 [ a = 1; c = 1]
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 363
Real Exponential Sequences (Continued 2)
x(n) x(n)
3 3
2 2
··· ···
1 1
··· ···
n n
−4 −3 −2 −1 1 2 3 4 −4 −3 −2 −1 1 2 3 4
1 1
2 2
3 3
x(n)
1
··· ···
n
−4 −3 −2 −1 1 2 3 4
−1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 364
Complex Sinusoidal Sequences
x(n) = ce jΩn ,
··· ···
n
−7 7
−1
2π
Im{e j(2π/7)n } = sin 7 n
··· ···
n
−7 7
−1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 366
Oscillation Rate of Complex Sinusoidal Sequences
x(n) = ce jΩn ,
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 367
General Complex Exponential Sequences
In the most general case of a complex exponential sequence x(n) = can ,
c and a are both complex.
Letting c = |c| e jθ and a = |a| e jΩ where θ and Ω are real, and using
Euler’s relation, we can rewrite x(n) as
Thus, Re{x} and Im{x} are each the product of a real exponential and
real sinusoid.
One of several distinct modes of behavior is exhibited by x, depending on
the value of a.
If |a| = 1, Re{x} and Im{x} are real sinusoids.
If |a| > 1, Re{x} and Im{x} are each the product of a real sinusoid and
a growing real exponential.
If |a| < 1, Re{x} and Im{x} are each the product of a real sinusoid and
a decaying real exponential.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 368
General Complex Exponential Sequences (Continued)
The various modes of behavior for Re{x} and Im{x} are illustrated
below.
|a| = 1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 369
Relationship Between Complex Exponentials and Real
Sinusoids
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 370
Unit-Step Sequence
u(n)
1
···
n
−3 −2 −1 0 1 2 3
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 371
Unit Rectangular Pulses
p(n)
1
···
n
a−3 a−2 a−1 a a+1 a+2 a+3 b−1 b b+1 b+2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 372
Unit-Impulse Sequence
The unit-impulse sequence (also known as the delta sequence), denoted
δ, is defined as
(
1 n=0
δ(n) =
0 otherwise.
The first-order difference of u is δ. That is,
δ(n) = u(n) − u(n − 1).
The running sum of δ is u. That is,
n
u(n) = ∑ δ(k).
k=−∞
A plot of δ is shown below.
δ(n)
n
−3 −2 −1 0 1 2 3
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 373
Properties of the Unit-Impulse Sequence
For any sequence x and any integer constant n0 , the following identity
holds:
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 374
Representing Rectangular Pulses (Using Unit-Step Sequences)
The idea from the previous slide can be extended to handle any sequence
that is defined in a piecewise manner (i.e., via an expression involving
multiple cases).
That is, by using unit-step sequences, we can always collapse a formula
involving multiple cases into a single expression.
Often, simplifying a formula in this way can be quite beneficial.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 376
Section 8.4
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 377
DT Systems
A system with input x and output y can be described by the equation
y = Hx,
x → y.
Note that the symbols “→” and “=” have very different meanings.
The symbol “→” should be read as “produces” (not as “equals”).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 378
Block Diagram Representations
Often, a system defined by the operator H and having the input x and
output y is represented in the form of a block diagram as shown below.
Input Output
x System y
H
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 379
Interconnection of Systems
Two basic ways in which systems can be interconnected are shown below.
x y
H1 +
x y
H1 H2
H2
Series
Parallel
A series (or cascade) connection ties the output of one system to the input
of the other.
The overall series-connected system is described by the equation
y = H2 H1 x.
A parallel connection ties the inputs of both systems together and sums
their outputs.
The overall parallel-connected system is described by the equation
y = H1 x + H2 x.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 380
Section 8.5
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 381
Memory
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 382
Memory (Continued)
n
−∞ n0 ∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 383
Causality
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 384
Causality (Continued)
n ≤ n0
n
−∞ n0 ∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 385
Invertibility
The inverse of a system H is another system H−1 such that, for every
sequence x,
H−1 Hx = x
(i.e., the system formed by the cascade interconnection of H followed by
H−1 is a system whose input and output are equal).
A system is said to be invertible if it has a corresponding inverse system
(i.e., its inverse exists).
Equivalently, a system is invertible if its input x can always be uniquely
determined from its output y.
An invertible system will always produce distinct outputs from any two
distinct inputs.
To show that a system is invertible, we simply find the inverse system.
To show that a system is not invertible, we find two distinct inputs that
result in identical outputs.
In practical terms, invertible systems are “nice” in the sense that their
effects can be undone.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 386
Invertibility (Continued)
A system H−1 being the inverse of H means that the following two
systems are equivalent (i.e., H−1 H is an identity):
x y x y
H H−1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 387
Bounded-Input Bounded-Output (BIBO) Stability
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 388
Time Invariance (TI)
A system H is said to be time invariant (TI) (or shift invariant (SI)) if,
for every sequence x and every integer n0 , the following condition holds:
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 389
Time Invariance (Continued)
x y
Sn0 H x y
H Sn0
System 1: y = HSn0 x
System
2: y = Sn0 Hx
y(n) = Hx0 (n)
y(n) = Hx(n − n0 )
x0 (n) = Sn0 x(n) = x(n − n0 )
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 390
Additivity, Homogeneity, and Linearity
A system H is said to be additive if, for all sequences x1 and x2 , the
following condition holds:
H(x1 + x2 ) = Hx1 + Hx2
(i.e., H commutes with sums).
A system H is said to be homogeneous if, for every sequence x and every
complex constant a, the following condition holds:
H(ax) = aHx
(i.e., H commutes with multiplication by a constant).
A system that is both additive and homogeneous is said to be linear.
In other words, a system H is linear, if for all sequences x1 and x2 and all
complex constants a1 and a2 , the following condition holds:
H(a1 x1 + a2 x2 ) = a1 Hx1 + a2 Hx2
(i.e., H commutes with linear combinations).
The linearity property is also referred to as the superposition property.
Practically speaking, linear systems are much easier to design and
analyze than nonlinear systems.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 391
Additivity, Homogeneity, and Linearity (Continued 1)
The system H is additive if and only if the following two systems are
equivalent (i.e., H commutes with addition):
x1 y x1 y
+ H H +
x2 x2
H
System 1: y = H(x1 + x2 )
System 2: y = Hx1 + Hx2
x y x y
a H H a
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 392
Additivity, Homogeneity, and Linearity (Continued 2)
The system H is linear if and only if the following two systems are
equivalent (i.e., H commutes with linear combinations):
x1 y x1 y
a1 + H H a1 +
x2 x2
a2 H a2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 393
Eigensequences of Systems
Hx = λx,
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 394
Part 9
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 395
Why Linear Time-Invariant (LTI) Systems?
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 396
Section 9.1
Convolution
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 397
DT Convolution
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 398
Practical Convolution Computation
we proceed as follows:
1 Plot x(k) and h(n − k) as a function of k .
2 Initially, consider an arbitrarily large negative value for n. This will result in
h(n − k) being shifted very far to the left on the time axis.
3 Write the mathematical expression for x ∗ h(n).
4 Increase n gradually until the expression for x ∗ h(n) changes form. Record
the interval over which the expression for x ∗ h(n) was valid.
5 Repeat steps 3 and 4 until n is an arbitrarily large positive value. This
corresponds to h(n − k) being shifted very far to the right on the time axis.
6 The results for the various intervals can be combined in order to obtain an
expression for x ∗ h(n) for all n.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 399
Properties of Convolution
The convolution operation is commutative. That is, for any two sequences
x and h,
x ∗ h = h ∗ x.
The convolution operation is associative. That is, for any sequences x, h1 ,
and h2 ,
(x ∗ h1 ) ∗ h2 = x ∗ (h1 ∗ h2 ).
The convolution operation is distributive with respect to addition. That is,
for any sequences x, h1 , and h2 ,
x ∗ (h1 + h2 ) = x ∗ h1 + x ∗ h2 .
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 400
Representation of Sequences Using Impulses
x ∗ δ = x.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 401
Circular Convolution
The convolution of two periodic sequences is usually not well defined.
This motivates an alternative notion of convolution for periodic sequences
known as circular convolution.
The circular convolution (also known as the DT periodic convolution) of
the N -periodic sequences x and h, denoted x ~ h, is defined as
N−1
x ~ h(n) = ∑ x(k)h(n − k) = ∑ x(k)h(mod(n − k, N)),
k=hNi k=0
(i.e., x0 (n) equals x(n) over a single period of x and is zero elsewhere).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 402
Section 9.2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 403
Impulse Response
y = x ∗ h.
In other words, a LTI system simply computes a convolution.
Furthermore, a LTI system is completely characterized by its impulse
response.
That is, if the impulse response of a LTI system is known, we can
determine the response of the system to any input.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 404
Step Response
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 405
Block Diagram of LTI Systems
x y
h
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 406
Interconnection of LTI Systems
x y x y
h1 h2 ≡ h1 ∗ h2
x y x y
h1 h2 ≡ h2 h1
x y
h1 +
x y
≡ h1 + h2
h2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 407
Section 9.3
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 408
Memory
h(n) = Kδ(n),
y = x ∗ (Kδ) = Kx
a causal sequence.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 410
Invertibility
h ∗ hinv = δ.
Consequently, a LTI system with impulse response h is invertible if and
only if there exists a sequence hinv such that
h ∗ hinv = δ.
Except in simple cases, the above condition is often quite difficult to test.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 411
BIBO Stability
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 412
Eigensequences of LTI Systems
As it turns out, every complex exponential is an eigensequence of all LTI
systems.
For a LTI system H with impulse response h,
zn
That is, is an eigensequence of a LTI system and H(z) is the
corresponding eigenvalue.
We refer to H as the system function (or transfer function) of the
system H.
From above, we can see that the response of a LTI system to a complex
exponential is the same complex exponential multiplied by the complex
factor H(z).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 413
Representation of Sequences Using Eigensequences
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 414
Part 10
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 415
Introduction
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 416
Section 10.1
Fourier Series
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 417
Harmonically-Related Complex Sinusoids
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 418
DT Fourier Series (DTFS)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 419
DT Fourier Series (DTFS) (Continued)
ak = 1
N ∑ x(n)e− j(2π/N)kn .
n=hNi
(The summation can be taken over any N consecutive integers due to the
N -periodic nature of x and e− j(2π/N)kn .)
The above equation for ak is often referred to as the Fourier series
analysis equation.
Due to the N -periodic nature of x and e− j(2π/N)kn , the sequence a is also
N -periodic.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 420
Convergence of Fourier Series
Since the analysis and synthesis equations for (DT) Fourier series involve
only finite sums (as opposed to infinite series), convergence is not a
significant issue of concern.
If an N -periodic sequence is bounded (i.e., is finite in value), its Fourier
series coefficient sequence will exist and be bounded and the Fourier
series analysis and synthesis equations must converge.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 421
Section 10.2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 422
Properties of (DT) Fourier Series
DTFS DTFS
x(n) ←→ ak and y(n) ←→ bk
Property
1 2 2
Parseval’s Relation N ∑n=hNi |x(n)| = ∑k=hNi |ak |
Even Symmetry x is even ⇔ a is even
Odd Symmetry x is odd ⇔ a is odd
Real / Conjugate Symmetry x is real ⇔ a is conjugate symmetric
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 423
Linearity
DTFS DTFS
Let x and y be N -periodic sequences. If x(n) ←→ ak and y(n) ←→ bk ,
then
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 424
Translation (Time Shifting)
DTFS
Let x denote a periodic sequence with period N . If x(n) ←→ ck , then
x(n − n0 ) ←→ e− jk(2π/N)n0 ck ,
DTFS
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 425
Modulation (Frequency Shifting)
DTFS
Let x denote a periodic sequence with period N . If x(n) ←→ ck , then
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 426
Reflection (Time Reversal)
DTFS
Let x denote a periodic sequence with period N . If x(n) ←→ ck , then
DTFS
x(−n) ←→ c−k .
That is, time reversing a sequence results in a time reversal of the
corresponding Fourier series coefficient sequence.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 427
Conjugation
DTFS
Let x denote a periodic sequence with period N . If x(n) ←→ ck , then
x∗ (n) ←→ c∗−k .
DTFS
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 428
Duality
DTFS
Let x denote a periodic sequence with period N . If x(n) ←→ a(k), then
a(n) ←→ N1 x(−k).
DTFS
That is, the analysis and synthesis equations are identical except for a
factor of N and different sign in the parameter for the exponential
function.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 429
Periodic Convolution
DTFS DTFS
Let x and y be N -periodic sequences. If x(n) ←→ ak and y(n) ←→ bk ,
then
x ~ y(n) ←→ Nak bk .
DTFS
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 430
Multiplication
DTFS DTFS
Let x and y be N -periodic sequences. If x(n) ←→ ak and y(n) ←→ bk ,
then
x(n)y(n) ←→ a ~ b(k).
DTFS
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 431
Parseval’s Relation
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 432
Even/Odd Symmetry
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 433
Real Sequences
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 434
Trigonometric Form of a Fourier Series
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 435
Other Properties of Fourier Series
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 436
Section 10.3
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 437
Prelude to the Discrete Fourier Transform (DFT)
Letting a0k = Nak , we can rewrite the Fourier series synthesis and analysis
equations, respectively, as
N−1 N−1
x(n) = 1
N ∑ a0k e j(2π/N)kn and a0k = ∑ x(n)e− j(2π/N)kn .
k=0 n=0
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 439
Properties of Discrete Fourier Transform (DFT)
Property
2 2
Parseval’s Relation ∑N−1 1 N−1
n=0 |x(n)| = N ∑k=0 |X(k)|
Even Symmetry x is even ⇔ X is even
Odd Symmetry x is odd ⇔ X is odd
Real / Conjugate Symmetry x is real ⇔ X is conjugate symmetric
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 440
Section 10.4
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 441
A New Perspective on Sequences: The Frequency Domain
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 442
Fourier Series and Frequency Spectra
To gain further insight into the role played by the Fourier series
coefficients ak in the context of the frequency spectrum of the N -periodic
sequence x, it is helpful to write the Fourier series with the ak expressed in
polar form as
N−1 N−1
x(n) = ∑ ak e j(2π/N)kn = ∑ |ak | e j([2π/N]kn+arg a ) .
k
k=0 k=0
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 444
Fourier Series and Frequency Spectra (Continued 2)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 445
Frequency Spectra of Real Sequences
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 446
Section 10.5
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 447
Frequency Response
Recall that a LTI system H with impulse response h is such that
H{zn }(n) = HZ (z)zn , where HZ (z) = ∑∞ −n
n=−∞ h(n)z . (That is, complex
exponentials are eigensequences of LTI systems.)
Since a complex sinusoid is a special case of a complex exponential, we
can reuse the above result for the special case of complex sinusoids.
For a LTI system H with impulse response h,
H e jΩn (n) = H(Ω)e jΩn ,
where Ω is real and
∞
H(Ω) = ∑ h(n)e− jΩn .
n=−∞
k=0
Thus, if the input x to a LTI system is a Fourier series, the output y is also a
DTFS DTFS
Fourier series. More specifically, if x(n) ←→ ak then y(n) ←→ H(kΩ0 )ak .
The above formula can be used to determine the output of a LTI system
from its input in a way that does not require convolution.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 449
Filtering
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 450
Ideal Lowpass Filter
An ideal lowpass filter eliminates all baseband frequency components
with a frequency whose magnitude is greater than some cutoff frequency,
while leaving the remaining baseband frequency components unaffected.
Such a filter has a frequency response of the form
(
1 |Ω| ≤ Ωc
H(Ω) =
0 Ωc < |Ω| ≤ π,
H(Ω)
1
Ω
−π −Ωc Ωc π
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 451
Ideal Highpass Filter
An ideal highpass filter eliminates all baseband frequency components
with a frequency whose magnitude is less than some cutoff frequency,
while leaving the remaining baseband frequency components unaffected.
Such a filter has a frequency response of the form
(
1 Ωc < |Ω| ≤ π
H(Ω) =
0 |Ω| ≤ Ωc ,
H(Ω)
Ω
−π −Ωc Ωc π
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 452
Ideal Bandpass Filter
An ideal bandpass filter eliminates all baseband frequency components
with a frequency whose magnitude does not lie in a particular range, while
leaving the remaining baseband frequency components unaffected.
Such a filter has a frequency response of the form
(
1 Ωc1 ≤ |Ω| ≤ Ωc2
H(Ω) =
0 |Ω| < Ωc1 or Ωc2 < |Ω| < π,
H(Ω)
Ω
−π −Ωc2 −Ωc1 Ωc1 Ωc2 π
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 453
Part 11
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 454
Motivation for the Fourier Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 455
Section 11.1
Fourier Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 456
Development of the Fourier Transform [Aperiodic Case]
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 457
Development of the Fourier Transform [Aperiodic Case] (Continued)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 458
Generalized Fourier Transform
The classical Fourier transform for aperiodic sequences does not exist
(i.e., ∑∞
n=−∞ x(n)e
− jΩn fails to converge) for some sequences of great
practical interest, such as:
2 a nonzero constant sequence;
2 a periodic sequence (e.g., a real or complex sinusoid); and
2 the unit-step sequence (i.e., u).
Fortunately, the Fourier transform can be extended to handle such
sequences, resulting in what is known as the generalized Fourier
transform.
For our purposes, we can think of the classical and generalized Fourier
transforms as being defined by the same formulas.
Therefore, in what follows, we will not typically make a distinction between
the classical and generalized Fourier transforms.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 459
DT Fourier Transform (DTFT)
The Fourier transform of the sequence x, denoted Fx or X , is given by
∞
Fx(Ω) = X(Ω) = ∑ x(n)e− jΩn .
n=−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 461
Convergence of the Fourier Transform
∞
∑ |x(k)| < ∞
k=−∞
∞
∑ |x(k)|2 < ∞
k=−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 462
Section 11.3
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 463
Properties of the (DT) Fourier Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 464
Properties of the (DT) Fourier Transform (Continued)
Property
Periodicity X(Ω) = X(Ω + 2π)
2 2
Parseval’s Relation ∑∞ 1 R
n=−∞ |x(n)| = 2π 2π |X(Ω)| dΩ
Even Symmetry x is even ⇔ X is even
Odd Symmetry x is odd ⇔ X is odd
Real / Conjugate Symmetry x is real ⇔ X is conjugate symmetric
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 465
(DT) Fourier Transform Pairs
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 466
Periodicity
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 467
Linearity
DTFT DTFT
If x1 (n) ←→ X1 (Ω) and x2 (n) ←→ X2 (Ω), then
DTFT
a1 x1 (n) + a2 x2 (n) ←→ a1 X1 (Ω) + a2 X2 (Ω),
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 468
Translation
DTFT
If x(n) ←→ X(Ω), then
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 469
Modulation
DTFT
If x(n) ←→ X(Ω), then
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 470
Conjugation
DTFT
If x(n) ←→ X(Ω), then
x∗ (n) ←→ X ∗ (−Ω).
DTFT
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 471
Time Reversal
DTFT
If x(n) ←→ X(Ω), then
DTFT
x(−n) ←→ X(−Ω).
This is known as the time-reversal property of the Fourier transform.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 472
Upsampling
DTFT
If x(n) ←→ X(Ω), then
DTFT
(↑ M)x(n) ←→ X(MΩ).
This is known as the upsampling property of the Fourier transform.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 473
Downsampling
DTFT
If x(n) ←→ X(Ω), then
M−1
1
∑X Ω−2πk
DTFT
(↓ M)x(n) ←→ M M .
k=0
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 474
Convolution
DTFT DTFT
If x1 (n) ←→ X1 (Ω) and x2 (n) ←→ X2 (Ω), then
DTFT
x1 ∗ x2 (n) ←→ X1 (Ω)X2 (Ω).
This is known as the convolution (or time-domain convolution)
property of the Fourier transform.
In other words, a convolution in the time domain becomes a multiplication
in the frequency domain.
This suggests that the Fourier transform can be used to avoid having to
deal with convolution operations.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 475
Multiplication
DTFT DTFT
If x1 (n) ←→ X1 (Ω) and x2 (n) ←→ X2 (Ω), then
Z
1
X1 (θ)X2 (Ω − θ)dθ.
DTFT
x1 (n)x2 (n) ←→ 2π
2π
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 476
Frequency-Domain Differentiation
DTFT
If x(n) ←→ X(Ω), then
d DTFT
nx(n) ←→ j dΩ X(Ω).
This is known as the frequency-domain differentiation property of the
Fourier transform.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 477
Differencing
DTFT
If x(n) ←→ X(Ω), then
x(n) − x(n − 1) ←→ 1 − e− jΩ X(Ω).
DTFT
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 478
Accumulation
DTFT
If x(n) ←→ X(Ω), then
n ∞
e jΩ
∑ πX(0) ∑ δ(Ω − 2πk).
DTFT
x(k) ←→ X(Ω) +
k=−∞ e jΩ − 1 k=−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 479
Parseval’s Relation
DTFT
If x(n) ←→ X(Ω), then
∞ Z
∑ |x(n)|2 = 1
2π
2π
|X(Ω)|2 dΩ
n=−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 480
Even and Odd Symmetry
2 x is odd ⇔ X is odd.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 481
Real Sequences
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 482
Section 11.4
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 483
Fourier Transform of Periodic Sequences
The Fourier transform can be generalized to also handle periodic
sequences.
Consider an N -periodic sequence x.
Define the sequence xN as
(
x(n) 0 ≤ n < N
xN (n) =
0 otherwise.
(i.e., xN (n) is equal to x(n) over a single period and zero elsewhere).
Let a denote the Fourier series coefficient sequence of x.
Let X and XN denote the Fourier transforms of x and xN , respectively.
The following relationships can be shown to hold:
∞
X(Ω) = 2π
N ∑ XN 2πk
N δ Ω − 2πk
N ,
k=−∞
∞
ak = N1 XN 2πk
N , and X(Ω) = 2π ∑ ak δ Ω − 2πk
N .
k=−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 484
Fourier Transform of Periodic Sequences (Continued)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 485
Section 11.5
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 486
Frequency Spectra of Sequences
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 487
Fourier Transform and Frequency Spectra
To gain further insight into the role played by the Fourier transform X in
the context of the frequency spectrum of x, it is helpful to write the Fourier
transform representation of x with X(Ω) expressed in polar form as
follows:
Z Z
x(n) = 1
2π X(Ω)e jΩn dΩ = 1
2π |X(Ω)| e j[Ωn+arg X(Ω)] dΩ.
2π 2π
In effect, the quantity |X(Ω)| is a weight that determines how much the
complex sinusoid at frequency Ω contributes to the integration result x(n).
Perhaps, this can be more easily seen if we express the above integral as
the limit of a sum, derived from an approximation of the integral using the
area of rectangles, as shown on the next slide. [Recall that
Rb n b−a
a f (x)dx = limn→∞ ∑k=1 f (xk )∆x where ∆x = n and xk = a + k∆x.]
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 488
Fourier Transform and Frequency Spectra (Continued 1)
Expressing the integral (from the previous slide) as the limit of a sum, we
obtain
`
∑ ∆Ω
0 0
x(n) = lim 1
X(Ω0 ) e j[Ω n+arg X(Ω )] ,
`→∞ 2π
k=1
where ∆Ω = 2π
` and Ω = k∆Ω.
0
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 490
Frequency Spectra of Real Sequences
Recall that, for a real sequence x, the Fourier transform X of x satisfies
X(Ω) = X ∗ (−Ω)
Ω
−π −B B π
where
Ex (Ω) = |X(Ω)|2 .
We refer to Ex as the energy-density spectrum of the sequence x.
The function Ex indicates how the energy in x is distributed with respect to
frequency.
For example, the energy contributed by frequencies in the range [Ω1 , Ω2 ]
is given by
Z Ω2
1
2π Ex (Ω)dΩ.
Ω1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 493
Section 11.6
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 494
Frequency Response of LTI Systems
Consider a LTI system with input x, output y, and impulse response h, and
let X , Y , and H denote the Fourier transforms of x, y, and h, respectively.
Since y(n) = x ∗ h(n), we have that
Y (Ω) = X(Ω)H(Ω).
The function H is called the frequency response of the system.
A LTI system is completely characterized by its frequency response H .
The above equation provides an alternative way of viewing the behavior of
a LTI system. That is, we can view the system as operating in the
frequency domain on the Fourier transforms of the input and output
signals.
The frequency spectrum of the output is the product of the frequency
spectrum of the input and the frequency response of the system.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 495
Frequency Response of LTI Systems (Continued 1)
|Y (Ω)| = |X(Ω)| |H(Ω)| and argY (Ω) = arg X(Ω) + arg H(Ω).
The magnitude spectrum of the output equals the magnitude spectrum of
the input times the magnitude response of the system.
The phase spectrum of the output equals the phase spectrum of the input
plus the phase response of the system.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 496
Frequency Response of LTI Systems (Continued 2)
(i.e., the magnitude response |H(Ω)| is even and the phase response
arg H(Ω) is odd).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 497
Unwrapped Phase
For many types of analysis, restricting the range of a phase function to an
interval of length 2π (such as (−π, π]), often unnecessarily introduces
discontinuities into the function.
This motivates the notion of unwrapped phase.
The unwrapped phase is simply the phase defined in such a way so as
not to restrict the phase to an interval of length 2π and to keep the phase
function continuous to the greatest extent possible.
For example, the function H(Ω) = e j3Ω has the unwrapped phase
Θ(Ω) = 3Ω.
Arg H(Ω) Θ(Ω) = 3Ω
3π 3π
2π 2π
π π
Ω Ω
−π − π2 π π −π − π2 π π
−π 2 −π 2
−2π −2π
−3π −3π
2 translating by −
arg H(Ω) arg H(Ω)
Ω (using bandlimited interpolation if − Ω 6∈ Z).
Therefore, the magnitude response determines how different complex
sinusoids are scaled (in amplitude) by the system.
Similarly, the phase response determines how different complex sinusoids
are translated (i.e., delayed/advanced) by the system.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 499
Magnitude Distortion
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 500
Phase Distortion
Recall that a LTI system H with frequency response H is such that
H{e jΩn }(n) = |H(Ω)| e jΩ[n−τp (Ω)] where τp (Ω) = − arg H(Ω)
Ω .
The function τp is known as the phase delay of the system.
If τp (Ω) = nd (where nd is a constant), the system shifts all complex
sinusoids by the same amount nd .
Since τp (Ω) = nd is equivalent to the (unwrapped) phase response being
of the form arg H(Ω) = −nd Ω (which is a linear function with a zero
constant term), a system with a constant phase delay is said to have
linear phase.
In the case that τp (Ω) = 0, the system is said to have zero phase.
If τp (Ω) is not a constant, different complex sinusoids are shifted by
different amounts, resulting in what is known as phase distortion.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 501
Distortionless Transmission
Consider a LTI system H with input x and output y given by
y(n) = x(n − n0 ),
where n0 is an integer constant.
That is, the output of the system is simply the input delayed by n0 .
This type of behavior is the ideal for which we strive in real-world
communication systems (i.e., the received signal y equals a delayed
version of the transmitted signal x).
Taking the Fourier transform of the preceding equation, we have
H(Ω) = e− jΩn0 .
Since the phase delay of the system is τp (Ω) = − −Ωn
Ω
0
= n0 , the
phase delay is constant and the system has linear phase.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 502
Block Diagram Representations of LTI Systems
Consider a LTI system with input x, output y, and impulse response h, and
let X , Y , and H denote the Fourier transforms of x, y, and h, respectively.
Often, it is convenient to represent such a system in block diagram form in
the frequency domain as shown below.
X Y
H
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 503
Interconnection of LTI Systems
X Y X Y
H1 H2 ≡ H2 H1
X Y
≡ H1 + H2
H2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 504
LTI Systems and Difference Equations
Let h denote the impulse response of the system, and let X , Y , and H
denote the Fourier transforms of x, y, and h, respectively.
One can show that H(Ω) is given by
Y (Ω) ∑M ak (e jΩ )−k ∑M
k=0 ak e
− jkΩ
H(Ω) = = Nk=0 = N .
X(Ω) ∑k=0 bk (e jΩ )−k ∑k=0 bk e− jkΩ
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 506
Duality Between DTFT and CTFS
The DTFT analysis and synthesis equations are, respectively, given by
∞ Z
X(Ω) = ∑
x(k)e− jkΩ and x(n) = 2π 1
X(Ω)e jnΩ dΩ.
2π
k=−∞
The CTFS synthesis and analysis equations are, respectively, given by
∞ Z
xc (t) = a(k)e ∑
jk(2π/T )t
and a(n) = T1
xc (t)e− jn(2π/T )t dt,
T
k=−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 508
Relationship Between DTFT and DFT
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 509
Spectral Sampling Example
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 510
Spectral Sampling Example: N = 4
X̃(k)
4
3
|X(Ω)|
2
1
Ω = π2 k
−π − 3π − 2π − π4 π 2π 3π π
4 4 4 4 4
Magnitude Spectrum
arg X̃(k)
arg X(Ω) π
π
2
Ω = π2 k
−π − 3π − 2π − π4 π 2π 3π π
4 4 − π2 4 4 4
−π
Phase Spectrum
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 511
Spectral Sampling Example: N = 8
X̃(k)
4
3
|X(Ω)|
2
1
Ω = π4 k
−π − 3π − 2π − π4 π 2π 3π π
4 4 4 4 4
Magnitude Spectrum
arg X̃(k)
arg X(Ω) π
π
2
Ω = π4 k
−π − 3π − 2π − π4 π 2π 3π π
4 4 − π2 4 4 4
−π
Phase Spectrum
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 512
Spectral Sampling Example: N = 16
X̃(k)
4
3
|X(Ω)|
2
1
Ω = π8 k
−π − 3π − 2π − π4 π 2π 3π π
4 4 4 4 4
Magnitude Spectrum
arg X̃(k)
arg X(Ω) π
π
2
Ω = π8 k
−π − 3π − 2π − π4 π 2π 3π π
4 4 − π2 4 4 4
−π
Phase Spectrum
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 513
Spectral Sampling Example: N = 64
X̃(k)
4
3
|X(Ω)|
2
1
π
Ω= 32 k
−π − 3π − 2π − π4 π 2π 3π π
4 4 4 4 4
Magnitude Spectrum
arg X̃(k)
arg X(Ω) π
π
2
π
Ω= 32 k
−π − 3π − 2π − π4 π 2π 3π π
4 4 − π2 4 4 4
−π
Phase Spectrum
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 514
Section 11.8
Application: Filtering
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 515
Filtering
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 516
Ideal Lowpass Filter
An ideal lowpass filter eliminates all baseband frequency components
with a frequency whose magnitude is greater than some cutoff frequency,
while leaving the remaining baseband frequency components unaffected.
Such a filter has a frequency response H of the form
(
1 |Ω| ≤ Ωc
H(Ω) =
0 Ωc < |Ω| ≤ π,
H(Ω)
1
Ω
−π −Ωc Ωc π
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 517
Ideal Highpass Filter
An ideal highpass filter eliminates all baseband frequency components
with a frequency whose magnitude is less than some cutoff frequency,
while leaving the remaining baseband frequency components unaffected.
Such a filter has a frequency response H of the form
(
1 Ωc < |Ω| ≤ π
H(Ω) =
0 |Ω| ≤ Ωc ,
H(Ω)
Ω
−π −Ωc Ωc π
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 518
Ideal Bandpass Filter
An ideal bandpass filter eliminates all baseband frequency components
with a frequency whose magnitude does not lie in a particular range, while
leaving the remaining baseband frequency components unaffected.
Such a filter has a frequency response H of the form
(
1 Ωc1 ≤ |Ω| ≤ Ωc2
H(Ω) =
0 |Ω| < Ωc1 or Ωc2 < |Ω| < π,
H(Ω)
Ω
−π −Ωc2 −Ωc1 Ωc1 Ωc2 π
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 519
Part 12
z Transform (ZT)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 520
Motivation Behind the z Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 521
Motivation Behind the z Transform (Continued)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 522
Section 12.1
z Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 523
(Bilateral) z Transform
The (bilateral) z transform of the sequence x, denoted Zx or X , is
defined as
∞
Zx(z) = X(z) = ∑ x(n)z−n .
n=−∞
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 525
Relationship Between Z and Fourier Transforms
Let X and XF denote the z and (DT) Fourier transforms of x, respectively.
The function X(z) evaluated at z = e jΩ (where Ω is real) yields XF (Ω).
That is,
X(e jΩ ) = XF (Ω).
Due to the preceding relationship, the Fourier transform of x is sometimes
written as X(e jΩ ).
The function X(z) evaluated at an arbitrary complex value z = re jΩ (where
r = |z| and Ω = arg z) can also be expressed in terms of a Fourier
transform involving x. In particular, we have
X(re jΩ ) = XF0 (Ω),
where XF0 is the (DT) Fourier transform of x0 (n) = r−n x(n).
So, in general, the z transform of x is the Fourier transform of an
exponentially-weighted version of x.
Due to this weighting, the z transform of a sequence may exist when the
Fourier transform of the same sequence does not.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 526
z Transform Examples
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 527
Section 12.2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 528
Disk
A disk with center 0 and radius r is the set of all complex numbers z
satisfying
|z| < r,
Im
r
Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 529
Annulus
An annulus with center 0, inner radius r0 , and outer radius r1 is the set of
all complex numbers z satisfying
Im
r1
Re
r0
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 530
Circle Exterior
The exterior of a circle with center 0 and radius r is the set of all complex
numbers z satisfying
|z| > r,
Im
r
Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 531
Example: Set Intersection
Im Im
Re Re
3/4 5/4
includes ∞
R1 R2
Im
Re
3/4 5/4
R1 ∩ R2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 532
Example: Scalar Multiple of a Set
Im Im
Re Re
1 2 2 4
R 2R
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 533
Example: Reciprocal of a Set
Im Im
Re Re
3/4 4/3
includes ∞
R R−1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 534
Region of Convergence (ROC)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 535
Property 1: General Form
Re Re Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 536
Property 2: Rational z Transforms
Re Re Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 537
Property 3: Finite-Duration Sequences
Re Re Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 538
Property 4: Right-Sided Sequences
Re Re Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 539
Property 5: Left-Sided Sequences
Re Re Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 540
Property 6: Two-Sided Sequences
Re Re Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 541
Property 7: More on Rational z Transforms
If a sequence x has a rational z transform X (with at least one pole), then:
1 If x is right sided, then the ROC of X is the region outside the circle of
radius equal to the largest magnitude of the poles of X (i.e., outside the
outermost pole), possibly including ∞.
2 If x is left sided, then the ROC of X is the region inside the circle of radius
Re Re Re Re
Thus, we can infer that, if X exists, the ROC can only be of one of the
forms listed above.
For example, the sets shown below would not be valid as ROCs.
Im Im
Re Re
Invalid Invalid
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 543
Section 12.3
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 544
Properties of the z Transform
Property
Initial Value Theorem x(0) = lim X(z)
z→∞
Final Value Theorem lim x(n) = lim[(z − 1)X(z)]
n→∞ z→1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 545
z Transform Pairs
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 546
z Transform Pairs (Continued)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 547
Linearity
ZT ZT
If x1 (n) ←→ X1 (z) with ROC R1 and x2 (n) ←→ X2 (z) with ROC R2 , then
ZT
a1 x1 (n) + a2 x2 (n) ←→ a1 X1 (z) + a2 X2 (z) with ROC R containing R1 ∩ R2 ,
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 548
Translation (Time Shifting)
ZT
If x(n) ←→ X(z) with ROC R, then
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 549
Z-Domain Scaling
ZT
If x(n) ←→ X(z) with ROC R, then
Re Re
r0 r1 |a| r0 |a| r1
R |a| R
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 550
Time Reversal
ZT
If x(n) ←→ X(z) with ROC R, then
ZT
x(−n) ←→ X(1/z) with ROC 1/R.
This is known as the time-reversal property of the z transform.
As illustrated below, the ROC R is reciprocated.
Im Im
Re 1 1 Re
r0 r1 r1 r0
R 1/R
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 551
Upsampling
Define (↑ M)x(n) as
(
x(n/M) n/M is an integer
(↑ M)x(n) =
0 otherwise.
ZT
If x(n) ←→ X(z) with ROC R, then
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 552
Downsampling
ZT
If x(n) ←→ X(z) with ROC R, then
M−1
1
∑ X e− j2πk/M z1/M with ROC RM .
ZT
(↓ M)x(n) ←→ M
k=0
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 553
Conjugation
ZT
If x(n) ←→ X(z) with ROC R, then
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 554
Convolution
ZT ZT
If x1 (n) ←→ X1 (z) with ROC R1 and x2 (n) ←→ X2 (z) with ROC R2 , then
ZT
x1 ∗ x2 (n) ←→ X1 (z)X2 (z) with ROC containing R1 ∩ R2 .
This is known that the convolution (or time-domain convolution)
property of the z transform.
The ROC always contains the intersection but can be larger than the
intersection (if pole-zero cancellation occurs).
Convolution in the time domain becomes multiplication in the z domain.
This can make dealing with LTI systems much easier in the z domain than
in the time domain.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 555
Z-Domain Differentiation
ZT
If x(n) ←→ X(z) with ROC R, then
ZTd
nx(n) ←→ −z dz X(z) with ROC R.
This is known as the z-domain differentiation property of the z
transform.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 556
Differencing
ZT
If x(n) ←→ X(z) with ROC R, then
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 557
Accumulation
ZT
If x(n) ←→ X(z) with ROC R, then
n
z
∑
ZT
x(k) ←→ X(z) for ROC containing R ∩ |z| > 1.
k=−∞ z−1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 558
Initial Value Theorem
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 559
Final Value Theorem
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 560
More z Transform Examples
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 561
Section 12.4
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 562
Finding the Inverse z Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 563
Section 12.5
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 564
System Function of LTI Systems
Consider a LTI system with input x, output y, and impulse response h, and
let X , Y , and H denote the z transforms of x, y, and h, respectively.
Since y(n) = x ∗ h(n), the system is characterized in the z domain by
Y (z) = X(z)H(z).
As a matter of terminology, we refer to H as the system function (or
transfer function) of the system (i.e., the system function is the z
transform of the impulse response).
When viewed in the z domain, a LTI system forms its output by multiplying
its input with its system function.
A LTI system is completely characterized by its system function H .
If the ROC of H includes the unit circle |z| = 1, then H(e jΩ ) is the
frequency response of the LTI system.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 565
Block Diagram Representation of LTI Systems
Consider a LTI system with input x, output y, and impulse response h, and
let X , Y , and H denote the z transforms of x, y, and h, respectively.
Often, it is convenient to represent such a system in block diagram form in
the z domain as shown below.
X Y
H
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 566
Interconnection of LTI Systems
X Y X Y
H1 H2 ≡ H2 H1
X Y
≡ H1 + H2
H2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 567
Causality
outermost pole of H or, if H has no poles, the entire complex plane; and
2 H is proper (i.e., when H(z) is expressed as a ratio of polynomials in z, the
order of the numerator polynomial does not exceed the order of the
denominator polynomial).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 568
BIBO Stability
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 569
Invertibility
H(z)Hinv (z) = 1,
1
Hinv (z) = .
H(z)
Since distinct systems can have identical system functions (but with
differing ROCs), the inverse of a LTI system is not necessarily unique.
In practice, however, we often desire a stable and/or causal system. So,
although multiple inverse systems may exist, we are frequently only
interested in one specific choice of inverse system (due to these
additional constraints of stability and/or causality).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 570
LTI Systems and Difference Equations
Let h denote the impulse response of the system, and let X , Y , and H
denote the z transforms of x, y, and h, respectively.
One can show that H(z) is given by
Y (z) ∑M ak z−k
H(z) = = k=0 .
X(z) ∑Nk=0 bk z−k
Observe that, for a system of the form considered above, the system
function is always rational.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 571
Section 12.6
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 572
Feedback Control Systems
Reference
Input Error Output
+ Controller Plant
−
Sensor
Feedback
Signal
input: desired value of the quantity to be controlled
output: actual value of the quantity to be controlled
error: difference between the desired and actual values
plant: system to be controlled
sensor: device used to measure the actual output
controller: device that monitors the error and changes the input of the
plant with the goal of forcing the error to zero
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 573
Stability Analysis of Feedback Control Systems
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 574
Section 12.7
Unilateral z Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 575
Unilateral z Transform
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 576
Unilateral z Transform (Continued 1)
Z−1 −1
u {Zu {x}}(n) = Zu {Z{xu}}(n)
= Z−1 {Z{xu}}(n)
= x(n)u(n)
(
x(n) n ≥ 0
=
0 otherwise.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 577
Unilateral z Transform (Continued 2)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 578
Properties of the Unilateral z Transform
Property
Initial Value Theorem x(0) = lim X(z)
z→∞
Final Value Theorem lim x(n) = lim[(z − 1)X(z)]
n→∞ z→1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 579
Unilateral z Transform Pairs
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 580
Solving Difference Equations [Using the Unilateral z Transform]
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 581
Part 13
Complex Analysis
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 582
Complex Numbers
z = x + jy,
where x and y are real numbers. The quantities x and y are called the real
part and imaginary part of z, and are denoted as Re z and Im z,
respectively.
The polar form of the complex number z expresses z in the form
where r and θ are real numbers and r ≥ 0. The quantities r and θ are
called the magnitude and argument of z, and are denoted as |z| and
arg z, respectively. [Note: e jθ = cos θ + j sin θ.]
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 583
Complex Numbers (Continued)
Since e jθ = e j(θ+2πk) for all real θ and all integer k, the argument of a
complex number is only uniquely determined to within an additive multiple
of 2π.
The principal argument of a complex number z, denoted Arg z, is the
particular value θ of arg z that satisfies −π < θ ≤ π.
The principal argument of a complex number (excluding zero) is unique.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 584
Geometric Interpretation of Cartesian and Polar Forms
Im Im
z z
y
θ
Re Re
x
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 585
The arctan Function
Im Im
(1, 1)
1 1
π + arctan( −1
−1
) −1
arctan( −1 )
arctan( 11 )
Re Re
−1 1 −1 1
−1 −1
(−1, −1)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 586
The atan2 Function
The angle θ that a vector from the origin to the point (x, y) makes with the
positive x axis is given by θ = atan2(y, x), where
arctan(y/x) x>0
π/2
x = 0 and y > 0
atan2(y, x) , −π/2 x = 0 and y < 0
arctan(y/x) + π x < 0 and y ≥ 0
arctan(y/x) − π x < 0 and y < 0.
The range of the atan2 function is from −π (exclusive) to π (inclusive).
For the complex number z expressed in Cartesian form x + jy,
Arg z = atan2(y, x).
Although the atan2 function is quite useful for computing the principal
argument (or argument) of a complex number, it is not advisable to
memorize the definition of this function. It is better to simply understand
what this function is doing (namely, intelligently applying the arctan
function).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 587
Conversion Between Cartesian and Polar Form
z = x + jy and z = re jθ .
To convert from polar to Cartesian form, we use the following identities:
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 588
Properties of Complex Numbers
z1 + z2 = z2 + z1 and
z1 z2 = z2 z1 .
For complex numbers, addition and multiplication are associative. That is,
for any three complex numbers z1 , z2 , and z3 ,
z1 (z2 + z3 ) = z1 z2 + z1 z3 .
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 589
Conjugation
z∗ = x − jy.
Geometrically, the conjugation operation reflects a point in the complex
plane about the real axis.
The geometric interpretation of the conjugate is illustrated below.
Im
z = x + jy
Re
z∗ = x − jy
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 590
Properties of Conjugation
|z∗ | = |z| ,
arg z∗ = − arg z,
zz∗ = |z|2 ,
Re z = 12 (z + z∗ ), and
1 ∗
Im z = 2 j (z − z ).
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 591
Addition
Cartesian form: Let z1 = x1 + jy1 and z2 = x2 + jy2 . Then,
z1 + z2 = (x1 + jy1 ) + (x2 + jy2 )
= (x1 + x2 ) + j(y1 + y2 ).
That is, to add complex numbers expressed in Cartesian form, we simply
add their real parts and add their imaginary parts.
Polar form: Let z1 = r1 e jθ1 and z2 = r2 e jθ2 . Then,
z1 + z2 = r1 e jθ1 + r2 e jθ2
= (r1 cos θ1 + jr1 sin θ1 ) + (r2 cos θ2 + jr2 sin θ2 )
= (r1 cos θ1 + r2 cos θ2 ) + j(r1 sin θ1 + r2 sin θ2 ).
That is, to add complex numbers expressed in polar form, we first rewrite
them in Cartesian form, and then add their real parts and add their
imaginary parts.
For the purposes of addition, it is easier to work with complex numbers
expressed in Cartesian form.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 592
Multiplication
That is, to multiply two complex numbers expressed in polar form, we use
exponent rules.
For the purposes of multiplication, it is easier to work with complex
numbers expressed in polar form.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 593
Division
Cartesian form: Let z1 = x1 + jy1 and z2 = x2 + jy2 . Then,
z1 z1 z∗2 z1 z∗2 (x1 + jy1 )(x2 − jy2 )
= = 2
=
z2 z2 z2 |z2 |
∗
x22 + y22
x1 x2 − jx1 y2 + jx2 y1 + y1 y2 x1 x2 + y1 y2 + j(x2 y1 − x1 y2 )
= = .
x22 + y22 x22 + y22
That is, to compute the quotient of two complex numbers expressed in
Cartesian form, we convert the problem into one of division by a real
number.
Polar form: Let z1 = r1 e jθ1 and z2 = r2 e jθ2 . Then,
z1 r1 e jθ1 r1
= jθ
= e j(θ1 −θ2 ) .
z2 r2 e 2 r2
That is, to compute the quotient of two complex numbers expressed in
polar form, we use exponent rules.
For the purposes of division, it is easier to work with complex numbers
expressed in polar form.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 594
Properties of the Magnitude and Argument
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 595
Euler’s Relation and De Moivre’s Theorem
e jθ = cos θ + j sin θ.
From Euler’s relation, we can deduce the following useful identities:
cos θ = 12 (e jθ + e− jθ ) and
1 jθ − jθ
sin θ = 2 j (e − e ).
De Moivre’s theorem. For all real θ and all integer n,
n
e jnθ = e jθ .
[Note: This relationship does not necessarily hold for real n.]
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 596
Roots of Complex Numbers
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 597
Quadratic Formula
az2 + bz + c = 0,
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 598
Complex Functions
F(z) = a0 + a1 z + a2 z2 + · · · + an zn ,
a0 + a1 z + a2 z2 + . . . + an zn
F(z) = ,
b0 + b1 z + b2 z2 + . . . + bm zm
where a0 , a1 , . . . , an , b0 , b1 , . . . , bm and z are complex.
Observe that a polynomial function is a special case of a rational function.
Herein, we will mostly focus our attention on polynomial and rational
functions.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 599
Continuity
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 600
Differentiability
A function F is said to be differentiable at a point z = z0 if the limit
F(z)−F(z0 )
F 0 (z0 ) = limz→z0 z−z0
An open disk in the complex plane with center z0 and radius r is the set of
complex numbers z satisfying
|z − z0 | < r,
Im
z0
Re
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 602
Analyticity
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 603
Zeros and Singularities
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 604
Zeros and Poles of a Rational Function
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 605
Part 14
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 606
Motivation for PFEs
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 607
Strictly-Proper Rational Functions
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 608
Section 14.1
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 609
Partial Fraction Expansions (PFEs) [CT and DT Contexts]
am vm + am−1 vm−1 + . . . + a0
F(v) = .
vn + bn−1 vn−1 + . . . + b0
Furthermore, the denominator polynomial D(v) = vn + bn−1 vn−1 + . . . + b0
in the above expression for F(v) can be factored to obtain
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 610
Simple-Pole Case [CT and DT Contexts]
D(v) = (v − p1 )(v − p2 ) · · · (v − pn ),
where the pk are distinct.
In this case, F has a partial fraction expansion of the form
A1 A2 An−1 An
F(v) = + +...+ + ,
v − p1 v − p2 v − pn−1 v − pn
where
Ak = (v − pk )F(v)|v=pk .
Note that the (simple) pole pk contributes a single term to the partial
fraction expansion.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 611
Repeated-Pole Case [CT and DT Contexts]
Suppose that the (rational) function F has at least one repeated pole.
In this case, F has a partial fraction expansion of the form
A1,1 A1,2 A1,q1
F(v) = + +...+
v − p1 (v − p1 )2 (v − p1 )q1
A2,1 A2,q2
+ +...+
v − p2 (v − p2 )q2
AP,1 AP,qP
+...+ +...+ ,
v − pP (v − pP )qP
where
1 h q −` i
d k qk
Ak,` = [(v − p k ) F(v)] .
(qk − `)! dv v=pk
Note that the qk th-order pole pk contributes qk terms to the partial fraction
expansion.
Note that n! = (n)(n − 1)(n − 2) · · · (1) and 0! = 1.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 612
Section 14.2
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 613
Partial Fraction Expansions (PFEs) [DT Context]
Any rational function F can be expressed in the form of
am vm + am−1 vm−1 + . . . + a1 v + a0
F(v) = .
bn vn + bn−1 vn−1 + . . . + b1 v + 1
Furthermore, the denominator polynomial
D(v) = bn vn + bn−1 vn−1 + . . . + b1 v + 1 in the above expression for F(v)
can be factored to obtain
q1 −1 q2 −1 qn
D(v) = (1 − p−1
1 v) (1 − p2 v) · · · (1 − pn v) ,
D(v) = (1 − p−1 −1 −1
1 v)(1 − p2 v) · · · (1 − pn v),
where
Ak = (1 − p−1
k v)F(v) v=pk
.
Note that the (simple) pole pk contributes a single term to the partial
fraction expansion.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 615
Repeated-Pole Case [DT Context]
Suppose that the (rational) function F has at least one repeated pole.
In this case, F has a partial fraction expansion of the form
" #
A1,1 A1,2 A1,q1
F(v) = −1
+ −1 2
+...+
1 − p1 v (1 − p1 v) (1 − p−1
1 v)
q1
" #
A2,1 A2,q2
+ −1
+...+
1 − p2 v (1 − p−1
2 v)
q2
AP,1 AP,qP
+...+ +...+ ,
1 − p−1
P v (1 − p−1P v)
qP
where
1 h q −` i
Ak,` = (−pk )qk −` dv
d k
[(1 − p−1
k v)qk
F(v)] .
(qk − `)! v=pk
Note that the qk th-order pole pk contributes qk terms to the partial fraction
expansion.
Note that n! = (n)(n − 1)(n − 2) · · · (1) and 0! = 1.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 616
Part 15
Miscellany
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 617
Sum of Arithmetic and Geometric Sequences
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 618
Part 16
Epilogue
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 619
Other Courses Offered by the Author of These Lecture
Slides
Computing
For further information about the above courses (including the URLs for
web sites of these courses), please refer to the slides that follow.
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 620
1
0.8
0.6
0.4
0.2
0
80
60 70
60
40 50
40
30
20 20
10
0 0
ECE 486/586:
Multiresolution Signal and Geometry Processing with C++
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 621
SENG 475:
Advanced Programming Techniques for Robust Efficient
Computing (With C++)
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 622
Part 17
References
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 623
Online Resources I
Copyright © 2013–2020 Michael D. Adams Signals and Systems Edition 3.0 625