$$
β Back to LINEAR SYSTEMS AND SIGNALS Overview
(l)
- 5.1-5 Showing all work, evaluate $β n=0 n(β3/2)βn.
- 5.1-6 Using only the z-transforms of Table 5.1, determine the z-transform of each of the following signals.
(a) u[n] βu[nβ2] (b) Ξ³ nβ2u[nβ 2] (c) 2n+1u[nβ 1] +enβ1u[n] (d) 2βn cosΟ 3 n u[nβ1] (e) nΞ³ nu[nβ 1] (f) n(nβ1)(nβ2)2nβ3u[nβm] for m = 0, 1, 2, 3 (g) (β1)nnu[n]
(h)
5.1-7 Find the inverse unilateral z-transform of each of the following: (a) z(z β4)
(a)
\n(b)
\n(c)
\n(d)
\n(e)
\n(f)
\n(g)
(h)
\n(i)
\n(j)
\n(k)
\n(k)
5.1-8 (a) Expanding X[z] as a power series in zβ1, find the first three terms of x[n] if
(b) Extend the procedure used in part (a) to find the first four terms of x[n] if
- 5.1-9 A right-sided signal x[n] has z-transform given by X[z] = z6+2z5+3z4+4z3 *z*4β1 . Using a power series expansion of X[z], determine x[n] over β5 β€ n β€ 5.
- 5.1-10 Find x[n] by expanding
as a power series in zβ1.
- 5.1-11 (a) In Table 5.1, if the numerator and the denominator powers of X[z] are M and N, respectively, explain why in some cases N β M = 0, while in others N β M = 1 or N β M = m (m any positive integer).
- (b) Without actually finding the z-transform, state what is N β M for X[z] corresponding to x[n] = Ξ³ nu[nβ4].
- 5.2-1 For a discrete-time signal shown in Fig. P5.2-1, show that
Find your answer by using the definition in Eq. (5.1) and by using Table 5.1 and an appropriate property of the z-transform.
Figure P5.2-1
- 5.2-2 Determine the unilateral z-transform of signal x[n] = (1βn) cos Ο 2 (nβ1) u[nβ1].
- 5.2-3 Suppose a DT signal x[n] = 2(u[nβ10] βu[nβ 6]) has a transform X(z). Define Y(z) = 1 2zβ3 d dzX(2z). Using graphic plot or vector notation, determine the corresponding signal y[n].
- 5.2-4 Suppose a DT signal x[n] = 3(u[n] βu[nβ5]) has a transform X(z). Define Y(z) = 2zβ4 d dzX z 2 . Using graphic plot or vector notation, determine the corresponding signal y[n].
- 5.2-5 Find the z-transform of the signal illustrated in Fig. P5.2-5. Solve this problem in two ways, as in Exs. 5.2d and 5.4. Verify that the two answers are equivalent.
- 5.2-6 Using z-transform techniques and properties (no time-domain convolution sum!), determine the convolution y[n] =( 1 2 )nu[nβ3]β( 1 3 )nβ6u[nβ4]. Express your answer in the form y[n] = c1Ξ³ nβN1 1 u[n β *N*1] + *c*2Ξ³ nβN2 2 u[n β N2], making sure to clearly identify the constants c1, c2, Ξ³1, Ξ³2, N1, and N2.
Figure P5.2-5
5.2-7 Determine the inverse unilateral z-transform x[n] of the signal
- 5.2-8 Using only the fact that Ξ³ nu[n]ββz/(zβΞ³ ) and properties of the z-transform, find the z-transform of each of the following:
- (a) n2u[n]
- (b) *n*2Ξ³ nu[n]
- (c) n3u[n]
- (d) an[u[n] βu[nβm]]
- (e) neβ2nu[nβm]
- (f) (nβ2)(0.5)nβ3 u[nβ4]
- 5.2-9 Using only pair 1 in Table 5.1 and appropriate properties of the z-transform, derive iteratively pairs 2 through 9. In other words, first derive pair 2. Then use pair 2 (and pair 1, if needed) to derive pair 3, and so on.
- 5.2-10 Find the z-transform of cos(Οn/4)u[n] using only pairs 1 and 11b in Table 5.1 and a suitable property of the z-transform.
- 5.2-11 Apply the time-reversal property to pair 6 of Table 5.1 to show that Ξ³ nu[β(n + 1)] ββ βz/(zβΞ³ ) and the ROC is given by |z| < |Ξ³ |.
- 5.2-12 (a) If x[n] ββ X[z], then show that (β1)nx[n] ββ X[βz].
- (b) Use this result to show that (βΞ³ )nu[n] ββ z/(z +Ξ³ ).
- (c) Use these results to find the z-transforms of xi[n]=[2nβ1 β(β2)nβ1]u[n] and xβii[n] = Ξ³ n cosΟn u[n]
- 5.2-13 (a) If x[n] ββ X[z], then show that
578 CHAPTER 5 DISCRETE-TIME SYSTEM ANALYSIS USING THE Z-TRANSFORM
- (b) Use this result to derive pair 2 from pair 1 in Table 5.1.
- 5.2-14 A number of causal time-domain functions are shown in Fig. P5.2-14. List the function of time that corresponds to each of the following functions of z. Few or no calculations are necessary! Be careful, the graphs may be scaled differently.
(a)
\n(b)
\n(c)
\n(d)
\n
(e)
(c) 0.75z
(f)
\n(g)
(g)
\n(h)
\n(i)
(j)
- 5.2-15 Suppose we upsample a causal signal x[n] by factor N to produce signal y[n]. Express Y(z) in terms of X(z), taking care to mathematically justify your result.
- 5.3-1 Using z-transform techniques, find the output y[n] of an LTID system specified by the equation y[n]β 1 3 y[nβ1] = x[nβ1] if the initial condition is y[β1] = 2 and the input is x[n]=βu[n].
- 5.3-2 Consider an LTID system y[n] β y[n β 2] = x[n] with y[β1] = 0, y[β2] = 1, and x[n] = u[n].
- (a) Determine Y[z], expressed as a rational function in standard factored form.
- (b) Use Y[z] from part (a) to solve for the system output y[n].
- 5.3-3 Solve Prob. 3.8-23 by the z-transform method.
- 5.3-4 Consider a DT system with transfer function H[z] = 2zβ2 zβ0.5 . Assuming the system is both controllable and observable, determine the ZIR yzir[n] given y[β1] = 1.
- 5.3-5 (a) Solve
when
and
-
(b) Find the zero-input and the zero-state components of the response.
-
5.3-6 Consider a LTID system that is described by the difference equation y[n] β 1 4 y[n β 2] = x[nβ1].
-
(a) Use transform-domain techniques to determine the zero-state response yzsr[n] to input x[n] = 3u[nβ5].
-
(b) Use transform-domain techniques to determine the zero-input response yzir[n] given yzir[β2] = yzir[β1] = 1.
-
5.3-7 (a) Find the output y[n] of an LTID system specified by the equation
=
for input x[n] = (4)βnu[n] and initial conditions y[β1] = 0 and y[β2] = 1.
- (b) Find the zero-input and the zero-state components of the response.
- (c) Find the transient and the steady-state components of the response.
- 5.3-8 Solve Prob. 5.3-7 if initial conditions y[β1] and y[β2] are instead replaced with auxiliary conditions y[0] = 3/2 and y[1] = 35/4.
- 5.3-9 (a) Solve
with y[β1] = 0, y[β2] = 1, and x[n] = u[n].
- (b) Find the zero-input and the zero-state components of the response.
- (c) Find the transient and the steady-state components of the response.
- 5.3-10 Solve
if y[β1] = 2, y[β2] = 3, and x[n] = (3)nu[n].
5.3-11 Solve
with y[β1] = 1, y[β2] = 0, and x[n] = u[n].
-
5.3-12 Consider a causal LTID system described as H(z) = 21(z2+1) 16(z2+ 1 4 zβ 3 8 ) .
- (a) Determine the standard delay-form difference equation description of this system.
- (b) Using transform-domain techniques, determine the system impulse response h[n].
- (c) Using transform-domain techniques, determine yzir[n] given y[β1] = 16 and y[β2] = 8.
-
5.3-13 Consider a causal LTID system described as y[n] β 5 6 y[n β 1] + 1 6 y[n β 2] = 3 2 x[n β 1] + 3 2 x[nβ2].
- (a) Determine the (standard-form) system transfer function H(z) and sketch the system pole-zero plot.
- (b) Using transform-domain techniques, determine yzir[n] given y[β1] = 2 and y[β2]=β2.
-
5.3-14 Solve
=
with y[0] = 0, y[1] = 1, and x[n] = enu[n].
- 5.3-15 A system with impulse response h[n] = 2(1/3)nu[n β 1] produces an output y[n] = (β2)nu[n β 1]. Determine the corresponding input x[n].
- 5.3-16 A professor recently received an unexpected $10 (a futile bribe attached to a test). Being the savvy investor that she is, the professor decides to invest the $10 into a savings account that earns 0.5% interest compounded monthly (6.17% APY). Furthermore, she decides to supplement this initial investment with an additional $5 deposit made every month, beginning the month immediately following her initial investment.
- (a) Model the professorβs savings account as a constant coefficient linear difference equation. Designate y[n] as the account balance at month n, where n = 0 corresponds to the first month that interest is awarded (and that her $5 deposits begin).
- (b) Determine a closed-form solution for y[n]. That is, you should express y[n] as a function only of n.
- (c) If we consider the professorβs bank account as a system, what is the system impulse response h[n]? What is the system transfer function H[z]?
- (d) Explain this fact: if the input to the professorβs bank account is the everlasting exponential x[n] = 1n = 1, then the output is not y[n] = 1nH[1] = H[1].
- 5.3-17 Sally deposits $100 into her savings account on the first day of every month except for each December, when she uses her money to buy
holiday gifts. Define b[m] as the balance in Sallyβs account on the first day of month m. Assume Sally opens her account in January (m = 0), continues making monthly payments forever (except each December!), and that her monthly interest rate is 1%. Sallyβs account balance satisfies a simple difference equation b[m] = (1.01)b[m β 1] + p[m], where p[m] designates Sallyβs monthly deposits. Determine a closed-form expression for b[m] that is only a function of the month m.
- 5.3-18 For each impulse response, determine the number of system poles, whether the poles are real or complex, and whether the system is BIBO-stable.
- (a) *h*1[n] = (β1+(0.5)n)u[n]
- (b) *h*2[n] = (j)n(u[n] βu[nβ10])
- 5.3-19 Find the following sums:
(a)
(b)
k=0 [Hint: Consider a system whose output y[n] is the desired sum. Examine the relationship between y[n] and y[n β 1]. Note also that y[0] = 0.]
5.3-20 Find the following sum:
[Hint: See the hint for Prob. 5.3-19.]
5.3-21 Find the following sum:
[Hint: See the hint for Prob. 5.3-19.]
- 5.3-22 Redo Prob. 5.3-19 using the result in Prob. 5.2-13a.
- 5.3-23 Redo Prob. 5.3-20 using the result in Prob. 5.2-13a.
- 5.3-24 Redo Prob. 5.3-21 using the result in Prob. 5.2-13a.
5.3-25 (a) Find the zero-state response of an LTID system with transfer function
and the input x[n] = e(n+1) u[n].
(b) Write the difference equation relating the output y[n] to input x[n].
5.3-26 Repeat Prob. 5.3-25 for x[n] = u[n] and
5.3-27 Repeat Prob. 5.3-25 for
- and the input x[n] is (a) (4)βnu[n] (b) (4)β(nβ2) u[nβ2] (c) (4)β(nβ2) u[n] (d) (4)βnu[nβ2]
- 5.3-28 Repeat Prob. 5.3-25 for x[n] = u[n] and
- 5.3-29 Find the transfer functions corresponding to each of the systems specified by difference equations in Probs. 5.3-5, 5.3-7, 5.3-9, and 5.3-14.
- 5.3-30 Find h[n], the unit impulse response of the systems described by the following equations:
- (a) y[n] +3y[nβ1] +2y[nβ2] = x[n] +3x[nβ 1] +3x[nβ2]
- (b) y[n + 2] + 2y[n + 1] + y[n] = 2x[n + 2] β x[n+1]
- (c) y[n]βy[nβ1]+0.5y[nβ2] = x[n]+2x[nβ 1]
- 5.3-31 Find h[n], the unit impulse response of the systems in Probs. 5.3-25, 5.3-26, and 5.3-28.
- 5.3-32 A system has impulse response h[n] = u[n β 3].
- (a) Determine the impulse response of the inverse system hβ1[n].
- (b) Is the inverse stable? Is the inverse causal?
- (c) Your boss asks you to implement hβ1[n] to the best of your ability. Describe your realizable design, taking care to identify any deficiencies.
5.4-1 A system has impulse response given by
This system can be implemented according to Fig. P5.4-1.
Figure P5.4-1
- (a) Determine the coefficients A1 and A2 to implement h[n] using the structure shown in Fig. P5.4-1.
- (b) What is the zero-state response y0[n] of this system, given a shifted unit step input x[n] = u[n+3]?
- 5.4-2 (a) Show the canonic direct form, a cascade, and a parallel realization of
- (b) Find the transpose of the realizations obtained in part (a).
- 5.4-3 Repeat Prob. 5.4-2 for
5.4-4 Repeat Prob. 5.4-2 for
5.4-5 Repeat Prob. 5.4-2 for
5.4-6 Repeat Prob. 5.4-2 for
5.4-7 Consider the LTID system shown in Fig. P5.4-7.
- Figure P5.4-7
- (a) Determine the standard delay-form difference equation description of this system.
- (b) Determine the impulse response h[n] of this system.
- (c) Is this realization canonical? Explain.
- (d) Is this system stable? Explain.
- (e) Is this system causal? Explain.
- 5.4-8 Realize a system whose transfer function is
5.4-9 Realize a system whose transfer function is given by
5.4-10 Consider the LTID system shown in Fig. P5.4-10, where parameter c is an arbitrary, real constant.
Figure P5.4-10
- (a) Determine the system transfer function H[z], expressed in standard rational form.
- (b) Determine all system poles and all system zeros.
- (c) Is the system of Fig. P5.4-10 canonical? Explain.
- (d) What constraints, if any, exist on parameter c to ensure that the system is stable?
- 5.4-11 This problem demonstrates the enormous number of ways of implementing even a relatively low-order transfer function. A second-order transfer function has two real zeros and two real poles. Discuss various ways of realizing
such a transfer function. Consider canonic direct, cascade, parallel, and the corresponding transposed forms. Note also that interchange of cascaded sections yields a different realization.
5.4-12 Consider a digital audio system: an input analog-to-digital converter (ADC) is used to collect input samples at a CD-quality rate Fs = 44 kHz. Input samples are processed with a digital filter to generate output samples, which are sent to a digital-to-analog converter (DAC) at the same rate Fs. Every sample interval T = 1/Fs, the digital processor executes the following MATLAB-compatible code:
% Read input sample from the ADCx = read_ADC;% Process input and...mem(1) = x - mem(3)*9/16;% ...compute output sampley = mem(1)*7/16 - mem(3)*7/16;% Send output sample to the DACwrite_DAC = y;% Update memory for next iterationmem(3) = mem(2);mem(2) = mem(1);- (a) Does the code implement DFI, DFII, TDFI, or TDFII? Support your answer by drawing the appropriate block diagram labeled in a manner that is consistent with the code.
- (b) Determine the transfer function H[z] of this system.
- (c) What is the basic filtering function of this system: LP, HP, BP, or BS? Justify your answer.
- (d) Determine the transfer function Hβ1[z] of the inverse system to H[z] and draw its DFI block implementation. How well will the inverse system operate?
- 5.4-13 Repeat Prob. 5.4-12 but instead use the code:
% Read input sample from the ADCx = read_ADC;% Compute output sampley = x*7/32+mem(1);% Send output sample to the DACwrite_DAC = y;% Update memory for next iterationmem(1) = mem(2);mem(2) = x*7/32 + y*9/16;- 5.5-1 A CT sinusoid x(t) = cos(Οt) is sampled at a greater-than-Nyquist rate Fs = 1000 Hz to produce a DT sinusoid x(t) = cos(n). Determine the analog frequency Ο if (a) = Ο 4
- (b) = 2Ο 3
- (c) = 7 8
- 5.5-2 Find the amplitude and phase response of the digital filters depicted in Fig. P5.5-2.
- 5.5-3 A causal LTID system H(z) = 21(zβj)(z+j) 16(zβ 1 2 )(z+ 3 4 ) has a periodic input x[n] that toggles between the values 1 and 2. That is, x[n]=[β¦, 1, 2, 1, β 2 , 1, 2, 1, β¦], where x[0] = 2.
- (a) Plot the magnitude response |H(ej)| over β2Ο β€ β€ 2Ο.
- (b) Plot the phase response H(ej) over β2Ο β€ β€ 2Ο.
- (c) Determine the system output y[n] in response to the periodic input x[n].
5.5-4 A causal LTID system H(z) = β7(z+1) 32(zβj 3 4 )(z+j 3 4 ) has a periodic input x[n] that cycles through the 4 values 3, 2, 1, and 2. That is, x[n] = [β¦, 3, 2, 1, 2, β 3, 2, 1, 2, β¦], where x[0] = 3.
- (a) Plot the magnitude response |H(ej)| over β2Ο β€ β€ 2Ο.
- (b) Plot the phase response H(ej) over β2Ο β€ β€ 2Ο.
- (c) Determine the system output y[n] in response to the periodic input x[n].
- 5.5-5 Find the amplitude and the phase response of the filters shown in Fig. P5.5-5. [Hint: Express H[ej] as eβj2.5Ha[ej].]
- 5.5-6 Find the frequency response for the moving-average system in Prob. 3.4-3. The inputβoutput equation of this system is given by
5.5-7 (a) Inputβoutput relationships of two filters are described by
(i) y[n]=β0.9y[nβ1] +x[n]
(ii) y[n] = 0.9y[nβ1] +x[n]
For each case, find the transfer function, the amplitude response, and the phase response. Sketch the amplitude response, and state the type (highpass, lowpass, etc.) of each filter.
- (b) Find the response of each of these filters to a sinusoid x[n] = cosn for = 0.01Ο and 0.99Ο. In general, show that the gain (amplitude response) of filter (i) at frequency 0 is the same as the gain of filter (ii) at frequency Ο β0.
- 5.5-8 For an LTID system specified by the equation
(a) Find the amplitude and the phase response.
- (b) Find the system response y[n] for the input x[n] = cos(0.5k β(Ο/3)).
- 5.5-9 For an asymptotically stable LTID system, show that the steady-state response to input ejnu[n] is H[ej]ejnu[n]. The steady-state response is that part of the response which does not decay with time and persists forever.
- 5.5-10 Express the following signals in terms of apparent frequencies:
- (a) cos(0.8Οn+ΞΈ )
- (b) sin(1.2Οn+ΞΈ )
- (c) cos(6.9n+ΞΈ )
Figure P5.5-2
(b)
Figure P5.5-5
- (d) cos(2.8Οn+ΞΈ ) +2 sin(3.7Οn+ΞΈ )
- (e) sinc(Οn/2)
- (f) sinc(3Οn/2)
- (g) sinc(2Οn)
- 5.5-11 Show that cos β (0.6Οn + (Ο/6)) + 3 cos(1.4Οn + (Ο/3)) = 2 cos(0.6Οn β (Ο/6)).
- 5.5-12 (a) A digital filter has the sampling interval T = 50Β΅s. Determine the highest frequency that can be processed by this filter without aliasing.
- (b) If the highest frequency to be processed is 50 kHz, determine the minimum value of the sampling frequency Fs and the maximum value of the sampling interval T that can be used.
- 5.5-13 Consider the discrete-time system represented by
- (a) Determine and plot the magnitude response |H[ej]| of the system.
- (b) Determine and plot the phase response H[ej] of the system.
- (c) Find an efficient block representation that implements this system.
- 5.6-1 Pole-zero configurations of certain filters are shown in Fig. P5.6-1. Sketch roughly the amplitude response of these filters.
- 5.6-2 Figure P5.6-2 displays the pole-zero plot of a second-order real, causal LTID system that has H[β1]=β1.
Figure P5.6-2
- (a) Determine the five constants b0, b1, b2, a1, and a2 that specify the transfer function H[z] = b0z2+b1z+b2 z2+a1z+a2 .
- (b) Using the techniques of Sec. 5.6, accurately hand-sketch the system magnitude response |H[ej]| over the range (β2Ο β€ β€ 0).
- (c) Determine the output y[n] of this system if the input is x[n] = sin Οn 2 .
- 5.6-3 Repeat Prob. 5.6-2 if the zero at z = 1 is moved to z = β1 and H[1]=β1 is specified rather than H[β1]=β1.
- 5.6-4 Figure P5.6-4 displays the pole-zero plot of a second-order real, causal LTID system that has H[β1] = 1.
- (a) Determine the five constants k, b1, b2, a1, and a2 that specify the transfer function H[z] = k z2+b1z+b2 z2+a1z+a2 .
Figure P5.6-4
- (b) Using the techniques of Sec. 5.6, accurately hand-sketch the system magnitude response |H[ej]| over the range (βΟ β€ β€ Ο ).
- (c) A signal x(t) = cos(100Οt) + cos(500Οt) is sampled at a greater than Nyquist rate Fs Hz and then input into the above LTID system to produce DT output y[n] = Ξ² cos(0n + ΞΈ ). Determine Fs and 0. You do not need to find constants Ξ² and ΞΈ.
- (d) Is the inpulse response h[n] of this system absolutely summable? Justify your answer.
- 5.6-5 Figure P5.6-5 displays the pole-zero plot of a second-order real, causal LTID system that has H[1]=β1.
(a) Determine the five constants k, b1, b2, a1, and a2 that specify the transfer function H[z] = k z2+b1z+b2 z2+a1z+a2 .
Figure P5.6-5
- (b) Using the techniques of Sec. 5.6, accurately hand-sketch the system magnitude response |H[ej]| over the range (βΟ β€ β€ Ο ).
- (c) A signal x(t) = cos(2Οft) is sampled at a rate Fs = 1 kHz and then input into the above LTID system to produce DT output y[n]. Determine, if possible, the frequency or frequencies f that will produce zero output, y[n] = 0.
- 5.6-6 The system y[n] β y[n β 1] = x[n] β x[n β 1] is an all-pass system that has zero phase response. Is there any difference between this system and the system y[n] = x[n]? Justify your answer.
- 5.6-7 Figure P5.6-7 displays the pole-zero plot of a second-order real, causal LTID system that has a repeated zero and H[1] = 4. The solid circle is the unit circle.
Figure P5.6-7
- (a) Determine the five constants k, b1, b2, a1, and a2 that specify the transfer function H[z] = k z2+b1z+b2 z2+a1z+a2 .
- (b) Using the techniques of Sec. 5.6, accurately hand-sketch the system magnitude response |H[ej]| over the range (βΟ β€ β€ Ο ).
- (c) Determine the steady-state output yss[n] of this system if the input is x[n] = cos( 3Ο*n* 4 )u[n].
- (d) State whether this system is LP, HP, BP, BS, or other. If the digital system operates at Fs = 8 kHz, what is the approximate hertzian cutoff frequency (or frequencies) of this system?
- 5.6-8 The magnitude and phase responses of a real, stable, LTI system are shown in Fig. P5.6-8.
- (a) What type of system is this: lowpass, highpass, bandpass, or bandstop?
- (b) What is the output of this system in response to
(c) What is the output of this system in response to
-
5.6-9 Consider an LTID system with system function H[z] = b0 *z*2+1 *z*2β9/16 .
- (a) Determine the constant b0 so that the system frequency response at = βΟ is β1.
- (b) Accurately sketch the system poles and zeros.
- (c) Using the locations of the system poles and zeros, sketch |H[ej]| over 0 β€ β€ 2Ο.
-
(d) Determine the response y[n] to the input x[n] = (β1+j)+j n +(1βj)sin(Οn+1).
-
(e) Draw an appropriate block diagram representation of this system.
-
5.6-10 Do Prob. 5.10-3 by graphical procedure. Do the sketches approximately, without using MATLAB.
-
5.6-11 Do Prob. 5.10-8 by graphical procedure. Do the sketches approximately, without using MATLAB.
-
5.6-12 (a) Realize a digital filter whose transfer function is given by
- (b) Sketch the amplitude response of this filter, assuming |a| < 1.
- (c) The amplitude response of this lowpass filter is maximum at = 0. The 3 dB bandwidth is the frequency at which the amplitude response drops to 0.707 (or 1/ β2) times its maximum value. Determine the 3 dB bandwidth of this filter when a = 0.2.
- 5.6-13 Design a digital notch filter to reject frequency 5000 Hz completely and to have a sharp recovery on either side of 5000 Hz to a gain of unity. The highest frequency to be processed is 20 kHz (Fh = 20,000). [Hint: See Ex. 5.15. The zeros should be at eΒ±jΟ*T* for Ο corresponding to 5000 Hz, and the poles are at aeΒ±jΟ*T* with a < 1. Leave your answer in terms of a. Realize this filter using the canonical form. Find the amplitude response of the filter.]
5.6-14 Consider the desired DT system magnitude response |H[ej]| in Fig. P5.6-14.
Figure P5.6-14
- (a) Is the filter LP, HP, BP, BS, or other? Explain.
- (b) Sketch the pole-zero plot of a 2nd-order system that behaves as a reasonable approximation of Fig. P5.6-14. What is the coefficient b0?
- 5.6-15 Show that a first-order LTID system with a pole at z = r and a zero at z = 1/r (r β€ 1) is an allpass filter. In other words, show that the amplitude response |H[ej]| of a system with the transfer function
is constant with frequency. This is a first-order allpass filter. [Hint: Show that the ratio of the distances of any point on the unit circle from the zero (at z = 1/r) and the pole (at z = r) is a constant 1/r.]
Generalize this result to show that an LTID system with two poles at z = reΒ±jΞΈ and two zeros at z = (1/r)eΒ±jΞΈ (r β€ 1) is an allpass filter. In other words, show that the amplitude response of a system with the transfer function
is constant with frequency.
5.6-16 (a) If h1[n] and h2[n], the impulse responses of two LTID systems are related by h2[n] = (β1)nh1[n], then show that
How is the frequency response spectrum *H*2[ej] related to the *H*1[ej]?
- (b) If H1[z] represents an ideal lowpass filter with cutoff frequency c, sketch *H*2[ej]. What type of filter is *H*2[ej]?
- 5.6-17 Mappings such as the bilinear transformation are useful in the conversion of continuous-time filters to discrete-time filters. Another useful type of transformation is one that converts a discrete-time filter into a different type of discrete-time filter. Consider a transformation that replaces z with βz.
- (a) Show that this transformation converts lowpass filters into highpass filters and highpass filters into lowpass filters.
- (b) If the original filter is an FIR filter with impulse response h[n], what is the impulse response of the transformed filter?
- 5.6-18 The bilinear transformation is defined by the rule s = 2(1βzβ1)/T(1+zβ1).
- (a) Show that this transformation maps the Ο axis in the s plane to the unit circle z = ej in the z plane.
- (b) Show that this transformation maps to 2 arctan(ΟT/2).
- 5.7-1 In Ch. 3, we used another approximation to find a digital system to realize an analog system. We showed that an analog system specified by Eq. (3.12) can be realized by using the digital system specified by Eq. (3.13). Compare that solution with the one resulting from the
impulse-invariance method. Show that one result is a close approximation of the other and that the approximation improves as T β 0.
- 5.7-2 A CT system has impulse response hct(t) = eβt u(t). Draw the DFI realization of the corresponding DT system designed by the impulse-invariance method with T = 0.1.
- 5.7-3 (a) Using the impulse-invariance criterion, design a digital filter to realize an analog filter with transfer function
- (b) Show a canonical and a parallel realization of the filter. Use a 1% criterion for the choice of T.
- 5.7-4 Use the impulse-invariance criterion to design a digital filter to realize the second-order analog Butterworth filter with transfer function
Use a 1% criterion for the choice of T.
-
5.7-5 Design a digital integrator using the impulse-invariance method. Find and give a rough sketch of the amplitude response, and compare it with that of the ideal integrator. If this integrator is used primarily for integrating audio signals (whose bandwidth is 20 kHz), determine a suitable value for T.
-
5.7-6 An oscillator by definition is a source (no input) that generates a sinusoid of a certain frequency Ο0. Therefore, an oscillator is a system whose zero-input response is a sinusoid of the desired frequency. Find the transfer function of a digital oscillator to oscillate at 10 kHz by the methods described in parts (a) and (b). In both methods, select T so that there are 10 samples in each cycle of the sinusoid.
- (a) Choose H[z] directly so that its zero-input response is a discrete-time sinusoid of frequency =ΟT corresponding to 10 kHz.
- (b) Choose Ha(s) whose zero-input response is an analog sinusoid of 10 kHz. Now use the impulse invariance method to determine H[z].
- (c) Show a canonical realization of the oscillator.
-
5.7-7 A variant of the impulse invariance method is the step-invariance method of digital filter synthesis. In this method, for a given Ha(s), we design H[z] in Fig. 5.24a such that y(nT) in Fig. 5.24b is identical to y[n] in Fig. 5.24a when x(t) = u(t).
- (a) Show that, in general,
(b) Use this method to design H[z] for
- (c) Use the step-invariance method to synthesize a discrete-time integrator and compare its amplitude response with that of the ideal integrator.
- 5.7-8 Use the ramp-invariance method to synthesize a discrete-time differentiator and integrator. In this method, for a given Ha(s), we design H[z] such that y(nT) in Fig. 5.24b is identical to y[n] in Fig. 5.24a when x(t) = tu(t).
- 5.7-9 In an impulse-invariance design, show that if Ha(s) is a transfer function of a stable system, the corresponding H[z] is also a transfer function of a stable system.
- 5.7-10 First-order backward differences provide the transformation rule s = (1βzβ1)/T.
- (a) Show that this transformation maps the Ο axis in the s plane to a circle of radius 1/2 centered at (1/2, 0) in the z plane.
- (b) Show that this transformation maps the left-half s plane to the interior of the unit circle in the z plane, which ensures that stability is preserved.
- 5.8-1 Find the z-transform (if it exists) and the corresponding ROC for each of the following signals:
- (a) (0.8)nu[n] +2nu[β(n+1)]
(b)
(c)
(d)
(e)
(f) [(0.8)n +3(0.4)n]u[n]
(g)
(h)
- 5.8-2 Using the definition, compute the bilateral z-transform X(z) of (a) x[n] = 3nu[βn] (b) x[n] = ( 1 3 )nu[n] Express your answers in standard rational form.
- 5.8-3 Determine the inverse z-transform x[n] of X[z] = *z*2β 1 3 z (zβ1)(z+2) with ROC 1 < |z| < 2.
5.8-4 Find the inverse
-transform of
when the ROC is,
\n(a)
\n(b)
\n(c)
5.8-5 Use partial fraction expansions, z-transform tables, and a region of convergence (|z| < 1/2) to determine the inverse z-transform of
- 5.8-6 Using z-transform techniques and properties (no time-domain convolution sum!), determine the convolution y[n] = ( 1 3 )nβ3u[n β 2] β (2)nu[βn]. Express your answer in the form y[n] = c1Ξ³ n 1 u[n + *N*1] + *c*2Ξ³ n 2 u[βn + N2], making sure to clearly identify the constants c1, c2, Ξ³1, Ξ³2, N1, and N2.
- 5.8-7 Using partial fraction expansions, z-transform tables, and the fact that h[n] is stable, determine the inverse z-transform of
5.8-8 Consider the system
- (a) Draw the pole-zero diagram for H[z] and identify all possible regions of convergence.
- (b) Draw the pole-zero diagram for Hβ1[z] and identify all possible regions of convergence.
- 5.8-9 A discrete-time signal x[n] has a rational z-transform that contains a pole at z = 0.5. Given *x*1[n] = (1/3)nx[n] is absolutely summable and
Problems 589
*x*2[n] = (1/4)nx[n] is not absolutely summable, determine whether x[n] is left-sided, right-sided, or two-sided. Justify your answer!
- 5.8-10 Let x[n] be an absolutely summable signal with rational z-transform X[z]. X[z] is known to have a pole at z = (0.75+0.75j), and other poles may be present. Recall that an absolutely summable signal satisfies %β ββ |x[n]| < β.
- (a) Can x[n] be left-sided? Explain.
- (b) Can x[n] be right-sided? Explain.
- (c) Can x[n] be two-sided? Explain.
- (d) Can x[n] be of finite duration? Explain.
- 5.8-11 Consider a causal system that has transfer function
When appropriate, assume initial conditions of zero.
- (a) Determine the output y1[n] of this system in response to *x*1[n] = (3/4)nu[n].
- (b) Determine the output y2[n] of this system in response to *x*2[n] = (3/4)n.
- (c) Determine the output y3[n] of this system in response to *x*3[n] = (3/4)nu[βnβ1].
- 5.8-12 Let x[n] =(β1)nu[nβn0]+Ξ±nu[βn]. Determine the constraints on the complex number Ξ± and the integer n0 so that the z-transform X[z] exists with region of convergence 1 < |z| < 2.
- 5.8-13 Using the definition, compute the bilateral z-transform, including the region of convergence (ROC), of the following complex-valued functions:
- (a) *x*1[n] = (βj)βnu[βn] +Ξ΄[βn]
- (b) x2[n] = (j)n cos(n+1)u[n]
- (c) x3[n] = jsinh(n)u[βn+1]
- (d) x4[n] = %0 k=ββ(2j)nΞ΄[nβ2k]
- 5.8-14 Use partial fraction expansions, z-transform tables, and a region of convergence (0.5 < |z| < 2) to determine the inverse z-transform of
- (a) X1[z] = 1 1+ 13 6 zβ1 + 1 6 zβ2 β 1 3 zβ3 1
(b)
5.8-15 Use partial fraction expansions, z-transform tables, and the fact that the systems are stable to determine the inverse z-transform of
(a)
(b)
5.8-16 By inserting N β 1 zeros between every sample of a unit step, we obtain a signal
Determine H[z], the bilateral z-transform of h[n]. Identify the number and location(s) of the poles of H[z].
- 5.8-17 Using transform-domain techniques, determine the zero-state response yzsr[n] of LTID system y[n] β 1 4 y[n β 2] = x[n] to the noncausal input x[n] = 2nu[2βn].
- 5.8-18 Determine the zero-state response of a system having a transfer function
and an input x[n] given by
- (a) x[n] = enu[n]
- (b) x[n] = 2nu[β(n+1)]
- (c) x[n] = enu[n] +2nu[β(n+1)]
- (d) x[n] = 2nu[n] +u[β(n+1)]
- (e) x[n] = eβ2nu[β(n+1)]
- 5.8-19 The discrete cross-correlation between real signal x[n] and real signal y[n] is
Let signal x[n] have z-transform X[z] with ROC Rx, and let signal y[n] have z-transform Y[z] with ROC Ry. Determine CXY [z] (the bilateral z-transform of cxy[n]) in terms of the z-transforms of x[n] and y[n].
5.8-20 Transform properties can be very useful. The accumulation property states
This property is the DT version of the Laplace transform βintegration in timeβ property. Given x[n] ββ X[z], prove the accumulation property. [Hint: Polynomial long division can be helpful.] 5.10-1 Use MATLAB to generate pole-zero plots for the causal systems with the following transfer functions:
(a)
(b)
3+ 3 2 zβ2+2zβ4 In each case, determine whether the system is
- stable.
- 5.10-2 For each of the following stable LTID systems, use MATLAB to generate magnitude and phase response plots over βΟ β€ <Ο. (a) *H*a[z] = cos(z) zβ0.5
- (b) Hb[z] = z3 sin(zβ1)
- 5.10-3 Consider an LTID system described by the difference equation 4y[n + 2] β y[n] = x[n + 2] +x[n].
- (a) Plot the pole-zero diagram for this system.
- (b) Plot the systemβs magnitude response |H[ej]| over βΟ β€ β€ Ο.
- (c) What type of system is this: lowpass, highpass, bandpass, or bandstop?
- (d) Is this system stable? Justify your answer.
- (e) Is this system real? Justify your answer.
- (f) If the system input is of the form x[n] = cos(n), what is the greatest possible amplitude of the output? Justify your answer.
- (g) Draw an efficient, causal implementation of this system using only add, scale, and delay blocks.
- 5.10-4 Consider an LTID system with system function H(z) = b0 *z*2+1 *z*2β4/9 .
- (a) Determine the constant b0 so that the system frequency response at = βΟ is β1.
- (b) Plot the pole-zero diagram for this system.
- (c) Plot the systemβs magnitude response |H[ej]| over βΟ β€ β€ Ο.
- (d) Plot the systemβs phase response H[ej] over βΟ β€ β€ Ο.
- (e) What type of system is this: lowpass, highpass, bandpass, or bandstop?
- (f) Determine the response y[n] to the input x[n] = (β1 β j) + (βj)n + (1 β j) cos(Ο*n* + 1 3 ).
- (g) Draw a TDFII block diagram representation of this system.
5.10-5 Consider the LTID system shown in Fig. P5.10-5, where parameters c1 and c2 are arbitrary constants.
Figure P5.10-5
- (a) Determine the system function H[z], expressed in standard rational form.
- (b) What is the order N of this system?
- (c) Determine the N poles and N zeros of this system. Use MATLAB to create the corresponding pole-zero plot.
- (d) What constraints, if any, exist on parameters c1 and c2 to ensure that the system is stable?
- (e) Determine c1 and c2 so that this system functions as an LPF with narrow passband. Use MATLAB to generate the corresponding magnitude response |H[ej]| over βΟ β€ β€ Ο.
- (f) Ms. Zeroine, the heroine of DT systems, believes that if x[n] = 0, then y[n] = 0 also. Is Ms. Zeroine correct? Fully justify your answer.
- (g) Dr. Strange suggests that by setting c1 = β1 and c2 = β2, the system will act as a highpass filter. Is Dr. Strange right? Fully justify your answer.
- 5.10-6 One interesting and useful application of discrete systems is the implementation of complex (rather than real) systems. A complex system is one in which a real-valued input can produce a complex-valued output. Complex systems that are described by constant coefficient difference equations require at least one complex-valued coefficient, and they are capable of operating on complex-valued inputs. Consider the complex discrete-time system
(a) Determine and plot the system zeros and poles.
- (b) Sketch the magnitude response |H[ejΟ]| of this system over β2Ο β€ Ο β€ 2Ο. Comment on the systemβs behavior.
- 5.10-7 Consider the complex system
Refer to Prob. 5.10-6 for an introduction to complex systems.
- (a) Plot the pole-zero diagram for H[z].
- (b) Plot the systemβs magnitude response |H[ej]| over βΟ β€ β€ Ο.
- (c) Explain why H[z] is a noncausal system. Do not give a general definition of causality; specifically identify what makes this system noncausal.
- (d) One way to make this system causal is to add two poles to H[z]. That is,
Find poles a and b such that |Hcausal[ej]| = |H[ej]|.
- (e) Draw an efficient block implementation of Hcausal[z].
- 5.10-8 A discrete-time LTI system is shown in Fig. P5.10-8.
- (a) Determine the difference equation that describes this system.
- (b) Determine the magnitude response |H[ej]| for this system and simplify your answer. Plot the magnitude response over βΟ β€ β€ Ο. What type of standard filter (lowpass, highpass, bandpass, or bandstop) best describes this system?
- (c) Determine the impulse response h[n] of this system.
Figure P5.10-8
- 5.10-9 Determine the impulse response h[n] for the system shown in Fig. P5.10-9. Is the system stable? Is the system causal?
- 5.10-10 An LTID filter has an impulse response function given by h[n] = Ξ΄[n β 1] + Ξ΄[n + 1]. Determine and carefully sketch the magnitude response |H[ej]| over the range βΟ β€ β€ Ο. For this range of frequencies, is this filter lowpass, highpass, bandpass, or bandstop?
- 5.10-11 A causal, stable discrete system has the rather strange transfer function H[z] = cos(zβ1).
- (a) Write MATLAB code that will compute and plot the magnitude response of this system over an appropriate range of digital frequencies . Comment on the system.
- (b) Determine the impulse response h[n]. Plot h[n] over (0 β€ n β€ 10).
- (c) Determine a difference equation description for an FIR filter that closely approximates the system H[z] = cos(zβ1). To verify proper behavior, plot the FIR filterβs magnitude response and compare it with the magnitude response computed in part (a).
- 5.10-12 The MATLAB signal-processing toolbox function butter helps design digital Butterworth filters. Use MATLAB help to learn how butter works. For each of the following cases, design the filter, plot the filterβs poles and zeros in the complex z plane, and plot the decibel magnitude response 20log10 |H[ej]|.
- (a) Design an eighth-order digital lowpass filter with c = Ο/3.
- (b) Design an eighth-order digital highpass filter with c = Ο/3.
- (c) Design an eighth-order digital bandpass filter with passband between 5Ο/24 and 11Ο/24.
- (d) Design an eighth-order digital bandstop filter with stopband between 5Ο/24 and 11Ο/24.
- 5.10-13 The MATLAB signal-processing toolbox function cheby1 helps design digital Chebyshev
type I filters. A Chebyshev type I filter has passband ripple and smooth stopband. Setting the passband ripple to Rp = 3 dB, repeat Prob. 5.10-12 using the cheby1 command. With all other parameters held constant, what is the general effect of reducing Rp, the allowable passband ripple?
5.10-14 The MATLAB signal-processing toolbox function cheby2 helps design digital Chebyshev type II filters. A Chebyshev type II filter has smooth passband and ripple in the stopband. Setting the stopband ripple Rs = 20 dB down, repeat Prob. 5.10-12 using the cheby2 command. With all other parameters held constant, what is the general effect of increasing Rs, the minimum stopband attenuation?
5.10-15 The MATLAB signal-processing toolbox function ellip helps design digital elliptic filters. An elliptic filter has ripple in both the passband and the stopband. Setting the passband ripple to Rp = 3 dB and the stopband ripple Rs = 20 dB down, repeat Prob. 5.10-12 using the ellip command.