Consider the multivariable feedback control system shown in the following figure. In order to quantify the multivariable stability margins and performance of such systems, you can use the singular values of the closed-loop transfer function matrices from r to each of the three outputs e, u, and y, viz.
where the L(s) is the loop transfer function matrix
Block Diagram of the Multivariable Feedback Control System
The two matrices S(s) and T(s) are known as the sensitivity function and complementary sensitivity function, respectively. The matrix R(s) has no common name. The singular value Bode plots of each of the three transfer function matrices S(s), R(s), and T(s) play an important role in robust multivariable control system design. The singular values of the loop transfer function matrix L(s) are important because L(s) determines the matrices S(s) and T(s).
The singular values of S(jω) determine the disturbance attenuation, because S(s) is in fact the closed-loop transfer function from disturbance d to plant output y — see Block Diagram of the Multivariable Feedback Control System. Thus a disturbance attenuation performance specification can be written as
where is the desired disturbance attenuation factor. Allowing to depend on frequency ω enables you to specify a different attenuation factor for each frequency ω.
The singular value Bode plots of R(s) and of T(s) are used to measure the stability margins of multivariable feedback designs in the face of additive plant perturbations ΔA and multiplicative plant perturbations ΔM, respectively. See the following figure.
Consider how the singular value Bode plot of complementary sensitivity T(s) determines the stability margin for multiplicative perturbations ΔM. The multiplicative stability margin is, by definition, the "size" of the smallest stable ΔM(s) that destabilizes the system in the figure below when ΔA = 0.
Taking to be the definition of the "size" of ΔM(jω), you have the following useful characterization of "multiplicative" stability robustness:
The size of the smallest destabilizing multiplicative uncertainty ΔM(s) is:
The smaller is , the greater will be the size of the smallest destabilizing multiplicative perturbation, and hence the greater will be the stability margin.
A similar result is available for relating the stability margin in the face of additive plant perturbations ΔA(s) to R(s) if you take to be the definition of the "size" of ΔA(jω) at frequency ω.
The size of the smallest destabilizing additive uncertainty ΔA is:
As a consequence of robustness theorems 1 and 2, it is common to specify the stability margins of control systems via singular value inequalities such as
where |W2(jω)| and |W3(jω)| are the respective sizes of the largest anticipated additive and multiplicative plant perturbations.
It is common practice to lump the effects of all plant uncertainty into a single fictitious multiplicative perturbation ΔM, so that the control design requirements can be written
as shown in Singular Value Specifications on L, S, and T.
It is interesting to note that in the upper half of the figure (above the 0 dB line),
while in the lower half of Singular Value Specifications on L, S, and T (below the 0 dB line),
This results from the fact that
if , and
Singular Value Specifications on L, S, and T
Thus, it is not uncommon to see specifications on disturbance attenuation and multiplicative stability margin expressed directly in terms of forbidden regions for the Bode plots of σi(L(jω)) as "singular value loop shaping" requirements, either as specified upper/lower bounds or as a target desired loop shape — see the preceding figure.
For those who are more comfortable with classical single-loop concepts, there are the important connections between the multiplicative stability margins predicted by and those predicted by classical M-circles, as found on the Nichols chart. Indeed in the single-input/single-output case,
which is precisely the quantity you obtain from Nichols chart M-circles. Thus, is a multiloop generalization of the closed-loop resonant peak magnitude which, as classical control experts will recognize, is closely related to the damping ratio of the dominant closed-loop poles. Also, it turns out that you can relate , to the classical gain margin GM and phase margin θM in each feedback loop of the multivariable feedback system of Block Diagram of the Multivariable Feedback Control System via the formulas:
(See .) These formulas are valid provided and are larger than 1, as is normally the case. The margins apply even when the gain perturbations or phase perturbations occur simultaneously in several feedback channels.
The infinity norms of S and T also yield gain reduction tolerances. The gain reduction tolerance gm is defined to be the minimal amount by which the gains in each loop would have to be decreased in order to destabilize the system. Upper bounds on gm are as follows:
loopsyn lets you design a stabilizing feedback
controller to optimally shape the open loop frequency response of a MIMO feedback control system
to match as closely as possible a desired loop shape Gd. The basic syntax of the
loopsyn loop-shaping controller synthesis command is:
K = loopsyn(G,Gd)
G is the LTI transfer function matrix of a MIMO plant model,
Gd is the target desired loop shape for the loop transfer function
K is the optimal loop-shaping controller. The
K has the property that it shapes the loop
L=G*K so that it matches the frequency response of
closely as possible, subject to the constraint that the compensator must stabilize the plant