Skip to content

$$

← Back to LINEAR SYSTEMS AND SIGNALS Overview

ThisisthetrigonometricformoftheDTFT.9.22Asignalx[n]canbeexpressedasthesumofevenandoddcomponents(Sec.1.51): This is the trigonometric form of the DTFT. **9.2-2** A signal *x*[*n*] can be expressed as the sum of even and odd components (Sec. 1.5-1):

x[n] = x_e[n] + x_o[n]

(a)Ifx[n]X(),showthatforrealx[n], (a) If *x*[*n*] ⇐⇒ *X*(), show that for real *x*[*n*],

x_e[n] \Longleftrightarrow \operatorname{Re}[X(\Omega)]

and and

x_o[n] \Longleftrightarrow j \operatorname{Im}[X(\Omega)]

(b)VerifytheseresultsbyfindingtheDTFToftheevenandoddcomponentsofthesignal(0.8)nu[n].9.23Forthefollowingsignals,findtheDTFTdirectly,usingthedefinitioninEq.(9.19).Assumeγ<1.(a)δ[n](b) - (b) Verify these results by finding the DTFT of the even and odd components of the signal (0.8)*nu*[*n*]. - **9.2-3** For the following signals, find the DTFT directly, using the definition in Eq. (9.19). Assume |γ | < 1. - (a) δ[*n*] (b)

\delta[n-k]

(c)<sup>γ</sup>nu[n1](d)<sup>γ</sup>nu[n+1](e)(γ)nu[n](f)γ<sup></sup>n<sup></sup>9.24UseEq.(9.18)tofindtheinverseDTFTforthefollowingspectra,givenonlyovertheintervalπ.Assume<sup>c</sup>and<sup>0</sup><π.(a)ejkintegerk(b)coskintegerk(c)cos2(/2)(d) - (c) <sup>γ</sup> *nu*[*n*−1] - (d) <sup>γ</sup> *nu*[*n*+1] - (e) (−γ )*nu*[*n*] - (f) γ <sup>|</sup>*n*<sup>|</sup> - **9.2-4** Use Eq. (9.18) to find the inverse DTFT for the following spectra, given only over the interval || ≤ π. Assume *<sup>c</sup>* and <sup>0</sup> < π. (a) *ejk* integer *k* - (b) cos*k* integer *k* - (c) cos2(/2) (d)

\Delta \left( \frac{1}{2\Omega_c} \right)

(e)2πδ(0)(f)π[δ(0)+δ(+0)]9.25(a)DetermineandplottheDTFTX()ofthetriangularsignalx[n]showninFig.P9.25.(b)UsingEx.9.5asaguide,useMATLABandtheFFTtovalidatetheDTFTcalculationsandplotofpart(a).FigureP9.259.26UsingEq.(9.18),showthattheinverseDTFTofrect((π/4)/π)is0.5sinc(πn/2)ejπn/4.9.27UsingEq.(9.19),findtheDTFTofthesignalsx[n]inFig.P9.27.9.28UsingEq.(9.19),findtheDTFTofthesignalsdepictedinFig.P9.28.9.29UseEq.(9.18)tofindtheinverseDTFTofthespectra(shownonlyforπ)inFig.P9.29.FigureP9.299.210UseEq.(9.18)tofindtheinverseDTFTofthespectra(shownonlyforπ)inFig.P9.210.9.211FindtheDTFTforthesignalsshowninFig.P9.211.FigureP9.210(c)(d)FigureP9.2119.212FindtheinverseDTFTofX()(shownonlyforπ)forthespectraillustratedinFig.P9.212.[Hint:X()=X()ej<sup>X</sup>().Thisproblemillustrateshowdifferentphasespectra(bothwiththesameamplitudespectrum)represententirelydifferentsignals.]9.213(a)Showthattimeexpandedsignalxe[n]inEq.(3.2)canalsobeexpressedas - (e) 2πδ(−0) - (f) π[δ(−0) +δ(+0)] - **9.2-5** (a) Determine and plot the DTFT *X*() of the triangular signal *x*[*n*] shown in Fig. P9.2-5. - (b) Using Ex. 9.5 as a guide, use MATLAB and the FFT to validate the DTFT calculations and plot of part (a). **Figure P9.2-5** - **9.2-6** Using Eq. (9.18), show that the inverse DTFT of rect((−π/4)/π ) is 0.5 sinc(π*n*/2) *ej*π*n*/4. - **9.2-7** Using Eq. (9.19), find the DTFT of the signals *x*[*n*] in Fig. P9.2-7. - **9.2-8** Using Eq. (9.19), find the DTFT of the signals depicted in Fig. P9.2-8. - **9.2-9** Use Eq. (9.18) to find the inverse DTFT of the spectra (shown only for || ≤ π) in Fig. P9.2-9. **Figure P9.2-9** **9.2-10** Use Eq. (9.18) to find the inverse DTFT of the spectra (shown only for || ≤ π) in Fig. P9.2-10. **9.2-11** Find the DTFT for the signals shown in Fig. P9.2-11. **Figure P9.2-10** (c) (d) **Figure P9.2-11** - **9.2-12** Find the inverse DTFT of *X*() (shown only for ||≤π) for the spectra illustrated in Fig. P9.2-12. [*Hint: X*() = |*X*()|*ej <sup>X</sup>*(). This problem illustrates how different phase spectra (both with the same amplitude spectrum) represent entirely different signals.] - **9.2-13** (a) Show that time-expanded signal *xe*[*n*] in Eq. (3.2) can also be expressed as

x_e[n] = \sum_{k=-\infty}^{\infty} x[k]\delta[n - Lk]

(b)FindtheDTFTofxe[n]byfindingtheDTFToftherighthandsideoftheequationinpart(a).(c)Usetheresultinpart(b)andTable9.1tofindtheDTFTofz[n],showninFig.P9.213.9.214(a)AglanceatEq.(9.18)showsthattheinverseDTFTequationisidenticaltotheinverse(continuoustime)FouriertransformEq.(7.10)forasignalx(t)bandlimitedtoπrad/s.Hence,weshouldbeabletousethecontinuoustimeFouriertransformTable7.1tofindDTFTpairsthatcorrespondtocontinuoustimetransformpairsforbandlimitedsignals.UsethisfacttoderiveDTFTpairs8,9,11,12,13,and14inTable9.1bymeansoftheappropriatepairsinTable7.1.(b)Canthismethodbeusedtoderivepairs2,3,4,5,6,7,10,15,and16inTable9.1?Justifyyouranswerwithspecificreason(s).9.215ArethefollowingfrequencydomainsignalsvalidDTFTs?Answeryesorno,andjustifyyouranswers.(a)X()=+π(b)X()=j+π(c)X()=sin(10)(d)X()=sin(/10)(e)X()=δ()9.31Usingonlypairs2and5(Table9.1)andthetimeshiftingpropertyofEq.(9.31),findtheDTFTofthefollowingsignals,assuminga<1.(a)u[n]u[n9](b)anmu[nm](c)an3(u[n]u[n10])1(b)FigureP9.213(d)<sup>a</sup>nmu[n](e)anu[nm](f)(nm)anmu[nm](g)(nm)anu[n](h)nanmu[nm]9.32Thetriangularpulsex[n]showninFig.P9.32aisgivenby - (b) Find the DTFT of *xe*[*n*] by finding the DTFT of the right-hand side of the equation in part (a). - (c) Use the result in part (b) and Table 9.1 to find the DTFT of *z*[*n*], shown in Fig. P9.2-13. - **9.2-14** (a) A glance at Eq. (9.18) shows that the inverse DTFT equation is identical to the inverse (continuous-time) Fourier transform Eq. (7.10) for a signal *x*(*t*) bandlimited to π rad/s. Hence, we should be able to use the continuous-time Fourier transform Table 7.1 to find DTFT pairs that correspond to continuous-time transform pairs for bandlimited signals. Use this fact to derive DTFT pairs 8, 9, 11, 12, 13, and 14 in Table 9.1 by means of the appropriate pairs in Table 7.1. - (b) Can this method be used to derive pairs 2, 3, 4, 5, 6, 7, 10, 15, and 16 in Table 9.1? Justify your answer with specific reason(s). - **9.2-15** Are the following frequency-domain signals valid DTFT's? Answer yes or no, and justify your answers. - (a) *X*() = +π - (b) *X*() = *j*+π - (c) *X*() = sin(10) - (d) *X*() = sin(/10) - (e) *X*() = δ() - **9.3-1** Using only pairs 2 and 5 (Table 9.1) and the time-shifting property of Eq. (9.31), find the DTFT of the following signals, assuming |*a*| < 1. - (a) *u*[*n*] −*u*[*n*−9] - (b) *an*−*mu*[*n*−*m*] - (c) *an*−3(*u*[*n*] −*u*[*n*−10]) 1 (b) **Figure P9.2-13** - (d) *<sup>a</sup>n*−*mu*[*n*] - (e) *anu*[*n*−*m*] - (f) (*n*−*m*)*an*−*mu*[*n*−*m*] - (g) (*n*−*m*)*anu*[*n*] - (h) *nan*−*mu*[*n*−*m*] - **9.3-2** The triangular pulse *x*[*n*] shown in Fig. P9.3-2a is given by

X(\Omega) = \frac{4e^{j6\Omega} - 5e^{j5\Omega} + e^{j\Omega}}{(e^{j\Omega} - 1)^2}

UsethisinformationandtheDTFTpropertiestofindtheDTFTofthesignalsx1[n],x2[n],x3[n],andx4[n]showninFigs.P9.32b,P9.32c,P9.32d,andP9.32e,respectively.9.33Supposesignal<sup>x</sup>[n]=sinc2(πn/2)modulatesacarriercos(cn)toproducesignaly[n]=x[n]cos(cn).FindandsketchtheDTFTof:(a)x[n](b)y[n]for<sup>c</sup>=π/2(c)y[n]for<sup>c</sup>=3π/4(d)y[n]for<sup>c</sup>=π9.34ShowthatperiodicconvolutionX()Y()=2πX()if Use this information and the DTFT properties to find the DTFT of the signals *x*1[*n*], *x*2[*n*], *x*3[*n*], and *x*4[*n*] shown in Figs. P9.3-2b, P9.3-2c, P9.3-2d, and P9.3-2e, respectively. **9.3-3** Suppose signal *<sup>x</sup>*[*n*] = sinc2(π*n*/2) modulates a carrier cos(c*n*) to produce signal *y*[*n*] = *x*[*n*] cos(c*n*). Find and sketch the DTFT of: (a) *x*[*n*] - (b) *y*[*n*] for <sup>c</sup> = π/2 - (c) *y*[*n*] for <sup>c</sup> = 3π/4 - (d) *y*[*n*] for <sup>c</sup> = π **9.3-4** Show that periodic convolution *X*()-∗ *Y*() = 2π*X*() if

X(\Omega) = \sum_{k=0}^{4} a_k e^{-jk\Omega}

and and

Y(\Omega) = \frac{\sin(5\Omega/2)}{\sin(\Omega/2)} e^{-j2\Omega}

whereakisasetofarbitraryconstants.9.35Usingonlypair2(Table9.1)andpropertiesofDTFT,findtheDTFTofthefollowingsignals,assuminga<1and<sup>0</sup><π.(a)ancos0nu[n](b)<sup>n</sup><sup>2</sup>anu[n](c)(nk)a<sup>2</sup>nu[nm]9.36Usepair10inTable9.1,andsuitablepropertiesoftheDTFT,toderivepairs11,12,13,14,15,and16.9.37Usethetimeshiftingpropertytoshowthat where *ak* is a set of arbitrary constants. - **9.3-5** Using only pair 2 (Table 9.1) and properties of DTFT, find the DTFT of the following signals, assuming |*a*| < 1 and <sup>0</sup> < π. - (a) *an* cos0*nu*[*n*] - (b) *<sup>n</sup>*<sup>2</sup>*anu*[*n*] - (c) (*n*−*k*)*a*<sup>2</sup>*nu*[*n*−*m*] - **9.3-6** Use pair 10 in Table 9.1, and suitable properties of the DTFT, to derive pairs 11, 12, 13, 14, 15, and 16. - **9.3-7** Use the time-shifting property to show that

x[n+k]+x[n-k] \Longleftrightarrow 2X(\Omega)\cos k\Omega

(e)FigureP9.32FigureP9.37UsethisresulttofindtheDTFTofthesignalsshowninFig.P9.37.9.38Usethetimeshiftingpropertytoshowthat (e) **Figure P9.3-2** **Figure P9.3-7** Use this result to find the DTFT of the signals shown in Fig. P9.3-7. **9.3-8** Use the time-shifting property to show that

x[n+k] - x[n-k] \Longleftrightarrow 2jX(\Omega) \sin k\Omega

UsethisresulttofindtheDTFTofthesignalshowninFig.P9.38.9.39Supposesignalx[n]hasspectrumX()thatisbandlimitedtoπ/2rad/sample.Next,definesignaly[n]as Use this result to find the DTFT of the signal shown in Fig. P9.3-8. **9.3-9** Suppose signal *x*[*n*] has spectrum *X*() that is bandlimited to π/2 rad/sample. Next, define signal *y*[*n*] as

y[n] = \begin{cases} x[n] & n \text{ even} \ 0 & n \text{ odd} \end{cases}

DeterminethespectrumofY()intermsofX().SketchY()if,overππ, Determine the spectrum of *Y*() in terms of *X*(). Sketch *Y*() if, over −π ≤ ≤ π,

Y(\Omega) = \begin{cases} |2\Omega/\pi| & -\pi/2 \le \Omega \le \pi/2\ 0 & \text{otherwise} \end{cases}

9.310RepeatProb.9.39ify[n]isinsteaddefinedas **9.3-10** Repeat Prob. 9.3-9 if *y*[*n*] is instead defined as

y[n] = \begin{cases} x[n] & n \text{ odd} \ 0 & n \text{ even} \end{cases}

9.311Usingonlypair2inTable9.1andtheconvolutionproperty,findtheinverseDTFTofX()=<sup>e</sup>2<sup>j</sup>/(ej<sup></sup>γ)2.9.312InTable9.1,youaregivenpair1.FromthisinformationandusingsuitablepropertiesoftheDTFT,derivepairs2,3,4,5,6,and7ofTable9.1.Forexample,startingwithpair1,derivepair2.Frompair2,usesuitablepropertiesoftheDTFTtoderivepair3.Frompairs2and3,derivepair4,andsoon.9.313Fromthepairej(0/2)<sup>n</sup><sup>2</sup>πδ(<sup></sup>(0/2))overthefundamentalband,andthefrequencyconvolutionproperty,findtheDTFTofej0n.Assume<sup>0</sup><π/2.9.314FromthedefinitionandpropertiesoftheDTFT,showthat(a) - **9.3-11** Using only pair 2 in Table 9.1 and the convolution property, find the inverse DTFT of *X*() = *<sup>e</sup>*2*<sup>j</sup>*/(*ej* <sup>−</sup>γ )2. - **9.3-12** In Table 9.1, you are given pair 1. From this information and using suitable properties of the DTFT, derive pairs 2, 3, 4, 5, 6, and 7 of Table 9.1. For example, starting with pair 1, derive pair 2. From pair 2, use suitable properties of the DTFT to derive pair 3. From pairs 2 and 3, derive pair 4, and so on. - **9.3-13** From the pair *ej*(0/2)*<sup>n</sup>* ⇐⇒ <sup>2</sup>πδ( <sup>−</sup> (0/2)) over the fundamental band, and the frequency-convolution property, find the DTFT of *ej*0*n*. Assume <sup>0</sup> <π/2. - **9.3-14** From the definition and properties of the DTFT, show that (a)

\sum_{n=-\infty}^{\infty} \operatorname{sinc}(\Omega_c n) = \frac{\pi}{\Omega_c} \quad \Omega_c < \pi

\n(b) \n(b)

\sum_{n=-\infty}^{\infty} (-1)^n \operatorname{sinc}(\Omega_c n) = 0 \quad \Omega_c < \pi

\n(c) \n(c)

\sum_{n=-\infty}^{\infty} \operatorname{sinc}^2(\Omega_c n) = \frac{\pi}{\Omega_c} \quad \Omega_c < \pi/2

\n(d) \n(d)

\sum_{n=-\infty}^{\infty} (-1)^n \operatorname{sinc}^2(\Omega_c n) = 0 \quad \Omega_c < \pi/2

\n(e) \n(e)

\int_{-\pi}^{\pi} \frac{\sin(M\Omega/2)}{\sin(\Omega/2)} = 2\pi \quad \text{odd } M

\int_{-\pi}^{\infty} \frac{\sin(\alpha z/2)}{\sin(\alpha z/2)} \sin(\alpha z/2)

(f) (f)

\sum_{n=-\infty}^{\infty} |\sin(\alpha z/2)|^4 = 2\pi/3\Omega_c \quad \Omega_c < \pi/2

9.315Showthattheenergyofsignalxc(t)specifiedinEq.(9.41)isidenticaltoTtimestheenergyofthediscretetimesignalx[n],assumingxc(t)isbandlimitedtoB1/2THz.[Hint:Recallthat **9.3-15** Show that the energy of signal *xc*(*t*) specified in Eq. (9.41) is identical to *T* times the energy of the discrete-time signal *x*[*n*], assuming *xc*(*t*) is bandlimited to *B* ≤ 1/2*T* Hz. [*Hint:* Recall that

\int_{-\infty}^{\infty} \operatorname{sinc} [\pi(t-m)] \operatorname{sinc} [\pi(t-n)] dt

= \begin{cases} 0 & m \neq n \ 1 & m = n \end{cases}

Thatis,sincfunctionsareorthogonal.]9.41UsetheDTFTmethodtofindthezerostateresponsey[n]ofacausalsystemwithfrequencyresponse That is, sinc functions are orthogonal.] **9.4-1** Use the DTFT method to find the zero-state response *y*[*n*] of a causal system with frequency response

H(\Omega) = \frac{e^{i\Omega} + 0.32}{e^{i2\Omega} + e^{i\Omega} + 0.16}

andtheinput<sup>x</sup>[n]=(0.5)nu[n].9.42RepeatProb.9.41for and the input *<sup>x</sup>*[*n*] = (−0.5)*nu*[*n*]. **9.4-2** Repeat Prob. 9.4-1 for

H(\Omega) = \frac{e^{i\Omega} + 0.32}{e^{i2\Omega} + e^{i\Omega} + 0.16}

andinputx[n]=u[n].9.43RepeatProb.9.41for and input *x*[*n*] = *u*[*n*]. **9.4-3** Repeat Prob. 9.4-1 for

H(\Omega) = \frac{e^{i\Omega}}{e^{i\Omega} - 0.5}

and and

x[n] = 0.8nu[n] + 2(2)nu[-(n+1)]

9.44DetermineandsketchthemagnitudeandphaseresponseforanLTIDsystemspecifiedbytheequation **9.4-4** Determine and sketch the magnitude and phase response for an LTID system specified by the equation

y[n] + 0.5y[n-1] = x[n] - 0.9x[n-1]

Determinethesystemoutputy[n]fortheinput<sup>x</sup>[n]=cos(<sup>π</sup><sup>n</sup><sup>3</sup>+0.5).9.45RepeatProb.9.44iftheLTIDsystemisinsteadspecifiedbytheequation Determine the system output *y*[*n*] for the input *<sup>x</sup>*[*n*] = cos( <sup>π</sup>*<sup>n</sup>* <sup>3</sup> +0.5). **9.4-5** Repeat Prob. 9.4-4 if the LTID system is instead specified by the equation

y[n] - 0.5y[n-1] = x[n] + 0.9x[n-1]

9.46Anaccumulatorsystemhasthepropertythataninputx[n]resultsintheoutput **9.4-6** An accumulator system has the property that an input *x*[*n*] results in the output

y[n] = \sum_{k=-\infty}^{n} x[k]

(a)Findtheunitimpulseresponseh[n]andthefrequencyresponseH()fortheaccumulator.(b)Usetheresultsofpart(a)tofindtheDTFTofu[n].9.47Anoncausal7pointmovingaverageisdescribedbytheequation - (a) Find the unit impulse response *h*[*n*] and the frequency response *H*() for the accumulator. - (b) Use the results of part (a) to find the DTFT of *u*[*n*]. - **9.4-7** A noncausal 7-point moving average is described by the equation

y[n] = \frac{1}{7} \sum_{k=-3}^{3} x[n-k]

(a)Findandsketchthemagnitudeandphaseresponsesofthesystem.(b)Howcanthissystembemadecausal?Plotthemagnitudeandphaseresponsesofthecausalsystem,andcommentonanydifferencesfrompart(a).9.48AnLTIDsystemfrequencyresponseoverπis - (a) Find and sketch the magnitude and phase responses of the system. - (b) How can this system be made causal? Plot the magnitude and phase responses of the causal system, and comment on any differences from part (a). - **9.4-8** An LTID system frequency response over || ≤ π is

H(\Omega) = \text{rect}\bigg(\frac{\Omega}{\pi}\bigg)e^{-j2\Omega}

Findtheoutputy[n]ofthissystem,iftheinputx[n]isgivenby(a)sinc(πn/2)(b)sinc(πn)(c)sinc2(πn/4)9.49(a)Ifx[n]X(),then,showthat(1)nx[n]<sup>X</sup>(π).(b)Sketch<sup>γ</sup>nu[n]and(γ)nu[n]for<sup>γ</sup><sup>=</sup>0.8;seethespectrafor<sup>γ</sup>nu[n]inFigs.9.5band9.5c.Fromthesespectra,sketchthespectrafor(γ)nu[n].(c)Anideallowpassfilterofcutofffrequency<sup>c</sup>isspecifiedbythefrequencyresponseH()=rect(/2c).Finditsimpulseresponseh[n].Findthefrequencyresponseofafilterwhoseimpulseresponseis(1)nh[n].Sketchthefrequencyresponseofthisfilter.Whatkindoffilteristhis?9.410Ananalogdifferentiator<sup>y</sup>(t)<sup>=</sup><sup>d</sup>dtx(t)canbeapproximatedusingabackwarddifferencesystemdescribedas Find the output *y*[*n*] of this system, if the input *x*[*n*] is given by - (a) sinc (π*n*/2) - (b) sinc(π*n*) - (c) sinc2 (π*n*/4) - **9.4-9** (a) If *x*[*n*] ⇐⇒ *X*(), then, show that (−1)*nx*[*n*] ⇐⇒ *<sup>X</sup>*(−π ). - (b) Sketch <sup>γ</sup> *nu*[*n*] and (−γ )*nu*[*n*] for <sup>γ</sup> <sup>=</sup> 0.8; see the spectra for <sup>γ</sup> *nu*[*n*] in Figs. 9.5b and 9.5c. From these spectra, sketch the spectra for (−γ )*nu*[*n*]. - (c) An ideal lowpass filter of cutoff frequency *<sup>c</sup>* is specified by the frequency response *H*() = rect(/2*c*). Find its impulse response *h*[*n*]. Find the frequency response of a filter whose impulse response is (−1)*nh*[*n*]. Sketch the frequency response of this filter. What kind of filter is this? - **9.4-10** An analog differentiator *<sup>y</sup>*(*t*) <sup>=</sup> *<sup>d</sup> dt x*(*t*) can be approximated using a backward difference system described as

y[n] = \frac{x[n] - x[n-1]}{T}

FindandsketchthemagnitudeandphaseresponsesofthisDTsystem.Forwhatfrequenciesdoesthesystemmostbehaveasadifferentiator?Forwhatfrequenciesdoesthesystemleastbehaveasadifferentiator?9.411Afilterwithimpulseresponseh[n]ismodifiedasshowninFig.P9.411.Determinetheresultingfilterimpulseresponseh1[n].FindalsotheresultingfilterfrequencyresponseH1()intermsofthefrequencyresponseH().HowareH()andH1()related?9.412(a)ConsideranLTIDsystemS1,specifiedbyadifferenceequationoftheformofEqs.(3.15)or(3.16)or(3.20)inCh.3.WeconstructanothersystemS<sup>2</sup>byreplacingcoefficientsai(i=0,1,2,...,N)bycoefficients(1)<sup>i</sup>aiandreplacingallcoefficientsbi(i=0,1,2,...,N)withcoefficients(1)<sup>i</sup>bi.Howarethefrequencyresponsesofthetwosystemsrelated?(b)IfS<sup>1</sup>representsalowpassfilter,whatkindoffilterisspecifiedbyS2?(c)Whattypeoffilter(lowpass,highpass,etc.)isspecifiedbythedifferenceequation Find and sketch the magnitude and phase responses of this DT system. For what frequencies does the system most behave as a differentiator? For what frequencies does the system least behave as a differentiator? **9.4-11** A filter with impulse response *h*[*n*] is modified as shown in Fig. P9.4-11. Determine the resulting filter impulse response *h*1[*n*]. Find also the resulting filter frequency response *H*1() in terms of the frequency response *H*(). How are *H*() and *H*1() related? - **9.4-12** (a) Consider an LTID system *S*1, specified by a difference equation of the form of Eqs. (3.15) or (3.16) or (3.20) in Ch. 3. We construct another system *S*<sup>2</sup> by replacing coefficients *ai* (*i*=0, 1, 2,...,*N*) by coefficients (−1)*<sup>i</sup> ai* and replacing all coefficients *bi* (*i* = 0, 1, 2,...,*N*) with coefficients (−1)*<sup>i</sup> bi*. How are the frequency responses of the two systems related? - (b) If *S*<sup>1</sup> represents a lowpass filter, what kind of filter is specified by *S*2? - (c) What type of filter (lowpass, highpass, etc.) is specified by the difference equation

y[n] - 0.8y[n-1] = x[n]

Whatkindoffilterisspecifiedbythefollowingdifferenceequation? What kind of filter is specified by the following difference equation?

y[n] + 0.8y[n-1] = x[n]

9.413(a)ThesystemshowninFig.P9.413containstwoidenticalLTIDfilterswithfrequencyresponseH0()andcorrespondingimpulseresponseh0[n].Itiseasytoseethatthesystemislinear.Showthatthissystemisalsotimeinvariant.DothisbyfindingtheFigureP9.411responseofthesystemtoinputδ[nk]intermsofh0[n].(b)IfH0()=rect(/2W)overthefundamentalband,and<sup>c</sup>+Wπ,findH(),thefrequencyresponseofthissystem.Whatkindoffilteristhis?9.51DeterminetheDTFTofx[n]=sin(0n)fromtheCTFTofxc(t)=sin(ω0t).9.52ACTsignalx(t),bandlimitedto25kHz,issampledat50kHztoproduce **9.4-13** (a) The system shown in Fig. P9.4-13 contains two identical LTID filters with frequency response *H*0() and corresponding impulse response *h*0[*n*]. It is easy to see that the system is linear. Show that this system is also time-invariant. Do this by finding the **Figure P9.4-11** response of the system to input δ[*n* − *k*] in terms of *h*0[*n*]. - (b) If *H*0() = rect(/2*W*) over the fundamental band, and *<sup>c</sup>* + *W* ≤ π, find *H*(), the frequency response of this system. What kind of filter is this? - **9.5-1** Determine the DTFT of *x*[*n*] = sin(0*n*) from the CTFT of *x*c(*t*) = sin(ω0*t*). - **9.5-2** A CT signal *x*(*t*), bandlimited to 25 kHz, is sampled at 50 kHz to produce

x[n] = \delta[n+4] - 2\delta[n+2] + \delta[n+1] - 3\delta[n] - \delta[n-1] - 2\delta[n-2] - \delta[n-4]

DeterminetheCTFTX(ω).9.71ThisproblemusesamatrixbasedapproachtoinvestigatethecomputationoftheinverseDTFS.(a)ImplementEq.(9.3),theinverseDTFS,usingamatrixbasedapproach.(b)ComparetheexecutionspeedofthematrixbasedapproachtotheIFFTbasedapproachforinputvectorsofsizes10,100,and1000.FigureP9.413(c)WhatistheresultofmultiplyingtheDFTmatrixWN<sup>0</sup>bytheinverseDTFSmatrix?Discussyourresult.9.72Astable,firstorderhighpassIIRdigitalfilterhastransferfunction Determine the CTFT *X*(ω). - **9.7-1** This problem uses a matrix-based approach to investigate the computation of the inverse DTFS. - (a) Implement Eq. (9.3), the inverse DTFS, using a matrix-based approach. - (b) Compare the execution speed of the matrix-based approach to the IFFT-based approach for input vectors of sizes 10, 100, and 1000. **Figure P9.4-13** - (c) What is the result of multiplying the DFT matrix **W***N*<sup>0</sup> by the inverse DTFS matrix? Discuss your result. - **9.7-2** A stable, first-order highpass IIR digital filter has transfer function

H[z] = \left(\frac{1+\alpha}{2}\right) \left(\frac{1-z^{-1}}{1-\alpha z^{-1}}\right)

(a)Deriveanexpressionrelatingαtothe3dBcutofffrequencyc.(b)Testyourexpressionfrompart(a)inthefollowingmanner.First,computeαtoachievea3dBcutofffrequencyof1kHz,assumingasamplingrateofF<sup>s</sup>=5kHz.Determineadifferenceequationdescriptionofthesystem,andverifythatthesystemisstable.Next,computeandplotthemagnituderesponseoftheresultingfilter.Verifythatthefilterishighpassandhasthecorrectcutofffrequency.(c)Holdingαconstant,whathappenstothecutofffrequency<sup>c</sup>asF<sup>s</sup>isincreasedto50kHz?WhathappenstothecutofffrequencyfcasF<sup>s</sup>isincreasedto50kHz?(d)IsthereawellbehavedinversefiltertoH[z]?Explain.(e)Determineαfor<sup>c</sup>=π/2.Commentontheresultingfilter,particularlyh[n].9.73Usingthefrequencysamplingmethod,designalength35linearphaseFIRhighstopfilterthathascutofffrequency<sup>c</sup>=2π/3.Plottheresultingfiltersimpulseresponseh[n]andmagnituderesponseH().9.74Usingthefrequencysamplingmethod,designalength71linearphaseFIRbandstopfilterthathasstopband(π/3<<π/2).Plottheresultingfiltersimpulseresponseh[n]andmagnituderesponseH().9.75FigureP9.75providesthedesiredmagnituderesponseH()ofarealfilter.Mathematically, - (a) Derive an expression relating α to the 3 dB cutoff frequency *c*. - (b) Test your expression from part (a) in the following manner. First, compute α to achieve a 3 dB cutoff frequency of 1 kHz, assuming a sampling rate of *F<sup>s</sup>* = 5 kHz. Determine a difference equation description of the system, and verify that the system is stable. Next, compute and plot the magnitude response of the resulting filter. Verify that the filter is highpass and has the correct cutoff frequency. - (c) Holding α constant, what happens to the cutoff frequency *<sup>c</sup>* as *F<sup>s</sup>* is increased to 50 kHz? What happens to the cutoff frequency *fc* as *F<sup>s</sup>* is increased to 50 kHz? - (d) Is there a well-behaved inverse filter to *H*[*z*]? Explain. - (e) Determine α for *<sup>c</sup>* = π/2. Comment on the resulting filter, particularly *h*[*n*]. - **9.7-3** Using the frequency sampling method, design a length-35 linear phase FIR highstop filter that has cutoff frequency <sup>c</sup> = 2π/3. Plot the resulting filter's impulse response *h*[*n*] and magnitude response |*H*()|. - **9.7-4** Using the frequency-sampling method, design a length-71 linear phase FIR bandstop filter that has stopband (π/3 < || < π/2). Plot the resulting filter's impulse response *h*[*n*] and magnitude response |*H*()|. - **9.7-5** Figure P9.7-5 provides the desired magnitude response |*H*()| of a real filter. Mathematically,

|H(\Omega)| = \begin{cases} 2\frac{4\Omega}{\pi} & 0 \leq \Omega < \frac{\pi}{4} \ 2 - \frac{4\Omega}{\pi} & \frac{\pi}{4} \leq \Omega < \frac{\pi}{2} \ 0 & \frac{\pi}{2} \leq \Omega \leq \pi \end{cases}

Sincethedigitalfilterisreal,H()=H()andH()=H(+2π)forall.(a)Canarealizablefilterhavethisexactmagnituderesponse?Explainyouranswer.(b)UsethefrequencysamplingmethodtodesignanFIRfilterwiththismagnituderesponse(orareasonableapproximation).UseMATLABtoplotthemagnituderesponseofyourfilter.9.76ArealFIRcombfilterisneededthathasmagnituderesponseH()=[0,3,0,3,0,3,0,3]for=[0,π/4,π/2,3π/4,π,5π/4,3π/2,7π/4],respectively.Providetheimpulseresponseh[n]ofafilterthataccomplishesthesespecifications.9.77ApermutationmatrixPhasasingleoneineachrowandcolumnwiththeremainingelementsallzero.Permutationmatricesareusefulforreorderingtheelementsofavector;theoperationPxreorderstheelementsofacolumnvectorxbasedontheformofP.(a)FullydescribeanN<sup>0</sup>×N<sup>0</sup>permutationmatrixnamedRN<sup>0</sup>thatreversestheorderoftheelementsofacolumnvectorx.(b)GivenDFTmatrixWN<sup>0</sup>,verifythat(WN<sup>0</sup>)(WN<sup>0</sup>)<sup>=</sup><sup>W</sup><sup>2</sup><sup>N</sup><sup>0</sup>producesascaledpermutationmatrix.HowdoesW<sup>2</sup>N0xreordertheelementsofx?(c)Whatistheresultof(W<sup>2</sup>N0)(W<sup>2</sup>N0)x=W<sup>4</sup>N0x? Since the digital filter is real, |*H*()|=|*H*(−)| and |*H*()|=|*H*(+2π )| for all . - (a) Can a realizable filter have this exact magnitude response? Explain your answer. - (b) Use the frequency-sampling method to design an FIR filter with this magnitude response (or a reasonable approximation). Use MATLAB to plot the magnitude response of your filter. - **9.7-6** A real FIR comb filter is needed that has magnitude response |*H*()|=[0, 3, 0, 3, 0, 3, 0, 3] for = [0,π/4,π/2, 3π/4,π, 5π/4, 3π/2, 7π/4], respectively. Provide the impulse response *h*[*n*] of a filter that accomplishes these specifications. - **9.7-7** A permutation matrix **P** has a single one in each row and column with the remaining elements all zero. Permutation matrices are useful for reordering the elements of a vector; the operation **Px** reorders the elements of a column vector **x** based on the form of **P**. - (a) Fully describe an *N*<sup>0</sup> × *N*<sup>0</sup> permutation matrix named **R***N*<sup>0</sup> that reverses the order of the elements of a column vector **x**. - (b) Given DFT matrix **W***N*<sup>0</sup> , verify that (**W***N*<sup>0</sup> )(**W***N*<sup>0</sup> ) <sup>=</sup> **<sup>W</sup>**<sup>2</sup> *<sup>N</sup>*<sup>0</sup> produces a scaled permutation matrix. How does **W**<sup>2</sup> *N*0 **x** reorder the elements of **x**? - (c) What is the result of (**W**<sup>2</sup> *N*0 )(**W**<sup>2</sup> *N*0 )**x** = **W**<sup>4</sup> *N*0 **x**?