Autovalores do Laplaciano - Departamento de Matemática - UFMG
Autovalores do Laplaciano - Departamento de Matemática - UFMG
Autovalores do Laplaciano - Departamento de Matemática - UFMG
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