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.

10.3. THE KRONECKER CANONICAL FORM OF A MATRIX PENCIL 173<br />

Working out the relationships for the coefficients of the powers of ρ 1 in (10.29),<br />

yields the following in matrix-vec<strong>to</strong>r format:<br />

⎛<br />

⎞<br />

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

⎛<br />

C B<br />

.<br />

0 C . .<br />

. . ..<br />

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

. ..<br />

⎜<br />

⎝<br />

⎜ . . .<br />

⎝ . . .. ⎟ B ⎠<br />

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

z 0<br />

−z 1<br />

z 2<br />

.<br />

(−1) ε z ε<br />

⎞ ⎛<br />

= M ε ⎟ ⎜<br />

⎠ ⎝<br />

z 0<br />

−z 1<br />

z 2<br />

.<br />

(−1) ε z ε<br />

⎞<br />

=0, (10.30)<br />

⎟<br />

⎠<br />

where the block matrix M ε has dimensions (ε +2)m × (ε +1)n and has rank <<br />

(ε +1)n. Note that any non-zero vec<strong>to</strong>r in the kernel of the matrix M ε in (10.30)<br />

yields a corresponding solution z(ρ 1 ) of degree ε in (10.27). Note furthermore that<br />

the matrices M ε only depend on the known matrices B and C.<br />

Now let the (algebraic) rank of a matrix M i be denoted by r i = rank(M i ). Because<br />

ε was chosen as the least possible degree of z(ρ 1 ) for which a solution exists, it holds<br />

for the matrices:<br />

( B<br />

M 0 =<br />

C<br />

⎛<br />

⎞<br />

B 0 ... 0<br />

⎛ ⎞<br />

)<br />

B 0<br />

C B<br />

.<br />

,M 1 = ⎝ C B ⎠ ,...,M ε−1 =<br />

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

, (10.31)<br />

0 C<br />

⎜<br />

⎝<br />

.<br />

.<br />

. .. ⎟ B ⎠<br />

0 0 ... C<br />

that r 0 = n, r 1 =2n, ..., r ε−1 = εn.<br />

Theorem 4 in Chapter 12 of [42] shows that if (10.27) has a solution z(ρ 1 )of<br />

minimal degree ε and ε>0, then (B + ρ 1 C) is strictly equivalent <strong>to</strong>:<br />

(<br />

Lε 0<br />

0 D + ρ 1 E<br />

)<br />

, (10.32)<br />

where the submatrix L ε of dimension ε×(ε+1) exhibits the structure given in (10.18).<br />

Furthermore, it holds that (D + ρ 1 E)z(ρ 1 ) = 0 has no solution z(ρ 1 ) of degree less<br />

then ε.<br />

The eigenvalue problem (B + ρ 1 C)ṽ = 0 is transformed using transformation<br />

matrices V and W as follows:<br />

W −1 (B + ρ 1 C) V (V −1 ṽ)=<br />

(<br />

Lε 0<br />

0 D + ρ 1 E<br />

) ( 0<br />

w =<br />

0<br />

)<br />

. (10.33)

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

Saved successfully!

Ooh no, something went wrong!