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.

178 CHAPTER 10. H 2 MODEL-ORDER REDUCTION FROM ORDER N TO N-2<br />

Formally, when the regular pencil D r + ρ 1 E r in (10.17) attains the canonical form<br />

(10.20), then the Kronecker canonical form of the pencil B + ρ 1 C is obtained.<br />

The Jordan blocks N j in (10.20) are square blocks of possibly different dimensions<br />

k j × k j with a structure as in (10.21). The equations resulting from N j w 1 = 0, where<br />

w 1 is the corresponding part of the eigenvec<strong>to</strong>r w, are:<br />

⎧<br />

w 1,1 + ρ 1 w 1,2 =0<br />

⎪⎨<br />

w 1,2 + ρ 1 w 1,3 =0<br />

.<br />

w 1,kj−1 + ρ 1 w 1,kj =0<br />

⎪⎩<br />

w 1,kj =0<br />

These equations only admit the trivial solution w 1 = 0. Therefore, ρ 1 = ∞ does<br />

not play a role in our problem. The blocks N j are said <strong>to</strong> generate infinite eigenvalues.<br />

The number of infinite eigenvalues amounts <strong>to</strong> the sum of the associated dimensions<br />

k j for these blocks.<br />

The Jordan blocks J j in (10.20) are square blocks of possibly different dimensions<br />

l j × l j with a structure as in (10.22). The corresponding equations following from<br />

J j w 1 = 0, are: ⎧⎪ ⎨<br />

⎪ ⎩<br />

ρ 1 w 1,1 + w 1,2 =0<br />

ρ 1 w 1,2 + w 1,3 =0<br />

.<br />

ρ 1 w 1,lj−1 + w 1,lj =0<br />

ρ 1 w 1,lj =0<br />

In this case there are two possible situations: if ρ 1 0, then J j w 1 only admits the<br />

trivial solution. If ρ 1 = 0 it admits the non-trivial solution w 1 =(1, 0,...,0) T ,orany<br />

vec<strong>to</strong>r proportional <strong>to</strong> it. This block corresponds <strong>to</strong> a zero eigenvalue with algebraic<br />

multiplicity l j .<br />

The Jordan blocks corresponding <strong>to</strong> non-zero finite eigenvalues are denoted by<br />

the R j blocks of dimension m j × m j in Equation (10.20). They exhibit the structure<br />

as in (10.23). The finite eigenvalues have the values ρ 1 = α j with multiplicity equal<br />

<strong>to</strong> the sum of the associated dimensions m j , for each value of α j . The corresponding<br />

equations resulting from R j w 1 = 0 are:<br />

⎧<br />

(ρ 1 − α j )w 1,1 + w 1,2 =0<br />

⎪⎨<br />

(ρ 1 − α j )w 1,2 + w 1,3 =0<br />

.<br />

(ρ 1 − α j )w 1,mj−1 + w 1,mj =0<br />

⎪⎩<br />

(ρ 1 − α j )w 1,mj =0

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

Saved successfully!

Ooh no, something went wrong!