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