| Robust Control Toolbox™ | |
| Provide feedback about this page |
Extract muinfo structure returned by mussv
Syntax
Description
A structured singular value computation of the form
returns detailed information in the structure muinfo. mussvextract is used to extract the compressed information within muinfo into a readable form.
The most general call to mussvextract extracts three usable quantities: VDelta, VSigma, and VLmi. VDelta is used to verify the lower bound. VSigma is used to verify the Newlin/Young upper bound and has fields DLeft, DRight, GLeft, GMiddle, and GRight. VLmi is used to verify the LMI upper bound and has fields Dr, Dc, Grc, and Gcr. The relation/interpretation of these quantities with the numerical results in bounds is described below.
Upper Bound Information
The upper bound is based on a proof that det(I - M*Delta) is nonzero for all block-structured matrices Delta with norm smaller than 1/bounds(1). The Newlin/Young method consists of finding a scalar 
and matrices D and G, consistent with BlockStructure, such that
Because some uncertainty blocks and M need not be square, the matrices D and G have a few different manifestations. In fact, in the formula above, there are a left and right D and G, as well as a middle G. Any such is an upper bound of mussv(M,BlockStructure).
It is true that if BlockStructure consists only of complex blocks, then all G matrices will be zero, and the expression above simplifies to
The LMI method consists of finding a scalar and matrices D and G, consistent with BlockStructure, such that
is negative semidefinite. Again, D and G have a few different manifestations to match the row and column dimensions of M. Any such is an upper bound of mussv(M,BlockStructure). If BlockStructure consists only of complex blocks, then all G matrices will be zero, and negative semidefiniteness of M´Dr M-2Dc is sufficient to derive an upper bound.
Lower Bound Information
The lower bound of mussv(M,BlockStructure) is based on finding a "small" (hopefully the smallest) block-structured matrix VDelta that causes det(I - M*VDelta) to equal 0. Equivalently, the matrix M*VDelta has an eigenvalue equal to 1. It will always be true that the lower bound (bounds(2)) will be the reciprocal of norm(VDelta).
Example
Suppose M is a 4-by-4 complex matrix. Take the block structure to be two 1-by-1 complex blocks and one 2-by-2 complex block.
You can calculate bounds on the structured singular value using the mussv command and extract the scaling matrices using mussvextract.
You can first verify the Newlin/Young upper bound with the information extracted from muinfo. The corresponding scalings are Dl and Dr.
Dl = VSigma.DLeft Dl = 1.0000e+000 0 0 0 0 9.9190e-001 0 0 0 0 1.1255e+000 0 0 0 0 1.1255e+000 Dr = VSigma.DRight Dr = 1.0000e+000 0 0 0 0 9.9190e-001 0 0 0 0 1.1255e+000 0 0 0 0 1.1255e+000 [norm(Dl*M/Dr) bounds(1)] ans = 4.3420e+000 4.3420e+000
You can first verify the LMI upper bound with the information extracted from muinfo. The corresponding scalings are Dr and Dc.
Dr = VLmi.Dr; Dc = VLmi.Dc; eig(M'*Dr*M - bounds(1)^2*Dc) ans = -2.0045e-005 +6.1649e-016i -1.4688e+001 -2.4975e-016i -2.0436e+001 -4.7583e-016i -1.9100e+001 +1.4136e-015i
Note that VDelta matches the structure defined by BlockStructure, and the norm of VDelta agrees with the lower bound,
VDelta VDelta = 1.0698e-001 -2.0405e-001i 0 0 0 0 1.4920e-001 +1.7556e-001i 0 0 0 0 -5.4173e-002 -1.0932e-002i -5.7140e-002 +9.6497e-002i 0 0 2.8071e-002 -8.0807e-002i -1.3608e-001 -1.0777e-001i [norm(VDelta) 1/bounds(2)] ans = 0.2304 0.2304
and that M*VDelta has an eigenvalue exactly at 1.
eig(M*VDelta) ans = 1.0000e+000 -8.3267e-017i -6.1108e-002 +2.5748e-001i 4.1427e-018 -5.8578e-018i -1.9637e-001 -5.6540e-002i
Keep the matrix the same, but change BlockStructure to be a 2-by-2 repeated, real scalar block and two complex 1-by-1 blocks. Run mussv with the 'C' option to tighten the upper bound.
You can compare the computed bounds. Note that bounds2 should be smaller than bounds, because the uncertainty set defined by BlockStructure2 is a proper subset of that defined by BlockStructure.
You can extract the D, G and Delta from muinfo2 using mussvextract.
As before, you can first verify the Newlin/Young upper bound with the information extracted from muinfo. The corresponding scalings are Dl, Dr, Gl, Gm and Gr.
Dl = VSigma2.DLeft; Dr = VSigma2.DRight; Gl = VSigma2.GLeft; Gm = VSigma2.GMiddle; Gr = VSigma2.GRight; dmd = Dl*M/Dr/bounds2(1) - sqrt(-1)*Gm; SL = (eye(4)+Gl*Gl)^-0.25; SR = (eye(4)+Gr*Gr)^-0.25; norm(SL*dmd*SR) ans = 1.0000
You can first verify the LMI upper bound with the information extracted from muinfo. The corresponding scalings are Dr, Dc, Grc and Gcr.
Dr = VLmi2.Dr; Dc = VLmi2.Dc; Grc = VLmi2.Grc; Gcr = VLmi2.Gcr; eig(M'*Dr*M - bounds(1)^2 *Dc + j*(Gcr*M-M'*Grc)) ans = -4.4665e-002 -4.4823e-019i -5.2486e-004 +1.5623e-018i -1.8028e-003 +3.2493e-019i -1.2558e-003 +1.2973e-019i
VDelta2 matches the structure defined by BlockStructure, and the norm of VDelta2 agrees with the lower bound,
VDelta2 VDelta2 = 0.2882 0 0 0 0 0.2882 0 0 0 0 -0.152 - 0.2448i 0 0 0 0 -0.0395 -0.2855i [norm(VDelta2) 1/bounds2(2)] ans = 0.2882 0.2882
and that M*VDelta2 has an eigenvalue exactly at 1.
eig(M*VDelta2) ans = -3.3623e-001 +2.1885e-001i -3.6805e-001 -1.5645e-001i 1.0000e+000 -1.4169e-016i 4.5066e-001 -3.4481e-001i
See Also
mussv (µ)
| Provide feedback about this page |
| mussv | ncfmargin | |
| © 1984-2008- The MathWorks, Inc. - Site Help - Patents - Trademarks - Privacy Policy - Preventing Piracy - RSS |
