Exercise 5

HINT: For problems 1 and 2, use a generalized plant description and the feedback or LFT command in matlab.

Discuss Handin 3

Leftovers from Exercise 4:

Solve Problem 4 and Problem 6

Problem 1

Consider the following midranging architecture based on error feedback where

LaTeX: P_{1} = \frac{e^{-s}}{s + 1}P1=ess+1, and  LaTeX: P_2 = \frac{e^{-5s}}{10s + 1}P2=e5s10s+1. The controllers LaTeX: C_1 = \frac{1 + 2/(3s) + s/2}{s/2 + 1}C1=1+2/(3s)+s/2s/2+1 and LaTeX: C_2 = \frac{1 + 1/(10s) + 3s}{2s + 1}C2=1+1/(10s)+3s2s+1 have been tuned to give reasonable nice performance in feedback with LaTeX: P_1P1 and LaTeX: P_2P2 respectively.

.midrange.png

  1. Why do we use positive feedback with LaTeX: C_2C2?, i.e. why LaTeX: u_2 = C_2(u - u_r)u2=C2(uur) and not LaTeX: C_2(u_r - u)C2(uru)?
  2. Simulate the responses to LaTeX: yy, LaTeX: u_1u1 and LaTeX: u_2u2 for reference step changes to LaTeX: y_ryr and LaTeX: u_rur and input load disturbances
  3. Generate bode plots of the closed loop responses from LaTeX: (y_r, d_1, d_2, n) \mapsto (y, u_1, u_2)(yr,d1,d2,n)(y,u1,u2), where LaTeX: d_1d1 and LaTeX: d_2d2 are the corresponding input load disturbances and LaTeX: nn is measurement noise acting on LaTeX: yy. Generate the responses for LaTeX: k_\text{ff} = 0kff=0 and LaTeX: k_\text{ff} = -P_2(0) / P_1(0)kff=P2(0)/P1(0).

Problem 2

With LaTeX: P_1 = \frac{e^{-5s}}{10s + 1}P1=e5s10s+1 and LaTeX: P_{2} = \frac{e^{-s}}{s + 1}P2=ess+1, consider a Cascade architecture as in the block diagram below:

cascade-2.png

Design reasonable controllers for the cases

LaTeX: C_1C1 LaTeX: C_2C2
PD PID
PID PD
PID PID

 

  1. Discuss which closed-loop transfer functions are interesting, and plot their frequency responses. What's Gang of Four here?
  2. How would you implement anti-windup in these three cases. Introduce actuator limitations and simulate. To help you get started you can use the following simulink implementation: get_started_p2.zip Download get_started_p2.zip

Problem 3

Explain why the standard Smith predictor does not work for processes with integration or unstable dynamics

Problem 4

The standard Smith predictor for a process LaTeX: P(s) = P_0(s)e^{-sL}P(s)=P0(s)esL with time delay can be factored into

LaTeX: C(s) = C_0(s)C_\text{pred}(s), \qquad C_\text{pred}(s) = \frac{1}{1 + \hat P_0(s) C_0(s)(1-e^{-sL})}C(s)=C0(s)Cpred(s),Cpred(s)=11+ˆP0(s)C0(s)(1esL) where LaTeX: C_0C0 is the nominal controller for the process LaTeX: P_0P0 without delay and LaTeX: LL is the time delay. The transfer function LaTeX: C_\text{pred}(s)Cpred(s) is actually a good predictor that also can be used for loop shaping.

  1. Explore the Bode plot of the predictor LaTeX: C_\text{pred}(s)Cpred(s) for LaTeX: P_0(s) = \frac{1}{s + 1}, \quad C_0(s) = 2 + \frac{4}{s}P0(s)=1s+1,C0(s)=2+4s. Demonstrate that LaTeX: C_\text{pred}(s)Cpred(s) is indeed a good phase lead compensator
  2. The properties of the compensator LaTeX: C_\text{pred}(s)Cpred(s) changes qualitatively with LaTeX: LL, explain the nature of the changes and find the smallest value of LaTeX: LL for which the change occurs.
  3. Use the insight from B) to discuss fundamental limitations of the Smith predictor LaTeX: C_\text{pred}Cpred
  4. Smith's predictor can be used a as lead compensator without reference to time delays. Explore if it can be used as a lead compensator for the oscillatory system LaTeX: P(s) = \frac{1}{s^2 + 0.02s + 1}P(s)=1s2+0.02s+1

Problem 5

In the same set up as Problem 1 construct an integrating LQG controller LaTeX: \begin{bmatrix}u_1 \\ u_2 \end{bmatrix} = C\begin{bmatrix} y_1 \\ y_2 \end{bmatrix}[u1u2]=C[y1y2], where LaTeX: y_1 = y + n_1y1=y+n1 and LaTeX: y_2 = u_1 + n_2y2=u1+n2. Can you tune the controller to get midranging behavior?

Smith Predictor

The control structure is based on the idea to use a model LaTeX: \hat P_0(s)ˆP0(s)to predict the output of a delayed system, LaTeX: P(s) = P_0(s) e^{-sL}P(s)=P0(s)esL. The hat is introduced to differentiate between model and reality. otto_smith.png

The controller becomes LaTeX: C(s) = \frac{C_0}{1 + C_0 \hat P_0(1 - e^{-sL})}C(s)=C01+C0ˆP0(1esL). If LaTeX: \hat P_0 = P_0ˆP0=P0, the closed-loop becomes LaTeX: T = \frac{P_0 C_0}{1 + P_0 C_0} e^{-sL}T=P0C01+P0C0esL.