Chapter 2 Additive Error Reduction
© National Instruments Corporation 2-19 Xmath Model Reduction Module
and finally:
These four matrices are the constituents of the system matrix of ,
where:
Digression:
This choice is related to the ideas of [Glo84] in the following way;
in [ Glo8 4], the complete set is identified of satisfying
with having a stable part of order ni–1
. The set is parameterized in
terms of a stable transfer function matrix K(s) which has to satisfy
with C2, B2 being two matrices appearing in the course of the algorithm
of [Glo84], and enjoying the property . The particular
choice
in the algorithm of [Glo84] and flagged in corollary 7.3 of [Glo84] is
equivalent to the previous construction, in the sense of yielding the
same , though the actual formulas used here and in [Glo84] for the
construction procedure are quite different. In a number of situations,
including the case of scalar (SISO)G(s), this is the only choice.
The next step of the algorithm is to call stable( ), to separate into
its stable and unstable parts, call them and , stable( ) will
always assign the matrix to , and the final step of the algorithm is
A
˜SB
1–A11 A12–A22
#A21
–()=
B
˜SB
1–B1A12–A22
#B2
–()=
C
˜C1C2A22
#A21
#
–=
D
˜DC
2A22
#B2
–=
G
˜s()
G
˜s() Grs() Gus()+=
G
˜s()
Gjω()G
˜jω()–∞σni
=
G
˜
C2Ks()B2
′
+0=
IK′jω–()Kjω()–0forallω≤
C2
′C2B2B2
′
=
Ks() C2C′2C2
()
#B2
–=
G
˜s
G
˜s()
G
˜s() G
˜us()
D
˜G
˜rs()