I am currently working on a project where I need to decompose my system into observable and unobservable subsystems in an efficient way, so I was looking for a function that could help me with that.
PS: I know about this function and it is not was I am looking for:
s = co.ss(A, B, C, D)
sys, T = co.observable_form(s)
In this case, the system needs to be fully observable.
Thank you!
In case someone is looking for an answer to this, I found a great library that has a function that does just this. It is called "harold" and you can find it here. The said function is called "staircase" and it is equivalent to the obsvf or ctrbf functions from MATLAB; it implements the Staircase Algorithm.