Cálculo de la matriz de transición de estado a partir del modelo de espacio de estado de CCF

Dadas las matrices de coeficientes de estado para la ecuación de estado

$$ \ dot {x} = Axe + Bu $$

$$ A = \ begin {bmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -1 & -3 & -3 \\ \ end {bmatrix} B = \ begin {bmatrix} 0 \\ 0 \\ 1 \\ \ end {bmatrix} $$

Para calcular la matriz de transición de estado, ¿existe un método más fácil que la técnica de Laplace inversa cuando las matrices de coeficientes se dan en forma canónica controlable?