View metadata, citation and similar papers at core.ac.uk brought to you by CORE
provided by Institute of Transport Research:Publications
Computational Techniques Based on the Blo ck-Diagonal
Form for Solving Large Systems Mo deling Problems
A. Varga
German Aerospace Research Establishment DLR
Institute for Rob otics and System Dynamics
Ob erpfa enhofen, 82230 Wessling, Germany.
Abstract in which each matrix A is square of order n . In ap-
i i
plications it is generally desirable that the diagonal
The reduction of the state-matrix of a linear time-
blo cks have orders as small as p ossible. Theoretically
invariant state-space mo del to a blo ck-diagonal form
the blo cks cannot b e smaller than the blo cks in the Jor-
by using a state co ordinate transformation is equiva-
dan canonical form of A. However, the computation of
lent with an additive decomp osition of the corresp on-
the Jordan form is generally a dicult and often, very
ding transfer-function matrix. Computationally invol-
ill-conditioned problem. From numerical p oint of view
ved and large storage demanding algorithms for solving
a sound approachisto determine blo cks with the lo-
several systems mo deling problems can b e conveniently
west achievable dimensions by using a well-conditioned
reformulated such that they p erform exclusively on the
transformation matrix. In this case, the application of
low order subsystems corresp onding to the individual
transformation to other matrices will not cause signi -
terms of suitable additive decomp ositions. Imp ortant
cant accuracy losses. Avery ecient and numerically
reductions of b oth the computational e ort and requi-
reliable algorithm for computing the BDF was prop o-
red memory usually by using the reformulated algo-
sed byBavely and Stewart [2]. Besides simplicity and
rithms and thus, their applicability can be extended
economy, the main advantage of this algorithm is its
to handle higher order systems. The pap er presents
abilitytokeep under control the condition number of
several algorithms suitable to p erform eciently on
the transformation matrix. Therefore, with this algo-
additively decomp osed systems. The e ectiveness of
rithm, the BDF can be computed with a prescrib ed
these algorithms for solving large order systems mo de-
accuracy loss. The computational e ort to compute
3 3
ling problems relies on a reliable numerical algorithm
the BDF is approximately 11n 15n op erations and
2
to compute the blo ck-diagonal form of a matrix.
all computations can b e p erformed by using only 2n
storage lo cations.
Consider now a state-space mo del of the form
1. Intro duction
x_ t = Axt+But
3
The use of condensed forms of systems matrices in
y t = Cxt
algorithms for analysis and design of control systems
was also considered by other authors. An overview of
denoted also fA; B ; C g, where x is the n-dimensional
computational techniqes based on the use of conden-
state vector, u is the m-dimensional input vector and
sed forms obtainable by employing orthogonal system
y is the p-dimensional output vector. The reduction
similarity transformations is presented in [10]. The use
of the state matrix A to a BDF by a system similarity
of the Hessenb erg form to enhance the eciency of se-
transformation of the form
veral algorithms for systems analysis and mo deling is