H-infinity methods in control theory

H (i.e. "H-infinity") methods are used in control theory to synthesize controllers achieving stabilization with guaranteed performance. To use H methods, a control designer expresses the control problem as a mathematical optimization problem and then finds the controller that solves this optimization. H techniques have the advantage over classical control techniques in that they are readily applicable to problems involving multivariate systems with cross-coupling between channels; disadvantages of H techniques include the level of mathematical understanding needed to apply them successfully and the need for a reasonably good model of the system to be controlled. It is important to keep in mind that the resulting controller is only optimal with respect to the prescribed cost function and does not necessarily represent the best controller in terms of the usual performance measures used to evaluate controllers such as settling time, energy expended, etc. Also, non-linear constraints such as saturation are generally not well-handled. These methods were introduced into control theory in the late 1970s-early 1980s by George Zames (sensitivity minimization),[1] J. William Helton (broadband matching),[2] and Allen Tannenbaum (gain margin optimization).[3]

The phrase H control comes from the name of the mathematical space over which the optimization takes place: H is the space of matrix-valued functions that are analytic and bounded in the open right-half of the complex plane defined by Re(s) > 0; the H norm is the maximum singular value of the function over that space. (This can be interpreted as a maximum gain in any direction and at any frequency; for SISO systems, this is effectively the maximum magnitude of the frequency response.) H techniques can be used to minimize the closed loop impact of a perturbation: depending on the problem formulation, the impact will either be measured in terms of stabilization or performance.

Simultaneously optimizing robust performance and robust stabilization is difficult. One method that comes close to achieving this is H loop-shaping, which allows the control designer to apply classical loop-shaping concepts to the multivariable frequency response to get good robust performance, and then optimizes the response near the system bandwidth to achieve good robust stabilization.

Commercial software is available to support H controller synthesis.

Problem formulation

First, the process has to be represented according to the following standard configuration:

The plant P has two inputs, the exogenous input w, that includes reference signal and disturbances, and the manipulated variables u. There are two outputs, the error signals z that we want to minimize, and the measured variables v, that we use to control the system. v is used in K to calculate the manipulated variable u. Notice that all these are generally vectors, whereas P and K are matrices.

In formulae, the system is:

\begin{bmatrix} z\\ v \end{bmatrix} = \mathbf{P}(s)\, \begin{bmatrix} w\\ u\end{bmatrix} = \begin{bmatrix}P_{11}(s) & P_{12}(s)\\P_{21}(s) & P_{22}(s)\end{bmatrix} \, \begin{bmatrix} w\\ u\end{bmatrix}
u = \mathbf{K}(s) \, v

It is therefore possible to express the dependency of z on w as:

z=F_\ell(\mathbf{P},\mathbf{K})\,w

Called the lower linear fractional transformation, F_\ell is defined (the subscript comes from lower):

F_\ell(\mathbf{P},\mathbf{K}) = P_{11} + P_{12}\,\mathbf{K}\,(I-P_{22}\,\mathbf{K})^{-1}\,P_{21}

Therefore, the objective of \mathcal{H}_\infty control design is to find a controller \mathbf{K} such that F_\ell(\mathbf{P},\mathbf{K}) is minimised according to the \mathcal{H}_\infty norm. The same definition applies to \mathcal{H}_2 control design. The infinity norm of the transfer function matrix F_\ell(\mathbf{P},\mathbf{K}) is defined as:

||F_\ell(\mathbf{P},\mathbf{K})||_\infty = \sup_\omega \bar{\sigma}(F_\ell(\mathbf{P},\mathbf{K})(j\omega))

where \bar{\sigma} is the maximum singular value of the matrix F_\ell(\mathbf{P},\mathbf{K})(j\omega).

The achievable H norm of the closed loop system is mainly given through the matrix D11 (when the system P is given in the form (A, B1, B2, C1, C2, D11, D12, D22, D21)). There are several ways to come to an H controller:

See also

References

  1. George Zames (1981). "Feedback and optimal sensitivity: Model reference transformations, multiplicative seminorms, and approximate inverses". IEEE Trans. Automatic Control 26 (2): 301–320. doi:10.1109/tac.1981.1102603.
  2. J. William Helton (1978). "Orbit structure of the Mobius transformation semigroup action on H-infinity (broadband matching)". Adv. in Math. Suppl. Stud. 3: 129–197.
  3. Allen Tannenbaum (1980). "Feedback stabilization of linear dynamical plants with uncertainty in the gain factor". International Journal of Control 32 (1): 1–16. doi:10.1080/00207178008922838.

Bibliography