MIT OpenCourseWare
http://ocw.mit.edu
16.323 Principles of Optimal Control
Spring 2008
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
Topic #14
16.31 Feedback Control Systems
MIMO Systems
• Singular Value Decomposition
Spr 2008 16.323 14–1
Multivariable Frequency Response
• In the MIMO case, the system G(s) is described by a p × m transfer
function matrix (TFM)
– Still have that G(s) = C(sI − A)−1B + D
– But G(s) → A, B, C, D MUCH less obvious than in SISO case.
– Also seen that the discussion of poles and zeros of MIMO systems
is much more complicated.
• In SISO case we use the Bode plot to develop a measure of the system
“size”.
– Given z = Gw, where G(jω) = |G(jω)|ejφ(w)
– Then w = |w|ej(ω1t+ψ) applied to |G(jω)|ejφ(w) yields
|w||G(jω1)|e j(ω1t+ψ+φ(ω1)) = |z |ej(ω1t+ψo) ≡ z
– Amplification and phase shift of the input signal obvious in the SISO
case.
• MIMO extension?
– Is the response of the system
� large or small?�
103/s 0
G(s) =
0 10−3/s
June 18, 2008
Spr 2008 16.323 14–2
• For MIMO systems, cannot just plot all of the gij elements of G
– Ignores the coupling that might exist between them.
– So not enlightening.
• Basic MIMO frequency response:
– Restrict all inputs to be at the same frequency
– Determine how the system responds at that frequency
– See how this response changes with frequency
• So inputs are w = wcejωt, where wc ∈ Cm
– Then we get z = G(s)|s=jω w, ⇒ z = zcejωt and zc ∈ Cp
– We need only analyze zc = G(jω)wc
• As in the SISO case, we need a way to establish if the system response
is large or small.
– How much amplification we can get with a bounded input.
�
• Consider zc = G(jω)wc and set �wc�2 = wcH wc ≤ 1. What can
we say about the �zc�2?
– Answer depends on ω and on the direction of the input wc
– Best found using singular values.
June 18, 2008
Spr 2008 16.323 14–3
Singular Value Decomposition
• Must perform SVD of the matrix G(s) at each frequency s = jω
G(jω) ∈ Cp×m U ∈ Cp×p Σ ∈ Rp×m V ∈ Cm×m
G = U ΣV H
– U H U = I, U U H = I, V H V = I, V V H = I, and Σ is diagonal.
– Diagonal elements σk ≥ 0 of Σ are the singular values of G.
� �
σi = λi(GH G) or σi = λi(GGH )
the positive ones are the same from both formulas.
– Columns of matrices U and V (ui and vj ) are the associated eigen
vectors
GH Gvj = σj2vj
GGH ui = σi2ui
Gvi = σiui
• If the rank(G) = r ≤ min(p, m), then
– σk > 0, k = 1, . . . , r
– σk = 0, k = r + 1, . . . , min(p, m)
– Singular values are sorted so that σ1 ≥ σ2 ≥ . . . ≥ σr
• An SVD gives a very detailed description of how a matrix (the
system G) acts on a vector (the input w) at a particular frequency.
June 18, 2008
Spr 2008 16.323 14–4
• So how can we use this result?
– Fix the size �wc�2 = 1 of the input, and see how large we can make
the output.
– Since we are working at a single frequency, we just analyze the
relation
zc = Gw wc, Gw ≡ G(s = jω)
• Define the maximum and minimum amplifications as:
σ ≡ max �zc�2
�wc �2 =1
σ ≡ min �zc�2
�wc �2 =1
• Then we have that (let q = min(p, m))
σ = �
σ1
σq p ≥ m “tall”
σ =
0 p < m “wide”
• Can use σ and σ to determine the possible amplification and attenu
ation of the input signals.
• Since G(s) changes with frequency, so will σ and σ
June 18, 2008
Spr 2008 16.323 14–5
SVD Example
• Consider (wide case)
⎡ ⎤
� � � �� � 1 0 0
5 0 0 1 0 5 0 0 ⎣
Gw = = 0 1 0⎦
0 0.5 0 0 1 0 0.5 0
0 0 1
= U ΣV H
so that σ1 = 5 and σ2 = 0.5
σ ≡ max �Gw wc�2 = 5 = σ1
�wc �2 =1
σ ≡ min �Gw wc�2 = 0 �= σ2
�wc �2 =1
• But now consider (tall case)
⎡ ⎤ ⎡ ⎤⎡ ⎤
5 0 1 0 0 5 0 � �
1 0
G̃w = ⎣ 0 0.5 ⎦ = ⎣ 0 1 0 ⎦ ⎣ 0 0.5 ⎦
0 1
0 0 0 0 1 0 0
= U ΣV H
so that σ1 = 5 and σ2 = 0.5 still.
σ ≡ max �Gw wc�2 = 5 = σ1
�wc �2 =1
σ ≡ min �Gw wc�2 = 0.5 = σ2
�wc �2 =1
June 18, 2008
Spr 2008 16.323 14–6
• For MIMO systems, the gains (or σ’s) are only part of the story, as
we must also consider the input direction.
• To analyze this point further, note that we can rewrite
⎡ ⎤
σ1
⎡
H ⎤
⎢ . . .
⎥ v1
⎥ ⎣
.. ⎦
Gw
=
U ΣV H
= u1 . . . up
⎢
� �
.
⎢
σm
⎦
H
⎥
vm
⎣
� m
= σiuiviH
i=1
– Assumed tall case for simplicity, so p > m and q = m
• Can now analyze impact of various alternatives for the input
– Only looking at one frequency, so the basic signal is harmonic.
– But, we are free to pick the relative sizes and phases of each of
the components of the input vector wc.
� These define the input direction
June 18, 2008
Spr 2008 16.323 14–7
• For example, we could pick wc = v1, then
� m �
�
zc = Gw wc = σiuiviH v1 = σ1u1
i=1
since viH vj = δij .
– Output amplified by σ1. The relative sizes and phases of each of
the components of the output are given by the vector zc.
• By selecting other input directions (at the same frequency), we can
get quite different amplifications of the input signal
�Gw wc�2
σ≤ ≤σ
�wc�2
• Thus we say that
– Gw is large if σ(Gw ) � 1
– Gw is small if σ(Gw ) � 1
• MIMO frequency response are plots of σ(jω) and σ(jω).
– Then use the singular value vectors to analyze the response at a
particular frequency.
June 18, 2008
