ctrbf
Compute controllability staircase form
Syntax
(Abar,Bbar,Cbar,T,k] = ctrbf(A,B,C)
ctrbf(A,B,C,tol)
Description
If the controllability matrix of (A,B) has rankr≤n, wherenis the size ofA, then there exists a similarity transformation such that
whereTis unitary, and the transformed system has astaircaseform, in which the uncontrollable modes, if there are any, are in the upper left corner.
where (Ac,Bc) is controllable, all eigenvalues ofAucare uncontrollable, and .
(Abar,Bbar,Cbar,T,k] = ctrbf(A,B,C)
decomposes the state-space system represented byA
,B
, andC
into the controllability staircase form,Abar
,Bbar
, andCbar
, described above.T
is the similarity transformation matrix andk
is a vector of lengthn, wherenis the order of the system represented byA
. Each entry ofk
represents the number of controllable states factored out during each step of the transformation matrix calculation. The number of nonzero elements ink
indicates how many iterations were necessary to calculateT
, andsum(k)
is the number of states inAc, the controllable portion ofAbar
.
ctrbf(A,B,C,tol)
uses the tolerancetol
when calculating the controllable/uncontrollable subspaces. When the tolerance is not specified, it defaults to10*n*norm(A,1)*eps
.
Examples
Compute the controllability staircase form for
A = 1 1 4 -2 B = 1 -1 1 -1 C = 1 0 0 1
and locate the uncontrollable mode.
(Abar,Bbar,Cbar,T,k]=ctrbf(A,B,C) Abar = -3.0000 0 -3.0000 2.0000 Bbar = 0.0000 0.0000 1.4142 -1.4142 Cbar = -0.7071 0.7071 0.7071 0.7071 T = -0.7071 0.7071 0.7071 0.7071 k = 1 0
The decomposed systemAbar
shows an uncontrollable mode located at -3 and a controllable mode located at 2.
Algorithms
ctrbf
implements the Staircase Algorithm of(1].
References
(1] Rosenbrock, M.M.,State-Space and Multivariable Theory, John Wiley, 1970.