- Page 1 and 2: P.C. Chau © 2001 Table of Contents
- Page 3 and 4: P.C. Chau © 2001 10. Multiloop Sys
- Page 5 and 6: Stephanopoulos (1984), of collectin
- Page 7 and 8: 1 - 2 Desired pH + - Error pH Contr
- Page 9 and 10: P.C. Chau © 2001 ❖ 2. Mathematic
- Page 11 and 12: 2 - 3 2.2 Laplace transform Let us
- Page 13 and 14: 2 - 5 on complex variables, our Web
- Page 15 and 16: 2 - 7 2. Dead time function (Fig. 2
- Page 17 and 18: 2 - 9 lim s → 0 0 ∞ df(t) dt e
- Page 19 and 20: 2 - 11 α 3 = 6s 2 - 12 (s + 1) (s
- Page 21 and 22: 2 - 13 f(t) = 1 2 (1 - j) e( - 2+3j
- Page 23 and 24: 2 - 15 2.6 Transfer function, pole,
- Page 25: 2 - 17 final change in y(t) relativ
- Page 29 and 30: 2 - 21 At steady state, Eq. (2-35)
- Page 31 and 32: 2 - 23 We now raise a second questi
- Page 33 and 34: 2 - 25 where Laplace transform give
- Page 35 and 36: 2 - 27 steady state and reformulati
- Page 37 and 38: 2 - 29 V dC' dt + Q in,s C' = C in,
- Page 39 and 40: 2 - 31 Y R = G 1+GH (2-63) The RHS
- Page 41 and 42: 2 - 33 ✎ Example 2.16. Derive the
- Page 43 and 44: 2 - 35 We now do the expansion: y(t
- Page 45 and 46: P.C Chau © 2001 ❖ 3. Dynamic Res
- Page 47 and 48: 3 - 3 The Laplace transform of Eq.
- Page 49 and 50: 3 - 5 3.2 Second order differential
- Page 51 and 52: 3 - 7 (1) ζ > 1, overdamped. The r
- Page 53 and 54: 3 - 9 3.3 Processes with dead time
- Page 55 and 56: 3 - 11 q o c o c1 V 1 V 2 c n-1 q n
- Page 57 and 58: 3 - 13 approximation is all we use
- Page 59 and 60: 3 - 15 We do not need to carry the
- Page 61 and 62: 3 - 17 where (a) τ 1 is roughly th
- Page 63 and 64: 3 - 19 p2=conv([1 1],[3 1]); % Seco
- Page 65 and 66: 4-2 ✎ Example 4.1: Derive the sta
- Page 67 and 68: 4-4 With Example 4.1 as the hint, w
- Page 69 and 70: 4-6 now take the Laplace transform
- Page 71 and 72: 4-8 In this case where U and Y are
- Page 73 and 74: 4-10 ✎ Example 4.7A: Repeat Examp
- Page 75 and 76: 4-12 We first take the inlet glucos
- Page 77 and 78:
4-14 state space representation. 4.
- Page 79 and 80:
4-16 (also phase variable canonical
- Page 81 and 82:
4-18 The time domain solution vecto
- Page 83 and 84:
5 - 2 which is the negative of the
- Page 85 and 86:
5 - 4 5.1.2 Proportional-Integral (
- Page 87 and 88:
5 - 6 I action General qualitative
- Page 89 and 90:
5 - 8 point T sp (or reference R)
- Page 91 and 92:
5 - 10 5.2.3 Synthesis of a single-
- Page 93 and 94:
5 - 12 (2) We use a valve with line
- Page 95 and 96:
5 - 14 We now take a formal look at
- Page 97 and 98:
5 - 16 There are two noteworthy ite
- Page 99 and 100:
5 - 18 transmitted to a controller,
- Page 101 and 102:
5 - 20 point, we now want to reduce
- Page 103 and 104:
5 - 22 6. When we developed the mod
- Page 105 and 106:
6 - 2 Process L = 0 Step input P =
- Page 107 and 108:
6 - 4 ITAE = ∞ t e'(t) dt 0 (6-5)
- Page 109 and 110:
6 - 6 While the calculations in the
- Page 111 and 112:
6 - 8 For set point change: K c = a
- Page 113 and 114:
6 - 10 where τ c is the system tim
- Page 115 and 116:
6 - 12 guideline suggests that we n
- Page 117 and 118:
6 - 14 6.2.3 Internal model control
- Page 119 and 120:
6 - 16 Substitution of (E6-5) and (
- Page 121 and 122:
6 - 18 ❐ Review Problems 4. Repea
- Page 123 and 124:
Table 6.3. Summary of methods to se
- Page 125 and 126:
P.C. Chau © 2001 ❖ 7. Stability
- Page 127 and 128:
7 - 3 Once again, if all three pole
- Page 129 and 130:
7 - 5 The two additional constraint
- Page 131 and 132:
7 - 7 ✎ Example 7.2A: Apply direc
- Page 133 and 134:
7 - 9 function y=f(x) y = 5*x + tan
- Page 135 and 136:
7 - 11 parts of the controller func
- Page 137 and 138:
7 - 13 1+K 1 (s +3) (s + 2) (s + 1)
- Page 139 and 140:
7 - 15 7.5 Root Locus Design In ter
- Page 141 and 142:
7 - 17 out, and the response in thi
- Page 143 and 144:
8 - 2 y(t) = AK p τ p ω τ p 2 ω
- Page 145 and 146:
8 - 4 transfer function can be "bro
- Page 147 and 148:
8 - 6 other hand, is easier to inte
- Page 149 and 150:
8 - 8 is obvious that the 70.7% com
- Page 151 and 152:
8 - 10 figure(1), bode(G); figure(2
- Page 153 and 154:
8 - 12 polar(phase,mag) ✎ Example
- Page 155 and 156:
P.C. Chau © 2001 8 - 14 8.3 Stabil
- Page 157 and 158:
8 - 16 margin between 30° and 45°
- Page 159 and 160:
8 - 18 ✎ Example 8.12. Derive the
- Page 161 and 162:
8 - 20 for a PI controller. From a
- Page 163 and 164:
8 - 22 closed-loop system apply onl
- Page 165 and 166:
8 - 24 underdamped time response cu
- Page 167 and 168:
8 - 26 angle and the Nyquist criter
- Page 169 and 170:
8 - 28 taui=3; taud=0.5; gc=tf([tau
- Page 171 and 172:
P.C. Chau © 2001 ❖ 9. Design of
- Page 173 and 174:
9 - 3 With the same reasoning that
- Page 175 and 176:
9 - 5 9.2 Pole Placement Design ✑
- Page 177 and 178:
9 - 7 where K r is some gain associ
- Page 179 and 180:
9 - 9 To evaluate the matrix polyno
- Page 181 and 182:
9 - 11 A note of caution is necessa
- Page 183 and 184:
9 - 13 tools of pole-placement for
- Page 185 and 186:
9 - 15 x e = (A ee - K er A 1e )x ~
- Page 187 and 188:
9 - 17 x = Fx ~ + Gu + Hy Find Eq.
- Page 189 and 190:
P.C. Chau © 2001 ❖ 10. Multiloop
- Page 191 and 192:
10-3 The remaining task to derive t
- Page 193 and 194:
10-5 τ = 0.05 s I τ = 0.5 s I τ
- Page 195 and 196:
10-7 expectation that we'll introdu
- Page 197 and 198:
10-9 10.3 Feedforward-feedback cont
- Page 199 and 200:
10-11 A more sophisticated implemen
- Page 201 and 202:
10-13 10.6 Multiple-input Multiple-
- Page 203 and 204:
10-15 If both references change sim
- Page 205 and 206:
10-17 ✑ 10.6.3 Relative gain arra
- Page 207 and 208:
10-19 other loops is in opposition
- Page 209 and 210:
10-21 G 11 G 12 G 21 G 22 d 11 d 12
- Page 211 and 212:
10-23 K = 0.07 - 0.05 0.1 - 0.15 Wi
- Page 213 and 214:
10-25 µ 2 = m 1 Find the gain matr
- Page 215 and 216:
10-27 expansion as we have done wit
- Page 217 and 218:
M1 - 2 It is important to know that
- Page 219 and 220:
M1 - 4 Of course, we can solve the
- Page 221 and 222:
So you say wow! But MATLAB can do m
- Page 223 and 224:
P.C. Chau © 2000 MATLAB Session 2
- Page 225 and 226:
F(s) = 6s 2 -12 6(s- 2)(s+ 2) (s 3
- Page 227 and 228:
M2 - 5 This example is simple enoug
- Page 229 and 230:
M3 - 2 The functions also handle mu
- Page 231 and 232:
M3.2 LTI Viewer M3 - 4 Features cov
- Page 233 and 234:
M4 - 2 From X 2 U = 1 (s 2 +2ζω n
- Page 235 and 236:
M4 - 4 If we examine the values of
- Page 237 and 238:
poly(A) % Check the characteristic
- Page 239 and 240:
P.C. Chau © 2001 MATLAB Session 5
- Page 241 and 242:
M5 - 3 M5.2 Control toolbox functio
- Page 243 and 244:
eig(ssm.a) M5 - 5 For fun, we can r
- Page 245 and 246:
M6 - 2 The point of the last two ca
- Page 247 and 248:
M6 - 4 A graphics window with pull-
- Page 249 and 250:
op_zero=[-0.5 -1.8]; % one zero in
- Page 251 and 252:
P.C. Chau © 2001 MATLAB Session 7
- Page 253 and 254:
G=tf(1,[1 0.4 1]); bode(G) % Done!
- Page 255:
M7 - 5 There is no magic in the fun