20.05.2014 Views

link to my thesis

link to my thesis

link to my thesis

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

204 CHAPTER 11. H 2 MODEL-ORDER REDUCTION FROM ORDER N TO N-3<br />

with z 0 (ρ 2 ),z 1 (ρ 2 ),...,z η1 (ρ 2 ) all polynomial in ρ 2 with a maximal degree of η 2 .The<br />

idea is <strong>to</strong> substitute this expression for z(ρ 1 ,ρ 2 ) in<strong>to</strong> (11.22), which leads <strong>to</strong>:<br />

(<br />

) (z0<br />

P (ρ 2 ) T + ρ 1 Q(ρ 2 ) T (ρ 2 ) − ρ 1 z 1 (ρ 2 )+ρ 2 1z 2 (ρ 2 )+...+(−1) η1 ρ η1<br />

1 z η 1<br />

(ρ 2 ) ) =<br />

1<br />

(<br />

·<br />

(<br />

+P (ρ 2 ) T z 0 (ρ 2 )<br />

)<br />

)<br />

+<br />

ρ 1 ·<br />

(<br />

−P (ρ 2 ) T z 1 (ρ 2 ) +Q(ρ 2 ) T z 0 (ρ 2 )<br />

)<br />

+<br />

ρ 2 1 · +P (ρ 2 ) T z 2 (ρ 2 ) −Q(ρ 2 ) T z 1 (ρ 2 ) +<br />

.<br />

.<br />

.<br />

(<br />

)<br />

.<br />

ρ η1<br />

1 ·<br />

(<br />

(−1) η1 P (ρ 2 ) T z η1 (ρ 2 ) (−1) η1−1 Q(ρ 2 ) T z η1−1(ρ 2 )<br />

)<br />

+<br />

ρ η1+1<br />

1 ·<br />

(−1) η1 Q(ρ 2 ) T z η1 (ρ 2 ) =0<br />

(11.24)<br />

When the relationships for the coefficients are worked out, this generates a system<br />

of equations which can be expressed in matrix-vec<strong>to</strong>r form as follows:<br />

⎛<br />

P (ρ 2 ) T 0 ... ... 0<br />

⎞<br />

. ⎛<br />

Q(ρ 2 ) T P (ρ 2 ) T ..<br />

0 Q(ρ 2 ) T . .. . ..<br />

. 0 0 .. . .. . ..<br />

⎜<br />

⎝<br />

⎜<br />

⎝<br />

.<br />

.<br />

.<br />

.. P (ρ2 ) T ⎟<br />

⎠<br />

0 0 ... ... Q(ρ 2 ) T<br />

⎛<br />

M η1,η 2<br />

(ρ 2 )<br />

⎜<br />

⎝<br />

z 0 (ρ 2 )<br />

−z 1 (ρ 2 )<br />

z 2 (ρ 2 )<br />

.<br />

(−1) η1 z η1 (ρ 2 )<br />

z 0 (ρ 2 )<br />

−z 1 (ρ 2 )<br />

z 2 (ρ 2 )<br />

.<br />

.<br />

(−1) η1 z η1 (ρ 2 )<br />

⎞<br />

=<br />

⎟<br />

⎠<br />

⎞<br />

=0<br />

⎟<br />

⎠<br />

(11.25)<br />

where the block matrix M η1,η 2<br />

(ρ 2 ) has dimensions (η 1 +2) m × (η 1 +1) n. The rank<br />

of the matrix M η1,η 2<br />

(ρ 2 ) is denoted by r η1,η 2<br />

and it holds that r η1,η 2<br />

< (η 1 +1)n.<br />

Since the variable η 2 denotes the degree of ρ 2 in z(ρ 1 ,ρ 2 ), z i (ρ 2 ) can be written<br />

as:<br />

z i (ρ 2 )=z i,0 + ρ 2 z i,1 + ρ 2 2 z i,2 + ...+ ρ η2<br />

2 z i,η 2<br />

, (11.26)<br />

for i =0,...,η 1 .<br />

When (11.26) is substituted in<strong>to</strong> the equations of (11.25) and the relationships<br />

for the coefficients are worked out, a system of equations can be set up which only<br />

involves constant coefficients. When P (ρ 2 ) T = B T + ρ 2 D T and Q(ρ 2 ) T = C T are<br />

plugged in<strong>to</strong> the expressions, the resulting system of equations can be written in

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!