Estimate parameters of Output-Error model from SISO data in Simulink software returning idpoly object
The OE block estimates the parameters of an Output-Error model, and returns the estimated model as an idpoly object.
Note: The OE Estimator block will be removed in a future release. There is no replacement for this block.
For information about the default algorithm settings used for model estimation, see oeOptions.
Each estimation generates a figure with the following plots:
Actual (measured) output versus the simulated or predicted model output.
Error in simulated model, which is the difference between the measured output and the model output.
The output-error model is defined, as follows:
w is the undisturbed output.
y(t) is the output at time t.
and are the parameters to be estimated.
is the number of poles of the transfer function from the input to the undisturbed output.
is the number of zeros of the transfer function from the input to the undisturbed output.
is the number of input samples that occur before the inputs that affect the current output.
are the previous inputs on which the current output depends.
e(t) is a white-noise disturbance value.
The OE model can also be written in a compact way using the following notation:
and is the backward shift operator, defined by .
The following block diagram shows the ARX model structure.
The block accepts two inputs, corresponding to the measured input-output data for estimating the model.
First input: Input signal.
Second input: Output signal.
The OE Estimator block outputs a sequence of multiple models (idpoly), estimated at regular intervals during the simulation.
The Length of Data window field in the block parameter dialog box specifies the number of data samples to use for estimation, as the simulation progresses.
The output format depends on whether you specify the Model Name in the block parameter dialog box.
Integers nb, nf, and nk specify the number of B and F model parameters and nk is the input-output delay, respectively.
Number of input data samples that specify the interval after which to estimate a new model.
Sampling time for the model.
Number of past data samples used to estimate each model. A longer data window should be used for higher-order models. Too small a value might cause poor estimation results, and too large a value leads to slower computation.
Name of the model.
Whether you specify the model name determines the output format of the resulting models, as follows:
If you do not specify a model name, the estimated models display in the MATLAB® Command Window in a transfer-function format.
If you specify a model name, the resulting models are output to the MATLAB workspace as a cell array.
Simulation: The algorithm uses only measured input data to simulate the response of the model.
Prediction: Specifies the forward-prediction horizon for computing the response K steps in the future, where K is 1, 5, or 10.
This example shows how you can use the OE Estimator block in a Simulink® model.
Specify the data from iddata1.mat for estimation:
load iddata1; IODATA = z1;
Create a new Simulink model, as follows:
Add the IDDATA Source block and specify IODATA in the Iddata object field of the IDDATA Source block parameters dialog box.
Add the OE Estimator block to the model. Set sample time in the block to 0.1 seconds and the simulation end time to 30 seconds.
Connect the Input and Output ports of the IDDATA Source block to the u and y ports of the OE Estimator block, respectively.
Run the simulation.
The estimated models appear in the MATLAB Command Window every 25 samples.