Considering a trivial LTI:

x' = Ax + Bu

y = Cx + Du

I was told that from this representation I can always derive a model without the direct feedthourgh matrix D in order to solve Lyapunov equation.

I was trying to do this by including the input vector u inside the state vector x, but i can't understand how.

What is the classical approach to this problem?

Similar questions and discussions