idnlarx/findop

[X,U] = findop(sys,'steady',InputLevel,OutputLevel) returns the operating-point state values, X, and input values, U, for theidnlarx model, sys, using steady-state input and output specifications.

[X,U] = findop(sys,spec) returns the steady-state operating point for sys using the operating-point specification, spec.

[X,U] = findop(___,Options) specifies optimization search options for all of the previous syntaxes.

[X,U,Report] = findop(___) returns a summary report on the optimization search results for all of the previous syntaxes.

[X,U] = findop(sys,'snapshot',T,Uin) returns the operating point for sys at a simulation snapshot at time, T, using the specified input, Uin. The initial states of sys are assumed to be zero.

[X,U] = findop(sys,'snapshot',T,Uin,X0) specifies the initial states of the simulation.

Examples

Find Steady-State Nonlinear ARX Operating Point Using Default Specifications

Estimate a nonlinear ARX model.

load iddata6;
M = nlarx(z6,[4 3 1]);

Find the steady-state operating point where the input level is fixed to 1 and the output is unknown.

[X,U] = findop(M,'steady',1,NaN);

Find Nonlinear ARX Operating Point Using Additional Specifications

Estimate a nonlinear ARX model.

load iddata7;
M = nlarx(z7,[4 3*ones(1,2) 2*ones(1,2)]);

Create a default operating point specification object.

spec = operspec(M);

Set the values for the input signals.

spec.Input.Value(1) = -1;
spec.Input.Value(2) = 1;

Set the maximum and minimum values for the output signal.

spec.Output.Max = 10;
spec.Output.Min = -10;

Find the steady-state operating point using the given specifications.

[X,U] = findop(M,spec);

Find Nonlinear ARX Operating Point Using Custom Options

Estimate a nonlinear ARX model.

load iddata6;
M = nlarx(z6,[4 3 2]);

Create a default findopOptions option set.

opt = findopOptions(M);

Modify the option set to specify a steepest descent gradient search method with a maximum of 50 iterations.

opt.SearchMethod = 'grad';
opt.SearchOptions.MaxIterations = 50;

Find the steady-state operating point using the specified options.

[X,U] = findop(M,'steady',1,1,opt);

Retrieve Nonlinear ARX Operating Point Search Report

Estimate a nonlinear ARX model.

load iddata7;
M = nlarx(z7,[4 3*ones(1,2) 2*ones(1,2)]);

Find the steady-state operating point where input 1 is set to 1 and input 2 is unrestricted. The initial guess for the output value is 2.

[X,U,R] = findop(M,'steady',[1 NaN],2);

Display the summary report.

disp(R);

            SearchMethod: 'auto'
                 WhyStop: 'Near (local) minimum, (norm(g) < tol).'
              Iterations: 12
               FinalCost: 0
    FirstOrderOptimality: 0
            SignalLevels: [1x1 struct]

Find Nonlinear ARX Simulation Snapshot Using Default Initial States

Load the estimation data and estimate a nonlinear ARX model.

load twotankdata;
z = iddata(y,u,1);
M = nlarx(z,[4 3 1]);

Find the simulation snapshot after 10 seconds, assuming initial states of zero.

[X,U] = findop(M,'snapshot',10,z);

Find Nonlinear ARX Simulation Snapshot Using Initial State Specifications

Load the estimation data and estimate a nonlinear ARX model.

load twotankdata;
z = iddata(y,u,1);
M = nlarx(z,[4 3 1]);

Create an initial state vector. The first four states correspond to delayed output values and the final three states correspond to delayed inputs.

X0 = [2;2;2;2;5;5;5];

Find the simulation snapshot after 10 seconds using the specified initial states.

[X,U] = findop(M,'snapshot',10,z,X0);

Input Arguments

`sys` — Nonlinear ARX model
`idnlarx` object

Nonlinear ARX model, specified as an idnlarx object.

`InputLevel` — Steady-state input level
vector

Steady-state input level for computing the operating point, specified as a vector. The length of InputLevel must equal the number of inputs specified in sys.

The optimization algorithm assumes that finite values in InputLevel are fixed input values. Use NaN to specify unknown input signals with initial guesses of 0. The minimum and maximum bounds for all inputs have default values of -Inf and +Inf respectively.

`OutputLevel` — Steady-state output level
vector

Steady-state output level for computing the operating point, specified as a vector. The length of OutputLevel must equal the number of outputs specified in sys.

The values in OutputLevel indicate initial guesses for the optimization algorithm. Use NaN to specify unknown output signals with initial guesses of 0. The minimum and maximum bounds for all outputs have default values of -Inf and +Inf respectively.

`spec` — Operating-point specifications
`operspec` object

Operating-point specifications, such as minimum and maximum input/output constraints and known inputs, specified as anoperspec object.

`T` — Operating point snapshot time
positive scalar

Operating point snapshot time, specified as a positive scalar. The value of T must be in the range [T₀, N*T_s], where N is the number of input samples, T_s is the sample time and T₀ is the input start time (Uin.Tstart).

`Uin` — Snapshot simulation input
`iddata` object | matrix

Snapshot simulation input, specified as one of the following:

Time-domain iddata object with a sample time and input size that matches sys.
Matrix with as many columns as there are input channels. If the matrix has N rows, the input data is assumed to correspond to the time vector (1:N)*sys.Ts.

`X0` — Initial states
column vector

Initial states of the simulation, specified as a column vector with size equal to the number of states in sys. X0 provides the initial conditions at the time corresponding to the first input sample (Uin.Start, if Uin is an iddata object, or sys.Ts if Uin is a double matrix).

For more information about the states of an idnlarx model, see Definition of idnlarx States.

`Options` — Operating point search options
`findopOptions` option set

Operating point search options, specified as a findopOptions option set.

Output Arguments

`X` — Operating point state values
column vector

Operating point state values, returned as a column vector of length equal to the number of model states.

`U` — Operating point input values
column vector

Operating point input values, returned as a column vector of length equal to the number of inputs.

`Report` — Search result summary
structure

Search result summary report, returned as a structure with the following fields:

Field	Description
`SearchMethod`	Search method used for iterative parameter estimation. See `SearchMethod` in `findopOptions` for more information.
`WhyStop`	Search algorithm termination condition.
`Iterations`	Number of estimation iterations performed.
`FinalCost`	Final value of the minimization objective function (sum of the squared errors).
`FirstOrderOptimality`	$\infty$ -norm of the search gradient vector when the search algorithm terminates.
`SignalLevels`	Structure containing the fields `Input` and `Output`, which are the operating point input and output signal levels respectively.

Algorithms