oe: Estimate Output-Error Models

In sysid: System Identification in R

Description

Fit an output-error model of the specified order given the input-output data

Usage

oe(x, order = c(1, 1, 0), init_sys = NULL, options = optimOptions())

Arguments

x	an object of class idframe
order	Specification of the orders: the four integer components (nb,nf,nk) are order of polynomial B + 1, order of the polynomial F, and the input-output delay respectively
init_sys	Linear polynomial model that configures the initial parameterization. Must be an OE model. Overrules the order argument
options	Estimation Options, setup using optimOptions

Details

SISO OE models are of the form

y[k] + f_1 y[k-1] + … + f_{nf} y[k-nf] = b_{nk} u[k-nk] + … + b_{nk+nb} u[k-nk-nb] + f_{1} e[k-1] + … f_{nf} e[k-nf] + e[k]

The function estimates the coefficients using non-linear least squares (Levenberg-Marquardt Algorithm)
The data is expected to have no offsets or trends. They can be removed using the detrend function.

Value

An object of class estpoly containing the following elements:

sys	an idpoly object containing the fitted OE coefficients
fitted.values	the predicted response
residuals	the residuals
input	the input data used
call	the matched call
stats	A list containing the following fields: vcov - the covariance matrix of the fitted coefficients sigma - the standard deviation of the innovations
options	Option set used for estimation. If no custom options were configured, this is a set of default options
termination	Termination conditions for the iterative search used for prediction error minimization: WhyStop - Reason for termination iter - Number of Iterations iter - Number of Function Evaluations

References

Arun K. Tangirala (2015), Principles of System Identification: Theory and Practice, CRC Press, Boca Raton. Sections 14.4.1, 17.5.2, 21.6.3

Examples

data(oesim)
z <- dataSlice(oesim,end=1533) # training set
mod_oe <- oe(z,c(2,1,2))
mod_oe 
plot(mod_oe) # plot the predicted and actual responses

sysid documentation built on May 2, 2019, 4:18 a.m.