Oppenheim Assignment Solutions
Oppenheim Assignment Solutions
Solutions to
Recommended Problems
S1.1
(b) Similarly,
Im{z}= Im 2 2 4
(c) The magnitude of z is the product of the magnitudes of 2 and eJT/4. However,
I| = i, while Iei"| = 1 for all 0. Thus,
IzI = e =I
- ie 4
| - 1 =
(d) The argument of z is the sum of the arguments of i and e'"/4 . Since 4 = 0 and
4 ei" = 0 for all 0,
= 4 e = + +es"=
(e) The complex conjugate of z is the product of the complex conjugates of 1 and
ei"/4 . Since * = and (ej")* = e ~i' for all 0,
z* =* * 4 =i -i'/4
(f) z + z* is given by
1r/4 1 . 4 ej"/4 + e /4 Vr2\/
z + z* =e + e-/4 = = cos
2 2 2 4 2
Alternatively,
z+ z*
Refz} 2, or z + z* = 2Re{z}= 2 4 2
S1.2
(a) Express z as z = a + jQ, where Re{z} = a and Im{z} = Q. Recall that z* is the
complex conjugate of z, or z* = a - jQ. Then
z + z* (a+jQ)+ (o -jQ) 2a + 0
2 2 2
(b) Similarly,
z - z* (a +fj) - (a -jQ) _ + 2jQ _.
2 2 2
S1-1
Signals and Systems
S1-2
S1.3
(a) Euler's relation states that e'" = cos 0 + j sin 0. Therefore, e-o = cos(-0) +
j sin (-0). But, cos (-0) = cos 0 and sin (-0) = -sin 0. Thus, e-j" = cos 0
j sin 0. Substituting,
e" + e - _ (cos 0 + j sin0) + (cos 0 - j sin0) 2 cos 0
2 2 2
(b) Similarly,
ej"- e -_ (cos 0 + j sin 0) - (cos 0 - j sin0) 2j sin0 s
= =-sin 6
2j 2j 2j
S1.4
(a) (i) We first find the complex conjugate of z = re'". From Euler's relation,
rej" = r cos 0 + Jr sin 0 = z. Thus,
z* = r cos 0 - jr sin 0 = r cos 0 + jr (-sin 0)
But cos 0 = cos (-0) and -sin 0 = sin (-0). Thus,
z* = r cos (-0) + jr sin (-0) = re -j"
2 20
z 2 = (re'")2 = r 2 (ejo) = r ei
2
(ii)
(iii) jz = er/ 2reo = reO+<,/2)
(iv) zz* = (re")(re-jo) = r 2em- 0 )>= r2 1
z 0
rej" e o+)e20
() z* re ~j"
(vi) = j__= e - e
z re 0 r
(b) From part (a), we directly plot the result in Figure S1.4-1, noting that for z =
re", r is the radial distance to the origin and 0 is the angle counterclockwise
subtended by the vector with the positive real axis.
Im{z}
3 jr/6
3
Re(z}
Figure S1.4-1
Introduction / Solutions
S1-3
(i)
(ii)
Im z)
z2 = 4 ej/3
4 9
9
Rez)
Figure S1.4-3
(iii)
Im Iz
jz= 2 ej27/3
3
2 21r
33
Refz)
Figure S1.4-4
(iv)
Indz)
4
9
Figure S1.4-5
Signals and Systems
S1-4
(v)
z*
Itnzl l=
z Rejzj
Figure S1.4-6
(vi)
7T ReWjz
-- 6
3
2
1 3 -jn/6
z 2
Figure S1.4-7
S1.5
This problem shows a useful manipulation. Multiply by e +a/ 2
e-ja/2 = 1, yielding
e+ja/2e -ja/2( 1 - eja) ja/
2
(e -ja/ 2
_ eia/2 )
Now we note that 2j sin (-x) = - 2j sin x = e -x - ex. Therefore,
S1.6
There are three things a linear scaling of the form x(at + b) can do: (i) reverse
direction => a is negative; (ii) stretch or compress the time axis = a I # 1; (iii) time
shifting => b # 0.
Introduction / Solutions
S1-5
x(t + 2)
t
-2
Figure S1.6-2
x(2t + 2)
t
-1 1
2
Figure S1.6-3
x(1 -3t)
t
2 1
3 3
Figure S1.6-4
Signals and Systems
S1-6
S1.7
= 1 -2a
Therefore,
fe 3- dt = 0 + 1e 6 56
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
2 Signals and Systems: Part I
Solutions to
Recommended Problems
S2.1
(a) We need to use the relations w = 21rf, where f is frequency in hertz, and
T = 2w/w, where T is the fundamental period. Thus, T = 1/f.
w /3 1 1
(i) f=-= =-Hz, T- 6s
2r 21 6 f
3r/4 3 8
(ii) f = =-Hz, T=-s
2w 8 3
3/4 3 8-x
(iii) f = = Hz, T = -s
3
2r 87r
Note that the frequency and period are independent of the delay r, and the
phase 0_.
(b) We first simplify:
cos(w(t + r) + 0) = cos(wt + wr + 0)
Note that wT + 0 could also be considered a phase term for a delay of zero.
Thus, if w, = w, and wX, + 0, = wor, + 6, + 2xk for any integer k, y(t) = x(t)
for all t.
(i) Wx = o ~ mr+Ox=
W,
2
W, Wyr, +O ,= 3 1 - 3 =O0+ 2wk
S2.2
(a) To find the period of a discrete-time signal is more complicated. We need the
smallest N such that UN = 21k for some integer k > 0.
(i) N = 2wk =* N = 6, k = 1
3
(ii) N = 2rk =o N = 8, k = 2
4
(iii) 2N = 2wk => There is no N such that aN = 2wk, so x[n] is not periodic.
(b) For discrete-time signals, if Ox = Q, + 2rk and Qxrx + Ox = Qr, + O, + 2k, then
x[n] = y[n].
S2-1
Signals and Systems
S2-2
S2.3
x[n-2]
-1 0 1 2 3 4 5 6
Figure S2.3-1
(ii) x[4 - n] = x[-(n - 4)], so we flip about the n = 0 axis and then shift
to the right by 4.
x[4-n]
0
0 2~l
-1 0 1
111
2 3 4 5 6
n
Figure S2.3-2
(iii) x[2n] generates a new signal with x[n] for even values of n.
x [2n]
0123
Figure S2.3-3
2 0 n
(b) The difficulty arises when we try to evaluate x[n/2] at n = 1, for example (or
generally for n an odd integer). Since x[i] is not defined, the signal x[n/2] does
not exist.
S2.4
By definition a signal is even if and only if x(t) = x(-t) or x[n] = x[-n], while a
signal is odd if and only if x(t) = -x(- t) or x[n] = -x[-n].
(a) Since x(t) is symmetric about t = 0, x(t) is even.
(b) It is readily seen that x(t) # x(- t) for all t, and x(t) K -x(- t) for all t; thus
x(t) is neither even nor odd.
(c) Since x(t) = -x(- t), x(t) is odd in this case.
Signals and Systems: Part I / Solutions
S2-3
(d) Here x[n] seems like an odd signal at first glance. However, note that x[n] =
-x[-n] evaluated at n = 0 implies that x[O] = -x[O] or x[O] = 0. The analo
gous result applies to continuous-time signals. The signal is therefore neither
even nor odd.
(e) In similar manner to part (a), we deduce that x[n] is even.
(f) x[n] is odd.
S2.5
(a) Let Ev{x[n]} = x[n] and Od{x[n]} = x[n]. Since xe[n] = y[n] for n >_ 0 and
xe[n] = x[ -fn], x,[n] must be as shown in Figure S2.5-1.
Xe[n] 2
41 4
n
-5 -4 -3 -2 -1 0 1 2 3 4 5
Figure S2.5-1
Since x[n] = y[n] for n < 0 and x[n] = -x,[-n], along with the property that
x0 [O] = 0, x[n] is as shown in Figure S2.5-2.
Xo [n]
1 2 3 4
0 n
-4 -3 -2 -1 0
Figure S2.5-2
Finally, from the definition of Ev{x[n]} and Od{x[n]}, we see that x[n] = x,[n] +
x[n]. Thus, x[n] is as shown in Figure S2.5-3.
Signals and Systems
S2-4
(b) In order for w[n] to equal 0 for n < 0, Od{w[n]} must be given as in Figure
S2.5-4.
Odfw[n]} 14
-4 - 3 --2 --
n
0 1 2 3 4
Figure S2.5-4
w[n] 21
p e - 0 n
-3 -2 -1 0 1 2 3 4
Figure S2.5-5
S2.6
x[n] 1
Figure S2.6-1
(b) We need to find a # such that e#' = (-e-')". Expressing -1 as ei", we find
e on = (ej'e -)T or 0 = -1 + jr
Note that any # = -1 + jfr + j27rk for k an integer will also satisfy the preced
ing equation.
Signals and Systems: Part I / Solutions
S2-5
Since cos 7rn = (- 1)' and sin 7rn = 0, Re{x(t)) and Im{y(t)} for t an integer are
shown in Figures S2.6-2 and S2.6-3, respectively.
1
Refe (- + ir)n
-1 1In n
0 -e
-e
Figure S2.6-2
Imfe (- +j7T)n I
0 0 0 0 0 -nf
-2 -1 0 1 2
Figure S2.6-3
S2.7
First we use the relation (1 + j) = \/Tej"!' to yield
x(t ) = \/ \/ej'4eJ"/4e(-I+j2*)t = 2eir/2e(-1+j2 )t
Re{x (t)}
envelope is 2e-r
2- -
Figure S2.7-1
Signals and Systems
S2-6
ImIx (t)}
envelope is 2et
2 -
Figure S2.7-2
(c) Note that x(t + 2) + x*(t + 2) = 2Re{x(t + 2)}. So the signal is a shifted
version of the signal in part (a).
x(t + 2) + x*(t + 2)
, 4e 2
Figure S2.7-3
S2.8
(a) We just need to recognize that a = 3/a and C = 2 and use the formula for SN,
N = 6.
6
= 2 =32
a /3\ _
-a)
(3)a
since - 1
- = o-4
S2.9
(a) The sum x(t) + y(t) will be periodic if there exist integers n and k such that
nT1 = kT 2, that is, if x(t) and y(t) have a common (possibly not fundamental)
period. The fundamental period of the combined signal will be nT1 for the small
est allowable n.
(b) Similarly, x[n] + y[n] will be periodic if there exist integers n and k such that
nN = kN 2. But such integers always exist, a trivial example being n = N 2 and
k = N1 . So the sum is always periodic with period nN, for n the smallest allow
able integer.
(c) We first decompose x(t) and y(t) into sums of exponentials. Thus,
13 1 ei(1 6Tt/ 3) -j(16irt/3)
x(t) = 1 ej( 1t/ ) + -e j(21rt/3) + e
2
e
__A6r3 ____
2 2 j
Y( jrt -e-jrt
23 2j
Multiplying x(t) and y(t), we get
3
z(t) - e ( -+/3 - e-j5w/s)
7 )t( + / )t -j 1/3
4j 4j 4j 4j
1 ej(19r/3 )t + 1 ej(13 r/3) t + 1 e j(13r/3)t - -j(19/3)t
2 2 2 2
We see that all complex exponentials are powers of e j(/3). Thus, the funda
mental period is 2 7r/(7r/3) = 6 s.
S2.10
M=Z-QO
x[mI
M= Sx
-OO
[-m] = -
since x[m] is odd. But the preceding sum equals -S. Thus, S = -S, or S = 0.
(b) Let y[n] = x 1[n]x2[n]. Then y[-n] = x1[-n]x 2[-n]. But x1[-nj = -xl[n] and
x2[-n] = x 2[n]. Thus, y[-n] = -x1[n]x2[n] = -y[n]. So y[n] is odd.
(c) Recall that x[n] = x[fn] + x[n]. Then
But from part (b), x,[nxo[n] is an odd signal. Thus, using part (a) we find that
the second sum is zero, proving the assertion.
Signals and Systems
S2-8
while 2r x(t)xo(t) dt = 0
S2.11
(a) x[n] = ewonT - ei2 nTTo. For x[n] = x[n + N], we need
x[n + N] = ej 2x(n +N)T/To eji[ 2 n(T/To) + 2rN(T/To)] = ej2nT/To
The two sides of the equation will be equal only if 27rN(T/TO) = 27rk for some
integer k. Therefore, TITO must be a rational number.
(b) The fundamental period of x[n] is the smallest N such that N(T/TO) = N(p/q)
= k. The smallest N such that Np has a divisor q is the least common multiple
(LCM) of p and q, divided by p. Thus,
N LCM(p, q);
note that k = LCM(p, q)
p q
The fundamental frequency is 2ir/N, but n = (kT 0)/T. Thus,
Q= 2 = = WT = q T
N kT,, k LCM(p, q) 0
(c) We need to find a value of m such that x[n + N] = x(nT + mT). Therefore,
N = m(T./T), where m(T./T) must be an integer, or m(q/p) must be an integer.
Thus, mq = LCM(p, q), m = LCM(p, q)/q.
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
3 Signals and Systems: Part II
Solutions to
Recommended Problems
S3.1
(a)
x[n]= 8[n] + 8 [n - 3]
n
0 1 2 3
Figure S3.1-1
(b)
x[n] = u[n]-u[n-- 5]
1111
0041T
-1 0 1 2 3 4 5
Figure S3.1-2
(c)
x [n]=S[n] + -1 [n -- +(_1)2S[n] +(_1)3 5[n -3]
-3 -2 -1
1
0 1 2 3 4 5 6 7
n
Figure S3.1-3
(d)
x(t)= u(t +3) - u(t - 3)
t
-3 0 3
Figure S3.1-4
S3-1
Signals and Systems
S3-2
(e)
x(t) =6(t + 2)
t
-2 0
Figure S3.1-5
(f)
S3.2
(1) h
(2) d
(3) b
(4) e
(5) a, f
(6) None
S3.3
(a) x[n] = b[n - 1] - 26[n - 2] + 36[n - 3] - 26[n 4] + b[n - 5]
(b) s[n] = -u[n + 3] + 4u[n + 1] - 4u[n 2] + u[n - 4]
S3.4
x (-t)
-12
Figure S3.4-2
x(1--t)
x1-0
-11 1
Figure S3.4-3
-1 0 1 2
t
Figure S3.4-4
6
t
-1 1
Figure S3.4-5
Signals and Systems
S3-4
2
3
t
Figure S3.4-6
2
3
Figure S3.4-7
S3.5
(a) and (b) are equivalent. (e) and (f) are equivalent.
S3.6
Memoryless:
(a) y(t) = (2 + sin t)x(t) is memoryless because y(t) depends only on x(t) and not
on prior values of x(t).
(d) y[n] = Ek=. x[n] is not memoryless because y[n] does depend on values of
x[-] before the time instant n.
(f) y[n] = max{x[n], x[n - 1], ... , x[-oo]} is clearly not memoryless.
Linear:
(a) y(t) = (2 + sin t)x(t) = T[x(t),
T[ax 1 (t) + bx 2 (t)] = (2 + sin t)[axi(t) + bx 2 (tt)
= a(2 + sin t)x 1(t) + b(2 + sin t)x 2 (t)
= aT[x 1 (t)] + bT[x 2 (t)]
Therefore, y(t) = (2 + sin t)x(t) is linear.
(b) y(t) = x(2t) = T[x(t)],
T[ax 1(t) + bx 2 (t)] = ax1 (2t) + bx 2 (t)
= aT[x1 (t)) + bT[x 2 (t)]
Therefore, y(t) = x(2t) is linear.
dxt
d
T[ax 1(t) + bx 2(t)] = -[ax 1 (t)+±bx
2 (t)]
dx 1 (t) dx (t)
= a dt + b dt2 = aT[x 1 (t)] + bT[x 2 (t)]
Time-invariant:
(a) y(t) = (2 + sin t)x(t) = T[x(t)],
T[x(t - T 0 )] = (2 + sin t)x(t - T0 )
9 y(t - T0 ) = (2 + sin (t - T0 ))x(t - T0 )
Therefore, y(t) = (2 + sin t)x(t) is not time-invariant.
(b) y(t) = x(2t) = T[x(t)],
T[x(t - # x(2t - TO) = y(t - T 0 )
T0 )] = x(2t - 2T0)
Therefore, y(t) = x(2t) is not time-invariant.
(d) y[n] = E
k=
x[k] = T[x[n]],
n n-NO
Yes, y[n] = E .x[k] is causal because the value of y[-] at any instant n
depends only on the previous (past) values of x[-].
Invertible:
(b) y(t) = x(2t) is invertible; x(t) = y(t/2).
(c) y[n] = E _.x[k] is not invertible. Summation is not generally an invertible
operation.
(e) y(t) = dx(t)/dt is invertible to within a constant.
Stable:
(a) If Ix(t) I < M, Iy(t) I < (2 + sin t)M. Therefore, y(t) = (2 + sin t)x(t) is stable.
(b) If |x(t)| < M, |x(2t)I < M and ly(t)| < M. Therefore, y(t) = x(2t) is stable.
(d) If |x[k]| 5 M, ly[n]j 5 M - E_,, which is unbounded. Therefore, y[n] =
E"Lx[k] is not stable.
Signals and Systems: Part II / Solutions
S3-7
S3.7
H~':
dx(t)
y(t) = d
dt
G: y(t) = x(2t)
G 1: y(t) = x(t/2)
(b)
Solutions to
Optional Problems
S3.8
(a) x 2(t) = xi(t) - xi(t - 2)
y)(t - 2)
2(t) = 1 (t)
Figure S3.8-1
Signals and Systems
S3-8
y 3(t)=y (t)+y 1 (t + 1)
-1 10 1 2
Figure S3.8-2
u(t
-2)
-u(1
-t)
y(t)=e-(t-1)u
-1)+u(-t)+
(c) x(t) = u(t - 1) - u(t - 2)
Figure S3.8-3
3 3
p p14
p n
-3 -2 -1 0 1 2 3
-4
Figure S3.8-4
Signals and Systems: Part 11/ Solutions
S3-9
y2 [n] 2 __
0 1 2 3 4 5
Figure S3.8-5
Y3[n] = y 1 [n + 1]
-1 0 234
Figure S3.8-6
S3.9
(a) (i) The system is linear because
Tlaxi(t) + bx 2 (t)] = 1
n=
23 [ax (t) + bx (t )](t - nT)
x2 (t) = sin (
+r =
cos (2)
Now the output
+00
= E cos(2,rt)b(t - nT)
cos(27rt)
0< t
Figure S3.9-1
T=1
0 tt
I ~11
2
Figure S3.9-2
Signals and Systems: Part II / Solutions
S3-11
y (t )
T=
T~__ t t
-1
Figure S3.9-5
y (t)
Figure S3.9-6
etcos(2nrt)
et
t
22
Figure S3.9-7
y (t)
T= I
T=1 tt
1 3
2
Figure S3.9-8
Signals and Systems
S3-12
y (t)
T= Ie3
2 2
2 I 2- t
-e1/2 I I 3/2
Figure S3.9-9
y (t)
T= YWe3
4t
F r2 2
-el2 -e 3 12
Figure S3.9-10
T= Y3 1 y~t)
12 2
1 2
Figure S3.9-12
S3.10
S3.11
(a) y[n] = x[n] + x[n - 11 = T[x[n]]. The system is linear because
T[ax 1 [n] + bx 2 [n]| = ax1[n] + ax1[n - 1] + bx 2[n] + bx 2[n - 1]
= aT[x1 [n]] + bT[x 2[n - 1]]
The system is time-invariant because
y[n] = x[n] + x[n - 1] = Tjx[n]],
T[x[n - N]] = x[n - N] + x[n 1 - N]
= y[n - N]
(b) The system is linear, shown by similar steps to those in part (a). It is not
time-invariant because
T[x[n N]] = x[n - N] + x[n - N - 1] + x[O]
# y[n - N] = x[n N] + x[n - N - 1] + x[-NJ
S3.12
(a) To show that causality implies the statement, suppose
x1 (t) - yl(t) (input x 1(t) results in output y 1(t)),
x 2(t) - y2),
where y 1(t) and y 2(t) depend on x 1(t) and x 2(t) for t < to. By linearity,
xI(t) - x 2(t) -+ y 1 (t) - y 2(t)
Signals and Systems
S3-14
If Xi(t) = x 2 (t) for t < to, then y 1 (t) = y 2 (t) for t < to. Hence, if x(t) = 0 for
t < to, y(t) = 0 for t < to.
(b) y(t) = x(t)x(t + 1),
x(t) = 0 for t < to =* y(t) = 0, for t < to
This is a nonlinear, noncausal system.
(c) y(t) = x(t) + 1 is a nonlinear, causal system.
(d) We want to show the equivalence of the following two statements:
Statement 1 (S1): The system is invertible.
Statement 2 (S2): The only input that produces the output y[n] = 0 for all n is
x[n] = 0 for all n.
To show the equivalence, we will show that
S2 false S1 false and
S1 false S2 false
S2 false == S1 false: Let x[n] # 0 produce y[n] = 0. Then cx[n] == y[n] = 0.
S1 false S2 false: Let xi => yi and x 2 =* Y2. If x 1 # X 2 but y1 = Y2, then
X1 - X2 0 0 but yi - yi = 0.
(e) y[n] = x 2 [n] is nonlinear and not invertible.
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
4 Convolution
Solutions to
Recommended Problems
S4.1
The given input in Figure S4.1-1 can be expressed as linear combinations of xi[n],
x 2[n], X3 [n].
x,[ n]
0 2
Figure S4.1-1
(b) Using superposition, y 4[n] = 2yi[n] - 2y 2[n] + y3 [n], shown in Figure S4.1-2.
-1 0 1
Figure S4.1-2
(c) The system is not time-invariant because an input xi[n] + xi[n - 1] does not
produce an output yi[n] + yi[n - 1]. The input x,[n] + xi[n - 11 is xi[n] +
xi[n - 1] = x2[n] (shown in Figure S4.1-3), which we are told produces y 2[n].
Since y 2[n] # yi[n] + yi[n - 1], this system is not time-invariant.
n
0 1
Figure S4.1-3
S4-1
Signals and Systems
S4-2
S4.2
The required convolutions are most easily done graphically by reflecting x[n] about
the origin and shifting the reflected signal.
(a) By reflecting x[n] about the origin, shifting, multiplying, and adding, we see
that y[n] = x[n] * h[n] is as shown in Figure S4.2-1.
(b) By reflecting x[n] about the origin, shifting, multiplying, and adding, we see
that y[n] = x[n] * h[n] is as shown in Figure S4.2-2.
y[n]
3
2
0 1 2 3 4 5 6
Figure S4.2-2
S4.3
(a) It is easiest to perform this convolution graphically. The result is shown in Fig
ure S4.3-1.
Convolution / Solutions
S4-3
4-
| t
4 8
Figure S4.3-1
(b) The convolution can be evaluated by using the convolution formula. The limits
can be verified by graphically visualizing the convolution.
1
t+
e (- dr, t> 0,
-0, t < 0,
Let r' = T - 1. Then
y( ) e- d r -e t > 0,
0 0 , t < 0
y(t)
It- I / \, t
1 3 5
Figure S4.3-2
Signals and Systems
S4-4
S4.4
y[n] = b[m
6 - no]h[n - m] = h[n - no]
m= -oO
y [n] = h1[12 I
ae|41 8 n
(n 1) no (n1+ 1)
Figure S4.4-1
For n > 0: 1 +
y[n] = = 2( 1
,
y[n] = 2 - (i)"
y[n]
2--
140
0 1 2
Figure S4.4-2
(c) Reversing the role of the system and the input has no effect on the output
because
S4.5
(a) (i) Using the formula for convolution, we have
y 1 (t) = x(r)h(t - r) dr
= r(-)-2 u(t - r) dr
t
e -( 2
= - dr, t > 0,
2e 10 = 2(1 e t > 0,
y(t) = 0, t< 0
y 1 (t)
-
2 -
t
0
Figure S4.5-1
=4(1 - e-t/2), 3 t : 0,
y 2 (t) = { 2e-(- /2
3
d-, t >_3,
= 4e- t/ 2 (e'/ 2 t 3,
- 1),
y 2(t) = 0, t 0
y 2 (t)
0 1 2 3 4 5 6
Figure S4.5-2
Signals and Systems
S4-6
(b) Since x 2 (t) = 2[xl(t) - xl(t - 3)] and the system is linear and time-invariant,
y 2(t) = 2[yi(t) - y1(t - 3)].
Solutions to
Optional Problems
S4.6
(a)
x(T)
T
0 1
Figure S4.6-1
h(-1 -r)
2,
-2 T
--1 0
--1 '
Figure S4.6-2
Convolution / Solutions
S4-7
h(0-r)
Figure S4.6-3
h(1 -- )
Figure S4.6-4
Signals and Systems
S4-8
Using these curves, we see that since y(t) = x(t) * h(t), y(t) is as shown in
Figure S4.6-6.
y(t)
1 W
--1 2
t
0 1
--- 4
Figure S4.6-6
h(r)
0 1 2
Figure S4.6-7
x(t -r)
Figure S4.6-8
x(t )
Figure S4.6-9
Convolution / Solutions
S4-9
x(t- r)
Figure S4.6-10
x(t t)
Figure S4.6-11
y(t)
3 -
I-
U' 1 2 3 4
Figure S4.6-12
S4.7
y[n] = x[n] * h[n]
=1 x[n - m]h[m]
nO
=- anmm,
- n > 0,
M=0
Signals and Systems
S4-10
y[n] = 0, n < 0
S4.8
(a) x(t) = E_= - kT) is a series of impulses spaced T apart.
x(t)
t
-2T -T 0 T 2T
Figure S4.8-1
y(t)
.. t
-2 3 -1 0 1 1 3 2
2 2 2 2
Figure S4.8-2
y(t)
-2 -3 -1 -j 0 1 1 3 2
22 2 2
Figure S4.8-3
Convolution / Solutions
S4-11
S4.9
(a) False. Counterexample: Let g[n] = b[n]. Then
x[n] * {h[n]g[n]} = x[n] h[0],
{x[n]* h[n]}g[n] = b[n] [x[n] *h[n]]
n=0
2
Let r' = T/ . Then
(c) True.
y(t) = x(t) * h(t)
y(- t) = x(-t) * h(-t)
S4.10
(a)
9(t) = Jro TO2 1 (r)22(t - r) dr,
O
(b) a+TO
TO+b
=
FTo
T1()- 2 t - r) dr +
T7 0 +b
1&)2(t - r) di
b TO
r) d1 2- ) di
=
=
T)TO Tb
21
t(-r
()A
)12( t T )
-
)
dr =q t)
(c) For 0 s t - I
9(t ) = ft
0e
e- di +
+
Ti±e-'di I
e/2+t1/2+1
=(-e-' t + (-e-T
0 11/2+t)
For s t < 1:
= ft
t(t) e-' di = e- ( 1/2) - e
- 1/2 )e
2 X1[n] *x2[n]
19
0 1 3 4 5
Figure S4.10-1
Convolution / Solutions
S4-13
(e)
S4.11
(a) Since y(t) = x(t) *h (t) and x(t) = g(t) *y(t), then g(t) *h(t) = 6(t). But
n F, n = 0,
k=0 n O
Therefore,
go = 1/h 0 ,
gi = -- hi/ho,
1 (-hl h2)
2 o ho ho
Therefore,
h(t) = (t -kT )
k=O
<
k=O
a f -w
_(T - kT)Ix(t - T)I dr
y(t) = ( aku(t-kT)
k=O
y(NT) = Z ak
k=O
Delay T
Figure S4.11
Convolution / Solutions
S4-15
S4.12
then
y[n] = Z akr[n - k]
k= -w
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
5 Properties of Linear,
Time-Invariant Systems
Solutions to
Recommended Problems
S5.1
The inverse system for a continuous-time accumulation (or integration) is a differ
entiator. This can be verified because
d[ x(r) dr =x(t)
Therefore, the input-output relation for the inverse system in Figure S5.1 is
dy(t)
x(t ) = dt
Figure S5.1
S5.2
(a) We want to show that
h[n] - ah[n - 11 = b[n]
Substituting h[n] = anu[n], we have
anu[n] - aa" -u[n - 1] = a"(u[n] - u[n - 1])
But
u[n] - u[n - 1] = b[n] and an6[n] = ab[n] = b[n]
(b) (i) The system is not memoryless since h[n] # kb[n].
(ii) The system is causal since h[n] = 0 for n < 0.
(iii) The system is stable for Ia I < 1 since
1
|al" 1 - |a
is bounded.
(c) The system is not stable for |a l > 1 since E"o 1a|" is not finite.
S5.3
(a) Consider x(t) = 6(t) -+ y(t) = h(t). We want to verify that h(t) = e -2u(t), so
dy(t) - -2e- 2
u(t) + e -' 6(t), or
dt
dy(t) + 2y(t) =e -21 6(t),
dt
S5-1
Signals and Systems
S5-2
but ed' b(t) = 6(t) because both functions have the same effect on a test func
tion within an integral. Therefore, the impulse response is verified to be correct.
(b) (i) The system is not memoryless since h(t) # kb(t).
(ii) The system is causal since h(t) = 0 for t < 0.
(iii) The system is stable since h(t) is absolutely integrable.
2t
IhtI dt = e dt = -le -2t
S5.4
By using the commutative property of convolution we can exchange the two systems
to yield the system in Figure S5.4.
Figure S5.4
h(t) = ds(t)
dt
S5.5
(a) By definition, an inverse system cascaded with the original system is the iden
tity system, which has an impulse response h(t) = 6(t). Therefore, if the cas
caded system has an input of b(t), the output w(t) = h(t) = 6(t).
(b) Because the system is an identity system, an input of x(t) produces an output
w(t) = x(t).
Solutions to
Optional Problems
S5.6
(a) If y(t) = ayi(t) + by 2(t), we know that since system A is linear, x(t) = ax,(t)
+ bx 2 (t). Since the cascaded system is an identity system, the output w(t) =
ax 1(t) + bx 2 (t).
Properties of Linear, Time-Invariant Systems / Solutions
S5-3
Figure S5.6-1
(b) If y(t) = y 1 (t - r), then since system A is time-invariant, x(t) = x,(t - -) and
also w(t) = xi(t - r).
Figure S5.6-2
(c) From the solutions to parts (a) and (b), we see that system B is linear and time-
invariant.
S5.7
(a) The following signals are obtained by addition and graphical convolution:
(x[n] + w[n]) * y[n] (see Figure S5.7-1)
x[n] * y[n] + w[n] * y[n] (see Figure S5.7-2)
x[n] +w[n]
yInI
1?
0 n -
0 1 -1 6
-1
(x[n] +w[n]) *y[n]
39
1 2
-l 0
-3 6
Figure S5.7-1
Signals and Systems
S5-4
x[n] *y[n]
2 2
1 0 n -0 1 9 n
-1 0 -1 0
-2
1 2
-I 0
-2<
-3
Figure S5.7-2
2 (x[n] *y[n])-w[n]
1 1 2
n 1 0 n
-1 0 0 3
-1 -1
-2
Figure S5.7-3
S5.8
Consider
= x(r)h(t - r) dr
S5.9
(a) True.
f0 |h(t)| dt =
f T |h(t )|
(b) False. If h(t) = 3(t - to) for to > 0, then the inverse system impulse response
dt = oo
Z
n= -oo
Ih[n]| = [
n= -oo
u[n] = c
n=-co
Zh[n]I = E
n=-K
Ih[n] = M (a number)
(e) False. h(t) = u(t) implies causality, but J u(t) dt = oo implies that the sys
tem is not stable.
Signals and Systems
S5-6
(f) False.
hi(t) = 6(t - ti), ti > 0 Causal
h 2(t) = b(t + 2 ), t2 > 0 Noncausal
h(t) = hI(t) * h 2(t) = 6(t + t 2 - ti), t2 !- t 1 Causal
(g) False. Suppose h(t) = e-'u(t). Then
(1 - e -') dt =t + e-'1 = oo
(h) True. We know that u[n] = E=O b[n - k] and, from superposition, s[n] =
Ef=0 h[n - k]. If s[n] # 0 for some n < 0, there exists some value of h[k] # 0
for some k < 0. If s[n] = 0 for all n < 0, h[k] = 0 for all k < 0.
S5.10
dx(t - r) dx(t)
dt r=O dt
(b) - [gt)ft)]_
g(t)f(t)uj(t) dt dt =
= -[g'(t)ftt) ± g(t)f'(t)]__
- IgMAO+ Xt 1(l It =0
= -[g'(O)f(O) + g(0)f'(0)],
(c)
-- o
x(r)u 2(r) dr = x(r))
_Y6
- f _.dr
-u (r)dr
*~dx dx/ d2 X
-J -uir) dr = 7 ) + d2U0(r) dr
2
d x
dr2 r=0
S5.11
(a) h(t) * g(t) = J'. h(t - r)g(r) dr = f' h(t - )g(r) dr since h(t) = 0 for t < 0
and g(t) = 0 for t < 0. But if t < 0, this integral is obviously zero. Therefore,
the cascaded system is causal.
(b) By the definition of stability we know that for any bounded input to H, the out
put of H is also bounded. This output is also the input to system G. Since the
input to G is bounded and G is stable, the output of G is bounded. Therefore, a
bounded input to the cascaded system produces a bounded output. Hence, this
system is stable.
S5.12
We have a total system response of
h = {[(h, * h 2 ) + (h 2 * h 2 ) - (h 2 * hl)] * hi + h-i} * hl
h = (h 2 * h) + (hi 1 * h 2 ')
S5.13
We are given that y[n] = x[n] * h[n].
Iy[n]I = E
k=-Oo
x[n - k]h[k]
k=-o
=maxm{x n]l} z:
k= -o
lh[k]|
that E Ih[k]| 5 1 max {Iy[n]|} s max {Ix[n]|}. This means that E Ih[k]I
- 1 is a sufficient condition. It is necessary because some x[n] always exists that
yields y[n] = E=_ Ih[k]1. (x[n] consists of a sequence of +1's and -l's.) There
fore, since max {x[n]} = 1, it is necessary that E= h[k] 1 5 1 to ensure that y[n]
S max {Ix[n]I} = 1.
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
6 Systems Represented by
Differential and Difference
Equations
Solutions to
Recommended Problems
S6.1
We substitute ya(t) = ay1(t) + Oy2(t) into the homogeneous differential equation
dya(t)
dt + ay3(t)= dd [ay 1(t) + #y2(t)] + a[ayi(t) + #y 2(t)]
dt dt
Since differentiation is distributive, we can express the preceding equation as
ady t) +0dy2t) + aay1(t) + afly2(t)
ddty ddty
aa _ di(
dt + ay1(t) I + dt
dt + ay2(t)
ay()
However, since both y 1(t) and y 2(t) satisfy the homogeneous differential equation,
the right side of the equation is zero. Therefore,
dy3 (t) + ay3(t)
= 0
cit
S6.2
(a) We are assuming that y(t) = es'. Substituting in the differential equation yields
d2 d
d (es') + 3 (e t ) + 2e" = 0
so that
s2e"t + 3se' + 2e" = es'(s 2 + 3s + 2) = 0
For any finite s, e" is not zero. Therefore, s must satisfy
0 = s2 + 3s + 2 = (s + 1)(s + 2), s = -1, -2
(b) From the answer to part (a), we know that both yi(t) = e-' and y 2 (t) = e-2
satisfy the homogeneous LCCDE. Therefore,
ys(t) = Kie~' + K 2e 2,
S6.3
(a) Assuming y(t) of the form
y(t) = Kest ,
we substitute into the LCCDE, setting x[n] = 0:
dy(t) 1 1 s
0= + y(t) =Ksel+K-7es = Ket±
dt 2 22
S6-1
Signals and Systems
S6-2
Since K # 0 and e" # 0, s must equal --. K then becomes arbitrary, so the
family of y(t) that satisfies the homogeneous equation is
y(t) =Ke-12
(b) Substituting into eq. (P6.3-1) y 1 (t) = Ae' for t > 0, we find
dy1 (t) 1 1 ~ t
t >0
dy I 1y t)= -Ae-'+ Ae'=e-'
Since e-' never equals zero, we can divide it out. This gives us an equation
for A,
A
asA = -2
2
(c) For y 1(t) = (2e -t/ 2 - 2e- t)u(t),
dyi(t) [2 (- -) e -1/2
- 2(-1)e -'], t > 0
dt t0, t <0,
dyi(t) 1 F(-e 1
t 2
+ 2e -t) + 1 (2e-' 2
- 2e- t) = e-', t> 0
dt 2 0, t < 0
= x(t)
S6.4
(a) Note that since y[n] is delayed by one sample by the delay element, we can label
the block diagram as shown in Figure S6.4.
x [n] - 1y [n[-1
_1
2
Figure S6.4
Thus y[n] = x[n] - 1],
- +y[n or y[n] + ly[n - 1] = x[n].
(b) Since the system is assumed to be causal, y[n] must be zero before a nonzero
input is applied. Therefore, x[n] = 0 for n < 0, and consequently y[n] must be
zero for n <0. Thus, y[-5] = 0.
(c) Since x[n] = [n] = 0 for n < 0, y[n] must also equal zero for n <0. For n = 0,
we have y[0] + ly[-1] = 1 or, substituting for y[n],
Ka~u[0] + Ka-lu[-1] = 1,
K + i- 0 = 1, or K=1
For n > 0, we have
y[n] + ly[n - 1] = 0 or an + ia"-1 = 0
Systems Represented by Differential and Difference Equations / Solutions
S6-3
since K = 1. Thus, a must equal -- for a" + -a1 to equal 0 for all n > 0.
Therefore, y[n] = (-1)"u[n]. Substituting into the left side of the difference
equation, we have
(-2)" u[n] + (1)"-1 u[n-1] = (-2)" u[n] - (-2)" u[n-1]
n=0
1,
-0, otherwise
(d) We can successively calculate y[n] by noting that y[ -1] = 0 and that
y[n] = -1y[n - 1] + b[n]
n=0, y[0] = -- - 0 + 1 = 1
n = 1, y[1] = - 1 + 0 = -1
n = 2, y[2] = -i (-)+0 =i
S6.5
(a) Performing the manipulations in inverse order to that done in the lecture (see
Figure S6.5-1) yields the system shown in Figure S6.5-2.
D
11
3
Figure S6.5-1
x [n] y [n]
-l__
_ - - - - --- L _ -
A B
Figure S6.5-2
Signals and Systems
S6-4
Since the system is linear and time-invariant, we can exchange the order of the
two boxes A and B, yielding the direct form I shown in Figure S6.5-3.
x [n] + y [n]
D | |D
x [n -1 -2 1y[n-1
B A
Figure S6.5-3
(b) From the direct form I, we see that the intermediate variable q[n] is related to
x[n] by
q[n] = x[n] - 2x[n - 1]
The signal y[n] can be described in terms of q[n] and y[n - 11 as
y[n] = q[n] + iy[n - 1]
Combining the two equations yields
y[n] = iy[n - 1] + x[n] - 2x[n - 11, or
y[n] - iy[n - 11 = x[n] - 2x[n - 1]
(c) (i) Figure S6.5-4 shows that if we concentrate on the right half of the dia
gram of direct form II given in Figure P6.5, we see the relation
y~n] = r[n] - 2r[n - 11
r [n] y[n]
r[n- 1] -2
Figure S6.5-4
(ii) Similarly, Figure S6.5-5 shows that if we concentrate on the first half of
the diagram, we obtain the relation
r[n] = x[n] + ir[n - 1], or x[n] = r[n] - ir[n - 1]
Systems Represented by Differential and Difference Equations / Solutions
S6-5
r [n]
x [n] - +
r[n- 1]
Figure S6.5-5
(iii) From the two equations obtained in parts (i) and (ii),
x[n] = r[n] - ir[n - 11 (S6.5-1)
and
y[n] = r[n] - 2r[n - 1], (S6.5-2)
we solve for r[n],obtaining
r[n] = fx[nJ - ky[n]
Substituting r[n] into eq. (S6.5-1), we have
x[n] = ,x[n] - y[n] - }{5xfn - 1] - iy[n - 1]},
which simplifies to
y[n] - iy[n - 1] = x[n] - 2x[n - 1]
S6.6
(a) Integrating both sides of eq. (P6.6-1) yields
y(t) + afy(t) dt = bx(t) + c fx(t) dt, or
y(t) = -a fy(t) dt + bx(t) + c fx(t) dt
Thus, we set up the direct form I in Figure S6.6-1.
x(t) P. y(t)
+ +
C -a
Figure S6.6-1
Signals and Systems
S6-6
(b) Since we are told that the system is linear and time-invariant, we can inter
change boxes A and B, as shown in Figure S6.6-2.
-'-4
b
-a c
L-_- -- ___._
B
Figure S6.6-2
Combining the two integrators yields the final answer, shown in Figure S6.6-3.
x(t) b + -1y t)
-a C
Figure S6.6-3
Solutions to
Optional Problems
S6.7
(a) In Figure S6.7 we convert the block diagram from Figure P6.7 to direct form I.
D D
-4
Figure S6.7
Systems Represented by Differential and Difference Equations / Solutions
S6-7
q[n] is given by
q[n] = x[n] + x[n - 1]
while
y[n] = q[n] - 4y[n - 1]
Substituting for q[n] yields
y[n] + 4y[n - 11 = x[n] + x[n - 1]
(b) The relation between x[n] and r[n] is r[n] = -4r[n 1] + x[n]. For such a -
simple equation, we solve it recursively when a[n] = x[n].
n b[n] r[n - 1] r[n]
<0 0 0 0
0 1 0 1
1 0 1 -4
2 0 -4 16
3 0 16 -64
We see that r[n] = (-4)"u[n].
(c) y[n] is related to r[n] by
y[n] = r[n] + r[n - 1]
Now y[n] = h[n], the impulse response, when x[n] = 6[n], and
h[n] = (-4)"u[n] + (-4)"-Iu[n - 1]
This expression for h[n] can be further simplified:
h[n] = (-4)"u[n] + (-4)"-Iu[n - 1]
or
0 0'
h[n] =
1, n =0
For n > 0,
h[n] = (-4)" + (-4)"-1
-3(-4)"-1
Thus,
h[n] = b[n] - 3(-4)"1 u[n - 1]
S6.8
Note that the system in Figure P6.8 is not in any standard form. Relating r(t) to x(t)
first, we have
x(t) r(t)
x(t)+r(t) f
Figure S6.8
S6.9
(a) Substituting y[n] = Az" into the homogeneous LCCDE, we have
Az" - -Azo ' = 0
Thus
K K
Y= 1 =\/ - cos 0 e
1
j tan- 1(sin no)/(2-cos go)1 or
cos 0 \/ o
KK
Therefore,
0
y[n] = Re[Ye' u[n]] = a_ Ko~ KK 'On
Retei tan l[(sin Oo)/(2-o
CoSjU~f
K
= B cos( 0 n + 0), where B =
sin g0
0 = -tan-'( o o
S6.10
The important observation to make is that if [dtr(t)]/dti is the input to the system
H, then [dis(t)]/dt' will be the output. Suppose that we construct a signal
M d'r(t)
q(t) = (a
i= 1
S6.11
(a) Substituting y(t) = Aeso' into the homogeneous LCCDE, we have
N dky(t) N dk
ak d =(a k(Ae)=0
k=0 k=0
= aks Aet =0
\k= 0)/
p(so) = ( aks = 0
k=0
d dp(s)
= - [p(s)Ae"] = Atp(s) + A est
ds ds
For s = so, p(so) = 0. Also, since p(s) is of the form
p(s) = (s - so) 2q(s),
we have
dp(s) 0
ds =S=
Therefore, AteO' satisfies the homogeneous LCCDE.
(c) Substituting y(t) = est , we get the characteristic equation
s 2 + 2s + 1 = 0, or sO = -1
Thus, y(t) = Ke~' + K 2te~'. For y(0) = 1 and y'(0) = 1, we need K, = 1 and
K 2 - Ki = 1, or K2 = 2. Thus,
y(t) = e~' + 2te
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
7 Continuous-Time Fourier Series
Solutions to
Recommended Problems
S7.1
(a) For the LTI system indicated in Figure S7.1, the output y(t) is expressed as
LTI
Figure S7.1
y(t) = f h()ew(- T) dr
= ejwt { h(r)e-j' dr
= e+iwtH(w)
(b) We are given that the first-order differential equation is of the form
dy(t) + ay(t) = x(t)
dt
From part (a), when x(t) = eiw', then y(t) = ej"'H(w). Also, by differentiating
y(t), we have
dy(t) = jcoej"'H(co)
dt)
Substituting, we get
jwewtH(w) + aejw'H(w) ewt
Hence,
jwH(o) + aH(w) = 1, or
1
H(w) =
a + jw.
S7.2
(a) The output of a discrete-time LTI system is given by the discrete-time convo
lution sum
y[n] = ( h[k]x[n - k]
k
S7-1
Signals and Systems
S7-2
y[n] = ( h[k]z"-k
k
= z" h[k]z -k
k
= z"H(z)
(b) We are given that the first-order difference equation is of the form
y[n] + ay[n - 1] = x[n]
From part (a), if x[n] = z", then y[n] = z"H(z). Hence,
y[n - 1] = z" -1H(z).
By substitution,
z"H(z) + az"- H(z) =z"
which implies
(1 + az-')H(z) = 1,
1
H(z) = _
1 + az-1
S7.3
2
-_ . _}15_ e -j2xts + ej t6 _ e-j2xt4
2j 2j 4j 4j
+ e j21rt4 ee -j2t6
4j 4j
Continuous-Time Fourier Series / Solutions
S7-3
Therefore,
x(t) = ( ake*4'
2
where wo = 1r.
ejr/6
a4 =- . a-4
43j 4j.
/6
e jr/6 _ e-ir
a5 = -, a-5, - I
2j 2j
e jr/6 _e-jr16
a6 - a_6 = .
4j 4j
All other ak's = 0.
S7.4
(a)
ak = - x(t)e -jk" 0 t dt
6 -jk(7r/3) 2 6 -jk(7r/3) 1
and
cos(27r/3)k - cos(r/3)k
ak = jrk
Signals and Systems
S7-4
ao = (1|T0 )
f x(t)dt.
(b)
x(t)
Figure S7.4-2
3/2
[6(t) 26(t - 1)]eikO' dt
2 1/2 -
1 . 1.
e-k-o = - (e-i-)k
2 2
Therefore,
a o - 2, ak = - ()k
It is instructive to plot ak, which we have done in Figure S7.4-3.
Continuous-Time Fourier Series / Solutions
S7-5
S7.5
(a) (i) and (ii)
From Problem 4.12 of the text (page 260), we have
x {t - T)= -x(t),
which means odd harmonics. Since x(t) is real and even, the waveform has real
coefficients.
(b) (i) and (iii)
-x(t) = x t
which means odd harmonics. Since x(t) is real and odd, the waveform has imag
inary coefficients.
(c) (i)
W~ T)
-x(t) = x t 4 )'
which means odd harmonics. Also, x(t) is neither even nor odd.
Solutions to
Optional Problems
S7.6
x(t) is specified in the interval 0 < t < T/4, as shown in Figure S7.6-1.
x(t)
| t
T T
8 4
Figure S7.6-1
(a) Since x(t) is even, we can extend Figure S7.6-1 as indicated in Figure S7.6-2.
x (t)
t
_T _T T T
4 8 8 4
Figure S7.6-2
Signals and Systems
S7-6
Since x(t) has only odd harmonics, it must have the property that x(t - T/2)
= -x(t), as shown in Figure S7.6-3.
x(t) in
x(t - T) = -x(t)
2 Figure S7.6-2 shifted
T T T 37
4 4 2 4
-x(t) in
Figure S7.6-2
/
Figure S7.6-3
x(t)
_T 0 T T 3T T
4 4 2 4
-1-
Figure S7.6-4
Since x(t) is odd, for - T/4 < t < T/4 it must be as indicated in Figure S7.6-6.
x(t)
T T T
4 8 4
Figure S7.6-6
Since x(t) has odd harmonics, x[t - (T/2)] = -x(t). Consequently x(t) is as
shown in Figure S7.6-7.
x(t)
T T T
8 4
-1
Figure S7.6-7
S7.7
1
ak = - x(t)e -jkw 0t dt
TO TO
Substituting r = - t, we have
1
ak = - = a-k
T0OfT 0 x(T)e kwOdr
Signals and Systems
S7-8
x*(t)e -jk 0t dt
50TO 0
(c) = -
= T0 x(t)ekwot dt = ak,
6k = ak
ak = - x(r)e -jk(2"/To)' dr = ak
TO To
Therefore,
S7.8
(a) Since O4(t) are eigenfunctions and the system is linear, the output is
y(t) = L
k=-00
XkCk'k(t).
2
d 2x(t)
2
dx(t)
(b) y(t)=t dt +t
dt2 dt
kk(t) = tk,
d0k(t) ktk1_
dt
d2dtk(t) "k-Jt2 '
2 ktk 1)tk-2
dt = -
S7.9
(a) 9(t) = 2 1(t) * 2 2(t)
= fT I 1(r)A2 (t - r) dr
= --
1
J 1 (r)e -<k(2./T 0)Tr
F
- r)e jk(2r T 0)(t-) dt
Toakbk
Continuous-Time Fourier Series / Solutions
S7-9
(b) Since z(t) * z(t) =x(t), as shown in Figure S7.9-1, then 2(t) is shown in Figure
S7.9-2.
x ( t)
2
z(t) z(t)
t t t
-1 1 -1 1 -2 2
Figure S7.9-1
(c) Without explicitly carrying out the convolutions, we can argue that the aperi
odic convolution of x 1(t) and x 2 (t) will be symmetric about the origin and is
nonzero from t = -2T to t = 2T. Now, if 1 1(t) and 2 2(t) are periodic with period
To, then the periodic convolution, P(t), will be periodic with period To. If To is
large enough, then q(t) is the periodic version of y(t) with period To. Hence, to
recover y(t) from P(t) we should extract only one period of P(t) from t =
-TO/2 to t = TO/2 and set y(t) = 0 for ItI > TO/2, where TO/2 2T, or To -
4T.
S7.10
(a) The approximation is
N
N(t)= akk(t
k=-N
Hence,
k = 1,
a 1,
otherwise,
we get b=
- T ak X*(tk(t) dt + T |akI
k a k
Since ai = bi + jc,
=E r OC x*(t)Oi(t) dt + 2bi
an d x(t)k7(t) dt
and
-= 0 and = 0,
we can multiply the second equation byj and add the two equations to get
(E aE
ab + - = 0
By substitution, we get
bi + jci = ax(t)O(t) dt
= a,
(b) If {4(t)} are orthogonal but not orthonormal, then the only thing that changes
from the result of part (a) is
Tab a
k1
qk*~c(t) dt Z Iak IAkk 2 k
ab
a = fbx(t)4it) dt
Continuous-Time Fourier Series / Solutions
S7-11
(c) Since
for all values of a, using parts (a) and (b) we can write
ai = a x(t)e-"nwo' dt
_ 1
=- x(t)e-nwot dt
T0 fTo
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
8 Continuous-Time
Fourier Transform
Solutions to
Recommended Problems
S8.1
(a)
x(t)
t
Tj Tj
2 2
Figure S8.1-1
i(t)
t
3T 1 -- T1 To T1 T1 To Tl 3
T 1 =O
2 2 2 2 2 2
Figure S8.1-2
X(w) = ox(t)e -j dt =
Ti/2 le-j" dt since x(t) = 0 for ItI>
S8-1
Signals and Systems
S8-2
ak = T f X(t)e ~jk 0
t dt,
S8.2
0 - (a+jw)t dt -1 e- (a+jo)t
a+jw o
X(w) + X*(O) a
Re{X(w)} = 2 a +
2
(2 i
= X(W) - X*(W) -W
Im{X(W)} 2 a2 + w 2 '>
Im{X(W)) _ w
-X(O) = tan- = -tan a
IX(w)I
-a a
4X(W)
2
Tro
-aa
2T
Figure S8.2-2
1 0
[ 1-+j(2-)t
1 + j(2 - )0
X(w) + X*(W) 1
Re{X()} =
2 1 + (w - 2)2
X(w) - X*(co) -(o - 2)
Im{X(W)} =
+ ((A -2)2
2 1
X(w)I
1 2 3
4X(W)
Figure S8.2-3
S8.3
Therefore,
Thus,
2
1- 7rb(w - wo)ei-' dw = ei-o'
2,r -"oo
Note that the integral relating 2-7rS(w - wo) and e-o' is exactly of the form
where x(t) = e jot and X(w) = 2-7rb(co - wo). Thus, we can think of ei-o' as the
inverse Fourier transform of 2 7rb(w - wo). Therefore, 27rb(o - wo) is the Fourier
transform of eiwo'.
(c) Using the result of part (a), we have
Therefore,
Z(w) = 27rab (W
k = -00
(d)
X(w)
2 sin 3 2 sin 2
To
3 3
Figure S8.3
S8.4
Xa(f) = X(w)
W=2,wf
We know that
x(t) =- f X(w)e'j t do
Continuous-Time Fourier Transform / Solutions
S8-7
IXa(f)I
TI
2 1 1 2
TI TI TF T,
Figure S8.4
(b) Comparing
we see that
=1
Xb(v)e jvt dv,
where we have substituted v for w. Thus, the factor of 1/2-x has been distributed
among the forward and inverse transforms.
S8.5
(a) By inspection, To = 6.
2
(b) ak = 1 TO )e -k( /To)t dt
TO J T0
We integrate from -3 to 3:
1 1 1
ak= 2 (ti + 1)+ 6(t) + -6(t
2 - + c 6
J-
= 1 +l 1 + 'je (1
-Jlk6.. Cos 2k
Signals and Systems
S8-8
(c) (i) X 1(w) = x 1 (t)e -"' dt - [6(t + 1) + b(t) + -1(t - 1)]e -j't dt
- ie * + 1 + ie -"1 + cos w
= 1 + ie -j + le -i5 w
x 1 (t)
2
t
-1 0 1
x 1 (t - 6)
0 5 6 7
xl(t + 6)
2
-7
t
-6
t
-5
|tt)
0
t
x(t)
-7 -6 -5 -1 0 1 5 6 7
Figure S8.5-1
Similarly, we can periodically repeat x 2(t) to get 1(t). Thus T2 = 6. See Figure
S8.5-2.
Continuous-Time Fourier Transform / Solutions
S8-9
x 2 (t)
tilt
t
01
X2 (t -6)
2
t
t 0
t2~ x 2 (t ± 6)
6
1 7
1
2
t
-6-5 0 7
X(t)
-6 -5
IL -1 0
t2 1
t It
5 6 7
t
Figure S8.5-2
(e) Since I(t) is a periodic repetition of x 1 (t) or X 2 (0), the Fourier series coefficients
of t(t) should be expressible as scaled samples of X,(co). Evaluate X 1 ( 0 ) at w=
21rk/6. Then
6
co-= k*k2rk 1 (27rk\
6
X2 = 1 + cos - = a.
A X6
Although X,(w) # X 2(w), they are equal for w =27rk/6.
Signals and Systems
S8-10
S8.6
(a) By inspection,
7 1
e -atu(t) 9 .
a + jw
Thus,
7I 1
e- 7tu(t) 1
7 + jw
Direct inversion using the inverse Fourier transform formula is very difficult.
(b) Xb(w) = 26(w + 7) + 26(o - 7),
Xb(t) = - Xb(w)ej do = - 2
[6(w + 7) + 6(w - 7)]ei-' dw
1 1 . 2
e-
It + e = cos 7t
(c) From Example 4.8 of the text (page 191), we see that
37 2a
e alti 9 _2a
a 2 + W2
However, note that
ax(t) aX(w)
since
Thus,
1 - 1 1 7 1
-e -" l or - e -31t
2a a2 + (02 9 + .2 6
(d) Xa(w)Xb(W) = Xa(w)[26(o + 7) + 26(w - 7)]
= 2X(-7)(co + 7) + 2Xa(7)6(w - 7)
2 2
XdCO)= - (w +7) + b(o7)
7 - j7 7 +j7
1f*[2 ±7+ 2
b(co - 7) ej" do
xaOt = f- .7 6(o + 7) +
27r-- 7 -37 7 +j7
1 1 . 1 1 .
XA(t) = - e 7' + - ej_ '
ir7-j7 xr7+j7
Note that
=_ 1=( ,
1 -jr/4 1 +jir/4
7 +j7 7 2 7 -j7 7 2
Thus
we ~j3 w, 0 5 w < 1,
3
(e) Xe(w) = we -j w, -1 ! (O 0,
t0, elsewhere,
x(t) = f X(w)e-' dw = i [J' we -jamei' do - we ~i3-ejwl dw
Note that
e ax
f dx
xe" = , -(ax
ax -1)
x(t) = -
21r
1
L(
(j(t
j(t - 3)w
- 3))
(j(t - 3)w - 1)
I
0
0
egj(t -
(j(t - 3)) 2
3)w
(j )) (j(t - 3)w - 1)
0
J
which can be simplified to yield
1 [ cos (t - 3) - 1 sin (t - 3)
x(t) = 7r (t - 3)2 + (t - 3)_
Solutions to
Optional Problems
S8.7
= X(w)H(w)
S8.8
(a) Using the analysis equation, we obtain
1 rT/2_ 1
Ilk = I T/ e(t)e -jk(21r/T) t dt-T
T -T/2 T
Thus all the Fourier series coefficients are equal to 1/T.
(b) For periodic signals, the Fourier transform can be calculated from ak as
X(w) = 21 ak w-T
k=-00
Signals and Systems
S8-12
In this case,
21 2rk
PM=T 6 T
P(wj)
21Tr
T
_47 _ 2n 0 2_ 4
T T T T
Figure S8.8
2xk
=jI -- X(w)5(c
Since each summation term is nonzero only at w = 27rk/T,
From this expression we see that the Fourier series coefficients of t(t) are
ak = - X ,
Ti T
which is consistent with our previous discussions.
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
9 Fourier Transform Properties
Solutions to
Recommended Problems
S9.1
The Fourier transform of x(t) is
(a)
IX(w)I
2
Figure S9.1-1
(b)
4X(O)
1T
____________________4
21
Figure S9.1-2
S9-1
Signals and Systems
S9-2
(c)
RelX(w)t
Figure S9.1-3
(d)
Im X(o)
Figure S9.1-4
S9.2
IX(co)|
w i W
Figure S9.2-1
<X(w) = tan-' X
'\XRGA))/
=wtan-'(1), IWI <W
Fourier Transform Properties / Solutions
S9-3
4X(w)
iT
| Co
-W W
Figure S9.2-2
(b) +
1M j wIl<W
X(- C) = 0,L l otherwise
I0
0,
1 + j, otherwise
|w l < W
0, otherwise
S9.3
For x(t) to be real-valued, X(w) is conjugate symmetric:
X(-w) = X*(W)
(a) X(w) = X(w)|eiX(W)
= IX(w) Icos(4IX(w)) + j IX(w) Isin(4X(x))
Therefore,
X(-o) = |X(-w)|cos(4X(-w)) + j|X(-w)|sin(4X(-w))
= |X(w) Icos(4X(w)) - j IX(w) Isin(4X(w))
= X*(W)
S9.4
X*(w) = f-x*(t)ej-' dt
Signals and Systems
S9-4
X*(W) = x(t)ei-t dt
Therefore,
X*(-o)= f x(t)e - dt
= X(W)
Therefore,
x*(-) =0 X*(w)ej"" do
- X(co)ei" dw = -f X*(w)ejw" dw
S9.5
(a) 5{e~""'} = 5{e-"'u(t) + e"'u(-t)}
1_ 1
a +jw a- jo
2a
a2 + W2
(b) Duality states that
g(t) -G(w)
G(t) -21rg(-w)
Since
ea 2a
e a2 + W2,
Fourier Transform Properties / Solutions
S9-5
we have
1 c7
1 + t2
1 s
-
1 X
x
(C) 1 + (3t) 2 - -e since x(at)
3 |al a
(d) We are given Figure S9.5-1.
x (t)
t
-T T
Figure S9.5-1
rT A
X(w) = A e--'' dt - . (e -jwT - e )wT
-r -Jw
- 2j sin coT
=A w
sin(wT)
T
= 2TA
eT
Sketches of y(t), Y(w), and X(w) are given in Figure S9.5-2.
X(<o)
/ N 10---.--CO
T T
y(t) 4TA
2T 22T
Figure S9.5-2
Signals and Systems
S9-6
w,2T = nr, or w = n
S9.6
x(0) = 1 X( ) de
27r --
X(0) =x(t) dt
S9.7
(a) We are given the differential equation
dy(t) + 2y(t) =
x(t) (S9.7-1)
dt
4 S2 ) 4 +2
IH(W)1 2
= (4 + c2)2 +
(4 + W2)2 (4 + W2)2>
IH(w)I =
\4 + W2
(b) We are given x(t) = e~'u(t). Taking the Fourier transform, we obtain
1
X(W)= ., Hx)= 2 +jW
1+j
Hence,
( 1 1 1 1
-(1 + j)(2 + j) 1 + jo 2 + jo
S9.8
Since each of the rectangular pulses on the right has a Fourier transform given
by (2 sin w)/w, the convolution property tells us that the triangular function will
have a Fourier transform given by the square of (2 sin w)/w:
2
sin w
X(()) = 4 (0).)2
Solutions to
Optional Problems
S9.9
We can compute the function x(t) by taking the inverse Fourier transform of X(w)
(ei-
2 t~(j) e
sin oot
t
Therefore,
[sin (wot) 1
y(t) = cos (wet) s
I t
II
Fourier Transform Properties / Solutions
S9-9
I T2C Iw
0 0
cj - W() Wc Wc ± W()ic C Cj C u cj( ± (W)
Figure S9.9
S9.10
(a) x(t) = e - cos wotu(t), a > 0
= e -"'u(t)cos(wot)
Therefore,
Therefore,
= - (6
X(w)X()=1 2
2)*
7
-[S(o -
.
wo) - b(W + Wo)]
21r 9 + w2 3
j3 j3
9 + (w + 2)2 9 + (w - 2)2
= l
IwI < 27,
X2(W)
10, otherwise
Signals and Systems
S9-10
XI) X 2 (W)
-w- o---I
T 7T -21T 27T
X(W)
-3 -2r -7 F r 2r 3
Figure S9.10
S9.11
We are given the LCCDE
dy(t)
=Acos oot
dt +±2ytt)
1
H(w) = . and x(t) = A cos oot
2 +j
We have already seen that for LTI systems,
1
= \4 -+w
A cos(wot + $)
Fourier Transform Properties / Solutions
S9-11
4 0+w
9
Therefore, wo = ± 5
S9.12
d2y(t) 2dy(t)
(a) (a
51 dt2 + dt + 3y(t) = -w 2Y(W) + 2jwY(w) + 3Y(w)
= (-W2 + j2w + 3)Y(o),
A(w) = -o 2 +j2w + 3
(4dx(t) x(t) I = 4jwX(w) - X(w)
(b)~ 3dt
= (j4w - 1)X(w),
B(o) = j4w - 1,
A(w)Y(w) = B(w)X(o),
B(o)
Y(W) =
A(w) X(x)
= H(o)X(w)
Therefore,
-1 + j4w
H(o) - B(w)
A(w) -W + 3 +j2w
2
1 - j4w
2 _ 3 -j2w
S9.13
sin W t
-W W
Figure S9.13-1
X(w)
sin 0t
it
-W0 0
Figure S9.13-2
Signals and Systems
S9-12
H(w)
-2w0 2w0
Figure S9.13-3
Y(w)
2
7r
Y(w) = X(W)H(w) =
-wo ww
-o0 "0
Figure S9.13-4
S9.14
(b)
1 1
-1 1 -2 2
Figure S9.14-2
sin t sin 2t
irt rt
S9.15
Given that
yi(t) = 2rX(-w)|,
we have
2
yi(t) = 1 x(u)ejtu du
Similarly, let y 2 (t) be the output due to passing x(t) through F twice.
= (2_)2 f -0 x(u)(27r)b(t + u) du
= (21)' X(-t )
Finally, let y(t) be the output due to passing x(t) through F three times.
S9.16
We are given
Let n = 1:
x(t) = e~"'u(t), a > 0,
1
X(w) = .
a + jw
Let n = 2:
x(t) = te~"u(t),
x(t) = e~"'
1-atu(t),
(n - 1)!
1
X(w) =
(a + jw)"
We consider the case for n + 1:
x(t) = - e- at)
n!
j d 1 1
X(w) = d (a +jw)
n do [(a + j,)"n
_j d
n dw
= (-n)(a + jo)-"-Ij
n
1
(a + jo)"+I
Therefore, it is true for all n.
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
10 Discrete-Time Fourier Series
Solutions to
Recommended Problems
S10.1
The output of a discrete-time linear, time-invariant system is given by
where h[n] is the impulse response and x[n] is the input. By substitution, we have
the following.
(a) y[n] =
)k k=0
e jwn 2
1 - i- 3
k 2 ]
(b) y[n] = (1ke j,(n-
k)4
ej Cr 4
k=0 k=0
1 - fjir/4)
-o
- O
4 s "4
S1O.2
To find the period of t 1[n], we set 2 1 [n] = 1 1[n + N] and determine N. Thus
27r
1 + sin 21rn 1 + sin [1 0 ( N
= 1 + sin n + N)
Since
2s 27n
sin (0n + 2r) =sin 10n
S10-1
Signals and Systems
S10-2
12
207
1+ sin (2 n +-2
ns
= 1+ sin
20ir (i
= t
2(n +N)
2 [n + N], we have
+-2
7+ 20-)
-1+ sin 20-x n + - + -N
(12 2 12
Hence, for -irN to be an integer multiple of 21r, N must be 6.
u rn
(b) 21[n] = 1 + sin (21
2j
Therefore, N = 12.
1 1
a 2 0 = 1, a 2 -1 e -fir/2) 1 a 21 - - ear/2> = -, and
2j 2' 23 2
a2 ± 2, . . . , a2±= 0
(c) The sequence alk is periodic with period 10 and a 2k is periodic with period 12.
S10.3
The Fourier series coefficients can be expressed as the samples of the envelope
1 sin[(2Ni + 1)9/2]
ak - where N, = 1 (see Example 5.3 on
N sin (9/2) 0=2wk/N
page 302 of the text)
sin (39/2)
N sin (0/2) 0=2irk/N
(a) For N = 6,
ak =-
sin 3(2k)]
1 sin [26/)
6
1
sin (2)
1 2r i (k
sin sin -- )
26
Discrete-Time Fourier Series / Solutions
S10-3
.i [3 (27rk{
1 sin 2 / 1
ak = - ___________
12
[112
12
(c) For N = 60,
(27rk
.i [ sin. -grk
k
1
60
2n 60)] 1 s 20)
[1 (2ck 60
sin 60 sin
S10.4
X(E2)
IT IT
Figure S10.4
Hence
21 0" 2,rk
YI(Q) = X(Q) - 3 ,
ak = 5- 1 1 + cos (5
, for all k
This result can also be obtained by using the fact that the Fourier series coeffi
cients are proportional to equally spaced samples of the discrete-time Fourier
transform of one period (see Section 5.4.1 of the text, page 314).
S10.5
(a) The given relation
3
x[n] = ake k(2r'')"
k=o
S10.6
(a) ak = ak+10 for all k is true since t[n] is periodic with period 10.
(b) ak = a_, for all k is false since I[n] is not even.
(c) ake ik(21/) is real. This statement is true because it would correspond to the Four
ier series of t[n + 2], which is a purely real and even sequence.
(d) ao = 0 is true since the sum of the values of I[n] over one period is zero.
Discrete-Time Fourier Series / Solutions
S10-5
Solutions to
Optional Problems
S10.7
The Fourier series coefficients of x[n], which is periodic with period N, are given by
ak = - T x[n]e -jk(2w/N)n
NI j(wNf
n=O
ForN = 8,
ak = E x[n]e-k(/
8n=0
4
)n (S10.7-1)
ei(kr/3)
- 1e-j(kr/3)] jk(/4)n
k=O II
k=O
7
1 [e(wn/ 4)+(7r/3) 1 1 [(7n/4)--(7r/3)
2j 1 -
4
)+(-/ )]
[eJ(rn/
3
2 1 j[(rn/4)-(/3)
7 4
S[(7n/4)+(7r/3)] {( -n/ )-(7-/3)j
2j 1 ee[(rn/4 )+(v/ 3
) 1 j[(n/4)-(r/3)j
7 7
Sx[n] k jk(2/8) akjk(w/4)n
k=O k=O
+ ejrn + ej(51/ 4 )n + ej( r/ )n
3 4 7 4
1 + ej(/ 4 )n + ej( -/ )n
S10.8
The impulse response of the LTI system is
h[n] = (i)'"
Signals and Systems
S10-6
H(Q) =
n=O
(- e -j"" +
n=-o
( ~"
~) -- 1
1 1
e~"
1+ 1 -i -1
3
5 - 4 cos 0
sin ( n) = sin - (n + N)
Thus
i 3
sin 3r n + -rN)
3)
si 4r
(4
We set 37rN/4 = 21rm to get N = 8 (m = 3). Hence, the period is 8.
7
x[n] T akejk(2w/8)n
k=O
Therefore,
1
a3 = =- a5
All other coefficients a, are zero. By the convolution property, the Fourier
series representation of y[n] is given by bk, where
bk = aJH(C
T s0=(2k)/8
Thus
1 3
b 2j 5 - 4 cos(3r/4)
=b*
bk akH(Q)I
1 3 3
f or all k
4 -=-
=
5 - 4 cos[(2r/2)k] 20
Discrete-Time Fourier Series / Solutions
S10-7
ak = 1 + 2 cos (k), 0 s k - 5
bk = akH(Q)
1 + 2 cos Irk 3
(3)5 - 4 cos[(7r/3)k]
(iv) x[n] = j" + (-1)"
The period of k[n] is 4. x[n] can be rewritten as
2
x[n] = [ej(r/ )n + (e'")"
3
T k jk(2/4)n
k=O
Hence,
ao = 0, ai = 1,
a2 = 1, a3 = 0
Therefore, bo = b = 0 and
3
bi
5 - 4 cos(7r/2) E
3 3
b2 =
5 - 4 cos r 9
(b) h[n] is sketched in Figure S10.8.
h [n ]
-2 -1
012
Figure S10.8
(ii) bk = H)
1 j. rk-j 1 j .7rk
sin- - sin k =- - -sin
42 2 2 4 2 2
(iv) b0 = 0,
b 1 = H(Q) =1 - 2j,
b 2 = H(Q) =1,
b= 0
S10.9
x[n] ak
N
(c) x[n] - x [n 2 . ak(l - ei-*),
-jr N even
0, k even,
2ak, k odd
2
a* = N1 N-i1 jk( /N)n
1-N+1
=- ~Lx[n]e -jk(2/N)n ak
Therefore, ak = atk.
S10.10
(a) i'i[n] = I[n] + 9[n],
i[n + NM]= .[n + NM] + P[n + NM]
= t[n] + 9[n]
= f>[n]
Hence, ?ID[n] is periodic with period NM.
Discrete-Time Fourier Series / Solutions
S10-9
(b) Ck = Z
NMn=0
NM-I
n jk(2W/NM)n
NM-1
E [N[n] + P[n]]e jk(2T/NM)n
=
1 NM-1 NM=1
= ZNMn4N
[nje n=0
-jk(2/NM)n + 1 N
>3
n=0 :P4nje-j(/N
jk(2/NM)n
NM [n] E
n 1=0
e -jk(2./NM)(n+LN) + N
NMnO 1=0
ik(2/NM)(n+IM)
'1 1
N akIM + 1bk/N, for k a multiple of M and N,
1
akIM, for k a multiple of M,
1
1 bk/N, for k a multiple of N,
0, otherwise
S10.11
2j 2j
Therefore,
e -j(i/4)
2j 2j
All other coefficients ak are zero, in the range 0 : k - 7. The magnitude and
phase of ak are plotted in Figure S10.11-1.
lak I
0 1 2 3 4 5 6 7
4ak 3
I k
0 2 3 4 5 6 7
- 31r1 1
Figure S10.11-1
Signals and Systems
S10-10
|ak
k
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
4ak
2
3
-k
0 1 2 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
2
Figure S10.11-2
1 1
akk
0 1234567
4ak
3
3
0 k
0 1 2 4 5 6 7
Fr
3
Figure S10.11-3
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
11 Discrete-Time Fourier Transform
Solutions to
Recommended Problems
S11.1
= (i" [nje -a
= (le -j")n
n=O
1
1 -ie-
Here we have used the fact that
1
a 1-a for |a l < 1
periodically repeated.
1 -e-4
1 - e~i"
16 ( "u[n] )" 7 16
16 ,j
2
16 16)1\+2 u[n~+ 21 7 16e' 0
1- je'
S11-1
Signals and Systems
S11-2
S11.2
(a) The difference equation y[n] - ly[n - 1] = x[n], which is initially at rest, has
a system transfer function that can be obtained by taking the Fourier transform
of both sides of the equation. This yields
Y(g)(1 - le -ia) =X(Q),
so
H(Q) = Y(Q) 1
X(Q) 1 (2)-il
(b) (i) If x[n] = b[n], then X(Q) = 1 and
1
Y(Q) = H(Q)X(Q) =
so
y[n] = (I)"U[n]
Y(Q) = e _ln
12e
AMi= (1 1 ) ( ) = -2 3
so
y[n] = -2(i)"u[nj + 3(2)"u[n]
S11.3
(a) We are given a system with impulse response
1 2
Using the modulation theorem, we have
1[ 1 1
H(Q) = 2 _ - -r/2 > + 1 -
(b) We expect the system output to be a sinusoid modified in amplitude and phase.
Using the results in part (a) and the fact that
2
x[n] = le(in/ ) + le -j(in/2)
Discrete-Time Fourier Transform / Solutions
S11-3
we have
H(Q)
H()
=
('1 1 1
Q=r/2 2 1-i + 3
2 2 4jin2
so
y[n] = 2ej(1n/ 2) +
3 3
4 -ir
= - -2n
3 2
S11.4
(a) The use of the Fourier transform simplifies the analysis of the difference
equation.
1 1
y[n] + 1y[n - 1]- y[n - 2] = x[n] - x[n - 1],
1 1
Y(Q)(1 + 1e ~i" 8 -j21) = X(Q)(1 - e -jf),
Y(H) 1 - e ~iu
= H(Q) =.
X(Q) (1 + le -o)(1 - le -i")
We want to put this in a form that is easily invertible to get the impulse
response h[n]. Using a partial fraction expansion, we see that
2 -1
H() 1 + ie j + 1 - 1 e jfl)
-" I
so
h[n] = 2(-i)"u[n] -()"u[n]
(b) At R = 0, H(Q) = 0. At 9 w/4, H(Q) = 0.65ex1 22 >.Since h[n] is real, H(Q) =
H*(-Q), so H(-Q) = H*(Q) and H(-7r/4) = 0.65e -j122>. Since H(Q) is periodic
in 21,
S11.5
(a) x[n] is an aperiodic signal with extent [0, N - 1]. The periodic signal
ak = T x[n]e -jk(21/N)n
n =0
Signals and Systems
S11-4
N-1
= Z x[nje -"
n=0
Therefore,
1 (27rk
-Xi,-=ak
N N
S11.6
Solutions to
Optional Problems
S11.7
S11.8
(a) X(Q - QO) is a shift in frequency of the spectrum X(Q). We will see later that
this is the result of modulating x[n] with an exponential carrier. To derive the
modification xm[n], we use the synthesis equation:
(c)(C --
1 f f mXgeidQ
Im{X(Q)}e'"" dQ 1 f2"[X(Q) - .X*(Q) ]ejnd
= -- ei" dG
1 1(1 *
-x[n = X(Q)e -"" dQ)
2j 2j \2 2 ,
1
= -{x[n] - x*[-n])
2j
(d) Since IX() 12 = X(Q)X*(Q), we see that the inverse transform will be in the form
of a convolution. Since
I X*(G)ein dQ = 21
X(Q)e -j"" d)
27 2, 7r
=x*[ -n],
then
21
7 2,
\X(Q)| 2
ein d2 = x[n] * x*[-n]
Signals and Systems
S11-6
S11.9
We are given an LTI system with impulse response
h[n] = sin(-xn/3)
irn
(a) We know from duality that H(Q) is a pulse sequence that is periodic with period
21r. Suppose we assume this and adjust the parameters of the pulse so that
1-
21r
J H()e'""dU = h[n]
Let a be the pulse amplitude and let 2W be the pulse width. Then
a
21r
fw
-w
gd.
2-7
a (ei"w-
jn
e iflW)
a 2 sin Wn
21r n '
so a = 1 and W = 7r/3, as indicated in Figure S11.9-1.
H(Q2)
-7T _IT A T
Figure S11.9-1
periodically repeated, and that multiplication by (- 1)' shifts the periodic spec
trum by 7r, so the spectrum Y(Q) is as shown in Figure S11.9-2.
Y(72)
IT IT
4 4
Figure S11.9-2
Therefore,
S11.10
Here
dX(Q)
Y(Q) = 2X(Q) + e-i"X(Q)
x[n - 1] e -i"X(Q)
The corresponding Y(Q) is
S11.11
2[n] = akejk(21/N)n
k=(N)
(a) If we multiply both sides of this equation by e -jl(2w/N)n and sum over (N), we
obtain
t! nje -jl(2./N)n =k j(k--)(2/Nn
n= (N) k=(N) n=(N)
If k is held fixed, the summation over (N) is zero unless k = 1,which yields Nal.
Thus
a, = - ne jl(21/N)n
n=(N)
Signals and Systems
S11-8
and therefore
- ( x[n]e -jk(2r/N)n
ak =
N n=(N)
(b) We are given that x[n] is an aperiodic signal
x[n] = - X(Q)e'j" d2
(i) By multiplying both sides by e -join and summing over all n, we have
x[n]e -jIn =
2
X(Q) e"(-")n di
n = oo 1r f2 _n 0
(ii) En ej' -"")n needs to be evaluated. We can recognize that this summa
tion is a Fourier series representation
7 e'j -n =T an e j[(2(QQ-i))/T]n
exa-oon = 27 (( - Q, + 27rn)
n= -oo n =-oo
(iii) Only a single impulse in the train appears in the integration interval of
one period. So
S11.12
(a) The Fouriertransform of ejk(21/N)n can be performed by inspection using the syn
thesis formula
ake
aejk(21/N)n k( -o aN
k
Q -
2 + 2rm
k =(N) M= -0 k =(N) \
Discrete-Time Fourier Transform / Solutions
S11-9
(c) We can change the double summation to a single summation since ak is periodic:
27k 027k
2,r1( akb Q N + 27rn =27r akb Q N
- k=(N) k=-w
So we have established the Fourier transform of a periodic signal via the use of
a Fourier series:
[n] = ake(21/N)n 1 k( 2)
k=(N) k=-w
(d) We have
Therefore,
1 2xk
=k = 1 X(-)
N 10=(2rk)/N
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
12 Filtering
Solutions to
Recommended Problems
S12.1
(a) The impulse response is real because
= H(-w)e' t dw
21-o
Since H(-w) = H(w),
h(-t) = H(w)ej' dw
= h(t)
x(t) = ( aejitsl2ir
k= -00
an = M~-j[(2"k*)ITI
ak = 9odt x(t)e
Here T = 9, s0
2rk
ak = - and 5J{ei(2*kt)!TI) =
9 T)
Consequently, the Fourier transform of the filter input is as shown in Figure
S12.1-1.
X(W)
27r
9
S12-1
Signals and Systems
S12-2
Since Y(w) = H(w)X(w), the Fourier transform of the filter output is as shown
in Figure S12.1-2.
Y(W)
2ni
9
2n 0 2r
9 9
Figure S12.1-2
S12.2
S12.3
dvo
(a) RC- +v = v,
dt
Taking the Fourier transform of this equation, we have
(RCjo + 1)Ve(w) = V,(w)
We now define
H, ( Ve(w) 1
V,(w) 1 + jwRC
We can see from this expression that vc(t) is a lowpass version of v,(t).
Filtering / Solutions
S12-3
OdB
-20
-40
-60
0.1/r 1/r 1O/T 100/r
r = RC
-37r/4 L
0.1/r 1/7- 10/k 100/r
r = RC
Figure S12.3-1
V,(w) 1 + jwRC
Signals and Systems
S12-4
3
0
0
0
0
CN
1 10
R C RC
Figure S12.3-2
S12.4
Consider 0 go !5 r. In this range, the gain of the filter IH(Q) I is go. The phase shift
for the positive frequency component is +7r/2 and the shift for the negative fre
quency component is -7r/2. Since
1
x[n] = cos (gon + 0) = -[e(On+0) + e -j6on+Oj,
2
2
y[n] = 0 [eI1On+O+(/ ) + e -j[Uon+O+(ir/2)
2
=j (on+") - e ~j(fl0 n+')
[e
2
y[n] = -9 0 sin (gon + 0)
It is apparent from this expression that H(Q) is a discrete-time differentiator. A sim
ilar result holds for -7r ! go s 0.
Filtering / Solutions
S12-5
x[n] = cos (n + )
= cos - n + ,
y[n] = 7 sin - n + )
S12.5
(a) We see by examining yi[n] and y 2[n] that y[n] averages x[n] and thus tends to
suppress changes while y 2[n] tends to suppress components that have not var
ied from x[n - 1] to x[n]. Therefore, the yi[n] system is lowpass and y 2[n] is
highpass.
(b) Taking the Fourier transforms yields
1+ e-
Yi(Q) = X(Q)
2'
Hi(9) = -(1 + e -ja
2
Figure S12.5-1
Y2(Q) = X(Q) 2e
H 2(0) = -(1 - e )
2
Signals and Systems
S12-6
S12.6
(a) By inspection we see that the impulse response is given by
[ NN
h,[n] = E b[n - k]
2N + 1 k= -N
sin Q2N + 1)
(b) H 2 () = 1
2N + 1
1
Ssin(/2)
(c)
|H1 (2)l
Figure S12.6-1
IH2WA
I I
2w 2tr2
2N+1 2N+1
Figure S12.6-2
S12.7
(a) From the specification that H(O) = 1, we know that
H(w) = a
a + jw
H()| = a
(b)
(h) I1~) 1 (a 2 + Wa2)1/2)
( a = H()I
(a2 + 4)'/2 1,2
Filtering / Solutions
S12-7
The low end specification is satisfied for a > 4, as shown in Figure S12.7-1.
IH(w)I
w=2
Figure S12.7-1
IH( o)l
w =6
-f2
Figure S12.7-2
Solutions to
Optional Problems
S12.8
The easiest method for solving this problem is to recognize that passing x(t) through
H(w) is equivalent to performing
-2 dx(t)
dt
This is easily seen since
dx(t )
-2 lt -2jwX(w)
dt
Signals and Systems
S12-8
1 _ - 6
(C) XMx = =6J + ". ,
jw(6 + jw) jw 6 +jw
1
(d) X(M) =
2 + jo
x(t) = e -2 1u(t)
-2 dx(t) - 2[-2e 2
'u(t) + e -2'tt)]= 4e -2 1 u(t) - 26(t)
dt
S12.9
(a) H(Q) = H,(Q)e ~m"
(i) H,(Q) is real and even:
h4n] - H,.()
From Table 5.1 of the text (page 335), we see that the even part of h4n]
has a Fourier transform that is the real part of H,(Q). This result is easily
verified:
Sh4-nle-j = ),h[n]e(G
so
s(h,[n] + h,[-) - [H,(Q) + H*(n)],
Ev~h,[n]} Re{H,.(Q)}
Now since
Re{H,.(Q)} = H,(Q),
we have that Ev{hn]} = h,.[n], i.e., h,[n] is even, and therefore
h,[n] = h,-n]
Filtering / Solutions
S12-9
so
H,(Q)e -j" h,[n -l ,
h[n] = hn - M]
(b) h4n] = h4-n]
Since h[n] = hAn - MI,
h[n + M] = hn],
h[M - n] = h,(M - n) - M] = h,[-n],
but
hn] = h,-n] =+ h[M - n| = h[M + n]
(c) h[M + n] = h[M - n] from part (b). Since h[n] is causal, h[M - n] = 0 for
n > M. But if h[M + n] = h[M - n], then
h[M+ n] = 0 for n > M,
so
h[n] = 0 for n > 2M
Summarizing, we have
h[n] = 0 for n < 0, n > 2M
S12.10
(a)
Hi (n)
1-
i --4- Q
- I _ ( T
c 2i( r + Qc 7T
Hv( N
Figure S12.10-1
Signals and Systems
S12-10
(b) If the cutoff frequency Q. = ir/N, the total system is an identity system.
N-1 N-1
(c) h[n] = L hk[n] = ej(2 7rnk/N)ho[n]
k=O k=O
[1- e j2rn
(2nN)Ih[n],
h[n] Nho[n], n = an integer multiple of N,
0, n # an integer multiple of N,
so r[n] is as shown in Figure S12.10-2.
(d) ho[n] = -, n = 0,
N)
ho[n] = 0, n
an integer multiple of N,
=
are the necessary and sufficient conditions.
S12.11
(b)
|H(w)|
"1 -a
I I
01 W 2
4H(w)
0 OI 2
Figure S12.11-2
|H(w)I
CI i oj
Fg1 S2
Figure S12.11-3 2
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
13 Continuous-Time Modulation
Solutions to
Recommended Problems
S13.1
(a) By the shifting property,
IY(W)I
2 4w.
2w~ 3w~
Figure S13.1-1
(b) Since ej 3 wc+ji2 e r/2 e 3 wt, we are modulating the same carrier as in part (a)
except that we multiply the result by eij/ 2 . Thus
Y(w) = ej' 2X(w - 3c)
Note in Figure S13.1-2 that the magnitude of Y(w) is unaffected and that the
phase is shifted by ir/2.
|Y(w)|
I w()
4 Y(w)
7r -
2wC 4wc
Figure S13.1-2
(c) Since
ej3wct g -j3wct
cos 3 xeot = 2 +
2 2'
S13-1
Signals and Systems
S13-2
jY(w)I
-3wc 3ic
4 Y(w)
-3w 3w
2
c
Figure S13.1-3
2 and 2
Thus, the magnitude and phase of Y(w) are as given in Figure S13.1-4. Note the
scaling by 1in the magnitude.
IY(w)I
-3oic 3oc
4 Y(W)
-4c -2wc 2c 4
Figure S13.1-4
(e) Since the phase terms are different in parts (c) and (d), we cannot just add spec
tra. We need to convert cos 3wct + sin 3wet into the form A cos(3wet + 0). Note
Continuous-Time Modulation / Solutions
S13-3
that
cos(a - #) = cos a cos # + sin a sin #
4
Let a = 3wet and # = 7r/ . Then
Thus
-r
cos 3Uet + sin 3wet = \/2 cos (3wet
2e
3
wct
- ("!) + 2 e j3wt -(r/4)J
2 2
Modulating by each exponential separately and then adding yields the magni
tude and phase given in Figure S13.1-5. (Note the scaling in the magnitude.)
Y(W)|
3co~
-3co~
-3wc 3 Ct
4 Y(W)
2wc
Figure S13.1-5
S13.2
In Figure S13.2-1 we redraw the system with some auxiliary signals labeled.
Figure S13.2-1
Signals and Systems
S13-4
R 1 (w)
1
2ir
Figure S13.2-2
(a) Since m(t) = d(t) = 1, y(t) is ri(t) filtered twice by the same ideal lowpass filter
with cutoff at we. Thus, comparing the resulting Fourier transform of y(t),
shown in Figure S13.2-3, we see that y(t) = 1/(27r)x(t), which is nonzero.
Y(W)
1
2ir
WCo
Figure S13.2-3
R 1 (o - wc)
1
27
S.
Figure S13.2-4
Continuous-Time Modulation / Solutions
S13-5
R I(w + wc)
1
2n
-2wc -c
Figure S13.2-5
Since cos wot = (ewct + e -iwct)/ 2 , modulating ri(t) by cos wct yields a Fourier
transform of r 2(t) given by
Ri(w - wc) + Ri(w + wo)
2
Thus, R 2 (W) is as given in Figure S13.2-6.
R 2(CJ)
2-t
-wc II 2
co
-2Fc -- oc Wc 23c.
Figure S13.2-6
R 3(W)
-- c Wc
Figure S13.2-7
Signals and Systems
S13-6
R 4 (w) is given by shifting R 3 (W) up and down by we and dividing by 2. See Figure
S13.2-8.
R 4 (W)
1
4w
Figure S13.2-8
Y(W)
1
4w
Figure S13.2-9
sin wet = 2j
then
R2(W)
I/
-4I
4w1
Figure S13.2-10
(d) In this case, it is not necessary to know r 3 (t) exactly. Suppose r 3 (t) is nonzero,
with RA(W) given as in Figure S13.2-11.
R 4 (w)
-2wc -~oc c 2o
Figure S13.2-12
After filtering, y(t) = 0 since R 4(w) has no energy from -we, to we,.
(e) For this part, let us calculate R 2(w) explicitly.
RI(w - 2wc) + R,(w + 2wc)
R2(W) = n
R 2 (W)
1
4ff
-2wc Wc Wc 2wc
Figure S13.2-13
Signals and Systems
S13-8
R 3 (w)
47r
-coc Wc
Figure S13.2-14
R 4 (W)
c S c
Figure S13.2-15
Finally, filtering R 4 (w) gives the Fourier transform of y(t), shown in Figure
S13.2-16.
Y(o)
8ir
Figure S13.2-16
Thus,
y(t) = 1- x(t)
87r
S13.3
_|
I
-lAwi Iol
Figure S13.3-1
1
14
Figure S13.3-2
S13.4
(a) In this case,
y(t) = [A + cos wL)Mt cos(Wct + Oc)
But
Cos oMt cos(Wct + O') = 1[cos((WM ~ cjt ~ Oc)+ cOs((WM + Oc)t + cj]
Thus,
1 1
y(t) = A cos(wct + Oc) + 1cos((wM w)t O + 1cos((WM + W')t + Oc
2 2
= eiwct + e -Jwct + -e ijceJ(WM- C)t
2 2 4
t
+ 1 ejoce -i(WM-Wc)t + 1 ejOCe(WM+)dt + ! e -joce -j(wm±(cc)
4 4 4
Signals and Systems
S13-10
TO k=
Thus,
P,=2(A) + 42
A2
2
Since
max Ix(t)| 1
A A
then
1 1
9 2M 2
4,
as shown in Figure S13.4-1.
Py
1 4M
1 1
overmodulation
Figure S13.4-1
(b) The power in the sidebands is found from P, when A = 0. Thus, P,y = and the
efficiency is
E= =
2
1/(2M ) + } 2 + M2 >
Solutions to
Optional Problems
S13.5
-WC C
Cos coCt
Figure S13.5-1
z(t) rX eyt
sin coc t
Figure S13.5-2
Therefore, A(t) = x 2 (t) + y 2(t). The block diagram in Figure S13.5-3 sum
marizes how to recover A(t) from z(t).
z(t) A (t)
sin ct
Figure S13.5-3
Note that to be able to recover A(t) in this way, the Fourier transform of A(t)
must be zero for o > oc| I and A(t) > 0. Also note that we are implicitly assum
ing that A(t) is a real signal.
S13.6
From Figures P13.6-1 to P13.6-3, we can relate the Fourier transforms of all the
signals concerned.
SI(M) = - [X (W 010--
2)
X(
2j
S2(W) = [ X
x '
-O + X o
SI (w)
S2(w)
1
2
0 0
2 2
Figure S13.6-1
S3 (w) S4 (Co)
Figure S13.6-2
[
S5(W)= 1 S3 ( W- We
-) - S3 (W + W, +2
Ss (C)
1
4
-- oc c
4
Figure S13.6-3
Y(W)
2 -
~-oc
I______
oc
CO
Figure S13.6-5
S13.7
Note that
qI(t) = [si(t)cos wot + s 2 (t)sin wot]cos wot
= s 1 (t)cos 2xot + s 2(t)sin wot cos wot
H(w)
Figure S13.7
Similarly,
q 2(t) = si(t)cos wet sin wot + s 2(t)sin 2Wot
SIM sin 2w t + s 2 (t) s 2 (t) 2wot
2 2 2
Using the same filter and imposing the same restrictions on s 2 (t), we obtain y 2(t) =
s2 (t).
S13.8
(a) X(w) is given as in Figure S13.8-1.
X(w)
-1
-ojM WM
Figure S13.8-1
For Y(w), the spectrum of the scrambled signal is as shown in Figure S13.8-2.
Y(W)
Figure S13.8-2
(b) Suppose we multiply x(t) by cos Wot. Denoting z(t) = x(t)cos wmt, we find that
Z(w) is composed of scaled versions of X(w) centered at ± wM. See Figure
S13.8-3.
Z(W)
-2wm M WM 2
WM
Figure S13.8-3
Filtering z(t) with an ideal lowpass filter with a gain of 2 yields y(t), as shown
in Figure S13.8-4.
x (t) MX 2. iy(t)
_M
- IjMJ
Cos WMt
Figure S13.8-4
(c) Suppose we use the same system to recover x(t). Let y(t)cos Wmt = r(t). Then
R(w) is as given in Figure S13.8-5.
R(o)
1
2
WM WM
Figure S13.8-5
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
14 Demonstration of
Amplitude Modulation
Solutions to
Recommended Problems
S14.1
(a) We see in Figure S14.1-1 that the modulating cosine wave has a peak amplitude
of 2K = 2, so that K = 1. At the point in time when the modulating cosine wave
is zero, the total signal is A = 2, so K/A = 0.5. Therefore, the signal has 50%
modulation. See Figure S14.1-1.
S14-1
Signals and Systems
S14-2
2K
.2 .4 .6 .8
Figure S14.1-2
2K
T-0
(till 11 zm "V
irri~irur~v
UW TU
Figure S14.1-3
Demonstration of Amplitude Modulation / Solutions
S14-3
S14.2
(a) (i)
y(t)
3-
2--
1--
10 12
Figure S14.2-1
(ii)
J, (t0
3\ /1
0-- t
1 2 3 4 5 6
-2
-3 -
Figure S14.2-2
Signals and Systems
S14-4
(iii)
= 2 X(2w)
where cos irt 5 H(w), and H(w) is as shown in Figure S14.2-5. There-
fore, Y(w) is as given in Figure S14.2-6.
H(w)
7T iT
-7T 7T
Figure S14.2-5
Y(W)
37r -n _7 3n
2 2 2 2
Figure S14.2-6
P(w)
-31r I i
3
-5 7 -7r T S5.
Figure S14.2-7
Signals and Systems
S14-6
S14.3
(a) ii
(b) i
(c) iii
(d) vi
(e) v
(f) iv
(g) vii
(h) x
(i) ix
(j) viii
S14.4
(a) We are considering
N-I
X(Q) = Z x[n]e'-j,
n=O
From the convolution theorem we can compute the Fourier transform of the
product of these two sequences:
Therefore,
IX(G2)I
0
-o 0 T 007 7
Figure S14.4-1
N-I
(i) For woT = 27r(2) and N = 5, the first term is zero for
k =. . . -3, 2, 7, . . .
However, when k = 2 we have the ratio of
1 ei 21(2 / 5 -k/ 5 ) 5 0
2 1 e j21(2/5--k/5) 0
Signals and Systems
S14-8
X( 27rk)
5
2.5 --
k
0 1 2 3 4
Figure S14.4-2
Now woT = 2r-3, and the numerator and denominator are nonzero for all
k. Evaluating the preceding expression yields X(k) as shown in Figure
S14.4-3.
IX(k)|
3--
k
0 1 2 3 4
Figure S14.4-3
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
15 Discrete-Time Modulation
Solutions to
Recommended Problems
S15.1
Recall that the Fourier transform of a train of impulses p(t) is P(w), as shown in
Figure S15.1-1.
P(w)
4T 27 0 27 4in
T T T T
Figure S15.1-1
X,(w) = -2 KX(O)P(w - 0) dO
X,(W)
1
T
27r 7T
_7' 27T
T 2T 2T T
Figure S15.1-2
Thus
1 1
X,(w) = X(wo) or x, = x(t)
S15-1
Signals and Systems
S15-2
S15.2
For go = ir/2, C(Q) is given as in Figure S15.2-1.
C(n)
-27r li31
T 7lI r 1 7
2
2 2 -~ I y f
Figure S15.2-1
Thus, Y(Q) is X(Q) centered on each impulse in Figure S15.2-1 and scaled by J, as
shown in Figure S15.2-2.
Y(R)
II I I I I
3If _T int 3nt 7T 3'r 21r
2 2 4 2 4 2
Figure S15.2-2
C(92)
7T T 7IT
7_it it itf7
4 4 4 4
Figure S15.2-3
Y(u)
1
2
2
-2-r 21T
4 4 4 2 4
Figure S15.2-4
S15.3
From the lecture we know that the system in Figure S15.3-1 is equivalent to a filter
with response centered at 9 = 1r, as shown in Figure S15.3-2.
(-1)"
Figure S15.3-1
H'(92)
S1
-
4 r 14
Figure S15.3-2
Therefore, the total response is the sum of H'(Q) and H(Q), shown in Figure
S15.3-3.
I
I
a
I
II I
I
1' 3iw _i
4 4
Figure S15.3-3
Signals and Systems
S15-4
X(92)
-7r if
Figure S15.3-4
Then, after multiplication by (-1)", the resulting signal has the Fourier transform
given in Figure S15.3-5.
F eIi
Figure S15.3-5
After filtering by H(Q), the resulting signal has the spectrum given in Figure
S15.3-6.
Figure S15.3-6
-IT
A i
E
Figure S15.3-7
Discrete-Time Modulation / Solutions
S15-5
Thus, the spectrum of y[n] is given by the sum of the spectrum in Figure S15.3-8
and X(Q), as shown in Figure S15.3-8.
Y(n)
FI S
Figure S15.3-8
S15.4
(a) P(Q) is composed of impulses spaced at 21r/N, where N is the period of the
sequence. In this case N = 2. The amplitude is 2 1ra,:
1
ak 2 T p[n]e j2kn2
n=O
= 1[le -j(2wkO/2) + Oe-j(2 kl/ 2 )I 2
We now perform the periodic convolution of X(Q) with P(Q) and scale by
1/(21r) to obtain the spectrum in Figure S15.4-2.
Signals and Systems
S15-6
(b) To recover x[n] from y[n], we can filter y[n] with H(Q) given as in Figure
S15.4-3.
(c) Using p[n] we can send only every other sample of x 1 [n]. Similarly, we can send
every other sample of x2[n] and interleave them over one channel. Note, how
ever, that we can do this only because X(Q) is bandlimited to less than 7r/2.
S15.5
We note that s(t) is a periodic signal. Therefore, S(w) is composed of impulses cen
tered at (21rk)/T for integer k. The impulse at w = 0 has area given by 2rao, where
ao is the zeroth Fourier series coefficient of s(t):
ao =
T r
s(t) dt = j 1dt
A/2
-sis
=
A
T
Thus, S(w) is as shown in Figure S15.5-1.
2na_1
2 irA
T
1 2ra
2 27
TrT
Figure S15.5-1
If X(w) = 0 for I I > w/T, then R(w) will equal (A/T)X(w) in the region I co I < w/T.
Therefore, for H(co) as in Figure S15.5-2, the signal y(t) = x(t).
H(w)
T T
Figure S15.5-2
Solutions to
Optional Problems
S15.6
(a) Consider the labeling of the system in Figure S15.6.
[n] r1 $2 [n]
Figure S15.6
r[n] = $1[n]x[n]
v[n] = E r[kjh[n - k] = $[kjx[kjh[n
0 - k]
k=- k=
Letp = k - m, k = p + m. Then
= z m h[(n - m) - p]z-x[p]
p=-w0
= y[n - m]
S15.7
In general, w(t) is recoverable from w,(t) if W,(w) contains repeated versions of
W(w) that do not overlap, i.e., that have no aliasing, as shown in Figure S15.7.
W,(C)
27r c _r 27r
T T T
W(W)
Figure S15.7
Since W(c) is repeated with period 27r/T, the largest frequency component of W(W),
Co,must be less than or equal to ir/T. From the modulation property,
1
W(W) = - X(W) * X 2(W)
27r
Discrete-Time Modulation / Solutions
S15-9
Thus, since the length of a convolution of two signals is the sum of the individual
lengths,
c= W1 ± W2
S15.8
(a) If a, = -12/21r, then the portion of X(Q) around Qj will be modulated down to
about 0 = 0 and then filtered by H(Q). We now need to reshift the spectrum
back to its original position. Therefore, we need to modulate by e2'i', or # =
+Qj/27r.
(b) Consider i = 0, 1. Then the corresponding filters are as given in Figure S15.8.
-20 0 E20 2r . 2 2 + E2
0
N N N 0
Figure S15.8
S15.9
(a) Since s(t) is periodic in T, S(w) will consist of impulses located at 27rk/T. See
Figure S15.9-1.
Signals and Systems
S15-10
S(W)
_ ___ t wt
6n _4r _2n 0 27r 47r 67r
T T T T T T
Figure S15.9-2
an X
where an is the nth Fourier series coefficient of one period of s(t). For some
region Y(w) to be zero, successive terms in the sum in eq. (S15.9-1) cannot over
lap. Thus, the maximum T is such that /irT= w,, or T = 7r/we.
(c) In general, we need to find some n such that an # 0. Then we use an ideal real
bandpass filter to isolate the nth term of the sum in eq. (S15.9-1). The resulting
signal r(t) has Fourier transform R(w) given by
r(t) = x(t) + On
2r COS
I,.o (T
(remember the effect of modulating by a cosine signal). Suppose we multiply
r(t) by
I cos + O
Then
If we now use a lowpass filter with cutoff wi/T, we get x(t). If we had picked the
smallest n such that an # 0, we could have avoided the bandpass filtering
because higher harmonics are eliminated by the lowpass filter.
Discrete-Time Modulation / Solutions
S15-11
S15.10
(a) Y(Q) will consist of repeated versions of X(Q) centered at (21r/5) + 21rk and
scaled by 1. Thus, Y(Q) is as shown in Figure S15.10-1.
Y(n)
1
2
-iT 2nr_ 2 7
5 5
Figure S15.10-1
(b) Z(Q) will consist, in turn, of repeated versions of Y(Q), centered at (47r/5) +
21rk and scaled by 2, as shown in Figure S15.10-2.
Z(W)
-ote thatthe 2i 2 4
5 5 5 5
Figure S15.10-2
Note that the version of Y(Q) centered at 67r/5 contributes to the spectrum
between -3.7r/5 and ir.
(c) Two possible choices are given in Figures S15.10-3 and S15.10-4.
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
16 Sampling
Solutions to
Recommended Problems
S16.1
If wo = 7r X 10', then
cos(won X 10-3) = cos(irn) =
S16.2
The sampling function
= 67r ( (w - 61rk),
P(W)
67]
Figure S16.2-1
cos(wot) has a spectrum given by rb(o - wo) + rb(o + wo), shown in Figure
S16.2-2.
cos (CO0 t)
F r S 1. 2
Figure S16.2-2
S16-1
Signals and Systems
S16-2
37r
Figure S16.2-3
37r
6 4
-81r - 7T - 7r -2In 2
7T
4
1r 67r 87r
Figure S16.2-4
Figure S16.2-6
(b) From part (a), it is clear that (i) and (iv) are identical.
S16.3
The signal x(t) = cos(wet + 0), where wo = 2wfo, can be written as
0
x(t) = lei'eiot + le -' e -jWOt
and the spectrum of x(t) is given by
X(w) = rej'"(w - wo) + re ~j'0 (w + wo)
The spectrum of p(t) is given by
2rk)
P(W)= 2
T k -0 T,
Therefore, the spectrum of x,(t) is
X,(W) =
2Tk -
Wo+ e -job (W -
2 rk+ )
/ L e6e
and the spectrum of X,(W) is given by
X,(w) = H(w)X,(w)
wr
(a) wo = 2r X 250, 0 =-) T = 10- ,
4,
2
+ e -j"S(w -
1 X 103k + 27 X 250)]
Hence, only the k = 0 term is passed by the filter:
X,(w) = w[ej'5(w - 2w X 250) + e -i"S(w + 21 X 250)]
and
25
X,(t) = 1 ei'ei2 x ot + 1 e -'e -22rX25o
2 2
= cos (2w X 250t + 0)
2
X,(w) = >" [eb(w - 2w X 103 k - r X 750)
2 2
+ e -j'b(w - r X 103 k + x X 750)]
Only the k = 1 term has nonzero contribution:
Hence,
X,(t) = cos (2w X 250t - 0)
= cos (27 X 250t -
S16.4
1 2 * 2-7
X(w) * - T -n -- ) e
27 2An _o
= X(W)* -
w,,
12(An)fl)
- n7
X,(w)
_3ff -r 7r 3L7r
3n2n n + , 2n 3n
Figure S16.4-1
Y(W)
Fi S1.4ff
Figure S16.4-2
(b)
xP(t) X x(t)
-WM CJM
cos("f)
Figure S16.4-3
(c)
2A
yt X x(t)
-CJM liM
cos(Z)
Figure S16.4-4
(d) A is maximum when ir/A is minimum. From part (a) we see that aliasing is
avoided in X,(w) if wm zr/A. Hence, Amax = */Wu.
Signals and Systems
S16-6
S16.5
(a) The transform of the sampled function appears as in Figure S16.5-1.
Passband
-- IW
2n __n_ 0 x21T
T 2T 2T T
Figure S16.5-1
Passband
Figure S16.5-2
M MKY7 2T 2T
M
Figure S16.5-3
S16.6
Since the input x,(t) cannot be distinguished for certain values of w,the output also
should not be distinguishable for certain values of w. Hence, Q(w) must be periodic
in w. Therefore, Figure P16.6-3 is a possible candidate, but Figure P16.6-2 is not.
Sampling / Solutions
S16-7
Solutions to
Optional Problems
S16.7
P(w)
... ..
4-A _ 21T 0 2w 4r
T T T T
Xp(w)
4
NA 0
NA
NAa N A
T NA 1 W2
2n
2T
T W2
Figure S16.7-1
- = x2 Tmax = 21,
Tax W2
which is sampling at half the Nyquist rate. X,(w) for this case is given in Figure
S16.7-2.
XP(W)
_ 2 0 2T 2T
T -T 2 T
Figure S16.7-2
Hence,
A = T, wb = 27r/T, Wa = WI
Signals and Systems
S16-8
S16.8
xix(t) X2 (t)
x(t) - )(- H, (c) X, (t)
eIwot p(t)
Figure S16.8-1
(a) X(o) and Xi(o) are as shown in Figure S16.8-2. X 2(w) is as shown in Figure
S16.8-3, and X,(o) is therefore as given in Figure S16.8-4.
X(W) X1(Go)
I W- ________W,
2 Wl
W w2 (W1 - W2) (W 2 - 1 )
2 2
Figure S16.8-2
X 2 (o)
Figure S16.8-3
Xp(w)
1
T
T T'
Figure 516.8-4
Sampling / Solutions
S16-9
Hence,
27
Tmax
(W2 - wi)
(c)
xP(t) x Re g X(t)
x2TH1(o)
e -'Ot
Figure S16.8-5
S16.9
We can alias the noise region to get maximum T. This corresponds to the aliased
spectrum, shown in Figure S16.9-2.
-2W -W 0 W 2W 3W 4W 5W
Figure S16.9-2
S16.10
The spectra of x, 2(t), where T = ir/W, given in Figures S16.10-1 and S16.10-2, could
have generated x,(t):
X 1 (O)
T-
-w 6
W W
2
W
6
W
2
6w
6
Figure S16.10-1
X 2 (O)
22 2 2
Figure S16.10-2
S16.11
X,(W) = Tk X k ,
T.ax
TmT
= I , A = T, W< W < - W
W T
(ii) X(w) = 0 for coI > 2W. Hence,
T. = , A = T, 2W< W < - 2W
2W T
Sampling / Solutions
S16-11
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
17 Interpolation
Solutions to
Recommended Problems
S17.1
It is more convenient to solve this problem in the time domain than in the frequency
domain. Since x,(t) = x(t)p(t) and p(t) is an impulse train, x,(t) is a sampled ver
sion of x(t), as shown in Figure S17.1-1.
x,(t)
xP\ t
-4T 0 4T
Figure S17.1-1
= Z
n= -oO
x(nT) 5(t - nT)
Since y(t) = x,(t) * h(t) and x,(t) is impulsive, the convolution carried out in the
time domain is as shown in Figure S17.1-2.
y(t)
\t
-4T 0 5T
Figure S17.1-2
S17-1
Signals and Systems
S17-2
Since h(t) * (1/T)h(t), shown in Figure S17.1-3, is wider than the sampling
period T, the resultant w(t) is not a triangularly sampled version of x(t). w(t) con
sists of the superposition of waveforms shown in Figure S17.1-4.
1-
\ t
0 T 2T
Figure S17.1-3
w(t)
--4T i 6T
Figure S17.1-4
We note that this superposition is actually a linear interpolation between the sam
ples of x(t). For example, Figure S17.1-5 convolved with Figure S17.1-6 equals
Figure S17.1-7.
x,(t)
1 Ji1 t
0 1 2
Figure S17.1-5
Interpolation / Solutions
S17-3
t t
1 1 0 1 2 3
Now adding the shifted and scaled triangles yields Figure S17.1-8, which we see is
the linear interpolation between samples of x,(t).
0 3
Figure S17.1-8
Now since (1/T)h(t) * h(t) is Figure S17.1-9, we expect that w(t) is the linear inter
polation of x,(t) shifted right by T, shown in Figure S17.1-10.
w (t)
/000 1'0
-4T 6T
2T
S17.2
Figure S17.2-1
(b)
Y 2 (t)
- - - - -_ --
0 2 3
Figure S17.2-2
(c)
y 3 (t)
sinin(t 1)
T (t -1)
sin n (t - 2)
2 (t - 2)
Figure S17.2-3
Interpolation / Solutions
S17-5
S17.3
(a) |
P(W)
4
XP(wj)
\ A'-Z
2Ir WCm
T TT
~+ WO
Figure S17.3-1
We have
22r
27 - COM> Win,
T
2
CO,> wm,
Xp(o)
2T 2T T7- T
Figure S17.3-2
HI(o)
2T 2T
Figure S17.3-3
H2(w)
7T 7 ITT7
T 2T 2T T
Figure S17.3-4
H3(w)
T 2T 2T T
Figure S17.3-5
Interpolation / Solutions
S17-7
(c) If
and H(w) is an ideal lowpass filter whose Fourier transform is shown in Figure
S17.3-6, then
H(w)
-Coc 'oc
Figure S17.3-6
S17.4
p(t)
Figure S17.4
We want to choose H(w) so that the cascade of the two filters is an ideal lowpass
filter. In this example,
so that
S17.5
X,(Go)
1
T
-27r 0 7T 27r
Y(2)
1
T
-27r 4 0 2r
Yc()
0 X104
Figure S17.5
Solutions to
Optional Problems
S17.6
(a) We want ho[n] such that x,[n] * ho[n] = x 0[n]. We sketch ho[n] in Figure
S17.6-1.
ho [n]
0 N-I
Figure S17.6-1
Interpolation / Solutions
S17-9
N/HO(Q), |A|
H(Q) =
0, 101 -<Iir
N(ejU(N -- 1)/2
sinW/2' ~N
0, N< |
h, [n]
-N
..?TII ii 0
I it.. N
Figure S17.6-2
H () = N
N
)1
[
|HO(U)|2 =H( o
2
2
(sin 9/2
sin N/2)'2 N'
S17.7
1 - e (T+jQ)
Y(Q)H(Q) = W(Q) = 1,
H(Q) = 1 - e -(T+jQ)
Now since H(Q) = 2=_oo h[n]e -", we see by inspection that
h[O] = 1,
h[1] =-e
h[n]= 0, n # 0, 1
S17.8
(a) Since X,(w) = X(w)P(w), we conclude that X,(w) is as shown in Figure S17.8-1.
X,(w)
K s
HwH
Figure S17.8-1
{pxt)}= 6 - nw)
71= -OD
y (t)
I
WS
0
'// a
27 T
WS
Figure S17.8-2
(d) x(t) may be recovered from y(t), assuming that no time-domain aliasing
has occurred, by low-time filtering y(t) from t = - T to T and applying a gain
of oS.
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
18 Discrete-Time Processing of
Continuous-Time Signals
Solutions to
Recommended Problems
S18.1
(a) Since x,(t) = xc(t)p(t), then X,(w) is just a replication of Xe(w) centered at mul
tiples of the sampling frequency, namely 8 kHz or 27r8 X 10' rad/s. The sam
pling period is T = 1/8000.
X,(W)
T
8000
Figure S18.1-1
(b) X(Q) is just a rescaling of the frequency axis, where 21r8 X 103 becomes 2 1r.
X(Q) is shown in Figure S18.1-2.
(c) Y(Q) is the product G(Q)X(Q). Therefore, Y(Q) appears as in Figure S18.1-3.
Y(2)
8000
6000
1 T
-2ff -- 4 u 48 2f
Figure S18.1-3
S18-1
Signals and Systems
S18-2
(d) Y,(w) is a frequency-scaled version of Y(w) but only in the range 0 = --w to 7,
as shown in Figure S18.1-4. Also note the gain of T.
Ye(W)
_ - 8000 -8000
r S 4
Figure S18.1-4
S18.2
(a) The maximum nonzero frequency component of H(w) is 5 00w. Therefore, this
frequency can correspond to, at most, the maximum digital frequency before
folding, i.e., Q = x. From the relation wT = Q,we get
T. = - =2 ms
500w
(b) Since w = 500w maps to Q = 7r, the discrete-time filter G(Q) is as shown in Figure
S18.2-1.
_' .- .04
Figure S18.2-1
(c) The complete system is given by Figure S18.2-2. Note the need for an anti-alias
ing filter.
Discrete-Time Processing of Continuous-Time Signals / Solutions
S18-3
S18.3
(a) Recall that Xc(w) is as given by Figure S18.3-1.
Xc(w)
7r IT
-10000r 100007T
Figure S18.3-1
Xc(w)
co
-54000r 540007r
Figure S18.3-4
2000OXc(20000(2 - 2n))
n= 0
200007r
54T
20
n= -1
n =+ 1
iT 27r
Figure S18.3-5
X(2 )
200007 -
F41T
i4gr IS
20 20
Figure S18.3-6
Discrete-Time Processing of Continuous-Time Signals / Solutions
S18-5
YJ(w)
-14000n 14000T
Figure S18.3-7
Note aliasing since 27000 Hz is above half the sampling rate of 20000 Hz.
(c) Xe(w) is as given in Figure S18.3-8.
tT74
Xc(w)
-34000ir 340007r
Figure S18.3-8
Yc(w)
-6000i 600r
Figure S18.3-11
S18.4
It is required that we sample at a rate such that the discrete-time frequency 7r/ 2 will
correspond to c. The relation between 9, and c is Q, = ocTo. Thus, we require
= TO
W"C
or = 4(c
(ir/2)/(Ac '
we need an anti-aliasing filter that will remove power at frequencies higher than
half the sampling rate; therefore wa = 2wc. Alternatively, we note that the "folding
frequency," or the frequency at which aliasing begins, is 0 = 7r. Since 9 = r/2 cor
responds to wc, then r must correspond to 2w,.
Discrete-Time Processing of Continuous-Time Signals / Solutions
S18-7
S18.5
2
(a) We sketch X(Q) by stretching the frequency axis so that 7 corresponds to the
sampling frequency with a gain of 1/TO. We then repeat the spectrum, as shown
in Figure S18.5-1.
X(O2)
1
7T 7T
Figure S18.5-1
Y(E2)
TO
3~ 3
Figure S18.5-2
(b) We see that Y(Q) looks like X(w) filtered and then sampled. The discrete-time
frequency is -/3. Again, 2-x corresponds to 27r/To, so 7/3 corresponds to ir/3T.
Thus, if x(t) is filtered by G(w) as given in Figure S18.5-3, then y[n] = z[n].
G(c)
3 To 3To
Figure S18.5-3
Signals and Systems
S18-8
Solutions to
Optional Problems
S18.6
(a) Since we are allowing all frequencies less than 1007r through the anti-aliasing
filter, we need to sample at least twice 100r, or 2 00r. Thus, 2 007r = 2-x/T or
To = 10 ms. To find K, recall that impulse sampling introduces a gain of 1/TO.
To account for this, K must equal To, or K = 0.01.
(b) (i) Since X(w) is bandlimited to 100-r, the anti-aliasing filter has no effect. The
Fourier transform of x,(t), the modulated pulse train, is given in Figure
S18.6-1.
Xp(2)
1 - 00
To
Figure S18.6-1
X(W)
200
-2w7r 2w
2 2
Figure S18.6-2
Y(n)
200
3 3~
Figure S18.6-3
Discrete-Time Processing of Continuous-Time Signals / Solutions
S18-9
(ii) There are three effects to note in D/C conversion: (1) a gain of To, (2) a
frequency scaling by a factor of T,, and (3) the removal of repeated spec
tra. Thus, Y(w) is as shown in Figure S18.6-4.
Y(W)
IT IT
3To 3To
Figure S18.6-4
S18.7
After the initial shock, you should realize that this problem is not as difficult as it
seems. If instead of h[n] we had been given the frequency response H(Q), then
He(w) would be just a scaled version of H(Q) bandlimited to ir/T. Let us find, then,
H(Q). Using properties of the Fourier transform, we have
1
Y(Q) =- ei"YG) + X(Q)
2
Y(Q) 1
H(0) X(Q) 1 - ie-j"
Thus,
1
IH(Q)I =
\ - cos 0
( sin g
-<H() = -tan' 1 - cos
Therefore, the magnitude and phase of He(c) are as shown in Figure S18.7.
1
He(w)| = \/ -coswT' T'
0, elsewhere
-tan-'
( isinwT
I,icoswT
1 - T
0, elsewhere
Signals and Systems
S18-10
IHc(w) I
r iT
T T
4Hc(w)
/00000* I /.A
7T
T
V 7V
Figure S18.7
S18.8
T T
( 8(t-nT)
n=
Figure S18.8
From our previous study, we know that Xe(o) in the range ±ir/Tlooks just like
X(Q) in the range ±ir. Similarly, Y,(w) between -- r/T and +r/T looks like Y(Q) in
the range -ir to r. Although there is a factor of T, we can disregard it in analyzing
this system because it is accounted for in the H(w) filter. The transformation of xc(t)
to yc(t) will correspond to filtering x[n], yielding y[n]. In fact, the equivalent system
will have a system function H(Q) given by
where He(w) is the Fourier transform of h(t). Thus, we need to find Hc(w). The rela
tion between yc(t) and xc(t) is governed by the following differential equation:
d td + 4'(dyc
+y4t dyt ~ 3yc(tt)=t)
+ ) = xc(t )
Discrete-Time Processing of Continuous-Time Signals / Solutions
S18-11
S18.9
(a) It is instructive to sketch a typical y,(t), which we have done in Figure S18.9-1.
Let us suppose that T is changed by being reduced. Then the envelope of y,(t)
seems to correspond to a higher-frequency cosine. At time kt,
2wk 2w(kT) 2wt
y,(t) = cos
N
- b(t - kT) = cos NT
NT 6(t - kT) = cos NT
6(t - kT),
-
where wo = 21/NT.
If the minimum wo is wi, and since T = 21/N(o,
2T
Tmax=
A) I
Similarly,
T.in =
NW2
Signals and Systems
S18-12
(b) Recall that sampling with an impulse train repeats the spectrum with a period
of 2r/T and a gain factor of 1/T. Since 5([cos(2t/NT)J is as given by Figure
S18.9-2, Y,(w) is then given by Figure S18.9-3.
7T 7T
2,n 27r
NT NT
Figure S18.9-2
T - NT T T
7rr 2n Nr i
2| r
27r 27/N - 1 ~-WO WO= 2rN - 1\ 21r
T T N T N T
Figure S18.9-3
S18.10
(a) By sampling sc(t), we get
s[n] = setnT) = x(nT) + ax(nT - TO) = x(nT) + ax[(n - 1)T]
since T = To. Let x[n] = x(nT). Then
s[n] = x[n] + ax[n - 1]
Therefore
x[n] = -ax[n - 11 + s[n]
Discrete-Time Processing of Continuous-Time Signals / Solutions
S18-13
This is a first-order difference equation, so given s[n], we can find x[n]. Since
x(t) is appropriately bandlimited, we can then set
y[n] = -ay[n - 1] + s[n]
which will make
A
yc(t) = -x(t)
T
(b) From part (a) we see that T = A will make y(t) = x(t).
(c) Since we do not want to alias, we still need T < 7r/wM. Now
s(t) = x(t) + ax(t - TO)
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
19 Discrete-Time Sampling
Solutions to
Recommended Problems
S19.1
x[n] is given by
x[n] = (-1)" = e'"
Hence, the Fourier transform of x[n] is
S19.2
(a) x,[n] is x[n] "stretched" by interspersing with zeros, as indicated in Figure
S19.2-1.
x [n] x,[n ]
0
n 0
0
1.
Figure S19.2-1
= ( x[nje -j(21)n + 0
= = -X
=X(29)
S19-1
Signals and Systems
S19-2
= ( x[2n]e-jn
XS(92)
Figr S IT
-7r 8 -9 IT
Figure S19.2-2
Xd(Q) = X + x
X - ,r)
1 X(2)
2 2
1
2
7rI IT 7T
Figure S19.2-3
Xd(E2 )
-27r -I - 0 T 27r
Figure S19.2-4
Discrete-Time Sampling / Solutions
S19-3
S19.3
(a) For N = 1, p[n] = 1. Hence
P(Q)
214
-27 0 2-n
Figure S19.3-1
For N = 2,
Hence
P(Q) = ( O - 7rk),
k= -o
-IT 0 iT
Figure S19.3-2
ForN =L
Hence
-2rk
P(Q) = k=
(b) X,(Q), the spectrum of x,[n], is proportional to the periodic convolution of P(Q)
and X(Q). Consequently, with P(Q) as indicated in Figure S19.3-3 and X(Q) as
indicated in Figure S19.3-4, X,(Q) is shown in Figure S19.3-5. In order that x[n]
be reconstructible from x,[n] using an ideal lowpass filter, aliasing must be
avoided, which requires that
X(92)
-21 -M M 2T
Figure S19.3-4
P(92)
FiM u M 2.3-5
N
Figure S19.3-5
3 10 3)10
5
(ii) Q, = 37r/ . To avoid aliasing,
> 3w
or N<
N 5 3
Since N must be a positive integer, this requires that N = 1, i.e., x[n]
cannot be sampled.
Discrete-Time Sampling / Solutions
S19-5
S19.4
(a) The sampling period Ti is 3 ms for the system in Figure P19.4-2 to be equivalent
to the one in Figure P19.4-1.
(b) X(Q) is sketched in Figure S19.4-1.
IX(92)I
-2T2 2T
Figure S19.4-1
|Y(G)|
| l- ' e I
Figure S19.4-2
S19.5
(a) Consider xd,[n] and Xdfn], and let
Xd3[n] = Xd,[n] + aXd in]
Then
± aXdn/N], n = 0, +N,
xP3 [n] = XdXnl/N]
0, otherwise
But
x,,[n] =
0, otherwise
And
=I d[/ ,
ax,2[n]
otherwise
Signals and Systems
S19-6
Hence,
Xd[f/N] + aXd2 [n/N], n = 0, ±N,...,
x,[n] + axP2[n] =
0, otherwise
and
xP3 [n] = xp,[n] + ax 2 [n]
So system A is linear.
(b) Take xd,[n] as shown in Figure S19.5-1, with N = 4.
-1 0 1 2
Figure S19.5-1
0 4 8
Figure S19.5-2
-4 0 4
Figure S19.5-3
(c) x,[n] =
xd[n/N], n = 0, N,... ,
0, otherwise
Hence
X (E2)
7T IT
3 3
Figure S19.5-4
X(Q)
3 F-
Figure S19.5-5
S19.6
(a) x,(t) is sketched in Figure S19.6-1, and Y,(w) is sketched in Figure S19.6-2.
-- - .....x~(t)
-2To -TO To 2 To
Figure S19.6-1
Y,(w)
-H-
-2wo -- 0 CO 2coo
Figure S19.6-2
Signals and Systems
S19-8
(b) X,(w) is sketched in Figure S19.6-3, and y,(t) is sketched in Figure S19.6-4.
X,(o)
1
To
2
17 7T 7T 7r
2
TO 2TO To
Figure S19.6-3
(c) Yes, y,(t) is periodic and this is reflected in Y,(w), which contains impulses.
Solutions to
Optional Problems
S19.7
n =0, 2, +4,. . .,
(a) x,[n] =x[n,
0, n =±1, 35,...
This is sketched in Figure S19.7-1.
Ip[nn
0
Figure S19.7-1
Figure S19.7-2
and
1 -
Xa(Q) = X- x- 1 ),
which are shown in Figures S19.7-3 and S19.7-4. See Problem P19.2(b).
Xp(E2)
4 4 4 4
Figure S19.7-3
Xd(n)= Xp ()
2 2
Figure S19.7-4
S19.8
(a) We know that the Fourier transform of p[n] is given by
27rk)
k= -0
Signals and Systems
S19-10
Aliasing will be just avoided when the sampled spectra will look as shown in
Figure S19.8-1.
N ='1A
121
N1 XP(2)
=1227
1
N=2 2
-27r
1
I7
37r_ 31
i
F - I
27r
N 3
3
-21T 47r -27r - 7r 3In 2r 4_ 21T
3 3 11 0 11 3 3
Figure S19.8-2
Discrete-Time Sampling / Solutions
S19-11
H(2)
3nr 37r
Figure S19.8-3
S19.9
4
3
-4 -3 -2 -1 0 1 2 3
Figure S19.9
(b) If
z[n] = [x[nj + x[n + 11 + x[n + 2]], for all n
Signals and Systems
S19-12
and
y[n] = z[3n],
we have expressed the processing as a combination of filtering and decimati
S19.10
If h[O] = 1 and h[n] = 0 for n = kN, k # 0, it is easy to see that the samples
xo[n] that came from x[n] will be unaffected. Hence,
y[kN] = x[k], for all k
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
20 The Laplace Transform
Solutions to
Recommended Problems
S20.1
(a) The Fourier transform of the signal does not exist because of the presence of
growing exponentials. In other words, x(t) is not absolutely integrable.
(b) (i) For the case a = 1, we have that
x(t)e -' = 3e'u(t) + 4e 2 u(t)
Although the growth rate has been slowed, the Fourier transform still
does not converge.
(ii) For the case a = 2.5, we have that
x(t)e-'' = 3e-0-5 t u(t) + 4e 0 5 tu(t)
The first term has now been sufficiently weighted that it decays to 0 as t
goes to infinity. However, since the second term is still growing exponen
tially, the Fourier transform does not converge.
(iii) For the case a = 3.5, we have that
x(t)e -"' = 3e -' 5 t u(t) + 4e - 0 "'u(t)
Both terms do decay as t goes to infinity, and the Fourier transform con
verges. We note that for any value of a > 3.0, the signal x(t)e -' decays
exponentially, and the Fourier transform converges.
(c) The Laplace transform of x(t) is
3 4 7(s -7)
X( s-2 s-3 (s - 2)(s - 3)'
and its pole-zero plot and ROC are as shown in Figure S20.1.
Im
s plane
--- -F/ Re
2 17 3
7
Figure S20.1
S20-1
Signals and Systems
S20-2
S20.2
Im
s plane
-4w 0 Re
Figure S20.2-1
(b) X(s) =
The Laplace transform converges for Re~s) + a > 0, so
o + a>0,0 or o>-a,
as shown in Figure S20.2-2.
Im
s plane
Re
Figure S20.2-2
-e (s+a)t dt = e(s+a) 0
(c) X(s) = -- e atu(-t)e - dt =
JW- s + a
s + a
if Re~s) + a < 0, o + a < 0, , < -a.
The Laplace Transform / Solutions
S20-3
Figure S20.2-3
S20.3
(a) (i) Since the Fourier transform of x(t)e ~'exists, a = 1 must be in the ROC.
Therefore only one possible ROC exists, shown in Figure S20.3-1.
Im
s plane
XRe
-2 2
Figure S20.3-2
Signals and Systems
S20-4
Im
s plane
)X Re
-2 2
Figure S20.3-3
(b) Since there are no poles present, the ROC exists everywhere in the s plane.
(c) (i) a = 1 must be in the ROC. Therefore, the only possible ROC is that shown
in Figure S20.3-4.
Im
s plane
0 Re
-2 2
Figure S20.3-5
The Laplace Transform / Solutions
S20-5
(d) (i) a = 1 must be in the ROC. Therefore, the only possible ROC is as shown
in Figure S20.3-7.
s plane
- Re
Figure S20.3-7
Im
2 s plane
Re
-2
Figure S20.3-8
Signals and Systems
S20-6
s plane
Figure S20.3-9
Table S20.3
S20.4
(a) For x(t) right-sided, the ROC is to the right of the rightmost pole, as shown in
Figure S20.4-1.
Im
s plane
Re
Figure S20.4-1
The Laplace Transform / Solutions
S20-7
Im
s plane
2u 1- Re
Figure S20.4-2
Since
- 1
X(S) =
s +1 s+ 2
we conclude that
x(t) = -e - t u(-t) - (-e 2- u(-t))
(c) For the two-sided assumption, we know that x(t) will have the form
fi(t)u(- t) + fA(t)u(t)
We know the inverse Laplace transforms of the following:
1 e ~'u(t), assuming right-sided,
S+ 1 I-e'u(-t), assuming left-sided,
1 e- u(t), assuming right-sided,
s + 2 -e -2'U(-t), assuming left-sided
Which of the combinations should we choose for the two-sided case? Suppose
we choose
x(t) = e -u(t) + (-e -2')u(-t)
We ask, For what values of a does x(t)e -0' have a Fourier transform? And we
see that there are no values. That is, suppose we choose a > -1, so that the
first term has a Fourier transform. For a > -1, e - 2 te -'' is a growing exponen
tial as t goes to negative infinity, so the second term does not have a Fourier
transform. If we increase a, the first term decays faster as t goes to infinity, but
Signals and Systems
S20-8
the second term grows faster as t goes to negative infinity. Therefore, choosing
a > -1 will not yield a Fourier transform of x(t)e -'. If we choose a 5 -1, we
note that the first term will not have a Fourier transform. Therefore, we con
clude that our choice of the two-sided sequence was wrong. It corresponds to
the invalid region of convergence shown in Figure S20.4-3.
Im
s plane
Re
-2 -- -0 '/
Figure S20.4-3
Im
splane
-, Re
-2 1 0
Figure S20.4-4
S20.5
There are two ways to solve this problem.
Method 1
This method is based on recognizing that the system input is a superposition of
eigenfunctions. Specifically, the eigenfunction property follows from the convolu
tion integral
y(t) = h(r)x(t - r) dr
so that
y(t) = 2e -t/2 + 3e ~'/' for all t.
Method 2
We consider the solution of this problem as the superposition of the response
to two signals x 1 (t), x 2(t), where x 1 (t) is the noncausal part of x(t) and x 2 (t) is the
causal part of x(t). That is,
xi(t) = e -t/2U(-t) + 2e t /3 U( -),
x 2(t) = e - t 2u(t) + 2e - t 3u(t)
This allows us to use Laplace transforms, but we must be careful about the ROCs.
Now consider L{xi(t)}, where C{-} denotes the Laplace transform:
1 2 1
-Lxi(t)} = Xi(s) - _ - , Re{s) < -1
s+ s+ 3 2
The pole-zero plot and associated ROC for Yi(s) is shown in Figure S20.5-1.
Im
s plane
Re
-1__u1
Figure S20.5-1
Y2 (s) = X 2(s)H(s) = 1 + 2
(s + 2)(s + 1) (s + 1)(s+ 1)'
+ 2 3 -3
Y2(s) =
s+i s+1 s+i s+1'
t 3
y 2 (t) = -5e -'u(t) + 2e ~1 /2U(t) + 3e - / u(t)
The pole-zero plot and associated ROC for Y2(s) is shown in Figure S20.5-2.
s plane
3
2
FY
Figure S20.5-2
S20.6
(a) Since
X(s) = x(t)e -s dt
The Laplace Transform / Solutions
S20-11
X(s)
s=,+j
= J
x(t)e -"'e-" dt
We see that the Laplace transform is the Fourier transform of x(t)e -'' from the
definition of the Fourier analysis formula.
x(t) = -2 rj
X(s)e' ds
_-joo
Solutions to
Optional Problems
S20.7
1
(a) X(s) = , Refs) > -1
s + 1
Therefore, x(t) is right-sided, and specifically
x(t) = e -u(t)
1
(b) X(s) = , Re{s) < -1
+ 1s
Therefore,
x(t) = -e- t u(-t)
s
(c) X(s) = Refs) > 0
S2 + 4'
Since
. 1
eisot
S -j W
L 1
s + jw 0
1. 1(1 1=+
L{cos(wot )u(t)} = £-e*ot + - e-"=- . +
12 22 Jo s+jo
s
L{cos(wot)u(t)) 2 + W2
s+1 _ s+1 -1 2
(d) X(s) = 2=- + ,so
s 2 + 5s + 6 (s + 2)(s + 3) s + 2 s + 3
x(t) = -e - 2 tu(t) + 2e - 3 u(t)
s+1 -1 2
(e) X(s)
(s + 2)(s + 3) s + 2 s + 3'
x(t) e - 2 tu(-t) - 2e - 3 u(-t)
s - s +1
(f) X(s) = 0 < Re{s} < 1
s2(s - 1)
t
x(t) = -e u(-t) - tu(t)
S2 _ S + 1
(g) X(s) = (S + 1)2 > < Res}
S(s + 1) 2 - 3s 3s
(S + 1)2 (S + 1)2
3(s + 1) 3
(S + 1)2 (S + 1)2>
x(t) = b(t) - 3e -'u(t) + 3te-'u(t)
s + 1
(h) X(s) =
(s + 1)2 + 4
Consider
Now
S20.8
The Laplace transform of an impulse ab(t) is a.Therefore, if we expand a rational
Laplace transform by dividing the denominator into the numerator, we require a
constant term in the expansion. This will occur only if the numerator has order
greater than or equal to the order of the denominator. Therefore, a necessary con
dition on the number of zeros is that it be greater than or equal to the number of
poles.
This is only a necessary and not a sufficient condition as it is possible to con
struct a rational Laplace transform that has a numerator order greater than the
The Laplace Transform / Solutions
S20-13
denominator order and that does not yield a constant term in the expansion. For
example,
s2 + 1 1
X(s) - ,
s s
which does not have a constant term. Therefore a necessary condition is that the
number of zeros equal or exceed the number of poles.
S20.9
(a) x(t) = e -a'u(t), a < 0,
X(s) = ,
s + a'
and the ROC is shown in Figure S20.9-1.
= 1,
X(S)
s -- a
and the ROC is shown in Figure S20.9-2.
Figure S20.9-2
Signals and Systems
S20-14
s plane
Figure S20.9-3
t
(d) x(t) = e -a" I a > 0,
t
= e -a u(t) + e a'u(- t),
1 -1
X(S) = + -,
s+a s-a'
and the ROC is shown in Figure S20.9-4.
Im
s plane
-a
Figure S2 0.9-4
X(s) = fe -S dt=-,
0 s
and the ROC is shown in Figure S20.9-5.
The Laplace Transform / Solutions
S20-15
s plane
Figure S20.9-5
X(s) = ( ak t
6(t - kT)e -s dt
k=0
= Lake -skT = 1
k=0
1
a2 -2sr < 1 - 2 log a - 2sT < 0 - s > - log a
T
we have that
x(t) = cos b costwot)u(t) - sin b sin(wot)u(t)
Wo
sin(wot) 2 + W2 ,
s2 0
Signals and Systems
S20-16
we have
[s - (tan b)wo]
X(s) =cosb
s2 + W2 '
and the ROC is shown in Figure S20.9-6.
s plane
- Re
--o0
Figure S20.9-6
(i) Consider
Using linearity and the preceding sin w0t, cos wot pairs, we have
[s + (cot b)wo]
X 1 (s) = sin bW
we have
Im
-a + jo
span~e
Re
-(a + wo cot b)'
-a -jcoo
Figure S20.9-7
S20.10
=X(-s),
but X 1(s) = X(s) since x(t) = x(-t). Therefore, X(s) = X(-s).
=-X(S),
but X 1(s) = X(s) since x(t) = -x(--t). Therefore, X(s) = -X(-s).
(c) We note that if X(s) has poles, then it must be two-sided in order for x(t) =
x(-t).
Ks
(i) X(s) = ( s
(s + 1)(S - 1),
-Ks -Ks
X(-s) = = # X(s),
(-s + 1)(-s - 1) (s - 1)(s + 1)
so x(t) # x(- t).
Signals and Systems
S20-18
Im
s plane
Re
Figure S20.10-2
Parts (c)(ii) and (c)(iv) do not have any possible two-sided ROCs. Part (c)(iii) is
even, as previously shown, and therefore cannot be odd.
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
21 Continuous-Time
Second-Order Systems
Solutions to
Recommended Problems
S21.1
a + s
(b) (i) hi(t) has a pole at -a and no zeros. Furthermore, since a > 0, the pole
must be in the left half-plane. Since hi(t) is causal, the ROC must be to
the right of the rightmost pole, as given in D, Figure P21.1-4.
(ii) h 2 (t) is left-sided; hence the ROC is to the left of the leftmost pole. Since
a is positive, the pole is in the left half-plane, as shown in A, Figure
P21.1-1.
(iii) h 3 (t) is right-sided and has a pole in the right half-plane, as given in E,
Figure P21.1-5.
(iv) h 4 (t) is left-sided and has a pole in the right half-plane, as shown in C,
Figure P21.1-3.
For a signal to be stable, its ROC must include the fw axis. Thus, C, D, and F
qualify. B is an ROC that includes a pole, which is impossible; hence it corre
sponds to no signal.
S21.2
(a) By definition,
X(s) = fl x(t)e -t dt
= e -te -S dt
S21-1
Signals and Systems
S21-2
Thus
= 0(-1) 1(-1) 1
X(s) Refs) > -1
1+s 1+s 1+s
The condition on Re{s} is the ROC and basically indicates the region for which
1/(1 + s) is equal to the integral defined originally. Similarly,
H(s) = e 2tu(t)e - dt = J e -( 2 ** dt = 1
Refs) > -2
0 s + 2'
(b) By the convolution property of the Laplace transform, Y(s) = H(s)X(s) in a
manner similar to the property of the Fourier transform. Thus,
1
Refs) > -1,
(s + 1)(s + 2)
where the ROC is the intersection of individual ROCs.
(c) Here we can use partial fractions:
1 A B
(s + 1)(s + 2) s + 1 s + 2'
A = Y(s)(s + 1) = 1,
s= -1
B = Y(s)(s + 2) = -1
S= -2
Thus,
1 1
Y(s)+-,2 Refs) > -1
s+1I s+ 2
Recognizing the individual Laplace transforms, we have
y(t) = e -u(t) - e ~2 t u(t)
S21.3
(a) The property to be derived is
Since we are not integrating over s or to, we can remove the e S'o term,
Note that wherever X(s) converges, the integral defining Y(s) also converges;
thus the ROC of X(s) is the same as the ROC of Y(s).
(b) Now we study one of the most useful properties of the Laplace transform.
y(t) = x 1 (ir)x 2 (t - r) dr
Then
Substituting, we have
We can associate this last integral with Xi(s) if we are also in the ROC of xi(t).
Thus Y(s) = X 2(s)Xi(s) for s inside at least the region R, n R 2 . It could happen
that the ROC is larger,but it must contain R, n R 2
S21.4
(a) From the properties of the Laplace transform,
Y(s) = X(s)H(s)
and using the linearity property of the Laplace transform, we can take the
Laplace transform of both sides of the differential equation, yielding
s2Y(s) - sY(s) - 2Y(s) = X(s).
Therefore,
Y(s) _ 1 1
H(s) S - =
X(s) S2 - s - 2 (s - 2)(s + 1)
Signals and Systems
S21-4
Im
s plane
X Re
-1 2
Figure S21.4-1
(b) (i) For a stable system, the ROC must include thejw axis. Thus the ROC must
be as drawn in Figure S21.4-2.
Figure S21.4-2
(ii) For a causal system, the ROC must be to the right of the rightmost pole,
as shown in Figure S21.4-3.
s plane
Figure S21.4-3
Continuous-Time Second-Order Systems/Solutions
S21-5
(iii) For a system that is not causal or stable, we are left with an ROC that is
to the left of s = -1, as shown in Figure S21.4-4.
Im
s plane
iX Re
--l 2
A
Figure S21.4-4
(c) To take the inverse Laplace transform, we use the partial fraction expansion:
1 1
H(s)
H'-= 1 =
A +
B _
= +
(s +1)(s - 2) s + 1 s - 2 s+ 1 s - 2
We now take the inverse Laplace transform of each term in the partial fraction
expansion. Since the system is causal, we choose right-sided signals in both
cases. Thus,
h(t) = -ie-u(t) + le 2
'u(t)
S21.5
o = 0: Since there is a zero at s = 0, 1H(jO)I = 0. You may think that the phase is
also zero, but if we move slightly on the jw axis, 4 H(jw) becomes
IH(jl)I - - 2
/ZF 2
The phase is
(H(jw)l
4H(jo)
Figure S21.5
S21.6
The pole-zero plot is shown in Figure S21.6.
Im
X -2
s plane
Re
-5 -0.1
X -2
Figure S21.6
Because the zero at s = -5 is so far away from thejw axis, it will have virtually no
effect on IH(jw)I. Since there is a zero at o = 0 and poles near o = 2, we estimate a
valley (actually a null) at w = 0 and a peak at o i 2.
±
Continuous-Time Second-Order Systems/Solutions
S21-7
Solutions to
Optional Problems
S21.7
(a) Let y(t) be the system response to the excitation x(t). Then the differential
equation relating y(t) to x(t) is
d yt) + 2 w. dy(t) + Woy(t) = w2x(t)
dt2 dt
y(t) = -2co,
t
y(r) dr - f.
tr
f
Direct Form II
In
u .7
Figure S21.7-2
(b) (i) For a constant w,, and 0 s < 1, H(s) has a conjugate pole pair on a circle
centered at the origin of radius wn. As changes from 0 to 1, the poles
move from close to the jo axis to -CO, as shown in Figures S21.7-3,
S21.7-4, and S21.7-5.
Figure S21.7-3 shows that for s~ 0 the pole is close to thejw axis, so
|H(jw) I has a peak very near w..
s plane
f ~0 /
Re
\ AfWn
|H(ico)
I
-on II-2 2 (On d 1 -2[2
Figure S21.7-3
Continuous-Time Second-Order Systems/Solutions
S21-9
Figure S21.7-4 shows that the peaks are closer together and more spread
out at { = 0.5.
s plane
0O.5/
Re
IH(jw)
-Cn -V 1 -2 2 Wn v1 -2 2
Figure S21.7-4
Figure S21.7-5 shows that at s~ 1 the poles are so close together and far
from the jw axis that IH(jw) I has a single peak.
Signals and Systems
S21-10
(ii) For constant between 0 and 1, the poles are located on two straight
lines. As w, increases, the peak frequency increases as well as the band
width, as indicated in Figures S21.7-6 and S21.7-7.
Im
s plane
Wn2 0
Re
IH(jco)|
Figure S21.7-6
Im
s plane
wn > 0
IH(jo)|
-on V 1 -2 2 Wnon 1 -2 2
Figure S21.7-7
Continuous-Time Second-Order Systems/Solutions
S21-11
S21.8
(a) (i) The parallel implementation of H(s), shown in Figure S21.8-1, can be
drawn directly from the form for H(s) given in the problem statement.
The corresponding differential equations for each section are as follows:
d 2y 1(t) dyi(t) dx(t)
+ + y i(t)
dt2 dt dt
d'y 2 (t) 2dy 2 (t)
+ + 2y(t) =xMt)
dt dt
y(t) = yi(t) + y 2(t)
Y1 (t)
y (t)
x(t)
Figure S21.8-1
r(t)
y(t)
Figure S21.8-2
H(s) = (s2 + 2s + 1 s+ 1
(S2 + S+1 )S 2+ 2s +2)
Thus, the cascade implementation is not unique.
S21.9
(a) Decompose sin wot as
gjwot _ -jwot
2j
Then
Xi(s) = -
2js-jwO S + jAOo
1 2jwo _ CO
2j S2 + wo S2 + w'
e 2 ' sin(wot)u(t) +
(s + 2)2 + w0
and the ROC is Re{s) > -2. Here we have used our answer to part (a).
(c) Since
£ dX(s)
tx(t))
te -2'u(t)
£
-1
d [ 1 1
ds (s + 2)'
Thus
S21.10
(a) (1), (2): An impulse has a constant Fourier transform whose magnitude is unaf
fected by a time shift. Hence, the Fourier transform magnitudes of (1) and (2)
are shown in (c).
(3), (5): A decaying exponential corresponds to a lowpass filter; hence, (3)
could be (a) or (d). By comparing it with (5), we see that (5) corresponds to
kte-"'u(t), which has a double pole at -a. Thus, (5) is a steeper lowpass filter
than (3). Hence, (3) corresponds to (d) and (5) corresponds to (a).
Signals and Systems
S21-14
(4), (7): These signals are of the form e-'" cos(wot)u(t). For larger a, the
poles are farther to the left. Hence H(jw) I for larger a is less peaky. Thus, (4)
corresponds to (f) and (7) corresponds to (g).
(6): If we convolve x(t) = 1 with h(t) given in (6), we find that the output
is zero. Thus (6) corresponds to a null at w = 0, either (b) or (h). Note that (6)
can be thought of as an h(t) given by (1) minus an h(t) given by (3). Thus, the
Fourier transform is the difference between a constant and a lowpass filter.
Therefore, (6) is a highpass filter, or (b).
(b) (a), (d): These are simple lowpass filters that correspond to (i) or (ii). Since (a)
is a steeper lowpass filter, we associate (a) with (ii) and (d) with (i).
(b), (h): These require a null at zero, and thus could correspond to (iii) or
(viii). In the case of (iii), as w increases, one pole-zero pair is canceled so that
for large w,H(s) looks like a lowpass filter. Hence, (b) corresponds to (viii) and
(h) corresponds to (iii).
(c): Here we need a pole-zero plot that is an all-pass system. The only pos
sible pole-zero plot is (vi).
(e): Here we need a null on thejw axis, but not at w = 0. The only possibility
is (v).
(f), (g): These are resonant second-order systems that could correspond to
(iv) or (vii). Since poles closer to thejw axis lead to peakier Fourier transforms,
(f) must correspond to (iv) and (g) to (vii).
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
22 The z-Transform
Solutions to
Recommended Problems
S22.1
H(z) = z
z -2
Setting the numerator equal to zero to obtain the zeros, we find a zero at z = 0.
Setting the denominator equal to zero to get the poles, we find a pole at z = 1.
The pole-zero pattern is shown in Figure S22.1.
z plane
Figure S22.1
(b) Since H(z) is the eigenvalue of the input z' and the system is linear, the output
is given by
y[n] = 1 3) +3 1 (2)"
= 3()" + 4(2)"
S22.2
(a) To see if x[n] is absolutely summable, we form the sum
N-I N-i
SN = x[nhl 2' = 1 2
n=o n=O
N(N- "
(C) SN
n=o(2
S22-1
Signals and Systems
S22-2
limN-oN is finite for Ir > 2. Therefore, the Fourier transform of r-"x[n] con
verges for Ir| > 2.
__ 1
= 2z- for 12z~'I < 1
Therefore, the ROC is IzI > 2.
1
(e) X 1(e") =
1 e -s
Therefore, x 1 [n] = (2)"u[n].
S22.3
(a) Since x[n] is right-sided, the ROC is given by Iz I > a. Since the ROC cannot
include poles, for this case the ROC is given by Iz I > 2.
(b) The statement implies that the ROC includes the unit circle Iz = 1. Since the
ROC is a connected region and bounded by poles, the ROC must be
5 < IzI < 2
(c) For this situation there are three possibilities:
(i) Izi < 3
(ii) -1 < I z| < 2
S22.4
1
1 - '
X 2(z) = - l()--m'
m= 1
GO 2z
= - L (2z) m
= - 1-2z
m1
1 - z '
with an ROC of 12z| < 1, or Iz < i.
The z-Transform / Solutions
S22-3
(b) (i)
Im (z
ROC
Re lz
Figure S22.4-1
(ii)
ImZl
ROC
Re (z)
lzl= 1
Figure S22.4-2
2z
(c) (i) X 3(z) 2 z- = 2 (1 > 1, as shown
=
n=O
-1z) -Z-1 .The ROC is Izi
in Figure S22.4-3.
-1i
z/2 z
n=Al
1 -(z/2) z -2'
with an ROC of Iz/21 < 1, or IzI < 2, shown in Figure S22.4-4.
Signals and Systems
S22-4
Im
/j
A. Re Iz)
Izl=2
Figure S22.4-4
(d) For the Fourier transform to converge, the ROC of the z-transform must include
the unit circle. Therefore, for x 1 [n] and x 4[n], the corresponding Fourier trans
forms converge.
S22.5
Consider the pole-zero plot of H(z) given in Figure S22.5-1, where H(a/2) = 1.
z plane
K ' 1 to 5zeros
Figure S22.5-1
(a) When H(z) = z/(z - a), i.e., the number of zeros is 1, we have
IH(en2)I
Figure S22.5-2
|H(es")| = 1 + a2 - 2a cos 0
Hence, we see that the magnitude of H(eia) does not change as the number of
zeros increases.
(b) For one zero at z = 0, we have
H(z)= z-a'
H(e") = e
e- a
We can calculate the phase of H(e ") by [Q - 4 (denominator)]. For two zeros
at 0, the phase of H(ein) is [29 - 4 (denominator)]. Hence, the phase changes
by a linear factor with the number of zeros.
(c) The region of the z plane where IH(z) I = 1 is indicated in Figure S22.5-3.
A z plane
+- Re z =
Figure S22.5-3
Signals and Systems
S22-6
S22.6
(a) (3)"[n]
n=O
= (3z)- 1 z
n=O i e z Z
Therefore, there is a zero at z = 0 and a pole at z = 1, and the ROC is
-<1 or IzI> ,
3,
as shown in Figure S22.6-1.
z plane
WZ
/
A, i - - - -
7
37
I Figure S22.6-1
Im z plane
Re
Figure S22.6-2
S22.7
We recognize that
x[n] = (-i)"u[n]
The z-Transform / Solutions
S22-7
(c) X(z) = -a
-az
S_ a2)
(1
a 1-az
1 a2)z
1 a2
a 1-a-IZ-1
Therefore,
a-[n] -
x[n) = -
aa I u[n
Solutions to
Optional Problems
S22.8
(a) 1{r-"x[n]} = 5{r-u[n]}
= r~"e -jan
n=O
= Z (re+yu-n
n=O
< 1
Thus, Ir| > 1.
(b) (i)
Im
r =1z plane
Re
Figure S22.8-1
Signals and Systems
S22-8
(ii)
Im
r= z plane
Re
Figure S22.8-2
(iii)
r= 3 z plane
Figure S22.8-3
(c)
Figure S22.8-4
S22.9
(a) The inverse transform of
1
1 -i -
The z-Transform / Solutions
S22-9
x[n] = (I)"u[n]
3 + 2z- 1 3 + 2z-'
(b) 2 + 3z-' + z-2
(2 + z-1)(1 + z-1)
1 1
2+z- 1 + z-'
2 1
1+ iz-' 1 + z- 1
= -1)u[n] + (-1)"u[n]
S22.10
Az - a)
(a) H(z) = ( - with A a constant
(1 - azTe)
Therefore,
H A(e -j' - a)
(e' ) 1 - ae j
I A 2(e ' a)(e" - a)
=
|H~es")2 a) = A2
(1 - ae-3)(1 - aej")
and thus,
|H(e'")| = Al
(b) (i) |v 112 = 1 + a2 - 2a cos 9
1 2
(ii) |v 2 12 = 1 + -2 - - cos 9
a a
1
= - (a 2 + 1 - 2a cos Q)
1
= -| vil2
S22.11
In all the parts of this problem, draw the vectors from the poles or zeros to the unit
circle. Then estimate the frequency response from the magnitudes of these vectors,
as was done in the lecture. The following rough association can be made:
(a) (i)
(b) (ii)
(c) (iv)
(d) (iii)
(e) (iv)
Signals and Systems
S22-10
S22.12
(a) x[n] = (I)"[u[n] - u[n - 101]. Therefore,
9
X(z) = E(i)"Z--"
n=o
1 0
-1(zY
= (2z)" = 1 (2z)-
n=O 1 (2z)-1
10 11
= |z| > 0,
zl
z9(z- (1)10I">0
- 2)
and
(1)-n
u[-n - ]- z < 2
2 ~z - 2'
Summing the two z-transforms, we have
-2z
X(z) = 2 1 < |zI < 2
(z 2)D(z - 2) 2
(See Figure S22.12-2.) The Fourier transform exists.
Therefore,
X(z) = 76 cos + 4z
7 00 n j(2r/6)n+(,r/4)] j(2w/6)n+(-/4)]
Z
- l + e n
2n=0 3
7 e _ e
6__ jw/4
le(2r/6)Z- (2v/6) 1
7z e e -j4
,r
2z cosir
(2 / 6)
2 cos (2,r
z
Z ire-
,r)
(e j(21/6)
7
7z 4 3 6 4 1
z- I (2r/6>)(z - j(2/6) where |z| >
The pole-zero plot and ROC are shown in Figure S22.12-3. Clearly, the Fourier
transform exists.
The z-Transform / Solutions
S22-11
pole -zero
cancel
Figure S22.12-1
z plane
Figure S22.12-2
Figure S22.12-3
Signals and Systems
S22-12
1 - z' 0 z10 -1
(d) X(z) = z- - - 1 = z
.=, 1 - Z' z"(z - 1
The ROC is all z except z = 0, shown in Figure S22.12-4. The Fourier transform
exists.
9th-order pole
.1.0 xPC 1 -1 1
pole -zero
cancel
Figure S22.12-4
S22.13
(a) From
2-n
x[n] = s
0, n >- 0
(b) We solve this similarly to the way we solved part (a).
(Oz-')'f
log 1 12 - 1-1 Z11<1
1
z-< 1
n 2 2
01
n > 0,
x[n]
=t
0, n s 0
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
23 Mapping Continuous-Time Filters
to Discrete-Time Filters
Solutions to
Recommended Problems
S23.1
Therefore,
1
X 2(z) = 1+ 1z3
S23-1
Signals and Systems
S23-2
z plane
-3
Figure S23.1
x[n - n 0] z "OX(z),
we have
X 4(z) = z- 5X 1(z) = -Z
Delaying the sequence does not affect the ROC of the corresponding z-trans
form, so the ROC is Izi > i.
(e) Using the time-shifting property, we have
X 5(z) = z
(f) X 6(z) = =
n=0 3 1
and the ROC is I z-'| < 1, or Izi >i.
(g) Using the convolution property, we have
X7(z) =X1(z)X
6(z)= (1-1i) 1-z
and the ROC is Iz I > i, corresponding to ROC1 n ROC.
S23.2
(a) We have
y[n] - 3y[n - 1] + 2y[n - 21 = x[n]
Taking the z-transform of both sides, we obtain
Y(z)[1 - 3z-' + 2z -2 =X(z),
and the ROC is outside the outermost pole for the causal (and therefore right-
sided) system, as shown in Figure S23.2.
z plane
Figure S23.2
The system is not stable because the ROC does not include the unit circle. We
can also conclude this from the fact that
E
n=O
12n+ 1 = o
1 1
Y(z) = H(z)X(z) = _
(1 - 3z 1) (1 - 2z-1)(1 - z 1)
corresponding to different ROCs. For the ROC Iz < 1 the system impulse
response is left-sided. Therefore, since
2 -1
H(z) = 1- 2z 1 z-'
then
h[n] = - (2)2"u[-n - 11 + (1)u[-n - 1]
= -2" 1u[-n - 1] + u[-n - 1]
For the ROC 1 < Iz I < 2, which yields a two-sided impulse response, we have
h[n] = -2"*lu[-n - 1] - u[n]
since the second term corresponding to -1/(1 - z-) has the ROC 1 < Iz|.
Neither system is stable since the ROCs do not include the unit circle.
S23.3
(a) Consider
= z-"i' ( x [m]z-m
M = -_0
= z-oX(z)
It is clear that the ROC of Xi(z) is identical to that of X(z) since both require
that IU _0 x[n]z-" converge in the ROC.
(b) Property 10.5.3 corresponds to multiplication of x[n] by a real or complex expo
nential. There are three cases listed in the text, which we consider separately
here.
- : x[n] (ze -. 7 )- n
= X(ze --io),
X 2(z)= zox[n]z
x[n]
=X
Mapping Continuous-Time Filters to Discrete-Time Filters / Solutions
S23-5
Letting z' = z/zo, we see that the ROC for X 2(z) are those values of z such
that z' is in the ROC of X(z'). If the ROC of X(z) is RO < Iz| < R 1, then
the ROC of X 2(z) isRozol < |zi < RjIzo|.
(iii) This proof is the same as that for part (ii), with a = zo.
(c) We want to show that
z dX(z)
nx[n] dz
Consider
X(z) = ( x[njz-"
n = -oo
Then
dX(z)
dz)= 700n
-nx[niz-*
dz n=_-o
= -z- n 1= -o ( nx[n]z-",
so
dX(z) = n
-Z dz nx[niz-',
dz n=_-o
which is what we wanted to show. The ROC is the same as for X(z) except for
possible trouble due to the presence of the z' term.
S23.4
(a)
IH(ein)|
r = 0.9
-r= 0.75
r = 0.5
0 iT 7 21r
4 2
Figure S23.4-1
Signals and Systems
S23-6
(b)
S23.5
(a)
Im z plane
()F X )( Re
3
Figure S23.5-1
Mapping Continuous-Time Filters to Discrete-Time Filters / Solutions
S23-7
(b)
The ROC is IzI > 2. The system is not stable because the ROC does not include
the unit circle.
(c)
z plane
Figure S23.5-3
The ROC is 2 > Iz I > i, which for this case includes the unit circle. The corre
sponding impulse response is two-sided because the ROC is annular. Therefore,
the system is not causal.
(d)
Im z plane
Re
Figure S23.5-4
Signals and Systems
S23-8
The remaining ROC does not include the unit circle and is not outside the out
ermost pole. Therefore, the system is not stable and not causal.
S23.6
daytt) dytt)
(a) + 5 + 6y(t) = x(t) + 2 dxt)
dx(t)
dt2 dt dt
is the system differential equation. Taking Laplace transforms of both sides, we
have
Y(s)(s 2 + 5s + 6) = X(s)(1 + 2s),
so
Y(s) 1 + 2s
2
X(s) s + 5s + 6
1+2s 5 -3
(s + 3)(s + 2) s + 3 s + 2
Assuming the system is causal, we obtain by inspection
h,(t) = 5e - 3 u(t) - 3e -2'u(t)
(b) Using the fact that the continuous-time system function Ak/(s - s) maps to the
discrete-time system function Ak/(l - e*kz- 1) (see page 662 of the text), we
have
5 3
Hd(z) 1 e "z- - e -2TZI
5 3
Hd(Z)
1 - e -0 0 3
z-'1- e ~o 02z- 1
0 03
Letting a = e - . , b = e- 002, we have
5 _ -
33
Hd(z) =
1 - az' 1- bz-
Solutions to
Optional Problems
S23.7
(a) The differential equation is
dy(t) + 0.5y(t) x(t)
dt
Taking the Laplace transform yields
Y(s)[s + 0.5] = X(s),
H(s) = Y(S) = 1
X(s) s+ 0.5 '
H(w) = 0
jF + 0.5'
which is sketched in Figure S23.7-1.
20 logH(w)
20lgI H(0)
slope
0- -20 dB/decade
-20-. i W
.5 5 50 ...
Figure S23.7-1
Letting T = 2 yields
Y(z) 2
= Hd(z) = -
Z(z) z
Now since
|Hd(ej0 )I = |Hd(z)I I
Signals and Systems
S23-10
we have
|Hd(e j)| = 2, for all Q,
which is an all-pass filter and is sketched in Figure S23.7-2.
(c) HA(z) =
1 1
0.5- + - z
T T
(0.5T - 1) + z
The pole is located at zo = -(0.5T - 1) and, since we assume causality, we
require that the ROC be outside this pole. When the pole moves onto or outside
the unit circle, stability does not exist. The filter is unstable for
Izo l 1 or I-(0.5T - 1)1 1,
|0.5T - 11 1,
T 4
Therefore, for T > 4, the system is not stable.
S23.8
X(z-1) = Z x[nZ"
X(z 1) =
M= -0
x[-m]z- = ( x[mjz-"' = X(z)
f (z - ak)
(b) X(z) = A k
J (z - bk)
f(z - bk)
Now if a term such as (z - ak) appears in Y(z), a term such as (z* - ak)
must also appear in Y(z). For example,
Y(z) = (z - ak)(z - ak),
*(z*) = [(z* - ak)(z - ak)*]
= (z - a*)(z - ak) = Yz)
So if a pole (or zero) appears at z = ak, a pole (or zero) must also appear
at z = a* because
(z - ak) = 0 = z = ak
(e) Both conditions discussed in parts (b) and (d) hold, i.e., a real, even sequence is
considered. A pole at z = z, implies a pole at 1/z, from part (b). The poles at
z = z, and z = 1/z, imply poles at z = z* and z = (1/zp)* from part (d). There
fore, if z, = pej", poles exist at
pei"
1 =1
p
(pej")* = pe -',
( *1pei p
S23.9
-
k=-
Y x1[k] [
( x2[n
n=--oo
- k]z-"
oO
= x1[kli2(z),
k= -oo
Signals and Systems
S23-12
where
= Z{x2[n - k]}
(b) Z{x 2 [n - k) = z-kX 2(z) from the time-shifting property of the z-transform, so
XA(z) = ( x 1 [k]z-kX2(z)
k=-oo
S23.10
Consider x[n] to be composed of a causal and an anticausal part:
x[n] = x[n]u[-n - 1] + x[n]u[n]
Let
x 1 [n] = x[n]u[-n - 1],
X2[n] = x[n]u[n],
so that
x[n] = x 1 [n] + x 2[n]
and
X(z) = X 1(z) + X 2(z)
It is clear that every pole of X 2(z) is also a pole of X(z). The only way for this not
to be true is by pole cancellation from X1 (z). But pole cancellation cannot happen
because a pole ak that appears in X 2(z) yields a contribution (ak)"u[n], which cannot
be canceled by terms of x 1[n] that are of the form (bk )U[-n - 11.
From the linearity property of z-transforms, if
y[n] = y 1 [n] + y 2[n]
then
Y(z) = Y1(z) + Y 2(z),
with the ROC of Y(z) being at least the intersection of the ROC of Yi(z) and the ROC
of Y 2(z). The "at least" specification is required because of possible pole cancella
tion. In our case, pole cancellation cannot occur, so the ROC of X(z) is exactly the
intersection of the ROC of X 1(z) and the ROC of X 2(z).
Now suppose X 2(z) has a pole outside the unit circle. Since x2[n] is causal, the
ROC of X 2(z) must be outside the unit circle, which implies that the ROC of X(z)
must be outside the unit circle. This is a contradiction, however, because x[n] is
assumed to be absolutely summable, which implies that X(z) has an ROC that
includes the unit circle.
Therefore, all poles of the z-transform of x[n]u[n] must be within the unit
circle.
Mapping Continuous-Time Filters to Discrete-Time Filters / Solutions
S23-13
S23.11
(a) If h[n] = hc(nT), then
Sd[n] = E hc(kT)
k= -o
= E hd[k],
k= -o
Sd[n] = E h,(kT)
k= -00
(b) If Sd[n] = se(nT), then hd[n] does not necessarily equal he(nT). For example,
hc(t) = e ~4'U(t),
sc(t) = f e -au(r)u(t -r) d
= t - "'dr = - (1 - e -a'), t : 0
o 1a
Sd[n] = Z
k= -oo
hd[k],
so
Sd[n] - sd[n - 1] = hd[n]
S23.12
(a) From the differential equation
we have
M
Y) bksk
Ys) = He(s) = k=O
X(s) N k
( aks"
k=0
Now consider
z - z
b'
Hd(z) = Yz
X(z) N _(zk
(ak
k=0
=He,(s)
Hd(ejR)
7T T7T iT
2 -2
-1
Figure S23.12
Hd(z) = Q'
] N,(z - zpi)Ni
Therefore, if a term such as (z - zo,) appears, (z-' - zo,) must also appear. If
Hd(z) has a pole within the unit circle, it must also have a pole outside the unit
circle. If the ROC includes the unit circle, it is therefore not outside the outer
most pole (which lies outside the unit circle) and, therefore, Hd(z) does not cor
respond to a causal filter.
Consider
1
He(s) =
S23.13
(a) We are given that
A
He(s) = (S 2
From Table 9.2 of the text (page 604), we see that hc(t) = Ate'Ou(t).
To verify, consider
S O= f .e'O'utt)e -"dt,
d
1 d * sotu(t)e Stdt]
ds (s s) ds u
=SS) teso'ult)e~"dt
Therefore,
.Cr
tesotu(t) ( 2
(s -so
1 _ = a u[n]z-"
d 1) d *0
-z = C-z E a u[n]z-"
dz -- 1 z dz =
-az-2
= (-n)anu[n]z-"
(1 - az )
= natu[n]z"
(1 -z) ,= _
(d) Hc(s) = = + +
(s+ 1)(s + 2)2 s + 1 s + 2 (s + 2)2
Using the first-order pole result for 1/(s + 1) and - 1/(s + 2) and the second-
order pole result for - 1/(s + 2)2, we have
1 1 Te 2Tz
Hd(z) = 1 - e -Tz 1 - e 2 Tz- e-2rz 1)2
1 e-T- (1
After some algebra, we obtain
z[z( -e -2T + e- - Te-2) + e -4T - e -3 T+ Te-3T]
Hd(z) = (z - e-T)(z - e -2T)2
Mapping Continuous-Time Filters to Discrete-Time Filters / Solutions
S23-17
Im z plane
double pole
cl 1 C41 vRe
e-2T e-T
e-4T _ e-3 T (1 - T)
zero at z =
2
e-T e T (1+T)
Figure S23.13
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
24 Butterworth Filters
Solutions to
Recommended Problems
S24.1
(a) For N = 5 and w = ( 2 7)1 kHz, IB(jw)1 2 is given by
1
IB(jw)1 2 = 10
1 + (2000 r)
Im
poles of B(s)B(-s)
20007r
Re
x /
Figure S24.1-1
(c) For B(s) to be stable and causal, its poles must be in the left half-plane, as
shown in Figure S24.1-2.
Im
poles of B(s)
Re
-20007T
Figure S24.1-2
S24-1
Signals and Systems
S24-2
(d) Since the total number of poles must be as shown in part (b), the poles of
B(-s) must be given as in Figure S24.1-3.
Im
X
x poles of B(-s)
2000n
-Re
Figure S24.1-3
S24.2
(a) When there is no aliasing, the relation in the frequency domain between the
continuous-time filter and the discrete-time filter corresponding to impulse
invariance is
H(eJ.) = 1HH.
T \T
- Y |1 1
| 7r
2 = wT
I I
TH(e n) | |
Hai)(j)
T+0
Figure S24.2-1
Butterworth Filters / Solutions
S24-3
Since the relation between 0 and w is linear, the shape of the frequency response
is preserved.
(b) For the bilinear transformation, there is no amplitude scaling of the frequency
response; however, there is the following frequency transformation:
9=2arctan (-2)
wT
As in part (a), we can find H(jw) by reflecting H(ej") through the preceding
frequency transformation, shown in Figure S24.2-2.
Because of the nonlinear relation between Q and w, H,(jo) does not exhibit a
linear slope as H(e'") does.
(c) We redraw the transformation of part (a) for the new H(e'0 ) in Figure S24.2-3.
As in part (a), the shape of the frequency response is preserved.
Signals and Systems
S24-4
2i
3
E2 wT
TH(ei2) I
Ha(jw)
T -
I2 7
iT 2ir
3T 3T
Figure S24.2-3
We redraw the transformation of part (b) for the new H(ej") in Figure S24.2-4.
Unlike part (b), the general shape of H(eju) is preserved because of the piece
wise-constant nature of H(ej").
2nT
3 =2 arctan w
2
7T
3
I |
H(e'*)
Ha(j) |
1 I-
I - I
2tan - 2 tan T
T 3 T 3
Figure S24.2-4
Butterworth Filters / Solutions
S24-5
S24.3
(a) Using the bilinear transformation, we get
1 -
(b) Since H(s) has a pole at -a, we need a > 0 for H(s) to be stable and causal.
(c) Figure S24.3 contains a plot of (1 - a)/(1 + a), the pole location of H(z), ver
sus a.
1- a
1+ a
-1-
Figure S24.3
We see that for a > 0, (1 - a)/(1 + a) is between -- 1 and 1. Since the only pole
of H(z) occurs at z = (1 - a)/(1 + a), H(z) must be stable whenever H(s) is
stable, assuming that H(z) represents a causal h[n].
S24.4
(a) For T = 1 and the impulse invariance method, B(jw) must satisfy
B (j =_ 1)2N
B (ij) 1 + j/4
4
_ 2N= (0. 2)
1j3r/4
1+ .'W
Signals and Systems
S24-6
1 2= (0.8)2,
j2 tan (w/8)
1+
SJwc
1 2= (0.2)2
1+ j2 tan (3-/8)
S24.5
(a) The relation between 0 and w is given by Q = wT, where T = 1/15000. Thus,
1 - IH(e 0 )|
1 0.9 for 0 : - 2wr
5
0.1 - IH(e'")|1 0 for 3 -x s7
5
Note that while Hd(jw) was restricted to be between 0.1 and 0 for all o larger
than 2r( 4 500), we can specify H(eja) only up to Q = x. For values higher than
w, we rely on some anti-aliasing filter to do the attenuation for us.
(b) Assuming no aliasing,
H(e) = G j
Therefore,
27
3 1|G(jw)|1 2.7, 0 s o
15,
7w T
0.3 1
IG(jw) 1 0, s o <
5 3
(c) The relation between o and Qis given by Q = 2 arctan (w). Thus,
T
1 IG(jo)I - 0.9, 0 s s tan ,
5)
37r
0.1 - IG(jw)I 0, tan -5 o < oo
10
(d) If T changes, then the specifications for G(jw) will change for either the impulse
variance method or the bilinear transformation. However, they will change in
such a way that the resulting discrete-time filter H(e'") will not change. Thus,
He(jo) will also not change.
Butterworth Filters / Solutions
S24-7
Solutions to
Optional Problems
S24.6
(a) We first assume that a B(s) exists such that the filter specifications are met
exactly. Since
2
1
IB(jw)1
1+
we require that
1
IB(j27r)1 2 = = (10-0.05)2 = 10-0,
1
IB(j3ir)1 2
= = 10-15
Substituting N = 5.88 and co, = 7.047, we see that the preceding equations are
satisfied.
(b) Since we know that N = 6, we use the first equation to solve for we:
10- = 1 1
1 +
Solving for co,, we find that w, = 7.032. The frequency response at ( = 0.31 is
given by
_1
IB(j3)1 2
- 12 = 0.02890,
1+ 7.032
20logo0 B(j3w)| = -15.4dB
(c) If we picked N = 5, there would be no value of we that would lead to a Butter-
worth filter that would meet the filter specifications.
S24.7
We require an Hd(z) such that
0 20 logio|Hd(eO)j -0.75, 0 s Q s 0.2613w,
-20 dB 20 logio
1 Hd(e)|), 0.4018w : Q ! w
We will for the moment assume that the specifications can be met exactly. Let 9, be
the frequency where
20 logI|Hd(eja,)| = -0.75, or |Hd(es0 ")1 2 = 10-0075
where
1+ ( s
The poles are drawn in Figure S24.7.
We associate with Ha(s) the poles that are on the left half-plane, as follows:
s, = -0.9805, S2 = 0.9805e"/ 14, S3 = s*,
0
S4 = 0.9805e' "'/ S5 = S*, S6 = 0.9805e 2' 14
/" S 7 = s*
Ha(s) is given by
Ha(S) = (0.9805)7
7
f
i= 1
(s - s)
Butterworth Filters / Solutions
S24-9
S24.8
(a) Assuming no aliasing, Hd(e'u) is related to Hb(jw) by
Hdtej") = $o j , T = 2
But
H 3a
=.2r
Thus
0.27r 23a=2
ft(3 2 33a=2a
Similarly,
(.3r)
j =2b
f(e") = - 1 [j (11 r)
S24.9
(a) Using properties of the Laplace transform, we have
1
sY(s) = X(s), or H(s)
s
Signals and Systems
S24-10
(b) Here h is given by T, a is given by x[(n - 1)T], and b is given by x(nT). There
fore, the area is given by
(a + b h = T[x((n - 1)T) + x(nT)] = An
Therefore,
P[n] = y[n -1] + An.
(d) From the answer to part (a), we substitute for An, yielding
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
25 Feedback
Solutions to
Recommended Problems
S25.1
(a)
YS)
Figure S25.1-1
We have
V(s) = X(s) - Y(s)K(s) (S25.1-1)
and
Y(s) = V(s)H(s) (S25.1-2)
ForK = 0,
2 2
Q(S) = and Q(z) = z -2'
s- 2
as shown in Figures S25.1-2 and S25.1-3, respectively.
S25-1
Signals and Systems
S25-2
Im
s plane z plane
)( Re
2 / 2
IzI 1
Figure S25.1-2
Figure S25.1-3
For K = -1,
2 2
Q(s) = - and Q(z) = ,
s - 4 z -4
as shown in Figures S25.1-4 and S25.1-5, respectively.
Im
Im s plane z plane
X Re
4 Re
1 4
/ge 2
Figure S25.1-4 Figure S25.1-5
For K = 1,
Im
s plane
Re
Figure S25.1-6
2
(c) Q(s) =
s - 2(1 - K)
The pole is located at s = 2(1 -K), as shown in Figure S25.1-8.
Feedback / Solutions
S25-3
Im s plane
K-*+00 2 Ko-
Re
K=O
Figure S25.1-8
Hence, the locus of the pole is the line Re{s} = 0. Similarly, for
2
Q(z) z - 2(1 - K)'
the locus of the pole is also the line Re{z} = 0, shown in Figure S25.1-9.
Im z plane
K++o 2 K--oo
Re
K=0
Figure S25.1-9
S25.2
We use Problem P25.1.
(iii) E = H(s)
E(s)
1
R(s) G(s)
Thus,
Y(z) = H 1 (z) + Ho(z)
X(z) 1 + G(z)H 1 (z)
G2 (S):
Figure S25.2
H,(s)H2(s)
Y(s) + G1 (s)H1 (s)
1
X(s) 1 + G2(s)H 1(s)H 2(s)
1 + G1(s)H1 (s)
H1 (s)H 2(s)
1 + G1(s)H1 (s) + G2 (s)H 1 (s)H 2 (s)
S25.3
(a)
Feedback / Solutions
S25-5
G(jw) H(jw)
100 - ~
---.. 5- 40-~10 J
Figure S25.3-2
(b) From the frequency response in part (a), clearly system 1 tends to make the
response more constant and system 2 tends to resemble the inverse of G(jw).
S25.4
For the system in Figure S25.4-1, we denote the closed-loop system function by
H
V +GH
x(t) + H y(t)
Figure S25.4-1
1
(s + 1)(s + 3) 1
(a) V(s) =
1 (s + 1)(s + 3) + 1
1 + (s + 1)(s + 3)
s2 + 4s + 4 (s + 2)2
Therefore,
v(t) = te 2'u(t)
s + 3
(b) V(s)
(s + 3) + (s + 1)
1+ (s + 1)
2s + 4 2s + 2
In this case,
v(t) = le -2 t u(t)
(c) The system function G(s) = e-113 corresponds to a delay of 1, i.e., the feedback
system of Figure P25.4(a) becomes that shown in Figure S25.4-2.
Signals and Systems
S25-6
Delay
Figure S25.4-2
= O
- 1j ( t - )
2z
(d) V(z) =
1+ 1 -(i' zz
z1 z
Z-1
(1-i ~)+ (!z~- Iz- 2)
z1
1+ z-1 -z-2
-1
(1 -'-)(1 + i2z-I)
6 6
5 5
1
1-z- 1+ z-
Therefore,
v[n] = [)"u[n] -- (-)u[n]]
H(z) 2 - iz-1
(e) V(z) =
1 + H(z)G(z)
1 Z-) z-1
+2
1 +36 z-'
1 -iz
(1 - iz-1) + (0 - iz-1)-1
2 - 2z- + nz-2
1-5z 1l~2
1+ 1z-1 - Iz-2
Thus,
v[n] = 3ib[n + 1] - 23[n] + 14[n - 1],
where f[n] is v[n] in part (d).
Feedback / Solutions
S25-7
Solutions to
Optional Problems
S25.5
y(t) = K 2w(t) + KK 2v(t) (S25.5-1)
By taking the transform of eq. (S25.5-1), we have
Y(s) = K 2W(s) + KK 2Vs)
Also
Therefore,
and
Y K2W(s) + K 1K 2X(s)
Y
=
- ) K 1K 2s
s + a
(s + a)[K 2W(s) + K1K 2X(s)]
(1 - K 1K2)s + a
S25.6
(a) The system function of the system given in Figure P25.6 must be determined
first. So we write down the difference equation
y[n] = x[n] + y[n - 1] + 4y[n - 2]
z2 z - 4=0, or z 2
2 - 2
Since Iz I > 1 for at least one pole the system is unstable.
(b) With closed-loop feedback, the difference equation is
y[n] = x,[n] - Ky[n - 1] + y[n - 1] + 4y[n - 2]
Thus,
z2
H(z) = z2 + (K - 1)z - 4
Signals and Systems
S25-8
2 2 - 2
K+ 1 \/(K - 1)2 + 16
Thus,
K 2 + 2K + 1 = K 2 - 2K + 17,
4K = 16, or K= 4
We can also calculate z 2:
Z2= -4
Similarly, zi = -1, z 2 = 4 for K = -2. Observe the root locus in Figure
S25.6-1.
Im Im
K> 0 z plane K <O z plane
K =-2
.00.9 K=4
Re Re
-1.56 2.56 -1.56 \ /12.56
Figure S25.6-1
Observe that if one of the poles is inside zI 1, the other is outside. Hence,
the system is unstable for all values of K.
(c) The difference equation can be written as
y[n] = x,[n] + y[n - 1] + (4 - K)y[n - 2]
Therefore,
z2
H(z) =
z- z + (K - 4)
In this case, the poles are located at
1 +/17 -4K
2 2
For a stable system, we want
izi < 1,
1 1/17 - 4K
|z| = ->+ 2
If we set 17 - 4K > 0, then
+ 17 - 4K <1 1
-<
Feedback / Solutions
S25-9
/17 -4K 1
2 2'
17 - 4K < 1,
K> 4
Now suppose 17 - 4K < 0. Then
K>O K=5
-1.56 2.56
Figure S25.6-2
S25.7
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
26 Feedback Example:
The Inverted Pendulum
Solutions to
Recommended Problems
S26.1
Ld 2 0(t)
(a) Ldz6(t) = g0(t) a(t) + Lx(t),
Ld2 (t)
dt2 - gA(t) = Lx(t)
Taking the Laplace transform of both sides yields
szLO(s) - gO(s) = LX(s),
O(s) = 2X(s)
s2 - g|L
0(s)
X(s) s- g/L (s + \g/L)(s - \ )
The pole at \g/L is in the right half-plane and therefore the system is unstable.
(b) We are given that a(t) = K0(t). See Figure S26.1-1.
Figure S26.1-1
0(s) H
X(s) 1 + GH'
so, with
and G=
H = s g/L
(s)/X(s) is given by
0(s) 1
X(s) s' - (g/L) + (K/L)
s=+K-g L '
which implies that the system is unstable. Any K < g will cause the system
poles to be pure imaginary, thereby causing an oscillatory impulse response.
S26-1
Signals and Systems
S26-2
1
H(s) =1
2 g K1 K2
L L L
1
2K2
s+- s ±
K1 -- g
L L
The poles are at
S26.2
(a) Here
2
H(s) =2 +2
2+ os + W2
G(s) = K
Feedback Example: The Inverted Pendulum / Solutions
S26-3
2
S +2 (n Wn
,, Z~
where w, = Wn(1 + K)"/2
± 2/(i) +
(s ( W,
here - = {
s2 + F,+ )2 w
(A.n
Therefore,
= wn(l + K)" 2 ,
2
(1 + K) ,
2
A = =
An 1 +K'
He(s) =
s2 + 2(,,s + (o2(1 + K)
(c)
K = +oo Im
s plane
K=_xo=
1 - ' Re
K = +oo
Figure S26.2
The poles start out at ± 0o, approach each other and touch at K = 2 - 1, and
then proceed to -fo, ±joo.
Signals and Systems
S26-4
S26.3
s= - r + KK 2a)
r + K1 K a
2 0
>0
If /3> 0, then r/a > -KK 2; if # < 0, then r/a < -K 1 K 2.
S26.4
K(s + 100)
H(s) = K
+ 1) s + 100 + Ks + K
1 + K(s
S + 100
K(s + 100)
(K + 1) s +100 + K)
K + 1]
(a) K = 0.01,
H(s) = 0.01(s + 100)
1.01(s + 99.0198)
The zero is at s = -100, and the pole is at s = -99.0198.
(b) K= 1,
S + 100
H(s) =
2 s+ 2
The zero is at s = -100; the pole is at s = -50.5.
(c) K =10,
H(s) = 10(s + 100)
11 s+ 110)
110
The zero is at s = -100; the pole is at s = -10.
(d) K = 100,
100(s + 100)
H(s) =
0200)
101 S
The zero is at s = -100; the pole is at s = -1.9802.
Feedback Example: The Inverted Pendulum / Solutions
S26-5
S26.5
s-t+ 1 1
(a) H(s) = K +=1
+ K s+ 1+ K
1++1
The pole is at s = -1 - K, as shown in Figure S26.5-1.
Im Im
K>O s plane K<O s plane
- . Re Re
-1 -1
Figure S26.5-1
The pole moves from infinity to negative infinity as K changes from negative
infinity to infinity.
1
s-1i s+3
(b) H(s) =
1 (s + 3)(s - 1) + K
+1 (K
+ s
s + 3
2
s + 2s + K - 3
The poles are at s, = -1 ± \1 - (K - 3), as shown in Figure S26.5-2.
Im Im
K>O KK 0
s plane s plane
K =4
Re
-3 -1 +1 3
Figure S26.5-2
The poles start at ± oo when K = -oo, move toward -1, touch when K = 4,
and proceed to -1 ± joo as K approaches positive infinity.
Signals and Systems
S26-6
Solutions to
Optional Problems
S26.6
(a) The poles for the closed-loop system are determined by the denominator of the
closed-loop transfer function
Kz
1 + = 0
(z - 1)(Z - i')'
so
(z - -)(z - 1) + Kz = 0
Since we are told a pole occurs when z = -1, we want to solve the equation
for K:
K (z - 2)(z-) _ 15
z z=-1 8
(b) In a similar manner to that in part (a),
K (z - 2)(z -) -3
z 1z=1 8
(c) From the root locus diagram in Figure P26.6, we see that for K > 0 when K
exceeds a critical value of K = -', as determined in part (a), one root remains
outside the unit circle. Similarly, when K < -8, one root is outside the unit
circle. Therefore, to ensure stability, we need
- < K < '
S26.7
(a) The closed-loop transfer function is
Y(s) _ H(s) _ He(s)H,(s)
X(s) 1 + G(s)H(s) 1 + He(s)H,(s)
and, therefore, from the given He(s) and H,(s), we have
Ka
Y(s) S+ a Ka Ka
X(s) + Ka s + a + Ka s + (K + 1)a
s + a
The system is stable for denominator roots in the left half of the s plane; there
fore -(K + 1)a < 0 implies that the system is stable.
Now since E(s)He(s)H,(s) = Y(s), we have
E(s) 1 s + a s + a
X(s) 1 + He(s)H,(s) s + a + Ka s + (K+ 1)a
The final value theorem, lim , e(t) = lim,_. sE(s), shows that
()Y(S_ He(s)H,(s)
X(s) 1 + He(s)H,(s)
K1+ -2 a
s )s + a
1 + K, + s2 a
SS + K2 Kia
(sK + K 2)a \ K 1
s(s + a) + (Kis + K 2)a S2 + sa(K + 1) + K 2a
The poles for this system occur at
- a(K + 1) (a(K+ 1) - K 2a
Note that if a(K + 1) > 0 and if K 2a > 0, we are assured that both poles are
in the left half-plane. Therefore, a(K + 1) > 0 and K 2a > 0 are the conditions
for stability. Now since
1
E(s) = X(s) 1 + He(s)H,(s)
1 s(s + a)
S s2 + a(K1 + 1)s + K 2a'
then
lim sE(s) = 0 implies that lim e(t) = 0,
s-0 t-CO
for a(KI + 1) > 0 and K 2a > 0, so we can track a step with this stable system.
S26.8
(a) = H(s)C(s)
X(s)
(s + 1)(s - 2)(s +
We can see from this expression that the overall transfer function for the sys
tem is
Y(s) 1
X(s) (s + 1)(s + 3)'
a stable system. In effect, the system was made stable by canceling a pole of
H(s) with a zero of C(s). In practice, if this is not done exactly, i.e., if any com
Signals and Systems
S26-8
ponent tolerances cause the zero to be slightly off from s = 2, the resultant
system will still be unstable.
2-v - (K - 2)
We see from this that at least one pole is in the right half-plane, i.e., there is
instability for all values of K.
1
(c)()Y(S) = 1+K(s + a) (s +
X() ____+_a)_(s+_1)(s_-_2 1 -
1)(s 2)
X(s) 1____
l+~sa(±+ )(- 2 )
K(s + a)
(s + 1)(s - 2) + K(s + a)
K(s + a) K(s + a)
s2 - s - 2 + Ks + Ka s 2 + (K - 1)s + (Ka - 2)
The poles are at
(K - 1) K_
- 2 -V -2- - (Ka - 2)
Now, if Ka - 2 > 0, the system is stable. K > 2/a because a > 0 is assumed.
This is true for 1 > a > 0 and 2 > a > 1. For a > 2, the system is stable for
K > 1.
Y(S) K(s + a)
(d) X(s) 2
s 2 + (K - 1)s + (Ka - 2)'
We want K - 1 = w,, 2K - 2 = W. So
(K - 1)2 = 2K - 2,
K = 3 or K= 1
IfK= 1, thenw, = 0, so K= 3 implies thatw, = 2.
S26.9
E(s) 1 S1
(a) -_____-where
X(s) 1 + H(s) s' + G(s)'
K f (s - #K
k=1
G(s) n-1
17 (s - aK)
k=1
S-I
(b) E(s) = for 1 = 1, x(t) = u- 2(t)
s + G(s)
So
1 1
lim sE(s) = Constant
S-0 s + G(s) 1,=o g
5 1-k 2-k
5
(c) E(s) = s + G(s), sE(s) = s +
G(s)
For k > 2,
lim sE(s) = co, lim e(t) = oo
s-o
S 1-k
(d) (i) E(s) = ,G sE(s)
S' + G(s) s' + G(s)
If k :5 1, then
1-k+1
lim sE(s) = lim 0 0,= 0,
s-O s-o s' + G(s) 0+ g
so lim,_., e(t) = 0.
(ii) If k = 1 + 1 and since
1-k
E(s) = s' +
G(s)
then
1 1
lim sE(s) = lim = Constant
S-0 S-O S' + G(s) g
Thus, lim,me(t) = Constant.
(iii) If k > 1 + 1, then since
S1-k gl-k+1
E(s) = S' + G(s) sE(s)
s' + G(s)
lims- 0 sE(s) = oc implies lim_. e(t) = 00.
S26.10
E(z) 1
(a) X(z) 1 + H(z)
z
X(z) z- 1 z(z + i)
E(z) = 1 + H(z) 1 (z - 1)(z + i) + 1
1+(z - 1)(z + 1)
z 2 + .z - 2
1±2+
z2 2Z + -2 Z'2
- Z + 2
The poles are at i + - . These poles are inside the unit circle and therefore
yield stable inverse z-transforms, so e[n] = b[n] + (2 stable sequences). So
lim-,.e[n] = 0.
Signals and Systems
S26-10
)(z - 1)B(z)
1I
-z(z - 1)B(z) for x[n] = u[n]
+ A(z)
zB(z)
(z - 1)B(z) + A(z)
Furthermore, we know that
Y(z) H(z) (z - 1)B(z)
X(z) 1 + H(z) (z - 1)B(z) + A(z)
There are no poles for Iz I > 1 because h[n] is stable. Therefore,
zB(z)
E(z) =
(z - 1)B(z) + A(z)
has no poles for IzI > 1, and lim,-e[n] = 0.
(c) H(z) =
- z 1'
E(z) 1 z - 1
X(z) 1 + H(z) z
z-1 fz-1\
E(z) = z1 X(z)
zz (z for x[n] = u[n]
= 1 =* e[n] = b[n],
so e[n] = 0, n -1
2z 1 + 1z- 2
(d) H(z)
(1+ -Iz ')(1 - z-1)'
E(z)
()X(z)
= 1H H(z) X(z 1
1 + H(z) ' E(z)
For x[n] = u[n], we have
1
X(z) =1 -1
Feedback Example: The Inverted Pendulum / Solutions
S26-11
We would like
N-1
e[n] = a[n - k],
k=O
N-1
E(z) = akz-'
k=O
Therefore,
N-1
1 z1)
Z1- T akz- k)
(k=0
H(z) = N-
(1-- 1) T akz -k)
(k=0
z-I + z- 2 - z- E(z) 1
(f) H(z) =
(1 + z-1)(1 - z)2 X(z) 1 + H(z)
Now x[n] = (n + 1)u[n] and
X(z) = (1 - z-1)2>
1
(1 + z-1)(1 - z-1)2 (1 - z-1)2
E(z) = 2 1 2
(1 + Z- 1)(1 - z-1) + z- + z- -z
1 + z-'
and
e[n] = b[n] + b[n - 1]
=0, n >- 2
MIT OpenCourseWare
http://ocw.mit.edu
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.