20.05.2014 Views

link to my thesis

link to my thesis

link to my thesis

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

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

there is a polynomial solution z(ρ 1 ,ρ 2 )ofA(ρ 1 ,ρ 2 ) T z = 0 where the degree of ρ 1 is<br />

η 1 = 1 and the degree of ρ 2 is η 2 =2.<br />

Using (11.40), the expression A(ρ 1 ,ρ 2 ) T z(ρ 1 ,ρ 2 ) = 0 can be written as:<br />

(<br />

)<br />

M 0 (ρ 2 )+ρ 1 M 1 (ρ 2 )+ρ 2 1 M 2 (ρ 2 )+ρ 3 1 M 3 (ρ 2 ) z(ρ 1 ,ρ 2 ) = 0 (11.47)<br />

where:<br />

M 0 (ρ 2 ) = M 0,0 + ρ 2 M 0,1 + ρ 2 2 M 0,2 + ρ 3 2 M 0,3<br />

M 1 (ρ 2 ) = M 1,0 + ρ 2 M 1,1 + ρ 2 2 M 1,2<br />

M 2 (ρ 2 ) = M 2,0 + ρ 2 M 2,1<br />

M 3 (ρ 2 ) = M 3,0 ,<br />

(11.48)<br />

and where every matrix M i,j is of dimension m × n. Recall again that the dimensions<br />

m and n are smaller than in the previous subsection because no linearization step is<br />

involved here.<br />

Every solution z(ρ 1 ,ρ 2 ) should take the form z 0 (ρ 2 ) − ρ 1 z 1 (ρ 2 ) because η 1 =1. If<br />

such a solution is substituted in<strong>to</strong> (11.47), and the relationships for the coefficients<br />

of the powers of ρ 1 are worked out, we find:<br />

⎛<br />

⎜<br />

⎝<br />

M 0 (ρ 2 ) 0<br />

M 1 (ρ 2 ) M 0 (ρ 2 )<br />

M 2 (ρ 2 ) M 1 (ρ 2 )<br />

M 3 (ρ 2 ) M 2 (ρ 2 )<br />

0 M 3 (ρ 2 )<br />

⎞<br />

⎟<br />

⎠<br />

(<br />

z0 (ρ 2 )<br />

−z 1 (ρ 2 )<br />

)<br />

(<br />

z0 (ρ<br />

= M 1,2 (ρ 2 )<br />

2 )<br />

−z 1 (ρ 2 )<br />

)<br />

= 0 (11.49)<br />

where the dimension of the block matrix M 1,2 (ρ 2 )is(η 1 + Ñ +1)m × (η 1 +1)n =<br />

5m × 2n.<br />

The degree of ρ 2 is η 2 = 2 and therefore z i (ρ 2 ) can be written as z i (ρ 2 ) =<br />

z i,0 + ρ 2 z i,1 + ρ 2 2 z i,2 for i =0, 1. This can be substituted in<strong>to</strong> (11.49). When<br />

the relationships for the coefficients of the powers of ρ 2 are worked out again, the<br />

following system of equations in matrix-vec<strong>to</strong>r form occurs:<br />

⎛<br />

M 1,2 ⎜<br />

⎝<br />

z 0,0<br />

z 0,1<br />

z 0,2<br />

z 1,0<br />

z 1,1<br />

z 1,2<br />

⎞<br />

= 0 (11.50)<br />

⎟<br />

⎠<br />

where the dimension of the matrix M 1,2 is (η 2 + Ñ + 1)(η 1 + Ñ +1) m × (η 1 + 1)(η 2 +

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

Saved successfully!

Ooh no, something went wrong!