Set or randomize initial parameter values
m = init(m0)
m = init(m0,R,pars,sp)
m = init(m0) randomizes initial parameter estimates for model
m0 for any linear or nonlinear identified model. It does
m is the same model
m0, but with a different nominal parameter vector. This
vector is used as the initial estimate by
m = init(m0,R,pars,sp) randomizes parameters around
pars with variances given by the row vector
Parameter number k is randomized as
e is a normal random
variable with zero mean and a variance of 1. The default value of
is all ones, and the default value of
pars is the nominal parameter
Only models that give stable predictors are accepted. If
sp = 'b',
only models that are both stable and have stable predictors are accepted.
sp = 's' requires stability only of the model, and
'p' requires stability only of the predictor.
sp = 'p'
is the default.
Sufficiently free parameterizations can be stabilized by direct means without any
random search. To just stabilize such an initial model, set
R = 0.
R > 0, randomization is also done.
For model structures where a random search is necessary to find a stable
model/predictor, a maximum of 100 trials is made by
init. It can be
difficult to find a stable predictor for high-order systems by trial and error.