29.04.2014 Views

PDF file - Johannes Kepler University, Linz - JKU

PDF file - Johannes Kepler University, Linz - JKU

PDF file - Johannes Kepler University, Linz - JKU

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.

CHAPTER 2. PRELIMINARIES 32<br />

Lemma 2.15. Assume that A is positive definite, C positive semi-definite, and that if B<br />

does not possess full row-rank then C is positive definite. Then the generalized eigenvalues<br />

λ of ( ( ) ( ( )<br />

A B<br />

T u A 0 u<br />

= λ<br />

B −C)<br />

p 0 S)<br />

p<br />

are contained in the set [<br />

Proof. To solve the problem<br />

−1; 1 − √ 5<br />

2<br />

]<br />

∪<br />

[<br />

1; 1 + √ 5<br />

2<br />

λAu = Au + B T p<br />

λSp = Bu − Cp<br />

]<br />

.<br />

(2.40a)<br />

(2.40b)<br />

we distinguish two cases.<br />

First if λ = 1 then (2.40a) can be fulfilled with p = 0 and (2.40a) by any u with Bu = 0.<br />

Now for λ ≠ 1 we get u = 1<br />

λ−1 A−1 B T p from (2.40a), insert this in (2.40b) and get<br />

(<br />

λ − 1 + √ 5<br />

2<br />

) (<br />

λ − 1 − √ 5<br />

2<br />

)<br />

BA −1 B T p + (λ + 1)(λ − 1)Cp = 0. (2.41)<br />

As both BA −1 B T and C are positive semi-definite and at least one of them is positive<br />

definite, (2.41) can only be fulfilled for<br />

[<br />

λ ∈ −1; 1 − √ ] (<br />

5<br />

∪ 1; 1 + √ ]<br />

5<br />

,<br />

2<br />

2<br />

which completes the proof.<br />

Remark 2.16. If C ≡ 0 then one can even show that λ ∈ {1, (1 ± √ 5/2}. Similar and more<br />

general results can be found e.g. in [RW92], [IRT93], or [SW94].<br />

and<br />

Now we can easily calculate that<br />

and deduce from Lemma 2.15 that<br />

and<br />

‖ ˆK −1 K‖ ˆK = ‖K1/2 ˆK−1 K 1/2 ‖ l2<br />

‖K −1 ˆK‖ ˆK = ‖ ˆK 1/2 K −1 ˆK1/2 ‖ l2 ,<br />

‖K 1/2 ˆK−1 K 1/2 ‖ l2 ≤ 1 + √ 5<br />

,<br />

2<br />

‖ ˆK 1/2 K −1 ˆK1/2 ‖ l2 ≤<br />

2<br />

√<br />

5 − 1<br />

.<br />

Thus we can estimate<br />

( )<br />

cond ˆK−1 ˆK<br />

K ≤ 3 + √ 5<br />

,<br />

2<br />

where the upper bound is clearly independent of h.

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

Saved successfully!

Ooh no, something went wrong!