Information Technology Reference
In-Depth Information
Table 3 Iterative SVDDRK-
SLS Algorithm
Iterative SVDDRK-SLS Algorithm:
(input:I q ; A q ; B q ; C q ; D q ; S q ; tol; output:V r q ; Z r q Þ
Switch q
1/*Choose the initial interpolation points*/
s i q for i q ¼ 1tor q
2/*Construction of the matrices V q by the dual rational-Krylov
based method knowing that*/
(a):V r q ¼ Span
A q Þ 1 B q ; ...; ð
A q Þ 1 B q
ð
s 1 q I q
s r q I q
With V r q V r q ¼ I r q
3/*Calculate the garamian matrix of observability for each
subsystem*/
(a)g o q ¼ R 1
0
e tA q C q C q e tA dt
4/*Construction of the matrix Z q by the SVD based method
knowing that*/
(b):Z r q ¼ Q q V r q ð
V r q Þ
5/*While (the relative change in si i
: ðð s i þ 1 s i Þ= s i Þ tol
(a):A r q ¼ Z r q A q V r q
(b):s i q ¼ k i q ð A r q Þ for i q ¼ 1:r q
(c):Construction of the matrix V q by the rational-Krylov based
method knowing that:
(d):V r q ¼ Span
A q Þ 1 B q ; ...; ð
A q Þ 1 B q
ð
s 1 q I q
s r q I q
With Z r q V r q ¼ I r q
(e):*Construction of the matrix Z r q by the SVD based method
knowing that*/
Z r q ¼ Q q V r q ð
V r q Þ
6/*Generate real V r q and Z r q for complex interpolation point*/
if there exist any is i q is not a real number
then V r q ð real Þ ¼ real ð V r q Þ , V r q ð imaginairy Þ ¼ imag ð V r q Þ ,
Z r q ð
real
Þ ¼ real
ð
Z r q Þ
, Z r q ð
imaginairy
Þ ¼ imag
ð
Z r q Þ
,
½V r q ;
rr q ¼QR½V r q ð
real
Þ;
V r q ð
imaginairy
Þ
,
½Z r q ; rr q ¼QR½Z r q ð real Þ; Z r q ð imaginairy Þ
,
end if
67*Parameters of reduced model*/
A r q ¼ Z r q A q V r q ;
B r q ¼ Z r q B q ;
C r q ¼ C q V r q ;
D r q ¼ D q
End Switch
Search WWH ::




Custom Search