How to hack tfest/oe identification algorithm in order to set constraints during estimation
2 views (last 30 days)
Good morning community,
I would like to hack into the estimation algorithm for system identification using the commands tfest(...) or oe(...). From the documentation it seems they use the same algorithm under the hood, specifically it is said:
"For discrete-time data, tfest uses the same algorithm as oe to determine the numerator and denominator polynomial coefficients. In this algorithm, the initialization is performed using arx, followed by nonlinear least-squares search-based updates to minimize a weighted prediction error norm."
I'm actually using discrete time-domain data for a MISO systems, specifically with two inputs. The structure of the model is the following
I would like to impose the constraints
How can I do this?
Khaled Aljanaideh on 24 Feb 2021
If the coloring of the noise term (e) in the model structure is not important, then this can be done by using the ARMAX or ARX model structures. Note that the arx model structure, for example, can be written as , which is equivalent to , where the model structure shown in the question is