Método do Ponto Fixo (MPF)

Ver Todos os Exercícios

Teoria

Introdução

Opa, chega mais!! Vamos começar mais um tópico da disciplina de Cálculo Numérico.

Se eu te pedisse pra encontrar uma raíz da função f x = x ³ - x - 1 , no intervalo [ 1,2 ] você saberia me dizer com facilidade? Caso não saiba, fique tranquilo!! Iremos mostrar aqui um novo método para encontrar esse zero da função: o Método do Ponto Fixo.

Nem sempre é fácil achar os zeros de uma função, às vezes as técnicas que aprendemos são ineficientes e acabamos dando voltas e mais voltas sem conseguir resolver o problema e encontrar as raízes da função.

Muitas vezes quando nos deparamos com um problema difícil de resolver nós substituímos ele por outro equivalente mais fácil de resolver. E é exatamente isso que vamos fazer neste tópico! Substituir o problema de encontrar os zeros da função por outro problema: Encontrar o ponto fixo de uma outra função! Calma que já já vamos ver como funciona esse macete!

Ponto Fixo de uma função qualquer

Primeiro vamos entender o que é o tal do Ponto Fixo de uma função. O ponto fixo p de uma função g x é simplesmente um ponto tal que g p = p . Ou seja, é um ponto onde a função g x encontra a reta x = y , de modo que o valor da função em um ponto vale ele mesmo.

Por exemplo:

A função g ( x ) = 3 x tem como ponto fixo o ponto x = 0 porque g ( 0 ) = 3.0 = 0

A função g ( x ) = 2 x - 1 tem como ponto fixo x = 1   porque g ( 1 ) = 2.1 - 1 = 1

AAAAh, e nada impede de uma função ter mais de 1 ponto fixo. Por exemplo:

A função g ( x ) = 4 x ² tem como ponto fixo x = 0   e x = 1 4

Ok, mas o que isso tem a ver com os zeros de uma função? Vem comigo que já vou explicar...

Método do Ponto Fixo

Imagina que queremos encontrar as raízes de uma função f ( x ) , mas não conseguimos encontrar a resposta igualando a função a zero…

Para resolver esse problema podemos aplicar o Método do Ponto Fixo. O método simplesmente quer transformar a equação f x = 0   em uma equação equivalente x = g x .

Ok ... um pouco mais devagar... O que está acontecendo aqui? Vamos tentar entender aos poucos, beleza?!

A partir da expressão f ( x ) = 0 (que é nosso problema inicial) podemos determinar uma nova expressão g ( x ) = x e encontrar as raízes da equação.

Note que o x   que é raiz de f ( x ) e ponto fixo da g ( x ) !

Essa função g ( x ) é chamada de função de iteração! E ela tem essa cara:

x k + 1 = g x k

Ou seja, o valor do x atual ( x k + 1 ) depende do valor da função anterior g x k .

Mas como funciona na prática essa g ( x ) ? Vamos ver um exemplo:

Exemplo

Começando com x 0 = 1 , calcule a função g x = x + 1 3 e encontre os próximos valores aplicados no ponto.

Vamos começar calculando o x 0 = 1 da função:

g x = x + 1 3

x 1 = g x 0 = 1 + 1 3 = 1,25992105

Agora, utilize o valor que você encontrou para aplicar novamente na função

x 2 = g x 1 = x 1 + 1 3 = 1,312293873

x 3 = g x 2 = x 2 + 1 3 = 1,32235381

x 4 = g x 3 = x 3 + 1 3 = 1,324268745

Percebe que os valores cada vez mais estão se aproximando? Convergindo para um número? Você está observando o processo de convergência. Isso porque x = 1,324   é o ponto fixo dessa função.

Aí você pode me perguntar: Beleza, mas que que isso tudo tem haver com os zeros da função f x   que queríamos calcular desde o começo?

A resposta é que quanto mais próximo o valor de x em f ( x ) chegar no zero da função, mas próximo o valor de x chega do ponto fixo de g ( x ) !

Ou seja, quanto mais próximo a gente estiver de encontrar o ponto fixo de g ( x ) também estaremos mais próximos de encontrar o valor do zero da função.

Lindo, não?! hahaha

Tá, mas e como encontrar a função g ( x ) a partir da f ( x ) ?

Bom, existem infinitas funções de iteração g x para a equação f x = 0 . O que elas possuem em comum é a forma geral:

g x = x + A x f x

onde, seja x * raiz de f , f x * = 0 e A x * ≠ 0 .

Vamos tentar um exemplo prático?!

Exemplo

Seja f x = x ³ - x - 1 . Queremos encontrar uma raiz no intervalo [ 1,2 ] . Tome x 0 = 1 .

Vamos lá! Precisamos encontrar nossa g x . Temos que f x = 0 , pois queremos uma raíz. Podemos isolar o x nessa equação :

x ³ - x - 1 = 0

x ³ = x + 1

x = x + 1 3

Lembra que a carinha da função de iteração é x k + 1 = g x k ?! Então conseguimos encontrar a nossa g x !!

Já calculamos seu ponto fixo anteriormente:

x 1 = g x 0 = 1 + 1 3 = 1,25992105

x 2 = g x 1 = x 1 + 1 3 = 1,312293873

x 3 = g x 2 = x 2 + 1 3 = 1,32235381

x 4 = g x 3 = x 3 + 1 3 = 1,324268745

Se o ponto fixo então é x = 1,324 , então encontramos a raíz procurada! Show de bola?!

Assim, a raiz de f x = x ³ - x - 1 no intervalo 1,2   é x = 1,324 .

Estudo da Convergência

Opa, ali em cima tava escrito infinitas funções?

Seria ótimo se todas as funções de interação funcionassem e convergissem para um ponto fixo né? BUT, infelizmente isso não acontece sempre.

O problema aqui é que nem toda g x vai convergir para a raíz, por mais que ela tenha a raíz como ponto fixo! Vou te dar um exemplo:

A equação x 3 + x - 10 = 0 tem como raiz x = 2 . Há duas funções de iteração possíveis para ela:

g 1 x = 10 - x 3           e             g 2 x = 10 - x 3

Vamos analisar primeiro g 1 x = 10 - x 3 , usando chute inicial x 0 = 1.5 :

x 1 = g x 0 = 6.625

x 2 = g x 1 = - 280.775

x 3 = g x 2 = 22134887.375

Bom, definitivamente não está convergindo para o 2 né?!

Agora vamos ver com g 2 x = 10 - x 3 e novamente x 0 = 1.5 :

x 1 = g x 0 = 2.04083

x 2 = g x 1 = 1.99659

x 3 = g x 2 = 2.00028

x 4 = g x 3 = 1.99998

Agora sim né?! Estamos claramente convergindo para o 2 !

E agora? Como vamos saber se escolhemos certo? Temos que começar a testar? Não! Por sorte alguém inventou um teorema que nos dá condições suficientes para que o processo seja convergente. Precisamos do seguinte:

Seja x * uma raiz da equação f x = 0 , isolada num intervalo I centrado em x * e seja g x uma função de iteração para a equação f x = 0 .

Se

g x e g ' x são contínuas em I ;

g ' x ≤ k < 1 , para todo x ∈ I ,

então o processo iterativo x k + 1 = g x k converge para x * .

Outro detalhe interessante é que, satisfeitas as condições acima, quanto menor for g ' x no ponto fixo de g x mais rápido esse método irá convergir para a raíz! Então, a melhor função de iteração será a que tiver o menor g ' p , onde p é ponto fixo de g x .Vamos verificar esse teorema? Vamos testar a g 1 x do nosso exemplo:

g 1 x = 10 - x 3

g 1 ' x = - 3 x 2

Ambas são contínuas em R então não teríamos problemas com algum intervalo I .

g 1 ' x < 1 ⇒ 3 x 2 < 1 ⇒ x 2 < 1 3 ⇒ - 1 3 < x < 1 3

Opa, repare que é impossível termos um intervalo I centrado em 2 com g 1 ' x < 1 . Por isso ela não convergiu para a raíz!

Vamos verificar a g 2 x :

g 2 x = 10 - x 3

g 2 ' x = - 1 3 10 - x 2 3

Temos que g 2 x é contínua em R e g 2 ' x é contínua em todos os reais menos o 10 . Então ainda é possível termos um intervalo I centrado em 2 , basta que não englobe o 10 .

g 2 ' x < 1

⇒   1 3 10 - x 2 3 < 1

⇒ 1 < 3 10 - x 2 3

⇒ 1 3 3 2 < 10 - x

⇒ x < 9.8

Aqui nossa única restrição é que x < 9.8 e x ≠ 10 . Então existem vários intervalos centrados em 2 possíveis que contenham x 0 = 1.5 e satisfaçam as condições acima, como por exemplo 1,3 . E foi por isso que ela convergiu, belezinha?! ;)

Critério de Parada

Repare que esse método não usa mais um intervalo, mas um chute inicial! Ou seja, não faria mais sentido usarmos b - a < ε como critério de parada! O que usaremos agora será a variação do valor de x entre duas iterações. Dessa forma, nosso critério de parada será:

x k + 1 - x k < ε

Isso normalmente indica que estamos bem perto da raiz. No entanto em alguns casos esse critério pode ser “enganado” pelo comportamento da função, então é bom ficar de olho também no valor da função f x k . Esse inclusive é outro critério de parada usado nesse método:

f x k < ε

Algoritmo

Chegamos na parte do Algoritmo novamente, onde eu escrevo o procedimento da maneira mais simples possível para ser passada para alguma linguagem. Vamos lá então:

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

a. Aproximação inicial x 0

b. Precisão do resultado final ε 1 e ε 2

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

Passo 3. Se f x 0 < ε 1 , escolha x 0 como x _ e FIM. Caso contrário, continue para o Passo 4.

Passo 4. k = k + 1 , aumentamos uma iteração no contador.

Passo 5. x = g x 0

Passo 6. Se x - x 0 < ε 2 , escolha x como x _ e FIM. Caso contrário continue para o Passo 7.

Passo 7. x 0 = x

Passo 8. Voltamos para o Passo 3 .

Como sempre, um bom critério de parada para se adicionar também é um número máximo de iterações.

O chato aqui é que para cada função que formos encontrar o zero teremos que definí-la e definir também uma g x para ela no nosso programa, ou seja, mais trabalho pra gente... mas é a vida! ahahahaha

Agora vamos exercitar! Mas repito, o mais importante aqui é que o conceito tenha sido bem entendido!!! 😊

Método do Ponto Fixo

Exemplo

Exemplo

Estudo da Convergência

Critério de Parada

Algoritmo

Exercícios Resolvidos

Exercício Resolvido #1

Use o método do Ponto Fixo para determinar a solução com precisão de 10 - 3 para a equação x 3 - x - 1 = 0 no intervalo 1,2 . Use x 0 = 1

Passo 1

Primeiro precisamos encontrar uma função de iteração g x = x convergente. Se liga:

x 3 - x - 1 = 0

x 3 = x + 1

x = x + 1 3

Como g x = x ,

g x = x + 1 3

Beleza, achamos uma função de iteração. Agora vamos testar se ela é convergente:

g x = x + 1 3

g ' x = 1 3 x + 1 2 3

g x é contínua em todo R e g ' x é contínua em todo R menos o - 1 .

g ' x = 1 3 x + 1 2 3 < 1 ⇒ 1 < 3 x + 1 2 3 ⇒ 1 3 3 2 < x + 1 ⇒ x > - 0.81

Como nossa raiz está no intervalo 1,2 , não será problema definir um intervalo I centralizado na raiz e que contenha x 0 = 1 obedecendo às condições acima. Logo, essa função converge para a nossa raiz.

Passo 2

x 0 = 1

ε 1 = ε 2 = 0 . 0 0 1

Iteração 1 :

x 1 = g x 0 = 1.2599

Critério de Parada:

x 1 - x 0 = 0.2 599 > 0.0 0 1

f x 1 = 0.2 599 > 0.0 0 1

Como não alcançamos o critério de parada continuamos as iterações.

x 1 = 1.2 599

ε 1 = ε 2 = 0 . 0 01

Iteração 2 :

x 2 = g x 1 = 1.312 3

Critério de Parada:

x 2 - x 1 = 0.052 4 > 0.0 0 1

f x 2 = 0.052 4 > 0.0 0 1

Como não alcançamos o critério de parada continuamos as iterações.

x 2 = 1.312 3

ε 1 = ε 2 = 0 . 0 0 1

Iteração 3 :

x 3 = g x 2 = 1.322 4

Critério de Parada:

x 3 - x 2 = 0.01 01 > 0.0 0 1

f x 3 = 0.01 01 > 0.0 0 1

Como não alcançamos o critério de parada continuamos as iterações.

x 3 = 1.3224

ε 1 = ε 2 = 0.001

Iteração 4 :

x 4 = g x 3 = 1.3243

Critério de Parada:

x 4 - x 3 = 0.0019 > 0.001

f x 4 = 0.0019 > 0.001

Como não alcançamos o critério de parada continuamos as iterações.

x 4 = 1.3243

ε 1 = ε 2 = 0.001

Iteração 5 :

x 5 = g x 4 = 1.3246

Critério de Parada:

x 5 - x 4 = 0.0004 < 0.001

f x 5 = 0.0004 < 0.001

Beleza, alcançamos o critério de parada! Então nossa aproximação é x - = 1.3246 .

Resposta

x - = 1.3246 .

Exercício Resolvido #2

Seja f x = x 2 - 4 sin ⁡ x . Determine um intervalo onde se encontre uma raiz e encontre uma função de iteração convergente para qualquer x 0 no intervalo. Em seguida, aplique o método do ponto fixo para encontrar uma aproximação da raiz com precisão de 10 - 4 .

Obs.: Desconsidere a raiz x = 0 .

Passo 1

Poxa, x = 0 tava molinho pra gente né? Mas não tem problema, vamos achar a outra raiz:

f 0 = 0

f 1 = 1 - 4 sin ⁡ 1 ≅ - 2.37

f 2 = 4 - 4 sin ⁡ 2 ≅ 0.36

Beleza, sabemos que existe uma raiz no intervalo 1,2 !

Passo 2

Agora vamos encontrar uma função de iteração g x = x convergente.

x 2 - 4 sin ⁡ x = 0

x 2 = 4 sin ⁡ x

x = 4 sin ⁡ x

x = 2 sin ⁡ x

Como g x = x ,

g x = 2 sin ⁡ x

Beleza, achamos uma função de iteração. Agora vamos testar se ela é convergente:

g x = 2 sin ⁡ x

g ' x = cos ⁡ x sin ⁡ x

g x é contínua para sin ⁡ x ≥ 0 , ou seja, x ∈ 2 n π , 2 n + 1 π , n ∈ R . O que quer dizer que ela é contínua em 0 , π , o que engloba nosso intervalo 1,2 .

g ' x é contínua para sin ⁡ x > 0 , ou seja, x ∈   ] 2 n π , ( 2 n + 1 ) π ] ,   n ∈ R . O que quer dizer que ela é contínua em ] 0 , π ] , o que engloba nosso intervalo 1,2 .

g ' x = cos ⁡ x sin ⁡ x < 1

A função g ' x é decrescente no intervalo 0 , π . Como g ' 1 ≅ 0.589 e g ' 2 ≅ - 0.436 , é seguro afirmar que em nenhum momento seu módulo é maior que 1 no intervalo 1,2 .

Logo, podemos ver que essa função de iteração converge.

Passo 3

Vamos tomar como estimativa inicial x 0 = 2

x 0 = 2

ε 1 = ε 2 = 0.0001

Iteração 1 :

x 1 = g x 0 = 1.90714

Critério de Parada:

x 1 - x 0 = 0.09285 > 0.0 0 01

f x 1 = 0.13868 > 0.00 0 1

Como não alcançamos o critério de parada continuamos as iterações.

x 1 = 1.90714

ε 1 = ε 2 = 0.0 0 01

Iteração 2 :

x 2 = g x 1 = 1.94316

Critério de Parada:

x 2 - x 1 = 0.03602 > 0.0 0 01

f x 2 = 0.04999 > 0.0 0 01

Como não alcançamos o critério de parada continuamos as iterações.

x 2 = 1.94316

ε 1 = ε 2 = 0.0 0 01

Iteração 3 :

x 3 = g x 2 = 1.93025

Critério de Parada:

x 3 - x 2 = 0.01291 > 0.0 0 01

f x 3 = 0.01847 > 0.0 0 01

Como não alcançamos o critério de parada continuamos as iterações.

x 3 = 1.93025

ε 1 = ε 2 = 0.0 0 01

Iteração 4 :

x 4 = g x 3 = 1.93503

Critério de Parada:

x 4 - x 3 = 0.00478 > 0.0 0 01

f x 4 = 0.00677 > 0.0 0 01

Como não alcançamos o critério de parada continuamos as iterações.

x 4 = 1.93503

ε 1 = ε 2 = 0.0 0 01

Iteração 5 :

x 5 = g x 4 = 1.93328

Critério de Parada:

x 5 - x 4 = 0.00175 > 0.0 0 01

f x 5 = 0.00249 > 0.0 0 01

Como não alcançamos o critério de parada continuamos as iterações.

x 5 = 1.93328

ε 1 = ε 2 = 0.0001

Iteração 6 :

x 6 = g x 5 = 1.93393

Critério de Parada:

x 6 - x 5 = 0.00064 > 0.0001

f x 6 = 0.00091 > 0.0001

Como não alcançamos o critério de parada continuamos as iterações.

x 6 = 1.93393

ε 1 = ε 2 = 0.0001

Iteração 7 :

x 7 = g x 6 = 1.93369

Critério de Parada:

x 7 - x 6 = 0.00024 > 0.0001

f x 7 = 0.00034 > 0.0001

Como não alcançamos o critério de parada continuamos as iterações.

x 7 = 1.93369

ε 1 = ε 2 = 0.0001

Iteração 8 :

x 8 = g x 7 = 1.93378

Critério de Parada:

x 8 - x 7 = 0.00008 < 0.0001

f x 8 = 0.00012 > 0.0001

Beleza, alcançamos o critério de parada! Então nossa aproximação é x - = 1.93378 .

Resposta

x - = 1.93378 .

Exercício Resolvido #3

Seja f x = 3 x 2 - x - 2 . Determine um intervalo onde se encontre uma raiz e encontre uma função de iteração convergente para qualquer x 0 no intervalo. Em seguida, aplique o método do ponto fixo para encontrar uma aproximação da raiz com precisão de 10 - 4 .

Passo 1

Vamos encontrar o intervalo onde essa raiz ta se escondendo:

Dessa vez temos que tomar cuidado já que a função não está definida no 0 .

f 0.5 = 9.5

f 1 = 0

Eita, achamos a raiz. Bom, vamos fingir que não achamos e continuar pra praticar um pouco haha. Beleza, sabemos que existe uma raiz no intervalo 0.5,1 !

Passo 2

Agora vamos encontrar uma função de iteração g x = x convergente.

3 x 2 - x - 2 = 0

3 x 2 = x + 2

3 = x 3 + 2 x 2

x 2 x + 2 = 3

x 2 = 3 x + 2

x = 3 x + 2

Como g x = x ,

g x = 3 x + 2

Beleza, achamos uma função de iteração. Agora vamos testar se ela é convergente:

g x = 3 x + 2

g ' x = - 3 2 1 x + 2 3 2

Ambas g x e g ' x são contínuas para x > - 2 , o que faz elas contínuas em 0.5,1

g ' x = 3 2 1 x + 2 3 2 < 1 ⇒ 1 x + 2 < 2 3 2 3 ⇒ 3 2 2 3 < x + 2

x > 3 2 2 3 - 2 ≅ - 1.091

Então no final das contas precisamos que x > - 1.091 , o que não causa nenhum problema pra gente criar o intervalo centrado em 1 (que nós sabemos ser a raiz).

Logo, podemos ver que essa função de iteração converge.

Passo 3

Vamos tomar como estimativa inicial x 0 = 0.5

x 0 = 0.5

ε 1 = ε 2 = 0.0001

Iteração 1 :

x 1 = g x 0 = 1.0954 5

Critério de Parada:

x 1 - x 0 = 0.59545 > 0.0001

f x 1 = 0.59545 > 0.0001

Como não alcançamos o critério de parada continuamos as iterações.

x 1 = 1.0954 5

ε 1 = ε 2 = 0.0001

Iteração 2 :

x 2 = g x 1 = 0.98446

Critério de Parada:

x 2 - x 1 = 0.11098 > 0.0001

f x 2 = 0.11098 > 0.0001

Como não alcançamos o critério de parada continuamos as iterações.

x 2 = 0.98446

ε 1 = ε 2 = 0.0001

Iteração 3 :

x 3 = g x 2 = 1.0026

Critério de Parada:

x 3 - x 2 = 0.01814 > 0.0001

f x 3 = 0.01814 > 0.0001

Como não alcançamos o critério de parada continuamos as iterações.

x 3 = 1.0026

ε 1 = ε 2 = 0.0001

Iteração 4 :

x 4 = g x 3 = 0.99957

Critério de Parada:

x 4 - x 3 = 0.00303 > 0.0001

f x 4 = 0.00303 > 0.0001

Como não alcançamos o critério de parada continuamos as iterações.

x 4 = 0.99957

ε 1 = ε 2 = 0.0001

Iteração 5 :

x 5 = g x 4 = 1.00007

Critério de Parada:

x 5 - x 4 = 0.00051 > 0.0001

f x 5 = 0.00051 > 0.0001

Como não alcançamos o critério de parada continuamos as iterações.

x 5 = 1.00007

ε 1 = ε 2 = 0.0001

Iteração 6 :

x 6 = g x 5 = 0.99999

Critério de Parada:

x 6 - x 5 = 0.00008 < 0.0001

f x 6 = 0.00008 < 0.0001

Pronto, alcançamos o critério de parada! Então nossa aproximação é x - = 0.99999 , o que é bem próximo do 1 né? haha

Resposta

x - = 0.99999 .

Exercício Resolvido #4

Encontre o zero de f x = e x - x 2 + 4 com precisão ε = 10 - 6 usando o Método do ponto Fixo. Escolha uma função de iteração apropriada e um ponto inicial x 0 .

Passo 1

Primeiro vamos encontrar uma aproximação inicial:

f - 1 ≅ 3.37

f 0 = 5

f 1 ≅ 5.71

Eita, e agora? Vamos tentar outros valores né, fazer o que:

f - 2 ≅ 0.135

Opa, estamos mais perto do 0 , vamos tentar descer mais:

f - 3 ≅ - 4.95

Beleza, nosso intervalo é - 3 , - 2 .

Passo 2

Agora vamos encontrar uma função de iteração g x = x convergente.

e x - x 2 + 4 = 0

x 2 = e x + 4

x = - e x + 4

Precisamos daquele pois estamos trabalhando em um intervalo negativo.

Como g x = x ,

g x = - e x + 4

Beleza, achamos uma função de iteração. Agora vamos testar se ela é convergente:

g x = - e x + 4

g ' x = - e x 2 e x + 4

Ambas g x e g ' x são contínuas em R .

g ' x = e x 2 e x + 4 < 1

Temos que g ' x é decrescente e g ' - 2 ≅ - 0.0333 , g ' - 3 ≅ - 0.01237 . Logo, é seguro falarmos que g ' x < 1 no intervalo dado.

Logo, podemos ver que essa função de iteração converge.

Passo 3

Vamos tomar como estimativa inicial x 0 = - 2

x 0 = - 2

ε 1 = ε 2 = 10 - 6

Iteração 1 :

x 1 = g x 0 = - 2.0335524

Critério de Parada:

x 1 - x 0 = 0.0335524 > 10 - 6

f x 1 = 0.0044655 > 10 - 6

Como não alcançamos o critério de parada continuamos as iterações.

x 1 = - 2.0335524

ε 1 = ε 2 = 10 - 6

Iteração 2 :

x 2 = g x 1 = - 2.0324541

Critério de Parada:

x 2 - x 1 = 0.0010982 > 10 - 6

f x 2 = 0.0001438 > 10 - 6

Como não alcançamos o critério de parada continuamos as iterações.

x 2 = - 2.0324541

ε 1 = ε 2 = 10 - 6

Iteração 3 :

x 3 = g x 2 = - 2.0324895

Critério de Parada:

x 3 - x 2 = 3.54 ∙ 10 - 5 > 10 - 6

f x 3 = 4.63 ∙ 10 - 6 > 10 - 6

Como não alcançamos o critério de parada continuamos as iterações.

x 3 = - 2.0324895

ε 1 = ε 2 = 10 - 6

Iteração 4 :

x 4 = g x 3 = - 2.0324884

Critério de Parada:

x 4 - x 3 = 1.1 ∙ 10 - 6 > 10 - 6

f x 4 = 1.5 ∙ 10 - 7 < 10 - 6

Pronto, alcançamos o critério de parada! Então nossa aproximação é x - = - 2.0324884 .

Resposta

x - = - 2.0324884 .

Exercício Resolvido #5

Encontre o zero de f x = x 3 - 2 x - 5 no intervalo 2,3 com precisão ε = 10 - 5 usando o Método do ponto Fixo. Escolha uma função de iteração apropriada e um ponto inicial x 0 .

Passo 1

Primeiro vamos encontrar uma função de iteração g x = x convergente.

x 3 - 2 x - 5   = 0

x 3 = 2 x + 5

x = 2 x + 5 3

Como g x = x ,

g x = 2 x + 5 3

Beleza, achamos uma função de iteração. Ao invés de usarmos aquele teorema vamos direto para o método e testamos por ele ;)

Passo 2

Vamos tomar como estimativa inicial x 0 = 2

x 0 = 2

ε 1 = ε 2 = 10 - 5

Iteração 1 :

x 1 = g x 0 = 2.080084

Critério de Parada:

x 1 - x 0 = 0.080084 > 10 - 5

f x 1 = 0.160168 > 10 - 5

Como não alcançamos o critério de parada continuamos as iterações.

x 1 = 2.080084

ε 1 = ε 2 = 10 - 5

Iteração 2 :

x 2 = g x 1 = 2.092351

Critério de Parada:

x 2 - x 1 = 0.012267 > 10 - 5

f x 2 = 0.024534 > 10 - 5

Como não alcançamos o critério de parada continuamos as iterações.

x 2 = 2.092351

ε 1 = ε 2 = 10 - 5

Iteração 3 :

x 3 = g x 2 = 2.094217

Critério de Parada:

x 3 - x 2 = 0.001866 > 10 - 5

f x 3 = 0.003733 > 10 - 5

Como não alcançamos o critério de parada continuamos as iterações.

x 3 = 2.094217

ε 1 = ε 2 = 10 - 5

Iteração 4 :

x 4 = g x 3 = 2.094501

Critério de Parada:

x 4 - x 3 = 0.000284 > 10 - 5

f x 4 = 0.000567 > 10 - 5

Como não alcançamos o critério de parada continuamos as iterações.

x 4 = 2.094501

ε 1 = ε 2 = 10 - 5

Iteração 5 :

x 5 = g x 4 = 2.094544

Critério de Parada:

x 5 - x 4 = 4.3 ∙ 10 - 5 > 10 - 5

f x 5 = 8.6 ∙ 10 - 5 > 10 - 5

Como não alcançamos o critério de parada continuamos as iterações.

x 5 = 2.094544

ε 1 = ε 2 = 10 - 5

Iteração 6 :

x 6 = g x 5 = 2.09455

Critério de Parada:

x 6 - x 5 = 6 ∙ 10 - 6 < 10 - 5

f x 6 = 1.3 ∙ 10 - 5 > 10 - 5

Pronto, alcançamos o critério de parada! Então nossa aproximação é x - = 2.09455 .

Resposta

x - = 2.09455 .

Exercício Resolvido #6

Encontre o zero de f x = e - x 2 - cos ⁡ x no intervalo 1,2 com precisão ε = 10 - 4 usando o Método do Ponto Fixo. Escolha uma função de iteração apropriada e um ponto inicial x 0 .

Passo 1

Primeiro vamos encontrar uma função de iteração g x = x convergente.

e - x 2 - cos ⁡ x = 0

Aqui ta meio ruim de isolar o x né? Nesses casos podemos usar aquela fórmula geral que vimos na teoria:

g x = x + A x f x

Logo, uma possível função de iteração seria

g x = x + cos ⁡ x - e - x 2

Beleza, achamos uma função de iteração. Ao invés de usarmos aquele teorema vamos direto para o método e testamos por ele ;)

Passo 2

Vamos tomar como estimativa inicial x 0 = 1.5

x 0 = 1.5

ε 1 = ε 2 = 10 - 4

Iteração 1 :

x 1 = g x 0 = 1.46534

Critério de Parada:

x 1 - x 0 = 0.03466 > 10 - 4

f x 1 = 0.01155 > 10 - 4

Como não alcançamos o critério de parada continuamos as iterações.

x 1 = 1.46534

ε 1 = ε 2 = 10 - 4

Iteração 2 :

x 2 = g x 1 = 1.45379

Critério de Parada:

x 2 - x 1 = 0.01155 > 10 - 4

f x 2 = 0.00408 > 10 - 4

Como não alcançamos o critério de parada continuamos as iterações.

x 2 = 1.45379

ε 1 = ε 2 = 10 - 4

Iteração 3 :

x 3 = g x 2 = 1.44972

Critério de Parada:

x 3 - x 2 = 0.00408 > 10 - 4

f x 3 = 0.00147 > 10 - 4

Como não alcançamos o critério de parada continuamos as iterações.

x 3 = 1.44972

ε 1 = ε 2 = 10 - 4

Iteração 4 :

x 4 = g x 3 = 1.44825

Critério de Parada:

x 4 - x 3 = 0.00147 > 10 - 4

f x 4 = 0.00053 > 10 - 4

Como não alcançamos o critério de parada continuamos as iterações.

x 4 = 1.44825

ε 1 = ε 2 = 10 - 4

Iteração 5 :

x 5 = g x 4 = 1.44772

Critério de Parada:

x 5 - x 4 = 0.00053 > 10 - 4

f x 5 = 0.00019 > 10 - 4

Como não alcançamos o critério de parada continuamos as iterações.

x 5 = 1.44772

ε 1 = ε 2 = 10 - 4

Iteração 6 :

x 6 = g x 5 = 1.44752

Critério de Parada:

x 6 - x 5 = 0.00019 > 10 - 4

f x 6 = 7 ∙ 10 - 5 < 10 - 4

Beleza, alcançamos o critério de parada! Então nossa aproximação é x - = 1.44752 .

Resposta

x - = 1.44752 .

Exercício Resolvido #7

Suponha que a > 0 . Certifique-se que as seguintes funções de iteração:

a φ 1 x = a x

b φ 2 x = 1 2 x + a x

possuem o ponto fixo ξ = a . Qual dessas funções você utilizaria para obter uma aproximação de a pelo Método do Ponto Fixo? Explique sua escolha.

Passo 1

Se a é ponto fixo dessas funções então   φ 1 a =   φ 2 a = a :

φ 1 a = a a = a ∙ a a ∙ a = a

φ 2 a = 1 2 a + a a = 1 2 a + a = a

Beleza, ambas as funções de iteração tem a como ponto fixo.

Passo 2

-

Agora o problema é, qual escolheríamos? O lógico seria escolhermos a que converge para o ponto fixo mais rapidamente certo? E lembra que quanto menor o módulo da derivada da função de iteração mais rápido o método vai convergir? Então queremos a função que tenha a menor derivada no ponto fixo.

φ 1 ' x = - a x 2

φ 1 ' a = a a = 1

Agora para a segunda função:

φ 2 ' x = 1 2 1 - a x 2 = 1 2 - a 2 x 2

φ 2 ' a = 1 2 - a 2 a = 1 2 - 1 2 = 0

A segunda função de iteração tem módulo da derivada igual a zero no ponto fixo! Então ela é a melhor escolha para função de iteração já que converge mais rapidamente.

Resposta

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

Exercício Resolvido #8

Seja f x = 3 x 2 - e x . Determine uma função de iteração φ x e um intervalo a , b no qual as iterações do Método do Ponto Fixo convergirão para uma raiz positiva de f . Encontre as soluções com precisão de 10 - 5 .

Passo 1

Vamos usar como função de iteração φ x = x - x 2 + e x 3

Passo 2

Agora precisamos encontrar um intervalo que contenha alguma raiz:

f 0 = - 1

f 1 ≅ 0.281

Então um intervalo é 0,1 . Agora vamos na fé que vai dar certo haha.

Passo 3

Tomando x 0 = 1 :

x 0 = 1

ε 1 = ε 2 = 10 - 5

Iteração 1 :

x 1 = g x 0 = 0.906094

Critério de Parada:

x 1 - x 0 = 0.093906 > 10 - 5

f x 1 = 0.011619 > 10 - 5

Como não alcançamos o critério de parada continuamos as iterações.

x 1 = 0.906094

ε 1 = ε 2 = 10 - 5

Iteração 2 :

x 2 = g x 1 = 0.909967

Critério de Parada:

x 2 - x 1 = 0.003873 > 10 - 5

f x 2 = 0.000121 > 10 - 5

Como não alcançamos o critério de parada continuamos as iterações.

x 2 = 0.909967

ε 1 = ε 2 = 10 - 5

Iteração 3 :

x 3 = g x 2 = 0.910007

Critério de Parada:

x 3 - x 2 = 4 ∙ 10 - 5 > 10 - 5

f x 3 = 9.8 ∙ 10 - 7 < 10 - 5

Beleza, chegamos no nosso critério de parada! Tranquilo né? Esse é o poder que uma boa escolha de função de iteração tem, mas pra saber qual usar e quando usar só praticando mesmo =/

Então nossa raiz é x - = 0.910007 .

Resposta

x - = 0.910007

Exercício Resolvido #9

Seja f x = x - cos ⁡ x . Determine uma função de iteração φ x e um intervalo a , b no qual as iterações do Método do Ponto Fixo convergirão para uma raiz positiva de f . Encontre as soluções com precisão de 10 - 5 .

Passo 1

Vamos usar como função de iteração φ x = x - x 2 + cos ⁡ x 2 = x 2 + cos ⁡ x 2 .

Passo 2

Agora precisamos encontrar um intervalo que contenha alguma raiz:

f 0 = - 1

f 1 ≅ 0 . 46

Então um intervalo é 0,1 .

Temos que φ ' x = 1 2 - sin ⁡ x 2 < 1 para todos os reais.

Logo, nossa função de iteração vai convergir no intervalo 0,1 .

Passo 3

Tomando x 0 = 1 :

x 0 = 1

ε 1 = ε 2 = 10 - 5

Iteração 1 :

x 1 = g x 0 = 0.770151

Critério de Parada:

x 1 - x 0 = 0.229849 > 10 - 5

f x 1 = 0.052346 > 10 - 5

Como não alcançamos o critério de parada continuamos as iterações.

x 1 = 0.770151

ε 1 = ε 2 = 10 - 5

Iteração 2 :

x 2 = g x 1 = 0.743978

Critério de Parada:

x 2 - x 1 = 0.026173 > 10 - 5

f x 2 = 0.008198 > 10 - 5

Como não alcançamos o critério de parada continuamos as iterações.

x 2 = 0.743978

ε 1 = ε 2 = 10 - 5

Iteração 3 :

x 3 = g x 2 = 0.739879

Critério de Parada:

x 3 - x 2 = 0.004099 > 10 - 5

f x 3 = 0.001329 > 10 - 5

Como não alcançamos o critério de parada continuamos as iterações.

x 3 = 0.739879

ε 1 = ε 2 = 10 - 5

Iteração 4 :

x 4 = g x 3 = 0.73921 5

Critério de Parada:

x 4 - x 3 = 0.000665 > 10 - 5

f x 4 = 0.000217 > 10 - 5

Como não alcançamos o critério de parada continuamos as iterações.

x 4 = 0.73921 5

ε 1 = ε 2 = 10 - 5

Iteração 5 :

x 5 = g x 4 = 0.739106

Critério de Parada:

x 5 - x 4 = 0.000108 > 10 - 5

f x 5 = 3.5 ∙ 10 - 5 > 10 - 5

Como não alcançamos o critério de parada continuamos as iterações.

x 5 = 0.739106

ε 1 = ε 2 = 10 - 5

Iteração 6 :

x 6 = g x 5 = 0.739089

Critério de Parada:

x 6 - x 5 = 1.8 ∙ 10 - 5 > 10 - 5

f x 6 = 5.7 ∙ 10 - 6 < 10 - 5

Pronto, alcançamos o critério de parada! Então nossa aproximação é x - = 0.739089 .

Resposta

x - = 0.739089

Exercício Resolvido #10

Seja f x = x 3 - 3 x + 1 :

a Mostre que f x possui uma raiz em 0,1 .

b Mostre que φ x = x 3 + 1 3 é uma possível função de iteração obtida a partir de f x .

c Verifique se φ x satisfaz as condições i e i i do Teorema de Convergência.

d Encontre uma estimativa para a raiz de f x através do método do ponto fixo e usando a função φ x do item c tal que f x < 0.007 .

Passo 1

a Para mostrar que f x possui uma raiz nesse intervalo é só compararmos o valor dela nos extremos:

f 0 = 1

f 1 = - 1

f 0 ∙ f 1 < 0

Logo, existe uma raiz em 0,1 .

Passo 2

b Fazendo algumas manipulações algébricas temos:

x 3 - 3 x + 1 = 0

x 3 + 1 = 3 x

x = x 3 + 1 3

φ x = x

φ x = x 3 + 1 3

Pronto, φ x é uma possível função de iteração obtida a partir de f x .

Passo 3

c Vamos começar pela condição i :

φ x = x 3 + 1 3

φ ' x = x 2

Ambas são contínuas em R , logo a condição i é satisfeita.

Condição i i :

φ ' x = x 2 < 1

- 1 < x < 1

Sendo assim, no intervalo 0,1 a função só não satisfaz as condições para o extremo x = 1 .

Passo 4

d Tomando x 0 = 0 :

ε 1 = 0.007

Iteração 1 :

x 1 = g x 0 = 0.3333

Critério de Parada:

f x 1 = 0.037 > 0.007

Como não alcançamos o critério de parada continuamos as iterações.

x 1 = 0.3333

ε 1 = ε 2 = 0.007

Iteração 2 :

x 2 = g x 1 = 0.3457

Critério de Parada:

f x 2 = 0.0042 < 0.007

Fim de papo, alcançamos o critério de parada! Então nossa aproximação é x - = 0.3457 .

Resposta

x - = 0.3457

Exercícios de Livros Relacionados

A equação x 2 - b = 0 tem como raiz ξ = b . Considere o MPF com φ ( x ) = b x : Comprove que φ ' ( ξ ) =- 1 ; o que acontece com a sequência { x k } tal que x k + 1 = φ ( x k ) ? sua conclusão do item
Considere a função f(x) = x3 - x - 1 (do Exemplo 19). Resolva-a pelo MPF com função de iteração φ ( x ) = 1 x + 1 x 2 e x 0 = 1 . Justifique seus resultados.
Ao se aplicar o método do ponto fixo (MPF) à resolução de uma equação, obtivemos os seguintes resultados nas iterações indicadas:Escreva o que puder a respeito da raiz procurada.
9. Se a função de iteração do MPF for tal que as condições do Teorema 2 estão satisfeitas:a) mostre que ξ - x kM 1 - M x k - x k - 1 ;b) para que valores de M teremos então que ξ - x kε se x k - x k -