(Not recommended) Balanced model truncation via Schur method
GRED = schurmr(G) GRED = schurmr(G,order) [GRED,redinfo] = schurmr(G,key1,value1,...) [GRED,redinfo] = schurmr(G,order,key1,value1,...)
schurmr returns a reduced order model
GRED of G and a struct array redinfo containing the error
bound of the reduced model and Hankel singular values of the original system.
The error bound is computed based on Hankel singular values of
G. For a stable system
Hankel singular values indicate the respective state energy of the system. Hence,
reduced order can be directly determined by examining the system Hankel SV's,
With only one input argument
G, the function will show a Hankel
singular value plot of the original model and prompt for model order number to
This method guarantees an error bound on the infinity norm of the additive error ∥
∥∞ for well-conditioned model reduced problems :
This table describes input arguments for
LTI model to be reduced (without any other inputs will plot its Hankel singular values and prompt for reduced order).
(Optional) an integer for the desired order of the reduced model, or optionally a vector packed with desired orders for batch runs
A batch run of a serial of different reduced order models can be generated by
order = x:y, or a vector of integers. By default, all
the anti-stable part of a system is kept, because from control stability point of
view, getting rid of unstable state(s) is dangerous to model a system.
' can be
specified in the same fashion as an alternative for
'. In this case, reduced order will be determined when the sum of
the tails of the Hankel sv's reaches the
A real number or a vector of different errors
Reduce to achieve H∞ error.
Optimal 1x2 cell array of LTI weights
Display Hankel singular plots (default
Integer, vector or cell array
Order of reduced model. Use only if not specified as 2nd argument.
Weights on the original model input and/or output can make the model reduction algorithm focus on some frequency range of interests. But weights have to be stable, minimum phase and invertible.
This table describes output arguments.
LTI reduced order model. Becomes multi-dimensional array when input is a serial of different model order array.
A STRUCT array with 3 fields:
G can be stable or unstable.
GRED can be either continuous or discrete.
Given a continuous or discrete, stable or unstable system,
following commands can get a set of reduced order models based on your
rng(1234,'twister'); G = rss(30,5,4); [g1, redinfo1] = schurmr(G); % display Hankel SV plot % and prompt for order (try 15:20) [g2, redinfo2] = schurmr(G,20); [g3, redinfo3] = schurmr(G,[10:2:18]); [g4, redinfo4] = schurmr(G,'MaxError',[0.01, 0.05]); for i = 1:4 figure(i); eval(['sigma(G,g' num2str(i) ');']); end
Given a state space (A,B,C,D) of a system and k, the desired reduced order, the following steps will produce a similarity transformation to truncate the original state-space system to the kth order reduced model .
Find the controllability and observability Gramians P and Q.
Find the Schur decomposition for PQ in both ascending and descending order, respectively,
Find the left/right orthonormal eigen-bases of PQ associated with the kth big Hankel singular values.
Find the SVD of (VTL,BIG VR,BIG) = U Σ VT
Form the left/right transformation for the final kth order reduced model
SL,BIG = V L,BIG UΣ(1:k,1:k)–½
SR,BIG = VR,BIGVΣ(1:k,1:k)–½
The proof of the Schur balance truncation algorithm can be found in .
 K. Glover, “All Optimal Hankel Norm Approximation of Linear Multivariable Systems, and Their L∝- error Bounds,” Int. J. Control, vol. 39, no. 6, pp. 1145-1193, 1984.
 M. G. Safonov and R. Y. Chiang, “A Schur Method for Balanced Model Reduction,” IEEE Trans. on Automat. Contr., vol. 34, no. 7, July 1989, pp. 729-733.