Exercícios Resolvidos de Método de Householder

Ver Teoria

Enunciado

Aplique as transformações de Householder à matriz 4   × 4 para calcular A ( 2 ) (primeira iteração).

A = 1 0 0 0 2 3 2 3 0 1 3 6 0 4 5 0

Passo 1

Para a primeira aplicação de uma transformação de Householder temos que α e r , que são:

α = - s g n a k + 1 ,     k k ∑ j = k + 1 n a j k k 2 1 2

r = 1 2 α 2 - a k + 1 ,   k k α

No nosso caso, ainda estamos em k = 1 :

α = - s g n ( a 21 ) ∑ j = 2 n a j 1 2

r = 1 2 ( α 2 - a 21 α )

Onde:

α = - 2 ∑ j = 2 4 a j 1 2 1 2 = - 4

r = 1 2 - 4 2 - 2 - 4 1 / 2 = 12

Passo 2

Para escrever o vetor v 1 fazemos:

v 1 =   0

Para v 2 :

v k + 1 k = a k + 1 , k k - α 2 r

v 2 =   2 - ( - 4 ) 2 . 12 = 12 4

Para o v 3 e o v 4 :

v j k = a j k k 2 r , j = k + 2 ,   k + 3 , … n

v 3 = 0 2 . 12 = 0

v 4 = 0 2 . 12 = 0

Vetor:

v 1 =   0 12 4 0 0 T

Passo 3

Então calculamos:

P k = I - 2 v k v k T

P 1 = I -   2 v 1 v 1 T

P ( 1 ) = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 - 2 0 12 4 0 0 0 12 4 0 0

P ( 1 ) = 1 0 0 0 0 - 0,5 0 0 0 0 1 0 0 0 0 1

Passo 4

Então calculamos:

A ( 2 ) = P 1 A P 1

1 0 0 0 0 - 0,5 0 0 0 0 1 0 0 0 0 1 . 1 0 0 0 2 3 2 3 0 1 3 6 0 4 5 0 . 1 0 0 0 0 - 0,5 0 0 0 0 1 0 0 0 0 1

A ( 2 ) = 1 0 0 0 - 1 0,75 - 1 - 1 , 5 0 - 0 , 5 3 6 0 - 2 5 0

Resposta

A ( 2 ) = 1 0 0 0 - 1 0,75 - 1 - 1,5 0 - 0,5 3 6 0 - 2 5 0