Register to reply

Controllability Matrix [Control Theory]

Share this thread:
yaang
#1
Nov14-10, 09:26 AM
P: 22
I looked at several books as well as internet sources, none of them explain how the controllability matrix is formed.

Given the linear time invariant system

x'(t)=Ax(t)+Bu(t)

A is an nxn matrix,
B is an nx1 matrix, (assuming single input)

Then controllability matrix R is given by:

R= [B AB A^2B A^3B .....A^(n-1)B]

System is controllable if Det(R)=/=0 or rank(R)=n

Can someone explain me the logic behind how this matrix was formed ?
Phys.Org News Partner Science news on Phys.org
Physical constant is constant even in strong gravitational fields
Montreal VR headset team turns to crowdfunding for Totem
Researchers study vital 'on/off switches' that control when bacteria turn deadly
paulfr
#2
Nov15-10, 01:27 PM
P: 179
I would begin by taking a very simple system and trying to see how
the C Matrix guarantees that the closed loop poles will be
in the Left Half Plane.

An interesting question is whether the Matrix can insure not only
Global Stability, but also Relative Stability ..... the poles will be
near the Real Axis for a stable, non oscillatory, transient response.
chingkui
#3
Nov15-10, 03:05 PM
P: 197
http://teal.gmu.edu/ececourses/ece52...ap2/node1.html

trambolin
#4
Nov15-10, 04:36 PM
P: 341
Controllability Matrix [Control Theory]

A slightly clearer but somewhat less rigorous connection (only C'bility [itex]\Longrightarrow[/itex] rank condition!) can be made as follows: We can solve the diff. eq. system that you have provided and obtain
[tex]
x(t) = \int_0^{\infty}e^{A(t-\tau)}Bu(\tau)d\tau + e^{At}x(0)
[/tex]

Let's assume zero initial conditions for simplicity. Now, since the controllability means that I can reach any x(t), the integral converges to x(t) with some u(t). Let's use the Taylor series of exponential

[tex]
x(t) = \int_0^{\infty}\left(I+A(t-\tau) + \frac{A^2}{2!}(t-\tau)^2+\cdots \right)Bu(\tau)d\tau
[/tex]

You can take the constant terms out and obtain a matrix-vector multiplication (though infinite dimensional)

[tex]
x(t) = \begin{bmatrix}B &AB &A^2B &\cdots\end{bmatrix}\begin{pmatrix}\int_0^{\infty}u(\tau)d\tau \\\int_0^{\infty}(t-\tau)u(\tau)d\tau \\ \int_0^{\infty}\frac{1}{2!}(t-\tau)^2u(\tau)d\tau\\ \vdots\end{pmatrix} = \mathcal{C}_\infty \mathcal{U}
[/tex]

I would denote the matrix part as [itex]\mathcal{C}_\infty[/itex] . Now, since we assume controllability, we should be able to obtain any x(t), hence [itex]\mathcal{C}_\infty[/itex] must be full row rank. But from Cayley-Hamilton theorem we know that the powers of A with degree higher then n-1, can be rewritten by the powers of A up to the degree n-1. (This is a bad sentence but looking it up is easy so I skip that part.) This means that no extra information about the rank of this matrix can be included after the [itex]A^{n-1}B[/itex] since the remaining terms are linear combinations of the first n terms. Thus,

[tex]rank(\mathcal{C}_\infty) = rank(\mathcal{C}) = rank(\begin{bmatrix}B &AB &A^2B &\cdots &A^{n-1}B\end{bmatrix}[/tex]

In case of SISO systems, [itex]\mathcal{C}[/itex] happens to be square so the rank condition equals to the determinant being nonzero.
yaang
#5
Nov15-10, 04:56 PM
P: 22
Thanks for the help everybody , after studying the example 1.2.1 @ http://teal.gmu.edu/ececourses/ece52...00000000000000

and trambolin's solution it's crystal clear now.

(How can i add "Solved" to the thread title ?)


Register to reply

Related Discussions
Controllability of state space equation Engineering, Comp Sci, & Technology Homework 1
Control Theory..please advise if cont theory should be in system design threads Electrical Engineering 6
Control Theory Engineering, Comp Sci, & Technology Homework 5
Control theory Engineering, Comp Sci, & Technology Homework 1
Control theory Calculus 1