Model order reduction by frequency weighted optimal Hankel-norm (HNA) 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.
HNA is an absolute error method which tries to minimize
||G-Gr|| = min
H
||V (G-Gr) W|| = min
H
where V and W denote output and input weightings.
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 corresponding to the projection op(V)*G1*op(W),
where G1 denotes 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 info.nu and the number of
Hankel singular values greater than max(tol1, ns*eps*info.hsv(1);
Specifies the computational approach to be used. Valid values corresponding to this key are:
Use the inverse free descriptor system approach.
Use the inversion based standard approach.
Switch automatically to the inverse free descriptor approach in case of badly conditioned feedthrough matrices in V or W. Default method.
LTI model of the left/output frequency weighting. The weighting must be antistable.
||V (G-Gr) . || = min
H
LTI model of the right/input frequency weighting. The weighting must be antistable.
||. (G-Gr) W || = min
H
LTI model of the left/output frequency weighting. The weighting must have only antistable zeros.
||inv(V) (G-Gr) . || = min
H
LTI model of the right/input frequency weighting. The weighting must have only antistable zeros.
||. (G-Gr) inv(W) || = min
H
LTI model of the left/output frequency weighting. The weighting must be stable.
||V (G-Gr) . || = min
H
LTI model of the right/input frequency weighting. The weighting must be stable.
||. (G-Gr) W || = min
H
LTI model of the left/output frequency weighting. The weighting must be minimum-phase.
||V (G-Gr) . || = min
H
LTI model of the right/input frequency weighting. The weighting must be minimum-phase.
||. (G-Gr) W || = min
H
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.
If ’order’ is not specified, tol1 contains the tolerance for determining the order of the reduced model. For model reduction, the recommended value of tol1 is c*info.hsv(1), where c lies in the interval [0.00001, 0.001]. tol1 < 1. If ’order’ is specified, the value of tol1 is ignored.
The tolerance for determining the order of a minimal realization of the ALPHA-stable part of the given model. tol2 <= tol1 < 1. If not specified, ns*eps*info.hsv(1) is chosen.
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.
Approximation Properties:
Algorithm
Uses SLICOT AB09JD by courtesy of
NICONET e.V.
Package: control