09.05.2013 Views

Autovalores do Laplaciano - Departamento de Matemática - UFMG

Autovalores do Laplaciano - Departamento de Matemática - UFMG

Autovalores do Laplaciano - Departamento de Matemática - UFMG

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.

Rodney Josué Biezuner 115<br />

4.26 Proposição. Sejam A ∈ Mn (R) uma matriz simétrica positiva <strong>de</strong>finida, v ∈ R n e W um subsespaço<br />

<strong>de</strong> R n . Então existe um único w ∈ W tal que<br />

v − w A = min<br />

z∈W v − z A .<br />

O vetor w é caracteriza<strong>do</strong> pela condição v − w ⊥A W .<br />

Segue <strong>de</strong>ste resulta<strong>do</strong> que e j A é minimiza<strong>do</strong> quan<strong>do</strong> escolhemos p = α0p 0 + α1p 1 + . . . + αj−1p j−1 ∈ Wj<br />

tal que e j = e 0 − p satisfaz<br />

e j ⊥A p i para i = 1, . . . , j − 1. (4.72)<br />

Definição. Dois vetores y, z que são ortogonais com respeito ao produto interno 〈·, ·〉 A , isto é, tais que<br />

são chama<strong>do</strong>s conjuga<strong>do</strong>s.<br />

〈y, z〉 A = 0<br />

Nosso objetivo então é <strong>de</strong>senvolver um méto<strong>do</strong> em que o erro a cada passo é conjuga<strong>do</strong> com todas as direções<br />

<strong>de</strong> busca anteriores. O próximo resulta<strong>do</strong>, que é basicamente uma reafirmação da Proposição 4.25, mostra<br />

que em qualquer méto<strong>do</strong> <strong>de</strong> <strong>de</strong>scida em que a busca na reta é exata satisfaz automaticamente e j ⊥A p j−1 ,<br />

isto é, (4.72) é váli<strong>do</strong> para a última iteração (o erro da iteração presente é A-ortogonal à direção <strong>de</strong> busca<br />

da iteração anterior).<br />

4.27 Proposição. Seja x k+1 = x k + αkp k obti<strong>do</strong> através <strong>de</strong> uma busca na reta exata. Então<br />

e<br />

Prova: Temos<br />

r k+1 ⊥ p k<br />

e k+1 ⊥A p k .<br />

b − Ax k+1 = b − Ax k − αkAp k ,<br />

<strong>de</strong> mo<strong>do</strong> que a seqüência <strong>do</strong>s resíduos é dada pela fórmula<br />

Logo,<br />

Além disso, como<br />

k+1 k<br />

r , p = r k+1 , p k k k<br />

− αk Ap , p = r k , p k −<br />

r k+1 = r k − αkAp k . (4.73)<br />

Ae k+1 = r k+1 ,<br />

p k , r k <br />

〈p k , Ap k 〉<br />

Ap k , p k = 0.<br />

segue que k+1 k<br />

e , p <br />

A = Ae k+1 , p k = r k+1 , p k = 0.<br />

<br />

O significa<strong>do</strong> geométrico <strong>de</strong>ste resulta<strong>do</strong> é que o mínimo <strong>do</strong> funcional f na reta xk + αkpk ocorre quan<strong>do</strong> a<br />

<strong>de</strong>rivada direcional <strong>de</strong> f na direção <strong>de</strong> busca é zero, ou seja,<br />

0 = ∂f k+1<br />

x<br />

∂pk<br />

= ∇f x k+1 k+1<br />

, pk = r , pk .<br />

De acor<strong>do</strong> com a Proposição 4.27, <strong>de</strong>pois <strong>do</strong> primeiro passo temos e 1 ⊥A p 0 . Para manter os erros<br />

subseqüentes conjuga<strong>do</strong>s a p 0 , como<br />

e k+1 = x − x k+1 = x − x k − αkp k

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

Saved successfully!

Ooh no, something went wrong!