O método das potências

Ver Todos os Exercícios

Teoria

Nós revimos até agora, alguns métodos importantes para encontrar autovalores e autovetores de uma matriz. Mas, vimos também que demanda um tempo maior pra resolver tudo!

Mas relaxa, existem métodos (talvez nem tão triviais assim) que vão permitir encontrar autovalor em módulo e um autovalor associado, de uma matriz A, SEM PRECISAR DESCOBRIR O POLINÔMIO CARACTERÍSTICO!

Um curiosidade sobre essse método, é que ele é responsável pelo ranqueamento de páginas usadas pelo Google.

Maaaaaaaaaaaaaaaaaaaaas, o método serve pra determinar autovalores de módulos grandes ou separados dos demais. É um método da potência é uma técnica iterativa (sim, meus caros... novamente iterações), usada para determinar o autovalor dominante de uma matriz, que basicamente é um autovalor com maior módulo.

Mas a diferença desse método é que ele não encontra SÓ o autovalor, maaaaas, encontra também o autovetor associado.

“O método das potências deriva do fato que as iterações exageram o tamanho relativo dos módulos autovalores”.

Para que o método das potências encontre uma aproximação para o maior autovalor em móduelo de uma matriz A   ϵ   R n × n , podemos supor também que a matriz tem n autovalores λ 1 , λ 2 , … ,   λ n e um conjunto de autovetores associados { v 1 , v 2 , … , v n } e linearmente independente, devemos supor também que A só tenha valor maior em módulo:

| λ 1 | > λ 2 ≥ | λ n | ≥ 0

Vamos definir alguns termos:

y p 1 1 = | | y ( 1 ) | |

Desse modo, x ( 1 ) :

x ( 1 ) = 1 y p 1 1 y ( 1 ) = 1 y p 1 1 A x ( 0 )

Então:

x p 1 1 = 1 = | | x ( 1 ) | |

Agora defina:

y ( 2 ) = A x ( 1 ) = 1 y p 1 1 A 2 x ( 0 )

De modo análogo, defina a sequiencia de vetores { x ( m ) } m = 0 ∞   e { y ( m ) } m = 0 ∞ e uma sequencia de escalares { µ ( m ) } m = 0 ∞ da seguinte forma:

Então vamos para um exemplo para ficar mais claro!

Exemplo

Temos uma matriz A (2 X 2):

Para isso, o exercício dá um vetor inicial :

v ( 0 ) = [ x           y ] t

Então, fazemos o processo de normalização:

u ( x ) = x x = a

u ( y ) = y x = b

u = [ a           b ] t

Agora, multiplicamos a matriz A pela matriz de vetors normalizados:

a 11 a 12 a 21 a 22 a b = x ( 1 )

O teste de convergência :

x ( 1 ) - x ( 0 ) x ( 1 ) < ε

Exemplo

Exercícios Resolvidos

Exercício Resolvido #1

Nós vimos outros métodos de encontrar autovalores e autovetores, mas, que esses métodos demandam um tempo computacional além do necessário em alguns casos. Rascunhe quais sãos as entradas e saídas necessárias para implementar o algorítimo para o método das potências.

Passo 1

A primeira coisa que devemos ter em mente é que para obter uma aproximação o para o autovalor dominante e autovetor associado da matriz A , n × n de um determinado vetor não nulo x.

Passo 2

Como ENTRADA:

Devemos primeiramente colocar ou pedir para que o usuário determine as dimensões n da minha matriz A .

Outro valor importante, é pedir os valores que irão compor a matriz A .

Além disso, devemos pedir o vetor x a tolerância (TOL) e o número máximo de iterações que queremos atingir.

Passo 3

Para Saída: Colocaremos o autovalor aproximado µ, o autovetor aproximado x ou uma mensagem dizendo que o número máximo de iterações foi excedido.

Resposta

ENTRADA: n , A , x , TOL, N.

SAÍDA: µ, x , N excedido.

Exercício Resolvido #2

Verdadeiro e falso:

Marque verdadeiro e falso, e corrija se necessário. Considere A como uma matriz n × n diagonalizável com um autovalor dominante λ .

  1. Seja x 0 = y 0 for um vetor inicial R n o maior componente deve ser 1.
  2. Para a maioria das escolhas de x 0 , m k converge para o autovalor dominante λ e y não converge para um autovalor dominante.
  3. O nome do Método da Potência deriva do fato de que as iterações exageram o tamanho relativo dos módulos dos autovalores.

Passo 1

  1. Seja x 0 = y 0 for um vetor inicial R n o maior componente deve ser 1.
  2. Essa premissa é verdadeira, o vetor inicial é sempre normalizado e transformado em valores vetor

    Passo 2

  3. Para a maioria das escolhas de x 0 , m k converge para o autovalor dominante λ e y não converge para um autovalor dominante.
  4. Essa premissa é falsa.

    O método da potência é uma técnica iterativa, usada para determinar o autovalor DOMINANTE de uma matriz, e converge para um autovalor dominante.

    Passo 2

  5. O nome do Método da Potência deriva do fato de que as iterações exageram o tamanho relativo dos módulos dos autovalores.

Verdadeiro.

Resposta

  1. V
  2. F
  3. V

Exercício Resolvido #3

Resolva a primeira iteração da matriz abaixo, pelo método das potências, considerando o vetor inicial:

x ( 0 ) = 1 1 1 t

A = - 4 14 0 - 5 13 0 - 1 0 2

Passo 1

Primeiro passo é fazer:

A x = y

Se:

A = - 4 14 0 - 5 13 0 - 1 0 2

x ( 0 ) = 1 1 1 t

Então:

- 4 14 0 - 5 13 0 - 1 0 2 . 1 1 1 = - 4 + 14 - 5 + 13 - 1 + 2 = 1 0 8 1

Passo 2

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

10 8 1

Onde µ será o termo maior 10, que servirá para normalizar:

x ( 1 ) = y ( 1 ) 10

x ( 1.1 ) = 10 10 = 1

x ( 1 . 2 ) = 8 10 = 0,8

x ( 1 . 3 ) = 1 10 = 0,1

O vetor x ( 1 ) , passa a ser

x ( 1 ) = 1 0,8 0 , 1 t

Resposta

x ( 1 ) = 1 0,8 0,1 t

Exercício Resolvido #4

Resolva a segunda iteração da matriz abaixo, pelo método das potências, considerando o vetor inicial:

x ( 0 ) = 1 1 t

A = 3 4 2 - 1

Passo 1

Primeira iteração:

A x = y

Se:

A = 3 4 2 - 1

x ( 0 ) = 1 1 t

Então:

3 4 2 - 1 . 1 1 = 3 + 4 2 - 1 = 7 1

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

7 1

Onde µ será o termo maior 7, que servirá para normalizar:

x ( 1 ) = y ( 1 ) 7

x ( 1.1 ) = 7 7 = 1

x ( 1.2 ) = 1 7 = 0 , 1429

O vetor x ( 1 ) , passa a ser

x ( 1 ) = 1 0 , 1429 t

Passo 2

Segunda iteração:

A x = y

Se:

A = 3 4 2 - 1

x ( 1 ) = 1 0,1429 t

Então:

3 4 2 - 1 . 1 0,1429 = 3 + 0,5716 2 - 0,1429 = 3 , 5716 1,8571

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

3,5716 1,8571

Onde µ será o termo maior 3,5716 , que servirá para normalizar:

x ( 2 ) = y ( 2 ) 3,5716

x ( 2.1 ) = 3,5716 3,5716 = 1

x ( 2.2 ) = 1,8571 3,5716 = 0,5199

O vetor x ( 1 ) , passa a ser

x ( 2 ) = 1 0,5199 t

Resposta

Após 2 iterações:

A u t o v a l o r = λ = 3,5716

A u t o v e t o r = x ( 2 ) = 1 0,5199 t

Exercício Resolvido #5

Resolva a segunda iteração da matriz abaixo, pelo método das potências, considerando o vetor inicial e calcule o erro:

x ( 0 ) = 1 1 t

A = 34,4 12,5 2 5,9 67,3

Passo 1

Primeira iteração:

A x = y

Se:

A = 34,4 12,5 25,9 67,3

x ( 0 ) = 1 1 t

Então:

34,4 12,5 25,9 67,3 . 1 1 = 3 4,4 + 12,5 2 5,9 + 67,3 = 46,9 93,2

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

46,9 93,2

Onde µ será o termo maior 93,2 , que servirá para normalizar:

x ( 1 ) = y ( 1 ) 93,2

x ( 1.1 ) = 46,9 93,2 = 0,5032

x ( 1.2 ) = 93,2 93,2 = 1

O vetor x ( 1 ) , passa a ser

x ( 1 ) = 0,5032 1 t

Passo 2

Segunda iteração:

A x = y

Se:

A = 34,4 12,5 25,9 67,3

x ( 1 ) = 0,5032 1 t

Então:

34,4 12,5 25,9 67,3 . 0,5032 1 = 17,31008 + 12,5 1 3 , 0 3288 + 67,3 = 29,81008 80,33288

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

29,81008 80,33288

Onde µ será o termo maior 80,33288 , que servirá para normalizar:

x ( 2 ) = y ( 2 ) 80,33288

x ( 2.1 ) = 29,81008 80,33288 = 3,5774

x ( 2.2 ) = 80,33288 80,33288 = 1

O vetor x ( 1 ) , passa a ser

x ( 2 ) = 3,5774 1 t

O erro será:

ε = 80,33288     - 93,2 80,33288     = 0 , 1608

Resposta

ε = 0,1608

Exercício Resolvido #6

Resolva a primeira iteração da matriz abaixo, pelo método das potências, considerando o vetor inicial e erro ε < 0,1 :

x ( 0 ) = 0,9 0,8 t

A = 2 1 1 2

Passo 1

Primeira iteração:

A x = y

Se:

A = 2 1 1 2

x ( 0 ) = 0,9 0,8 t

Devemos normalizar o vetor:

x ( 1 ) = y ( 1 ) 0,9

x ( 1.1 ) = 0,9 0,9 = 1

x ( 1.2 ) = 0,8 0 , 9 = 0,8 889

Então:

2 1 1 2 . 1 0,8 889 = 2 + 0,8889 1 + 1,7778 = 2,8889 2,7779

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

2,8889 2,7779

Onde µ será o termo maior 2,8889, que servirá para normalizar:

x ( 1 ) = y ( 1 ) 2,8889

x ( 1.1 ) = 2,8889 2,8889 = 1

x ( 1.2 ) = 2,7779 2,8889 = 0 , 9615

O vetor x ( 1 ) , passa a ser

x ( 1 ) = 1 0 , 9615 t

E agora, calculamos o erro:

ε = 2,8889 - 1 2,8889 = 0,6538 > 0,1

Passo 2

Segunda iteração:

A x = y

Se:

A = 2 1 1 2

x ( 1 ) = 1 0,9615 t

Então:

2 1 1 2 . 1 0 , 9615 = 2 + 0 , 9615 1 + 1,923 = 2,9615 2,9231

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

2,9615 2,9231

Onde µ será o termo maior 2,9615 , que servirá para normalizar:

x ( 2 ) = y ( 2 ) 2,9615

x ( 2 . 1 ) = 2,9615 2,9615 = 1

x ( 2 . 2 ) = 2,9231 2,9615 = 0,9 870

O vetor x ( 1 ) , passa a ser

x ( 2 ) = 1 0,9 870 t

E agora, calculamos o erro:

ε = 2,9615 - 2,8889 2,9 615 = 0 , 0245 < 0,1

Resposta

Após 2 iterações:

A u t o v a l o r = λ = 2,9615

A u t o v e t o r = x ( 2 ) = 1 0,9870 t

Exercício Resolvido #7

Ache o autovetor e o autovalor da matriz abaixo, pelo método das potências, considerando o vetor inicial e erro ε < 0,1 :

x ( 0 ) = 1 1 t

A = 4 3 1 2

Passo 1

Primeira iteração:

A x = y

Se:

A = 4 3 1 2

x ( 0 ) = 1 1 t

Então:

4 3 1 2 . 1 1 = 4 + 3 1 + 2 = 7 3

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

7 3

Onde µ será o termo maior 7, que servirá para normalizar:

x ( 1 ) = y ( 1 ) 7

x ( 1.1 ) = 7 7 = 1

x ( 1.2 ) = 3 7 = 0 , 4286

O vetor x ( 1 ) , passa a ser

x ( 1 ) = 1 0,4286 t

E agora, calculamos o erro:

ε = 7 - 1 7 = 0 , 8 5 71 > 0,1

Passo 2

Segunda iteração:

A x = y

Se:

A = 4 3 1 2

x ( 1 ) = 1 0,4286 t

Então:

4 3 1 2 . 1 0 , 4286 = 4 + 1,2858 1 + 0,8572 = 5,2858 1,8572

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

5,2858 1,8572

Onde µ será o termo maior 5,2858 , que servirá para normalizar:

x ( 2 ) = y ( 2 ) 5,2858

x ( 2.1 ) = 5,2858 5,2858 = 1

x ( 2.2 ) = 1,8572 5,2858 = 0 , 3514

O vetor x ( 2 ) , passa a ser

x ( 2 ) = 1 0 , 3514 t

E agora, calculamos o erro:

ε = 5,2858 - 7 5,2858 = 0 , 3243 > 0,1

Passo 3

Segunda iteração:

A x = y

Se:

A = 4 3 1 2

x ( 2 ) = 1 0,3514 t

Então:

4 3 1 2 . 1 0,3514 = 4 + 1 , 0542 1 + 0 , 7028 = 5 , 0542 1 , 7028

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

5,0542 1,7028

Onde µ será o termo maior 5 , 0542 , que servirá para normalizar:

x ( 3 ) = y ( 2 ) 5,0542    

x ( 3 . 1 ) = 5,0542 5,0542 = 1

x ( 3 . 2 ) = 1 , 7028 5,0542 = 0,3 369

O vetor x ( 3 ) , passa a ser

x ( 3 ) = 1 0,3 369 t

E agora, calculamos o erro:

ε = 5,0542 - 5,2858 5,0542 = 0 , 0458 < 0,1

Resposta

Após 3 iterações:

A u t o v a l o r = λ = 5,0542

A u t o v e t o r = x ( 3 ) = 1 0,3369 t

Exercício Resolvido #8

Resolva até a segunda iteração da matriz abaixo, pelo método das potências, considerando o vetor inicial e encontre o erro na ultima iteração:

x ( 0 ) = 1 1 1 t

A = 3 0 1 2 2 2 4 2 5

Passo 1

Primeira iteração) Primeiro passo é fazer:

A x = y

Se:

A = 3 0 1 2 2 2 4 2 5

x ( 0 ) = 1 1 1 t

Então:

3 0 1 2 2 2 4 2 5 . 1 1 1 = 3 + 0 + 1 2 + 2 + 2 4 + 2 + 5 = 4 6 11

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

4 6 11

Onde µ será o termo maior 11 , que servirá para normalizar:

x ( 1 ) = y ( 1 ) 11

x ( 1.1 ) = 4 11 = 0,3637

x ( 1.2 ) = 6 11 = 0,5454

x ( 1.3 ) = 11 11 = 1

O vetor x ( 1 ) , passa a ser

x ( 1 ) = 0,3637 0,5 454 1 t

Passo 2

Primeiro passo é fazer:

A x = y

Se:

A = 3 0 1 2 2 2 4 2 5

x ( 1 ) = 0,3637 0,5454 1 t

Então:

3 0 1 2 2 2 4 2 5 . 0,3637 0,5454 1 = 1,0911 + 0 + 1 0,7274 + 1,0908 + 2 1,4548 + 1,0908 + 5 = 2,0911 3,8182 7,5456

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

2,0911 3,8182 7,5456

Onde µ será o termo maior 7,5456 , que servirá para normalizar:

x ( 2 ) = y ( 2 ) 7,5456

x ( 2 . 1 ) = 2,0911 7,5456 = 2,7712

x ( 2 . 2 ) = 3,8182 7,5456 = 0,506

x ( 2 . 3 ) = 7,5456 7,5456 = 1

O vetor x ( 2 ) , passa a ser:

x ( 2 ) = 2,7712 0,506 1 t

O erro será:

ε = 7,5456 - 11 7,5456 = 0 , 4578

Resposta

ε = 0,4578

Exercício Resolvido #9

Resolva até a segunda iteração da matriz abaixo, pelo método das potências, considerando o vetor inicial e encontre o erro na ultima iteração :

x ( 0 ) = 1 1 1 t

A = 528,2 547,2 156,4 273,8 312,8 98,0 78,2 98,0 39,0

Passo 1

Primeira iteração) Primeiro passo é fazer:

A x = y

Se:

A = 528,2 547,2 156,4 273,8 312,8 98,0 78,2 98,0 39,0

x ( 0 ) = 1 1 1 t

Então:

528,2 547,2 156,4 273,8 312,8 98,0 78,2 98,0 39,0 . 1 1 1 = 528,2 + 547,2 + 156,4 273,8 + 318,8 + 98 72,2 + 98 + 39 = 1231,8 684,6 215,3

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

1231,8 684,6 215,3

Onde µ será o termo maior 1231,8 , que servirá para normalizar:

x ( 1 ) = y ( 1 ) 1231,8

x ( 1.1 ) = 1231,8 1231,8 = 1

x ( 1.2 ) = 684,6 1231,8 = 0,5556

x ( 1.3 ) = 215,3 1231,8 = 0,17478

O vetor x ( 1 ) , passa a ser

x ( 1 ) = 1 0,5556 0,17478 t

Passo 2

Primeiro passo é fazer:

A x = y

Se:

A = 528,2 547,2 156,4 273,8 312,8 98,0 78,2 98,0 39,0

x ( 1 ) = 1 0,5556 0,17478 t

Então:

528,2 547,2 156,4 273,8 312,8 98,0 78,2 98,0 39,0 . 1 0,5556 0,17478 = 528,2 + 304,024 + 27,335 273,8 + 173,8 + 17,128 78,2 + 54,45 + 6,8164 = 859,7 464,7 139,46

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

859,7 464,7 139,46

Onde µ será o termo maior 859,7 , que servirá para normalizar:

x ( 2 ) = y ( 1 ) 859,7

x ( 2.1 ) = 859,7 859,7 = 1

x ( 2.2 ) = 464,7 859,7 = 0,5406

x ( 2.3 ) = 139,46 859,7 = 0,1622

O vetor x ( 2 ) , passa a ser:

x ( 2 ) = 1 0,5406 0,1622 t

Resposta

x ( 2 ) = 1 0,5406 0,1622 t

Exercício Resolvido #10

Resolva até a segunda iteração da matriz abaixo, pelo método das potências, considerando o vetor inicial e encontre o erro na ultima iteração:

x ( 0 ) = 1 1 1 t

A = 2 3 1 0 3 - 1 0 0 1

Passo 1

Primeira iteração) Primeiro passo é fazer:

A x = y

Se:

A = 2 3 1 0 3 - 1 0 0 1

x ( 0 ) = 1 1 1 t

Então:

2 3 1 0 3 - 1 0 0 1 . 1 1 1 = 2 + 3 + 1 0 + 3 - 1 0 + 0 + 1 = 6 2 1

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

6 2 1

Onde µ será o termo maior 6 , que servirá para normalizar:

x ( 1 ) = y ( 1 ) 6

x ( 1.1 ) = 6 6 = 1

x ( 1.2 ) = 2 6 = 0,3334

x ( 1.3 ) = 1 6 = 0 , 1 667

O vetor x ( 1 ) , passa a ser

x ( 1 ) = 1 0,3334 0,1667 t

Passo 2

Primeiro passo é fazer:

A x = y

Se:

A = 2 3 1 0 3 - 1 0 0 1

x ( 1 ) = 1 0,3334 0,1667 t

Então:

2 3 1 0 3 - 1 0 0 1 . 1 0,3334 0 , 1 667 = 2 + 1 , 0 002 + 0 , 1 667 0 + 1,0 002 - 0,1667 0 + 0 + 0,1667 = 3,1669 1,1669 0,1667

No próximo passo, devemos transformar a matriz calculada em uma matriz normalizada .

3,1669 1,1669 0,1667

Onde µ será o termo maior 3,1669 , que servirá para normalizar:

x ( 2 ) = y ( 2 ) 3,1669

x ( 2.1 ) = 3,1669 3,1669 = 1

x ( 2.2 ) = 1,1669 3,1669 = 0 , 3685

x ( 2.3 ) = 0,1667 3,1669 = 0,0526

O vetor x ( 2 ) , passa a ser:

x ( 2 ) = 1 0 , 3685 0,0526 t

O erro será:

ε = 3,1669 - 6 3,1669 = 0 , 8946

Resposta

ε = 0,8946

Exercícios de Livros Relacionados

Considere o seguinte sistema de equações diferenciais, com coeficientes constantes: d y 1 d x = f 1 ( x , y 1 ,y 2 … , y n ) d y 2 d x = f 2 ( x , y 1 ,y 2 … , y n ) d y 3 d x = f 3 ( x , y 1 ,y 2 … ,
A curvatura de uma coluna delgada sujeita a uma carga P pode ser modelada por: d 2 y d x 2 = M E I Onde d 2 y d x 2 especifica a curvatura, M é o momento de curvatura, E é o modulo de elasticidade, e
Seja A uma matriz quadrada de ordem n e sejam λ 1 , λ 2 ,···, λ n seus autovalores. Quais são os autovalores de A-q I onde q é uma constante e I é a matriz identidade?
Para cada uma das matrizes: A = - 2 5 1 - 3 , A = 1 4 3 0 3 1 0 2 - 1 encontre um polinômio que tenha a matriz como raiz.
Sabendo que uma matriz de ordem 3 tem como autovalores λ 1 = - 1 , λ 2 = 2 , λ 3 = 3 .a) Qual é o polinômio característico de A ?b) Quanto vale t r ( A 2 ) ? c) Quais são os autovalores de A - 1 ?d) A