Solutions. Number of Problems: 10

Transcription

1 Final Exam February 9th, 2 Signals & Systems (5-575-) Prof. R. D Andrea Solutions Exam Duration: 5 minutes Number of Problems: Permitted aids: One double-sided A4 sheet. Questions can be answered in English or German. Use only the prepared sheets for your solutions. Additional paper is available from the supervisors.

2 Problem A linear, time invariant system with input x[n] and output y[n] has the transfer function H(z) = Y(z) X(z) = 2z a 2 +2az +z 2 a) For what real values of a is the system causal and stable? (3 points) b) Write down a causal difference equation of the above system of the following form: ( y[n] = M ) N b k x[n k] a k y[n k] a k= k= (2 points) Solution a) The poles of the system are at = a 2 +2az +z 2 = (a+z ) 2 z = a z = a For the system to be causal and stable, the poles must lie within the unit circle. It follows that a > must hold. b) Y(z) X(z) = 2z a 2 +2az +z 2 Y(z)(a 2 +2az +z 2 ) = X(z)( 2z ) a 2 y[n]+2ay[n ]+y[n 2] = x[n] 2x[n ] y[n] = a 2(x[n] 2x[n ] 2ay[n ] y[n 2])

3 Problem 2 A system is governed by the difference equation ( y[n] = x[n]x[n ]+cos 3πn π ), 3 where x[n] is the input and y[n] is the output. Is the system a) linear? Yes No b) time invariant? Yes No c) bounded input bounded output (BIBO) stable? Yes No d) memoryless? Yes No e) causal? Yes No Please circle either Yes or No. If you change your mind, please cross out both Yes and No and write either Yes or No alongside, or leave it blank. It is not necessary to justify choices. You can get a maximum of and a minimum of points for this problem. For each subproblem you get: + for a correct answer, - for an incorrect answer and for no answer. Solution 2 a) No b) No c) Yes d) No e) Yes

4 Problem 3 The output of a system T is y[n] = T{x[n]}. T is a causal, linear, time invariant system. Given the input and output below, calculate the impulse response h[n] of the system T for all n. n < > 5 x[n] 2 y[n] Solution 3 Because the system is given to be causal, h[n] = for all n <. One can see from the input output behavior that the system has a finite impulse response (the output goes to zero a finite number of steps after the input goes to zero). The input is zero for n > 2, and the output is zero for n > 5; therefore the impulse response must be zero for n > 3. n = : n = : n = 2: y[] = h[]x[] h[] = y[] x[] = 3 y[] = h[]x[] + h[]x[] h[] = y[] h[]x[] x[] = y[2] = h[]x[2] + h[]x[] + h[2]x[] h[2] = y[2] h[]x[2] h[]x[] x[] = 5

5 n = 3: y[3] = h[]x[3]+h[]x[2]+h[2]x[]+h[3]x[] h[3] = y[3] h[]x[3] h[]x[2] h[2]x[] x[] = 2

6 Problem 4 The transfer function of a filter is H(z) = Y(z) X(z) = b +a z. Calculate the coefficients b and a such that the filter is stable and causal, and such that the frequency response H(Ω) of the filter fulfills the two criteria H(Ω = ) =, and ( H Ω = π ) = 2. 2 Solution 4 The first criterion yields = b +a e j = b +a, +a = b. From the second criterion, we obtain = b 2 ja = b. +a 2 Combining the two criteria, we obtain = +a. 2 +a 2 It follows that = (+a ) 2 +a 2 2. = a 2 +4a +.

7 Solving the quadratic equation for a, we obtain the two solutions a = 2± 3. Since the filter needs to be stable, causal, and has a pole at z = a, the only possible solution is a = 2+ 3, from which follows b = + 3.

8 Problem 5 Associate the following impulse responses ((a) to (e)) and frequency responses (() to (5)) with the corresponding difference equation by filling out the following table. It is not necessary to justify choices. Difference equation y[n] = 3 (x[n]+x[n ]+x[n 2]) y[n] = 2 (x[n] x[n ]) y[n] =.8x[n]+.2y[n ] y[n] = 3 (x[n+]+x[n]+x[n ]) y[n] =.2x[n]+.8y[n ] Impulse Response Frequency Response (a) to (e) () to (5) Impulse response h[n] Time index n Impulse response h[n] Time index n (a) (b) Impulse response h[n] Time index n Impulse response h[n] Time index n (c) (d) Impulse response h[n] Time index n (e)

9 Amplitude Amplitude Phase ( ) 9 9 Phase ( ) Frequency Ω Frequency Ω () (2) Amplitude Amplitude Phase ( ) 9 9 Phase ( ) Frequency Ω Frequency Ω (3) (4) Amplitude Phase ( ) Frequency Ω (5)

10 Solution 5 Difference equation Impulse Response Frequency Response (a) to (e) () to (5) y[n] = 3 (x[n]+x[n ]+x[n 2]) d y[n] = 2 (x[n] x[n ]) c 5 y[n] =.8x[n]+.2y[n ] e 4 y[n] = 3 (x[n+]+x[n]+x[n ]) b 3 y[n] =.2x[n]+.8y[n ] a 2

11 Problem 6 For a zero mean signal x[n], R xx [] = σ 2 x, where R xx is the autocorrelation function and σ 2 x is the variance of the signal x[n]. The autocorrelation function is defined as R xx [k] := E(x[n]x[n k]), where E( ) is the expected value operator. Let x[n] be white noise with varianceσ 2 x = 2. Thesignalx[n]isappliedtoalineartimeinvariantsystem T to generate the output y[n] = T{x[n]}. The frequency response H(Ω) of T can be approximated by { H(Ω) = 2 e jω, Ω Ω c, Ω c < Ω π. Calculate the frequency Ω c such that the variance of the output signal y[n] is σ 2 y =.4. Solution 6 The power spectral density function of y[n] is { S yy (Ω) = H(Ω) 2 S xx (Ω) = 4 σ2 x = 2, Ω Ω c, Ω c < Ω π. The system T is linear, therefore the zero mean input x[n] results in a zero mean output y[n]. We calculate the variance of y[n] using the inverse Fourier transform σy 2 = R yy [k = ] = π S yy (Ω)e jω dω 2π = π Ωc 2 dω = Ω c 2π. Solving for Ω c with σ 2 y =.4, we obtain π Ω c = 4π 5.

12 Problem 7 Consider the continuous time signal ( x(t) = 4cos 6πt+ π ) +2cos(8πt+π) 2. 3 a) Determine the largest possible sampling time T s in seconds to sample the signal without aliasing effects. b) Sample the given continuous time signal with the sampling time you obtained in a), starting at t = : ( point) ( point) x[] = x(), x[] = x(t s ),..., x[n] = x(nt s ). Determine the fundamental period N of the discrete time signal x[n]. c) State the Fourier series coefficients c k given by (3 points) c k = N N n= x[n]e j2πk N n, k =,,...,N. Solution 7 a) The largest possible sampling time is T s = /8s. b) The sampled signal is x[n] = 4cos ( 3 4 πn+ π ) +2cos(πn+π) 2. 3 The fundamental period of the signal is N = 8. c) The coefficients are k c k -2 2e jπ 3-2 2e jπ 3. For a detailed solution, please refer to the solution of problem 4 of the 29 final exam or the recitation notes posted on Nov. 5 on the class website.

13 Problem 8 A continuous time, linear, time invariant system with input x(t) and output y(t) is given by the state space description q(t) = Aq(t) + Bx(t) y(t) = Cq(t) + Dx(t), where the state vector has two elements q(t) = {q (t),q 2 (t)} T and the matrices are [ ] [ ] A = B = C = [ ] D = [.5] Assume the input x(t) is piece-wise constant x(t) = x[k] kt s t < (k +)T s with sampling time T s. Compute A d,b d,c d,d d of a discrete time state space description of the system when the output y[k] is defined as q[k +] = A d q[k]+b d x[k] y[k] = C d q[k]+d d x[k], y[k] := y(t = kt s ). Solution 8 We solve this problem using the matrix exponential (see lecture notes for derivation). We define: [ ] A B M := =. The discrete time matrices A d,b d are then obtained by calculating [ ] Ad B d = e MT s.

14 Since the matrix M is nilpotent, i.e. M 3 =, it is straightforward to calculate the matrix exponential: [ ] Ad B d = I +MT s + 2 M2 Ts, 2 where I is the identity matrix. The matrices are therefore [ ] [ ] Ts A d = B d = 2 T2 s T s. C d = C D d = D

15 Problem 9 You are given the state space description of a system: [ ] [ ] 2 q[n+] = 4 2 q[n]+ x[n] 25 3 y[n] = [ ] q[n]+ [ 4 ] x[n] a) Is the system bounded input bounded output stable? (2 points) b) Is the system controllable? ( point) c) Calculate the step response of the system for n =,,2. ( point) d) Consider the second system ( point) [ ] [ ] q[n+] = 4 2 q[n]+ x[n] y[n] = [ ] q[n]+ [ 3 ] x[n] is the input output behavior of this system identical to the first one? Solution 9 a) The eigenvalues of the matrix A are λ = λ 2 = 2. The system is therefore stable. [ ] 2 7 b) The controllability matrix is [B AB] = 4. The matrix has full 53 row rank and the system is therefore controllable. c) n = : [ ][ ] [ ] 2 q[] = [ ] 2 = y[] = [ ][ ] + [ 4 ] = 4

16 n = : [ ][ ] [ ] 2 q[2] = [ 25 ] = 4 52 y[] = [ ][ ] 2 + [ 4 ] = n = 2: y[2] = [ ][ ] [ 4 ] 52 = 99 4 d) No. An easy way to see this is [ the ] first value of the step response: The state of both systems is q[] =. The output of the first system is therefore y[] = 4, while the output of the second system is y[] = 3.

17 Problem Assume that you are identifying a linear time invariant system T with input x[n] and output y[n]: y[n] = T{x[n]} You applied the input x[n] and measured the output y[n] for samples, i.e. n, as shown below: Input signal x[n] Output signal y[n] Time Index n You want to use this data to identify a model of the form in the frequency domain. Ĥ(z) = b +b z +b 2 z 2 + +b M z M +a z +a 2 z 2 + +a N z N a) Let Y(Ω) be the Discrete Fourier Transform of the measured output y[n]. How many frequency points will the Transform have in the range Ω π? Now assume that you have calculated the system transfer function at each frequency point Ω k. H(Ω k ) = Y(Ω k) X(Ω k ) b) When fitting a model to the system transfer function, what are the highest possible model orders N and M that you can identify using only the given measurement data? ( point) (2 points)

18 c) A colleague has identified a model from the given measurement data. A plot of the frequency response of his estimated model Ĥ(z) is shown below. Is this result plausible, and why? (2 points) Amplitude Phase ( ) Frequency Ω Solution a) There will be. The Discrete Fourier Transform will result in frequency domain points in a 2π interval. This means that the spacing between frequency points will be Ω = 2π. There will be one point at Ω =. There will be 5 more points in the interval, because 5 Ω < π and 5 Ω > π. b) There are pieces of information (the point at Ω = is real, and all points at Ω are complex.). Therefore, N +M + N +M. c) No. The system has a phase of 9 at Ω =, which means that the system has a complex gain. This is clearly not the case, as the inputoutput data is real.