MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.003: Signals and Systems-Fall 2003 PROBLEM SET 11 SOLUTIONS Problem 1(O&W 1029(d)) In this problem we are asked to sketch the magnitude of the Fourier transform associated with the pole-zero diagram, Figure P10.29(d). In order to do so, we need to make some assumptions (a) The ROC includes the unit circle to ensure the existence of the Fourier transform (b)The gain factor is one so that the system we are dealing with has the following form H(2)= (z-z1)(x+1) where z1 is a positive real number whose magnitude is less than 1 To obtain the frequency response of a DT system, we need to evaluate the magnitude and phase of H(a) along the unit circle in z-plane, i.e., 2=eu for 0<w<2. First,we look at the magnitude plot. In general, for a fixed w we can think of H(e3u)l as tot zeros length of a vector connecting i th zero to eju) IH(eju) i=I( length of a vector connecting j th pole to eju If of zeros or poles is zero, then we define the product above to be 1. In are two poles and no zeros. Thus, the above expression can be s implied lo. our case. there HH(e)l where v and U2 are vectors shown in the figure below:
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.003: Signals and Systems—Fall 2003 Problem Set 11 Solutions Problem 1 (O&W 10.29 (d)) In this problem we are asked to sketch the magnitude of the Fourier transform associated with the pole-zero diagram, Figure P10.29 (d). In order to do so, we need to make some assumptions: (a) The ROC includes the unit circle to ensure the existence of the Fourier transform. (b) The gain factor is one so that the system we are dealing with has the following form: H(z) = 1 (z − z1)(z + z1) , where z1 is a positive real number whose magnitude is less than 1. To obtain the frequency response of a DT system, we need to evaluate the magnitude and phase of H(z) along the unit circle in z− plane, i.e., z = e jω for 0 ≤ ω < 2π. First, we look at the magnitude plot. In general, for a fixed ω we can think of |H(e jω )| as |H(e jω )| = Π # of zeros i=1 ( length of a vector connecting i th zero to e jω ) Π # of poles j=1 ( length of a vector connecting j th pole to e jω ) If # of zeros or poles is zero, then we define the product above to be 1. In our case, there are two poles and no zeros. Thus, the above expression can be s implied to: |H(e jω )| = 1 |v1||v2| , where v1 and v2 are vectors shown in the figure below: 1
∠t From the plot, we see that when w=2, the product of lun and lu2l becomes maximum Thus, we expect to have the minimum magnitude at the frequency. Also, when w=0 or T, the product of the lengths or the two vectors becomes minimum; thus the magnitude of H(eJu) becomes maximum. The magnitude plot of H(eu) for 0 <w<T is shown below t Although you are not asked to sketch the phase, here is a briefly outline on how to sketch
<e =m 1 1 −1 −1 ω × ∠v~1 v~1 × ∠v~2 v~2 −z1 z1 e jω From the plot, we see that when ω = π 2 , the product of |v1| and |v2| becomes maximum. Thus, we expect to have the minimum magnitude at the frequency. Also, when ω = 0 or π, the product of the lengths or the two vectors becomes minimum; thus the magnitude of H(e jω ) becomes maximum. The magnitude plot of H(e jω ) for 0 ≤ ω < π is shown below: |H(e jω )| ω π 2 π Although you are not asked to sketch the phase, here is a briefly outline on how to sketch the phase. 2
The phase ZH(e) can be described as: ∠H(e) angle of vector connecting i th zero to e of poles angle of vec ng j th pole to e) Again for our specific case, using the vectors u and u2 defined above we have ∠H(e-)=-(∠v1+∠v2), The phase starts off at 0 when w=0 and decreases to when w=3 symmetric poles, the phase keeps decreasing to -2T at w = T. The phase plot is bele ∠H(e) Problem 2(O&W 10.34) ym]=yn-1]+yn-2]+r{n-1] Taking the z-transform of this equation Y(z)=z-1Y()+x-2Y(2)+x-1Xx(x2) H(z) X 1+√5 1-√5 H(z) has a zero at z=0 and poles at z1 and Since the system is causal, the ROC of H(a)will be outside the circle containing its outermost pole:[=>|a1 The pzmap and ROC are depicted below
The phase ∠H(e jω ) can be described as: ∠H(e jω ) = # X of zeros i=1 ( angle of vector connecting i th zero to e jω ) − # X of poles j=1 ( angle of vector connecting j th pole to e jω ). Again for our specific case, using the vectors v1 and v2 defined above we have ∠H(e jω ) = −(∠v1 + ∠v2), The phase starts off at 0 when ω = 0 and decreases to −π when ω = π 2 . Because of the symmetric poles, the phase keeps decreasing to −2π at ω = π. The phase plot is shown below: ∠H(e jω ) ω π 2 π −π −2π Problem 2 (O&W 10.34) (a) y[n] = y[n − 1] + y[n − 2] + x[n − 1] Taking the z-transform of this equation: Y (z) = z −1Y (z) + z −2Y (z) + z −1X(z) H(z) = Y (z) X(z) = z −1 1 − z −1 − z −2 = z z 2 − z − 1 = z z − 1 + √ 5 2 ! z − 1 − √ 5 2 ! H(z) has a zero at z = 0 and poles at z1 = 1+√ 5 2 and z2 = 1− √ 5 2 . Since the system is causal, the ROC of H(z) will be outside the circle containing its outermost pole: |z| > |z1|. The pzmap and ROC are depicted below: 3
H(z) 1+√5 z-1+ B 1+√5 (2) 21=-4、l B H(2) 25 H
0 × 1+√ 5 2 × 1− √ 5 2 <e =m (b) H(z) = −z −1 z −1 + 1 + √ 5 2 ! z −1 + 1 − √ 5 2 ! = A z −1 + 1 + √ 5 2 + B z −1 + 1 − √ 5 2 A = z −1 + 1 + √ 5 2 ! H(z)| z−1=− 1+√ 5 2 = − 1 + √ 5 2 √ 5 B = z −1 + 1 − √ 5 2 ! H(z)| z−1=− 1− √ 5 2 = 1 − √ 5 2 √ 5 H(z) = − √ 1 5 1 + 2 1+√ 5 z −1 + √ 1 5 1 + 2 1− √ 5 z −1 , |z| > 1 + √ 5 2 4
Taking the inverse z-transform: h/n= √5(1+√3 an (c)The system is unstable, as its ROC does not contain the unit circle. The instability is 山m() term will grow indefinitely as In order to make the system stable, the rOC must contain the unit circle. For stabilit the ROC should be: ==v5<|z<=tv5 The inverse z-transform of H(a) with this ROC is 4间=(125)(-2)时n Problem 3(O&W 10.42 Because we are dealing with a system which has initial conditions, we may want to the unilateral z-transform. From the properties of the unilateral z-transform, we get the following relationships (x)+y{-1 In each part of this problem, the first step is to take the unilateral z-transform of both sides of the difference equation. To find the zero-input response(ZIR), set the input to 0 and solve for Y(z). To find the zero-state response(ZSR), set the initial conditions to zero and solve for Y(a) ym]+3n-1=xm,y-1=1 In/=/1 uny Taking the unilateral z-transform of both sides of the difference equation Solving for Y(a) gives X 1+ ZIR
Taking the inverse z-transform: h[n] = − 1 √ 5 − 2 1 + √ 5 n u[n] + 1 √ 5 − 2 1 − √ 5 n u[n] (c) The system is unstable, as its ROC does not contain the unit circle. The instability is also apparent in h[n], as the − 2 1− √ 5 n term will grow indefinitely as n → ∞. In order to make the system stable, the ROC must contain the unit circle. For stability, the ROC should be: 2 −1− √ 5 < |z| < 2 −1+√ 5 . The inverse z-transform of H(z) with this ROC is: h[n] = − 1 √ 5 − 2 1 + √ 5 n u[n] − 1 √ 5 − 2 1 − √ 5 n u[−n − 1] Problem 3 (O&W 10.42) Because we are dealing with a system which has initial conditions, we may want to use the unilateral z-transform. From the properties of the unilateral z-transform, we get the following relationships: y[n] ←→ Y (z) y[n − 1] ←→ z −1Y (z) + y[−1] In each part of this problem, the first step is to take the unilateral z-transform of both sides of the difference equation. To find the zero-input response (ZIR), set the input to 0 and solve for Y (z). To find the zero-state response (ZSR), set the initial conditions to zero and solve for Y (z). (a) y[n] + 3y[n − 1] = x[n], y[−1] = 1, x[n] = 1 2 n u[n] Taking the unilateral z-transform of both sides of the difference equation, Y (z) + 3z −1Y (z) + 3y[−1] = X(z). Solving for Y (z) gives, Y (z) = X(z) 1 + 3z −1 | {z } ZSR + −3y[−1] 1 + 3z −1 | {z } ZIR . 5