State-transition matrix explained

x

at an initial time

t0

gives

x

at a later time

t

. The state-transition matrix can be used to obtain the general solution of linear dynamical systems.

Linear systems solutions

The state-transition matrix is used to find the solution to a general state-space representation of a linear system in the following form

x

(t)=A(t)x(t)+B(t)u(t),x(t0)=x0

,where

x(t)

are the states of the system,

u(t)

is the input signal,

A(t)

and

B(t)

are matrix functions, and

x0

is the initial condition at

t0

. Using the state-transition matrix

\Phi(t,\tau)

, the solution is given by:[1] [2]

x(t)=\Phi(t,t0)x(t0)+\int

t
t0

\Phi(t,\tau)B(\tau)u(\tau)d\tau

The first term is known as the zero-input response and represents how the system's state would evolve in the absence of any input. The second term is known as the zero-state response and defines how the inputs impact the system.

Peano–Baker series

The most general transition matrix is given by a product integral, referred to as the Peano–Baker series

\begin{align} \Phi(t,\tau)=I&+

tA(\sigma
\int
1)d\sigma

1\\ &+

tA(\sigma
\int
1)\int
\sigma1
\tau

A(\sigma2)d\sigma2d\sigma1\\ &+

tA(\sigma
\int
1)\int
\sigma1
\tau

A(\sigma2)\int

\sigma2
\tau

A(\sigma3)d\sigma3d\sigma2d\sigma1\\ &+ \end{align}

where

I

is the identity matrix. This matrix converges uniformly and absolutely to a solution that exists and is unique. The series has a formal sum that can be written as

\Phi(t,\tau)=\exp

tA(\sigma)d\sigma
l{T}\int
\tau
where

l{T}

is the time-ordering operator, used to ensure that the repeated product integral is in proper order. The Magnus expansion provides a means for evaluating this product.

Other properties

The state transition matrix

\Phi

satisfies the following relationships. These relationships are generic to the product integral.

1. It is continuous and has continuous derivatives.

2, It is never singular; in fact

\Phi-1(t,\tau)=\Phi(\tau,t)

and

\Phi-1(t,\tau)\Phi(t,\tau)=I

, where

I

is the identity matrix.

3.

\Phi(t,t)=I

for all

t

.[3]

4.

\Phi(t2,t1)\Phi(t1,t0)=\Phi(t2,t0)

for all

t0\leqt1\leqt2

.

5. It satisfies the differential equation

\partial\Phi(t,t0)
\partialt

=A(t)\Phi(t,t0)

with initial conditions

\Phi(t0,t0)=I

.

6. The state-transition matrix

\Phi(t,\tau)

, given by

\Phi(t,\tau)\equivU(t)U-1(\tau)

where the

n x n

matrix

U(t)

is the fundamental solution matrix that satisfies
U

(t)=A(t)U(t)

with initial condition

U(t0)=I

.

7. Given the state

x(\tau)

at any time

\tau

, the state at any other time

t

is given by the mapping

x(t)=\Phi(t,\tau)x(\tau)

Estimation of the state-transition matrix

In the time-invariant case, we can define

\Phi

, using the matrix exponential, as

\Phi(t,t0)=

A(t-t0)
e
. [4]

In the time-variant case, the state-transition matrix

\Phi(t,t0)

can be estimated from the solutions of the differential equation
u

(t)=A(t)u(t)

with initial conditions

u(t0)

given by

[1, 0,\ldots, 0]T

,

[0, 1,\ldots, 0]T

, ...,

[0, 0,\ldots, 1]T

. The corresponding solutions provide the

n

columns of matrix

\Phi(t,t0)

. Now, from property 4,

\Phi(t,\tau)=\Phi(t,t0)\Phi(\tau,

-1
t
0)
for all

t0\leq\tau\leqt

. The state-transition matrix must be determined before analysis on the time-varying solution can continue.

See also

Further reading

Notes and References

  1. Baake. Michael. Schlaegel. Ulrike. The Peano Baker Series. Proceedings of the Steklov Institute of Mathematics. 2011. 275. 155–159. 10.1134/S0081543811080098. 119133539.
  2. Book: Rugh. Wilson. Linear System Theory. 1996. Prentice Hall. Upper Saddle River, NJ . 0-13-441205-2.
  3. Book: Brockett, Roger W.. Finite Dimensional Linear Systems. John Wiley & Sons. 1970. 978-0-471-10585-5.
  4. Reyneke . Pieter V. . Polynomial Filtering: To any degree on irregularly sampled data . Automatika . 2012 . 53 . 4 . 382–397. 10.7305/automatika.53-4.248 . 40282943 . free . 2263/21017 . free .