ELEC3114 Sample Paper Solutions Matthew Davis [email protected] 20th December 2013 Abstract These solutions are for the sample nal exam given out in S2 2013 for ELEC3114. These solutions are not aliated with the Electrical Engineering and Telecommunications, they are meant as a guide only. Other solutions and notes can be found at http://www.elsoc.net/notes.php. Question 1 This paper is the one which starts with G (s) = ( ) 10 s2 + 0.2s + 1 s (s2 + 1.2s + 9) and C (s) = a) s + 0.01 s (s + 0.1) ( ) 10 s2 + 0.2s + 1 (s + 0.01) H (s) = C (s) × G (s) = 2 2 s (s + 1.2s + 9) (s + 0.1) The open loop transfer function has a factor of s2 in the denominator, therefore it is type 2. b) The steady state error for r (t) = 1 and r (t) = 0 will be zero, since this is a type 2 system. Noting that the standard parabola input is 21 t2 (not t2 ) and the question asks for t2 , that means 1 e (∞) = because ka = 10×1×0.01 9×0.1 c) 2 2 = = 18 ka 0.111 = 0.111 · · · ( 2 )( s ) s + 0.2s + 1 0.01 +1 ( s ) H (s) = 2 2 s (s + 1.2s + 9) 0.1 +1 √ −0.2± 0.2 −4 The zeros of the transfer function of 0.01 rad s−1 and √ = 2 2 2 −0.1 ± 0.99j , which corresponds to a double break frequency of 0.1 + 0.99 ≈ 1rad s−1 . √ 2 −4×9 The poles of the transfer function are 0.1rad s−1 and −1.2± 1.2 = 2 −0.6 ± 2.939j 2 √ Which corresponds to a double break frequency at 0.62 + 2.9392 = 3rad s−1 Now putting these in ascending order we have • 0.01rad s−1 (root) • 0.1rad s−1 (pole) • 1rad s−1 (double zero) • 3rad s−1 (double root) d) I don't actually know what this question is asking. Does it want the y = mx + b equations for each asymptote? That seems like way too much eort. e) ±20 log √ 2 ± 3dB f) For s2 + 0.2s + 1 0.2 ωn = 1, ζ = 2×1 = 0.1 Therefore the bump corresponding to these break points is at MP = 1 1 √ √ = = 5.025 = 14dB 2 s × 0.1 × 1 − 0.12 2ζ 1 − ζ below the corresponding asymptote. √ For s2 + 1.2s + 9 we can work out ωn = 9 = 3, ζ = the bump corresponding to these break points is at 2 1.2 2×3 = 0.2 Therefore MP = 1 1 √ √ = = 2.55 = 8.14dB 2 2 × 0.2 × 1 − 0.22 2ζ 1 − ζ above the corresponding asymptote. g) 3 You can check this on Matlab with bode ( t f ( [ 1 0 . 0 1 ] , [ 1 0 . 1 ] ) ∗ t f ( 1 0 ∗ [ 1 0 . 2 1 ] , [ 1 1 . 2 9 0 0 ] ) ) The corrections due to resonance are the reason the Matlab output looks dierent in the phase plot. Question 2 a) Yes it is stable. The poles are all in the LHP, except one at the origin (but that's ok because it's only rst order). 4 5 b) From a quick bode plot, or by solving G (jω) ∈ ℜ we see ω = ±1 c) See above. d) The intercept on the negative real axis in the Nyquist plot was at -1.58. There1 fore the system is stable for K < 1.5 = 0.67 e) The gain margin for K=1 is 4dB (from the bode plot). K = margin is − (4 − 9.5) = 5.5dB f) The characteristic polynomial is s3 + 2s2 + s + 3K s3 s2 s1 s0 2−3K 2 1 2 =⇒ 2 − 3K 3K So we want 2 − 3K > 0 6 1 3K 0 =⇒ 0 0 1 3 so the gain 0 < K < 0.67 Which is consistent with part d. Question 3 Credit goes to Shannon `British racing' Green for the solutions to this question. a) ωn = sp = −ωn ζ ± ωn √ π √ = TP 1 − ζ 2 π 3 π √ = 3.015 1 − 0.12 1 − ζ 2 j = −0.3015 × 0.1 ± 3.015 × √ 1 − 0.12 j = −0.3015 ± 3j b) By plotting the closed and open loop poles and zeros on the complex plane, then summing the angles, we can nd z1 My advice is to buy a Casio fx_100au. It will make this so much simpler. ) ( ) ( ) 3 3 3 −1 −1 (2k + 1)×180 = tan −2 tan −tan z1 − 0.03015 −0.03015 2 − 0.03015 ( ) ( ) 3 3 −2 tan + tan 8 − 0.03015 10 − 0.03015 ◦ −1 ( ( ) ( ) ( ) 3 3 3 −1 −1 (2k + 1)×180 = tan +2 tan −tan z1 − 0.03015 0.03015 2 − 0.03015 ( ) ( ) 3 3 −2 tan + tan−1 8 − 0.03015 10 − 0.03015 ( ) 3 ◦ −1 (2k + 1) × 180 = tan + 103.9◦ z1 − 0.03015 ( ) 3 (2k + 1) × 180◦ − 103.9◦ = tan−1 z1 − 0.03015 ◦ −1 z1 = 1.0449 This solution means z1 < p1 , that's a good indication that we haven't gone horrible wrong. (−0.3+3j)2 (−0.3+3j+0.1)((−0.3+3j)2 +1.2(−0.3+3j)+9) 1 = 8.23 = K= |Clead (−0.3 + 3k) × G (−0.3 + 3j) | 10×((−0.3+3j)2 +0.2(−0.3+3j)+1)(−0.3+3j+0.01) 7 c) The root locus is only on the real axis when there is an odd number of roots and poles to the right of that interval. Therefore [−10, −8] and [−2, z1 ] lie on that interval. d) At s = 0 the angle contributions from all open loop poles and zeros not at the origin is 0◦ . ◦ = ±90◦ Therefore the departure angle is ±180 2 e) θa = σa = 180◦ = 90◦ 4−2 −8 − 2 − 0 − 0 − (−10 − 1.05) = 0.53 4−2 f) We can work out the solution is K=114 in Matlab using r l o c u s ( zpk ([ − 10 − 1 . 0 5 ] , [ 0 0 −8 − 2] ,1)) I can't work out the solution by hand. 8 g) h) Ka = 20 = lim s2 × G (s) Clag (s) Clead (s) = s→0 z2 = 8.23 × 10 × 1.05 × z2 8 × 2 × 0.01 20 × 8 × 2 × 0.01 = 0.037 8.23 × 10 × 1.05 Therefore Clag (s) = s + 0.037 s + 0.01 i) Clag (s) Clead (s) G (s) = K (s + 0.037) (s + 1.05) (s + 2) (s + 0.01) 0 = (s + 2) (s + 0.01) + K (s + 0.036) (s + 1.05) 0 = (1 + K) s2 + (2.01 + 1.086K) s + (0.02 + 0.0378K) Putting this into a Routh table we get: 9 s2 s1 s0 1+K 2.01+1.086K 0.02+0.0378K 0.02+0.0378K 0 0 For K>0 the system is stable. Question 4 a) [ ] 0 −4 A= −1 0 [ ] 0 B= 1 [ ] C= 1 2 D=0 Using the formula provided, −1 T (s) = C × (sI − A) [ T (s) = 1 ( [ ] [ ] 1 0 0 2 × s − 0 1 −1 ×B+D ])−1 [ ] 2 2s − 4 −4 0 = × = 2 0 1 s −4 s+2 b) 2 s+2 T (s) = 1 0 2 0 No sign changes in the rst column mean that this is stable for all K. s1 s0 c) [ C ] [ 1 2 ] The observer matrix is = CA −2 −4 The determinant of this is 0 so the system is not observable. I don't know where to go from here. 10 d) Using the 3 poles given to form the desired transfer function we have (s + 3 − 4j) (s + 3 + 4j) (s + 20) = 0 s3 + 26s2 + 145s + 500 = 0 (1) Then we use the matrix [ ] A − BK BKe −C 0 ] [ ]) [ ] ([ 0 −4 0 0 0 − K] 1 K2 Ke M = −1 0 [ − 1 2 0 0 −4 0 = − (1 + K1 ) −K2 Ke −1 −2 0 Therefore s sI − M = 1 + K1 1 4 s + K2 2 0 −Ke s Then det (sI − M) = s3 + K2 s2 + (2Ke − 4K1 − 4) s − 4Ke = 0 Matching coecients of equations 1 and 2 we get K1 = −99.75 K2 = 26.3 Ke = −125 11 (2)
© Copyright 2024