Método de Newton

Ver Todos os Exercícios

Teoria

Introdução

Nós nos acostumamos ao longo da nossa vida acadêmica a trabalhar com sistemas de equações lineares. Isso porque eles têm uma série de propriedades que facilitam muito nossas contas (e vidas :P). O problema é que no mundo real os processos, de forma geral, não ocorrem de forma linear!

Assim, quando estamos solucionando um problema prático, tentamos fazer aproximações de um modelo real para que este se comporte de maneira linear. No entanto, muitas vezes isso não é possível e acabamos tendo que resolver um sistema de equações não lineares.

Mas afinal que diabos é um sistema não linear? É simplesmente um sistema onde as equações não são retas! Por exemplo:

x 1 2 + x 2 2 - 4 = 0 x 1 2 - x 2 2 4 - 1 = 0

Graficamente:

Nesse caso, como podemos ver, temos 4 soluções, onde as equações se encontram!

E pra resolver problemas desse tipo nós vamos ver métodos iterativos para a solução de sistemas não lineares. Mas antes, precisamos de um pouco de base teórica.

Primeira coisa, vamos resolver problemas do tipo F x = 0 . Assim oh:

f 1 x 1 , x 2 , … , x n = 0 f 2 x 1 , x 2 , … , x n = 0 ⋮ f n x 1 , x 2 , … , x n = 0

E vamos usar essa notação aqui:

F x = f 1 x f 2 x ⋮ f n x ; x = x 1 x 2 ⋮ x n

Onde cada função f i x é uma função não linear. Repare que no caso de sistemas lineares temos f i x = a i x + b i .

Agora nós vamos definir o vetor gradiente de f i x denotado por ∇ f i x :

∇ f i x = ∂ f i x ∂ x 1 ∂ f i x ∂ x 2 ⋯ ∂ f i x ∂ x n

E pra definir isso temos que assumir alguns detalhes chatos tipo que F x está definida em um conjunto aberto D ⊂ R n e que tem derivadas contínuas nesse conjunto (é aquele tipo de coisa que muito provavelmente nunca vai interferir em algum exercício, mas temos que mencionar assim mesmo haha).

E ainda vamos que supor que existe pelo menos um ponto x * ∈ D , tal que F x * = 0 (na linguagem dos mortais, estamos assumindo que o sistema tem solução).

Aí a gente acha a matriz Jacobiana definida assim:

J x = ∇ f 1 x ∇ f 2 x ⋮ ∇ f n x = ∂ f 1 x ∂ x 1 ∂ f 1 x ∂ x 2 ⋯ ∂ f 1 x ∂ x n ∂ f 2 x ∂ x 1 ∂ f 2 x ∂ x 2 ⋯ ∂ f 2 x ∂ x n ⋮ ⋮ ⋱ ⋮ ∂ f n x ∂ x 1 ∂ f n x ∂ x 2 ⋯ ∂ f n x ∂ x n

No nosso exemplo temos:

∇ f 1 x = ∂ f 1 x ∂ x 1 ∂ f 1 x ∂ x 2 = 2 x 1 2 x 2

∇ f 2 x = ∂ f 2 x ∂ x 1 ∂ f 2 x ∂ x 2 = 2 x 1 - x 2 2

Logo,

J x = 2 x 1 2 x 2 2 x 1 - x 2 2

E é isso aí! Agora que aprendemos esses conceitos vamos ver como funciona um método desses na prática.

Método de Newton

Então, o método de Newton é o método que é mais usado para resolver sistemas de equações não lineares (é isso aí, se não tem o nome de Gauss vai ter o de Newton, parece que eles competiam haha).

Então, sem mais delongas, o que vamos fazer é o seguinte. Como todo método iterativo, nosso método de Newton vai precisar de um chute inicial x 0 . Nós vamos aplicar esse chute inicial na nossa matriz Jacobiana J x e na nossa F x e depois vamos resolver o sistema:

J x 0 s 0 = - F x 0

Aí o nosso x 1 vai ser obtido:

x 1 = x 0 + s 0

De forma geral, x k + 1 é dado por:

x k + 1 = x k + s k ,

onde s k é solução do sistema:

J x k s k = - F x k

Vamos ver um exemplo. Resolva o sistema não linear F x = 0 tal que:

F x = x 1 + x 2 - 5 x 1 2 + x 2 2 - 25

Usando como aproximação inicial x 0 = - 1 6 T e ε = 0.1

Passo 1: Encontramos a matriz Jacobiana.

J x = ∂ f 1 x ∂ x 1 ∂ f 1 x ∂ x 2 ∂ f 2 x ∂ x 1 ∂ f 2 x ∂ x 2 = 1 1 2 x 1 2 x 2

Passo 2: Calculamos J x 0 e F x 0 :

J x 0 = 1 1 - 2 12

F x 0 = 0 12

Passo 3: Resolvemos o sistema J x 0 s 0 = - F x 0 :

1 1 - 2 12 s 1 0 s 2 0 = 0 - 12

Calculando isso aí você vai encontrar:

s 0 = 0.86 - 0.86 T

Então,

x 1 = x 0 + s 0 = - 0.14 5.14

Passo 4: Critério de Parada

Opa, chegamos numa parte importantíssima que eu não tinha comentado ainda.

Tem dois critérios de parada principais usados nesse método. Eles são, dado um erro ε :

F x k + 1 < ε

ou,

max   ⁡ x k + 1 - x k < ε

Se a questão não impuser nenhum critério você pode usar um desses, ou usar os dois juntos se quiser. Então:

F x 1 = 0 1.44 = 0 2 + 1.44 2 = 1.44 > ε

x 1 1 - x 1 0 = 0.86 > ε

Passo 5: Calculamos J x 1 e F x 1 :

J x 1 = 1 1 - 0.28 10.28

F x 1 = 0 1.44

Passo 6: Resolvemos o sistema J x 1 s 1 = - F x 1 :

1 1 - 0.28 10.28 s 1 1 s 2 1 = 0 - 1.44

Calculando isso aí você vai encontrar:

s 1 = 0.14 - 0.14 T

Então,

x 2 = x 1 + s 1 = 0 5

Passo 7: Critério de Parada

F x 2 = 0 0 = 0 < ε

Que beleza, achamos o resultado exato! Então paramos por aqui!

Já deu pra pegar o espírito da coisa né? O bom de métodos iterativos é que não tem nenhuma surpresa, é sempre repetição!

Vamos dar uma olhada agora no algoritmo.

Algoritmo

Passo 1. Dados iniciais (o que seu programa vai receber de você):

a. Vetor aproximação inicial x = x 1 , x 2 , … x n

b. Erro ε

Passo 2. Definimos um contador de iterações k = 1

Passo 3. Definimos a variável n como o número de elementos em x

Passo 4. Calcule F x e J x onde J x i j = ∂ f i x ∂ x j para 1 ≤ i , j ≤ n

Passo 5. Resolva o sistema linear n × n , J x y = - F x

Passo 6. x = x + y

Passo 7. Se y < ε então retornar x e FIM. Caso contrário,

Passo 8. k = k + 1

Passo 9. Para i variando de 1 a n , defina x i = 1 a i i b i - ∑ j = 1 i - 1 a i j ∙ x j - ∑ j = i + 1 n a i j ∙ x j 0

Passo 10. Se max 1 ≤ i ≤ n ⁡ x i - x i 0 < ε , FIM e retorne x 1 , x 2 , … , x n . Caso contrário, continue para o Passo 5

Passo 11. k = k + 1

Passo 12. Para i variando de 1 a n , defina x i 0 = x i

Passo 13. Volte para o passo 3

Repare que nesse algoritmo o único critério de parada usado é a diferença x k + 1 - x k . Se você quiser pode adicionar o outro critério que comentamos aqui no capítulo também. Outro critério interessante é estipularmos um número máximo de iterações como em todos os outros métodos iterativos.

Observações

Então galera, aqui vão algumas observações interessantes sobre o método de Newton!

Primeiro que para garantir a convergência precisamos de uma aproximação inicial próxima da solução do sistema. Quanto mais próximo o chute, mais rápida vai ser nossa convergência!

Outro ponto fraco do método é a necessidade de resolvermos um sistema matricial a cada iteração. Se esse sistema tiver muitas equações será necessário muito poder computacional (ou manual no caso da prova haha) e o método acaba sendo muito custoso.

Quando J x k é singular (tem o determinante igual a zero), o método de Newton pode não estar definido.

Método de Newton

Algoritmo

Observações

Exercícios Resolvidos

Exercício Resolvido #1

Dado o sistema não linear

x + y = 1 x 2 + 2 y = 2

aplique uma iteração do Método de Newton para encontrar uma solução do sistema, a partir do ponto x 0 , y 0 T = 0,0 T . Calcule os erros absoluto e relativo. O ponto calculado é uma boa aproximação para a solução do sistema?

Passo 1

-

Primeira coisa que temos que fazer é passar as constantes para o outro lado! Esse método é capaz apenas de encontrar o zero de funções não lineares! Assim, temos:

x + y - 1 = 0 x 2 + 2 y - 2 = 0

Onde,

F x = x + y - 1 x 2 + 2 y - 2

A matriz Jacobiana fica:

J x = ∂ f 1 x , y ∂ x ∂ f 1 x , y ∂ y ∂ f 2 x , y ∂ x ∂ f 2 x , y ∂ y = 1 1 2 x 2

Iteração 1:

Agora calculamos J x 0 e F x 0 com x 0 = 0 :

J x 0 = 1 1 0 2

F x 0 = - 1 - 2

Resolvemos o sistema J x 0 s 0 = - F x 0 :

1 1 0 2 s 1 0 s 2 0 = 1 2

Chegando em:

s 0 = 0 1 T

Então,

x 1 = x 0 + s 0 = 0 0 + 0 1 = 0 1

Como podemos ver, o ponto calculado não é apenas uma boa aproximação, ele é uma solução do sistema! Logo, nossos erros absoluto e relativo serão iguais a zero.

Resposta

Ei, a resposta está no passo a passo :)

Exercício Resolvido #2

Use o método de Newton com x 0 = 0 para calcular x 2 nos sistemas não-lineares abaixo:

a 4 x 1 2 - 20 x 1 + 1 4 x 2 2 + 8 = 0 1 2 x 1 x 2 2 + 2 x 1 - 5 x 2 + 8 = 0

b sin ⁡ 4 π x 1 x 2 - 2 x 2 - x 1 = 0 4 π - 1 4 π e 2 x 1 - e + 4 e x 2 2 - 2 e x 1 = 0

c x 1 1 - x 1 + 4 x 2 = 12 x 1 - 2 2 + 2 x 2 - 3 2 = 25

d 5 x 1 2 - x 2 2 = 0 x 2 - 0.25 sin ⁡ x 1 + cos ⁡ x 2 = 0

Passo 1

a Beleza, temos que:

F x = 4 x 1 2 - 20 x 1 + 1 4 x 2 2 + 8 1 2 x 1 x 2 2 + 2 x 1 - 5 x 2 + 8

Primeiro precisamos encontrar a matriz Jacobiana!

J x = ∂ f 1 x ∂ x 1 ∂ f 1 x ∂ x 2 ∂ f 2 x ∂ x 1 ∂ f 2 x ∂ x 2 = 8 x 1 - 20 x 2 2 x 2 2 2 + 2 x 1 x 2 - 5

Iteração 1:

Agora calculamos J x 0 e F x 0 com x 0 = 0 :

J x 0 = - 20 0 2 - 5

F x 0 = 8 8

Resolvemos o sistema J x 0 s 0 = - F x 0 :

- 20 0 2 - 5 s 1 0 s 2 0 = - 8 - 8

Chegando em:

s 0 = 0.4 1.76 T

Então,

x 1 = x 0 + s 0 = 0 0 + 0.4 1.76 = 0.4 1.76

Normalmente agora seria a hora de testarmos o critério de parada, mas como a questão só pede para encontrarmos x 2 vamos continuar!

Iteração 2:

Calculando J x 1 e F x 1 :

J x 1 = - 16.8 0.88 3.55 - 4.3

F x 1 = 1.41 0.62

Resolvemos o sistema J x 1 s 1 = - F x 1 :

- 16.8 0.88 3.55 - 4.3 s 1 1 s 2 1 = - 1.41 - 0.62

Chegando em:

s 1 = 0.096 0.223 T

Então,

x 2 = x 1 + s 1 = 0.4 1.76 + 0.096 0.223 = 0.496 1.983

Passo 2

b Agora temos:

F x = sin ⁡ 4 π x 1 x 2 - 2 x 2 - x 1 4 π - 1 4 π e 2 x 1 - e + 4 e x 2 2 - 2 e x 1

Primeiro encontramos a matriz Jacobiana:

J x = ∂ f 1 x ∂ x 1 ∂ f 1 x ∂ x 2 ∂ f 2 x ∂ x 1 ∂ f 2 x ∂ x 2 = 4 π x 2 cos ⁡ 4 π x 1 x 2 - 1 4 π x 1 cos ⁡ 4 π x 1 x 2 - 2 4 π - 1 2 π e 2 x 1 - 2 e 8 e x 2

Iteração 1:

Agora calculamos J x 0 e F x 0 com x 0 = 0 :

J x 0 = - 1 - 2 - 3.596 0

F x 0 = 0 - 1.582

Resolvemos o sistema J x 0 s 0 = - F x 0 :

- 1 - 2 - 3.596 0 s 1 0 s 2 0 = 0 1.582

Chegando em:

s 0 = - 0.44 0.22 T

Então,

x 1 = x 0 + s 0 = 0 0 + - 0.44 0.22 = - 0.44 0.22

Iteração 2:

Calculando J x 1 e F x 1 :

J x 1 = - 0.041 - 3.919 - 4.673 4.784

F x 1 = - 0.938 0.798

Resolvemos o sistema J x 1 s 1 = - F x 1 :

- 0.041 - 3.919 - 4.673 4.784 s 1 1 s 2 1 = 0.938 - 0.798

Chegando em:

s 1 = - 0.073 - 0.239 T

Então,

x 2 = x 1 + s 1 = - 0.44 0.22 + - 0.073 - 0.239 = - 0.513 - 0.019

Passo 3

c Nessa letra vamos ter que jogar as constantes pro outro lado, já que só sabemos achar o zero dessas funções! Então ficamos com:

F x = x 1 1 - x 1 + 4 x 2 - 12 x 1 - 2 2 + 2 x 2 - 3 2 - 25

Agora encontramos a matriz Jacobiana:

J x = ∂ f 1 x ∂ x 1 ∂ f 1 x ∂ x 2 ∂ f 2 x ∂ x 1 ∂ f 2 x ∂ x 2 = 1 - 2 x 4 2 x - 2 4 2 y - 3

Iteração 1:

Agora calculamos J x 0 e F x 0 com x 0 = 0 :

J x 0 = 1 4 - 4 - 12

F x 0 = - 12 - 12

Resolvemos o sistema J x 0 s 0 = - F x 0 :

1 4 - 4 - 12 s 1 0 s 2 0 = 12 12

Chegando em:

s 0 = - 48 15 T

Então,

x 1 = x 0 + s 0 = 0 0 + - 48 15 = - 48 15

Iteração 2:

Calculando J x 1 e F x 1 :

J x 1 = 97 4 - 100 108

F x 1 = - 2304 3204

Resolvemos o sistema J x 1 s 1 = - F x 1 :

97 4 - 100 108 s 1 1 s 2 1 = 2304 - 3204

Chegando em:

s 1 = 24.057 - 7.391 T

Então,

x 2 = x 1 + s 1 = - 48 15 + 24.057 - 7.391 = - 23.943 7.609

Passo 4

d Temos:

F x = 5 x 1 2 - x 2 2 x 2 - 0.25 sin ⁡ x 1 + cos ⁡ x 2

A matriz Jacobiana será dada por:

J x = ∂ f 1 x ∂ x 1 ∂ f 1 x ∂ x 2 ∂ f 2 x ∂ x 1 ∂ f 2 x ∂ x 2 = 10 x 1 - 2 x 2 - 0.25 cos ⁡ x 1 0.25 sin ⁡ x 2 + 1

Iteração 1:

Agora calculamos J x 0 e F x 0 com x 0 = 0 :

J x 0 = 0 0 - 0.25 1

F x 0 = 0 - 0.25

Opa, temos um problema aqui né? Não é possível obtermos x 1 pois demos de cara com um sistema de uma equação e duas variáveis! Quando uma linha da matriz Jacobiana zera, dizemos que ela é singular. Logo, não podemos continuar o exercício, pois J 0 é singular.

Resposta

a x 2 = 0.496 1.983

b x 2 = - 0.513 - 0.019

c x 2 = - 23.943 7.609

d Não podemos encontrar x 1 pois J 0 é singular.

Exercício Resolvido #3

Use o método de Newton com x 0 = 0 para calcular x 2 nos sistemas não-lineares abaixo:

a x 1 2 + x 2 - 37 = 0 x 1 - x 2 2 - 5 = 0 x 1 + x 2 + x 3 - 3 = 0

b 15 x 1 + x 2 2 - 4 x 3 = 13 x 1 2 + 10 x 2 - x 3 = 11 x 3 3 - 25 x 3 = - 22

Passo 1

a Beleza, temos que:

F x = x 1 2 + x 2 - 37 x 1 - x 2 2 - 5 x 1 + x 2 + x 3 - 3

Primeiro precisamos encontrar a matriz Jacobiana!

J x = ∂ f 1 x ∂ x 1 ∂ f 1 x ∂ x 2 ∂ f 1 x ∂ x 3 ∂ f 2 x ∂ x 1 ∂ f 2 x ∂ x 2 ∂ f 2 x ∂ x 3 ∂ f 3 x ∂ x 1 ∂ f 3 x ∂ x 2 ∂ f 3 x ∂ x 3 = 2 x 1 1 0 1 - 2 x 2 0 1 1 1

Iteração 1:

Agora calculamos J x 0 e F x 0 com x 0 = 0 :

J x 0 = 0 1 0 1 0 0 1 1 1

F x 0 = - 37 - 5 - 3

Resolvemos o sistema J x 0 s 0 = - F x 0 :

0 1 0 1 0 0 1 1 1 s 1 0 s 2 0 s 3 0 = 37 5 3

Chegando em:

s 0 = 5 37 - 39

Então,

x 1 = x 0 + s 0 = 0 0 0 + 5 37 - 39 = 5 37 - 39

Normalmente agora seria a hora de testarmos o critério de parada, mas como a questão só pede para encontrarmos x 2 vamos continuar!

Iteração 2:

Calculando J x 1 e F x 1 :

J x 1 = 10 1 0 1 - 74 0 1 1 1

F x 1 = 25 - 1369 0

Resolvemos o sistema J x 1 s 1 = - F x 1 :

10 1 0 1 - 74 0 1 1 1 s 1 1 s 2 1 s 3 1 = - 25 1369 0

Chegando em:

s 1 = - 0.649 - 18.509 19.158

Então,

x 2 = x 1 + s 1 = 5 37 - 39 + - 0.649 - 18.509 19.158 = 4.351 18.491 - 19.842

Passo 2

c Nessa letra vamos ter que jogar as constantes pro outro lado, já que só sabemos achar o zero dessas funções! Então ficamos com:

F x = 15 x 1 + x 2 2 - 4 x 3 - 13 x 1 2 + 10 x 2 - x 3 - 11 x 3 3 - 25 x 3 + 22

Agora encontramos a matriz Jacobiana:

J x = ∂ f 1 x ∂ x 1 ∂ f 1 x ∂ x 2 ∂ f 1 x ∂ x 3 ∂ f 2 x ∂ x 1 ∂ f 2 x ∂ x 2 ∂ f 2 x ∂ x 3 ∂ f 3 x ∂ x 1 ∂ f 3 x ∂ x 2 ∂ f 3 x ∂ x 3 = 15 2 x 2 - 4 2 x 1 10 - 1 0 0 3 x 3 2 - 25

Iteração 1:

Agora calculamos J x 0 e F x 0 com x 0 = 0 :

J x 0 = 15 0 - 4 0 10 - 1 0 0 - 25

F x 0 = - 13 - 11 22

Resolvemos o sistema J x 0 s 0 = - F x 0 :

15 0 - 4 0 10 - 1 0 0 - 25 s 1 0 s 2 0 s 3 0 = 13 11 - 22

Chegando em:

s 0 = 1.101 1.188 0.880

Então,

x 1 = x 0 + s 0 = 0 0 0 + 1.101 1.188 0.880 = 1.101 1.188 0.880

Iteração 2:

Calculando J x 1 e F x 1 :

J x 1 = 15 2.376 - 4 2.202 10 - 1 0 0 - 22.677

F x 1 = 1.406 1.212 0.681

Resolvemos o sistema J x 1 s 1 = - F x 1 :

15 2.376 - 4 2.202 10 - 1 0 0 - 22.677 s 1 1 s 2 1 s 3 1 = - 1.406 - 1.212 - 0.681

Chegando em:

s 1 = - 0.069 - 0.103 0.030

Então,

x 2 = x 1 + s 1 = 1.101 1.188 0.880 + - 0.069 - 0.103 0.030 = 1.032 1.085 0.910

Resposta

a x 2 = 4.351 18.491 - 19.842

b x 2 = 1.032 1.085 0.910

Exercício Resolvido #4

Se aplicarmos o método de Newton para sistemas não lineares para resolver um sistema linear A x = b ,   A ∈ R n × n ,   x , b ∈ R n , em quantas iterações o método irá convergir? Justifique sua resposta.

Passo 1

-

Eaí, alguma ideia? É, eu também não :p

Bom, só tem um jeito de descobrir então né! Vamos aplicar o método em um sistema linear e ver no que dá:

α 11 α 12 ⋯ α 1 n α 21 α 22 ⋯ α 2 n ⋮ ⋮ ⋱ ⋮ α n 1 α n 2 ⋯ α n n x 1 x 2 ⋮ x n = b 1 b 2 ⋮ b n

Isso nos deixa com um sistema do tipo:

α 11 x 1 + α 12 x 2 + … + α 1 n x n - b 1 = 0 α 2 1 x 1 + α 2 2 x 2 + … + α 2 n x n - b 2 = 0 ⋮ α n 1 x 1 + α n 2 x 2 + … + α n n x n - b n = 0

Onde:

F x = α 11 x 1 + α 12 x 2 + … + α 1 n x n - b 1 α 21 x 1 + α 22 x 2 + … + α 2 n x n - b 2 ⋮ α n 1 x 1 + α n 2 x 2 + … + α n n x n - b n

A matriz Jacobiana ia ficar:

α 11 α 12 ⋯ α 1 n α 21 α 22 ⋯ α 2 n ⋮ ⋮ ⋱ ⋮ α n 1 α n 2 ⋯ α n n

E para qualquer aproximação x 0 que usarmos, teremos:

F x 0 = α 11 α 12 ⋯ α 1 n α 21 α 22 ⋯ α 2 n ⋮ ⋮ ⋱ ⋮ α n 1 α n 2 ⋯ α n n x 1 0 x 2 0 ⋮ x n 0 - b 1 b 2 ⋮ b n

Assim, o sistema J x 0 s 0 = - F x 0 vai ficar:

α 11 α 12 ⋯ α 1 n α 21 α 22 ⋯ α 2 n ⋮ ⋮ ⋱ ⋮ α n 1 α n 2 ⋯ α n n s 1 0 s 2 0 ⋮ s n 0 = b 1 b 2 ⋮ b n - α 11 α 12 ⋯ α 1 n α 21 α 22 ⋯ α 2 n ⋮ ⋮ ⋱ ⋮ α n 1 α n 2 ⋯ α n n x 1 0 x 2 0 ⋮ x n 0

Se multiplicarmos J - 1 x 0 nos dois lados:

s 1 0 s 2 0 ⋮ s n 0 = α 11 α 12 ⋯ α 1 n α 21 α 22 ⋯ α 2 n ⋮ ⋮ ⋱ ⋮ α n 1 α n 2 ⋯ α n n - 1 b 1 b 2 ⋮ b n - x 1 0 x 2 0 ⋮ x n 0

Lembrando que α 11 α 12 ⋯ α 1 n α 21 α 22 ⋯ α 2 n ⋮ ⋮ ⋱ ⋮ α n 1 α n 2 ⋯ α n n - 1 b 1 b 2 ⋮ b n é a própria solução do sistema linear! Então:

s 1 0 s 2 0 ⋮ s n 0 = x 1 x 2 ⋮ x n - x 1 0 x 2 0 ⋮ x n 0

Logo, teremos que:

x 1 = x 0 + s 0 = x 1 0 x 2 0 ⋮ x n 0 + x 1 x 2 ⋮ x n - x 1 0 x 2 0 ⋮ x n 0 = x 1 x 2 ⋮ x n

Ou seja, sempre encontraremos a solução na primeira iteração!

Resposta

Ei, a resposta está no passo a passo :)

Exercício Resolvido #5

Considere o sistema não linear

e x 1 - 1 = 0 e x 2 - 1 = 0

que tem como solução x 1 = x 2 = 0 .

a Verifique que a matriz Jacobiana é invertível em x = 0,0 T

b Resolva esse sistema através do método de Newton com x 0 = - 10 , - 10 T

Passo 1

a Calculando a matriz Jacobiana temos:

J x = ∂ f 1 x ∂ x 1 ∂ f 1 x ∂ x 2 ∂ f 2 x ∂ x 1 ∂ f 2 x ∂ x 2 = e x 1 0 0 e x 2

Em 0,0 T , temos:

J = 1 0 0 1

Que não só é invertível, mas tem como inversa ela mesma, já que é a matriz identidade. Então está verificado haha.

Passo 2

b Iteração 1:

Calculando J x 0 e F x 0 com x 0 = - 10 , - 10 T :

J x 0 = e - 10 0 0 e - 10

F x 0 = e - 10 - 1 e - 10 - 1

Resolvemos o sistema J x 0 s 0 = - F x 0 :

e - 10 0 0 e - 10 s 1 0 s 2 0 = 1 - e - 10 1 - e - 10

Chegando em:

s 0 = e 10 - 1 e 10 - 1 T

Então,

x 1 = x 0 + s 0 = - 10 - 10 + e 10 - 1 e 10 - 1 = e 10 - 11 e 10 - 11

Que é um número estupidamente grande né? Isso aconteceu porque a matriz Jacobiana avaliada no ponto inicial - 10 , - 10 era praticamente igual a matriz nula né? Eu diria que e - 10 é bem próximo de zero haha.

E lembra que quando a matriz Jacobiana se iguala a zero temos uma matriz singular e não podemos continuar com o método? Então, foi exatamente isso que aconteceu aqui ;)

Resposta

Não podemos resolver o problema, pois J - 10 , - 10 é muito próximo da matriz nula.

Exercícios de Livros Relacionados

Resolva os sistemas não lineares abaixo usando seus programas para os métodos de Newton e Newton Modificado com e = 1 o-4:
Resolva os sistemas não lineares abaixo usando seus programas pata os métodos de Newton e Newton Modificado com e = 1 o-4:
O método de Newton pode ser aplicado para a resolução de um sistema linear Ax = b. Neste caso, quantas iterações serão realizadas? Por quê?
Resolva os sistemas não lineares abaixo usando seus programas pata os métodos de Newton e Newton Modificado com e = 1 o-4:
Resolva os sistemas não lineares abaixo usando seus programas pata os métodos de Newton e Newton Modificado com e = 1 o-4: