Data Based Design of 3Term Controllers. Data Based Design of 3 Term Controllers p. 1/10
|
|
- Cynthia Strickland
- 5 years ago
- Views:
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 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 informationAnalysis 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 informationRobust 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 informationCOMPUTER 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 informationDESIGN 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 informationSTABILITY 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 informationSTABILITY 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 informationOutline. 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 informationFrequency 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 informationLecture 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 informationDidier 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 informationThe 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 informationECE 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 informationROOT 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 informationINTRODUCTION 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 informationHomework 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 informationECEN 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 informationECE317 : 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 informationMODERN 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 informationIntro 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 informationLecture 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 informationMEM 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 informationAN 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 informationRadar 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 informationFrequency 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 informationChapter 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 informationRaktim 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 informationPID 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 informationFEL3210 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 informationPlan 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 informationIntroduction. 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 informationChapter 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 informationSystems 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 information1 (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 informationDiscrete 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 informationControl 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 informationSystems 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 informationControl 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 informationLet 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 informationRobust 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 informationControl 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 informationSynthesis 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 informationDefinition 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 informationThe 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 informationFrequency 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 informationAutomatic 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 informationD(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 informationMEM 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 informationControl 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 informationControl 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 informationIntroduction 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 informationNyquist 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 informationTopic # 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 informationSystems 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 informationRobust 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!
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 informationCourse 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 informationMAS107 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 informationBoise 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 informationTopic # 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 informationLecture 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 informationDesign 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 informationSolutions 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 informationSystems 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 informationRobust 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 informationECSE 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 informationSchool 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 informationExercises 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 informationStability 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 informationELECTRONICS & 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 information1 (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 informationRichiami 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 informationProcess 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 informationDr 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 informationRoot 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 informationControl 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 informationEECE 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 informationPart 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 informationDIGITAL 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 informationMAE143a: 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 informationTheory 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 informationAPPLICATIONS 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 informationStep 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 informationRaktim 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 informationLecture 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 informationA 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 informationMASSACHUSETTS 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 informationCDS 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 informationDesign 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 informationControl 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 informationSystems 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 informationFREQUENCY-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 informationCDS 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 informationLecture 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 informationDigital 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 informationStability 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 informationRobust 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 informationSystems 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 informationClassify 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 informationTime 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