MOESP algorithm by using n4sid command for discrete time data
2 views (last 30 days)
I am trying to apply MOESP algorithm by using the n4sid command to disctete time data but the iddata object requires me to input the sampling time,which, is not possible for discrete time data. at the value of Ts = 1, I obtain the matrix A with abnormally high values, How do I correct this?
data = iddata(L',U',1)
sys = n4sid(data,2)
Command Window Output:
Discrete-time identified state-space model:
x(t+Ts) = A x(t) + B u(t) + K e(t)
y(t) = C x(t) + D u(t) + e(t)
x1 -2.517e+12 2.094e+12
x2 -3.026e+12 2.517e+12
Rajiv Singh on 9 Aug 2022
The issue is not with the sample time choice but with the system itself. It seems that the data from a system with undamped poles. In that case n4sid can run into numerical issues and deliver poor results.
Note that for discrete-time estimations using time-domain data, the sample time does not affect the numerical parameter values computed.