Board Reference

    Flight Controller
    Field Oriented Controller

    Motor Drive

    PMSM

    Motion Control

    Attitude Control

    Tutorial

    Mathematics
    Classical Control
    State Estimation
    Feedforward Control
    Disturbance Observer
    Motion Control

Robustness

Problem statement

A problem is to control an output y\bm{y} of following system using an input u\bm{u}.
y=Pu+d,\begin{align} \bm{y} &= \bm{P}\bm{u} + \bm{d}, \end{align}
where d\bm{d} and d\bm{d} denote a plant model and disturbance. An input should be designed so that it make an output track a desired trajectory and a disturbance be suppressed. For this purpose, feedforward and feedback control should be introduced. Now, let us consider the nonideality of input and observation:
u=uref+vz=y+w,\begin{align} \bm{u} &= \bm{u}^{\rm ref} + \bm{v} \\ \bm{z} &= \bm{y} + \bm{w}, \end{align}
where z\bm{z}, v\bm{v}, and w\bm{w} stand for a sensor output, process noise, and sensor noise. A feedback controller is constructed using a sensor output.

Simple feedback control

1-DOF control systems enable the design of one characteristic of the control system. The variables r\bm{r}, y\bm{y}, d\bm{d}, v\bm{v}, and w\bm{w} are the reference, the output, the disturbance, the process noise, and the observation noise. The noises v\bm{v} and w\bm{w} are the Gaussian noises, which are seeds of the noises that superimpose into the system. The blocks C\bm{C}, P\bm{P}, H1\bm{H}_1, and H2\bm{H}_2 denote the feedback controller, the plant, and the noise-shaping filters for the process noise and the observation noise, respectively. These noise-shaping filters determine the colors and magnitudes of the noises. The term ΔP\bm{\Delta}_P represents the fluctuation of the plant dynamics due to uncertainty or perturbation on the system. The multiplicative form P(I+ΔP)\bm{P}(\bm{I}+\bm{\Delta}_P) is useful in introducing a fluctuation rate into the controller design. For the simplicity of discussion, let the entire plant dynamics P(I+ΔP)\bm{P}(\bm{I}+\bm{\Delta}_P) be expressed as Pm\bm{P}_{\rm m}. Herein, the output of this controller is generated as
y=(I+PmC)1PmCr(I+PmC)1Pmd+(I+PmC)1PmH1v(I+PmC)1PmCH2w,\begin{align} \bm{y}=&(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}\bm{P}_{\rm m}\bm{C}\bm{r}-(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}\bm{P}_{\rm m}\bm{d} \notag\\ &+(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}\bm{P}_{\rm m}\bm{H}_1\bm{v}-(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}\bm{P}_{\rm m}\bm{CH}_2\bm{w}, \end{align}
where I\bm{I} is the unit matrix. The nominal output can be checked by replacing Pm\bm{P}_{\rm m} with P\bm{P}. In the controller design, the uncertainty is once removed from the plant model, and the nominal model P\bm{P} is mainly used. Herein, the condition that the fluctuation of the plant system does not destabilize the control system is provided by the small gain theorem;
(I+PC)1PCΔP1,\begin{align} \| (\bm{I}+\bm{PC})^{-1}\bm{PC}\bm{\Delta}_P \|_{\infty} \leq 1, \end{align}
where \|\cdot\|_{\infty} denotes the HH_{\infty} norm. In the equation above, the two coefficients are standing out;
S=(I+PC)1T=(I+PC)1PC,\begin{align} \bm{S}&=(\bm{I}+\bm{PC})^{-1}\\ \bm{T}&=(\bm{I}+\bm{PC})^{-1}\bm{PC}, \end{align}
where S\bm{S} is the sensitivity function, and T\bm{T} is the complementary-sensitivity function. As their names denote, these are in the complementary relationship, which is expressed as
S+T=I.\begin{align} \bm{S}+\bm{T}=\bm{I}. \end{align}
The sensitivity means susceptibility of the transfer function from the command to the output. Here, let the transfer function from the command to the output be represented as
Try(P)=(I+PC)1PC.\begin{align} \bm{T}_{ry}(\bm{P})=(\bm{I}+\bm{P}\bm{C})^{-1}\bm{P}\bm{C}. \end{align}
Herein, Try\bm{T}_{ry} is described as the function of P\bm{P}. When the plant system fluctuates from P\bm{P} to Pm\bm{P}_{\rm m}, the rates of the changes of the plant system RP\bm{R}_{P} and the transfer function RT\bm{R}_{T} can be described as
RP=(PmP)Pm1RT=(Try(Pm)Try(P))Try1(Pm).\begin{align} \bm{R}_{P}&=(\bm{P}_{\rm m}-\bm{P})\bm{P}_{\rm m}^{-1}\\ \bm{R}_{T}&=\left(\bm{T}_{ry}(\bm{P}_{\rm m})-\bm{T}_{ry}(\bm{P})\right)\bm{T}_{ry}^{-1}(\bm{P}_{\rm m}). \end{align}
With the simple algebraic operations, we obtain
RT=ITry(P)Try1(Pm)=I(I+PC)1PCC1Pm1(I+PmC)=(I+PC)1(IPPm1),\begin{align} \bm{R}_{T}&=\bm{I}-\bm{T}_{ry}(\bm{P})\bm{T}_{ry}^{-1}(\bm{P}_{\rm m}) \notag\\ &=\bm{I}-(\bm{I}+\bm{P}\bm{C})^{-1}\bm{P}\bm{C}\bm{C}^{-1}\bm{P}_{\rm m}^{-1}(\bm{I}+\bm{P}_{\rm m}\bm{C}) \notag\\ &=(\bm{I}+\bm{PC})^{-1}(\bm{I}-\bm{P}\bm{P}_{\rm m}^{-1}), \end{align}
and the following relationship can be found
RT=(I+PC)1RP=SRp.\begin{align} \bm{R}_{T}&=(\bm{I}+\bm{PC})^{-1}\bm{R}_{P}=\bm{S}\bm{R}_{\rm p}. \end{align}
Thus, the sensitivity function shows the robustness of the tracking performance against the fluctuation. On the other hand, the complementary-sensitivity function indicates the robustness of the stability against the fluctuation, as denoted in \eq{math:chap2_RobustStable}. These are called robust performance and robust stability. Besides, the index that represents the effect of the disturbance on the responses is called the quasi-sensitivity. The quasi-sensitivity function is defined as
Sq=SP.\begin{align} \bm{S}_{\rm q}&=\bm{S}\bm{P}. \end{align}
Using these functions, the output of the nominal system can be rewritten as
y=TrSqd+SqH1vTH2w.\begin{align} \bm{y}=&\bm{T}\bm{r}-\bm{S}_{\rm q}\bm{d}+\bm{S}_{\rm q}\bm{H}_1\bm{v}-\bm{T}\bm{H}_2\bm{w}. \end{align}
This equation shows five tradeoffs between
  1. tracking performance and disturbance suppression performance;
  2. tracking performance and noise sensitivity;
  3. tracking performance and robust stability;
  4. disturbance suppression performance and noise sensitivity;
  5. disturbance suppression performance and robust stability.
The 1-DOF controller can handle only one characteristic due to the lack of DOFs. The first one is a well-known problem referred to as the pole-zero assignment problem. On the 1-DOF controller, the transfer function from the command to the response Try\bm{T}_{ry} and that from the disturbance to the response Tdy\bm{T}_{dy} have the same poles. Following the internal model principle \cite{IMC_Francis_1976}, most controllers have an integrator to remove a steady-state error, but this integrator puts an unexpected zero on Try\bm{T}_{ry}. Adding poles to compensate the disturbance increases the order of the system, and generates unexpected zeros. Therefore, a conservative design that sacrifices one side imposes on engineers. It should be noted that the added zeros and poles also influence the system noise. Furthermore, the controller should have a high gain so that both a numerator and a denominator of Try\bm{T}_{ry} get almost the same value to improve the tracking performance. This approach causes to amplify the system noises and makes the response oscillatory. For those factors, it is hard to elicit the performance when using the 1-DOF controller.

2-DOF Control

The 1-DOF control architecture has been revised. Here, the blocks K\bm{K}, Nn\bm{N}_{\rm n}, and Dn\bm{D}_{\rm n} stand for the tracking controller and the irreducible numerator and denominator of the left coprime factorization of a nominal plant modelPn\bm{P}_{\rm n}, i.e., Pn=NnDn1\bm{P}_{\rm n}=\bm{N}_{\rm n}\bm{D}^{-1}_{\rm n}. This architecture enables the design of the tracking performance at will by introducing a model-based feedforward loop. Since the 1-DOF controller generates a control effort for tracking from trajectory error, a phase lag between the command and the response is inevitable. Furthermore, the tracking performance is affected by the controller since the control effort is the output of the controller. The solution which the architecture shows is to generate the control effort required for command tracking in a designed loop. This controller is referred to as the 2-DOF controller, which provides two available characteristics, which are tracking performance and disturbance suppression performance. In other words, the 2-DOF controller includes a feedback controller and a feedback controller. The feedforward controller determines the behavior of a plant system when a plant system is not disturbed from an external system. The feedback controller compensates tracking error which is caused by disturbance. Herein, the output is generated as
y=(I+PmC)1(PmPn1+PmC)NnKr(I+PmC)1Pmd+(I+PmC)1PmH1v(I+PmC)1PmCH2w=NnKr+(I+PmC)1ΔPNnKr(I+PmC)1Pmd+(I+PmC)1PmH1v(I+PmC)1PmCH2w.\begin{align} \bm{y}=&(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}(\bm{P}_{\rm m}\bm{P}_{\rm n}^{-1}+\bm{P}_{\rm m}\bm{C})\bm{N}_{\rm n}\bm{K}\bm{r}-(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}\bm{P}_{\rm m}\bm{d} \notag\\ &+(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}\bm{P}_{\rm m}\bm{H}_1\bm{v}-(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}\bm{P}_{\rm m}\bm{CH}_2\bm{w}\\ =&\bm{N}_{\rm n}\bm{K}\bm{r}+(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}\bm{\Delta}_{P}\bm{N}_{\rm n}\bm{K}\bm{r}-(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}\bm{P}_{\rm m}\bm{d} \notag\\ &+(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}\bm{P}_{\rm m}\bm{H}_1\bm{v}-(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}\bm{P}_{\rm m}\bm{CH}_2\bm{w}. \end{align}
Because tracking performance relies on model-based feedforward control, the mismatch of a model causes an error. The additional factor is only the feedforward loop, and hence, the stability against plant fluctuation T\bm{T} is the same as that of the 1-DOF control system. Regarding to the sensitivity of the system response against plant fluctuation, the following conditions are found
RP=(PmP)Pm1RT=ITry(P)Try1(Pm)=I(I+PC)1P(C+Pn1)NnKK1Nn1(C+Pn1)1Pm1(I+PmC)=(I+PC)1(IPPm1).\begin{align} \bm{R}_{P}&=(\bm{P}_{\rm m}-\bm{P})\bm{P}_{\rm m}^{-1}\\ \bm{R}_{T}&=\bm{I}-\bm{T}_{ry}(\bm{P})\bm{T}_{ry}^{-1}(\bm{P}_{\rm m}) \notag\\ &=\bm{I}-(\bm{I}+\bm{P}\bm{C})^{-1}\bm{P}\left(\bm{C}+\bm{P}_{\rm n}^{-1}\right)\bm{N}_{\rm n}\bm{K}\bm{K}^{-1}\bm{N}_{\rm n}^{-1}\left(\bm{C}+\bm{P}_{\rm n}^{-1}\right)^{-1}\bm{P}_{\rm m}^{-1}(\bm{I}+\bm{P}_{\rm m}\bm{C}) \notag\\ &=(\bm{I}+\bm{PC})^{-1}(\bm{I}-\bm{P}\bm{P}_{\rm m}^{-1}). \end{align}
Therefore, the sensitivity, the quasi-sensitivity, and the complementary-sensitivity function are the same as those of the 1-DOF controller. When the system model well describes the plant system or when the plant fluctuation is small, the output is rewritten as
y=NnKr(I+PmC)1Pmd+(I+PmC)1PmH1v(I+PmC)1PmCH2w=NnKrSqd+SqH1vTH2w.\begin{align} \bm{y}=&\bm{N}_{\rm n}\bm{K}\bm{r}-(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}\bm{P}_{\rm m}\bm{d}+(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}\bm{P}_{\rm m}\bm{H}_1\bm{v}-(\bm{I}+\bm{P}_{\rm m}\bm{C})^{-1}\bm{P}_{\rm m}\bm{CH}_2\bm{w} \notag\\ =&\bm{N}_{\rm n}\bm{K}\bm{r}-\bm{S}_{\rm q}\bm{d}+\bm{S}_{\rm q}\bm{H}_1\bm{v}-\bm{T}\bm{H}_2\bm{w}. \end{align}
In this equation, the coefficient of r\bm{r} does not have the characteristics of the controller C\bm{C}. It means that the tradeoffs relating to the tracking performance, including the pole-zero assignment problem, are removed. Therefore, the remaining problems are the two tradeoffs between
  1. disturbance suppression performance and noise sensitivity;
  2. disturbance suppression performance and robust stability.
It should be noted that the noise characteristics and the plant characteristics are not something that can be shaped by modifying a control architecture, and a design of the hardware is critical to reducing the constraints above. These should be handled by the integrated design of the control system, with the theory of control, power electronics, and computer science.

© DigitalServo