Model order reduction by Balanced Stochastic Truncation (BST) method. The aim of model reduction is to find an LTI system Gr of order nr (nr < n) such that the input-output behaviour of Gr approximates the one from original system G.
BST is a relative error method which tries to minimize
-1 ||G (G-Gr)|| = min inf
Inputs
LTI model to be reduced.
The desired order of the resulting reduced order system Gr. If not specified, nr is chosen automatically according to the description of key ’order’.
Optional pairs of keys and values. "key1", value1, "key2", value2
.
Optional struct with keys as field names.
Struct opt can be created directly or
by function options
. opt.key1 = value1, opt.key2 = value2
.
Outputs
Reduced order state-space model.
Struct containing additional information.
The order of the original system G.
The order of the alpha-stable subsystem of the original system G.
The Hankel singular values of the phase system corresponding to the alpha-stable part of the original system G. The ns Hankel singular values are ordered decreasingly.
The order of the alpha-unstable subsystem of both the original system G and the reduced-order system Gr.
The order of the obtained reduced order system Gr.
Option Keys and Values
The desired order of the resulting reduced order system Gr.
If not specified, nr is the sum of NU and the number of
Hankel singular values greater than MAX(TOL1,NS*EPS)
;
nr can be further reduced to ensure that
HSV(NR-NU) > HSV(NR+1-NU)
.
Approximation method for the H-infinity norm. Valid values corresponding to this key are:
Use the square-root Balance & Truncate method.
Use the balancing-free square-root Balance & Truncate method. Default method.
Use the square-root Singular Perturbation Approximation method.
Use the balancing-free square-root Singular Perturbation Approximation method.
Specifies the ALPHA-stability boundary for the eigenvalues of the state dynamics matrix G.A. For a continuous-time system, ALPHA <= 0 is the boundary value for the real parts of eigenvalues, while for a discrete-time system, 0 <= ALPHA <= 1 represents the boundary value for the moduli of eigenvalues. The ALPHA-stability domain does not include the boundary. Default value is 0 for continuous-time systems and 1 for discrete-time systems.
Use [G, beta*I]
as new system G to combine
absolute and relative error methods.
BETA > 0 specifies the absolute/relative error weighting
parameter. A large positive value of BETA favours the
minimization of the absolute approximation error, while a
small value of BETA is appropriate for the minimization
of the relative error.
BETA = 0 means a pure relative error method and can be
used only if rank(G.D) = rows(G.D) which means that
the feedthrough matrice must not be rank-deficient.
Default value is 0.
If ’order’ is not specified, tol1 contains the tolerance for determining the order of reduced system. For model reduction, the recommended value of tol1 lies in the interval [0.00001, 0.001]. tol1 < 1. If tol1 <= 0 on entry, the used default value is tol1 = NS*EPS, where NS is the number of ALPHA-stable eigenvalues of A and EPS is the machine precision. If ’order’ is specified, the value of tol1 is ignored.
The tolerance for determining the order of a minimal realization of the phase system (see METHOD) corresponding to the ALPHA-stable part of the given system. The recommended value is TOL2 = NS*EPS. TOL2 <= TOL1 < 1. This value is used by default if ’tol2’ is not specified or if TOL2 <= 0 on entry.
Boolean indicating whether equilibration (scaling) should be
performed on system G prior to order reduction.
Default value is true if G.scaled == false
and
false if G.scaled == true
.
Note that for MIMO models, proper scaling of both inputs and outputs
is of utmost importance. The input and output scaling can not
be done by the equilibration option or the prescale
function
because these functions perform state transformations only.
Furthermore, signals should not be scaled simply to a certain range.
For all inputs (or outputs), a certain change should be of the same
importance for the model.
BST is often suitable to perform model reduction in order to obtain low order design models for controller synthesis.
Approximation Properties:
Algorithm
Uses SLICOT AB09HD by courtesy of
NICONET e.V.
Package: control