control_gramian: Controllability Gramian
In netcontrol: Control Theory Methods for Networks

Description Usage Arguments Value Examples

View source: R/control_stats.R

Compute the (infinite time) controllability Gramian for the discrete linear time invariant system described by x(t+1) = Ax(t) + Bu(t). The infinite time controllability Gramian is the solution to the discrete Lyapunov equation AWA^\prime-W = -BB^\prime, while the finite time Gramian for time T is

W_t = ∑_{t = 0}^T A^tBB^\prime(A^\prime)^t

1	control_gramian(A, B, t = NA)

`A`	n x n matrix.
`B`	n x m matrix.
`t`	Either NA for infinite time Gramian, or a positive non-zero integer. Defaults to NA.

The infinite time or finite time controllability Gramian

A = matrix(c(0,-3,-2,2,-2,1,-1,2,-1), 3,3)
B = diag(3)

#Infinite time Gramian
W_inf = control_gramian(A, B)

#4 time Gramian
W_4 = control_gramian(A,B,4)