The basic state-space model in innovations form can be written
x(t+1) = A x(t) + B u(t) + K e(t)
y(t) = C x(t) + D u(t) + e(t)
The SITB supports two kinds of parametrizations of state-space models: Black box, free parametrizations and parametrizations tailor made to the application.
To estimate a state space model, select Parametric Models in the pop-up menu Estimate in the ident window. In the dialog that opens, choose State space as the model structure, and enter the desired order. The Order editor handles options to fix K, D and x0 to zero in the black box case.
Note that in the black box case, there is a special feature to select the model order by entering a vector (like 1:10) for the model order. You can push the Order selection button to fill out the model order field in this case.
By entering a vector (e.g. 1:10), all orders will be computed using a preliminary method, and you will have to choose order(s) based on information in a special graph. You can also use the button Order selection to fill out the Model order field with a default model order range.
There are two basic methods for the estimation: PEM and N4SID
There are a number of structure options that can be reached in the Order Editor:
By fixing the matrix K to zero, an Output Error method is obtained.
The initial state vector x(0) can either be estimated or fixed to zero.
The D matrix can either be estimated or fixed to zero.
The search for a minimum is controlled by a number of options. These are accessed from the Options button in the Iteration control... dialog.
The quality of the resulting estimates may significantly depend on an auxiliary order (like a prediction horizon). This auxiliary order can be given within parentheses after the order in the Orders: edit box (e.g., 4 (9)). If no auxiliary order is given, a default value is used.
If the auxiliary order is given as a vector (e.g. 4 (5:17)), models for all these auxiliary orders are computed and evaluated. The model that gives the best fit between model and measured output is chosen. A figure is opened that illustrates the fit as a function of the auxiliary order. The fit is either a simulation error fit if K=0 is chosen, or a one-step ahead prediction error fit.
See the commands PEM and N4SID in the manual for more information.
To use them in conjunction with ident, define the appropriate structure in the MATLAB command line and enter its variable name in the Orders: edit box of the Parametric Models dialog. If desired, select the appropriate iteration options for PEM by pressing the Iteration control... button.
(file iduiss.htm)