Chapter 2 Additive Error Reduction
© National Instruments Corporation 2-21 Xmath Model Reduction Module
2. Find a stable order ns – 2 approximation Gns–2
of Gns–1
(s), with
3. (Step ns–nr):
Find a stable order nsr approximation of Gnsr+1
,
with
Then, because for ,
for , ..., this being a property of the algorithm, there follows:
The only difference that arises when singular values have multiplicity in
excess of 1 is that the degree reduction at a given step is greater. Thus, if
σns(G) has multiplicity k, then G(s) is approximated by Gns–k(s) of degree
ns –k.
No separate optimization of the D matrix of Gnsr is required. The
approximation error bound is the same as for the first algorithm. The actual
approximation error may be different. Notice that this second algorithm
does not calculate an unstable Gu(s) such that
Discrete-Time SystemsNo special algorithm is presented for discrete-time systems. Any stable
discrete-time transfer-function matrix G(z) can be used to define a stable
continuous-time transfer-function matrix by a bilinear transformation, thus
when α is a positive constant.
Gns 1–jω()Gns 2–jω()–∞σns 1–Gns 1–
()=
.
.
.
Gnsr 1+jω()Gnsr jω()–∞σnsr 1+Gnsr 1+
()=
σiGns 1–
()σ
iG()≤ins<
σiGns 2–
()σ
iGns 1–
()≤
isi–≤
Gjω()Gnsr jω()–σnsr 1+Gnsr 1+
()... σns
++ G()≤
σiG()
insr1+=
ns
∑
≤
Gjω()Gnsr jω()–Gujω()–∞σnsr 1+
=
Hs() Gαs+
αs–
------------
⎝⎠
⎛⎞
=