Data Based Design of 3Term Controllers. Data Based Design of 3 Term Controllers p. 1/10

Size: px
Start display at page:

Download "Data Based Design of 3Term Controllers. Data Based Design of 3 Term Controllers p. 1/10"

Transcription

1 Data Based Design of 3Term Controllers Data Based Design of 3 Term Controllers p. 1/10

2 Data Based Design of 3 Term Controllers p. 2/10 History Classical Control - single controller (PID, lead/lag) is designed based on transfer function or frequency response data

3 Data Based Design of 3 Term Controllers p. 2/10 History Classical Control - single controller (PID, lead/lag) is designed based on transfer function or frequency response data Modern/Post Modern Control (H, H 2, l 1 ) - optimal controller of high order based on state space model

4 Data Based Design of 3 Term Controllers p. 2/10 History Classical Control - single controller (PID, lead/lag) is designed based on transfer function or frequency response data Modern/Post Modern Control (H, H 2, l 1 ) - optimal controller of high order based on state space model Fuzzy Neural Control- model free approach but no guarantee of stability or performance

5 Data Based Design of 3 Term Controllers p. 2/10 History Classical Control - single controller (PID, lead/lag) is designed based on transfer function or frequency response data Modern/Post Modern Control (H, H 2, l 1 ) - optimal controller of high order based on state space model Fuzzy Neural Control- model free approach but no guarantee of stability or performance Our Recent Results - complete set of PID controllers achieving stability and performance based on transfer function model

6 Data Based Design of 3 Term Controllers p. 2/10 History Classical Control - single controller (PID, lead/lag) is designed based on transfer function or frequency response data Modern/Post Modern Control (H, H 2, l 1 ) - optimal controller of high order based on state space model Fuzzy Neural Control- model free approach but no guarantee of stability or performance Our Recent Results - complete set of PID controllers achieving stability and performance based on transfer function model Present Paper - extend these results to the case where only data is available

7 Data Based Design of 3 Term Controllers p. 3/10 Drawbacks of Identification States have no physical significance

8 Data Based Design of 3 Term Controllers p. 3/10 Drawbacks of Identification States have no physical significance Parameters in the model have no physical significance

9 Data Based Design of 3 Term Controllers p. 3/10 Drawbacks of Identification States have no physical significance Parameters in the model have no physical significance State variables have no meaningful dimensions or units

10 Data Based Design of 3 Term Controllers p. 4/10 Main Idea: To synthesize the set of all stabilizing PID/First Order controllers from experimental data of the plant rather than from mathematical models.

11 Data Based Design of 3 Term Controllers p. 4/10 Main Idea: To synthesize the set of all stabilizing PID/First Order controllers from experimental data of the plant rather than from mathematical models. To synthesize the set of all PID/First Order controllers satisfying the given performance requirements from the experimental data of the plant.

12 Data Based Design of 3 Term Controllers p. 4/10 Main Idea: To synthesize the set of all stabilizing PID/First Order controllers from experimental data of the plant rather than from mathematical models. To synthesize the set of all PID/First Order controllers satisfying the given performance requirements from the experimental data of the plant. Knowledge of the frequency response magnitude and phase, equivalently, P(jω), ω [0, )

13 Data Based Design of 3 Term Controllers p. 4/10 Main Idea: To synthesize the set of all stabilizing PID/First Order controllers from experimental data of the plant rather than from mathematical models. To synthesize the set of all PID/First Order controllers satisfying the given performance requirements from the experimental data of the plant. Knowledge of the frequency response magnitude and phase, equivalently, P(jω), ω [0, ) Knowledge of the number of RHP poles p +.

14 Data Based Design of 3 Term Controllers p. 5/10 Assumption The plant is stabilizable (i.e., N(s) and D(s) are coprime). where P(s) = N(s) D(s)

15 Data Based Design of 3 Term Controllers p. 5/10 Assumption The plant is stabilizable (i.e., N(s) and D(s) are coprime). where P(s) = N(s) D(s) The plant has no jω poles.

16 Data Based Design of 3 Term Controllers p. 5/10 Assumption The plant is stabilizable (i.e., N(s) and D(s) are coprime). where P(s) = N(s) D(s) The plant has no jω poles. Only available information

17 Data Based Design of 3 Term Controllers p. 5/10 Assumption The plant is stabilizable (i.e., N(s) and D(s) are coprime). where P(s) = N(s) D(s) The plant has no jω poles. Only available information P(jω) for ω [0, )

18 Data Based Design of 3 Term Controllers p. 5/10 Assumption The plant is stabilizable (i.e., N(s) and D(s) are coprime). where P(s) = N(s) D(s) The plant has no jω poles. Only available information P(jω) for ω [0, ) Number of plant RHP poles, p +.

19 Data Based Design of 3 Term Controllers p. 6/10 Preliminaries Consider a real rational function R(s) = A(s) B(s) where A(s) and B(s) are polynomials of real coefficients with degrees m and n, respectively. Assume that A(s) and B(s) have no zeros on the jω axis.

20 Data Based Design of 3 Term Controllers p. 6/10 Preliminaries Consider a real rational function R(s) = A(s) B(s) where A(s) and B(s) are polynomials of real coefficients with degrees m and n, respectively. Assume that A(s) and B(s) have no zeros on the jω axis. Write R(jω) = R r (ω) + jr i (ω) where R r (ω) and R i (ω) are real rational functions in ω. Note that R r (ω) and R i (ω) have no real poles for ω (,+ ).

21 Data Based Design of 3 Term Controllers p. 6/10 Preliminaries Consider a real rational function R(s) = A(s) B(s) where A(s) and B(s) are polynomials of real coefficients with degrees m and n, respectively. Assume that A(s) and B(s) have no zeros on the jω axis. Write R(jω) = R r (ω) + jr i (ω) where R r (ω) and R i (ω) are real rational functions in ω. Note that R r (ω) and R i (ω) have no real poles for ω (,+ ). Let 0 = ω 0 < ω 1 < ω 2 < < ω l 1 and define ω l = denote the finite zeros of R i (ω) = 0

22 Data Based Design of 3 Term Controllers p. 7/10 Real Hurwitz Signature Lemma Define sgn[x] = 8 >< >: +1 if x > 0 0 if x = 0 1 if x < 0

23 Data Based Design of 3 Term Controllers p. 7/10 Real Hurwitz Signature Lemma Define If n m is even sgn[x] = 8 >< >: +1 if x > 0 0 if x = 0 1 if x < 0 σ(r) =! l 1 sgn[r r (ω + 0 )] + 2 X ( 1) j sgn[r r (ω j )] + ( 1) l sgn[r r (ω l )] j=1 ( 1) l 1 sgn[r i ( )] If n m is odd! Xl 1 σ(r) = sgn[r r (ω 0 )] + 2 ( 1) j sgn[r r (ω j )] ( 1) l 1 sgn[r i ( )] j=1

24 Data Based Design of 3 Term Controllers p. 8/10 Complex Hurwitz Signature Lemma Consider a complex rational function Q(s) = D(s) E(s), Q(jω) = Q r(ω) + jq i (ω) where Q r (ω) and Q i (ω) are real rational functions. Q r (ω) and Q i (ω) have no real poles for ω (,+ ).

25 Data Based Design of 3 Term Controllers p. 8/10 Complex Hurwitz Signature Lemma Consider a complex rational function Q(s) = D(s) E(s), Q(jω) = Q r(ω) + jq i (ω) where Q r (ω) and Q i (ω) are real rational functions. Q r (ω) and Q i (ω) have no real poles for ω (,+ ). Lemma σ(q) =! Xl 1 ( 1) l 1 j sgn[q r (ω j )] sgn[q i ( )] j=1

26 Data Based Design of 3 Term Controllers p. 9/10 PID controller design + C(s) P(s)

27 Data Based Design of 3 Term Controllers p. 9/10 PID controller design + C(s) P(s) Let the PID controller be of the form C(s) = K i + K p s + K d s 2 s(1 + st), T > 0

28 Data Based Design of 3 Term Controllers p. 10/10 Main Results The complete set of stabilizing PID gains for a given LTI plant can be found from the frequency response data P(jω) and the knowledge of the number of RHP poles.

29 Data Based Design of 3 Term Controllers p. 10/10 Main Results The complete set of stabilizing PID gains for a given LTI plant can be found from the frequency response data P(jω) and the knowledge of the number of RHP poles. Using the result above, the subset of the PID gains that satisfy the several given performance requirements.

30 Data Based Design of 3 Term Controllers p. 11/10 PID Controller Design Let us consider the plant and PID controller pair of the form: P(s) = N(s) D(s) C(s) = K i + K p s + K d s 2 s(1 + st) where deg[d(s)] = n and deg[n(s)] = m., T > 0

31 Data Based Design of 3 Term Controllers p. 11/10 PID Controller Design Let us consider the plant and PID controller pair of the form: P(s) = N(s) D(s) C(s) = K i + K p s + K d s 2 s(1 + st) where deg[d(s)] = n and deg[n(s)] = m., T > 0 Then consider the rational function F(s) = s(1 + st) + `K i + K p s + K d s 2 P(s)

32 Data Based Design of 3 Term Controllers p. 11/10 PID Controller Design Let us consider the plant and PID controller pair of the form: P(s) = N(s) D(s) C(s) = K i + K p s + K d s 2 s(1 + st) where deg[d(s)] = n and deg[n(s)] = m., T > 0 Then consider the rational function F(s) = s(1 + st) + `K i + K p s + K d s 2 P(s) Closed-loop stability is equivalent to the condition that zeros of F(s) lie in the LHP. This is also equivalent to the condition σ(f(s)) = n + 2 (p p + )

33 Data Based Design of 3 Term Controllers p. 12/10 PID Controller Design Continue... Now consider the rational function F(s) = F(s)P( s)

34 Data Based Design of 3 Term Controllers p. 12/10 PID Controller Design Continue... Now consider the rational function F(s) = F(s)P( s) Note that σ ` F(s) = σ(f(s)) + σ(p( s)) = n + 2 (p p + ) + (z + z ) (p + p ) = n + 2 z + + z = n m {z } 2z relative degree of P(s)

35 Data Based Design of 3 Term Controllers p. 13/10 PID Controller Design Continue... Write F(jω) = jω(1 + jωt)p( jω) + (K i + jωk p ω 2 K d )P(jω)P( jω) = F r (ω) + j F i (ω) where F r (ω) = (K i K d ω 2 ) P(jω) 2 ω 2 TP r (ω) + ωp i (ω) F i (ω) = ω `K p P(jω) 2 + P r (ω) + ωtp i (ω)

36 Data Based Design of 3 Term Controllers p. 13/10 PID Controller Design Continue... Write F(jω) = jω(1 + jωt)p( jω) + (K i + jωk p ω 2 K d )P(jω)P( jω) = F r (ω) + j F i (ω) where F r (ω) = (K i K d ω 2 ) P(jω) 2 ω 2 TP r (ω) + ωp i (ω) F i (ω) = ω `K p P(jω) 2 + P r (ω) + ωtp i (ω) K p appears only in F i (ω) K i, K d only in F r (ω)

37 Data Based Design of 3 Term Controllers p. 14/10 Lemma: Determining the required signature Relative Degree and Net Phase Change: A. In the Bode magnitude plot of the LTI system P(jω), the high frequency slope is (n m)20db/decade where n m is the relative degree of the plant P(s).

38 Data Based Design of 3 Term Controllers p. 14/10 Lemma: Determining the required signature Relative Degree and Net Phase Change: A. In the Bode magnitude plot of the LTI system P(jω), the high frequency slope is (n m)20db/decade where n m is the relative degree of the plant P(s). B. The net change of phase of P(jω), ω [0, ), denoted 0 (φ) is: 0 (φ) = ˆ(n m) 2(p + z + ) π 2 where p + and z + are numbers of RHP poles and zeros of P(s), respectively.

39 Data Based Design of 3 Term Controllers p. 15/10 Proof of Lemma The statement A is obvious from the property of Bode magnitude plot.

40 Data Based Design of 3 Term Controllers p. 15/10 Proof of Lemma The statement A is obvious from the property of Bode magnitude plot. Let p and z are number of LHP poles and zeros of P(s). Then the net phase change of P(jω) for ω [0, ) is 0 (φ) = ˆ(z z + ) (p p + ) π 2 = ˆ(n m) 2(p + z + ) π 2.

41 Data Based Design of 3 Term Controllers p. 15/10 Proof of Lemma The statement A is obvious from the property of Bode magnitude plot. Let p and z are number of LHP poles and zeros of P(s). Then the net phase change of P(jω) for ω [0, ) is 0 (φ) = ˆ(z z + ) (p p + ) π 2 = ˆ(n m) 2(p + z + ) π 2. From known n m from Bode plot, given p +, and measured 0 compute z +. (φ), we can

42 Data Based Design of 3 Term Controllers p. 16/10 Determining z + and p + Assume that a known feedback controller C(s) stabilizes the plant P and the closed-loop response can be measured and denoted by G(jω), for ω [0, )

43 Data Based Design of 3 Term Controllers p. 16/10 Determining z + and p + Assume that a known feedback controller C(s) stabilizes the plant P and the closed-loop response can be measured and denoted by G(jω), for ω [0, ) Then P(jω) is the computed frequency response of the unstable plant P(jω) = G(jω) C(jω)(1 G(jω))

44 Data Based Design of 3 Term Controllers p. 16/10 Determining z + and p + Assume that a known feedback controller C(s) stabilizes the plant P and the closed-loop response can be measured and denoted by G(jω), for ω [0, ) Then P(jω) is the computed frequency response of the unstable plant P(jω) = G(jω) C(jω)(1 G(jω)) Knowledge of C(s) and G(jω) is sufficient to determine z + and p +.

45 Data Based Design of 3 Term Controllers p. 17/10 Determining z + and p + (continue...) Lemma: z + = 1 2 ˆ rp r C 2z + c σ(g) p + = 1 2 [σ(p) σ(g) r C] z + c where z + c denotes the number of RHP zeros of C(s).

46 Data Based Design of 3 Term Controllers p. 18/10 Proof of Lemma G(s) = P(s)C(s) 1 + P(s)C(s) since G(s) is stable, σ(g) = `z + z c `z+ + z + c (n + nc ) = r P r C 2z + c 2z+ which implies z + = 1 2 ˆ rp r C 2z + c σ(g).

47 Data Based Design of 3 Term Controllers p. 18/10 Proof of Lemma G(s) = P(s)C(s) 1 + P(s)C(s) since G(s) is stable, σ(g) = `z + z c `z+ + z + c (n + nc ) = r P r C 2z + c 2z+ which implies z + = 1 2 ˆ rp r C 2z + c σ(g). From σ(p) applied to P(s), we have p + = z σ(p) r P. Then we now have p + = 1 2 [σ(p) σ(g) r C] z + c.

48 Data Based Design of 3 Term Controllers p. 19/10 Determining thezerosof F i (ω) = 0 Recall F i (ω) = ω(k p P(jω) 2 + P r (ω) + ωtp i (ω)) = 0

49 Data Based Design of 3 Term Controllers p. 19/10 Determining thezerosof F i (ω) = 0 Recall F i (ω) = ω(k p P(jω) 2 + P r (ω) + ωtp i (ω)) = 0 Then for ω 0, or K p = P r(ω) + ωtp i (ω) P(jω) 2 cos φ(ω) + ωt sin φ(ω) K p = P(jω)

50 Data Based Design of 3 Term Controllers p. 20/10 PID Controller Design Continue... The set of PID stabilizing controllers can be found as follows: Fix K = K p and set F i (ω, K p ) = 0.

51 Data Based Design of 3 Term Controllers p. 20/10 PID Controller Design Continue... The set of PID stabilizing controllers can be found as follows: Fix K = K p and set F i (ω, K p ) = 0. Let i t {+1,0, 1} and j {+1, 1}.

52 Data Based Design of 3 Term Controllers p. 20/10 PID Controller Design Continue... The set of PID stabilizing controllers can be found as follows: Fix K = K p and set F i (ω, K p ) = 0. Let i t {+1,0, 1} and j {+1, 1}. Determine strings of integers {i 0, i 1, } satisfying:

53 Data Based Design of 3 Term Controllers p. 20/10 PID Controller Design Continue... The set of PID stabilizing controllers can be found as follows: Fix K = K p and set F i (ω, K p ) = 0. Let i t {+1,0, 1} and j {+1, 1}. Determine strings of integers {i 0, i 1, } satisfying: For n m even: o ni 0 2i 1 + 2i ( 1) i 1 2i l 1 + ( 1) l i l ( 1) l 1 j = n m + 2z + + 2

54 Data Based Design of 3 Term Controllers p. 20/10 PID Controller Design Continue... The set of PID stabilizing controllers can be found as follows: Fix K = K p and set F i (ω, K p ) = 0. Let i t {+1,0, 1} and j {+1, 1}. Determine strings of integers {i 0, i 1, } satisfying: For n m even: o ni 0 2i 1 + 2i ( 1) i 1 2i l 1 + ( 1) l i l ( 1) l 1 j = n m + 2z For n m odd: i0 2i 1 + 2i ( 1) i 1 2i l 1 ( 1) l 1 j = n m + 2z

55 Data Based Design of 3 Term Controllers p. 21/10 PID Controller Design Continue... For each string (sign sequence of appropriate number of terms) satisfying the signature formula, the conditions for stability are: sgn [R(ω t, K i, K d )] i t > 0, for t = 0, 1, 2,,

56 Data Based Design of 3 Term Controllers p. 21/10 PID Controller Design Continue... For each string (sign sequence of appropriate number of terms) satisfying the signature formula, the conditions for stability are: sgn [R(ω t, K i, K d )] i t > 0, for t = 0, 1, 2,, Each valid string produces a set of linear inequalities in (K i, K d ) space.

57 Data Based Design of 3 Term Controllers p. 22/10 Example Consider the stable plant 0 Bode Diagram Magnitude (db) Phase (deg) Frequency (rad/sec) 6 π 2 = (n m) {z } =2 2( p + z + ) {z} =0! π 2 z+ = 2

58 Data Based Design of 3 Term Controllers p. 23/10 Example (continue...) The required signature for stability can now be determined and is σ( F) = (n m) + 2z = (2) + 2(2) + 2 = 8.

59 Data Based Design of 3 Term Controllers p. 23/10 Example (continue...) The required signature for stability can now be determined and is σ( F) = (n m) + 2z = (2) + 2(2) + 2 = 8. Since n m is even, we have i 0 2i 1 + 2i 2 2i 3 + 2i 4 + ( 1) l i l = 8, and it is clear that at least four terms are required to satisfy the above.

60 Data Based Design of 3 Term Controllers p. 23/10 Example (continue...) The required signature for stability can now be determined and is σ( F) = (n m) + 2z = (2) + 2(2) + 2 = 8. Since n m is even, we have i 0 2i 1 + 2i 2 2i 3 + 2i 4 + ( 1) l i l = 8, and it is clear that at least four terms are required to satisfy the above. In other words l 4.

61 Data Based Design of 3 Term Controllers p. 24/10 Example (continue...) From the figure it is easy to see that Kp solutions and therefore we have has at most three positive frequencies as i 0 2i 1 + 2i 2 2i 3 + i 4 = necessary range of ω 10 0 feasible range of K p ω

62 Data Based Design of 3 Term Controllers p. 25/10 Example (continue...) Fix K p = 1 and compute the set of ω s that satisfies cos φ(ω) + ωt sin φ(ω) P(jω) = ω 1 =0.742 ω 2 =1.865 ω 3 = K p 10 K p = ω

63 Data Based Design of 3 Term Controllers p. 26/10 Example (continue...) This leads to the requirement i 0 2i 1 + 2i 2 2i 3 = 7

64 Data Based Design of 3 Term Controllers p. 26/10 Example (continue...) This leads to the requirement i 0 2i 1 + 2i 2 2i 3 = 7 giving the feasible string F = {i 0, i 1, i 2, i 3 } = {1, 1, 1, 1}.

65 Data Based Design of 3 Term Controllers p. 26/10 Example (continue...) This leads to the requirement i 0 2i 1 + 2i 2 2i 3 = 7 giving the feasible string F = {i 0, i 1, i 2, i 3 } = {1, 1, 1, 1}. Thus, we have the following set of linear inequalities for stability: K i > K i K d < K i K d > K i K d < 0

66 Data Based Design of 3 Term Controllers p. 27/10 Example: Stabilizing PID Setfor K p = Kdbar Kibar

67 Data Based Design of 3 Term Controllers p. 28/10 Example: Entire Stabilizing PID Set 5 0 Kpbar Kdbar Kibar

68 Data Based Design of 3 Term Controllers p. 29/10 Performance Specifications Many performance attainment problems can be cast as stabilization of families of real and complex plants. For example, The problem of achieving a gain margin is equivalent to stabilizing the family of real plants P c (s) = {KP(s) : K [K min, K max ]}.

69 Data Based Design of 3 Term Controllers p. 29/10 Performance Specifications Many performance attainment problems can be cast as stabilization of families of real and complex plants. For example, The problem of achieving a gain margin is equivalent to stabilizing the family of real plants P c (s) = {KP(s) : K [K min, K max ]}. The problem of achieving prescribed phase margin θ m is equivalent to stabilizing the family of complex plants P c (s) = n o e jθ P(s) : θ [0, θ m ].

70 Data Based Design of 3 Term Controllers p. 30/10 The problem of achieving an H norm specification on the sensitivity function S(s), that is, W(s)S(s) < γ is equivalent to stabilizing the family of complex plants (" # ) P c 1 (s) = P(s) : θ [0,2π]. γ ejθ W(s)

71 Data Based Design of 3 Term Controllers p. 30/10 The problem of achieving an H norm specification on the sensitivity function S(s), that is, W(s)S(s) < γ is equivalent to stabilizing the family of complex plants (" # ) P c 1 (s) = P(s) : θ [0,2π]. γ ejθ W(s) The problem of achieving an H norm specification on the complementary sensitivity function T(s), that is, W(s)T(s) < γ is equivalent to stabilizing the family of complex plants P c (s) = j» P(s) γ ejθ W(s) ff : θ [0,2π].

72 Data Based Design of 3 Term Controllers p. 31/10 Assumption The only information available to the designer is: Knowledge of the frequency response magnitude and phase, equivalently, P c (jω), ω (, + ). Knowledge of the number of RHP poles, p +.

73 Data Based Design of 3 Term Controllers p. 32/10 Determining Performance Set The complete set of stabilizing PID gains for a given complex LTI plant can be found from the frequency response data P c (jω) and the knowledge of the number of RHP poles

74 Data Based Design of 3 Term Controllers p. 32/10 Determining Performance Set The complete set of stabilizing PID gains for a given complex LTI plant can be found from the frequency response data P c (jω) and the knowledge of the number of RHP poles The set of stabilizing PID gains can be computed by the following procedure:

75 Data Based Design of 3 Term Controllers p. 32/10 Determining Performance Set The complete set of stabilizing PID gains for a given complex LTI plant can be found from the frequency response data P c (jω) and the knowledge of the number of RHP poles The set of stabilizing PID gains can be computed by the following procedure: Determine the relative degree n c m c from the high frequency slope of the Bode magnitude plot where n c and m c are degrees of numerator and denominator of P c.

76 Data Based Design of 3 Term Controllers p. 32/10 Determining Performance Set The complete set of stabilizing PID gains for a given complex LTI plant can be found from the frequency response data P c (jω) and the knowledge of the number of RHP poles The set of stabilizing PID gains can be computed by the following procedure: Determine the relative degree n c m c from the high frequency slope of the Bode magnitude plot where n c and m c are degrees of numerator and denominator of P c. Fix K p = K p and solve Kp cos φ(ω) + ωt sin φ(ω) = P c (jω) and let ω 1 < ω 2 < < ω l 1 denote the distinct frequencies which are solutions of the above.

77 Data Based Design of 3 Term Controllers p. 33/10 Set ω 0 =, ω l = + and determine all strings of integers i t {+1,0, 1} and j { 1,+1} such that l 1 X ( 1) l 1 r i r j = n c m c + 2z c r=1 where n c and m c denote the numerator and denominator degrees of P c (s) and z + c the number of RHP zeros.

78 Data Based Design of 3 Term Controllers p. 33/10 Set ω 0 =, ω l = + and determine all strings of integers i t {+1,0, 1} and j { 1,+1} such that l 1 X ( 1) l 1 r i r j = n c m c + 2z c r=1 where n c and m c denote the numerator and denominator degrees of P c (s) and z + c the number of RHP zeros. For the fixed K p = K p chosen in Step 1, solve for the stabilizing (K i, K d ) from: for t = 0,1,.» K i K d ωt 2 + ω t sin φ(ω t ) ωt 2T cos φ(ω t) P c i t > 0 (jω)

79 Data Based Design of 3 Term Controllers p. 33/10 Set ω 0 =, ω l = + and determine all strings of integers i t {+1,0, 1} and j { 1,+1} such that l 1 X ( 1) l 1 r i r j = n c m c + 2z c r=1 where n c and m c denote the numerator and denominator degrees of P c (s) and z + c the number of RHP zeros. For the fixed K p = K p chosen in Step 1, solve for the stabilizing (K i, K d ) from: for t = 0,1,.» K i K d ωt 2 + ω t sin φ(ω t ) ωt 2T cos φ(ω t) P c i t > 0 (jω) Repeat the previous three steps by updating K p over prescribed ranges.

80 Data Based Design of 3 Term Controllers p. 34/10 Performance Example Taking the same frequency domain data set P(jω) used in the previous example, we consider the problem of achieving an H norm specification on the complementary sensitivity function T(s), that is, W(s)T(s) < 1 where W(s) = s s + 1. By solving the complex stabilization problem, we have the stabilizing PID controller parameter region that satisfies the given H norm specification.

81 Data Based Design of 3 Term Controllers p. 35/10 Thecomplete set ofstabilizingpidgainsfor H specification when K p = K d 0 2 Stabilizing region for P c (s) 4 Stabilizing region for P(s) K i

82 Data Based Design of 3 Term Controllers p. 36/10 Nyquist plot of W(s)T(s) By selecting a point, we verify that the point selected satisfied the given H specification imag real

83 Data Based Design of 3 Term Controllers p. 37/10 Entireset of stabilizingpidgainssatisfyingthe H specification By sweeping K p, we have the entire stabilizing PID gains that satisfy the given H specification as shown in Figure. 2 0 K p K d K i 2 3 4

84 Data Based Design of 3 Term Controllers p. 38/10 Example with Gain/Phase Margin Specification Consider the following nonminimum phase plant: P(s) =» s 3 4s 2 + s + 2 s 5 + 8s s s s + 17 e s

85 Data Based Design of 3 Term Controllers p. 38/10 Example with Gain/Phase Margin Specification Consider the following nonminimum phase plant: P(s) =» s 3 4s 2 + s + 2 s 5 + 8s s s s + 17 e s We wish to find the entire set of controllers that simultaneously satisfies the following specifications: PID controllers must stabilize the given plant with a delay The closed-loop system must guarantee the following gain and phase margins:

86 Data Based Design of 3 Term Controllers p. 38/10 Example with Gain/Phase Margin Specification Consider the following nonminimum phase plant: P(s) =» s 3 4s 2 + s + 2 s 5 + 8s s s s + 17 e s We wish to find the entire set of controllers that simultaneously satisfies the following specifications: PID controllers must stabilize the given plant with a delay The closed-loop system must guarantee the following gain and phase margins: Gain margin : K + 2 (about 6 [db])

87 Data Based Design of 3 Term Controllers p. 38/10 Example with Gain/Phase Margin Specification Consider the following nonminimum phase plant: P(s) =» s 3 4s 2 + s + 2 s 5 + 8s s s s + 17 e s We wish to find the entire set of controllers that simultaneously satisfies the following specifications: PID controllers must stabilize the given plant with a delay The closed-loop system must guarantee the following gain and phase margins: Gain margin : K + 2 (about 6 [db]) Phase margin : [θ, θ + ] = [ 10 o,60 o ]

88 Data Based Design of 3 Term Controllers p. 39/10 (A) All Stabilizing PID Controllers The feasible ranges of k p given are: [ 19.1, 8.5], [ 8.5,4.23], [4,23, ]. We can easily conclude that the only region that contains the solutions is [ 8.5,4.23]. The complete set of stabilizing PID controllers is: K p Kp= 6.49 Kp= Kp= Kp= 5.82 Kp= 5.15 Kp= 4.48 Kp= 3.80 Kp= 2.46 Kp= Kp= 1.12 Kp= Kp= 3.56 Kp= 2.89 Kp= 0.21 Kp= 0.46 Kp= 2.22 Kp= 1.55 Kp= k i k d

89 Data Based Design of 3 Term Controllers p. 40/10 (B)All PID Satisfying thegm and PM Requirement We obtained the final result Kp= 2.22 Kp= 1.55 Kp= 0.88 k p Kp= 0.21 Kp= 0.46 Kp= 1.12 Kp= Kp= 3.13 Kp= k d k i 2

90 S C k p 2-D Regions with k p = 0.88 S kp : stable set with a fixed k p S C k p : unstable set S kp (G) : set satisfying the GM requirement S kp (φ) : set satisfying the PM requiremen S kp (G) S kp (φ) set satisfying both gain and phase margin requirements 5 Kp= S kp (G) S kp Kd S kp (G) S kp (φ) S kp (φ) Ki Data Based Design of 3 Term Controllers p. 41/10

91 First Order Controller Design Data Based Design of 3 Term Controllers p. 42/10

92 Data Based Design of 3 Term Controllers p. 42/10 First Order Controller Design Consider the First order controller of the form for an LTI plant. C(s) = x 1s + x 2 s + x 3

93 Data Based Design of 3 Term Controllers p. 42/10 First Order Controller Design Consider the First order controller of the form for an LTI plant. C(s) = x 1s + x 2 s + x 3 Assumption The plant has no jω poles or zeros. Available frequency domain data P(jω) for ω [0, ). Knowledge of p +.

94 Data Based Design of 3 Term Controllers p. 43/10 Root Invariant Regions Let P(jω) = P r (ω) + jp i (ω) F(s) = (s + x 3 ) + (sx 1 + x 2 )P(s)

95 Data Based Design of 3 Term Controllers p. 43/10 Root Invariant Regions Let P(jω) = P r (ω) + jp i (ω) F(s) = (s + x 3 ) + (sx 1 + x 2 )P(s) For closed-loop stability, it is necessary and sufficient that σ(f(s)) = n + 1 (p p + )

96 Data Based Design of 3 Term Controllers p. 43/10 Root Invariant Regions Let P(jω) = P r (ω) + jp i (ω) F(s) = (s + x 3 ) + (sx 1 + x 2 )P(s) For closed-loop stability, it is necessary and sufficient that σ(f(s)) = n + 1 (p p + ) Let F(s) = F(s)P( s) and the stability condition is F(s) = (s + x 3 )P( s) + (sx 1 + x 2 )P(s)P( s)

97 Data Based Design of 3 Term Controllers p. 44/10 Root Invariant Region (continue...) In other words, F(jω, x 1, x 2, x 3 ) = x 2 P(jω) 2 + ωp i (ω) + x 3 P r (ω) {z } F r (ω,x 1,x 2,x 3 ) +jω `x 1 P(jω) 2 x 3 P i (ω) + P r (ω). {z } F i (ω,x 1,x 2,x 3 )

98 Data Based Design of 3 Term Controllers p. 44/10 Root Invariant Region (continue...) In other words, F(jω, x 1, x 2, x 3 ) = x 2 P(jω) 2 + ωp i (ω) + x 3 P r (ω) {z } F r (ω,x 1,x 2,x 3 ) +jω `x 1 P(jω) 2 x 3 P i (ω) + P r (ω). {z } F i (ω,x 1,x 2,x 3 ) The curves F r (ω, ) = 0 and F i (ω, ) = 0, 0 ω < along with the F(0, ) = 0 and F(, ) = 0 partition the (x 1, x 2, x 3 ) parameter space into signature invariant regions.

99 Data Based Design of 3 Term Controllers p. 44/10 Root Invariant Region (continue...) In other words, F(jω, x 1, x 2, x 3 ) = x 2 P(jω) 2 + ωp i (ω) + x 3 P r (ω) {z } F r (ω,x 1,x 2,x 3 ) +jω `x 1 P(jω) 2 x 3 P i (ω) + P r (ω). {z } F i (ω,x 1,x 2,x 3 ) The curves F r (ω, ) = 0 and F i (ω, ) = 0, 0 ω < along with the F(0, ) = 0 and F(, ) = 0 partition the (x 1, x 2, x 3 ) parameter space into signature invariant regions. By plotting these curves and selecting a test point from each of these regions we can determine the stability regions corresponding to those with signature equal to n m + 2z

100 Data Based Design of 3 Term Controllers p. 45/10 Procedure for First Order Controller Design 1. Determine the relative degree n m from the high frequency slope of the Bode magnitude plot.

101 Data Based Design of 3 Term Controllers p. 45/10 Procedure for First Order Controller Design 1. Determine the relative degree n m from the high frequency slope of the Bode magnitude plot. 2. Let 0 [φ(ω)] denote the net change of phase P(jω) for ω [0, ). Determine z + from knowledge of p + and 0 [φ(ω)] = ˆ(n m) 2(p + z + ) π 2.

102 Data Based Design of 3 Term Controllers p. 45/10 Procedure for First Order Controller Design 1. Determine the relative degree n m from the high frequency slope of the Bode magnitude plot. 2. Let 0 [φ(ω)] denote the net change of phase P(jω) for ω [0, ). Determine z + from knowledge of p + and 0 [φ(ω)] = ˆ(n m) 2(p + z + ) π Plot the curves below in the (x 1, x 2 ) plane for a fixed x 3. x 3 + x 2 P(0) = 0 8 < : x 1 (ω) = 1 P(jω) sin φ(ω) x ω 3 cos φ(ω), for 0 < ω < x 2 (ω) = 1 P(jω) (cos φ(ω)x 3 + ω sin φ(ω)), 1 + P( )x 2 = 0. for 0 < ω <

103 Data Based Design of 3 Term Controllers p. 46/10 Procedure for First Order Controller Design (continue...) 4. The curves x 1 (ω) and x 2 (ω) partition the (x 1, x 2 ) plane into disjoint signature invariant regions. The stabilizing regions correspond to those for which F(s) has a signature of n m + 2z

104 Data Based Design of 3 Term Controllers p. 47/10 Example For illustration, we have collected the frequency domain (Nyquist-Bode) data of a stable plant: P(jω) = {P(jω) : ω (0,10) sampled every 0.01}.

105 Data Based Design of 3 Term Controllers p. 47/10 Example For illustration, we have collected the frequency domain (Nyquist-Bode) data of a stable plant: P(jω) = {P(jω) : ω (0,10) sampled every 0.01}. The Nyquist plot of the plant obtained is: 50 Nyquist Diagram 5 Nyquist Diagram Imaginary Axis Magnification shown in Fig. 3 Imaginary Axis Real Axis Real Axis

106 Data Based Design of 3 Term Controllers p. 48/10 Example (continue...) From the data P(jω), we have P(0) = and P( ) = 0. Then it is easy to see that the straight line is not applicable.

107 Data Based Design of 3 Term Controllers p. 48/10 Example (continue...) From the data P(jω), we have P(0) = and P( ) = 0. Then it is easy to see that the straight line is not applicable. After fixing x 3 = 0.2, the data points representing the straight line and the curve are depicted.

108 Data Based Design of 3 Term Controllers p. 48/10 Example (continue...) From the data P(jω), we have P(0) = and P( ) = 0. Then it is easy to see that the straight line is not applicable. After fixing x 3 = 0.2, the data points representing the straight line and the curve are depicted. By testing a point for each root invariant region, we have the following x Stabilizing regions x 1

109 Data Based Design of 3 Term Controllers p. 49/10 Example (continue...) We now consider the problem of determining the entire set of first order stabilizing controllers satisfying the required closed-loop performance described by the requirement on the H norm of the weighted complementary sensitivity function: W(jω)T(jω) < γ, for all ω

110 Data Based Design of 3 Term Controllers p. 49/10 Example (continue...) We now consider the problem of determining the entire set of first order stabilizing controllers satisfying the required closed-loop performance described by the requirement on the H norm of the weighted complementary sensitivity function: W(jω)T(jω) < γ, for all ω This is equivalent to the problem of simultaneously stabilizing the specified complex family as well as the original plant P(s).

111 Data Based Design of 3 Term Controllers p. 50/10 Example (continue...) In this problem, we let γ = 1 and x 3 = 2.5. We superimpose on the top of the stabilizing region shown (left) for the real plant, the stabilizing sets for the complex plant families P c (jω, θ) for θ = 0, π 3, 2π 3, π, 4π 3, 5π, 2π are plotted (right) x x stabilizing first order set x x 1

112 Data Based Design of 3 Term Controllers p. 51/10 Example (continue...) To verify, a number of points inside the performance region, construct the corresponding controllers, and Nyquist plots of W(s)T(s) have been plotted as shown. These points are shown as * in (right). We observe from the Nyquist plots and every test set satisfies the H performance requirement imag real

113 Data Based Design of 3 Term Controllers p. 52/10 Example of FO Controller Design Consider an unstable plant with time delay. G(s) =» s + 1 s 4 + 8s s s 1 e s The design specifications are given as follows: the closed-loop system must satisfy the following gain and phase margin requirements:

114 (A) The Complete Set of Stabilizing FO Controllers The feasible range is found as: x 3 [ 6, ]. We chose to execute the algorithm for x 3 [ 6,70]. Then we obtain the the set of FO Controller parameters so that the each and every corresponding closed-loop system satisfies the given gain and phase margin requirements. X X3=70.00 X3=61.55 X3=53.11 X3=44.66 X3=36.22 X3=27.77 X3= X3= X3= Data Based Design of 3 Term Controllers p. 53/10

115 Data Based Design of 3 Term Controllers p. 54/10 (B) A Set of FOC that Meets the Desired Overshoot and Settling Time Requirements Computing a feasible range of the generalized time constant based on the CRA, the feasible range of the time constant is obtained as: τ [2.692, 6.258] x x 2 x

116 Data Based Design of 3 Term Controllers p. 55/10 Stepresponseswith FO Controllers in S Step Responses y(t) t[s]

117 2-D setof FO Controllers, S x 3,when x 3 = 53.1 To proceed, we pick a FO Controller from the controller set S and examine various performance of the corresponding closed-loop system. In this example, we first select x 3 = 53.1 and the corresponding 2-D set is depicted as in figure X3= S x Sx3 (G) S x3 (φ) S x3 X S x3 (G) S x3 (φ) S x3 (T) X1 Data Based Design of 3 Term Controllers p. 56/10

118 Data Based Design of 3 Term Controllers p. 57/10 Various characteristics of the CL system with the selected FOC A selection of the controller point from the set S x 3 with the fixed x 3 enables us to display the following six figures displaying various characteristics of the closed-loop system with the FO Controller chosen. X3=53.1 Step Reponse X2 800 y(t) X t[s] Roots Space Im 0 Imag Re Gain and Phase Margin Gm=( 24.8, 11.0) [db] Pm=( 180.0, 66.4) [deg] u(t) Real Control Input 1 t[s]

119 Data Based Design of 3 Term Controllers p. 58/10 The top left figure shows the controller set S x 3 and + indicates the selected controller of the parameter values: x 1 = 83.3, x 2 = 831, x 3 = 53.1.

120 Data Based Design of 3 Term Controllers p. 58/10 The top left figure shows the controller set S x 3 and + indicates the selected controller of the parameter values: x 1 = 83.3, x 2 = 831, x 3 = The top right figure shows the step response of the closed-loop system. Two figures in the middle show the Nyquist plot, and the closed-loop poles and zeros. The bottom two figures show the values of gain and phase margins, and the control input signal, respectively.

121 Data Based Design of 3 Term Controllers p. 59/10 Data-Robust Design (An Example) What we know: P(jω) ± 20% for all ω

122 Data Based Design of 3 Term Controllers p. 59/10 Data-Robust Design (An Example) What we know: P(jω) ± 20% for all ω Number of RHP zeros is z + = 2

123 Data Based Design of 3 Term Controllers p. 59/10 Data-Robust Design (An Example) What we know: P(jω) ± 20% for all ω Number of RHP zeros is z + = 2 Relative degree is r P = 1 due to high frequency slope (-20db/decade)

124 Data Based Design of 3 Term Controllers p. 59/10 Data-Robust Design (An Example) What we know: P(jω) ± 20% for all ω Number of RHP zeros is z + = 2 Relative degree is r P = 1 due to high frequency slope (-20db/decade) Signature required is r P + 2z = 1 + 2(2) + 2 = 7.

125 Data Based Design of 3 Term Controllers p. 59/10 Data-Robust Design (An Example) What we know: P(jω) ± 20% for all ω Number of RHP zeros is z + = 2 Relative degree is r P = 1 due to high frequency slope (-20db/decade) Signature required is r P + 2z = 1 + 2(2) + 2 = 7. It requires at least 3 frequencies to satisfy K p = g(ω). Notations: P(jω) = Family of response shown in Bode plot P max r (ω) = max P(jω) P(jω) P r(ω), P min r (ω) = min P(jω) P(jω) P r(ω), P max i (ω) = max P(jω) P(jω) P i(ω), P min i (ω) = min P(jω) P(jω) P i(ω), for every ω for every ω for every ω for every ω

126 Data Based Design of 3 Term Controllers p. 59/10 Data-Robust Design (An Example) What we know: P(jω) ± 20% for all ω Number of RHP zeros is z + = 2 Relative degree is r P = 1 due to high frequency slope (-20db/decade) Signature required is r P + 2z = 1 + 2(2) + 2 = 7. It requires at least 3 frequencies to satisfy K p = g(ω). Notations: P(jω) = Family of response shown in Bode plot P max r (ω) = max P(jω) P(jω) P r(ω), P min r (ω) = min P(jω) P(jω) P r(ω), P max i (ω) = max P(jω) P(jω) P i(ω), P min i (ω) = min P(jω) P(jω) P i(ω), for every ω for every ω for every ω for every ω

127 Data-Robust Design (An Example) Continue... Data Based Design of 3 Term Controllers p. 60/10

128 Data Based Design of 3 Term Controllers p. 60/10 Data-Robust Design (An Example) Continue... K = P r(ω) + ωtp i (ω) P(jω) 2 := g(ω)

129 Data Based Design of 3 Term Controllers p. 60/10 Data-Robust Design (An Example) Continue... K = P r(ω) + ωtp i (ω) P(jω) 2 := g(ω) We overbound the g(ω) family: 8 < g max (ω) = max : g(ω) : Pr max 8 < g min (ω) = min : g(ω) : Pr max (ω), P min r (ω), Pi max P(jω) max, P(jω) min (ω), Pi max P(jω) max, P(jω) min (ω), P min r (ω), Pi min (ω), (ω), Pi min (ω), 9 = ; 9 = ;

130 Data Based Design of 3 Term Controllers p. 61/10 g(ω) graph: K p = 5 20 ω ω 2 10 K p =5 g(ω) 5 0 ω 1 ω 3 ω ω 1 + ω

131 Data Based Design of 3 Term Controllers p. 62/10 Interval Linear Programming y mx c > 0, m [m, m + ], c [c, c + ].

132 Data Based Design of 3 Term Controllers p. 62/10 Interval Linear Programming y mx c > 0, m [m, m + ], c [c, c + ]. c c+ c m m+ m

133 Data Based Design of 3 Term Controllers p. 62/10 Interval Linear Programming y mx c > 0, m [m, m + ], c [c, c + ]. c c+ y (m +, c + ) (m +, c ) c m m+ m x (m, c + ) (m, c )

134 Data Based Design of 3 Term Controllers p. 62/10 Interval Linear Programming y mx c > 0, m [m, m + ], c [c, c + ]. c c+ y (m +, c + ) (m +, c ) c m m+ m x (m, c + ) (m, c ) 8 >< >: y m x c > 0 y m x c + > 0 y m + x c > 0 y m + x c + > 0

135 Data Based Design of 3 Term Controllers p. 63/10 A set of linear programmings b t = ω tp i (ω t ) + ω 2 t TP r(ω t ) P(jω t ) 2

136 Data Based Design of 3 Term Controllers p. 63/10 A set of linear programmings b t = ω tp i (ω t ) + ω 2 t TP r(ω t ) P(jω t ) 2 8 >< >: K i (ω 1 )2 K d b 1 > 0 K i (ω 2 )2 K d b 2 < 0 K i (ω 3 )2 K d b 3 > 0 K i (ω 4 )2 K d b 4 < 0

137 Data Based Design of 3 Term Controllers p. 63/10 A set of linear programmings b t = ω tp i (ω t ) + ω 2 t TP r(ω t ) P(jω t ) 2 8 >< K i (ω 1 )2 K d b 1 > 0 K i (ω 2 )2 K d b 2 < 0 8 >< K i (ω 1 )2 K d b + 1 > 0 K i (ω 2 )2 K d b + 2 < 0 >: K i (ω 3 )2 K d b 3 > 0 K i (ω 4 )2 K d b 4 < 0 >: K i (ω 3 )2 K d b + 3 > 0 K i (ω 4 )2 K d b + 4 < 0

138 Data Based Design of 3 Term Controllers p. 63/10 A set of linear programmings b t = ω tp i (ω t ) + ω 2 t TP r(ω t ) P(jω t ) 2 8 >< K i (ω 1 )2 K d b 1 > 0 K i (ω 2 )2 K d b 2 < 0 8 >< K i (ω 1 )2 K d b + 1 > 0 K i (ω 2 )2 K d b + 2 < 0 >: K i (ω 3 )2 K d b 3 > 0 K i (ω 4 )2 K d b 4 < 0 >: K i (ω 3 )2 K d b + 3 > 0 K i (ω 4 )2 K d b + 4 < 0 8 >< K i (ω + 1 )2 K d b 1 > 0 K i (ω + 2 )2 K d b 2 < 0 >: K i (ω + 3 )2 K d b 3 > 0 K i (ω + 4 )2 K d b 4 < 0

139 Data Based Design of 3 Term Controllers p. 63/10 A set of linear programmings b t = ω tp i (ω t ) + ω 2 t TP r(ω t ) P(jω t ) 2 8 >< K i (ω 1 )2 K d b 1 > 0 K i (ω 2 )2 K d b 2 < 0 8 >< K i (ω 1 )2 K d b + 1 > 0 K i (ω 2 )2 K d b + 2 < 0 >: K i (ω 3 )2 K d b 3 > 0 K i (ω 4 )2 K d b 4 < 0 >: K i (ω 3 )2 K d b + 3 > 0 K i (ω 4 )2 K d b + 4 < 0 8 >< K i (ω + 1 )2 K d b 1 > 0 K i (ω + 2 )2 K d b 2 < 0 8 >< K i (ω + 1 )2 K d b + 1 > 0 K i (ω + 2 )2 K d b + 2 < 0 >: K i (ω + 3 )2 K d b 3 > 0 K i (ω + 4 )2 K d b 4 < 0 >: K i (ω + 3 )2 K d b + 3 > 0 K i (ω + 4 )2 K d b + 4 < 0

140 Data Based Design of 3 Term Controllers p. 64/10 Stabilizing set for K p = 5 7 For k p = K d K i

141 Digital ControlDesign Data Based Design of 3 Term Controllers p. 65/10

142 Data Based Design of 3 Term Controllers p. 66/10 Digital Control Design: Objectives To synthesize the set of all stabilizing PID controllers from experimental data of the discrete-time plant rather than from mathematical models.

143 Data Based Design of 3 Term Controllers p. 66/10 Digital Control Design: Objectives To synthesize the set of all stabilizing PID controllers from experimental data of the discrete-time plant rather than from mathematical models. To synthesize the set of all stabilizing First Order controllers from experimental data of the discrete-time plant rather than from mathematical models.

144 Data Based Design of 3 Term Controllers p. 67/10 Preliminaries: Tchebyshev Decomposition Consider a real polynomial in z, P(z) = a n z n + a n 1 z n 1 + a n 2 z n a 1 z + a 0.

145 Data Based Design of 3 Term Controllers p. 67/10 Preliminaries: Tchebyshev Decomposition Consider a real polynomial in z, P(z) = a n z n + a n 1 z n 1 + a n 2 z n a 1 z + a 0. Since z k z=e jθ = coskθ + j sinkθ,

146 Data Based Design of 3 Term Controllers p. 67/10 Preliminaries: Tchebyshev Decomposition Consider a real polynomial in z, P(z) = a n z n + a n 1 z n 1 + a n 2 z n a 1 z + a 0. Since z k z=e jθ = coskθ + j sinkθ, we have P(e jθ ) = (a n cos nθ + + a 1 cos θ + a 0 ) + j (a n sinnθ + + a 1 sinθ)

147 Data Based Design of 3 Term Controllers p. 67/10 Preliminaries: Tchebyshev Decomposition Consider a real polynomial in z, P(z) = a n z n + a n 1 z n 1 + a n 2 z n a 1 z + a 0. Since z k z=e jθ = coskθ + j sinkθ, we have P(e jθ ) = (a n cos nθ + + a 1 cos θ + a 0 ) + j (a n sinnθ + + a 1 sinθ) Let then we have u = cosθ, e jθ = u + j 1 u 2.

148 Data Based Design of 3 Term Controllers p. 68/10 Tchebyshev Decomposition (continue...) We now have P(e jθ ) u= cos θ = R(u) + j p 1 u 2 T(u)

149 Data Based Design of 3 Term Controllers p. 68/10 Tchebyshev Decomposition (continue...) We now have P(e jθ ) u= cos θ = R(u) + j p 1 u 2 T(u) where R(u) = a n c n (u) + a n 1 c n 1 (u) + + a 1 c 1 (u) + a 0 T(u) = a n s n (u) + a n 1 s n 1 (u) + + a 1 s 1 (u)

150 Data Based Design of 3 Term Controllers p. 68/10 Tchebyshev Decomposition (continue...) We now have P(e jθ ) u= cos θ = R(u) + j p 1 u 2 T(u) where R(u) = a n c n (u) + a n 1 c n 1 (u) + + a 1 c 1 (u) + a 0 T(u) = a n s n (u) + a n 1 s n 1 (u) + + a 1 s 1 (u) k c k (u) s k (u) 1 u 1 2 2u 2 1 2u 3 4u 3 + 3u 4u u 4 8u u 3 + 4u 5 16u u 3 5u 16u 4 12u 2 + 1

151 Data Based Design of 3 Term Controllers p. 68/10 Tchebyshev Decomposition (continue...) We now have P(e jθ ) u= cos θ = R(u) + j p 1 u 2 T(u) where R(u) = a n c n (u) + a n 1 c n 1 (u) + + a 1 c 1 (u) + a 0 T(u) = a n s n (u) + a n 1 s n 1 (u) + + a 1 s 1 (u) k c k (u) s k (u) 1 u 1 2 2u 2 1 2u 3 4u 3 + 3u 4u u 4 8u u 3 + 4u 5 16u u 3 5u 16u 4 12u s k (u) = 1 K dc k(u) du, c k+1(u) = uc k (u) `1 u 2 s k (u), k = 1, 2,

152 Data Based Design of 3 Term Controllers p. 69/10 Preliminaries (continue...) Consider a real rational function Q(z) = P 1(z) P 2 (z) where P 1 (z) and P 2 (z) are polynomials of real coefficients with degrees m and n, respectively. Assume that P 1 (z) and P 2 (z) have no zeros on the unit circle.

153 Data Based Design of 3 Term Controllers p. 69/10 Preliminaries (continue...) Consider a real rational function Q(z) = P 1(z) P 2 (z) where P 1 (z) and P 2 (z) are polynomials of real coefficients with degrees m and n, respectively. Assume that P 1 (z) and P 2 (z) have no zeros on the unit circle. Write Q(z) z= u+j 1 u 2 = R q(u) + j p 1 u 2 T q (u)

154 Data Based Design of 3 Term Controllers p. 69/10 Preliminaries (continue...) Consider a real rational function Q(z) = P 1(z) P 2 (z) where P 1 (z) and P 2 (z) are polynomials of real coefficients with degrees m and n, respectively. Assume that P 1 (z) and P 2 (z) have no zeros on the unit circle. Write Q(z) z= u+j 1 u 2 = R q(u) + j p 1 u 2 T q (u) where R q (u) = R(u) D(u), T q(u) = T(u) D(u)

155 Data Based Design of 3 Term Controllers p. 69/10 Preliminaries (continue...) Consider a real rational function Q(z) = P 1(z) P 2 (z) where P 1 (z) and P 2 (z) are polynomials of real coefficients with degrees m and n, respectively. Assume that P 1 (z) and P 2 (z) have no zeros on the unit circle. Write Q(z) z= u+j 1 u 2 = R q(u) + j p 1 u 2 T q (u) where R q (u) = R(u) D(u), T q(u) = T(u) D(u) R(u) = R 1 (u)r 2 (u) + `1 u 2 T 1 (u)t 2 (u) T(u) = T 1 (u)r 2 (u) R 1 (u)t 2 (u) D(u) = R 2 2(u) + `1 u 2 T 2 2 (u)

156 Data Based Design of 3 Term Controllers p. 70/10 Preliminaries (continue...) Recall Q(z) z= u+j 1 u 2 = R(u) D(u) + j 1 u 2 T(u) D(u)

157 Data Based Design of 3 Term Controllers p. 70/10 Preliminaries (continue...) Recall Q(z) z= u+j 1 u 2 = R(u) D(u) + j 1 u 2 T(u) D(u) Since D(u) = R 2 2(u) + ( 1 u 2) T 2 2 (u) > 0 for all u [ 1, 1], the zeros of T q (u) are identical to those of T(u).

158 Data Based Design of 3 Term Controllers p. 70/10 Preliminaries (continue...) Recall Q(z) z= u+j 1 u 2 = R(u) D(u) + j 1 u 2 T(u) D(u) Since D(u) = R 2 2(u) + ( 1 u 2) T 2 2 (u) > 0 for all u [ 1, 1], the zeros of T q (u) are identical to those of T(u). Let t 1,,t k denote the real distinct zeros of T(u) of odd multiplicity ordered as 1 < t 1 < t 2 < < t k < +1.

159 Data Based Design of 3 Term Controllers p. 71/10 Schur Signature Lemma Let Q(z) be a real rational function with i z zeros and i p poles, respectively, inside the unit circle C and no zeros/poles on the unit circle and Q(z) z= u+ 1 u 2 = R(u) + j 1 u 2 T(u).

160 Data Based Design of 3 Term Controllers p. 71/10 Schur Signature Lemma Let Q(z) be a real rational function with i z zeros and i p poles, respectively, inside the unit circle C and no zeros/poles on the unit circle and Q(z) z= u+ 1 u 2 = R(u) + j 1 u 2 T(u). Then the net change in phase of Q(e jθ ) as θ runs from 0 to π: π 0 Q(e jθ ) = π (i z i p ) = πσ[q(z)]

161 Data Based Design of 3 Term Controllers p. 71/10 Schur Signature Lemma Let Q(z) be a real rational function with i z zeros and i p poles, respectively, inside the unit circle C and no zeros/poles on the unit circle and Q(z) z= u+ 1 u 2 = R(u) + j 1 u 2 T(u). Then the net change in phase of Q(e jθ ) as θ runs from 0 to π: The signature is π 0 Q(e jθ ) = π (i z i p ) = πσ[q(z)] σ[q(z)] = i z i p = 1 ht ( 1)i 0 1 kx 2 sgn 1)] + 2 ( 1) j sgn[r(t j )] + ( 1) k+1 sgn[r(+1)] A. j=1

162 Data Based Design of 3 Term Controllers p. 72/10 Assumption The plant is stabilizable (mathematically, N(z) and D(z) are coprime). where P(z) = N(z) D(z)

163 Data Based Design of 3 Term Controllers p. 72/10 Assumption The plant is stabilizable (mathematically, N(z) and D(z) are coprime). where P(z) = N(z) D(z) The plant has no poles on the unit circle.

164 Data Based Design of 3 Term Controllers p. 72/10 Assumption The plant is stabilizable (mathematically, N(z) and D(z) are coprime). where P(z) = N(z) D(z) The plant has no poles on the unit circle. Only available information:

165 Data Based Design of 3 Term Controllers p. 72/10 Assumption The plant is stabilizable (mathematically, N(z) and D(z) are coprime). where P(z) = N(z) D(z) The plant has no poles on the unit circle. Only available information: P(e jθ ) for θ [0, 2π]

166 Data Based Design of 3 Term Controllers p. 72/10 Assumption The plant is stabilizable (mathematically, N(z) and D(z) are coprime). where P(z) = N(z) D(z) The plant has no poles on the unit circle. Only available information: P(e jθ ) for θ [0, 2π] Number of unstable poles of the plant, o p

167 Data Based Design of 3 Term Controllers p. 72/10 Assumption The plant is stabilizable (mathematically, N(z) and D(z) are coprime). where P(z) = N(z) D(z) The plant has no poles on the unit circle. Only available information: P(e jθ ) for θ [0, 2π] Number of unstable poles of the plant, o p The relative degree r

168 Data Based Design of 3 Term Controllers p. 73/10 PID Controller Design for Discrete-time Systems + C(z) P(z)

169 Data Based Design of 3 Term Controllers p. 73/10 PID Controller Design for Discrete-time Systems + C(z) P(z) Let the PID controller be of the form C(z) = K P + K ITz z 1 + K D(z 1) Tz

170 Data Based Design of 3 Term Controllers p. 73/10 PID Controller Design for Discrete-time Systems + C(z) P(z) Let the PID controller be of the form equivalently, C(z) = K P + K ITz z 1 + K D(z 1) Tz C(z) = K 2z 2 + K 1 z + K 0 z(z 1)

171 Data Based Design of 3 Term Controllers p. 73/10 PID Controller Design for Discrete-time Systems + C(z) P(z) Let the PID controller be of the form equivalently, C(z) = K P + K ITz z 1 + K D(z 1) Tz C(z) = K 2z 2 + K 1 z + K 0 z(z 1) where K P = K 1 2K 0, K I = K 0 + K 1 + K 2 T, K D = K 0 T

172 Data Based Design of 3 Term Controllers p. 74/10 PID Controller Design (continue...) The closed-loop characteristic polynomial is δ(z) := z(z 1)D(z) + `K 2 z 2 + K 1 z + K 0 N(z).

173 Data Based Design of 3 Term Controllers p. 74/10 PID Controller Design (continue...) The closed-loop characteristic polynomial is δ(z) := z(z 1)D(z) + `K 2 z 2 + K 1 z + K 0 N(z). Closed-loop stability requires σ(δ) = n + 2.

174 Data Based Design of 3 Term Controllers p. 74/10 PID Controller Design (continue...) The closed-loop characteristic polynomial is δ(z) := z(z 1)D(z) + `K 2 z 2 + K 1 z + K 0 N(z). Closed-loop stability requires σ(δ) = n + 2. Consider Π(z) = δ(z) D(z) = z(z 1) + `K 2 z 2 + K 1 z + K 0 P(z)

175 Data Based Design of 3 Term Controllers p. 74/10 PID Controller Design (continue...) The closed-loop characteristic polynomial is δ(z) := z(z 1)D(z) + `K 2 z 2 + K 1 z + K 0 N(z). Closed-loop stability requires σ(δ) = n + 2. Consider Π(z) = δ(z) D(z) = z(z 1) + `K 2 z 2 + K 1 z + K 0 P(z) Consequently, the stability requires that σ(π) = n + 2 i p where i p is the number of poles of the plant located inside unit circle.

176 Data Based Design of 3 Term Controllers p. 75/10 Stability Condition with PID Controllers Let P(z) be the plant with relative degree r. Let the PID controller be C(z) = K 2z 2 + K 1 z + K 0 z(z 1).

177 Data Based Design of 3 Term Controllers p. 75/10 Stability Condition with PID Controllers Let P(z) be the plant with relative degree r. Let the PID controller be C(z) = K 2z 2 + K 1 z + K 0 z(z 1) Then the closed-loop system is stable if and only if. σ(ν) = r + o z + 1 where ν(z) = z 1 P(z 1 )Π(z) Π(z) = z(z 1) + `K 2 z 2 + K 1 z + K 0 P(z) and o z is the number of non-minimum phase zeros of P(z).

178 Data Based Design of 3 Term Controllers p. 76/10 PID Controller Design (continue...) Now ν(z) z= u+jv = z 1 P(z 1 )Π(z) z= u+jv = ˆz 1 P(z 1 ) + `K 0 z 1 + K 1 + K 0 z P(z)P(z 1 ) z= u+jv = ( u 1 + jv) (R p (u) jvt p (u)) + [K 0 ( u jv) + K 1 + K 2 ( u + jv)] m 2 (u) where v = p 1 u 2, m(u) = P(z) z= u+jv.

179 Data Based Design of 3 Term Controllers p. 76/10 PID Controller Design (continue...) Now ν(z) z= u+jv = z 1 P(z 1 )Π(z) z= u+jv = ˆz 1 P(z 1 ) + `K 0 z 1 + K 1 + K 0 z P(z)P(z 1 ) z= u+jv = ( u 1 + jv) (R p (u) jvt p (u)) + [K 0 ( u jv) + K 1 + K 2 ( u + jv)] m 2 (u) where Then v = p 1 u 2, m(u) = P(z) z= u+jv. ν(z) = R ν (u, K 0, K 1, K 3 ) + jvt ν (u, K 3 ) where K 3 := K 2 K 0

180 Data Based Design of 3 Term Controllers p. 77/10 PID Controller Design (continue...) R ν (u, K 0, K 1, K 3 ) = (u 1)R p (u) (1 u 2 )T p (u) + K 1 m 2 (u) u(2k 0 + K 3 )m 2 (u) T ν (u, K 3 ) = R p (u) (u + 1)T p (u) + K 3 m 2 (u). Note that P(z) z= u+j 1 u 2 = R p(u) + j p 1 u 2 T p (u)

181 Data Based Design of 3 Term Controllers p. 77/10 PID Controller Design (continue...) R ν (u, K 0, K 1, K 3 ) = (u 1)R p (u) (1 u 2 )T p (u) + K 1 m 2 (u) u(2k 0 + K 3 )m 2 (u) T ν (u, K 3 ) = R p (u) (u + 1)T p (u) + K 3 m 2 (u). Note that P(z) z= u+j 1 u 2 = R p(u) + j p 1 u 2 T p (u) R p (u) and T P (u) can be obtained directly from the experimental data P(e jθ ) u= cos θ,

182 Data Based Design of 3 Term Controllers p. 77/10 PID Controller Design (continue...) R ν (u, K 0, K 1, K 3 ) = (u 1)R p (u) (1 u 2 )T p (u) + K 1 m 2 (u) u(2k 0 + K 3 )m 2 (u) T ν (u, K 3 ) = R p (u) (u + 1)T p (u) + K 3 m 2 (u). Note that P(z) z= u+j 1 u 2 = R p(u) + j p 1 u 2 T p (u) R p (u) and T P (u) can be obtained directly from the experimental data P(e jθ ) u= cos θ, R ν (u, K 0, K 1, K 3 ) and T ν (u, K 3 ) can also be obtained from P(e jθ )

183 Data Based Design of 3 Term Controllers p. 78/10 Algorithm: PID Controller Design For fixed K 3 = K 3, solve K 3 = R p(u) + (u + 1)T p (u) m 2 (u) = g(u) determine the roots u 1, u 2, u 0 = 1 < u 1 < u 2 < < u l < u l+1 = +1

184 Data Based Design of 3 Term Controllers p. 78/10 Algorithm: PID Controller Design For fixed K 3 = K 3, solve K 3 = R p(u) + (u + 1)T p (u) m 2 (u) = g(u) determine the roots u 1, u 2, u 0 = 1 < u 1 < u 2 < < u l < u l+1 = +1 Develop linear inequalities corresponding to stability as follows. Let I j = n o i j 0, ij 1,, ij l, ij l+1 denote a string where i j i {0,1, 1} such that i j 0 2ij 1 + 2ij 2 + ( 1)l+1 i j l+1 = r + o z + 1.

185 Data Based Design of 3 Term Controllers p. 79/10 Algorithm: PID Controller Design (continue...) For each string I j satisfying the above, we have the set of inequalities sgn [R ν (u t, K 0, K 1, K 3)] i j t > 0 which is a set of linear inequalities in K 0, K 1 space for fixed K 3.

186 Data Based Design of 3 Term Controllers p. 79/10 Algorithm: PID Controller Design (continue...) For each string I j satisfying the above, we have the set of inequalities sgn [R ν (u t, K 0, K 1, K 3)] i j t > 0 which is a set of linear inequalities in K 0, K 1 space for fixed K 3. By constructing these inequalities for each string satisfying the above, we obtain the stabilizing set for K 3 = K 3.

187 Data Based Design of 3 Term Controllers p. 79/10 Algorithm: PID Controller Design (continue...) For each string I j satisfying the above, we have the set of inequalities sgn [R ν (u t, K 0, K 1, K 3)] i j t > 0 which is a set of linear inequalities in K 0, K 1 space for fixed K 3. By constructing these inequalities for each string satisfying the above, we obtain the stabilizing set for K 3 = K 3. By sweeping over K 3 we can generate the complete set.

188 Data Based Design of 3 Term Controllers p. 79/10 Algorithm: PID Controller Design (continue...) For each string I j satisfying the above, we have the set of inequalities sgn [R ν (u t, K 0, K 1, K 3)] i j t > 0 which is a set of linear inequalities in K 0, K 1 space for fixed K 3. By constructing these inequalities for each string satisfying the above, we obtain the stabilizing set for K 3 = K 3. By sweeping over K 3 we can generate the complete set. The range of K 3 to be swept is determined by the requirement that g(u) should have o z+1 2 roots at least.

189 Data Based Design of 3 Term Controllers p. 80/10 An Example: PID Design Available Information: Frequency domain data P(e jωt ) := j P(e jωt ), ω = 2π T ff sampled every T = 0.01.

190 Data Based Design of 3 Term Controllers p. 80/10 An Example: PID Design Available Information: Frequency domain data P(e jωt ) := j P(e jωt ), ω = 2π T ff sampled every T = The plant is stable. In other words, the number of unstable poles of the plant is 0, that is, o p = 0.

191 Data Based Design of 3 Term Controllers p. 80/10 An Example: PID Design Available Information: Frequency domain data P(e jωt ) := j P(e jωt ), ω = 2π T ff sampled every T = The plant is stable. In other words, the number of unstable poles of the plant is 0, that is, o p = 0. The relative degree of the plant is 2, that is, r = 2.

192 Data Based Design of 3 Term Controllers p. 81/10 An Example: PID Design (continue...) The Nyquist plot of the plant P(z): 1.5 Nyquist Diagram Imaginary Axis Real Axis

193 Data Based Design of 3 Term Controllers p. 82/10 An Example: PID Design (continue...) Net phase: π 0 P(e jθ ) = π [r + (o z o p )] := 2π. o z = 2 r +o p = = 0.

194 Data Based Design of 3 Term Controllers p. 82/10 An Example: PID Design (continue...) Net phase: π 0 P(e jθ ) = π [r + (o z o p )] := 2π. o z = 2 r +o p = = 0. The stability requirement is equivalent to σ [ν(z)] = = 3.

195 Data Based Design of 3 Term Controllers p. 82/10 An Example: PID Design (continue...) Net phase: π 0 P(e jθ ) = π [r + (o z o p )] := 2π. o z = 2 r +o p = = 0. The stability requirement is equivalent to σ [ν(z)] = = 3. Applying this to the signature lemma,! 1 2 sgn[t( 1)] sgn[r( 1)] 2sgn[R(t 1)] + 2sgn[R(t 2 )] sgn[r(1)] := 3 where t i are the zeros of g(u) in g(u) for fixed K 3.

196 Data Based Design of 3 Term Controllers p. 83/10 An Example: PID Design (continue...) It is easy to see that at least two zeros t i are required and also that the only feasible string of sign sequences is: sgn of T( 1) R( 1) R(t 1 ) R(t 2 ) R(1)

197 Data Based Design of 3 Term Controllers p. 83/10 An Example: PID Design (continue...) It is easy to see that at least two zeros t i are required and also that the only feasible string of sign sequences is: sgn of T( 1) R( 1) R(t 1 ) R(t 2 ) R(1) The feasible range of K 3 values is that corresponding to the requirement of two zeros in T(u).

198 Data Based Design of 3 Term Controllers p. 83/10 An Example: PID Design (continue...) It is easy to see that at least two zeros t i are required and also that the only feasible string of sign sequences is: sgn of T( 1) R( 1) R(t 1 ) R(t 2 ) R(1) The feasible range of K 3 values is that corresponding to the requirement of two zeros in T(u). Plot the right hand side of g(u): g(u) =! 1 P c (u) 2 P r (u) (1 + u)p i (u) = K 3 where P(e jωt ) ωt=θ = P(e jθ ) u= cos θ = P r (u) + j p 1 u 2 P i (u).

199 Data Based Design of 3 Term Controllers p. 84/10 An Example: PID Design (continue...) fesible region of K 3 values 1 K u

200 Data Based Design of 3 Term Controllers p. 85/10 An Example: PID Design (continue...) At K 3 = 1.3, it is found from the graph that u = := t 1, u = := t 2.

201 Data Based Design of 3 Term Controllers p. 85/10 An Example: PID Design (continue...) At K 3 = 1.3, it is found from the graph that u = := t 1, u = := t 2. Then the set of linear inequalities corresponding to K 3 = 1.3 is T( 1) = 1 R( 1) = K K 2 > 0 R( ) = K K 2 < 0 R( ) = K K 2 > 0 R(1) = K K 2 < 0

202 Data Based Design of 3 Term Controllers p. 86/10 An Example: PID Design (continue...) By sweeping K 3 over ( 0.7,1.4), we have the stabilizing PID parameter regions shown: 2 1 K K K

203 Data Based Design of 3 Term Controllers p. 87/10 1st Order Controllers for Discrete-time Systems Consider the frequency response of the discrete-time plant P : P(z) z= u+jv = R p (u) + jvt p (u).

204 Data Based Design of 3 Term Controllers p. 87/10 1st Order Controllers for Discrete-time Systems Consider the frequency response of the discrete-time plant P : P(z) z= u+jv = R p (u) + jvt p (u). Note that R P (u) and T p (u) for 1 u 1 are immediately available from the given frequency response data points provided by P(e jθ ) for θ [0, π].

205 Data Based Design of 3 Term Controllers p. 87/10 1st Order Controllers for Discrete-time Systems Consider the frequency response of the discrete-time plant P : P(z) z= u+jv = R p (u) + jvt p (u). Note that R P (u) and T p (u) for 1 u 1 are immediately available from the given frequency response data points provided by P(e jθ ) for θ [0, π]. Consider the real rational function F(z) = (z + x 3 ) + (zx 1 + x 2 ) P(z).

206 Data Based Design of 3 Term Controllers p. 88/10 Stability Condition with First Order Controllers Let P(z) be the plant with the number of unstable poles being o p. Let the first order controller be C(z) = x 1z + x 2 z + x 3.

207 Data Based Design of 3 Term Controllers p. 88/10 Stability Condition with First Order Controllers Let P(z) be the plant with the number of unstable poles being o p. Let the first order controller be C(z) = x 1z + x 2 z + x 3. Then the closed-loop system is stable if and only if σ(π) = o p + 1 where Π(z) = (z + x 3 ) + (zx 1 + x 2 ) P(z) and o p is the number of non-minimum phase poles of P(z).

208 Data Based Design of 3 Term Controllers p. 89/10 1st Order Controllers Design Consider Π(z) = (z + x 3 ) + (zx 1 + x 2 ) P(z) and Π(z) z= u+jv = ( u + x 3 + jv) + ux 1 + x 2 + jvx 1! R p (u) + jvt p (u)! = (x 3 u) + R p (u)x 2 ur p (u) + v 2 T p (u)! x 1 +jv "! # R p (u) ut p (u) x 1 + T p (u)x

209 Data Based Design of 3 Term Controllers p. 90/10 Stability Conditions For complex root crossing, we now have the expression of the curve in (x 1, x 2 ) space for every fixed x ur p (u) + v 2 T p (u)! v R p (u) ut p (u)! R p (u) vt p (u) x 1(u) 7 5 x 2 (u) 3 5 = 2 4 (x 3 u) v 3 5. {z } A(u)

210 Data Based Design of 3 Term Controllers p. 90/10 Stability Conditions For complex root crossing, we now have the expression of the curve in (x 1, x 2 ) space for every fixed x ur p (u) + v 2 T p (u)! v R p (u) ut p (u)! R p (u) vt p (u) x 1(u) 7 5 x 2 (u) 3 5 = 2 4 (x 3 u) v 3 5. {z } A(u) Since det[a(u)] = `R p(u) 2 + v 2 Tp 2 (u) 2 v = v P(e jθ ),

211 Data Based Design of 3 Term Controllers p. 90/10 Stability Conditions For complex root crossing, we now have the expression of the curve in (x 1, x 2 ) space for every fixed x ur p (u) + v 2 T p (u)! v R p (u) ut p (u)! R p (u) vt p (u) x 1(u) 7 5 x 2 (u) 3 5 = 2 4 (x 3 u) v 3 5. {z } A(u) Since det[a(u)] = `R p(u) 2 + v 2 Tp 2 (u) 2 v = v P(e jθ ), the solution of the above is 2 4 x 1(u) x 2 (u) = P(e jθ ) (u x 3 ) T p (u) + R p (u) (1 ux 3 ) T p (u) + x 3 R p (u) 3 5.

212 Data Based Design of 3 Term Controllers p. 91/10 Stability Conditions (continue...) The two straight lines representing the real root crossing can be obtained by letting u = 1 and u = 1, equivalently letting θ = 0 and θ = π. (x 3 1) + P(e j0 )(x 2 x 1 ) = 0 (x 3 1) + P(e jπ )(x 2 x 1 ) = 0.

213 Data Based Design of 3 Term Controllers p. 92/10 An Example: First Order Controller Design Available Information: Frequency domain data P(e jωt ) := np(e jωt ), ω = π o T sampled every T = 0.01.

214 Data Based Design of 3 Term Controllers p. 92/10 An Example: First Order Controller Design Available Information: Frequency domain data P(e jωt ) := np(e jωt ), ω = π o T sampled every T = The plant is stable, i.e., the number of poles outside the unit circle is 0, that is, o p = 0.

215 Data Based Design of 3 Term Controllers p. 93/10 An Example: First Order Controller Design (continue...) At x = 0.75, the following region is obtained Region 4 Region x Region Region x 1

216 Data Based Design of 3 Term Controllers p. 93/10 An Example: First Order Controller Design (continue...) At x = 0.75, the following region is obtained Region 4 Region x Region Region x 1 Each separated region represents a set of controller parameters that gives a fixed number of unstable poles of the closed-loop system.

217 Data Based Design of 3 Term Controllers p. 93/10 An Example: First Order Controller Design (continue...) At x = 0.75, the following region is obtained Region 4 Region x Region Region x 1 Each separated region represents a set of controller parameters that gives a fixed number of unstable poles of the closed-loop system. To identify the stabilizing region, we arbitrarily select a point from each region and plot the corresponding Nyquist plot

218 Data Based Design of 3 Term Controllers p. 93/10 An Example: First Order Controller Design (continue...) At x = 0.75, the following region is obtained Region 4 Region x Region Region x 1 Each separated region represents a set of controller parameters that gives a fixed number of unstable poles of the closed-loop system. To identify the stabilizing region, we arbitrarily select a point from each region and plot the corresponding Nyquist plot

219 Data Based Design of 3 Term Controllers p. 94/10 An Example: First Order Controller Design (continue...) The following figure shows the Nyquist plots with selected controllers from the four specified regions. 6 A controller from Region A controller from Region 2 4 Imag 2 Imag Real A controller from Region Real A controller from Region 4 5 Imag Imag Real Real

220 Data Based Design of 3 Term Controllers p. 95/10 An Example: First Order Controller Design (continue...) The Nyquist plot with a controller from Region 1 shows that the encirclement around 1 point is 2π, i.e., N = 2.

221 Data Based Design of 3 Term Controllers p. 95/10 An Example: First Order Controller Design (continue...) The Nyquist plot with a controller from Region 1 shows that the encirclement around 1 point is 2π, i.e., N = 2. Since o p = 0, the corresponding closed-loop system will have 2 poles outside the unit circle.

222 Data Based Design of 3 Term Controllers p. 95/10 An Example: First Order Controller Design (continue...) The Nyquist plot with a controller from Region 1 shows that the encirclement around 1 point is 2π, i.e., N = 2. Since o p = 0, the corresponding closed-loop system will have 2 poles outside the unit circle. Similarly, corresponding closed-loops system with controllers from Region 3 and 4 will have 2 and 3 poles outside the unit circle, respectively.

223 Data Based Design of 3 Term Controllers p. 95/10 An Example: First Order Controller Design (continue...) The Nyquist plot with a controller from Region 1 shows that the encirclement around 1 point is 2π, i.e., N = 2. Since o p = 0, the corresponding closed-loop system will have 2 poles outside the unit circle. Similarly, corresponding closed-loops system with controllers from Region 3 and 4 will have 2 and 3 poles outside the unit circle, respectively. This test led us to the conclusion that the region 2 is the only stabilizing controller parameter region.

224 Data Based Design of 3 Term Controllers p. 96/10 An Example: First Order Controller Design (continue...) By sweeping over x 3, we have the entire first order stabilizing controllers for the given plant x x x 1

225 Data Based Design of 3 Term Controllers p. 97/10 Data Measurement by Impulse Response Typically data is obtained from the measurement by a sinusoidal input.

226 Data Based Design of 3 Term Controllers p. 97/10 Data Measurement by Impulse Response Typically data is obtained from the measurement by a sinusoidal input. Alternative input - Impulse Output - Markov parameters y[k] = [m 0, m 1, m 2,, m k, ]

227 Data Based Design of 3 Term Controllers p. 97/10 Data Measurement by Impulse Response Typically data is obtained from the measurement by a sinusoidal input. Alternative input - Impulse Output - Markov parameters y[k] = [m 0, m 1, m 2,, m k, ] Y (z) = m 0 + m 1 z m k z k +

228 Data Based Design of 3 Term Controllers p. 97/10 Data Measurement by Impulse Response Typically data is obtained from the measurement by a sinusoidal input. Alternative input - Impulse Output - Markov parameters y[k] = [m 0, m 1, m 2,, m k, ] Y (z) = m 0 + m 1 z m k z k + By truncating the terms, we have approximation of the frequency response of the system. P(z) z= u+jv = Y n (z) z= u+jv where n is the number of the terms taken from Y (z).

229 Data Based Design of 3 Term Controllers p. 98/10 Example: Data Measurement by Impulse Input Markov parameter obtained y[k] = [0, 0, 1, 0, 0.25, 0, , 0, , 0, , ]

230 Data Based Design of 3 Term Controllers p. 98/10 Example: Data Measurement by Impulse Input Markov parameter obtained y[k] = [0, 0, 1, 0, 0.25, 0, , 0, , 0, , ] Stabilizing region with K 3 = 1.2 for n = 3,5,7,10, n=3 K n= n=5 n=10, K 1

231 Data Based Design of 3 Term Controllers p. 99/10 Data Measurement by Step Input Alternative input - Step Input Step Response parameters y s [k] = [y 0, y 1, y 2,, y k, ]

232 Data Based Design of 3 Term Controllers p. 99/10 Data Measurement by Step Input Alternative input - Step Input Step Response parameters y s [k] = [y 0, y 1, y 2,, y k, ] =» z 1 H(z) z= u+jv = Y s (z) z z= u+jv y 0 + y 1 z 1 + y 2 z y k z k + `1 z 1 z= u+jv = y 0 + (y 1 y 0 )z 1 + (y 2 y 1 )z (y k y k 1 )z k + = m 0 + m 1 z 1 + m 2 z m k z k + where m 0, m 1,, m k, are the Markov parameters.

233 Data Based Design of 3 Term Controllers p. 99/10 Data Measurement by Step Input Alternative input - Step Input Step Response parameters y s [k] = [y 0, y 1, y 2,, y k, ] =» z 1 H(z) z= u+jv = Y s (z) z z= u+jv y 0 + y 1 z 1 + y 2 z y k z k + `1 z 1 z= u+jv = y 0 + (y 1 y 0 )z 1 + (y 2 y 1 )z (y k y k 1 )z k + = m 0 + m 1 z 1 + m 2 z m k z k + where m 0, m 1,, m k, are the Markov parameters. By truncating the terms, we have approximation of the frequency response of the system. P(z) z= u+jv = H n (z) z= u+jv where n is the number of the terms taken from H(z).

234 Data Based Design of 3 Term Controllers p. 100/10 Example: Data Measurement by Step Input The step response and the Markov parameters computed y s [k] = [0, 0, 1, 1, 1.25, 1.25, , , , , , m(k) = [0, 0, 1, 0, 0.25, 0, , , 0, , ]

235 Example: Data Measurement by Step Input The step response and the Markov parameters computed y s [k] = [0, 0, 1, 1, 1.25, 1.25, , , , , , m(k) = [0, 0, 1, 0, 0.25, 0, , , 0, , ] Stabilizing resion with K 3 = 1.2 for n = For K 3 = K K 1 Data Based Design of 3 Term Controllers p. 100/10

236 Data Based Design of 3 Term Controllers p. 101/10 Computer Aided Design using Labview The algorithm for the design of a PID Controller from the frequency response data of the system has programmed in LabVIEW due to its user-friendly graphical environment.

237 Data Based Design of 3 Term Controllers p. 101/10 Computer Aided Design using Labview The algorithm for the design of a PID Controller from the frequency response data of the system has programmed in LabVIEW due to its user-friendly graphical environment. The Virtual Instrument (VI) in Labview have a front panel that is displayed to the user and a block diagram, where the computations are performed.

238 Data Based Design of 3 Term Controllers p. 101/10 Computer Aided Design using Labview The algorithm for the design of a PID Controller from the frequency response data of the system has programmed in LabVIEW due to its user-friendly graphical environment. The Virtual Instrument (VI) in Labview have a front panel that is displayed to the user and a block diagram, where the computations are performed. The inputs to the program are the frequency response data and the number of RHP poles of the system.

239 Data Based Design of 3 Term Controllers p. 101/10 Computer Aided Design using Labview The algorithm for the design of a PID Controller from the frequency response data of the system has programmed in LabVIEW due to its user-friendly graphical environment. The Virtual Instrument (VI) in Labview have a front panel that is displayed to the user and a block diagram, where the computations are performed. The inputs to the program are the frequency response data and the number of RHP poles of the system. Given these inputs, the entire range of K p that can stabilize the system is displayed,

240 Data Based Design of 3 Term Controllers p. 101/10 Computer Aided Design using Labview The algorithm for the design of a PID Controller from the frequency response data of the system has programmed in LabVIEW due to its user-friendly graphical environment. The Virtual Instrument (VI) in Labview have a front panel that is displayed to the user and a block diagram, where the computations are performed. The inputs to the program are the frequency response data and the number of RHP poles of the system. Given these inputs, the entire range of K p that can stabilize the system is displayed, and as the user scrolls through the stabilizing range of K p, entire stabilizing ranges of K i and K d are displayed.

241 Data Based Design of 3 Term Controllers p. 102/10 CAD using LABVIEW(continue...) Left: The file containing the frequency response data from a stable system (number of right hand plane poles equals zero) is fed into the program through the file path box located at the top. Right: The front panel shows stabilizing sets of K p, K i and K d.

242 Data Based Design of 3 Term Controllers p. 103/10 CAD using LABVIEW(continue...) Left: T is set to a very small number.

243 Data Based Design of 3 Term Controllers p. 103/10 CAD using LABVIEW(continue...) Left: T is set to a very small number. g(ω) versus frequency generated from the inputs is used to compute ω values associated with a particular value of K p.

244 Data Based Design of 3 Term Controllers p. 103/10 CAD using LABVIEW(continue...) Left: T is set to a very small number. g(ω) versus frequency generated from the inputs is used to compute ω values associated with a particular value of K p. Using the values of ω obtained above, linear equations are solved by the program to compute stabilizing sets of K i and K d for a fixed K p.

245 Data Based Design of 3 Term Controllers p. 103/10 CAD using LABVIEW(continue...) Left: T is set to a very small number. g(ω) versus frequency generated from the inputs is used to compute ω values associated with a particular value of K p. Using the values of ω obtained above, linear equations are solved by the program to compute stabilizing sets of K i and K d for a fixed K p. As the selected K p is changed, the K i, K d graph changes dynamically to show the new stabilizing ranges of K i and K d for the selected K p.

246 Data Based Design of 3 Term Controllers p. 103/10 CAD using LABVIEW(continue...) Left: T is set to a very small number. g(ω) versus frequency generated from the inputs is used to compute ω values associated with a particular value of K p. Using the values of ω obtained above, linear equations are solved by the program to compute stabilizing sets of K i and K d for a fixed K p. As the selected K p is changed, the K i, K d graph changes dynamically to show the new stabilizing ranges of K i and K d for the selected K p. Right: Gives the stabilizing range of K i and K d for K p = 1.

247 Data Based Design of 3 Term Controllers p. 104/10 CAD using LABVIEW: 3D graph of stabilizing sets The stabilizing set of K p, K i and K d for the given system is shown.

248 Data Based Design of 3 Term Controllers p. 105/10 CAD using LABVIEW(continue...) Front Panel of VI showing performance indices for specific K p, K i and K d.

249 Data Based Design of 3 Term Controllers p. 106/10 CAD using LABVIEW(continue...) Front Panel of VI that satisfies multiple performance index specifications

250 Data Based Design of 3 Term Controllers p. 106/10 CAD using LABVIEW(continue...) Front Panel of VI that satisfies multiple performance index specifications The generated set of points satisfies multiple performance indices simultaneously. Performance indices used in the program: Gain Margin= 3db, Phase Margin= 45 o, Overshoot= 30%.

PID CONTROLLER SYNTHESIS FREE OF ANALYTICAL MODELS

PID CONTROLLER SYNTHESIS FREE OF ANALYTICAL MODELS PID CONTROLLER SYNTHESIS FREE OF ANALYTICAL MODELS L.H. Keel,1 S.P. Bhattacharyya, Tennessee State University, Nashville, Tennessee, USA Texas A&M University, College Station, Texas, USA Abstract: In this

More information

Analysis of SISO Control Loops

Analysis of SISO Control Loops Chapter 5 Analysis of SISO Control Loops Topics to be covered For a given controller and plant connected in feedback we ask and answer the following questions: Is the loop stable? What are the sensitivities

More information

Robust Control 3 The Closed Loop

Robust Control 3 The Closed Loop Robust Control 3 The Closed Loop Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /2/2002 Outline Closed Loop Transfer Functions Traditional Performance Measures Time

More information

COMPUTER AIDED SYNTHESIS AND DESIGN OF PID CONTROLLERS. A Thesis SANDIPAN MITRA

COMPUTER AIDED SYNTHESIS AND DESIGN OF PID CONTROLLERS. A Thesis SANDIPAN MITRA COMPUTER AIDED SYNTHESIS AND DESIGN OF PID CONTROLLERS A Thesis by SANDIPAN MITRA Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements for the

More information

DESIGN USING TRANSFORMATION TECHNIQUE CLASSICAL METHOD

DESIGN USING TRANSFORMATION TECHNIQUE CLASSICAL METHOD 206 Spring Semester ELEC733 Digital Control System LECTURE 7: DESIGN USING TRANSFORMATION TECHNIQUE CLASSICAL METHOD For a unit ramp input Tz Ez ( ) 2 ( z ) D( z) G( z) Tz e( ) lim( z) z 2 ( z ) D( z)

More information

STABILITY ANALYSIS TECHNIQUES

STABILITY ANALYSIS TECHNIQUES ECE4540/5540: Digital Control Systems 4 1 STABILITY ANALYSIS TECHNIQUES 41: Bilinear transformation Three main aspects to control-system design: 1 Stability, 2 Steady-state response, 3 Transient response

More information

STABILITY ANALYSIS. Asystemmaybe stable, neutrallyormarginallystable, or unstable. This can be illustrated using cones: Stable Neutral Unstable

STABILITY ANALYSIS. Asystemmaybe stable, neutrallyormarginallystable, or unstable. This can be illustrated using cones: Stable Neutral Unstable ECE4510/5510: Feedback Control Systems. 5 1 STABILITY ANALYSIS 5.1: Bounded-input bounded-output (BIBO) stability Asystemmaybe stable, neutrallyormarginallystable, or unstable. This can be illustrated

More information

Outline. Classical Control. Lecture 1

Outline. Classical Control. Lecture 1 Outline Outline Outline 1 Introduction 2 Prerequisites Block diagram for system modeling Modeling Mechanical Electrical Outline Introduction Background Basic Systems Models/Transfers functions 1 Introduction

More information

Frequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability

Frequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability Lecture 6. Loop analysis of feedback systems 1. Motivation 2. Graphical representation of frequency response: Bode and Nyquist curves 3. Nyquist stability theorem 4. Stability margins Frequency methods

More information

Lecture 5: Frequency domain analysis: Nyquist, Bode Diagrams, second order systems, system types

Lecture 5: Frequency domain analysis: Nyquist, Bode Diagrams, second order systems, system types Lecture 5: Frequency domain analysis: Nyquist, Bode Diagrams, second order systems, system types Venkata Sonti Department of Mechanical Engineering Indian Institute of Science Bangalore, India, 562 This

More information

Didier HENRION henrion

Didier HENRION  henrion GRADUATE COURSE ON POLYNOMIAL METHODS FOR ROBUST CONTROL PART II.1 ROBUST STABILITY ANALYSIS: POLYTOPIC UNCERTAINTY Didier HENRION www.laas.fr/ henrion henrion@laas.fr Cité de l espace with Ariane launcher

More information

The Complete Set of PID Controllers with Guaranteed Gain and Phase Margins

The Complete Set of PID Controllers with Guaranteed Gain and Phase Margins Proceedings of the th IEEE onference on Decision and ontrol, and the European ontrol onference 5 Seville, Spain, December -5, 5 ThA7.5 The omplete Set of PID ontrollers with Guaranteed Gain and Phase Margins

More information

ECE 486 Control Systems

ECE 486 Control Systems ECE 486 Control Systems Spring 208 Midterm #2 Information Issued: April 5, 208 Updated: April 8, 208 ˆ This document is an info sheet about the second exam of ECE 486, Spring 208. ˆ Please read the following

More information

ROOT LOCUS. Consider the system. Root locus presents the poles of the closed-loop system when the gain K changes from 0 to. H(s) H ( s) = ( s)

ROOT LOCUS. Consider the system. Root locus presents the poles of the closed-loop system when the gain K changes from 0 to. H(s) H ( s) = ( s) C1 ROOT LOCUS Consider the system R(s) E(s) C(s) + K G(s) - H(s) C(s) R(s) = K G(s) 1 + K G(s) H(s) Root locus presents the poles of the closed-loop system when the gain K changes from 0 to 1+ K G ( s)

More information

INTRODUCTION TO DIGITAL CONTROL

INTRODUCTION TO DIGITAL CONTROL ECE4540/5540: Digital Control Systems INTRODUCTION TO DIGITAL CONTROL.: Introduction In ECE450/ECE550 Feedback Control Systems, welearnedhow to make an analog controller D(s) to control a linear-time-invariant

More information

Homework 7 - Solutions

Homework 7 - Solutions Homework 7 - Solutions Note: This homework is worth a total of 48 points. 1. Compensators (9 points) For a unity feedback system given below, with G(s) = K s(s + 5)(s + 11) do the following: (c) Find the

More information

ECEN 605 LINEAR SYSTEMS. Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability 1/27

ECEN 605 LINEAR SYSTEMS. Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability 1/27 1/27 ECEN 605 LINEAR SYSTEMS Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability Feedback System Consider the feedback system u + G ol (s) y Figure 1: A unity feedback system

More information

ECE317 : Feedback and Control

ECE317 : Feedback and Control ECE317 : Feedback and Control Lecture : Routh-Hurwitz stability criterion Examples Dr. Richard Tymerski Dept. of Electrical and Computer Engineering Portland State University 1 Course roadmap Modeling

More information

MODERN PID CONTROLLER PARAMETER DESIGN FOR DC MOTOR

MODERN PID CONTROLLER PARAMETER DESIGN FOR DC MOTOR MODERN PID CONTROLLER PARAMETER DESIGN FOR DC MOTOR Thesis by ZHANGXIN ZHOU Submitted to the Office of Graduate and Professional Studies of Texas A&M University in partial fulfillment of the requirements

More information

Intro to Frequency Domain Design

Intro to Frequency Domain Design Intro to Frequency Domain Design MEM 355 Performance Enhancement of Dynamical Systems Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Closed Loop Transfer Functions

More information

Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.

Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition

More information

MEM 355 Performance Enhancement of Dynamical Systems

MEM 355 Performance Enhancement of Dynamical Systems MEM 355 Performance Enhancement of Dynamical Systems Frequency Domain Design Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 5/8/25 Outline Closed Loop Transfer Functions

More information

AN AUTOMATED VIRTUAL TOOL TO COMPUTE THE ENTIRE SET OF PROPORTIONAL INTEGRAL DERIVATIVE CONTROLLERS FOR A CONTINUOUS LINEAR TIME INVARIANT SYSTEM

AN AUTOMATED VIRTUAL TOOL TO COMPUTE THE ENTIRE SET OF PROPORTIONAL INTEGRAL DERIVATIVE CONTROLLERS FOR A CONTINUOUS LINEAR TIME INVARIANT SYSTEM AN AUTOMATED VIRTUAL TOOL TO COMPUTE THE ENTIRE SET OF PROPORTIONAL INTEGRAL DERIVATIVE CONTROLLERS FOR A CONTINUOUS LINEAR TIME INVARIANT SYSTEM A Thesis by BHARAT NARASIMHAN Submitted to the Office of

More information

Radar Dish. Armature controlled dc motor. Inside. θ r input. Outside. θ D output. θ m. Gearbox. Control Transmitter. Control. θ D.

Radar Dish. Armature controlled dc motor. Inside. θ r input. Outside. θ D output. θ m. Gearbox. Control Transmitter. Control. θ D. Radar Dish ME 304 CONTROL SYSTEMS Mechanical Engineering Department, Middle East Technical University Armature controlled dc motor Outside θ D output Inside θ r input r θ m Gearbox Control Transmitter

More information

Frequency domain analysis

Frequency domain analysis Automatic Control 2 Frequency domain analysis Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011

More information

Chapter 2. Classical Control System Design. Dutch Institute of Systems and Control

Chapter 2. Classical Control System Design. Dutch Institute of Systems and Control Chapter 2 Classical Control System Design Overview Ch. 2. 2. Classical control system design Introduction Introduction Steady-state Steady-state errors errors Type Type k k systems systems Integral Integral

More information

Raktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Frequency Response-Design Method

Raktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Frequency Response-Design Method .. AERO 422: Active Controls for Aerospace Vehicles Frequency Response- Method Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. ... Response to

More information

PID Controllers for S ystems Systems with with Time- Time Delay

PID Controllers for S ystems Systems with with Time- Time Delay PID Controllers for Systems with Time-Delay 1 General Considerations CHARACTERISTIC EQUATIONS FOR DEALY SYSTEMS u(t) Delay T y(t) =u(t T ) u(t) _y(t) Integrator y(t) _y(t)+ay(t T )=u(t) a Delay T 2 u(t)

More information

FEL3210 Multivariable Feedback Control

FEL3210 Multivariable Feedback Control FEL3210 Multivariable Feedback Control Lecture 5: Uncertainty and Robustness in SISO Systems [Ch.7-(8)] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 5:Uncertainty and Robustness () FEL3210 MIMO

More information

Plan of the Lecture. Goal: wrap up lead and lag control; start looking at frequency response as an alternative methodology for control systems design.

Plan of the Lecture. Goal: wrap up lead and lag control; start looking at frequency response as an alternative methodology for control systems design. Plan of the Lecture Review: design using Root Locus; dynamic compensation; PD and lead control Today s topic: PI and lag control; introduction to frequency-response design method Goal: wrap up lead and

More information

Introduction. Performance and Robustness (Chapter 1) Advanced Control Systems Spring / 31

Introduction. Performance and Robustness (Chapter 1) Advanced Control Systems Spring / 31 Introduction Classical Control Robust Control u(t) y(t) G u(t) G + y(t) G : nominal model G = G + : plant uncertainty Uncertainty sources : Structured : parametric uncertainty, multimodel uncertainty Unstructured

More information

Chapter 6 - Solved Problems

Chapter 6 - Solved Problems Chapter 6 - Solved Problems Solved Problem 6.. Contributed by - James Welsh, University of Newcastle, Australia. Find suitable values for the PID parameters using the Z-N tuning strategy for the nominal

More information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 23: Drawing The Nyquist Plot Overview In this Lecture, you will learn: Review of Nyquist Drawing the Nyquist Plot Using

More information

1 (20 pts) Nyquist Exercise

1 (20 pts) Nyquist Exercise EE C128 / ME134 Problem Set 6 Solution Fall 2011 1 (20 pts) Nyquist Exercise Consider a close loop system with unity feedback. For each G(s), hand sketch the Nyquist diagram, determine Z = P N, algebraically

More information

Discrete Systems. Step response and pole locations. Mark Cannon. Hilary Term Lecture

Discrete Systems. Step response and pole locations. Mark Cannon. Hilary Term Lecture Discrete Systems Mark Cannon Hilary Term 22 - Lecture 4 Step response and pole locations 4 - Review Definition of -transform: U() = Z{u k } = u k k k= Discrete transfer function: Y () U() = G() = Z{g k},

More information

Control Systems I. Lecture 9: The Nyquist condition

Control Systems I. Lecture 9: The Nyquist condition Control Systems I Lecture 9: The Nyquist condition Readings: Åstrom and Murray, Chapter 9.1 4 www.cds.caltech.edu/~murray/amwiki/index.php/first_edition Jacopo Tani Institute for Dynamic Systems and Control

More information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 21: Stability Margins and Closing the Loop Overview In this Lecture, you will learn: Closing the Loop Effect on Bode Plot Effect

More information

Control Systems I. Lecture 9: The Nyquist condition

Control Systems I. Lecture 9: The Nyquist condition Control Systems I Lecture 9: The Nyquist condition adings: Guzzella, Chapter 9.4 6 Åstrom and Murray, Chapter 9.1 4 www.cds.caltech.edu/~murray/amwiki/index.php/first_edition Emilio Frazzoli Institute

More information

Let the plant and controller be described as:-

Let the plant and controller be described as:- Summary of Fundamental Limitations in Feedback Design (LTI SISO Systems) From Chapter 6 of A FIRST GRADUATE COURSE IN FEEDBACK CONTROL By J. S. Freudenberg (Winter 2008) Prepared by: Hammad Munawar (Institute

More information

Robust fixed-order H Controller Design for Spectral Models by Convex Optimization

Robust fixed-order H Controller Design for Spectral Models by Convex Optimization Robust fixed-order H Controller Design for Spectral Models by Convex Optimization Alireza Karimi, Gorka Galdos and Roland Longchamp Abstract A new approach for robust fixed-order H controller design by

More information

Control Systems. Frequency Method Nyquist Analysis.

Control Systems. Frequency Method Nyquist Analysis. Frequency Method Nyquist Analysis chibum@seoultech.ac.kr Outline Polar plots Nyquist plots Factors of polar plots PolarNyquist Plots Polar plot: he locus of the magnitude of ω vs. the phase of ω on polar

More information

Synthesis of Controllers for Non-minimum Phase and Unstable Systems Using Non-sequential MIMO Quantitative Feedback Theory

Synthesis of Controllers for Non-minimum Phase and Unstable Systems Using Non-sequential MIMO Quantitative Feedback Theory Synthesis of Controllers for Non-minimum Phase and Unstable Systems Using Non-sequential MIMO Quantitative Feedback Theory Chen-yang Lan Murray L Kerr Suhada Jayasuriya * Graduate Student Graduate Student

More information

Definition of Stability

Definition of Stability Definition of Stability Transfer function of a linear time-invariant (LTI) system Fs () = b 2 1 0+ b1s+ b2s + + b m m m 1s - - + bms a0 + a1s+ a2s2 + + an-1sn- 1+ ansn Characteristic equation and poles

More information

The stability of linear time-invariant feedback systems

The stability of linear time-invariant feedback systems The stability of linear time-invariant feedbac systems A. Theory The system is atrictly stable if The degree of the numerator of H(s) (H(z)) the degree of the denominator of H(s) (H(z)) and/or The poles

More information

Frequency Response Techniques

Frequency Response Techniques 4th Edition T E N Frequency Response Techniques SOLUTION TO CASE STUDY CHALLENGE Antenna Control: Stability Design and Transient Performance First find the forward transfer function, G(s). Pot: K 1 = 10

More information

Automatic Control 2. Loop shaping. Prof. Alberto Bemporad. University of Trento. Academic year

Automatic Control 2. Loop shaping. Prof. Alberto Bemporad. University of Trento. Academic year Automatic Control 2 Loop shaping Prof. Alberto Bemporad University of Trento Academic year 21-211 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 21-211 1 / 39 Feedback

More information

D(s) G(s) A control system design definition

D(s) G(s) A control system design definition R E Compensation D(s) U Plant G(s) Y Figure 7. A control system design definition x x x 2 x 2 U 2 s s 7 2 Y Figure 7.2 A block diagram representing Eq. (7.) in control form z U 2 s z Y 4 z 2 s z 2 3 Figure

More information

MEM 355 Performance Enhancement of Dynamical Systems

MEM 355 Performance Enhancement of Dynamical Systems MEM 355 Performance Enhancement of Dynamical Systems Frequency Domain Design Intro Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /5/27 Outline Closed Loop Transfer

More information

Control Systems 2. Lecture 4: Sensitivity function limits. Roy Smith

Control Systems 2. Lecture 4: Sensitivity function limits. Roy Smith Control Systems 2 Lecture 4: Sensitivity function limits Roy Smith 2017-3-14 4.1 Input-output controllability Control design questions: 1. How well can the plant be controlled? 2. What control structure

More information

Control System Design

Control System Design ELEC ENG 4CL4: Control System Design Notes for Lecture #11 Wednesday, January 28, 2004 Dr. Ian C. Bruce Room: CRL-229 Phone ext.: 26984 Email: ibruce@mail.ece.mcmaster.ca Relative Stability: Stability

More information

Introduction to Feedback Control

Introduction to Feedback Control Introduction to Feedback Control Control System Design Why Control? Open-Loop vs Closed-Loop (Feedback) Why Use Feedback Control? Closed-Loop Control System Structure Elements of a Feedback Control System

More information

Nyquist Criterion For Stability of Closed Loop System

Nyquist Criterion For Stability of Closed Loop System Nyquist Criterion For Stability of Closed Loop System Prof. N. Puri ECE Department, Rutgers University Nyquist Theorem Given a closed loop system: r(t) + KG(s) = K N(s) c(t) H(s) = KG(s) +KG(s) = KN(s)

More information

Topic # Feedback Control. State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback

Topic # Feedback Control. State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback Topic #17 16.31 Feedback Control State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback Back to reality Copyright 21 by Jonathan How. All Rights reserved 1 Fall

More information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 22: The Nyquist Criterion Overview In this Lecture, you will learn: Complex Analysis The Argument Principle The Contour

More information

Robust Loop Shaping Controller Design for Spectral Models by Quadratic Programming

Robust Loop Shaping Controller Design for Spectral Models by Quadratic Programming Robust Loop Shaping Controller Design for Spectral Models by Quadratic Programming Gorka Galdos, Alireza Karimi and Roland Longchamp Abstract A quadratic programming approach is proposed to tune fixed-order

More information

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 9. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -

More information

Course Outline. Closed Loop Stability. Stability. Amme 3500 : System Dynamics & Control. Nyquist Stability. Dr. Dunant Halim

Course Outline. Closed Loop Stability. Stability. Amme 3500 : System Dynamics & Control. Nyquist Stability. Dr. Dunant Halim Amme 3 : System Dynamics & Control Nyquist Stability Dr. Dunant Halim Course Outline Week Date Content Assignment Notes 1 5 Mar Introduction 2 12 Mar Frequency Domain Modelling 3 19 Mar System Response

More information

MAS107 Control Theory Exam Solutions 2008

MAS107 Control Theory Exam Solutions 2008 MAS07 CONTROL THEORY. HOVLAND: EXAM SOLUTION 2008 MAS07 Control Theory Exam Solutions 2008 Geir Hovland, Mechatronics Group, Grimstad, Norway June 30, 2008 C. Repeat question B, but plot the phase curve

More information

Boise State University Department of Electrical Engineering ECE461 Control Systems. Control System Design in the Frequency Domain

Boise State University Department of Electrical Engineering ECE461 Control Systems. Control System Design in the Frequency Domain Boise State University Department of Electrical Engineering ECE6 Control Systems Control System Design in the Frequency Domain Situation: Consider the following block diagram of a type- servomechanism:

More information

Topic # Feedback Control

Topic # Feedback Control Topic #4 16.31 Feedback Control Stability in the Frequency Domain Nyquist Stability Theorem Examples Appendix (details) This is the basis of future robustness tests. Fall 2007 16.31 4 2 Frequency Stability

More information

Lecture 7 (Weeks 13-14)

Lecture 7 (Weeks 13-14) Lecture 7 (Weeks 13-14) Introduction to Multivariable Control (SP - Chapters 3 & 4) Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 7 (Weeks 13-14) p.

More information

Design and Tuning of Fractional-order PID Controllers for Time-delayed Processes

Design and Tuning of Fractional-order PID Controllers for Time-delayed Processes Design and Tuning of Fractional-order PID Controllers for Time-delayed Processes Emmanuel Edet Technology and Innovation Centre University of Strathclyde 99 George Street Glasgow, United Kingdom emmanuel.edet@strath.ac.uk

More information

Solutions to Skill-Assessment Exercises

Solutions to Skill-Assessment Exercises Solutions to Skill-Assessment Exercises To Accompany Control Systems Engineering 4 th Edition By Norman S. Nise John Wiley & Sons Copyright 2004 by John Wiley & Sons, Inc. All rights reserved. No part

More information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 23: Drawing The Nyquist Plot Overview In this Lecture, you will learn: Review of Nyquist Drawing the Nyquist Plot Using the

More information

Robust Performance Example #1

Robust Performance Example #1 Robust Performance Example # The transfer function for a nominal system (plant) is given, along with the transfer function for one extreme system. These two transfer functions define a family of plants

More information

ECSE 4962 Control Systems Design. A Brief Tutorial on Control Design

ECSE 4962 Control Systems Design. A Brief Tutorial on Control Design ECSE 4962 Control Systems Design A Brief Tutorial on Control Design Instructor: Professor John T. Wen TA: Ben Potsaid http://www.cat.rpi.edu/~wen/ecse4962s04/ Don t Wait Until The Last Minute! You got

More information

School of Mechanical Engineering Purdue University. DC Motor Position Control The block diagram for position control of the servo table is given by:

School of Mechanical Engineering Purdue University. DC Motor Position Control The block diagram for position control of the servo table is given by: Root Locus Motivation Sketching Root Locus Examples ME375 Root Locus - 1 Servo Table Example DC Motor Position Control The block diagram for position control of the servo table is given by: θ D 0.09 See

More information

Exercises for lectures 13 Design using frequency methods

Exercises for lectures 13 Design using frequency methods Exercises for lectures 13 Design using frequency methods Michael Šebek Automatic control 2016 31-3-17 Setting of the closed loop bandwidth At the transition frequency in the open loop is (from definition)

More information

Stability Margin Based Design of Multivariable Controllers

Stability Margin Based Design of Multivariable Controllers Stability Margin Based Design of Multivariable Controllers Iván D. Díaz-Rodríguez Sangjin Han Shankar P. Bhattacharyya Dept. of Electrical and Computer Engineering Texas A&M University College Station,

More information

ELECTRONICS & COMMUNICATIONS DEP. 3rd YEAR, 2010/2011 CONTROL ENGINEERING SHEET 5 Lead-Lag Compensation Techniques

ELECTRONICS & COMMUNICATIONS DEP. 3rd YEAR, 2010/2011 CONTROL ENGINEERING SHEET 5 Lead-Lag Compensation Techniques CAIRO UNIVERSITY FACULTY OF ENGINEERING ELECTRONICS & COMMUNICATIONS DEP. 3rd YEAR, 00/0 CONTROL ENGINEERING SHEET 5 Lead-Lag Compensation Techniques [] For the following system, Design a compensator such

More information

1 (s + 3)(s + 2)(s + a) G(s) = C(s) = K P + K I

1 (s + 3)(s + 2)(s + a) G(s) = C(s) = K P + K I MAE 43B Linear Control Prof. M. Krstic FINAL June 9, Problem. ( points) Consider a plant in feedback with the PI controller G(s) = (s + 3)(s + )(s + a) C(s) = K P + K I s. (a) (4 points) For a given constant

More information

Richiami di Controlli Automatici

Richiami di Controlli Automatici Richiami di Controlli Automatici Gianmaria De Tommasi 1 1 Università degli Studi di Napoli Federico II detommas@unina.it Ottobre 2012 Corsi AnsaldoBreda G. De Tommasi (UNINA) Richiami di Controlli Automatici

More information

Process Control & Instrumentation (CH 3040)

Process Control & Instrumentation (CH 3040) First-order systems Process Control & Instrumentation (CH 3040) Arun K. Tangirala Department of Chemical Engineering, IIT Madras January - April 010 Lectures: Mon, Tue, Wed, Fri Extra class: Thu A first-order

More information

Dr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review

Dr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review Week Date Content Notes 1 6 Mar Introduction 2 13 Mar Frequency Domain Modelling 3 20 Mar Transient Performance and the s-plane 4 27 Mar Block Diagrams Assign 1 Due 5 3 Apr Feedback System Characteristics

More information

Root Locus Methods. The root locus procedure

Root Locus Methods. The root locus procedure Root Locus Methods Design of a position control system using the root locus method Design of a phase lag compensator using the root locus method The root locus procedure To determine the value of the gain

More information

Control Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

Control Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich Control Systems I Lecture 7: Feedback and the Root Locus method Readings: Jacopo Tani Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 2, 2018 J. Tani, E. Frazzoli (ETH) Lecture 7:

More information

EECE 460 : Control System Design

EECE 460 : Control System Design EECE 460 : Control System Design SISO Pole Placement Guy A. Dumont UBC EECE January 2011 Guy A. Dumont (UBC EECE) EECE 460: Pole Placement January 2011 1 / 29 Contents 1 Preview 2 Polynomial Pole Placement

More information

Part II. Advanced PID Design Methods

Part II. Advanced PID Design Methods Part II Advanced PID Design Methods 54 Controller transfer function C(s) = k p (1 + 1 T i s + T d s) (71) Many extensions known to the basic design methods introduced in RT I. Four advanced approaches

More information

DIGITAL CONTROLLER DESIGN

DIGITAL CONTROLLER DESIGN ECE4540/5540: Digital Control Systems 5 DIGITAL CONTROLLER DESIGN 5.: Direct digital design: Steady-state accuracy We have spent quite a bit of time discussing digital hybrid system analysis, and some

More information

MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions

MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions Question 1. SIGNALS: Design of a noise-cancelling headphone system. 1a. Based on the low-pass filter given, design a high-pass filter,

More information

Theory of Machines and Automatic Control Winter 2018/2019

Theory of Machines and Automatic Control Winter 2018/2019 Theory of Machines and Automatic Control Winter 2018/2019 Lecturer: Sebastian Korczak, PhD, Eng. Institute of Machine Design Fundamentals - Department of Mechanics http://www.ipbm.simr.pw.edu.pl/ Lecture

More information

APPLICATIONS FOR ROBOTICS

APPLICATIONS FOR ROBOTICS Version: 1 CONTROL APPLICATIONS FOR ROBOTICS TEX d: Feb. 17, 214 PREVIEW We show that the transfer function and conditions of stability for linear systems can be studied using Laplace transforms. Table

More information

Step input, ramp input, parabolic input and impulse input signals. 2. What is the initial slope of a step response of a first order system?

Step input, ramp input, parabolic input and impulse input signals. 2. What is the initial slope of a step response of a first order system? IC6501 CONTROL SYSTEM UNIT-II TIME RESPONSE PART-A 1. What are the standard test signals employed for time domain studies?(or) List the standard test signals used in analysis of control systems? (April

More information

Raktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries

Raktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries . AERO 632: of Advance Flight Control System. Preliminaries Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Preliminaries Signals & Systems Laplace

More information

Lecture 9 Infinite Impulse Response Filters

Lecture 9 Infinite Impulse Response Filters Lecture 9 Infinite Impulse Response Filters Outline 9 Infinite Impulse Response Filters 9 First-Order Low-Pass Filter 93 IIR Filter Design 5 93 CT Butterworth filter design 5 93 Bilinear transform 7 9

More information

A system that is both linear and time-invariant is called linear time-invariant (LTI).

A system that is both linear and time-invariant is called linear time-invariant (LTI). The Cooper Union Department of Electrical Engineering ECE111 Signal Processing & Systems Analysis Lecture Notes: Time, Frequency & Transform Domains February 28, 2012 Signals & Systems Signals are mapped

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall 2007

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall 2007 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering.4 Dynamics and Control II Fall 7 Problem Set #9 Solution Posted: Sunday, Dec., 7. The.4 Tower system. The system parameters are

More information

CDS 101/110: Lecture 10.3 Final Exam Review

CDS 101/110: Lecture 10.3 Final Exam Review CDS 11/11: Lecture 1.3 Final Exam Review December 2, 216 Schedule: (1) Posted on the web Monday, Dec. 5 by noon. (2) Due Friday, Dec. 9, at 5: pm. (3) Determines 3% of your grade Instructions on Front

More information

Design Methods for Control Systems

Design Methods for Control Systems Design Methods for Control Systems Maarten Steinbuch TU/e Gjerrit Meinsma UT Dutch Institute of Systems and Control Winter term 2002-2003 Schedule November 25 MSt December 2 MSt Homework # 1 December 9

More information

Control of Manufacturing Processes

Control of Manufacturing Processes Control of Manufacturing Processes Subject 2.830 Spring 2004 Lecture #19 Position Control and Root Locus Analysis" April 22, 2004 The Position Servo Problem, reference position NC Control Robots Injection

More information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 24: Compensation in the Frequency Domain Overview In this Lecture, you will learn: Lead Compensators Performance Specs Altering

More information

FREQUENCY-RESPONSE ANALYSIS

FREQUENCY-RESPONSE ANALYSIS ECE450/550: Feedback Control Systems. 8 FREQUENCY-RESPONSE ANALYSIS 8.: Motivation to study frequency-response methods Advantages and disadvantages to root-locus design approach: ADVANTAGES: Good indicator

More information

CDS 101/110a: Lecture 10-1 Robust Performance

CDS 101/110a: Lecture 10-1 Robust Performance CDS 11/11a: Lecture 1-1 Robust Performance Richard M. Murray 1 December 28 Goals: Describe how to represent uncertainty in process dynamics Describe how to analyze a system in the presence of uncertainty

More information

Lecture 6 Classical Control Overview IV. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Lecture 6 Classical Control Overview IV. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lecture 6 Classical Control Overview IV Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lead Lag Compensator Design Dr. Radhakant Padhi Asst.

More information

Digital Control Systems

Digital Control Systems Digital Control Systems Lecture Summary #4 This summary discussed some graphical methods their use to determine the stability the stability margins of closed loop systems. A. Nyquist criterion Nyquist

More information

Stability and Robustness 1

Stability and Robustness 1 Lecture 2 Stability and Robustness This lecture discusses the role of stability in feedback design. The emphasis is notonyes/notestsforstability,butratheronhowtomeasurethedistanceto instability. The small

More information

Robust and Optimal Control, Spring A: SISO Feedback Control A.1 Internal Stability and Youla Parameterization

Robust and Optimal Control, Spring A: SISO Feedback Control A.1 Internal Stability and Youla Parameterization Robust and Optimal Control, Spring 2015 Instructor: Prof. Masayuki Fujita (S5-303B) A: SISO Feedback Control A.1 Internal Stability and Youla Parameterization A.2 Sensitivity and Feedback Performance A.3

More information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 12: Overview In this Lecture, you will learn: Review of Feedback Closing the Loop Pole Locations Changing the Gain

More information

Classify a transfer function to see which order or ramp it can follow and with which expected error.

Classify a transfer function to see which order or ramp it can follow and with which expected error. Dr. J. Tani, Prof. Dr. E. Frazzoli 5-059-00 Control Systems I (Autumn 208) Exercise Set 0 Topic: Specifications for Feedback Systems Discussion: 30.. 208 Learning objectives: The student can grizzi@ethz.ch,

More information

Time Response Analysis (Part II)

Time Response Analysis (Part II) Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary

More information