Understanding Digital Signal Processing Solutions
Understanding Digital Signal Processing Solutions
Understanding Digital
Signal Processing
Third Edition
Richard G. Lyons
Antoine Trux
This work is protected by United States copyright laws and is provided solely for the use
of instructors in teaching their courses and assessing student learning. Dissemination or
sale of any part of this work (including the World Wide Web) will destroy the integrity of
the work and is not permitted. The work and materials from it should never be made
available to students except by instructors using the accompanying text in their classes.
All recipients of this work are expected to abide by these restrictions and to honor the
intended pedagogical purposes and the needs of their instructors who rely on these
materials.
ISBN-13: 978-0-13-218106-8
ISBN-10: 0-13-218106-1
Solution: 1.1
The solution is:
100
π ≈ 4 ⋅ ∑ ( −1) ⋅ 1 .
n
2n + 1
n =0
Solution: 1.2
There are many correct solutions to this problem. An example of a continuous
time-domain signal that has a finite number of amplitude values is a
squarewave, such as the x1(t) shown in Figure S1–2(a). Signal x1(t) has only two
possible amplitude values. Of course, any bi-level pulsed signal, like the x2(t)
shown in Figure S1–2(b) is also a correct solution to this problem. Figure S1–
2(c) shows a continuous signal having only three possible amplitude values.
Solution:
Three possible solutions:
x1(t)
1 ...
(a)
0
0 t
x2(t)
1 ...
(b)
0
0 t
x3(t)
2
(c) 1 ...
0
0 t
Figure S1–2
Solution: 1.3
The code: PI = 2*asin(1.0) correctly defines π
under the assumption that the arcsin's angle argument, "(1.0)",
is measured in radians.
Solution: 1.4
Property of Pearson Education. Not permissible for redistribution.
From one of the Laws of Exponents, we can
write:
xp p–q
xq = x .
Solve the problem by setting q = p, giving us
xp p–p 0
xp = 1 = x = x .
Solution: 1.5
The cosine sequences are as follows:
Solution: 1.6
Solution:
x1(n)
1 ...
1 3 5
(a) 0
0 2 4 6 n
–1
x2(n)
1 ...
2 6
(b) 0
0 1 3 4 5 n
–1
x3(n)
1 ...
(c) 0
0 1 2 3 4 5 6 n
Figure S1–6
Solution: 1.7
x2(n)
1 ...
3
(b) 0
0 1 2 4 5 6 n
–1
x3(n)
1 ...
(c) 0
0 1 2 3 4 5 6 n
Figure S1–7
Solution: 1.8
The desired xshift(n) = x(n+1) sequence is shown in Figure S1–8(b).
Solution:
x(n)
1 ...
2 6
(a) 0
0 1 3 4 5 7 n
–1
xshift(n) = x(n+1)
1 ...
1 5
(b) 0
0 2 3 4 6 7 n
–1
Figure S1–8
Solution: 1.9
In our text, we represent a sinusoidal sequence using the form
m(n) = sin(2πfonts)
Property of Pearson Education. Not permissible for redistribution.
where fo is the sinusoid's frequency measured in Hz. Setting that expression
equal to the problem's m(n) expression, we can write:
2πfonts = 0.8πn.
Recalling the definition that ts = 1/fs, solving the above expression for the
frequency fo, we have our solution of:
0.8πn 0.8
fo = = ⋅ f s = 0.4 ⋅ 2500 = 1000 Hz.
2πnts 2
Solution: 1.10
With N = 6 and n = 9, the computation needed to compute y(9) is
5
y(9) = ∑ x(9–p) = x(9) + x(8) + x(7) + x(6) + x(5) + x(4).
p=0
Solution: 1.11
(a) The block diagram implementing
n
y(n) =
∑
k=n–4
1
5 x(k)
Solution:
x(n) x(n–1) x(n–2) x(n–3) x(n–4)
Delay Delay Delay Delay
y(n)
y(n)
0.2 ...
(b)
0
0 1 2 3 4 5 6 7 8 9 n
Property of Pearson
FigureEducation.
S1–11 Not permissible for redistribution.
(c) Implementing Eq. (P1–1) is preferred over implementing Eq. (P1–2)
because Eq. (P1–1) requires fewer multiplications (lower
computational workload) to compute each y(n) output sample.
Solution: 1.12
On the musical scale, the dimension of the x-axis is time, and the dimension of
the y-axis is frequency.
Solution: 1.13
Using the trigonometric identity:
cos(α+β) + cos(α–β) = 2cos(α)cos(β) (1.13–1)
and the problem's original
x(n) = cos(2πfonts + φ) + cos(2πfonts)
expression, we can write two simultaneous equations as
α + β = 2πfonts + φ (1.13–2)
and
α – β = 2πfonts. (1.13–3)
Solving Eqs. (1.13–2) and (1.13–3) for α and β yields
α = 2πfonts + φ/2, and β = φ/2.
Substituting α and β into Eq. (1.13–1) we write
cos(2πfonts + φ) + cos(2πfonts)
Solution: 1.14
The x = α and y = sin(α) curves are shown in Figure S1–14. There we see that
over the range of roughly α = –0.1π to α = 0.1π the statement "For small α,
sin(α) = α" statement is valid.
–π/2
Figure S1–14
Solution: 1.15
The solutions are:
Solution: 1.16
First we determine how many discrete samples are required to represent a single
cycle of the analog sinewave. That number of samples per analog signal cycle is
found using
fs in samples/sec
samples/cycle = analog sinewave frequency in cycles/sec
100x106 samples/sec
= 25x106 cycles/sec = 4 samples/cycle
x(n)
1 ...
3 7
0
0 1 2 4 5 6 n
–1
Figure S1–16
Solution: 1.17
The proportionality characteristic of a linear system, in the text's Eq. (1–14),
states that if input sequence x(n) yields output y(n),
results in
x(n) y(n),
then a scaled input sequence cx(n), where c is some constant scalar value, yields
a scaled output cy(n),
results in
cx(n) cy(n).
(a) For system ya(n) = x(n–1)/6, the answer is Yes.
For example, if we consider a new x'(n) = 2x(n) input, then the new ya'(n)
output sequence is
ya'(n) = x'(n–1)/6 = 2x(n–1)/6 = 2ya(n).
Solution: 1.18
Decimation is not time-invariant.
An example of this, where yshift(m) ≠ y(m+1), is as follows:
Figure S1–18
Solution: 1.19
We prove the two networks in Figure S1–19 exhibit the commutative property
of linear time-invariant systems as follows:
For the network in Figure S1–19(a) we write output y1(n) as
⎡ B ⎤
y2 (n) = A ⎢ x(n) + y2 (n − 1) ⎥ = Ax(n) + By2 (n − 1)
⎣ A ⎦
which is identical in form to the above y1(n) expression, which is what we set
out to prove.
Delay Delay
A A
y1(n–1) y2(n–1)/A
B B
Property
(a) of Pearson Education. Not
(b) permissible for redistribution.
Figure S1–19
Solution: 1.20
The block diagram solutions to these problems are shown in Figure S1–20.
Solution:
4th-order comb filter
x(n) x(n–1) x(n–2) x(n–3) x(n–4)
Delay Delay Delay Delay
Notice the
(a)
minus sign
–
yC(n)
Delay Delay
A
yI(n–1) yLI(n–1)
(b) (c)
(1-A)
0.5 –0.5 + –
yD(n)
yD(n)
(d) 0.5
x(n) + yD(n)
x(n)
Delay Delay
Delay –
0.5 + – 0.5
yD(n) Delay
Figure S1–20
Solution: 1.21
The impulse response solutions to these problems are shown in Figure S1–21.
y(n) Integrator
1
(b) ...
0
0 1 2 3 4 5 6 7 8 9 n
y(n)
0.5 Leaky integrator
(c) 0.0312
0.25
0.0625
0.125 0.0156 ...
0
0 1 2 3 4 5 6 7 8 9 n
y(n)
Differentiator
0.5
...
(d) 0
0 1 2 3 4 5 6 7 8 9 n
–0.5
Figure S1–21
Solution: 1.22
The step response solutions to this problem are shown in Figure S1–22.
y(n) Integrator
9
7
(b) 5 ...
3
1
0 1 2 3 4 5 6 7 8 9 n
y(n) Differentiator
0.5
(d) ...
0
0 1 2 3 4 5 6 7 8 9 n
Figure S1–22
Solution: 1.23
(a) The original s(t) and the negative of the fundamental frequency
[(4A/π)sin(2πfot)] are shown in Figure S3–23(a). Adding those two
waveforms results in the interesting waveform in Figure S3–23(b) that is the
solution to this problem.
s(t) –4Asin(2πfot)/π
A
... ...
(a)
t=0 t (Time)
–A
Figure S1–23
(b) The operating frequency range of an amplifier needed to exactly double the
ideal s(t) squarewave's peak-peak amplitude would be infinitely wide!
Solution: 1.24
Step 5 is the illegal step because it is an incomplete square root operation.
The square root of q2 is equal to ±q. So following Step 4, Step 5 should have
been:
Solution: 2.1
(a) The photos would incorrectly indicate that the clock's minute
hand is rotating counterclockwise (anti-clockwise).
(b) With the idea of lowpass sampling in mind, because the minute hand rotates
at a frequency of one cycle/hour, to show the true clockwise minute hand
rotation:
Photos must be taken at rate of more than two photos/hour. That is, the
time between photos must be less than 1/2 hour.
Thinking a bit more deeply about this problem, it is true to say that if the
time T between photos is in the range K < T < (K + 0.5), where T is in
hours and K is a positive integer, then the sequence of photos would
indicate correct clockwise minute hand rotation. (The above Part (b)
solution is the case where K = 0.) For example, if T was one hour and
twenty minutes, clockwise minute hand rotation would be seen in the
sequence of photos.
Solution: 2.2
The important missing information is the fs sample rate.
Without knowing fs, we cannot use the x(n) samples to
characterize either the time-domain or frequency-
domain nature of the x(t) signal.
Solution: 2.3
(a) The ts period for a 2 GHz sample rate is:
1 1
ts = f = 2x109 = 0.5x10–9 seconds = 0.5 nanoseconds.
s
Note: 0.5 nanoseconds is a very short period of time—in that time light
travels only 5.9 inches (15 cm).
(b) With the time between samples being ts = 0.5x10–9 seconds, 256x106
samples represents a time interval of:
= cos[2πn/8 + π/7].
Solution: 2.5
The sample period ts must be less than half the sinewave's period of to = 1/fo.
to 1
Thus ts must be less than 2 = 2f , or
o
1
0 < ts < 2f .
o
Solution: 2.6
The answer is N = 2.
We must obtain no less than 2 discrete samples, in time, per analog
sinewave cycle. We verify this for an analog sinewave of frequency fo Hz
using the lowpass Nyquist criterion, fs ≥ 2fo. The periods of fs and 2fo must
satisfy
1 1 to
fs ≤ 2fo , or t s≤
2 , or 2ts ≤ to .
Solution: 2.7
Yes, a continuous sinewave signal whose frequency is fs/2.
For example, a continuous x(t) signal defined by:
x(t) = sin[2π(fs/2)t].
With the continuous t replaced with nts, and recalling that ts = 1/fs, the
discrete samples are:
x(n)
x(t) ts
1 ...
0
0 1 2 3 4 5 6 n
–1
Figure S2–7
Solution: 2.8
Ignoring the days when the Stock Exchange is closed, the ts period
for standard stock market charts is one day, 24 hours.
Solution: 2.9
We can write x(t) in the standard form of:
Equating the first and last angle arguments of the above cosine expressions
yields
nπ/2 = 2πfon/fs,
or
2πfon
fs = = 4fo = 4(2000) = 8000 Hz.
nπ/2
Figure S2–9
Solution: 2.10
= 0.5cos(3800πt) + 0.5cos(4200πt).
The highest-frequency spectral component in x(t) is 2100 Hz. So, the
minimum acceptable fs sample rate is twice that frequency, or
Solution: 2.11
The first step to solving this problem is finding frequency fo in terms of fs.
We can write x(n) in the standard form of:
Equating the first and last angle arguments of the above sin expressions
yields
nπ/4 = 2πfon/fs,
or
fsnπ/4 fs 160
fo = = 8 = 8 = 20 Hz.
2πn
Next, knowing this one possible value for fo (20 Hz), we write the text's Eq.
(2–5) as
x(n) = sin(nπ/4) = sin(2πfonts) = sin(2π(fo + kfs)nts),
reminding us that "aliases" of fo are fo+kfs, where k is any integer. Using
k = 1 and k = 2, we solve the problem by stating that three possible positive
frequency values for fo that would result in sequence x(n) are:
fo = 20 Hz,
fo = 20 + 1·160 = 180 Hz, and
fo = 20 + 2·160 = 340 Hz.
or any of:
fo = 20 + k·160 Hz, for integer k.
Solution: 2.12
By drawing the Xd(f) spectrum in Figure S2–12 showing spectral
replications as dashed lines, the problem solution is the following table.
Figure S2–12
Solution: 2.13
Based on the expression x(t) = sin(2π700t), the frequency of the analog
sinusoidal tone is 700 Hz and its spectrum is shown in Figure S2–13(a).
(a)
–700 0 700 Hz
Figure S2–13
Solution:
|X(f)|, spectrum of x(n)
... ...
(b)
–2000 –1000 –fs/2 0 fs/2 1000 2000 Hz
(–2fs) (–fs) (fs) (2fs)
–700 700
–1300 –300 1700
(700–2fs) (700–fs) (700+fs)
–1700 300 1300
(–700–fs) (–700+fs) (–700+2fs)
|A(f)|
(a)
–70 0 70 Freq
(MHz)
Solution:
|X(f)|
56 56
MHz MHz
(b)
Figure S2–14
Solution: 2.15
As shown in Figure S2–15, the minimum fs sample rate for lowpass
sampling such that no spectral overlap occurs in the frequency range of 2 -
to- 9 kHz in the spectrum of the discrete x(n) samples is:
fs,min = 19 kHz.
Sampled |X(f)|
7 kHz 7 kHz
–1 0 1 2 9 10 17 19 Freq
(fs) (kHz)
Solution: 2.16
The signal's bandwidth (928 Hz) is too wide relative to its center frequency
(711 Hz) to enable bandpass sampling.
|X(f)|
928 Hz
Figure S2–16
Solution: 2.17
(a) The spectrum of the u(t) output of the mixer will be the spectrum of signal
w(t), shown in Figure S2–17(a), translated up in frequency centered at 50
MHz + fLO and a copy of the spectrum of w(t) translated down in frequency
to be centered at 50 – fLO, as shown in Figure S2–17(b). To ensure that the
low-frequency image is centered at 15 MHz, then
fLO = 50 – 15 = 35 MHz.
(a)
–50 0 50 MHz
|U(f)|
(b)
0
–85 –15 15 85 MHz
(–fc) (fc)
|X(f)|
bandwidth B = 5 MHz
(c)
0
–85 –15 15 85 MHz
(–fc) (fc)
Figure S2–17
(b) Analog bandpass filter# 2 has two purposes. First, it filters out (attenuates)
the high-frequency image spectrum that's centered at 85 MHz in U(f).
Second, it reduces the bandwidth of the desired image signal centered at 15
MHz to a bandwidth of 5 MHz, as shown in Figure S2–17(c).
(c) The correct table entries are as follows (don't forget, B = 5 MHz for x(t)):
(d) The text's Eq. (2–11) tells us that to force the sampled spectrum to be
centered at fs/4, then fs must satisfy:
4fc
fs = 2k – 1 ,
As such, the fc center frequencies of the U(f) and X(f) spectra, and the
acceptable fLO local oscillator frequencies, can be list (in MHz) as:
k fc fLO
Property of Pearson Education. Not permissible for redistribution.
1 3 50 – 3 = 47 MHz
2 9 50 – 9 = 41 MHz
3 15 50 – 15 = 35 MHz
4 21 50 – 21 = 29 MHz
... ... ...
Solution: 2.18
If we use the variable f–60 dB to represent the frequency at which the anti-
aliasing filter must have an attenuation value of –60 dB, as shown in Figure
S2–18(a), then
f–60 dB = fs –B Hz.
|H(f)|
Analog filter
2B magnitude
response
(a) 60 dB
0 f–60 dB Freq
|X(m)|
B B
(b) 60 dB
–B 0 fs Freq
fs –2B
f–60 dB = B + (fs –2B) = fs –B
Figure S2–18
Solution: 2.19
The problem solution is the spectrum shown in Figure S2–19.
Figure S2–19
Solution: 2.20
(a) For bandpass sampling, the absolute minimum sampling rate is fs = 2B, as
shown in Figure S2–20(a). So the minimum fc center frequency, in terms of
x(t)'s bandwidth B, that enables bandpass sampling of x(t) is
(b) Here we force the student to find out, on their own, that the AM broadcast
radio band in North America covers the frequency range of roughly 530
kHz to 1710 kHz (depending on who you ask). This AM broadcast-band
spectrum has a bandwidth of BAM = 1180 kHz and a center frequency of
fc,AM = 1120 kHz, as shown in Figure S2–20(b).
From the above Part (a), because fc,AM is less than 1.5BAM (1770
kHz), the full spectrum of the commercial AM broadcast band
cannot by bandpass sampled.
|X(f)|
B B B B
(a)
–fc 0 Freq
fc fs
1.5B
(b)
Figure S2–20
Solution: 2.21
The correct table entries are as follows:
Solution: 2.22
The answer is no, the suggested algorithm does not compute the
absolute minimum fs bandpass sampling rate.
We show this by implementing the algorithm based on the text's Section 2.3
example of a B = 5 MHz-wide bandpass signal center at fc = 20 MHz. Doing
that, we have
4fc + 2B 80 + 10
Z=⎣ 4B ⎦ = ⎣ 20 ⎦ = ⎣4.5⎦ = 4.
Given Z = 4, the algorithm's fs,min,web is
4fc 80
fs,min,web = 2Z –1 = 8 –1 = 11.43 MHz.
From the text's Table 2-1, we see that the minimum fs sample rate (when
m = 3) is 11.25 MHz, which is less than the web site algorithm's fs,min,web.
This verifies that the suggested algorithm does not compute absolute
minimum fs bandpass sampling rates.
Figure S2–22(a) shows the spectrum of the analog bandpass signal. Figure
S2–22(b) and (c) show the spectra of the discrete bandpass signals sampled
at the two different sample rates.
(a)
–25 –20 –15 –10 –5 0 5 10 15 20 25 MHz
(b)
–25 –20 –15 –10 –5 0 5 10 15 20 25 MHz
11.43 11.43
(c)
–25 –20 –15 –10 –5 0 5 10 15 20 25 MHz
11.25 11.25
Figure S2–22
Solution: 3.2
(a) The frequency-domain sample spacing of the DFT is:
fs
DFT sample spacing = ,
N
so
fs 1000
N= = = 22.2222.... samples
DFT sample spacing 45
which is not an integer. It is impossible to have a discrete
signal sequence having a non-integer number of samples.
Solution: 3.3
(a) Because the frequency-domain sample spacing of the DFT's X(m) is:
fs 44100
X (m) sample spacing = = = 1 Hz,
N N
the solution is N = 44100/1 = 44100 samples.
(b) The time duration of the x(n) sequence is the number of samples minus one,
N – 1, times the time period between samples (ts = 1/fs). We show that
notion in Figure S3–3, where a seven-sample w(n) sequence has a time
duration of six ts sample periods.
w(n)
0.5
3 4 5
0
1 2 6 n
–0.5
0
ts = 1/fs
Figure S3–3
Solution: 3.4
(a) The frequency spacing of the X(m) DFT samples is
fs 3000
X (m) sample spacing = = = 6 Hz.
N 500
(b) The highest-frequency spectral component that can be present in x(t) with no
aliasing errors in x(n) is half the fs sample rate, ... in this case that frequency
is
highest freq = fs/2 = 1500 Hz.
(c) The spacing between spectral replications is always equal to the fs sample
rate, ... in this case
the spectral replication spacing is 3000 Hz.
Solution: 3.5
(a) For this x1(n) = 9, 9, 9, 9, 9, 9, 9, 9 time-domain sequence, we use x1(n) and
go through the correlation steps used in the DFT Example 1 in Section 3.1.1.
In computing |X1(m)| in this way, the student should see that the |X1(m)|
samples are all zero when m is greater than zero. In computing |X1(0)|, the
student should use Eq. (3–13') to obtain:
7 7
|X1(0)| = ∑ x1(n) = ∑ 9 = 72.
n=0 n=0
(d) The relationship of the |X2(m)| and |X3(m)| DFT samples is:
|X3(m)| = |X2(m)|.
Solution: 3.6
The positive-frequency spectral energy occurs at m = 3, so the cyclic
frequency of the x(n) sinusoid, with N = 8, is:
mfs 3(4000 Hz)
freq of x(n) in Hz, f = N = 8 = 1500 Hz.
Because x(n) is an integer number of cycles of a sinusoid, from the text's Eq.
(3–17) we find the sinusoid's peak amplitude be Ao = (2Mr)/N or:
2·|X(3)| 2·8
1500 Hz tone peak amplitude, Ao = N = 8 = 2.
The phase of X(3) is:
Imag. part of X(3) 5.657
X(3) phase, φ = tan–1[ Real part of X(3) ] = tan–1[ 5.657 ]
Solution: 3.8
(a) With m = N/2, the DFT equation becomes
N–1
X(N/2) = ∑x(n)[cos(2πn(N/2)/N) – jsin(2πn(N/2)/N)]
n=0
N–1 N–1
= ∑x(n)[cos(πn) – jsin(πn)] = ∑x(n)(–1)n.
n=0 n=0
= 0 + cos(πn)sin(θ) = (–1)nsin(θ).
Finally, applying our input x(n) expression to the above X(N/2) expression,
we have our solution of
N–1 N–1
X(N/2) = ∑(–1)nsin(θ)·(–1)n = ∑(1)2nsin(θ)
n=0 n=0
= N·sin(θ).
X(N/2) = N·sin(0) = 0.
X(N/2) = N·sin(π/2) = N.
Solution: 3.9
(a) The first challenge for the student is to correctly define the x(n) time-domain
expression for a complex sinusoid with magnitude Ao (i.e., x(n) = Aoej2πfnts),
having exactly three cycles over N samples. Substituting 1/fs for ts, and 3fs/N
for f, the x(n) expression becomes:
(b) The x(n) time-domain expression for a real-only sinusoid of peak amplitude
Ao with exactly three cycles over N samples is:
Property of Pearson Education. Not permissible for redistribution.
x(n) = Aosin(2πfnts) = Aosin[2πfn/fs] = Aosin[2π3fsn/(Nfs)]
= Aosin(6πn/N).
Using Euler's identity, sin(α) = (ejα – e–jα)/2j, we write x(n) as:
Ao –jAo
x(n) = 2j (ej6πn/N –e–j6πn/N) = 2 (ej6πn/N –e–j6πn/N) .
∑
N–1
–jAo j6πn/N –j6πn/N –j6πn/N
X(3) = ∑ x(n)e–j6πn/N = 2 (e –e )e
n=0
n=0
Solution: 3.10
The first output sample, X(0), of an N-point DFT is equal to the sum of the
input time-domain samples. If X2(0) = X1(0) + 20, then the sum of x2(n) must
be 20 greater then the sum of x1(n). Therefore
Q = 41.
Solution: 3.12
Given a DFT's X(m) samples, we can compute the x(n) time samples using the
inverse DFT defined as:
1 N–1
x(n) = N ∑ X(m)ej2πnm/N , for n = 0, 1, 2, ..., N–1.
m=0
The first sample of x(n) is found by evaluating the above expression at n = 0, or:
1 N–1 1 N–1
x(0) = N ∑ X(m)ej2π(0)m/N
= N ∑ X(m)e0
m=0 m=0
1 N–1
= N ∑ X(m) .
m=0
N–1
∑ X(m) = N·x(0).
m=0
Solution: 3.13
(a) The solution begins by conjugating both sides of the inverse DFT equation
as:
Property of Pearson Education. Not permissible for redistribution.
1 ⎡ N–1 ⎤* 1 N–1
x*(n) = N ⎢ ∑X(m)ej2πnm/N ⎥ = N ∑ [X(m)ej2πnm/N]* .
⎢ ⎥
⎣m=0 ⎦ m=0
1 N–1 1 N–1
x*(n) = N ∑X(m)*(ej2πnm/N)* = N ∑X(m)*e–j2πnm/N.
m=0 m=0
Next we recognize the right side of the above expression to be in the form of
a forward DFT.
Taking the conjugate of both sides of the above expression yields our
desired expression for x(n),
1⎡ N–1 ⎤
x(n) = N⎢ ∑X(m)*e–j2πnm/N ⎥*
⎢ ⎥
⎣m=0 ⎦
showing how to compute inverse DFTs using a forward DFT process.
Solution: 3.14
Windowing the x(n) time-domain samples will reduce DFT leakage.
However, to minimize spectral leakage as much as possible we force the
analog signal's fo frequency to be one of the DFT's analysis frequencies,
such as in the scenario in the text's Figure 3–10(a). Thus,
where m is an integer is in the range 1 < m < N/2–1 (to satisfy the
Nyquist sampling criterion.) The analog x(t) test signal becomes
x(t) = sin(2π[mfs/N]t).
Solution: 3.15
(a) In the X(m) DFT results the fundamental frequency spectral component is
located at m = 9. We find our Part (a) solution using
Property of Pearson Education. Not permissible for redistribution.
mfs 9(22.255 kHz)
fundamental freq, at X(9), = N = 902 = 222.06 Hz.
(b) The frequency of the highest nonzero spectral component of the guitar's
"A3" note is located at m = 54. As was done above,
mfs 54(22.255 kHz)
highest freq, at X(54), = N = 902 = 1.332 kHz.
Solution: 3.16
(a) The original h1(n) comprises a 0 Hz (DC) component of amplitude 0.5 plus
a single-cycle sinusoid whose peak amplitude is 0.5. The DC component
accounts for the non-zero H1(0) sample in H1(m). The sinusoid component
of h1(n) is what produces the non-zero H1(1) sample in H1(m).
(d) Considering the h1(n), h2(n), and h3(n) sequences, and their |H1(m)| ,
|H2(m)|, and |H3(m)| spectral magnitude samples, the statement that can be
made is:
5
0
0 1 2 3 m 4 5 15 16
1
h (n)
0.8 3
0.6
(c)
0.4
0.2
0
0 5 10 15 20 n 25 30 35 40 47
25
|H3(m)| |H3(0)| = 24,
20
|H3(3)| = 12
15
(d)
10
5
0
0 1 2 3 m 4 5 23 24
Figure S3–16
Solution: 3.17
We repeat the Problem's alternate Hanning window expression here as:
⎛ πn ⎞
whan,alt (n) = sin2⎜ N ⎟, for n = 0, 1, 2, . . ., N–1.
⎝ ⎠
1 –cos(2α)
Using the power relations trigonometry identity: sin2(α) = 2 ,
Solution: 3.18
An N-point DFT provides spectral samples whose cyclic frequency sample
spacingProperty
is fs/N. (fof
s isPearson
the sample rate in Hz of
Education. Notthepermissible
original x(n)for
sequence). A
redistribution.
Q-point DFT provides spectral samples whose frequency sample spacing is
fs/Q. Because Q > N, the DFT spectrum of the zero-padded sequence has
finer spectral granularity (smaller spectral sample spacing).
Solution: 3.19
If the N-point DFT has a sample spacing of fs/N = 100 Hz, then the sample
spacing of a 5N-point DFT would be
Solution: 3.20
The solution to this problem is found using the text's Eq. (3–33). The DFT
processing gain increase (in dB) achieved in using a one million-point DFT
compared to using 100-point DFT is:
⎛ 106 ⎞
SNRN=million = SNRN=100 + 20log10⎜ ⎟
⎝ 102 ⎠.
⎛ 106 ⎞
Improved gain = 20log10⎜ ⎟
⎝ 102 ⎠ = 20log10(100) = 40 dB.
Solution: 3.21
(a) Because N = 16 for X2(m), the desired x2(n) is the 16-point time sequence
shown in Figure S3–21. The issues the student should consider in
determining x2(n) are:
Solution:
0.5 x2(n)
0
x(t)
–0.5
0 2 4 6 n 8 10 12 14 15
Property of PearsonFigure
Education.
S3–21 Not permissible for redistribution.
(b) The x2(n) sequence is a sampled version of the analog x(t) signal
with a sampling rate that is twice the sampling rate of x1(n).
Sequence x2(n) is an interpolated-by-two version of x1(n).
Solution: 3.22
The problem's rectangular expression for S, given as
can be written as
Now if the DFT's x(n) input sequence is real only, then X(0) = a + jb is real
only, and the DFT samples have conjugate symmetry. Thus we can state:
• b = 0,
• c = e, and
• d = –g.
Solution: 3.23
The DTFT of x(n) is a continuous function that cannot be plotted using a
digital computer. The N-point DFT of x(n) produces discrete frequency-
domain samples of the DTFT of x(n), which can be plotted.
Solution: 3.24
(a) The Ximp(ω) spectrum is computed using the DTFT equation as:
∞
Ximp(ω) = ∑ximp(n)e–jωn .
n=–∞
Because only ximp(0) is nonzero, we write:
0
Ximp(ω) = ∑ximp(n)e–jωn = ximp(0)e–jω0 = 1.
n=0
The problem solution is the constant-valued Ximp(ω) curve shown in Figure
S3–24.Property of Pearson Education. Not permissible for redistribution.
Solution:
Ximp(ω)
0.5
0
0 π ω 2π
(fs/2) (fs)
Figure S3–24
(b) The answer to this question reflects on one of the most important principles
of signal analysis. That is, signals that are very narrow in the time domain
have spectra that are very wideband. (And vice versa. Bandlimited signals
that have their time durations expanded yield narrowed spectra.) So
lightening, which is narrow in time with high energy, has a very wideband
noise spectrum. Thus two AM radio receivers tuned to two very different
center frequencies both experience amplitude noise from lightening.
Solution: 3.25
The DTFT of x(n) is:
∞
X(ω) = ∑x(n)e–jωn .
n=–∞
Because only x(0) and x(1) are nonzero, we write:
1
X(ω) = ∑x(n)e–jωn = 1·e–jω(0) + 1·e–jω(1)
n=0
|X(ω)| = 2 + 2cos(ω) .
To aid in drawing |X(ω)|, we evaluate |X(ω)| at five convenient frequency
points as:
Solution:
2
1
0
–π –π/2 0 ω π/2 π
Figure S3–25
Solution: 4.1
(a) There is no difference in the results of performing an N-point DFT
and an N-point FFT. The primary difference is that the FFT requires
far fewer multiplication and addition operations than does a DFT.
(b) The radix-2 FFT restriction is that the FFT size N must be integer
power of two.
Solution: 4.2
Based on the property that the FFT's X(m) frequency-domain sample
spacing is:
fs
X (m) sample spacing = = 1 Hz,
N
then the minimum number of time samples is:
fs 44100
N= = = 44100.
1 Hz 1
N = 65536 samples.
Solution: 4.3
(a) The next highest integer power of two greater than 3800 is 4096. Thus the
number of zero-valued samples that must be appended to the original time-
domain sequence to extend its length to 4096 is
(b) The FFT sample spacing is the discrete signal's fs sample rate (in Hz)
divided by the length of the FFT. That is:
fs
FFT bin spacing = N .
The fs sample rate is the discrete signal sequence's length (in samples)
dividedProperty
by the signal's time duration,
of Pearson or Not permissible for redistribution.
Education.
3800 samples
fs = two seconds = 1900 samples/second (1900 Hz).
Solution: 4.4
There are N2 complex multiplies necessary to compute an N-point DFT.
Thus a 32768-point DFT requires
N 32768
MFFT = 2 log2(N) = 2 log2(32768)
The solution to this problem, the ratio of MDFT over MFFT, is:
MDFT 1.074x109
Ratio = M = 2.46x105 ≈ 4.37x103.
FFT
Solution: 4.5
(a) To satisfy the Nyquist criterion, fs must be greater than
2·50 = 100 Hz.
(b) If the average value of x(n) is 17, we use the property described in Section
3.4 stating the FFT's X(0) sample, which is always real, is the product of the
average value ("DC bias") of the x(n) sequence and the FFT size N. That is
Solution 4.6
Property of Pearson Education. Not permissible for redistribution.
If the spectrum analyzer that can perform a 1024-point FFT in 50
microseconds, its maximum throughput rate would occur if a memory bank
finishes loading 1024 x(n) samples every 50 microseconds. In that situation
the ts time period between x(n) samples from the analog-to-digital converter
is
50x10–6 –8
ts = 1024 = 4.88x10 seconds/sample.
Because the fs sample rate is 1/ts, we know that
1 1
fs = t = 4.88 x 10–8 = 204.8x105 samples/second = 20.48 MHz.
s
In order to satisfy the Nyquist sampling criterion the analog x(t) signal can
have a maximum one-sided bandwidth of half the sampling rate. So our
maximum x(t) bandwidth solution is
fs 20.48 MHz
Bmax = 2 = 2 = 10.24 MHz.
Solution 4.7
When the measured frequency of a spectral component changes by
an amount equal to a change in sample rate, that spectral component
in x(t) has a center frequency greater than both fs/2 (20 kHz) and f's/2
(19 kHz), and aliasing is taking place.
We illustrate this situation in Figure S4–7(a) where the solid curves are the
spectrum of the analog x(t) signal and the dashed curves are spectral
replications in X1(m). The x(t) analog signal has a ±15 kHz spectral noise
component (greater than fs/2) whose negative-frequency spectral replication
resides at +5 kHz.
When the sample rate is reduced to f's = 19 kHz, the spectral
replication of x(t)'s negative 15 kHz spectral noise component now aliases to
4 kHz as presented in Figure S4–7(b).
(a)
–20 –16 –12 –8 0 4 8 12 16 20 kHz
(fs)
–5 5 (fs/2)
Spectral replication
aliased to 4 kHz
|X2(m)| f's = 19 kHz
(b)
–20 –16 –12 –8 –4 0 4 8 12 16 20 kHz
–19 (f's/2) 19
(f's)
Figure S4–7
Solution: 4.8
(a) α = e–j2π(3)/16 = e–j6π/16 = e–j3π/8.
Solution: 4.9
(a) All eight x(n) input samples affect the value of the X(2) output sample. That
is, there is a signal flow path from all x(n) inputs to the X(2) output.
(b) Likewise, all eight x(n) input samples affect the value of the X(5) output
sample.
Solution: 4.10
The signal flow solution using optimized decimation-in-time butterflies is
provided in Figure S4–10(a). Because W40 = 1, the 4-point FFT can also be
drawn as shown in Figure S4–10(b).
x(0) X(0)
x(0) X(0)
x(1) –1 X(2)
(b)
x(2) –1 X(1)
Solution: 4.11
By inspection, in a standard N-point decimation-in-time (DIT) FFT, the
number of unique twiddle factors in the qth stage is 2q. So
R = 2q.
The interval of the twiddle factors' exponents for the qth stage decrease by a
factor of two relative to the interval of the twiddle factors' exponents for the
(q–1)th stage. So P must increase as q increases. By inspection we see that
P = 2q,
allowing us to express the values of the 2q unique twiddle factors in the qth
stage of a general N-point DIT FFT as:
q
kth twiddle factor of qth stage: = WNkN / 2 , for k = 0,1,2,...,2q − 1.
Solution: 4.12
(a) To show the signal flow diagram of a standard 8-point decimation-in-time
FFT with bit-reversed outputs with the butterfly full twiddle factors shown
in rectangular notation, first we show the twiddle factors in complex
exponential notation as in Figure S4–12–I.
Property of Pearson Education. Not permissible for redistribution.
stage 1 stage 2 stage 3
x(0) X(0)
–j2π0/8 e–j2π0/8
e–j2π0/8 e
x(4) X(1)
e–j2π1/8
e–j2π2/8
Figure S4–12–I
x(4) –1 X(1)
–j 0.707–j0.707
x(6) –1 +j X(3)
–0.707–j0.707
Figure S4–12–II
(b) What the student should notice about the twiddle factors in the FFT's 1st and
2nd stages in Figure S4–12–II is that those twiddle factors require no
multiplications. The computations within those two stages are only sign
changes or real/imaginary-part data assignments.
Solution: 4.13
We proceed through this grueling solution using the notation that
Wq = e–j2πq/8.
The samples at nodes A through F inside the FFT structure, in terms of the
x(n) input samples, are
A = x(0) + W0x(4) = x(0) + x(4),
B = x(2) + W0x(6) = x(2) + x(6),
C = x(1) + W0x(5) = x(1) + x(5),
D = x(3) + W0x(7) = x(3) + x(7),
E = A + W4B = x(0) + x(4) + W4x(2) + W4x(6),
F = C + W4D = x(1) + x(5) + W4x(3) + W4x(7).
showing that XDFT(2) = XFFT(2) which is what we set out to prove. (Whew!)
Solution: 4.14
For a 16-point decimation-in-time FFT, the X(m) output index order is the
bit reversal of the x(n) input index order. With in-order inputs, the output
samples will be ordered as shown at the FFT output in Figure S4–14.
Solution:
Correct indices
x(0) X(0)
x(1) X(8)
x(2) X(4)
x(3) X(12)
x(4) X(2)
x(5) 16-point X(10)
x(6) FFT X(6)
x(7) X(14)
x(8) X(1)
x(9) X(9)
x(10) X(5)
x(11) X(13)
x(12) X(3)
x(13) X(11)
x(14) X(7)
x(15) X(15)
e –j2πk/N
–j2πm/N
B e D B e –j2πm/N D
(a) (b)
Figure S4–15
Solution: 4.16
There are many ways to draw a real-valued block diagram of a decimation-
in-time butterfly. One way is shown in Figure S4–16–I. All correct solutions
to this problem must perform the necessary computations shown on the right
side of Figure S4–16–I.
Solution:
Necessary computations
Decimation-in-time butterfly
AR CR = AR + BRcos(θ) + BIsin(θ)
AI CI = AI + BIcos(θ) –BRsin(θ)
DR = AR –BRcos(θ) –BIsin(θ)
–
DI = AI + BRsin(θ) –BIcos(θ)
–
BR
–
BI
cos(θ) –sin(θ)
Solution:
Necessary computations
Decimation-in-time butterfly
AR CR = AR + BRcos(θ) + BIsin(θ)
BR DR = AR –BRcos(θ) –BIsin(θ)
–
cos(θ)
sin(θ)
–sin(θ)
AI CI = AI + BIcos(θ) –BRsin(θ)
BI DI = AI + BRsin(θ) –BIcos(θ)
–
cos(θ)
Figure S4–16-II
Solution: 4.17
(a) We determine the maximum possible decimal value of the real part of the
complex output sample C as follows: The expression for the C output of the
FFT butterfly is
C = Areal + jAimag + e–j2πk/N(Breal + jBimag)
= Areal + jAimag + [cos(2πk/N) –jsin(2πk/N)](Breal + jBimag)
= Areal + cos(2πk/N)Breal + sin(2πk/N)Bimag
+ j[Aimag + cos(2πk/N)Bimag –sin(2πk/N)Breal] .
So the real part of the C output sample, Creal , comprises three terms
Creal = Areal + cos(2πk/N)Breal + sin(2πk/N)Bimag .
Here's the solution to this problem:
If Areal, Breal, and Bimag are each equal to 127 when they enter the butterfly,
and if the twiddle factor angle 2πk/N happens to be π/4 = 45o, then output
Creal can have a decimal value as large as
Creal = 127 + cos(π/4)127 + sin(π/4)127 = 306.61
Property of Pearson Education. Not permissible for redistribution.
which is more than twice the maximum possible decimal value (127) for the
real or imaginary parts of input samples A and B.
(b) Because the maximum possible decimal value for the real part of output
sample C is greater than 255 (which would require 9 binary bits), we require
10 binary bits to store the real part of output sample C.
The real part of C is equally likely to be more negative than a decimal –255
after a single FFT butterfly if A and B are large negative numbers.) We
remind the student that without accounting for, or mitigating, this potential
binary-word growth in fixed-point FFT systems, overflow errors could
render FFT results useless.
Solution: 4.18
The number of butterfly operations in a single 16384-point FFT is:
N 16384
BFFT = 2 log2(N) = 2 log2(16384)
≈ 2x1010 butterflies/second.
A C
–1 A–B
B D
–j2πk/N
e
Figure S4–18
Solution: 5.1
(a) The h(k) impulse response of the system is:
h(k)
1 ...
0
0 1 2 3 4 5 k
Figure S5–1
x(n) y(n)
x(n) y(n)
or
1
h(k) = 1, 0, 0, 0, 0, ...
The student should recall that we can estimate the frequency response of the
system by taking the discrete Fourier transform (DFT) of h(k), the solution
of which has a magnitude of
|H(m)| = 1
For this h(k) impulse response, there are two methods the student might use to
determine the frequency response, i.e., the DFT magnitude samples.
Solution Method# 2: Recalling the text's Eq. (3–43) for the DFT of a
general rectangular time sequence, the magnitude of that expression is:
sin(πmK/N)
|H(m)| = .
sin(πm/N)
Solution: 5.2
(a) Based on the expression x(t) = cos(2π800t), the frequency of the analog
sinusoidal tone is 800 Hz and its spectrum is shown in Figure S5–2(a).
(b) The spectral magnitude of y(n) is shown in Figure S5–2(c), where we show
replications of the filter's frequency magnitude response.
(c) Because x(n)'s ±200 Hz sinusoid is within the filter's unity-gain passband,
(a)
–0.8 0 0.8 kHz
Spectral replications
Spectrum of x(n)
Figure S5–2
Solution: 5.3
The time delay through the linear-phase FIR filter is equal to the filter's
group delay. From the text's description of group delay (measured in
seconds), D/(2fs) where D is the number of filter delay elements, we can
write
Max number of delay elements Dmax
6x10–3 > 2fs = 2f ,
s
or
Dmax < 6x10–3 · 2fs = 6x10–3 · 96x103 = 576.
So the number of filter delay elements must be less than 576 which will be
satisfied so long as:
Solution: 5.4
The problem solution is the network shown in Figure S5–4, where the 'D'
symbol indicates a single delay element.
If input x(n) is a unit impulse sequence, the sequence at node L will be the
hLow(n) impulse response, and the sequence at node H will be the required
hHigh(n) impulse response.
Property of Pearson Education. Not permissible for redistribution.
Solution:
Woofer
L D/A
hLow(k)
converter
x(n)
– H D/A
D D D D D D D
+ converter
Tweeter
D = Delay
Figure S5–4
Solution: 5.5
The answer is simple. First we design the desired lowpass filter
to obtain filter coefficients h(k). Then we merely multiply the
h(k) coefficients by 2, yielding the FIR filter a gain of 2.
Solution: 5.6
The answer to this problem is "Nothing."
The lowpass filter's h(k) coefficients need not be changed in any way. With
the actual x(n) sample rate being 40 kHz, the X(f) spectral plot and the
filter's H(f) magnitude response plot remain unchanged with the exception
that their frequency axis values (Hz) are increased by a factor of two.
Solution: 5.7
(a) The characteristic of the filter's frequency response causing the filter's
output sequence to go to all zeros is:
Hopefully the student realizes that fact, and validates that notion by noticing
the filter input sequence has 8 samples/cycle making its frequency fs/8. An
8-point moving average FIR filter has a frequency magnitude response null
at fs/8.
Solution: 5.8
Abrupt changes in the amplitude of a signal are associated
with high frequency spectral content.
Property of Pearson Education. Not permissible for redistribution.
Likewise, a time-domain signal represented by only small amplitude
changes, as time advances, is a low frequency signal.
Solution: 5.9
(a) Because FIR digital filter performs time-domain convolution, we apply the
text's Eq. (5–29). With a filter impulse response h(k) being of length P = 3
and the filter's input sequence being of length Q = 3,
(b) The nonzero filter output sequence samples are [0.5, 1.5, 2.0, 1.5, 0.5], so
Solution: 5.10
To compute the filter's y(n) output sequence, we convolve the x(n) input
sequence with the h(k) impulse response. Using the text's Eq. (5–6) with a
filter impulse response length of M = 3, our convolution is described by
2
y(n) = ∑ h(k)x(n–k).
k=0
We flip the time order of the filter's x(n) input sequence, pass the flipped
sequence over the h(k) sequence, and perform the sums of products as
shown in Figure S5–10(a).
h(k) y(n)
1.0 2.25
0.75 2.0
0.5 1.75
0.25 1.5
0 1.25
x(n) flipped 0 1 2 k
1.0 1.0
Slide x(n) one
0.75 sample to the right
0.75
0.5 for each new y(n) 0.5
computation.
0.25 0.25
0 0
4 3 2 1 0 n 0 1 2 3 4 5 6 n
(a) (b)
Figure S5–10
Property of Pearson Education. Not permissible for redistribution.
Solution: 5.11
The important topic we learned in Chapter 5 that verifies the spectral
replications of discrete sequences is the Convolution Theorem. The
spectrum of a discrete version of the product y(t) = s(t)x(t) is the
convolution of the individual X(f) and S(f) spectra.
X(f)
S(f)
Convolution
Y(f)
... ...
Solution: 5.12
(a) Recalling the convolution theorem, which states that multiplication in the
time domain is equivalent to convolution in the frequency domain, the Y(m)
spectra is the convolution of X(m) with itself. The convolution of a
rectangular function with itself is a triangular function as shown in Figure
S5–12(b). (The "*" symbol means convolution.)
X(m)
Given:
(a)
–0.4fs 0 0.4fs Freq
Solution:
Y(m) Spectral
Y(m) = X(m) * X(m)
replication
(b)
–fs 0 fs Freq
–0.8fs 0.8fs
X(m)
Given:
Solution: 5.13
(a) Based on the interpretation in Figure S5–13(a), the difference equation for
linear interpolation is:
x(n) – x(n–1) x(n) x(n–1)
y(n) = x(n–1) + 2 = 2 + 2 .
y(n)
x(1)
x(0)
0
1 n
Figure S5–13
(b) The various (equivalent) block diagrams for a linear interpolation system
are shown in Figure S5–13(b).
0.5 0.5
y(n)
(b) y(n)
0.5
0.5 Delay
0.5
y(n)
0.5
Solution: 5.14
The DC gain of filter h1 is equal to the sum of filters' coefficients. That is:
6
Filter h1 DC gain = H1(0) = ∑h1(k) = 99.6.
k=0
Likewise, the DC gain of filter h2 is equal to:
6
Filter h2 DC gain = H2(0) = ∑h2(k) = 103.6.
k=0
So H2(0) and H1(0) differ by a value of four. Each of filter h1's two 25.5-
valued coefficients were changed to Q. Thus 2Q = 2·25.5 + 4, giving us our
solution of:
Q = 25.5 + 2 = 27.5.
Solution: 5.15
The filter's difference equation is:
Due to phase linearity, h(0) = h(4) and h(1) = h(3). So we may rewrite the
difference equations as:
Solution:
x(n–4) x(n–3)
z–1 z–1
x(n) x(n–1)
z–1 z –1
y(n)
Figure S5–15
Solution: 5.16
To derive an equation for the H1(ω) and H2(ω) frequency responses, we
write the discrete-time Fourier transform (DTFT) of the filters' coefficients
as:
∞ 2
H1(ω) = ∑h1(k)e-jkω = ∑h1(k)e-jkω = e–j0ω + 2e–jω + e–j2ω
k=–∞ k=0
= 1 + 2e–jω + e–j2ω.
and
∞ 1
H2(ω) = ∑h2(k)e -jkω
= ∑h2(k)e-jkω = e–j0ω + e–jω = 1 + e–jω.
k=–∞ k=0
Next we set the H1(ω) and H2(ω) expressions equal to each other to
determine frequency ωo as:
H1(ωo) = H2(ωo) → 1 + 2e–jωo + e–j2ωo = 1 + e–jωo.
With algebraic simplification the above equality can be written as:
e–jωo = –e–j2ωo.
or
1 = –e–jωo.
or
e–jωo = –1.
Because e–j±π = –1, the desired ωo frequency is:
ωo = ±π radians/sample
Property
which is theofsolution
Pearson Education.
to our problem. Not permissible for redistribution.
In terms of cyclic frequency, the ωo = ±π frequency is equivalent to ±fs/2 Hz
(half the sample rate). The student could have explicitly solved the 1 = –e–jωo
equation for the desired ωo frequency using natural logarithms as:
ln(1) = ln(–1)·(–jωo)
0 = ±jπ –jωo
ωo = ±π.
This problem deserves additional discussion. The frequency magnitude
responses of the two FIR filters is shown in Figure S5–16(a). Their phase
responses are provided in Figure S5–16(b). At a frequency of ω = 2.095
radians, the two filters have identical frequency magnitude responses.
However, at that frequency their phase responses differ, so their complex
frequency responses are not equal at ω = 2.095 radians.
4
|H1(ω)|
(a) 2
|H2(ω)|
0
–π –π/2 0 π/2 π
ω
ω = 2.095
90
(b) 0 H2(ω) phase
–90
–180
–π –π/2 0 π/2 π
ω
ω = 2.095
|H2(ω)|
0
–20
dB
(c) |H1(ω)|
–40
–60
–π –π/2 0 π/2 π
ω
Figure S5–16
Solution: 5.17
(a) The filter's center coefficient is:
h1 = 1.6617.
(b) The DC gain of the filter is the sum of its coefficients, or:
0
|H(ω)| Noise
–20
dB
–40
–5 –4 –3 –2 –1 0 1 2 3 4 5
Freq (MHz)
3.35 MHz
Property of PearsonFigure
Education.
S5–17 Not permissible for redistribution.
(c) The 3-tap FIR filter has linear phase because its
coefficients are symmetrical, h(k) = [1, h1, 1].
Solution: 5.18
As described in the text, to achieve a linear phase response the
h(k) coefficients of an N-tap tapped-delay line FIR filter must
be either symmetrical,
h(k) = h(N–k–1)
or anti-symmetrical,
h(k) = –h(N–k–1)
where 0≤k≤(N–1)/2 for odd N, and 0≤k≤(N/2)–1 for even N.
Solution: 5.19
A tapped-delay line FIR filter having N = 512 taps requires N – 1 = 511
addition operations to compute each filter output sample.
Solution: 5.20
From Section 5.10 we learned that the group delay, measured in samples, of
a tapped-delay line filter network having symmetrical (or antisymmetrical)
coefficients is found using
D
group delay = G = 2 samples,
where D is the number of unit-delay elements in the filter's delay line. With
the number of FIR filter delay elements, D, being one less than its number
of taps, the group delay of the filter is
255 – 1
group delay = 2 = 127 samples.
Solution: 5.21
Because the input signal's frequency is within the lowpass filter's passband,
the filter's description of "linear-phase" is the key here. That means we have
symmetrical coefficients and the time delay, equal to the filter's group delay,
through the 70-tap filter is independent of frequency. With the number of
FIR filter delay elements, D, being one less than its number of taps, from
Section 5.10 the group delay (measured in seconds) of the filter is
D 70 – 1 69
group delay = G = 2f = 2f = 4x106 = 17.25 microseconds.
s s
where D is the number of unit-delay elements in the filter's delay line. With
the number of FIR filter delay elements, D, being one less than its number
of taps, the group delay of the first filter is
17 – 1
group delay G1 = 2 = 8 samples.
(b) The group delay of the 17-tap H2(f) half-band FIR filter depends on how it
is implemented. Recall, the first and last coefficients of a 17-tap half-band
filter are zero-valued. If the filter is implemented with a tapped-delay line
structure with 17 multipliers, then its group delay will be
17 – 1
group delay G2,17 taps = 2 = 8 samples.
However, that filter can be implemented with only 15 taps by discarding the
first and last zero-valued coefficients with no change in its frequency
response. In that case the group delay of the 15-tap H2(f) half-band FIR filter
will be
15 – 1
group delay G2,15 taps = 2 = 7 samples.
Solution: 5.23
(a) The reverberator's h(n) impulse response is shown in Figure S5–23(a).
Solution:
1
h(n)
(a) 0
–1
0 5 10 15 n 20 25 30
8
Figure S5–23
n=–∞
= e–j0ω + e–j8ω = cos(0)–jsin(0) + cos(8ω)–jsin(8ω) = 1 + cos(8ω) –jsin(8ω).
So |H(ω)| is
|H(ω)| = |2cos(4ω)|.
(c) The plot of |H(ω)| = |2cos(4ω)| is given in Figure S5–23(b). A simple delay
line implementation of a digital reverberator does not have a flat frequency
response, and should not be used in practice.
Solution:
2 |H(ω)|
1
(b)
0
–π –π/2 0 π/4 π/2 π
(–fs/2) ω (fs/2)
(–fs/4) (fs/8) (fs/4)
π/8
(fs/16)
Solution: 5.24
The number of unit-delay elements in the upper path of the parallel-path
filter must be equal to the group delay of the tapped-delay line subfilter in
the bottom path. Because the bottom path subfilter has symmetrical
coefficients and four unit-delay elements, as shown in Figure S5–24, its
group delay is
D 4
G = 2 = 2 = 2 samples.
Delay Delay
y(n)
x(n) x(n–4) –
Delay Delay Delay Delay
w(n)
Figure S5–24
Solution: 5.25
The student should recall that the discrete (sampled) frequency response of
an FIR filter is the DFT of the filter's impulse response (the DFT of the
filter's h(k) coefficients). For an N-tap FIR filter (k = 0, 1, 2, ..., N–1) the N-
sample discrete frequency response is:
N–1
H(m) = ∑ h(k)e–j2πmk/N , for m = 0, 1, 2, ..., N–1.
k=0
Solution: 5.26
The magnitude response solution is:
(N–1)/2
|H(ω)| = |h(0) + 2 ∑ h(k)cos(ωk)| .
k=1
(N–1)/2
= h(0) + ∑h(k)[ejωk + e–jωk]
k=1
Recalling Euler's 2cos(α) = ejα + e–jα identity, we write H(ω) as
(N–1)/2
H(ω) = h(0) + 2 ∑ h(k)cos(ωk).
k=1
Note: The above H(ω) applies to odd-N FIR filters whose symmetrical
coefficients (impulse response) are defined by a symmetrical k index. (With
h(0) being the center coefficient.) Such filters are called "noncausal". For
the more standard, causal, filter indexing where h(0) is the first coefficient,
as shown by example in Figure S5–26, the correct frequency response
expression is
(N–1)/2
Hcausal(ω) = [h(0) + 2 ∑ h(k)cos(ωk)]e–jω(N–1)/2
k=1
with its linear phase shift factor. However, the magnitude of H(ω) and
Hcausal(ω) are equal.
Causal h(k)
0.2
0.1
–0.1
0 1 2 3 4 5 6 k
Figure S5–26
Solution: 5.27
The solution is simple. Just rearrange
Atten
Nfir ≈ 22(f – f )
stop pass
Solution: 5.28
(a) The time-domain difference equation for the filter can be found by way of
inspection of the filter's block diagram, or by redrawing the block diagram
to that shown in Figure S5–28.
b 256 –2b b
y(n)
Figure S5–28
(c) From Figure S5–28, the number of delay elements in the filter's tapped-
delay line is D = 4. So the group delay, measured in samples, is:
D 4
group delay = G = 2 = 2 = 2 samples.
Solution: 5.29
Computing the product of 24 and 13 by way of convolution is shown in
Figure S5–29.
0 0 2 4
3 1
t Sum the
products
2
0 0 2 4
3 1
t Sum the
products
10
0 0 2 4
3 t Sum the
products
12
Figure S5–29
A rough sketch of the Laplace s-plane showing a shaded area that corresponds
to the z-plane shaded band's characteristics is shown in Figure S6–1.
Solution:
jω
(Imag.)
s-plane
(s = σ + jω) fs/4
σ
(Real)
–fs/4
Figure S6–1
Solution: 6.2
The filter transfer functions are:
1
(a) H(z) = 1 + z–2 .
1 + 3z–1 + 2z–2
(b) H(z) = 1 + z–3 .
(a) Order is 2.
(b) Order is 3.
(c) Order is 4.
Solution: 6.4
The filter frequency response equations are those filters' z-domain transfer
functions with variable z replaced with ejω, or:
1
(a) Hpolar(ω) = .
1 + e–j2ω
1 1
(a) HRect(ω) = = .
1 + [cos(2ω) –jsin(2ω)] 1 + cos(2ω) –jsin(2ω)
1 + 3e–jω + 2e–j2ω
(b) Hpolar(ω) = .
1 + e–j3ω
1 + 3[cos(ω) –jsin(ω)] + 2[cos(2ω) –jsin(2ω)]
HRect(ω) =
1 + [cos(3ω) –jsin(3ω)]
Solution: 6.5
(a) One or more poles outside the unit circle means the filter is unstable.
Solution: 6.6
(a) Given H(z) as
z2 + 0.3z + 1
H(z) = z2 + 0.3z + 0.8 ,
(b) The filter is stable, because the filter's poles are within the unit circle. We
verify this as follows: Given the original form of H(z) as
z2 + 0.3z + 1
H(z) = z2 + 0.3z + 0.8 ,
we find the roots of the denominator, which are the locations of the filter's
poles, by first setting the denominator of H(z) to zero as
z2 + 0.3z + 0.8 = 0.
Using the quadratic factorization formula, we factor the above 2nd-order
polynomial to yield
z2 + 0.3z + 0.8 = (z + p1)(z + p2)
(c) The Direct Form I structure for H(z) are shown in Figure S6–6(a).
(d) The Direct Form II structure for H(z) are shown in Figure S6–6(b).
(a) OR
1 –0.8 –0.8
Solution:
Direct Form II Simplified Direct Form II
x(n) y(n) x(n) y(n)
z–1 1 z–1
(b) OR
–0.8 1 –0.8
Figure S6–6
Solution: 6.7
(a) The difference equation describing the filter is
(b) The z-transform of the y(n) difference equation from Part (a) is
Y(z)
H(z) = X(z) = h(0) + h(1)z–1 + h(2)z–2.
(d) The order of the filter is 2 due to the z–2 term in H(z).
Solution: 6.8
(a) The answer is "No."
(b) If we know the H(z) z-domain transfer function equation of a digital filter, to
determine the filter's frequency response we replace z with ejω in H(z),
creating H(ω), and evaluate the H(ω) frequency-domain expression over the
range of –π ≤ ω ≤ π.
Solution: 6.9
The filter block diagrams can be drawn by inspection of Hcas(z), rewritten as:
⎛ 1 − 4 z −1 + 2 z −2 ⎞
H cas ( z ) =
Y ( z) ⎛
=⎜
1 ⎞
( −1
⎟ ⋅ 1 − 4z + 2z
X ( z ) ⎝ 1 − 0.3 z −1 ⎠
−2
)
=⎜ −1 ⎟
.
⎝ 1 − 0.3z ⎠
Or if necessary, the block diagrams can be drawn based on the filter's time-
domain difference equation found by taking the inverse z-transform of Hcas(z).
Doing that, the filter's difference equation is:
y(n) = x(n) – 4x(n–1) + 2x(n–2) + 0.3y(n–1).
Based on either Hcas(z) or y(n), the problem's solutions are the Direct Form I and
Direct Form II block diagrams given in Figure S6–09.
Solution:
Direct Form I Direct Form II
x(n) y(n) x(n) y(n)
–4 0.3 0.3 –4
z–1 z–1
2 2
Figure S6–09
Solution: 6.10
Property of Pearson Education. Not permissible for redistribution.
The h1(k) impulse response of Filter H1 is given in Figure S6–10(a)
showing this filter to be an IIR filter.
The H1 filter's z-plane pole/zero plot is provided in Figure S6–10(b)
showing a pole at z = 1 on the unit circle.
Filter H1
h1(k) 1
Imaginary part
2
pole
1.5
0
1
0.5
–1
0 –1 0 1
0 5 k 10 15 Real part
(a) (b)
Filter H2
h2(k) 1
1 Imaginary part
pole
0
0.5 zero
–1
0 –1 0 1
0 5 k 10 15 Real part
(c) (d)
Figure S6–10
Solution: 6.11
Because the networks contain cascaded subnetworks, we start our solution
by examining the subnetworks. Subnetwork-A is shown in Figure S6–11(a).
Subnetwork-A
Network-A
x(n) w(n) y(n)
(a)
z–1 z–1 z–1
Subnetwork-B
Network-B
x(n) w(n) y(n)
(b)
z–1 z–1 z–1
|HNet-B(ω)| = 1·|HNet-A(ω)|.
Solution: 6.12
We prove that the z-plane pole locations for the two filters are identical by
examining their H(z) transfer functions. For Filter# 1 we write
Y(z) = X(z) + 2Acos(α)Y(z)z–1 – A2Y(z)z–2.
Rearranging the above, we have
Y(z) 1
H1(z) = X(z) = .
1 – 2Acos(α)z–1 + A2z–2
Next, analyzing Filter# 2 we write
U(z) = X(z) – Asin(α)Y(z)z–1 + Acos(α)W(z)z–1
Rearranging the above,
Property we have
of Pearson Education. Not permissible for redistribution.
X(z) – Asin(α)Y(z)z–1
U(z) = . (6.12–1)
1 – Acos(α)z–1
For Filter# 2 we can also write
Y(z) = Asin(α)U(z)z–1 + Acos(α)Y(z)z–1.
Rearranging the above, we have
Y(z)[1 – Acos(α)z–1]
U(z) = . (6.12–2)
Asin(α)z–1
Setting Eqs. (6.12–1) and (6.12–2) equal to each other, we write
Rearranging the above to solve for H2(z) = Y(z)/X(z), and recalling the trig
identity: cos2(α) + sin2(α) = 1, we write
Y(z) Asin(α)z–1
H2(z) = X(z) = .
1 – 2Acos(α)z–1 + A2z–2
Because H1(z) and H2(z) have identical denominator terms their z-plane
pole locations will be identical, which is what we set out to prove.
Solution: 6.13
(a) We can derive the DC bias removal filter's H(z), from Figure S6–13–I(a), by
writing
y(n) = x(n) – x(n–1) + Ay(n–1).
Next we convert that expression to the z-domain and write
Y(z) = X(z) – X(z)z–1 + A(z)Y(z)z–1.
Solving the above expression for Y(z)/X(z) we have our solution of
Y(z) 1 – z–1
H(z) = X(z) = 1 – Az–1 .
Alternatively, the student could find the above H(z) using the text's Eq. (6–25).
Imag.
z-plane
x(n) y(n)
zero at
z=1
z–1 z–1
Real
Unit pole at
circle z=A
–1 A
(a) (b)
Figure S6–13–I
Property
(b) We find of z-plane
the filter's Pearson Education.
zero Not permissible
by setting H(z)'s for redistribution.
numerator equal to zero as:
1 – z–1 = 0.
Solving the above expression for z we see that the filter indeed has a
z-plane zero at z = 1 (at the cyclic frequency of zero Hz), as shown
in S6–13–I(b).
(c) Two equivalent versions of the Direct Form II implementation of the filter
are provided in S6–13–II. This DC bias removal filter is both stable and
often used in practice.
Solution:
Direct Form II Alternate Direct Form II
x(n) y(n) x(n) y(n)
–
z–1 z–1
A –1 A
Figure S6–13–II
Solution: 6.14
(a) This part of the problem is simple because we gave the student the factored
form of the notch filter's transfer function as:
Y(z) 1 –2cos(ωc)z–1 + z–2 (1 –ejωcz–1)(1 –e–jωcz–1)
H(z) = X(z) = = ,
1 –2Rcos(ωc)z–1 + R2z–2 (1 –Rejωcz–1)(1 –Re–jωcz–1)
and R = 0.9. Finding the poles and zeros locations on the z-plane merely
means determining the values of z that make the H(z) numerator and
denominator factors equal to zero. Thus we write
(1 –ejωcz–1) = 0, and (1 –e–jωcz–1) = 0,
specifying two zeros at, z = z0 and z = z1, where
z0 = e+jωc, and z1 = e–jωc. [Both zeros are on the unit circle.]
The location of the filter's poles are found by writing
(1 –Rejωcz–1) = 0, and (1 –Re–jωcz–1) = 0,
specifying two poles, z = p0 and z = p1, where
p0 = 0.9e+jωc, and p1 = 0.9e–jωc.
The pole/zero diagram for this generic notch filter is provided in Figure S6–
14(a).
ωc
(a) –ωc
Real
–1 1
Imag.
z-plane ωc = 0.419
radians
(120 Hz)
ωc
(b) –ωc Real
–1 1
–20
–30
–800 –600 –400 –200 0 200 400 600 800
Freq (Hz) fs/2
120
Figure S6–14
(b) For a signal sample rate of fs = 1800 Hz, to position the notch filter's center
frequency at 120 Hz we set ωc to
120
ωc = 2π . 1800 = 0.419 radians
Solution: 6.15
Property of Pearson Education. Not permissible for redistribution.
To find the locations of the filter's two z-plane zeros, first we set the transfer
function polynomial equal to zero as:
1 + Bz–1 + z–2 = 0.
Next we multiply both sides of the above expression by z2 converting it to
z2 + Bz + 1 = 0
so we can use the following quadratic factorization formula to factor the
polynomial
⎛ b b2 –4ac ⎞ ⎛ b b2 –4ac ⎞
az2 + bz + c = ⎜z + 2a + 2 ⎟ · ⎜z +
2a –
⎟
4a2 ⎠ .
⎝ 4a ⎠ ⎝
Because a = 1, b = B, and c = 1 for our filter polynomial, we can write
⎛ B B2 –4 ⎞ ⎛ B B2 –4 ⎞
z2 + Bz + 1 = ⎜z + 2 + ⎟ · ⎜z +
2 –
⎟
⎝ 4 ⎠ ⎝ 4 ⎠
= (z – z0)(z – z1)
where z0 and z1 are the locations of the filter's two zeros on the z-plane. The
sum of z0 and z1 is
⎛ B B2 –4 ⎞ ⎛ B B2 –4 ⎞
z0 + z 1 = ⎜– 2 – ⎟ ⎜
4 ⎠ + ⎝– 2 +
⎟
⎝ 4 ⎠
or
B B
z0 + z1 = – 2 – 2 = –B
Solution: 6.16
(a) The z-domain transfer function of the filter is, by inspection,:
Y(z) 1 – 6z–1 + 8z–2
H(z) = X(z) = 1 – 2.5z–1 + z–2 .
(b) In preparation for polynomial factoring, to find the filter's poles and zeros,
we multiply H(z) by z2/z2 to the exponents of z positive. Doing so yields
z2 – 6z + 8 (z – 2)(z – 4)
H(z) = z–2 – 2.5z + z–2 = (z – 0.5)(z – 2) .
From the above factored form of H(z) we see that the filter has zeros at z = 1
and z = 4. In addition, the filter has poles at z = 0.5 and z = 2. So the
problem solution is the z-plane pole/zero diagram shown in Figure S6–16(a).
Imaginary part
0.5 cancelation
0
(a)
–0.5
–1
–1 0 1 2 3 4
Real part
Figure S6–16
The block diagram of the simpler equivalent filter (having fewer multipliers)
is shown in Figure S6–16(b).
Solutions:
X(z) Y(z)
z–1
(b)
0.5 –4
Solution: 6.17
The answer is "No." It is not possible to have a real-coefficient filter
whose frequency-domain phase response is of the form:
θ(ω) = C.
The above θ(ω) implies that at zero Hz (ω = 0) we have a non-zero phase
response of C radians, and this situation has no meaning. Stated in different
words, a phase of θ(0) = C implies that with a constant (zero Hz) input signal
applied to the filter, the constant (zero Hz) filter output will be shifted in phase
by C radians. However the notion of phase is not defined (has no meaning) for a
constant-amplitude, DC, signal. All digital filters (having real-valued
coefficients) have a phase response of zero radians at a frequency of zero Hz.
Solution: 6.18
The time-domain difference equations are.
Property of Pearson Education.
w(n) Not permissible for redistribution.
= x(n) + 0.25w(n–1)
and
y(n) = 3[2x(n) + w(n)] + 3w(n–1).
Converting those expressions to the z-domain we write
X(z)
W(z) = X(z) + 0.25W(z)z–1 = 1 – 0.25z–1 (6.18–1)
and
Y(z) = 3[2X(z) + W(z)] + 3W(z)z–1
X(z)
= 6X(z) + 3 1 – 0.25z–1 (1 + z–1). (6.18–2)
or
9 + 1.5z–1
H(z) = 1 – 0.25z–1 .
Solution: 6.19
(a) The network's time-domain difference equation is:
y(n) = x(n) + y(n–1) –Qy(n–1) = x(n) + (1–Q)y(n–1).
To simplify the notation we temporarily replace the (1–Q) factor with the
symbol β, and replace x(n) with the constant D, to write:
y(n) = D + βy(n–1).
At various values of time index n the network's outputs are:
n = 0, y(0) = D,
n = 1, y(1) = D + βD,
n = 2, y(2) = D + β(D + βD) = D + βD + β2D,
n = 3, y(3) = D + β(D + βD + β2D) = D + βD + β2D + β3D,
...
n = N, y(N) = D + βD + β2D + β3D + ... + βND
N
= D + D ∑ βn .
n=1
At time n = 100 we can write:
Property of Pearson Education. Not permissible for redistribution.
100
y(100) = D + D ∑ βn .
n=1
The above y(100) contains a geometric series which can, after reviewing
Appendix B, be written as:
β – β101
y(100) = D + D .
1–β
1–Q D D
y(100) = D + D Q = D + Q –D = Q
Input x(n) = D = 2
10 Imag.
Q = 0.2 z-plane Single pole at
8 zpole = 1–Q
y(n)
6
Q = 0.4
4 Q = 0.6 Real
2 Q = 0.8
0 -1 1
0 5 10 15 20 25
n
(a) (b)
Figure S6–19
1 1
H(z)|z=1 = 1– (1–Q) = Q
(c) To verify stability, the location of the network's poles is found by setting its
H(z) transfer function's denominator equal to zero, and solving for z. Doing
so we have
1– (1–Q)z–1 = 0.
zpole = 1–Q
so the filter is stable when |zpole| = |1–Q| < 1. With Q being a real
number, the |1–Q| < 1 condition means that 0 < Q < 2 guarantees
stability. Thus any Q in the range of 0 < Q ≤ 1 yields a stable system as
shown in Figure S6–19(b).
Solution: 6.20
(a) The pole/zero plot is found, of course, by determining the poles and zeros of
the system's transfer function. We find the desired H(z) transfer function
from the time-domain difference equation of
y(n) = x(n) –y(n–1), or
with a zero at z = 0, and a pole at z = –1. The desired pole/zero plot is shown
in Figure S6–20-I.
Solution:
Imag.
z-plane
pole zero
Real
–1 1
Figure S6–20-I
(c) The time-domain impulse response of the system is shown in Figure S6–20-II.
Solution:
h(n)
1
–1
0 1 3 5 7 9 11 13 15 17 19 21 Time
n
Figure S6–20-II
Solution: 6.21
The solutions to this problem are found using the following steps:
The poles of H4(z) are located at z = 0 + j0.9, and z = 0 – j0.9 as shown in Figure
S6–21. The angles of those poles are ±π/2 radians which correspond to cyclic
frequencies of ±fs/4. Thus the y4(n) filter is associated with the |HA(m)|
frequency magnitude response, and we indicate that by writing:
y4(n) → |HA(m)|.
Imag.
y1(n) filter Pole at
z = 0 +j0.9
Real
–1 1
Figure S6–21
Solution: 6.22
Finding the DC gain of the IIR filter starts with the filter's transfer function:
Solution: 6.23
(a) The allpass filter's Hap(ω) frequency response is the filters' z-domain transfer
function with variable z replaced with ejω. That is
–K + e–jω
Hap(ω) = .
1 – Ke–jω
Next we describe the frequency magnitude squared, |H(ω)|2, as
–K + e–jω [–K + cos(ω)]2 + [–sin(ω)]2
|Hap(ω)|2 = | |2 =
1 – Ke–jω [1 –Kcos(ω)]2 + [Ksin(ω)]2 ]
2 K2 –2Kcos(ω) + 1
|Hap(ω)| = = 1, for all ω.
1 –2Kcos(ω) + K2
(b) The Direct Form I and Direct Form II structures for Hap(z) are shown in
Figure S6–23.
Solution:
Direct Form I Direct Form II
x(n) y(n) x(n) y(n)
K K
Figure S6–23
Alternatively, the student could have replaced z in Hg(z) with ejω to obtain
the Hg(ω) frequency response expression for the filter. Evaluating that
expression at ω = 0 (0 Hz) yields an Hg(ω = 0) of 0.99. Evaluating Hg(ω) at
ω = π (fs/2 Hz) gives an Hg(ω = π) of 0.82. Again showing that the filter is a
lowpass filter.
1 z-plane
Imaginary Part
–1
–1 0 1
Real Part
Figure S6–24
Solution: 6.25
The derivation of the z-domain transfer function of the general 2nd-order
recursive network, in Figure S6–25, proceeds as follows:
b(2) a(2)
Figure S6–25
or
So one solution for the transfer function of the general recursive 2nd-order
network is:
Solution: 6.26
(a) The filter's time-domain difference equation is
y(n) = 2x(n) – e–0.88y(n–1).
Next we convert that expression to the z-domain and write
Property of PearsonY(z)Education. Not
= 2X(z) – e–0.88 permissible
Y(z)z–1
. for redistribution.
Solving the above expression for Y(z)/X(z) we have the z-domain transfer
function of
Y(z) 2
H(z) = X(z) = 1 + e–0.88z–1 .
Replacing H(z)'s variable z with ejω, we write our desired H(ω) frequency
response as:
2
H(ω) = –0.88 –jω .
1+e e
(b) Setting the denominator of H(z) to zero, we find the location of the single z-
plane pole as:
1 + e–0.88zpole–1 = 0.
Solving for zpole, the z-plane pole value is
2 2
Mπ = H(ω)|ω=π = –0.88 –jπ = –0.88
1+e e 1+e [cos(π) –jsin(π)]
1 1
Mπ = 1 + e–0.88[–1 – j0] = 1 – e–0.88 = 3.4175.
Solution: 6.27
(a) The difference equation for the averaging network, shown at the left side of
Figure S6–27(a), is:
x(n) n–1
y(n) = n + n y(n–1).
The right side of Figure S6–27(a) is provided to remind us that
multiplication by 1/n can be interpreted as a division operation.
Implementing the y(n) difference equation requires 2 divisions, 1 multiply,
and 1 addition per output sample. We can do better. (Division operations in
binary arithmetic are computation intensive and are to be avoided whenever
possible.)
z–1 _. z–1
1/n = .n
(a) 1/n _..
n
(n–1)/n (n–1)
Figure S6–27
Solution:
Solution 1
x(n) y(n) x(n) _.. n y(n)
– –
(b) 1/n z–1 z–1
y(n–1) y(n–1)
(n–1) (n–1)
Solution:
Solution 2
x(n) y(n) x(n) _.. n y(n)
(c) – –
1/n z–1 z–1
y(n–1) y(n–1)
Setting the H(z) denominator equal to zero to find the network's single pole
location gives us
1
zpole = 1 – n .
1
zpole = 1 – n [for all equivalent networks].
When n = 1 the networks' poles are located a z = 0 which is inside the unit
circle guaranteeing stability. As n increases, the magnitude of the pole
locations are always less than one revealing that the networks in Figure S6–
27 are always stable.
Solution: 6.28
(a) Given the two poles and single zero on the z-plane, defined by
p0 = 0.25 + j0.25, p1 = 0.25 – j0.25, and z0 = –1,
the associated H(z) transfer function is
G(z – z0) z +1
H(z) = (z – p )(z – p ) = G· (z –0.25 –j0.25)(z –0.25 + j0.25) .
0 1
Next, multiplying the numerator and denominator terms by z–2, we have our
solution of:
z–1 + z–2
H(z) = G · 1 –0.5z–1 + 0.125z–2 .
(b) The Direct Form I block diagram of the H(z) filter is given in Figure S6–
28(a). An alternate, but equivalent, block diagram depiction is given in
Figure S6–28(b).
G 0.125 0.125
(a) (b)
G z–1 G z–1
z–1 z–1
0.125 0.125
(c) (d)
Simplified Structure-III
x(n) y(n)
z–1
z–1 G
z–1
0.5 z–1
(e) 0.125
Figure S6–28
Solution: 6.29
To find the roots of polynomial P, we set P equal to zero, as
Roots of P → z2 + bz + c = 0,
and solve for z. To find the roots of polynomial Q, we set Q equal to zero, as
Roots of Q → Gz2 + Gbz + Gc = 0,
and solve for z. If we multiply both sides of the above 'Roots of Q'
expression by 1/G we have
Roots of Q → z2 + bz + c = 0
Property of Pearson Education. Not permissible for redistribution.
which is equal to the 'Roots of P' expression. Thus the roots of polynomial
Q are equal to the roots of polynomial P, which is what we set out to prove.
Solution: 6.30
(a) Given the pole/zero characteristics of z-plane A, the associated |HA(f)|
magnitude response is shown in Figure S6–30–I(a).
Solution:
Approximate |HA(f)|
∞
Note
(a)
0
–fs/2 –3fs/8 –fs/4 –fs/8 0 fs/8 fs/4 3fs/8 fs/2
Frequency
Figure S6–30–I
Because the poles are directly on the unit circle, the magnitude response is
infinity at f = ±fs/8 Hz.
Solution:
Approximate |HB(f)|
max |Hb(f)|
(< ∞)
(b)
|HB(–fs/2)| = 0 |HB(fs/2)| = 0
0
–fs/2 –3fs/8 –fs/4 –fs/8 0 fs/8 fs/4 3fs/8 fs/2
Frequency
Because the poles are inside the unit circle, the magnitude response is less
than infinity at f = ±fs/8 Hz. Because there is a zero at z = –1, the magnitude
response is zero at f = ±fs/2 Hz.
Imaginary Part
Imaginary Part
0 0
–1 –1
–1 0 1 –1 0 1
Real Part Real Part
(a) (b)
Solution:
Approximate |HC(f)| Passband (main lobe)
max |HC(f)|
Sidelobes Sidelobes
(c)
0
–fs/2 –3fs/8 –fs/4 –fs/8 0 fs/8 fs/4 3fs/8 fs/2
Frequency
Figure S6–30–II
(d) This problem is a bit "tricky". All of the zeros in z-plane D, shown in Figure
S6–30–III(a), affect the |HD(f)| magnitude response but to generate a rough
sketch of |HD(f)| we can ignore all the zeros not lying on the unit circle. As
such, we can use the approximate pole/zero plot shown in Figure S6–30–
III(b) providing the problem solution's |HD(f)| magnitude response shown in
Figure S6–30–III(c).
Imaginary Part
0 0
–1 –1
–1 0 1 –1 0 1
Real Part Real Part
(a) (b)
(c)
0
–fs/2 –3fs/8 –fs/4 –fs/8 0 fs/8 fs/4 3fs/8 fs/2
Frequency
Figure S6–30–III
Solution: 6.31
The solution begins by examining the filter's pole/zero plot. We find the
poles and zeros by factoring the filter's transfer function as
1 – z–2 (1 – z–1)(1 + z–1)
H(z) = 1 – z–1 = 1 – z–1 .
The two factors in the numerator of H(z), when each is set equal to zero,
produce the two z-plane zeros at z = ±1 as shown in Figure S6–31(a). The
single H(z) denominator factor, when set equal to zero, produces the pole at
z = 1. The pole-zero pair at z = 1 cancel each other so the filter's equivalent
pole/zero plot is that shown in Figure S6–31(b). A filter having the pole/zero
plot of Figure S6–31(b), having no poles (no feedback), with a zero at z = –1
will have the following transfer function:
HEquiv(z) = 1 + z–1.
0 0
–1 –1
–1 0 1 –1 0 1
Real part Real part
(a) (b)
Solution:
Simplified HEquiv(z)
x(n) y(n)
(c)
z–1
Solution: 6.32
First, we need expressions for the locations of the two poles in terms of
coefficients ap1 and ap2. The transfer function of the IIR notch filter is
Y(z) 1 + b(1)z–1 + z–2
H(z) = X(z) = 1 –a(1)z–1 –a(2)z–2 .
Setting those H(z) denominator factors equal to zero allows us to specify the
filter's pole locations in terms of coefficients ap1 and ap2. Doing so, we write
(1 + ap1z–1) = 0, and (1 + ap2z–1) = 0 ,
which specifies two poles, p0 and p1, located at
p0 = –ap1, and p1 = –ap2 .
Values –ap1 and –ap2 are complex numbers. OK, from the original problem
specification, we know that
p0 = Rejωn, and p1 = Re–jωn ,
defining coefficients ap1 and ap2 as
ap1 = –Rejωn and ap2 = –Re–jωn .
We can now write H(z) as
1 + b(1)z–1 + z–2
H(z) =
(1 –Rejωnz–1)(1 –Re–jωnz–1)
1 + b(1)z–1 + z–2
H(z) = .
1 –2Rcos(ωn)z–1 + R2z–2
H(z) is now in the standard form allowing us to identify the IIR notch filter's
feedback coefficients. So the solutions to this problem are:
a(1) = 2Rcos(ωn),
a(2) = –R2 .
The IIR notch filter's structure is shown in Figure S6–32.
Property of Pearson Education. Not permissible for redistribution.
IIR biquad, Direct Form II
x(n) y(n)
z–1
2Rcos(ωn) –2cos(ωn)
z–1
–R2
Figure S6–32
Solution: 6.33
(a) We determine the filter's stability by determining the locations of the filter's
H(z) transfer function zeros on the z-plane. H(z) is found from the filter's
difference equation of
y(n) = x(n) + 0.5y(n–1) – 0.81y(n–2).
Taking the z-transform of the above expression, rearranging terms, and
solving for H(z) = Y(z)/X(z), we have
Y(z) 1
H(z) = X(z) = 1 – 0.5z–1 + 0.81z–2 .
Multiplying H(z) by z2/z2, yielding the desired form with z having positive-
only exponents, we have
z2
H(z) = z2 – 0.5z1 + 0.81 .
The first factor is zero when z = p0 = 0.25 + j0.865, and the second pole
location is z = p1 = 0.25 – j0.8646. The poles are each located inside the z-
plane's unit circle at a radius of 0.9 as shown in Figure S6–33, therefore
Radii of poles
are 0.9.
ωr
–ωr Real
1
–1
p1
Figure S6–33
(b) We determine the range of negative values of A for which the filter will be
stable by examining H(z)'s denominator polynomial in terms of the A
coefficient. That expression is
z2 – 0.5z – A.
Using the quadratic factorization formula, we factor the above 2nd-order
polynomial showing the two filter poles located on the z-plane at
1 1
p0 = 0.25 + 16 + A and p1 = 0.25 – 16 + A .
We can write:
A p0 p1 Radii of Poles
0 0.5 0 0.5, 0
–1/32 0.4268 0.0732 0.4268, 0.0732
–1/16 0.25 0.25 0.25
–1/8 0.25 + j0.25 0.25 –j0.25 0.3536
–1/4 0.25 + j0.433 0.25 –j0.433 0.5
–1/2 0.25 + j0.6614 0.25 –j0.6614 0.7071
–3/4 0.25 + j0.8292 0.25 –j0.8292 0.866
–1 0.25 + j0.9682 0.25 –j0.9682 1.0
When A is in the range –1/16 ≤ A ≤ 0, the filter is stable because the two
poles lie on the z-plane's real axis at radii that do not exceed 0.5. When
A < –1/16 the poles move off the real axis and become conjugate poles with
equal radii.
1 1 1
p0 = 0.25 + 16 + A = 0.25 + –j2( 16 + A) = 0.25 + j –A – 16 .
(c) From the above analysis, the magnitude of p0 is unity when |A| is unity, and
we write:
Conditionally stable when: A = –1.
(d) The positive resonant frequency of the filter, in terms of fs, is found from the
angle (shown in Figure S6–33) of the p0 pole when A = –0.81. From the
above p0 = 0.25 + j0.8646 we find angle ωr using
⎛ 0.8646 ⎞
ωres = tan–1⎜ 0.25 ⎟ = 1.289 radians.
⎝ ⎠
Because an angle of 2π radians corresponds to fs Hz, the filter's cyclic
resonant frequency is found using
ωresfs 1.289fs
fres = = 6.283 = 0.205fs Hz.
2π
Solution: 6.34
Given one of the filter's zeros is z0 = 0.5657 + j0.5657 = 0.8ejπ/4, and the
filter has real-valued coefficients, there must be a complex conjugate zero
whose value is z1 = 0.5657 – j0.5657 = 0.8e–jπ/4. Because the filter is linear-
phase, the z0 and z1 zeros must be accompanied by two zeros having the
same angles with magnitudes equal to 1/0.8 = 1.25. Those two zeros are
shown as z2 and z3 in Figure S6–34.
The solutions to this problem are the following four z-domain zeros:
z0 = 0.5657 + j0.5657 = 0.8ejπ/4 (Given)
z1 = 0.5657 – j0.5657 = 0.8e–jπ/4
z2 = 0.8839 + j0.8839 = 1.25ejπ/4
z3 = 0.8839 – j0.8839 = 1.25e–jπ/4.
Imaginary part
z2 = 1.25ejπ/4
0 α
–α
z3 = 1.25e–jπ/4
z1 = 0.8e–jπ/4
–1
–1 0 1 2
Real part
Figure S6–34
Solution: 6.35
(a) The two-stage bandpass IIR filter will detect the B4 musical note (approximately 493
Hz), when the sample rate is 8000 samples/second. We determine this by finding the
pole locations of the two 2nd-order IIR filters by factoring the denominators of their
transfer functions. From the text's Eq. (6–25), the transfer functions are
0.1032 –0.1837z–1 + 0.1032z–2 0.1032z2 –0.1837z + 0.1032
H1(z) = 1 –1.8275z–1 + 0.9834z–2 = z2–1.8275z + 0.9834
and
0.3034 –0.5768z–1 + 0.3034z–2 0.3034z2 –0.5768z + 0.3034
H2(z) = –1
1 –1.8462z + 0.9843z –2 = z2–1.8462z + 0.9843 .
Using the quadratic factorization formula from the text's Eq. (6–15), we
have
First filter denominator polynomial: z2–1.8275z + 0.9834
α1 α2
Real Real
-1 1 -1 1
Figure S6–35
The angles of the filters' poles determine the resonant frequency of each
filter. The angle of the first filter's positive-frequency pole is
α1 = 0.3991 radians.
Multiplying α1 by 8000/2π yields a resonant frequency of 508.1 Hz. The
angle of the second filter's positive-frequency pole is
α2 = 0.3752 radians.
Multiplying α2 by 8000/2π yields a resonant frequency of 477.7 Hz. The
average of these two resonant frequencies is our desired fc resonant
frequency of the cascaded filter. So the solution is:
508.1 + 477.7
fc = 2 = 493 Hz. [Musical note B4]
(b) Yes, the two 2nd-order IIR filters are stable. The first filter's pole
magnitudes are 0.9917, and the second filter's pole magnitudes are 0.9921.
Because these pole magnitudes are less than unity (inside the unit circle),
both filters are stable.
Solution: 6.36
There are two correct solutions to this problem. The transfer function of the IIR
filter is:
B + Bz −1
H ( z) = ,
1 − Az −1
which can be re-written as:
1 + z −1
H ( z) = ⋅B.
1 − Az −1
The above H(z) gives us one solution shown in Figure S6–36(a) requiring only
two multiplies per filter output sample. Because they are linear, the factor B and
the ratio of polynomials factor in H(z) can be swapped in order to give us the
second correct solution, shown in Figure S6–36(b), requiring only two
multiplies Property of Pearson
per filter output sample.Education. Not permissible for redistribution.
Solution:
x(n) y(n) x(n) y(n)
z–1 B B z–1
A A
(a) (b)
Figure S6–36
Solution: 6.37
The transposed structure of a traditional 3-coefficient tapped-delay line FIR
filter is given Figure S6–37.
Solution:
x(n)
Figure S6–37
Solution: 6.38
The transposed structure of Network I is given Figure S6–38(a). An equivalent
depiction is provided in Figure S6–38(b). The transposed structure of Network
II is given Figure S6–38(c).
Solution:
Transposed Transposed Network I
Network I (alternate depiction)
y(n)
x(n)
z–1 A C B
z–1
x(n)
y(n)
C B A
(a) (b)
Figure S6–38
Property of Pearson Education. Not permissible for redistribution.
Solution: 6.39
(a) Starting at time index n = 0, we compute the values of the internal filter
nodes as shown in Table S6–1:
Given the above y(n) impulse response output, the problem solution is
shown in Figure S6–39.
Solution:
10
Quantized y(n) with rounding to nearest q = 0.05,
when x(n) = 1, 0, 0, 0, 0, ...
5
y(6) = 0.05
0
y(7) = –0.05
–5
0 5 n 10 15
Figure S6–39
(b) Comparing the problem's given impulse responses and the above Figure S6–
39, we can state:
The peak-to-peak amplitude of the limit cycles is equal to
twice the value of the rounding precision factor q.
Solution: 6.40
The impulse response of the hCas(k) cascaded combination filter is the
convolution of the h1(k) and h2(k) impulse responses of the two filters.
Solution: 6.41
Property of Pearson Education. Not permissible for redistribution.
(a) The first part of the problem is solved as follows: We define the maximum
and minimum passband gain with the (1 + R) and (1 – R) values,
respectively, as shown in Figure S6–41.
|H(ω)|
1+R
1 P = 2 dB
1–R
Linear
scale
0 fpass Freq
Figure S6–41
(b) The general equation that defines the linear R deviation parameter as a
function of the logarithmic peak-peak passband ripple parameter P is taken
from the above derivation, and is
10P/20 – 1
R(P) = 1 + 10P/20 .
Solution: 6.42
The overall frequency response of cascaded and parallel filters was
discussed in Section 6.8.1 of this chapter. The solution to this problem is:
Property of Pearson Education. Not permissible for redistribution.
Hcombination(ω) = H1(ω)H2(ω)[H3(ω) + H4(ω)]
as indicated in Figure S6–42 below.
H1(ω)H2(ω)[H3(ω) + H4(ω)]
H3(ω)
x(n) y(n)
H1(ω) H2(ω) +
H4(ω)
Figure S6–42
Solution: 6.43
(a) The feedback system is shown in Figure S6–43(a).
x(n) y(n)
+ A(z)
–
B(z)
(a)
C(z) + D(z)
D(z)
(b) (c)
Figure S6–43
By substituting ejω for z in the A(z) and B(z) polynomials the desired
expression for the H(ω) frequency response of the system is:
Y(ω) A(ω)B(ω)
H(ω) = = .
X(ω) 1 + A(ω)B(ω)[C(ω) + D(ω)]
Solution: 6.44
(a) As discussed in the text, the impulse response of a parallel combination of
subfilters is the sum of the individual subfilters' impulse responses. For this
problem we write
hPar(k) = h(k) + [–hHigh(k)].
Solving for h(k), we have
h(k) = hPar(k) + hHigh(k).
Given the following sample values for hPar(k) and hHigh(k),
hHigh(k) = hHigh(0), hHigh(1), ..., hHigh(5), hHigh(6), hHigh(7), ..., hHigh(11), hHigh(12)
h(k) = 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0.
(b) The parallel lowpass filter network, showing the h(k) subfilter being a delay
line whose length is six samples, is shown in Figure S6–44.
Solution:
z–6
x(n) y(n)
–hHigh(k)
Figure S6–44
Solution:
x(n) y(n)
–
G2 z–1 z–M G1
–
z–1 z–M
(a) –
A/G2 G1B
z–1 z–M
1/G2 G1
x(n) y(n)
–
G2 z–1 z–M G1
–
(b) z–1 z–M
–
A G1B
z–1 z–M
G1
Figure S6–45
Solution: 6.46
(a) The impulse invariance Method 2 integrator design proceeds as follows:
Step 3: The analog integrator's H(s) transfer function is already in the form
of individual single pole filters,
1
H (s) = where p1 = 0,
s + p1
so Step 3 can be skipped.
1 1
H ii ( z ) = = .
1 − e0 z −1 1 − z −1
Step 3: Substituting
2 ⎛ 1 − z −1 ⎞
⋅
ts ⎜⎝ 1 + z −1 ⎟⎠
for s in H(s), when ts = 2, yields the desired bilinear transform z-domain
transfer function of
1 1 + z −1
H blt ( z ) = = .
1 − z −1 1 − z −1
1 + z −1
1 1 1
H ii ( z ) = = = , unacceptable gain at fs/2,
z =−1 1− z −1
1 − (−1) 2
1 + z −1 1 −1
H blt ( z ) = =
−1 1 − (−1)
= 0 , acceptable gain at fs/2.
z =−1 1− z
Mag (linear)
|Hblt(f)|
2
0
–fs/2 0 Freq fs/2
Figure S6–46
Solution: 6.47
(a) Given the analog filter's transfer function
1
H (s) = ,
1 + RCs
using the impulse invariance Method 2 process we skip to Method 2 Step 3
of the design process. Because the filter is 1st-order having one path (no
parallel paths) and no partial fraction expansion is necessary, from Method 2
Step 3 we write
A1 1/ RC
H ( s) = = ,
s + p 1 s + 1/ RC
Solution:
Hii filter
x(n) y(n) x(n) y(n)
z–1 1 1 z–1
(a) RC RC
e–1/RC e–1/RC
Figure S6–47
1 1 + z −1
H bt ( z ) = ⋅ .
1 + 2 RC 1 + ⎛ 1 − 2 RC ⎞ z −1
⎜ ⎟
⎝ 1 + 2 RC ⎠
Three equivalent Direct Form II block diagrams of the Hbt(z) digital filter
are shown in Figures S6–47(b) and S6–47(c).
Solution:
Hbt filter
x(n) y(n) x(n) y(n)
z–1 1 z–1 1
1 + 2RC 1 + 2RC
(b)
2RC – 1 1 2RC – 1
1 + 2RC 1 + 2RC 1 + 2RC
Hbt filter
x(n) y(n)
1 z–1
(c) 1 + 2RC
2RC – 1
1 + 2RC
(c) The |Hbt(z)| response has large attenuation at high frequencies because
it has a z-plane zero at z = –1 (infinite attenuation at fs/2 Hz).
Solution: 6.48
Given the analog filter's Laplace-domain transfer function of
s
H ( s) = ,
s +
Property of Pearsonωo Education. Not permissible for redistribution.
using the bilinear transform we substitute
2 ⎛ 1 – z–1 ⎞
ts ⎜⎝ 1 + z–1 ⎟⎠
for s in H(s) yielding the z-domain transfer function
2 1 − z −1
⋅
ts 1 + z −1
H ( z) = .
2 1 − z −1
⋅ + ωo
ts 1+ z −1
Multiplying H(z)'s numerator and denominator by ts·(1 + z–1), we have
H ( z) =
( )
2 1 − z −1
.
2 (1 − z −1 ) + ts (1 + z −1 ) ω o
H ( z) =
1
1 + 0.5ts ωo (
⋅ 1 − z −1
.
)
1 − 0.5ts ωo −1
1− z
1 + 0.5ts ωo
0.761 − 0.761z −1
H ( z) = .
1 − 0.522z −1
Solution: 6.49
Given the analog filter's transfer function
5
H (s) = ,
s( s − 0.8)
using the bilinear transform we substitute
2 ⎛ 1 – z–1 ⎞ –1
⎛1–z ⎞
⎜ –1 ⎟ = a ⎜
ts ⎝ 1 + z ⎠ –1 ⎟, where a = 2/ts,
⎝1+z ⎠
for s in H(s) yielding the z-domain transfer function
5
= 1 – –1
1 – –1 –1 .
⎛ z ⎞⎛ z ⎞ ⎛ 0.8a(1 – z ) ⎞
a2⎜ 1 + z–1 ⎟ ⎜ 1 + z–1 ⎟ – ⎜ 1 + z–1 ⎟
⎝ ⎠⎝ ⎠ ⎝ ⎠
5 + 10z–1 + 5z–2
= (a2 – 0.8a) –2a2z–1 + (a2 + 0.8a)z–2 .
5 + 10z–1 + 5z–2
H(z) = 3.9984x106 –8x106z–1 + 4.0016x106z–2 .
Solution: 6.50
(a) Given that the analog lowpass filter's cutoff frequency is fa = 3.8 kHz, the
digital filter's fd cutoff frequency is found by mapping fa to fd using
YFd(z)
HFd(z) = X(z) = 1– z–1.
The differentiator's frequency magnitude response is the square root of the sum
of HFd(ω)'s real and imaginary parts squared, or
4 – 4cos(ω) 1 – cos(ω)
|HFd(ω)| = 2 – 2cos(ω) = 2 =2 2 .
YCd(z)
HCd(z) = X(z) = 0.5 – 0.5z–2.
The differentiator's frequency magnitude response is the square root of the sum
of HCd(ω)'s real and imaginary parts squared, or
2 2
|HCd(ω)| = [0.5
Property – 0.5cos(2ω)]
of Pearson + 0.25sin
Education. Not (2ω)
permissible for redistribution.
= 0.25 – 0.5cos(2ω) + 0.25cos2(2ω) + 0.25sin2(2ω) .
1 – cos(2ω)
|HCd(ω)| = 0.5 – 0.5cos(2ω) = 2 .
Solution: 7.2
Two block diagrams of a central-difference differentiator having only one
multiplier are shown in Figure S7–2.
Solution:
Central-difference differentiators
x(n) x(n)
z–1 z–1 z–1 z–1
– –
0.5
ycd(n) ycd(n)
0.5
Figure S7–2
Solution: 7.3
Rocky would be incorrect because the web page
statement is not true!
Solution Method# 1: The impulse response of two cascaded first-difference
differentiators is a first-difference differentiator's impulse response
[hFd(n) = [1,–1]) convolved with itself, or
hFd,cascaded(n) = [1,–2,1].
= 1 – 2e–jω + e–j2ω.
Next, obtained in a similar manner, the central-difference differentiator's
frequency response is
HCd(ω) = 0.5 – 0.5e–j2ω.
Because HFd,cascaded(ω) ≠ HCd(ω),
again, the original Internet statement is not true.
Solution: 7.4
We can estimate the acceleration of the motor shaft by first computing the
derivative of the Apos(n) shaft position signal to obtain the Avel(n) motor shaft
velocity signal as shown in Figure S7–4–I(a). Next, we compute the derivative
of Avel(n) to obtain the desired Aacc(n) acceleration signal. Thus our Acceleration
Measurement Network is the cascade of two digital differentiators as shown in
Figure S7–4(a).
Due to the high-frequency noise in the Apos(n) signal, we eliminate first-
difference differentiators as possible solutions. Due to the restriction of performing
no more than one multiplication per Apos(n) input sample both the text's Lanczos
and wideband differentiators are eliminated, leaving only central-difference
differentiators as possible solutions as shown in Figure S7–4(b).
Apos(n)
z–1 z–1 z–1 z–1
(b) + – + –
Aacc(n)
Avel(n)
0.5 0.5
Figure S7–4
Property of Pearson Education. Not permissible for redistribution.
The two multipliers in Figure S7–4(b) can be combined into one multiplier giving
us the two equivalent correct solutions shown in Figures S7–4(c) and (d).
Solutions:
Apos(n)
z–1 z–1 z–1 z–1
(c) + – + –
Aacc(n)
Avel(n)
0.25
Apos(n)
z–1 z–1 z–1 z–1
(d) + – + –
Aacc(n)
Avel(n)
0.25
Solution: 7.5
The student should use the text's Eq. (7–14),
(–1)k
h(k)ωc=π = k ,
to solve this wideband differentiator design problem. Doing so, with –3≤k≤3 and k ≠ 0,
using a table would produce
k hωc=π(k)
–3 (–1)–3
–3 = 0.333
–2 (–1)–2
–2 = –0.5
–1 (–1)–1
–1 = 1
0 0 (set to 0 by
design)
1 (–1)1
1 = –1
2 (–1)2
2 = 0.5
3 (–1)3
3 = –0.333
Note: If the student uses the text's more general Eq. (7–13),
ωccos(ωck) sin(ωck)
hgen(k) = – ,
πk πk2
that equation reduces to
πcos(πk) sin(πk) cos(πk)
hgen(k) = – 2 =
πk πk k
yielding the desired coefficients of
hgen(k) = 0.333, –0.5, 1, 0, –1, 0.5, and –0.333.
(Again, because N = 7 is odd we must set the center coefficient, hgen(3), to
zero.)
Solution: 7.6
From Section 5.10 we learned that the group delay, measured in samples, of a
tapped-delay line filter network having antisymmetrical coefficients is found
using
D
Gdiff = 2 samples,
where D is the number of unit-delay elements in the filter's delay line. The block
diagram (structure) of the ydiff(n) differentiator is shown in Figure S7–6(a).
Redrawing that figure to show the individual unit-delay elements (z–1) gives us
Figure S7–6(b) showing the differentiator to have D = 6 unit-delay elements.
Thus the solution to this problem is
6
group delay Gdiff = 2 = 3 samples, which
is an integer number of samples.
x(n)
z–1 z–1 z–1 z–1 z–1 z–1
–1/16 1/16
(b) –
ydiff(n)
So the real part of HRe(ω) is equal to 0.5 for all ω, which validates the
problem's statement.
Solution: 7.8
(a) Given the x(n) input, shown in Figure S7–8(a), the integrator's y(n) output
sequence is that shown in Figure S7–8(b).
x(n)
1 ...
(a) 3 5 7
0
0 1 2 4 6 8 n
Property
–1 of Pearson Education. Not permissible for redistribution.
Solution:
y(n)
(b) 1
...
0
0 1 2 3 4 5 6 7 8 n
Figure S7–8
(b) Comparing the peak-peak amplitudes of x(n) input and y(n) output
sequences we see that, for an input sinusoid whose frequency is fs/2 Hz,
the integrator has an amplitude loss by a factor of 0.5.
(c) The text's equation for the frequency response of the integrator is:
1
H Re (ω) = .
1 − e − jω
|HRe(ω)|
(c)
0.5
0
0 π 2π ω
(fs/2) (fs)
Solution: 7.9
Using a trapezoidal rule integrator to estimate the area under the x(t) curve will
produce integrator output samples equal to the areas of the shaded trapezoids
shown in Figure S7–9. Those integrator output samples (areas of the shaded
trapezoids) are:
x(n)
Continuous x(t)
A(1) A(3)
A(2)
0
0 1 2 3 n
(ts) (2ts) (3ts)
Figure S7–9
Solution: 7.10
(a) From Eq. (6–105) in the text's Section 6.11, to perform a bilinear
transformation we substitute:
2 1 – z–1
s = t 1 + z–1
s
Substituting e–jω for z in the above HBilin(z) yields the desired bilinear
transform-designed integrator's frequency response of
ts 1 + e–jω
HBilin(ω) = 2 .
1 –e–jω
Solution: 7.11
Solving for the roots of the four transfer functions' denominators, we find
that the four integrators have z-domain pole locations of:
Solution: 7.12
(a) The block diagram of the FIR matched filter is shown in Figure S7–12.
y(n)
Figure S7–12
Solution: 7.13
Given the problem's xs(n) signal-of-interest, and assuming the filter's
coefficients are a time-reversed version of xs(n), the maximum filter output
sample is the sum of the xs(n) samples squared. The algebra expression is:
Maximum y(n) = N–1xs(k)2 .
∑
k=0
Solution: 7.14
Using x1(n) as our signal-of-interest in a matched filter, the filter's output
sequence is:
y1(n) = [9, 30, 43, 30, 9].
The matched filter's y1(n) output SNR is:
Solution: 7.16
(a) The passband width of the lowpass prototype FIR filter is M = 5 times the
lowpass IFIR filter's desired passband width (4 kHz), or
Prototype filter passband width = 5·4 kHz = 20 kHz
as shown in Figure S7–16(a).
(b) Because the shaping subfilter's passband images are centered at integer
multiples of fs/M = fs/5 Hz,
the shaping subfilter will have two passband
images residing between zero and fs/2 Hz
as shown in Figure S7–16(b).
|Hp(f)|
(a)
0 20 fs/2 fs
kHz
|Hsh(f)| M=5
(b)
0 fs/5 2fs/5 3fs/5 4fs/5 fs
fs/2
Figure S7–16
Solution: 7.17
(a) The desired IFIR filter's transition bandwidth is 61.2 kHz minus 60 kHz, or
1.2 kHz. Thus the normalized transition region bandwidth, and the
normalized passband width, are:
ftrans = 1.2x103/3x106 = 0.0004, and fpass = 60x103/3x106 = 0.02.
Using those normalized ftrans and fpass frequencies in the text's Figure 7–
18(a), as shown in Figure S7–17 , yields an
optimum shaping subfilter expansion factor of M = 16.
(c) Using those normalized ftrans and fpass frequencies in the text's Figure 7–
18(b), the computational reduction we expect from the IFIR filter is
roughly
92%.
10 0.04
0.06
5 0.08
Figure S7–17
Solution: 7.18
(a) The time-domain difference equation of the complex multiplier output is
Solution:
+
vR(n) wR(n)
–
vI(n) w(n) =
wR(n) + jwI(n)
v(n) =
vR(n) + jvI(n) wI(n)
cos(ω) sin(ω)
e jω = cos(ω) + jsin(ω)
Figure S7–18
Property of Pearson Education. Not permissible for redistribution.
(c) A single complex multiply requires four real multiplies and two real
additions.
Solution: 7.19
(a) From the first comb filter we see that N = 8 and r = 1. The pole/zero
plots of the cascaded comb and the resonators are given in Figures S7–19–
I(a), (b) and (c). With their pole–zero cancellations, the desired pole/zero
plot solution is shown in Figure S7–19–I(d).
Imaginary part
Imaginary part
Two
0.5 zeros 0.5
2 2 0 2
0
Imaginary part
Figure S7–19–I
(c) This part of the problem is a bit tricky. For the three stages of the FSF to
have equal gains, the gain factors following the second and third resonators
must be twice the value of the gain factor following the first (k = 0, or zero
Hz) resonator. However all three gain factors, 1/16, are equal so the second
and third stages have half the gain as the first stage. Thus the FSF's
frequency magnitude response looks roughly as that shown in Figure S7–
19–II.
Figure S7–19–II
Solution: 7.20
The cascaded (series combination) comb filters are shown below in
Figure S7–20–I(a).
Solution:
Cascaded-comb impulse response
1
h(n)
(b) 0
–1
0 2 4 6 8 10 12 14 16
n
Figure S7–20–I
(b) Yes, the cascaded-comb combination filter is linear phase because the
nonzero samples of the h(n) impulse response are symmetrical as shown in
Figure S7–20–I(b). (The cascade of two linear-phase filters will always
have linear phase.)
(c) With damping factor r = 0.9, the cascaded combination of comb filters in
Figure S7–20–II(a) will not have linear phase because the nonzero samples
Property of Pearson Education. Not permissible for redistribution.
of their h(n) impulse response are not symmetrical. This is shown in Figure
S7–20–II(b).
Figure S7–20–II
Solution: 7.21
The purpose of this problem is to test the student's understanding of the
text's Figure 7–50. If we expand a portion of Figure 7–50, shown below as
Figure S7–21, we see that
0
attenuation (dB)
Min stopband
Figure S7–21
Solution: 7.22
Property of Pearson Education. Not permissible for redistribution.
Yes, it is possible to design an FSF that achieves 85 dB of stopband
attenuation.
From the text's Figure 7–50, shown here as Figure S7–22, we see that
Type-IV FSFs with three transition band coefficients can easily achieve 85
dB of stopband attenuation.
0
256 No transition coefficients R = passband
(R ≈ 1.6 dB)
Min stopband attenuation (dB)
peak-peak ripple
–20
64 32 16
One transition coefficient
(R ≈ 0.7 dB)
–40 128
256 64 32
–60 16 Two transition coefficients
128
(R ≈ 0.35 dB)
32
Three transition
–80 256 64 16 coefficients
–85 dB (R ≈ 0.16 dB)
–100 256 128 64
32
16
–120
0 0.05fs 0.1fs 0.15fs 0.2fs 0.25fs
Transition bandwidth
Figure S7–22
Solution: 7.23
(a) From Table H–4 (for odd N), in Appendix H, we use the values N = 23 and
BW = 3 to find the transition coefficient values of
H(3) = T1 = 0.64635467,
H(5) = T3 = 0.0077623.
(b) From the same row of Table H–4 as was used above, we can expect
an FSF stopband attenuation of roughly 95 dB.
(c) The structure (block diagram) of the six-section Type-IV FSF is shown in
Figure S7–23 below.
k=0
0.5
k=1
x(n) – y(n)
–
k=2
z–23 z–2 0.64635467 –
k=3
–r23 –r2 0.16260027
k=4
0.0077623
k=5
Figure S7–23
Instructor: The key here is to see that the student filled in all the details such as
the comb filters' coefficients, the 0.5 gain factor for the k = 0 (DC) section,
and the correct plus and minus signs for the final summation.
Solution: 7.24
There are two ways to solve this problem; a DFT analysis method, and a z-
domain transfer function analysis method. Both solution methods are provided
below.
(a) DFT Analysis Method: Given this problem, we expect the student to begin
their solution by drawing a block diagram like that shown below in either
Figure S7–24–I(a) or Figure S7–24–I(b).
1
H0(z) =
1 –e j2π0/16 z–1
H(0)/16
1
H1(z) =
1 –e j2π1/16 z–1
x(n) H(1)/16 X2(n)
1 –z–16 1
H2(z) =
1 –e j2π2/16 z–1
(a)
.. .. H(2)/16
. . X2(n)X2(n)*
1
H15(z) = |X2(n)|2
j2π15/16 –1
1 –e z
H(15)/16
z–1
H(0)/16
x(n) e j2π0/16
X2(n)
e j2π2/16
. .
. .
. .
z–1 H(15)/16
e j2π15/16
Figure S7–24–I
The FSF's H(k) gain terms at the outputs of the resonators are the DFT of
the FIR filter's coefficients in the problem's Figure P7–24–II(a). So a key
part of the solution is determining the values for the H(k)/16 gain factors in
the above Figure S7–24–I. Based on the text's Figure 7–20 discussion that
the FSF H(k) gain terms are the DFT of the equivalent nonrecursive FIR
filter's coefficients, we find the desired H(k) gain terms by performing the
DFT of the problem's Figure P7–24–II(a) filter coefficients when m = 2.
N–1 N–1
= ∑e –j2π[2N–2–i(2–k)]/N
= ∑e–j2π2N/Nej2π2/Nej2πi(k–2)/N .
i=0 i=0
Because e–j2π2N/N = 1, we can write the FSF's H(k) gain terms as
N–1
H(k) = ej2π(2)/N ∑e–j2πi(k–2)/N .
i=0
The summation in the above H(k) is a geometric series that can be
converted to a closed form expression (See Appendix B) as
N–1
1 – e–j2π(k–2)
H(k) = e j2π(2)/N
∑e –j2π(k–2)i/N
=e j2π(2)/N
1 – e–j2π(k–2)/N
i=0
sin[π(k–2)]
= e–jπ(k–2)(N–1)/N ej2π(2)/N .
sin[π(k–2)/N]
Substituting various values of k in H(k), to obtain the FSF's gain terms,
reveals that the sin[π(k–2)] terms are zero for all k except k = 2. So we are
left with a single-resonator FSF whose H(2) gain term is
sin[π(k–2)]
H(2) = ej2π(2)/N .
sin[π(k–2)/N]
When k = 2, the phase of H(2) is ej2π(2)/N = ej4π/16 but its magnitude is 0/0
(indeterminate). Using L'Hopital's Rule on
sin[π(k–2)]
|H(2)| = .
sin[π(k–2)/N]
yields
d{sin[π(k–2)]}/dk cos[π(k–2)] d[π(k–2)]/dk
|H(2)| = = ·
d{sin[π(k–2)/N]}/dk cos[π(k–2)/N] d[π(k–2)/N]/dk
cos[π(k–2)] πk
= · .
cos[π(k–2)/N] πk/N
When k = 2, the above |H(2)| becomes
cos(0) 2π
|H(2)| = cos(0) · = N = 16.
2π/N
e jπ/4
Figure S7–24–II
= 15e–j2π(2)(15–k)/16X(z)z–k
∑ .
k=0
k=0 k=0
= e–j2π(2)15/16 15ej2π(2)k/16z–k
∑ .
k=0
Using the material in Appendix B, we can write the above geometric series
in closed form as
1 – (ej2π(2)/16z–1)16 1 – ej2π(2)z–16
HFSF(z) = e–j2π(2)15/16 j2π(2)/16 –1 = e–j2π(2)15/16
1–e z 1 – ejπ/4z–1
1
= (1 – z–16) ejπ/4
1 – ejπ/4z–1
whose implementation is that shown in S7–24–II(b).
Note: The Figure S7–24–II(b) solution to this problem is called a "sliding
DFT", and is well known in the field of DSP. A detailed discussion of the
sliding DFT is provided in Section 13.18.
(b) There are two ways to solve this part of the problem. First, we could use our
knowledge of DFT bin center frequencies. If the fs sample rate of the x(n)
input is 200 kHz, the center frequency of the m = 2 DFT bin is
The second way to find the center frequency of the m = 2 DFT bin is to find
the resonant frequency of the single-section FSF resonator in Figure S7–
24–II(b) by determining the angle of the resonator's z-plane pole location.
From the text's Eq. (7–40), the resonator's Hres(z) transfer function is
1
Hres(z) = .
1 – ejπ/4z–1
Setting Hres(z)'s denominator equal to zero yields
1 – ejπ/4zpole–1 = 0
or
zpole = ejπ/4
defining the resonator's z-plane pole location to be on the unit circle at an
angle of π/4 radians. Because a resonant frequency of 2π radians
corresponds to fs Hz, the filter's cyclic resonant frequency (and the desired
single-bin DFT center frequency) is found using
π . fs 200 kHz
fr = f(m=2 center) = 4 = 8 = 25 kHz.
2π
Solution: 8.2
(a) Rectangular notation proof: Given a complex number C = a +jb,
multiplication by j is:
jC = ja – b = –b + ja.
Numbers C = a +jb and jC = –b + ja are shown in Figure S8–2.
Imag
jC a
θ3
θ2
b C
θ1
–b a Real
Figure S8–2
that anglesofθPearson
Notice Property 1 and θ2 are
Education.
equal as: Not permissible for redistribution.
⎛b⎞
θ1 = tan–1⎜ a ⎟ = θ2.
⎝ ⎠
The amount of rotation is angle θ3 minus angle θ1 or:
⎛π ⎞ π π
θ3 – θ1 = ⎜ 2 + θ2⎟ – θ1 = 2 + θ1 – θ1 = 2 = 90o.
⎝ ⎠
Polar notation proof: Given a complex number C = Mejθ, and using the text's
Eq. (8–12), j = ejπ/2, multiplication by j is:
jC = Mejθejπ/2 = Mej(θ + π/2).
Thus the angle of jC is θ + π/2 radians, or θ + 90o, which is what we set out
to prove. Again, using polar notation yields a simpler solution.
Solution: 8.3
(a) With a rectangular form of C as C = Mcos(φ) + jMsin(φ):
CC* = [Mcos(φ) + jMsin(φ)][Mcos(φ) – jMsin(φ)]
Solution: 8.4
If the sum of a complex number C plus its reciprocal (C + 1/C) is
real only, then the magnitude of C must be unity. That is, |C| = 1.
R I
= (R + R2 + I2 ) + j(I – R2 + I2 ).
leading to
I
I = R2 + I2 .
Solution: 8.5
This is a trick question. It is not valid to use Qefficient to compute the
desired Q value! In general, the real part of the ratio of two complex
numbers is not equal to the ratio of their individual real parts.
We show this as by writing the original ratio of complex numbers:
Ca Ra + jIa
Cb = Rb + jIb .
Using the principles given in Appendix A's Eq. (A–20) , we write:
Ca (RaRb + IaIb) + j(RbIa – RaIb)
Cb = Rb2 + Ib2 .
The correct value for our desired Q is the real part of the above ratio, or:
⎡ Ca ⎤ RaRb + IaIb
Q = real part of ⎢ C ⎥ = R 2 + I 2 .
⎣ b⎦ b b
Solution: 8.6
Property of Pearson Education. Not permissible for redistribution.
The solution proceeds as follows:
cos(α + β) = Re{ej(α + β)} = Re{ejαejβ}
= Re{[cos(α) + jsin(α)]·[cos(β) + jsin(β)]}
= Re{cos(α)cos(β) + jcos(α)sin(β) + jsin(α)cos(β) – sin(α)sin(β)}
= cos(α)cos(β) – sin(α)sin(β).
Solution: 8.7
Given the notation for the spectrum of a continuous cosine wave in the text's Figure
8–10(a), shown here as Figure S8–7(a), we describe the problem's given x(t) signal,
in the time domain, as shown in S8–7(b). Thus the problem solution is:
Imag
Real
–fo
cos(2πfot) =
e j2πfot + e–j2πfot
0
2 2
(a) fo
Freq
Imag
–100 Real
0 100
(b)
Freq (Hz)
e j(2π100t +– π) e –j(2π100t +– π)
2
+
2 – π)
= cos(2π100t +
Figure S8–7
Solution: 8.8
To add the two complex exponentials, we convert them to rectangular form for
addition, and convert the sum back to polar form, as follows:
(
Ae j (ωt +α) + Be j (ωt +β) = e jωt Ae jα + Be jβ )
= e jωt [ A cos(α) + B cos(β) + j ( A sin(α) + B sin(β) )]
⎡ j ⋅ tan −1 ⎡⎢
A sin(α) + B sin(β) ⎤
⎣ A cos(α) + B cos(β) ⎥⎦
⎤
⋅ ⎢ [ A cos(α) + B cos(β) ] + [ A sin(α) + B sin(β) ] ⋅ e
jωt 2 2
=e ⎥
⎢⎣ ⎥⎦
Property of Pearson Education. Not permissible for redistribution.
{
j ωt + tan −1 ⎡⎢
A sin(α) + B sin(β) ⎤
}
⎣ A cos(α) + B cos(β) ⎥⎦
[ A cos(α) + B cos(β) ] + [ A sin(α) + B sin(β) ] ⋅e
2 2
= .
and
⎡ A sin(α) + B sin(β) ⎤
θ = ωt + tan −1 ⎢ ⎥.
⎣ A cos(α) + B cos(β) ⎦
Solution: 8.9
The proof proceeds as follows:
sin(α)
tan(α) = = sin(α) · [1/cos(α)]
cos(α)
ejα – e–jα 2
= · jα
2j e + e–jα
ejα – e–jα
= .
j(ejα + e–jα)
Solution: 8.10
Given two complex numbers, C1 = a + jb and C2 = c + jd, either rectangular
or polar notation may be used to prove |C1|·|C2| = |C1C2|.
= [(a2 + bEducation.
Property of Pearson
2
) (c2 + d2)]1/2Not permissible(above).
= ProdOfMag for redistribution.
which is what we set out to prove. (Whew!)
Polar notation proof: Using C1 = M1ejθ and C2 = M2ejφ, the product of the
magnitudes is:
Hopefully the student is learning that polar notation is often easier to use
than rectangular notation for solving complex-valued algebra problems.
(M1ejθM2ejφ)* = M1e–jθM2e–jφ
(M1M2ej(θ+φ))* = M1M2e–j(θ+φ)
M1M2e–j(θ+φ) = M1M2e–j(θ+φ)
which is what we set out to prove.
Solution: 8.12
Because the signals' magnitudes are |c1(t)| = 1, and |c2(t)| = 0.75, the equation for
the sum of the two signals, having exponents measured in radians, is:
Solution: 8.13
Converting to polar notation makes the proof quite simple. That is,
Here we learn that polar notation is often easier to work with than
rectangular notation in the algebra of quadrature processing.
(c) The real part of q(n) begins with a unity-valued sample (when n = 0) and is
a cosine wave sequence decreasing in amplitude as time index n increases.
(d) The imaginary part of q(n) begins with a zero-valued sample (when n = 0)
and is a sinewave sequence decreasing in amplitude as time index n
increases.
Solution:
1
q(n) q(2)
0.8
q(1)
0.6
q(3)
0.4
Imaginary
0.2
(a) 0
q(4) q(0)
-0.2
q(5) q(7)
-0.4
-0.6 q(6)
-0.8
-1
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
Real
2
1.5
Imag Part
1
0.5
0
(b)
-0.5
-1
-2
0 2
2 1
Time 4 0 Real Part
6
8 -1
Figure S8–14
Solution: 8.15
Using Euler's equation, the cosine sequence x(n) may be written in complex
exponential form as
ej[2π(–1000)nts + π/4] e–j[2π(–1000)nts + π/4]
x(n) = 2 + 2 .
Solution: 8.16
The text's Figure 8–11 gives us the answer, but we can also use algebra to arrive at
the solution. If the highest frequency spectral component of |XC(ω)| is ωo
radians/sample, shown by the bold arrow in Figure S8–16(a), in the discrete time
domain we can express that spectral component as:
Figure S8–16
M·ejωon/fs M·e–jωon/fs
M[cos(ωon/fs) = 2 + 2 .
Solution: 8.17
The complex constant ejπ/2 is equal to the j operator,
ejπ/2 = j.
So we may write y(n) as
y(n) = j·x(n) = j·[xi(n) + jxq(n)] = –xq(n) + jxi(n)
leading to the simple problem's block diagram solution shown in Figure S8–
17.
Solution:
xi(n) yi(n)
–1
xq(n) yq(n)
Figure S8–17
Solution: 8.18
Using the function-product trigonometric identity:
cos(α – β) cos(α + β)
sin(α)sin(β) = 2 – 2
we write:
cos(2πft – 2πft – θ) – cos(2πft + 2πft + θ)
sin(2πft)sin(2πft + θ) = 2 .
Property
Because cos(–θ) of Pearson
= cos(θ), we Education. Not permissible
write our solution as: for redistribution.
cos(θ) cos(2π(2f)t+θ)
sin(2πft)sin(2πft + θ) = 2 – 2 .
1
0
sin(2πft)
–1
1 The average of
θ = π/4 sin(2πft)sin(2πft + π/4)
0 equals
sin(2πft + θ) cos(π/4)/2 = 0.354.
–1
1
θ = π/4
0
sin(2πft)sin(2πft + θ)
–1
0 0.2 0.4 0.6 0.8 1
t
Figure S8–18
Solution: 8.19
To determine the constants A, B, and C, we start with
x(n) = 5ej(2πn/20 + π/4) + 3ej(2πn/20 + π/6)
and factor out the ej2πn/20 term yielding
x(n) = ej2πn/20(5ejπ/4 + 3ejπ/6).
Next we convert the complex exponentials inside the parenthesis to
rectangular form as
x(n) = ej2πn/20[5cos(π/4) + j5sin(π/4) + 3cos(π/6) + j3sin(π/6)].
The terms inside the brackets are constants which we evaluate as
x(n) = ej2πn/20[3.54 + j3.54 + 2.598 + j1.5] = ej2πn/20[6.13 + j5.04]
Property of Pearson Education. Not permissible for redistribution.
= ej2πn/20[7.94ej0.687] = 7.94ej(2πn/20 + 0.687) .
So our desired A, B, and C constants are:
Solution: 8.20
(a) The real part of x(t) is 1 – cos(2π1t), as shown in Figure S8–20(a).
4 – 4cos(2π1t) 1 – cos(2π1t)
= 2 =2 2 .
|x(t)| = |2sin(πt)|.
(a) 1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
1
Imaginary part of x(t)
(b) 0
–1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
2
|x(t)|
(c) 1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Time (seconds)
Figure S8–20
Solution: 8.21
(a) The θ2 arctangent algorithm has the lowest average error magnitude.
(c) The required number of processor clock cycles for the three algorithms are:
Because of its lower computational overhead and its lower average error,
the θ3 arctangent algorithm should be chosen over the θ1 algorithm.
Solution: 8.22
Given that the algorithm's error is
X
E(X) = tan-1(X) – 1 + AX2 radians,
we can find the value Xmax error where the E(X) function is at its maximum by
taking the derivative of E(X) with respect to X, setting that derivative equal
to zero, and solving for Xmax error. Knowing Xmax error, we can compute the
desired angle θ associated with Xmax error. We do this as follows, Given the
above E(X), its derivative with respect to X is
d(1 + AX2)
(1 + AX2) –
dE(X) 1 dX 1 1 + AX2 – 2AX2
dX = 1 + X2 – 2 2
(1 + AX )
= 1 + X2 – 1 + 2AX2 + A2X4 .
With A = 0.28125,
–(0.84375 –1)
X2 = 0.36035 = 0.4336.
So X is
X = ± 0.4336 = ±0.6585
giving the two values for X at which the error is maximized. Our solution is
the arctangent of X = ±0.6585, or
Property of Pearson Education. Not permissible for redistribution.
θmax error = tan–1(±0.6585) = ±0.5823 radians = ±33.36o.
0.4
Error (degrees)
0.2
–0.2
–0.4
–40 –20 0 20 40
True angle θ (degrees)
–33.36o 33.36o
Figure S8–22
Solution: 8.23
In our text, we have represented a complex exponential sequence using the form
m(n) = ej2πfonts
where fo is the complex exponential's frequency measured in Hz. Setting that
expression equal to the problem's m(n) expression, we may write:
m(n) = ej2πfonts = ej0.8πn
Setting the above angle arguments equal to each other, we write:
2πfonts = 0.8πn.
Recalling the definition that ts = 1/fs, solving the above expression for the
frequency fo, we have our solution of:
0.8πn 0.8
fo = = ⋅ f s = 0.4 ⋅ 5000 = 2000 Hz.
2πnts 2
Solution: 8.24
(a)
xi(t) = cos[2π(fc+10)t]cos(2πfct).
= cos(2π10t)/2 + cos[2π(2fc+10)t]/2.
= sin(2π10t)/2 –sin[2π(2fc+10)t]/2.
(c) The minus sign on the second term for xq(t) represents a 180o (π radians)
phase shift between the low and high frequency components in xq(t).
(d) After lowpass filtering (removing the high frequency term in xi(t)),
(e) After lowpass filtering (removing the high frequency term in xq(t)),
Solution: 8.25
The first challenge for the student is to correctly define the xc(n) time-
domain expression for a discrete complex sinusoid, whose frequency is –fc
Hz, used to translate (mix) the x(n) = cos(2πfonts) signal down in frequency
by fc Hz. That discrete complex sinusoid expression is:
xc(n) = e–j2πfcnts.
Using Euler's identity, cos(α) = (ejα + e–jα)/2, we write x(n) as:
1
x(n) = 2 (ej2πfonts + e–j2πfonts).
The problem solution is the xp(n) product of x(n) and xc(n) expressed as:
1
xp(n) = 2 (ej2πfonts + e–j2πfonts)e–j2πfcnts
1
= 2 (ej2π(fo–fc)nts + e–j2π(fo+fc)nts)
Property of Pearson Education. Not permissible for redistribution.
containing spectral components at frequencies 2π(fo–fc) and –2π(fo+fc)
radians/second.
(b) With the spectrum of the original x(n) being that shown in Figure S8–25(a),
the solution is the two spectral components of xp(n) as shown in Figure S8–
25(b).
Original X(m)
1/2 1/2
–fo 0 fo Freq
(Hz)
Solution:
Xp(m)
1/2 1/2
–fo 0 fo Freq
–fo–fc fo–fc (Hz)
[–(fo+fc)]
Figure S8–25
(Complex signals need not have spectral magnitude symmetry around the
zero Hz point as do real signals.)
Solution: 8.26
Given the expression
⎡A A(1+ε) ⎤ –jωot
+ ⎢2 +
⎣ 2 ·[cos(α) – jsin(α)]⎥⎦e
⎡ –Aε Aα ⎤ ⎡ Aα ⎤
mimp(t) ≈ ⎢ 2 – j 2 ⎥ ejωot + ⎢A – j 2 ⎥ e–jωot
⎣ ⎦ ⎣ ⎦
Factoring A from within the square roots and dividing both sides by A we have:
ε2 α2 –3 α2
4 + 4 = 10 1+ 4 .
ε2 α2 –6 ⎡ α2 ⎤
+ = 10 ⎢ 1 + ⎥
4 4 ⎣ 4 ⎦.
Finally, because α<<1, (1+α2/4) ≈ 1 and we write the solution to this problem
as:
ε2 + α2 ≈ 4x10–6
telling us that the sum of errors ε2 and α2 must not be greater
than 4x10–6.
Solution: 8.27
(a) The original incorrect filter structure, repeated here in Figure S8–27(a), has
incorrect subtractions, and incorrect signs for the "B" coefficients.
Property of Pearson Education. Not permissible for redistribution.
Incorrect signs
yi(n)
xi(n) –
z–1
–B
(a)
+B
xq(n)
z–1 yq(n)
–
Figure S8–27
To prove this, we determine the real-valued difference equation for the quadrature
bandpass filter as:
yi(n) + jyq(n) = [xi(n) + jxq(n)] + [A + jB][yi(n–1) + jyq(n–1)]
–B B B
xq(n) xq(n) – xq(n) –
z–1 yq(n) z–1 yq(n) z–1 yq(n)
–
A A –A
(b) (c) (d)
Y(z) 1
H(z) = X(z) = .
1 –ej2πfr/fsz–1
(d) The filter is only conditionally stable. Setting the denominator of H(z) equal
to zero and solving for z gives:
1 –ej2πfr/fsz–1 = 0, or
z = ej2πfr/fs.
This results in a single filter pole at z = ej2πfr/fs directly on the z-plane's unit
circle, and that's what we call "conditionally stable".
1
H(f) = –j2π(f – fr/fs) ← Polar form
1 –e
1
= . ← Rectangular form
1 –cos[2π(f – fr/fs)] +jsin[2π(f – fr/fs)]
+
i(n) i'(n)
–
q(n)
(a) q'(n)
Solution:
i(n) i'(n)
xc'(n) = i'(n) + jq'(n)
(b) q'(n)
cos(2πfcnts)
sin(2πfcnts) up-conversion
–sin(2πfcnts) down-conversion
Figure S8–28
Solution:
+
i(n) i'(n)
–
q(n)
(c)
cos(2πfcnts) sin(2πfcnts) up-conversion
–sin(2πfcnts) down-conversion
Solution: 8.29
First we need to know the difference equation for our complex resonator,
which is:
Property of Pearson Education. Not permissible for redistribution.
y(n) = x(n) + ejωry(n–1)
+ j[cos(ωr)yimag(n–1) + sin(ωr)yreal(n–1)].
So we build the complex oscillator by implementing
yreal(n) = x(n) + cos(ωr)yreal(n–1) – sin(ωr)yimag(n–1)
and
yimag(n) = cos(ωr)yimag(n–1) + sin(ωr)yreal(n–1).
Solution:
Complex digital resonator (oscillator)
x(n) yreal(n)
z–1
yimag(n)
z–1
+
–
sin(ωr) cos(ωr)
Figure S8–29
Solution: 8.30
(a) Because the real-coefficient lowpass filter's transfer function is
Hreal(z) = 1 + z–1,
that filter has a z-plane zero at z = –1. We want a filter having a z-plane zero
at z = e–jπ/2, so the transfer function of our desired complex filter is:
HProperty of Pearson
cmplx(z) = 1 – e
–jπ/2 –1 Education. Not permissible for redistribution.
z .
Because –e–jπ/2 = ejπ/2, an alternate (and equivalent) Hcmplx(z) transfer
function is:
(b) Block diagrams of the two equivalent Hcmplx(z) filters are shown in Figures
S8–30(a) and (b).
Solutions:
Hcmplx(z) Alternate Hcmplx(z)
x(n) y(n) x(n) y(n)
–
z–1 z–1
e–jπ/2 ejπ/2
(a) (b)
Figure S8–30
Solution: 8.31
The operation of the communications system is shown in Figure S8–31.
I(t)cos(ωct)
I(t)
cos(ωct)
–sin(ωct) +
I(t)cos(ωct)
–Q(t)sin(ωct) – Q(t)sin(ωct)
Q(t)
Antenna
I(t)
LPF
I(t) [cos(0t) + cos(2ωct)]/2 2
–Q(t) [sin(2ωct) - sin(0t)]/2
cos(ωct)
Q(t)
LPF
–I(t)[sin(2ωct) – sin(0t)]/2 2
–sin(ωct) + Q(t) [cos(0t) – cos(2ωct)]/2
Figure S8–31
In the demodulator the lowpass filters (LPF) eliminate the signal components
whose frequencies are 2ωc. Recalling that cos(0) = 1, and sin(0) = 0,
the output of the top LPF is:
xr(t) = cos(ωt)
+
xa(t) = 2cos(ωt)
sin(ωt) –cos(ωt) –
HT HT
+
xb(t) = 0
cos(ωt) +
Figure S9–1
Solution: 9.2
(a) Given the 3-dimensional Xr(f) spectrum of xr(t) = Asin(2πfot) shown in
Figure S9–2(a), the three-dimensional spectrum of the Hilbert transform of
xr(t), Xi(f), is shown in Figure S9–2(b).
Solution:
Xr(f) Xi(f)
Imag Imag
A/2
Real Hilbert –fo Real
–fo transform
–A/2
0 fo 0 fo
–A/2
–A/2 Freq Freq
(a) (b)
Figure S9–2
(b) Representing the Hilbert transform of xr(t) as xi(t), from Figure S9–2(b) we
see that xi(t) = –Acos(2πfot). Thus the equation for the analytic signal, xa(t),
associated with xr(t) is
(e) The xa(t) analytic signal, on a complex plane, at time t = 0 is shown as the
dot in Figure S9–2(c).
Imag
jA
(c) –A A
xa(0) Real
t=0
–jA
Solution: 9.3
As shown by the table in Figure S9–3(a):
(a) y(t) = x(t).
(b) v(t) = –x(t).
(c) w(t) = HT–1[x(t)].
(d) The system that uses a single Hilbert transform operation to compute
inverse Hilbert transforms is shown in Figure S9–3(b).
Solution:
x(t) = x(t) =
HT3[x(t)]
cos(ωot) sin(ωot)
u(t), HT[(x(t)] sin(ωot) –cos(ωot) Inverse HT2[x(t)]
Hilbert
v(t), HT2[(x(t)] –cos(ωot) –sin(ωot) transform x(t)
HT
w(t), HT3[(x(t)] –sin(ωot) cos(ωot)
Inverse
y(t), HT4[(x(t)] cos(ωot) sin(ωot) –1 Hilbert
transform
(a) (b) of x(t).
Figure S9–3
Solution: 9.4
(a) We should expect their frequency magnitude responses to be zero (a
magnitude null) at zero Hz because the DC gain (gain at zero Hz) is
equal to the sum of the systems' impulse responses. For both Hilbert
transformers the sum of the impulse response samples is equal to zero.
(b) The z-domain transfer function for a 6-tap Hilbert transform filter is:
Solution: 9.5
There are two correct solutions to this problem. The equivalent structures of an
11-tap FIR Hilbert transformer are shown in Figure S9–5.
xr(n)
z–1 z–1 z–1 z–1 z–1 z–1 z–1 z–1 z–1 z–1
(a)
h(0) h(2) h(4) h(6) h(8) h(10)
+
xi(n)
+
xi(n)
Figure S9–5
Solution: 9.6
(a) There are several correct solutions. They are:
(I) H(z) = h(0) + h(1)z–1 + h(2)z–2 + h(3)z–3 + h(4)z–4 + h(5)z–5 + h(6)z–6
+ h(7)z–7 + h(8)z–8 + h(9)z–9 + h(10)z–10,
10
(II) H(z) = ∑ h(k)z–k.
k=0
2sin2(πn/2)
hfs=1(n) = = h'(n)
πn
which is what we set out to prove.
Solution: 9.8
The operations performed on X(m) to create an Xnew(m) sequence are shown
below:
N
Xnew(m) = –j·X(m), for positive frequencies (1≤m≤ 2 –1),
N
Xnew(m) = j·X(m), for negative frequencies ( 2 +1≤m≤N–1).
Solution: 9.9
The equation for the hhilb(k) coefficients when the half-band filter h(k)
coefficients' index k is defined using our standard notation of k = 0, 1, 2, ..., N–1
is found by replacing the original sin function's n index with
N–1
n=k– 2 .
This substitution shifts the sin function's indexing to correspond with our
traditional filter coefficient k indexing. The correct hhilb(k) expression is
Property of Pearson Education. Not permissible for redistribution.
π N–1
hhilb(k) = 2sin[ 2 (k – 2 )]h(k).
Solution: 9.10
Because h(0) = –h(6) and h(2) = –h(4), two equivalent alternatives of the
original 7-tap FIR Hilbert transformer, shown in Figure S9–10(a), implemented
to reduce the number multipliers by a factor of two are shown in Figure S9–
10(b) and Figure S9–10(c). The difference between the two structures are:
x'r(n)
xr(n) xr(n–2) xr(n–4)
z–1 z–1 z–1 z–1 z–1 z–1 xr(n–6)
+
xi(n)
Solutions:
Using the h(0) and h(2) coefficients
x(n–6) x(n–4)
z–1 z–1
x(n) x(n–2)
z–1 z–1 z–1 z–1
(b) – –
x'r(n)
h(0) h(2)
xi(n)
xi(n)
Figure S9–10
Solution:
This process is difficult to
implement in an accurate,
computationally-efficient, way Half-sample x'r(n)
delay
xr(n)
z–1 z–1 z–1 z–1 z–1
+
xi(n)
Figure S9–11
(b) The problem with an odd-order (even number of taps) Hilbert transformer
is that a non-integer (2.5 samples) delay is needed to obtain a time-aligned
x'r(n) sequence. Such non-integer, fractional, delay systems can be built, but
they require many arithmetic computations per output sample. From a
computational standpoint, even-order (odd number of taps) Hilbert
transformers are preferred because the x'r(n) output sequence is always
available at no computational cost.
Solution: 9.12
A complex-coefficient FIR filter's frequency magnitude response is periodic
because the filter's impulse response comprises discrete samples. All digital
filters have periodic frequency responses because the Fourier transforms of
discrete impulse response sequences are always periodic.
Solution: 9.13
(a) To determine an expression for the frequency response of a (K–1)th-order
FIR Hilbert transformer, first we write the transformer's time-domain
difference equation as
K–1
xi(n) = ∑ h(k) x (n–k).
r
k=0
Next we can describe the transfer function of the K-tap Hilbert transformer
defined as H(z) = Xi(z)/Xr(z). Rearranging the above expression yields
Κ–1
Xi(z)
H(z) = X (z) =
r ∑ h(k)z –k
.
k=0
Finally, substituting ejω for z in the H(z) transfer function gives us our
desired expression for the (K–1)th-order Hilbert transformer's H(ω)
frequency response, in radians per sample, as
Κ–1
H(ω) = H(z)|z = ejω = ∑ h(k)e –jkω
.
k=0
If the student substituted ej2πf for z in the H(z) transfer function, the desired
expression for the (K–1)th-order Hilbert transformer's H(f) frequency
response, in cyclic frequency Hz, would be
Κ–1
H(f) = HHT(z)|z = ej2πf = ∑ h(k)e –jk2πf
.
k=0
(d) For a 43rd-order Hilbert transformer, the above DFT equation will yield 44
frequency-domain
Property ofsamples.
PearsonTo increase theNot
Education. number of frequency-domain
permissible for redistribution.
samples, we would pad (append to) the h(k) impulse response sequence with
zero-values samples, and perform a larger-sized DFT of the lengthened
time-domain unit impulse response sequence.
Solution: 9.14
With K being an odd number, we can create a table whose entries are:
K Number of
multipliers
3 2
5* 2
7 4
9* 4
11 6
13* 6
15 8
17* 8
and so on.
* values where K – 1 is an integer multiple of 4.
K=5
xr(n)
z–1 z–1 z–1 z–1 xr(n)
z–1 z–1 z–1
0 h(1) 0 h(3) 0
h(1) h(3)
+
+ xi(n)
h(0) = h(2) = h(4) = 0 xi(n)
Figure S9–14
Solution: 9.15
Writing an expression for Xc(ω) in terms of Xr(ω) and H(ω), we have
Xc(ω) = XProperty
r(ω) + jXiof
(ω)Pearson
= Xr(ω) +Education.
jXr(ω)H(ω) Not permissible for redistribution.
= Xr(ω)[1 + jH(ω)].
⎧⎪ 1, for − π < ω ≤ 0
1 + jH (ω) = ⎨
⎪⎩ 0, for 0 < ω < π
⎧⎪ 0, for − π < ω ≤ 0
H (ω) = ⎨
⎪⎩ j, for 0 < ω < π.
Solution: 9.16
The problem's solution, found by determining the spectra at two intermediate
nodes, A(f) and B(f), is shown in Figure S9–16(d). The network eliminates all
negative-frequency spectral components in x(t).
A(f)
xr(t) yr(t) Imag
0.5 Real
xi(t) yi(t)
–5
HT 0 10
a(t) –1 b(t)
Hz
HT
1
(a) (b)
Solution:
B(f)
Imag Y(f) = X(f) + B(f)
Imag
Real Real
–5 1
0.5
0 2
10 0
10
Hz
Hz
(c) (d)
Figure S9–16
Solution: 9.17
Over the time interval 100 ≤ t < 150 milliseconds, for example, signal xa(t) is
xa(t)100≤t<150 = 3cos(ωt) + jHT[3cos(ωt)] = 3cos(ωt) + j3sin(ωt)
where the "HT[]" notation means Hilbert transform. The magnitude of xa(t)100≤t<150 is
2 2
|xa(t)|Property
100≤t<150 = of [3cos(ωt)] + [3sin(ωt)]Not
Pearson Education. 9[cos2(ωt) for
= permissible + sin 2
(ωt)] .
redistribution.
Using the trig identity: cos2(α) + sin2(α) = 1, we write
|xa(t)|100≤t<150 = 9 = 3.
Using that same analysis over the full signal time duration, the problem solution
is the |xa(t)| waveform shown in Figure S9–17. Signal |xa(t)| is called the
"envelope" of the original xr(t) signal.
Solution:
4
|xa(t)|
3
2
1
0.5
0
0 50 100 150 200 250
Time (milliseconds)
Figure S9–17
(b) Figure S10–1 shows the frequency magnitude response of an ideal (zero
transition region width) decimation-by-four lowpass filter.
Solution:
Filter mag. response
1
... ...
Figure S10–1
Solution: 10.2
There is no frequency difference between the x(n) and y(m) sequences. As
shown in Figure S10–2, the set of 36 samples in the original x(n) and y(m)
plots extend over two different intervals of time. Both sinusoids repeat every
18ts1 = 9ts2 = 18 milliseconds, so the x(n) and y(m) sinusoidal sequences
have identical frequencies.
Figure S10–2
Solution: 10.3
(a) Given that w(n)'s frequency is 128 Hz, the FFT size is N = 2048, and
fs = 1024 Hz, the positive-frequency value of index m for the FFT sample
having the largest magnitude in |W(m)| is
(128)N (128)2048
mmax = fs = 1024 = 256.
(b) Figure S10–3(a) shows the first 20 samples of the original w(n). Figure
S10–3(b) shows the first 20 samples of x(n). Now x(n)'s frequency is still
128 Hz, the FFT size is N = 1024, and fs = 512 Hz, the positive-frequency
value of index m for the FFT sample having the largest magnitude in |X(m)|
is
(128)N (128)1024
mmax,dec=2 = fs = 512 = 256.
(c) This is almost a trick question. Figure S10–3(c) shows the first 20 samples
of y(n). They are all zero-valued samples! Because y(n) are all zeros, the
Y(m) sequence is also all zeros. Therefore:
0
(a)
w(n)
–1
0 5 n 10 15
1
x(n)
(b) 0
–1
0 5 n 10 15
1
y(n)
(c) 0
–1
0 5 n 10 15
Figure S10–3
Solution: 10.4
(a) The sample rate of y(n) is fs2 = fs1/M = 1000/4 = 250 Hz.
(b) Because 100 Hz is within the unity-gain filter's passband, the peak
amplitude of the 100 Hz sinusoid in the w(n) sequence is P.
(c) Because decimation that does not violate the Nyquist criterion causes no
amplitude loss, the peak amplitude of the 100 Hz sinusoid in the y(m)
sequence is P.
(d) This is a little tricky for beginners in DSP. From Chapter 3, hopefully the
student will remember that DFT magnitudes are proportional to DFT size.
Because the x(n) 100 Hz sinusoid's 4N-point DFT magnitude was K,
the y(m) 100 Hz sinusoid's N-point DFT magnitude will be K/4.
(Magnitude loss by a factor of M = 4.)
The lessons to be learned from this problem, so far, are that decimation of a
finite-length sequence does not cause a change in time-domain amplitudes,
but decimation of a finite-length sequence does cause a loss in frequency-
domain magnitudes.
(e) The equation defining the downsampled y(m) sequence, in terms of w(n), is
y(m) = w(4n).
Solution: 10.5
Property of Pearson Education. Not permissible for redistribution.
The spectral magnitudes of the complex signals out of the filter pairs, and the problem's
solution of the |Xc(f)| spectrum, are provided in Figure S10–5.
wQ(n) xc(m)
(a) uQ(n)
hsin(k) hHP(k) 3 xQ(m)
|Xr(f)|
(b)
–12 –8 –4 0 4 8 12 kHz
(–fs/2) (fs/2)
|U(f)|
(c)
–12 –8 –4 0 4 8 12 kHz
(–fs/2) (fs/2)
|W(f)|
(d)
–12 –8 –4 0 4 8 12 kHz
(–fs/2) (fs/2)
Solution:
|Xc(f)|
(e)
kHz
–16 –12 –8 –4 0 4 8 12 16
(–2fs) (–fs) (fs) (2fs)
Figure S10–5
Solution: 10.6
(a) Using the text's Eq. (10–3), and Atten = 50, the number of taps in the LPF0
lowpass filter is
Atten 50
N LPF0 ≈ = ≈ 455 taps.
22( fstop − f pass ) 22(2.3/120 − 1.7 /120)
(b) With M = 30, B' = 1700 Hz, and the fstop frequency is 2300 Hz we compute
F = (fstop–B')/fstop = F = 600/2300 = 0.261.
From the text's Eq. (10–2), we estimate the optimum M1 decimation factor
as
1 − MF /(2 − F ) 1 − 30 ⋅ 0.261/(2 − 0.261)
M1,opt ≈ 2M = 60 = 11.05.
Property2 −ofFPearson
( M + 1) Education.
2 − 0.261(30 + 1)
Not permissible for redistribution.
The integer submultiple of 30 closest to 11.05 is 10, so we set
M1 = 10, and M2 = M/M1 = 30/10 = 3.
(d) The number of filter taps using the single-filter decimation system is 455
taps. The number of filter taps using the two-stage decimation system is
32+45 = 77 taps.
The reduction in the number of filter taps using the two-
stage decimation system is 455–77 = 378 taps.
LPF1 response
(a)
0 1.7 6 10.3 12 Freq
(fs,old /10) (kHz)
LPF2 response
(b) B'
0 1 1.7 2.3 3 4 Freq
(fs,new ) (kHz)
Figure S10–6
Solution: 10.7
We need 23 xold(n) input samples to fill the LPF1 filter with data. After that we
need M1 xold(n) input samples to LPF1 for each delay element in the LPF2 filter.
So the number of xold(n) input samples needed to fill both filters with input data
is:
Number of input samples = 23 + M1·(10–1) = 23 + 63 = 86 samples.
(b) Figure S10–8 shows the frequency magnitude response of an ideal (zero
transition region width) interpolation-by-three lowpass filter.
Solution:
Filter mag. response
3
... ...
–2π –π π 2π Rad./sample
(–fs) (–fs/2) –π/3 π/3 (fs/2) (fs) (Hz)
(–fs/6) (fs/6)
Figure S10–8
Solution: 10.9
(a) The equation defining the upsampled y(n) sequence, in terms of x(n), is
⎧⎪ x(n / 3), for n = 0,3,6,9,...
y ( p) = ⎨
⎪⎩ 0, otherwise.
(b) The 24-point Y(m) spectrum of the upsampled-by-3 y(n) sequence is shown
in Figure S10–9.
Solution:
|Y(m)|
4
0
0 7 9 15 17 23 m
Figure S10–9
The above Y(m) spectrum is obtained by merely repeating the original X(m)
spectrum three times over 24 frequency-domain samples.
Solution: 10.10
Property of Pearson Education. Not permissible for redistribution.
Upsampling (zero sample insertion) by a factor of two results in a spectrum that
is compressed by a factor of two. That is, given the original |Xold(f)| spectrum in
Figure S10–10(a), the spectral magnitude of xnew(m) will be the |Xnew(f)| shown
in Figure S10–10(b). The desired expression for the frequency points of spectral
symmetry of |Xnew(f)| is:
kfs
Points of spectral symmetry = ± 4
where k is an integer.
... ...
(a) |Xold(f)|
... ...
... ...
(b)
|Xnew(f)|
... ...
Figure S10–10
Solution: 10.11
The spectral magnitudes of the complex sequences at nodes A, B, C, and the
real part of the y(m) output sequence are shown in Figure S10–11.
|X(f)|
(a)
–4 –2 –fs/2 0 fs/2 2 4 Freq
(MHz)
(b)
–4 –2 –fs/2 0 fs/2 2 4 Freq
(MHz)
B
(c)
–4 –fs/2 –1 0 1 fs/2 4 Freq
(MHz)
C
(d)
–4 –fs/2 –1 0 1 fs/2 4 Freq
(MHz)
Output of 2nd
multiplier
(e)
–4 –fs/2 –1 0 1 fs/2 4 Freq
(MHz)
Solution:
|Y(f)|
(f)
–4 –fs/2 –1 0 1 fs/2 4 Freq
(MHz)
Figure S10–11
Solution: 10.12
(a) The solution to this problem is interpolation combined with appropriate
bandpass filtering as shown in Figure S10–12(a).
Solution:
4
x(t) u(n) v(m) Bandpass w(m) y(m)
A/D 4
(a) filter, h(k)
Figure S10–12
(b) Justification of this design is provided in the spectral plots in Figure S10–
12(b) through (e). Interpolation by four exhibits a gain loss of four as shown
by the K/4 value in Figure S10–12(c), so if the bandpass filter has a gain of
one, we need to multiply the w(m) sequence by four.(We can eliminate the
multiplier by designing a filter whose passband gain is four.)
Property of Pearson Education. Not permissible for redistribution.
Solutions:
|U(f)|
... K ...
(b)
–3 –2 –1 0 1 2 3 kHz
–4 4
(–fs1/2) (fs1/2)
|V(f)|
K/4
... ...
(c)
–12 –8 –4 0 4 8 12 kHz
–16 16
(–fs2/2) (fs2/2)
|W(f)|
... K/4 ...
(d)
–12 –8 –4 0 4 8 12 kHz
–16 16
(–fs2/2) (fs2/2)
|Y(f)|
... K ...
(e)
–12 –8 –4 0 4 8 12 kHz
–16 16
(–fs2/2) (fs2/2)
Solution: 10.13
(a) Based on the interpretation in Figure S10–13(a), the difference equation for
linear interpolation is
x(n) – x(n–1)
y(n) = x(n–1) + 2 = 0.5x(n) + 0.5x(n–1).
y(n) x(1)
x(0)
0
1 n
Property of Pearson Education. Not permissible for redistribution.
Figure S10–13
(b) The linear interpolator's frequency response is its H(z) transfer function
with variable z replaced with ejω, or:
H(ω) = 0.5 + 0.5e–jω = 0.5 + 0.5cos(ω) –j0.5sin(ω).
The frequency magnitude response of a linear interpolation filter is
|H(ω)| = |0.5 + 0.5cos(ω) –j0.5sin(ω)|
= [0.5 + 0.5cos(ω)]2 + [0.5sin(ω)]2
= 0.25 + 0.5cos(ω) + 0.25cos2(ω) + 0.25sin2(ω)
= 0.25 + 0.5cos(ω) + 0.25 = 0.5 + 0.5cos(ω) .
A sketch of this |H(ω)| over the frequency range of ω = ±π radians/sample
(±fs/2 Hz) is given as the solid curve in Figure S10–13(b).
Solution:
1
|H(ω)|
0.8
0.6
(b)
0.4
0.2
0
–π/2 –π/4 0 π/4 π/2
(–fs/2) (–fs/4) Frequency (fs/4) (fs/2)
Advantages:
• Simple to compute. It's merely a two-point moving averager.
• The multiplies by 0.5 can be performed by binary data right shifts
requiring no actual multiplication operations.
Disadvantage:
• Linear interpolation's frequency magnitude response is highly inferior
to a high-performance interpolation by two lowpass filter, shown as
the dashed lines in Figure S10–13(b).
• Least accurate of all methods of interpolation.
Solution: 10.14
Sample rate conversion by a rational factor establishes the following
relationships:
L ⋅ f s ,in L ⋅ f s ,CD
f s ,out = = f s ,DAT =
.
PropertyMof Pearson Education.
M Not permissible for redistribution.
or
L ⋅ f s ,CD 160 ⋅ 44100
M = = = 147.
f s ,DAT 48000
Solution: 10.15
The solutions to this problem are:
2
Length = 20
0
xo(n)
Decimation of an –2
even-length
0 5 n 10 15
sequence. Time ts = 1 msec. Time = 19ts = 19 msec.
durations are not
equal. 2 Decimation of an
Length = 10 odd-length
0 sequence. Time
xD(mD) durations are
–2 equal.
0 2 4 mD 6 8
tsD = 2 msec. Time = 9tsD = 18 msec.
2
Length = 39
0
xI(mI)
–2
0 5 10 15 20 mI 25 30 35
tsI = 0.5 msec. Time = 38tsI = 19 msec.
Figure S10–15
Solution: 10.16
The table should be filled in as follows:
Solution: 10.17
The desired spectral plots are as follows (dashed curves are spectral replications
at multiples of the sample rate):
Solutions:
|Q(f)| Spectral replications
K/4
(a)
-1600 -1200 -800 -400 0 400 800 1200 1600 Freq
(fs) (Hz)
-100 100
(b)
-1600 -1200 -800 -400 0 400 800 1200 1600 Freq
(fs) (Hz)
-100 100
Figure S10–17
Solution: 10.18
(a) The frequency responses of the L(f) lowpass and H(f) highpass filters are
shown in Figure S10–18(a). The system's first scrambling network is shown
in Figure S10–18(b), with signal nodes A through F appropriately marked.
The spectral magnitude of the x(n) input to the first scrambling network is
shown in Figure S10–18(c). The spectral magnitudes at the various signal
nodes are shown in Figures S10–18(d) through Figures S10–18(f).
The solution to this Part (a) of the problem, the spectrum of the first
scrambling network's output, is shown in Figure S10–18(g). (Compare that
spectrum with the original X(f) spectrum in Figure S10–18(c).)
|L(f)| |H(f)|
(a)
–8 –4 0 4 8 kHz –8 –4 0 4 8 kHz
(–fs) (fs) (–fs) (fs)
|X(f)|
(c)
–8 –4 0 kHz
4 8
(–fs) (fs)
Node A Node D
(d)
–8 –4 0 4 8 kHz –8 –4 0 4 8 kHz
(–fs) (fs)
Node B Node E
(e)
–8 –4 0 4 8 kHz –8 –4 0 4 8 kHz
(–fs) (fs)
Node C Node F
(f)
–8 –4 0 4 8 kHz –8 –4 0 4 8 kHz
(–fs) (fs)
(g)
–8 –4 0 4 8 kHz
(–fs) (fs)
Figure S10–18
(b) Because the scrambling networks swap the low- and high-frequency
spectral components of their inputs,
Solution: 10.19
(a) Polyphase filters are useful for interpolation because no unnecessary
multiplications are performed on zero-valued time samples.
Solution:
x(n) at
fs rate
H0(z) 4
z–1
H1(z) 4
z–1
H2(z) 4
y(m) at
z–1 4fs rate
H3(z) 4
Figure S10–19
Solution: 10.20
This problem can be a bit tricky for a DSP novice to contemplate. The solution
is to maintain the time relationships between the four time-domain sequences.
That is, make sure the sequences are delayed from each other by one sample-
time period, with the oldest (earlier in time) sequence being applied to the "A"
input port. Using standard "z–1" delay elements, the solution is shown in Figure
S10–20(a). (If desired, the system could be re-drawn as that in Figure S10–
20(a).)
(a) (b)
Figure S10–20
Solution: 10.21
(a) In the first system 12 multiplications must be performed for each x(n) input sample. With
30 input samples arriving per second, the total multiplication rate is:
= 360 multiplications/second.
(b) In the second system 12 multiplications must be performed for every three x(n) input
samples. With 30 input samples arriving per second, the total multiplication rate is:
Solution: 10.22
The decimation by four filter using the switches is shown in Figure S10–22(a),
and a decimation by four polyphase filter having 12 multipliers is shown in
Figure S10–22(b). Let's define ts as the time between xold(n) samples.
xnew(m)
z–1 z–1
h(3) h(7) h(11)
z–1 z–1
+
(b)
z–1 z–1
xnew(m)
h(1) h(5) h(9) +
z–1 z–1
h(2) h(6) h(10)
Figure S10–22
Solution: 10.23
(a) Polyphase filters are useful for decimation because they only
compute output samples that will be retained. No unnecessary
computations are performed.
Figure S10–23
Solution: 10.24
(a) Resampling by the rational factor 5/4 mandates that we perform interpolation by L = 5
followed by decimation by M = 4. With output y(m)'s index is m = 7, the commutating
switch's port position value (index) k is computed as
where the subscripted 5 means modulo-5. The index n of the most recent
x(n) input sample applied to the subfilters, when m = 7, is computed as
(b) Interpolation by L has an inherent gain loss by a factor of L. For the resampler to have a DC
(zero Hz) gain of unity the original prototype lowpass FIR filter must have a
DC gain of L = 5.
Solution: 10.25
Given the spectrum of the filter output a(n) sequence shown in Figure S10–
25(a), the spectra of the b(n), c(m), and y(p) sequences are those shown in
Figures S10–25(b), (c), and (d).
|A(f)|
... ...
(a)
–1600 –800 –400 0 400 800 1600 Freq
(fs)
... ...
(b)
–800 –400 0 400 800 Freq
(fs)
|C(f)|
... ...
(c)
–800 –400 –200 0 200 400 800 Freq
(fs)
|Y(f)|
... ...
(d)
–400 –200 0 200 400 Freq
(fs)
Figure S10–25
Solution: 10.26
(a) The answer is no, the shapes of a comb filter's passband curves are not a
function of f2 which would make them parabolic. We prove this statement
by finding the equation for the frequency response of a comb filter. That
frequency response is the Hcomb(z) transfer function evaluated on the unit
circle. We start by substituting ejω for z in Hcomb(z), because z = ejω defines
the unit circle, giving
Hcomb(ejω) = Hcomb(z)|z=ejω = (1 - e-j8ω).
Factoring out the half-angled exponential e-jω8/2, we have
Hcomb(ejω) = e-jω8/2 (ejω8/2 - e-jω8/2) = e-j4ω (ej4ω - e-j4ω).
Using Euler's identity 2jsin(α) = ejα - e-jα, we arrive at
Hcomb(ejω) = e-j4ω [2jsin(4ω)] = je-j4ω [2sin(4ω)]
or
Hcomb(f) = je-j4ω [2sin(8πf)].
Ignoring the phase shift term (complex exponential) above, the frequency-
domain magnitude response of a comb filter is
|Hcomb(ejω)| = |2sin(4ω)|
or
|Hcomb(f)| = |2sin(8πf)|.
Solution: 10.27
The impulse response of the decimation CIC filter, shown in Figure S10–27(a),
can be obtained from the y(n) column in the table in Figure S10–27(b). At time
index n = 5 sequence w(n–5) becomes all ones and the comb's output becomes
all zeros. The y(n) impulse response is plotted in Figure S10–27(c).
Solution:
n x(n) w(n) w(n–5) y(n)
Decimation CIC filter –1 0 0 0 0
Integrator Comb 0 1 1 0 1
1 0 1 0 1
x(n) y(n) 2 0 1 0 1
w(n)
3 0 1 0 1
–
4 0 1 0 1
z–1 z–5 5 0 1 1 0
6 0 1 1 0
w(n-1) w(n-5)
7 0 1 1 0
(a) (b)
1
(c)
0.5
0 ...
0 5 10
n
Figure S10–27
Solution: 10.28
(a) The solution for this part of the problem are the u(n) and y(n) sequences for
the D = 4 interpolation and decimation CIC filters shown in Figure S10–28.
1 ...
0.5
0
0 5 n 10
u(n) in decimation
.
CIC filter .. y(n) in decimation CIC filter
10
6
8 ...
4
6
2
4
0
2 0 5 n 10
0
0 5 n 10
Figure S10–28
The u(n) sequence in the decimation filter is unbounded, but this does not
mean that decimation CIC filters are not used in practice. (We never apply a
step sequence to CIC decimators in practice. In practice, we apply signals
that have both positive and negative sample values to CIC decimation
filters.)
(b) The number of binary bits needed to accommodate the u(n) and y(n) samples,
for each CIC filter, up to time index n = 500 are shown in Table S10–2. At
n = 500, u(500) = 501 requiring a 9-bit memory location (or hardware register)
for the decimation CIC filter.
Solution: 10.29
(a) This problem is solved by evaluating the 5th-order CIC filter's frequency
magnitude response equation at the frequency fo indicated by Atten. That
equation is
Property of Pearson Education. Not permissible for redistribution.
M
⎪ sin(πfD) ⎪
|Hcic,Mth-order(f)| = ⎪ ⎪
⎪ sin(πf) ⎪
with M = 5 and D = 6. Because the first CIC filter magnitude null is at
fs,in/D = 0.1667 as shown in Figure S10–29, our fo frequency of interest
(normalized to fs,in) is
fo = 1/D – B/2 = 0.1667 – 0.04/2 = 0.1467
where bandwidth B is also normalized to fs,in.
0
|HCIC(f)|
B Atten
dB
Figure S10–29
or 86 dB below the CIC filter's passband peak, which is what we set out to
find.
(b) The Gain loss is found in a similar manner by evaluating the 5th-order CIC
filter's frequency magnitude response equation at the normalized frequency
f = B/2 = 0.02. Doing so we have
5
⎪ sin(0.02π/6) ⎪
|Hcic,Mth-order(0.02)| = ⎪ ⎪ = (5.863)5 = 6.926x103.
⎪ sin(0.02π) ⎪
The magnitude response at f = 0.02 in decibels is
Gain Loss = 20 log10( 6.926x103/65) = –1 dB,
or 1 dB below the passband peak, which is what we set out to find.
Solution: 10.30
The number of unit-delay elements in the upper path of the two-path filter
must beProperty
equal to of
thePearson Education.
group delay of the CICNot permissible
filter. formay
The student redistribution.
recall
that a D = 9 CIC filter is equivalent to a D = 9-tap tapped-delay line moving
average filter (having symmetrical coefficients) as shown in Figure S10–30.
Because that moving average filter (and the CIC filter) has eight delay
elements, its group delay is:
the upper path of the parallel filter must contain four unit-delay elements.
x(n)
z–1 z–1 z–1 z–1 z–1 z–1 z–1 z–1
y(n)
Figure S10–30
Given that Hø(f), as shown in Appendix F, the group delay of this filter is:
–1 d(Hø(f)) –1 d(–8πf)
G(f) = · d(f) = · d(f) = 4 samples.
2π 2π
So, for a D = 9 CIC filter in the bottom path, the upper path of the parallel
filter must have four unit-delay elements.
1 N 1 4 1
xave = N ∑ x(n) = 4 ∑ x(n) = 4 [1 + 2 + 3 + 4] = 2.5.
n=1 n=1
n=1
1
xvar = 4 [(–1.5)2 + (–0.5)2 + (0.5)2 + (1.5)2] = 1.25.
Solution: 11.2
(a) By its inherent nature, xave is the single value such that the sum of
diff(n) = x(n) – xave is zero. That is,
6 6
∑ diff(n) = ∑ [x(n) – xave] = 0.
n=1 n=1
Solution: 11.3
The answer is yes, it is valid to state that zave = xave + yave.
We verify this as follows: Given that
1 N 1 N
zave = N ∑ z(n) = N ∑ [x(n) + y(n)]
Property
n=1 of Pearson
n=1 Education. Not permissible for redistribution.
[x(1) + y(1)] + [x(2) + y(2)] + [x(3) + y(3)] + ... + [x(N) + y(N)]
= N ,
we can write
x(1) + x(2) + x(3) + ... + x(N) + y(1) + y(2) + y(3) + ... + y(N)
zave = N
x(1) + x(2) + x(3) + ... + x(N) y(1) + y(2) + y(3) + ... + y(N)
= N + N
1 N 1 N
= N ∑ x(n) + N ∑ y(n) = xave + yave.
n=1 n=1
which is what we set out to prove.
Note: The statement that zave = xave + yave tells us that the average of a signal-
plus-additive-noise is equal to the average of the noise-free signal plus the
average of the noise signal. If the noise signal's average is zero (noise
samples that are equally likely to be plus or minus), we can measure the
average of a signal by computing the average of the signal-plus-additive-
noise. This is what signal averaging is all about!
Solution: 11.4
The average phase angle (in degrees) of π/4 radians, –3π/4 radians, and –π/4
radians is found by first finding the average of three complex numbers
having the appropriate arguments, as:
ejπ/4 + e–j3π/4 + e–jπ/4
Averagecomplex = ( 3 )
The argument of e–jπ/4 is –π/4 radians, so the average phase angle in radians
is:
Average phaseradians = tan–1(e–jπ/4) = –π/4 radians.
k = 20 FFTs.
Solution: 11.6
Using the geometric series-to-closed form expression description in Appendix
B, we can write Hma(z) as
1 1 – z–N
= N 1 – z–1 = Hrma(z)
Solution: 11.7
To determine the frequency magnitude responses of a recursive running sum
filter when N = 4, 8, and 16, we take the DFT of the filter's impulse response for
the various N values. Those impulse responses are rectangular sequences of all
ones, with lengths of 4, 8, and 16. The DFT of those impulse responses are
sin(x)/x-like functions whose mainlobe amplitudes are equal to N, with
magnitude nulls at multiples of fs/N, as shown in Figure S11–7.
Solutions:
Recursive running sum
frequency magnitude response
16
N = 16
12
N=8
8
N=4
4
0
0 fs/8 fs/4 fs/2
Solution: 11.9
The solution to this problem is the bold curve shown in Figure S11–9(b). Over
the positive frequency range, the location of the N = 3 averager's z-plane zero is
found using the text's Eq. (11–26) for k = 1. That zero location is
2π 2π
θzeros = k N = 3 radians
as shown in Figure S11–9(a). From the text's Eq. (11–26'), the cyclic frequency
associated with an angle of 2π/3 radians is fs/3 Hz as shown by magnitude null
at fs/3 Hz in the bold curve in Figure S11–9(b). The magnitude at fs/2 (half the
sample rate) is found using the text's Eq. (11–25) with f = 0.5 as
1 sin(πfN) 1 sin(3π/2) 1 –1 1
|Hma(f)| = N · | | = 3 ·| | = 3 ·| 1 |= 3 .
sin(πf) sin(π/2)
Solution:
θ = 2π/3 radians |Hma(f)|
= fs/3 Hz 1
1 0.8
Imaginary part
Figure S11–9
Solution: 11.10
(a) Cascading an N = 4 moving averager and an N = 2 moving averager results
in a cascaded network whose frequency magnitude response is the product
of the individual moving averagers' frequency magnitude responses as
shown in by the bold curve in Figure S11–10(a).
Figure S11–10
(b) The phase of the cascaded (two-stage) filter is linear because the product of
two linear-phase frequency responses is itself linear phase.
Another way to show that the phase of the cascaded (two-stage) filter is
linear is to remember that the combined hcasc(k) impulse response of two
cascaded filters is the convolution of their individual impulse responses. As
shown in Figure S11–10(b), convolving the two moving averagers' impulse
responses yields a cascaded hcasc(k) impulse response that is symmetrical,
meaning that the cascaded (two-stage) filter exhibits linear phase.
h1(k)
0.5
0.25
0 hcasc(k)
0 1 2 3 k 0.5
(b) Convolve 0.25
h2(k) 0
0.5 0 1 2 3 4 k
0.25
0
0 1 k
Solution: 11.12
When α = 0 the filter's feed forward coefficient is zero and the filter output
samples will be an all-zero sequence. (What filter could be more stable than
this?) Thus we would never actually use a value of α = 0.
Solution: 11.13
The alternate (feedback only) exponential averager's transfer function is:
Y(z) 1
H(z) = W(z) = .
1– (1–α)z–1
Replacing z with ejω yields the feedback loop's frequency response of:
1
H(ω) = .
1– (1–α)e–jω
At zero Hz, ω = 0, the feedback loop's frequency response (its DC gain) is:
1 1 1
H(ω)|ω=0 = –j0 = = .
1– (1–α)e 1– (1–α) α
Solution: 11.14
Given the averager's frequency response of
α
Hexp(ω) =
1 – (1–α) cos(ω) + j(1–α).sin(ω)
.
the desired frequency magnitude response is found by finding the magnitude
expression for the complex denominator of Hexp(ω) as:
⎪ α ⎪
|Hexp(ω)| = ⎪ . . ⎪
⎪ 1 – (1–α) cos(ω) + j(1–α) sin(ω) ⎪
α
=
[1 – (1–α).cos(ω)]2 + [(1–α).sin(ω)]2
Squaring the terms in brackets, our solution is:
α
|Hexp(ω)| = .
1 – 2.(1–α).cos(ω) + (1–α)2
Solution: 11.15
Property of Pearson Education. Not permissible for redistribution.
When a unit impulse sequence (1,0,0,0,0,...,) input is applied to the standard
integrator, the unity-valued sample's contribution to the y(n) output sequence
lasts forever as shown in Figure S11–15(a). When a unit impulse is applied to
an exponential averager, the unity-valued sample makes an immediate
contribution to the averager's output sequence, however that sample's
contribution to the output diminishes with time (it leaks away) as shown in
Figure S11–15(b).
0.1 0.1
0
0 1 2 3 4 5 6 7 8 9 0
n 0 1 2 3 4 5n6 7 8 9
(a) (b)
Figure S11–15
Solution: 11.16
(a) The exponential averager's response, to a unity-valued input sample applied
at time n = 0, is
h(0) = α,
h(1) = α(1–α),
h(2) = α(1–α)(1–α) = α(1–α)2,
h(3) = α(1–α)(1–α)(1–α) = α(1–α)3,
h(4) = α(1–α)(1–α)(1–α)(1–α) = α(1–α)4,
etc.
The student should see a pattern occurring here, and write:
h(0) = α,
h(n) = α(1–α)n, for n > 0.
The above h(0) and h(n) results are an acceptable solution to this problem.
However, a perceptive student may realize that h(0) = α = α(1–α)0, and
write the h(n) impulse response in the more concise form of:
h(n) = α(1–α)n, for n ≥ 0.
which is what we seek.
Note: The hyperperceptive student may notice that there is a problem with
the h(n) formula when α = 1, because then h(0) = 00. Now, 00 is sometimes
Property
regarded of Pearson
as indeterminate. AsEducation.
Antoine TruxNot permissible
points for redistribution.
out for example, the
Standard C library pow(x,y) (which computes xy) generates an error if
x = y = 0. However, most mathematicians now consider 00 to be equal to 1,
and this agrees well with the first output sample (h(0) = 1) of our impulse
response when α = 1.
(b) Here we hope the student remembers that a filter's gain at zero Hz (DC
gain) is the sum of the filter's impulse response samples. Using the h(n)
expression from Part (a), we can write the exponential averager's gain at 0
Hz (DC gain) as:
∞ ∞
Hexp(ω)|ω=0 = Hexp(0) = ∑ h(n) = ∑ α(1–α)n.
n=0 n=0
We could have, just as well, evaluated the text's Eq. (11–34) Hexp(ω)
frequency response, repeated here as
α
Hexp(ω) = ,
1 – (1–α)e–jω
at ω = 0 . Doing so we have
α α
DC gain = Hexp(0) = –j0 = = 1.
1 – (1–α)e α
which is what we set out to prove.
Note: We ignore the indeterminate case when α = 0 because the exponential
averager has all zero-valued output samples in that situation.
Solution: 12.2
(a) $A231 = 4152110.
(b) 0x71F = 182310.
Solution: 12.3
Hex numbers $07 and $E2 converted to binary format are:
$07 = 0000 0111 and $E2 = 1110 0010.
To perform the desired subtraction we add the two's complement of 1110 0010
($E2) to 0000 0111 ($07). Doing so yields:
710 ← ($07)
–(–3010) ← ($E2)
--------------
3710
Solution: 12.4
(a) Converting $45 to binary format results in: 0100 01012.
Sign-extending to sixteen bits yields: 0000 0000 0100 01012.
Converting back to hex format gives the solution: $0045.
Property of Pearson Education. Not permissible for redistribution.
(b) Converting $B3 to binary format results in: 1011 00112. (negative number!)
Sign-extending to sixteen bits yields: 1111 1111 1011 00112.
Converting back to hex format gives the solution: $FFB3.
Solution: 12.5
Using the symbol "◊" to represent a binary point:
Solution: 12.6
With 16-bit data words (b = 16), in the two's complement integer number format
the most positive representable decimal number is:
most positive = 2b–1–1 = 215 –1 = 32,76710
and the most negative representable decimal number is:
most negative = –2b–1 = –215 = –32,76810.
Solution: 12.7
The solution is given in the following tables:
Solution: 12.8
Multiplying decimal 0.165 by decimal 32768 gives 540.672. Rounding that
product to the nearest integer yields 541. Converting decimal 541 to a binary
integer proceeds as follows:
Placing the binary point to the right of the MSB of the 16-bit binary number
gives us our binary and hexadecimal solutions of
Tax rate = 0◊000 0010 0001 1101 = 021D16 ($021D).
Solution: 12.9
The number system allowing decimal 1/3 to be exactly represented with a
finite number of digits must have a base that is a multiple of three.
For example; using fixed-point numbering systems having bases of three, six,
and nine, decimal 1/3 can be represented by
1/310 = 0.13 = 0.26 = 0.39.
Solution: 12.10
This is a trick question! In a base 6 numbering system, the
digit "7" in 42736 has no meaning (is not defined).
Property of Pearson Education. Not permissible for redistribution.
Solution: 12.11
(a) To a sixteen decimal digit precision, the maximum positive decimal value is
2(# of integer bits – 1) – 2–(# of fraction bits) = 2(1 – 1) – 2–31 = 1 – 2–31
Solution: 12.12
The frequency resolution of the AD9958 chip, in Hz, is the signal
generator's maximum frequency (250 MHz) minus the generator's minimum
frequency (0 Hz) divided by two raised to the power of the number of
frequency-control bits. That is:
Solution: 12.13
The combined data output rate of the digital portion, measured in bytes (8-bit
binary words) per second, of a stereo CD player is:
Solution: 12.14
As stated in the text's Eq. (12–22) with regard to digital networks, such as
that shown in Figure S12–14, the hardware register containing the y(n)
sequence must have a word width that accommodates a value equal to the
network's DC (zero Hz) gain G times the input signal, i.e., G·x(n). That
number of y(n) bits is
y(n) bits = number of bits in x(n) + ⎡log2(G)⎤
z–1
0.85 z–1
–0.12
Figure S12–14
Next we set z = 1 giving us our desired DC (zero Hz) gain solution of:
1
G = H(z)|z=1 = 1 –0.85 + 0.12 = 3.704.
= 8 + ⎡1.89⎤ = 8 + 2 = 10 bits.
Solution: 12.15
The decimal values of all possible filter coefficient values using a four-bit
unsigned binary words in a 2.2 (two dot two) "integer plus fraction" format are
given in the following tables:
Solution: 12.17
As shown in Figure S12–17, the input to the summation operation is
x(n)/22 + x/24 + x/26 + x/28 + x/210 + x/212 + x/214
≈ x(n)·(0.3333) ≈ x(n)/3.
x(n) 2
2
2
x(n)/22 2
x(n)/24 2
x(n)/26 2
x(n)/28
x(n)/210 2
x(n)/212
x(n)/214
x(n)/3
Figure S12–17
Solution: 12.18
The solution to this problem is A = 3, and B = 1 as shown in Figure S12–18.
9 . x(n)
x(n)
Left shift
by 3 bits y(n) = 54.x(n)
Left shift
by 3 bits –
Left shift
by 1 bit
Multiplication
by 9, (23 + 1) Multiplication
by 6, (23 - 21)
Property of PearsonFigure
Education.
S12–18Not permissible for redistribution.
Solution: 12.19
The data-word bit manipulations necessary to eliminate the two multipliers in
Figure S12–19(a) are negation and binary bit shifting. The desired multiplier-
free implementation is shown in Figure S12–19(b).
• The two's complement negation comprises inverting the bits of the x(n–2)
sample and adding one. The bits of the negated result are arithmetically shifted
left by one bit before being added to x(n) and x(n–4).
• The magnitude bits of the sum, w(n), are then shifted right by two bits.
An "Arithmetic left shift" in Figure S12–19(b) means that a zero bit is inserted
as the new least significant bit (LSB). After that shift, the overflow and MSB
bits must be examined to maintain proper data word polarity. An "Arithmetic
right shift" means that the appropriate bit, to maintain proper polarity, is
inserted as the new most significant bit.
Two's complement
negation
x(n) x(n–2) x(n–4)
z–2 z–2
Arithmetic left
–2 shift by 1 bit
y(n)
w(n) w(n)
(a) (b)
Figure S12–19
Solution: 12.20
With a sample rate of fs = 2x109 samples/second, the time period between
x(n) samples is:
1 1
ts = f = 2x109 = 0.5x10–9 seconds.
s
128·5
= 511 ≈ 1.25 volts.
Solution: 12.22
(a) The A/D converters' quantization-level voltage, q, is the converter's full
peak-peak voltage range divided by two raised to a power equal to the
number of converter bits. Thus the quantization-level voltage for a 12-bit
A/D converter is:
10 volts 10 volts
q= 2 12 = 4096 = 0.0024 volts = 2.4 millivolts.
(b) The A/D converters' maximum positive and maximum negative quantization
error (noise) voltages are plus and minus half the converter's quantization
voltage:
±q ±0.0024 volts
max plus and minus quantization error = 2 = 2
(c) With a 7-volt peak-peak sinusoidal voltage at the converter's analog input,
we use the text's Eq. (12–13) to compute the converter's SNRA/D output
signal to quantization noise value, in dB. That SNRA/D value is:
SNRA/D = 6.02(12) + 4.77 + 20log10(LF) = 77.01 + 20log10(LF).
The loading factor, LF, is
rms voltage of the input (7/2)(0.707)
LF = peak converter input voltage = 5 = 0.495.
Plugging the LF = 0.495 value into the SNRA/D expression, we compute our
desired SNR result as:
SNRA/D-max
90
80
71 70
60 12-bit
dB
50 10-bit
40
8-bit
30
6-bit
20
Figure S12–22
Solution: 12.23
Given that a 12-bit A/D converter's signal to quantization noise level is 67
dB, when driven to its full-scale input voltage range, we compute the
desired effective bit value beff using the text's Eq. (12–16) as:
Solution: 12.24
(a) The quantization error as a function of the continuous x(t) input for the
truncating and rounding A/D converters are shown as the bold diagonal lines
in Figure S12–24.
Figure S12–24
(b) The table entries of quantization error properties, in terms of the A/D
converters' quantization-level voltage q, are provided in the following table.
Solution: 12.25
To determine the A/D converter's number of bits, b, to accommodate the
accuracy and operating range of the thermocouple, we first determine the
dynamic range of the measurement system's voltage v(t). That is done as
follows:
Temp. range
Voltage v(t) dynamic range = Temp. resolution
Solution:
# of occurrences
(a)
0
–1 –0.5 0 0.5 1
Sample value
Figure S12–26
That bathtub-like curve in Figure S12–26(a) shows that there are more
occurrences of A/D output samples whose values are near ±1 than
occurrences of sample values close to zero amplitude.
We can see that this is true from the sinewave A/D converter output in
Figure S12–26(b), where we see that there are more samples in the
amplitude Range-1 and Range-3 (near ±1) than there are in the amplitude
Range-2 (near zero). Histogram testing of A/D converters is described in
more detail in Chapter 13.
The precise histogram of an ideal A/D converter's output samples,
when the converter's input is an analog sinewave whose peak value is one, is
shown in Figure S12–26(c).
1
Range-1
# of occurrences
0.5
0 Range-2
–0.5
Range-3
–1
–1 –0.6 –0.2 0 0.2 0.6 1
0 Time Sample value
(b) (c)
Solution: 12.27
To determine the value of A of a uniform pdf whose variance is equal to
two, we use Appendix D's Figure D–4, repeated below in Figure S12–27,
and Eq. (D–12) from Appendix D. That equation, which applies to Figure
Property
S12–27, is of Pearson Education. Not permissible for redistribution.
(b + a)2
2
uniform pdf variance = σ = 12 .
p(f)
1
b+a
–a 0 b f
Figure S12–27
For this problem, the limits of the Figure S12–27 uniform pdf are
a = A, and b = A. So we can write the above σ2 expression as:
(A + A)2 4A2 A2
σ2 = 2 = 12 = 12 = 3 .
Solving the above for A we have:
A2 2
3 = 2, or A = 6,
so
A= 6 ≈ 2.4495.
Solution: 12.28
Given a single data sample value in binary floating point format, we can
perform a multiply-by-4 by:
(i) Shifting the sample's fraction bits to the left by two bits, or
(ii) Incrementing the sample's exponent bits two times.
Solution: 12.29
To solve this problem we use the text's expression
valueIEEE = (–1)s . 1◊ f . 2e–127 ,
where symbol ◊ denotes the binary point and f and e are decimal numbers.
To examine the floating-point number's individual bits, we convert the given
hex number to binary format as:
$C2ED0000 = 1100 0010 1110 1101 0000 0000 0000 0000
= –1.8515625 . 64 = –118.510.