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?