Zeros de funções
Prof. Rafael mesquita
[email protected]
Introdução
๏ฎ
๏ฎ
๏ฎ
Nas mais diversas áreas das ciências exatas
ocorrem situações que envolvem a solução de
uma equação do tipo f(x) = 0
Um número real ๐œบ é um zero da função f(x) ou uma
raíz da equação f(x)=0 se f(๐œบ) = 0
Os valores de x que anulam f(x) podem ser reais
ou complexos
โ€“ Nosso estudo interessa-se apenas pelos valores reais
Introdução
๏ฎ
Podemos resolver f(x) = 0 por dois caminhos
distintos
โ€“ Métodos diretos
๏‚ง Métodos analíticos
๏‚ง Numero finito de operações
๏‚ง Processos particulares
โ€“ Cada tipo de função deve possuir seu próprio caminho
para a solução
โ€“ Métodos iterativos
๏‚ง Partem de uma aproximação inicial da solução
๏‚ง A cada iteração, uma nova aproximação é gerada
โ€“ Até que uma solução satisfatória seja encontrada
Introdução
๏ฎ
Para algumas equações, como as polinômiais de
grau 2 existem fórmulas diretas para se encontrar
as raízes
โ€“ Métodos diretos!
๏ฎ
No entanto em polinômios de grau mais alto e no
caso de funções mais complicadas essa tarefa
não é trivial
โ€“ Métodos iterativos!
โ€“ Uso de aproximações
โ€“ Precisão prefixada
Introdução
๏ฎ
1.
Idéia central de métodos numéricos baseados em
aproximações para encontrar zeros de funções:
Localização ou isolamento das raízes
โ€ข
2.
Obtenção de um intervalo contendo a raíz
Refinamento
โ€ข
Melhorar a aproximação inicial do intervalo até que
seja obtida a precisão prefixada
Aproximação Inicial
๏ฎ
Graficamente, os zeros reais são representados
pelas abcissas dos pontos onde uma curva
intercepta o eixo x
Aproximação Inicial
๏ฎ
๏ฎ
๏ฎ
Análise teórica e gráfica da função f(x)
Sucesso da fase de refinamento depende do
resultado dessa fase
Teorema de Bolzano:
โ€“ Seja f(x) uma função continua no intervalo [a,b]
โ€“ Se f(a).f(b) < 0, então existe pelo menos um ponto x = ๐œบ
entre a e b que é zero de f(x)
Aproximação Inicial
๏ฎ
Graficamente
โ€“ Se f(a).f(b) < 0, então existe pelo menos um ponto x = ๐œบ
entre a e b que é zero de f(x)
Aproximação Inicial
๏ฎ
Graficamente
โ€“ Se f(a).f(b) < 0, então existe pelo menos um ponto x = ๐œบ
entre a e b que é zero de f(x)
Não significa que
exista exatamente
uma raíz !
Aproximação Inicial
๏ฎ
Graficamente
โ€“ Se fโ€™(x) existir e preservar o sinal em (a,b), então esse
intervalo contém um único zero de f(x)
Aproximação Inicial
๏ฎ
Conclusão:
โ€“ Se f(a).f(b) < 0, então existe pelo menos um ponto x = ๐œบ
entre a e b que é zero de f(x)
โ€“ Se fโ€™(x) existir e preservar o sinal em (a,b), então esse
intervalo contém um único zero de f(x)
๏‚ง Intervalo de separação!
โ€“ Além disso, f deve ser contínua no intervalo [a;b]
Aproximação Inicial
๏ฎ
Conclusão:
โ€“ Uma forma de isolar as raízes de f(x) é tabelar f(x) para
diversos valores de x e analisar as mudanças de sinal
de f(x) e o sinal da derivada nos intervalos em que f(x)
mudou de sinal
Aproximação Inicial
๏ฎ
Exemplo 1
Aproximação Inicial
๏ฎ
Exemplo 1
Intervalos contém pelo
menos um zero de f(x)
Aproximação Inicial
๏ฎ
Exemplo 1
Intervalos contém pelo
menos um zero de f(x)
Como o polinômio é de
grau 3, sabemos que
cada intervalo contem
exatamente uma raíz!
Aproximação Inicial
๏ฎ
Exemplo 2
Intervalo contém pelo
menos um zero de f(x)
Aproximação Inicial
๏ฎ
Para saber se existe um único zero no intervalo,
analisamos o sinal de fโ€™(x)
๏ฎ
f(x) é continua dentro do intervalo [1;2]
fโ€™(x) não muda de sinal dentro do intervalo [1;2]
Como f(1).f(2)<0, concluímos que existe
exatamente uma raíz no intervalo [1;2]
๏ฎ
๏ฎ
Aproximação Inicial
๏ฎ
1.
2.
Na análise gráfica da função f(x) podemos utilizar
um dos seguintes processos:
Esboçar o gráfico da função f(x) e localizar as
abcissas dos pontos onde a curva intercepta o
eixo x
A partir da equação f(x) = 0, obter a equação
equivalente g(x) = h(x), esboçar os gráficos de
g(x) e h(x) e localizar os pontos x onde as duas
curvas se interceptam
1. f(๐œบ) = 0 <=> g(๐œบ) = h(๐œบ)
Aproximação Inicial
๏ฎ
Exemplo: ๐’‡ ๐’™ = ๐’™๐Ÿ‘ โˆ’ ๐Ÿ—๐’™ + ๐Ÿ‘
Aproximação Inicial
โ€ข
Exemplo: ๐’‡ ๐’™ = ๐’™๐Ÿ‘ โˆ’ ๐Ÿ—๐’™ + ๐Ÿ‘
Esboço do gráfico
๐’‡ ๐’™ = ๐’™๐Ÿ‘ โˆ’ ๐Ÿ—๐’™ + ๐Ÿ‘
โ€ข
๐’‡โ€ฒ ๐’™ = ๐Ÿ‘๐’™๐Ÿ โˆ’ ๐Ÿ—
โ€ข
๐’‡โ€ฒ ๐’™ = ๐ŸŽ โ‡’ ๐’™ = ± ๐Ÿ‘
โ€ข
๐’‡โ€ฒโ€ฒ ๐’™ = ๐Ÿ”๐’™ โ‡’ ๐’‡โ€ฒโ€ฒ โˆ’ ๐Ÿ‘ < ๐ŸŽ โ‡’ ๐’‘๐’•๐’ ๐’Ž๐’‚๐’™
โ€ข
โ‡’ ๐’‡โ€ฒโ€ฒ + ๐Ÿ‘ > ๐ŸŽ โ‡’ ๐’‘๐’•๐’ ๐’Ž๐’Š๐’
๏ฎ
1.
Aproximação Inicial
๏ฎ
1.
Exemplo: ๐’‡ ๐’™ = ๐’™๐Ÿ‘ โˆ’ ๐Ÿ—๐’™ + ๐Ÿ‘
Esboço do gráfico
Aproximação Inicial
๏ฎ
๏ฎ
๏ฎ
๏ฎ
Exemplo:๐’™๐Ÿ‘ โˆ’ ๐Ÿ—๐’™ + ๐Ÿ‘ = ๐ŸŽ
๐’‡ ๐’™ โ‡’๐’ˆ ๐’™ โˆ’๐’‰ ๐’™
๐’™๐Ÿ‘ โˆ’ ๐Ÿ—๐’™ + ๐Ÿ‘ = ๐ŸŽ
โ€“ โ‡’ ๐’™๐Ÿ‘ = ๐Ÿ—๐’™ โˆ’ ๐Ÿ‘
๐’ˆ ๐’™ = ๐’™๐Ÿ‘ e ๐ก ๐ฑ = ๐Ÿ—๐’™ โˆ’ ๐Ÿ‘
Aproximação Inicial
๏ฎ
โ€ข
Exemplo: ๐’‡ ๐’™ = ๐’™ โˆ’ ๐Ÿ“๐’†โˆ’๐’™
๐’ˆ ๐’™ = ๐’™ e ๐ก ๐ฑ = ๐Ÿ“๐’†โˆ’๐’™
g(x) < h(x)
Aproximação Inicial
๏ฎ
โ€ข
Exemplo: ๐’‡ ๐’™ = ๐’™ โˆ’ ๐Ÿ“๐’†โˆ’๐’™
๐’ˆ ๐’™ = ๐’™ e ๐ก ๐ฑ = ๐Ÿ“๐’†โˆ’๐’™
g(x) < h(x)
Aproximação Inicial
๏ฎ
โ€ข
Exemplo: ๐’‡ ๐’™ = ๐’™ โˆ’ ๐Ÿ“๐’†โˆ’๐’™
๐’ˆ ๐’™ = ๐’™ e ๐ก ๐ฑ = ๐Ÿ“๐’†โˆ’๐’™
g(x) < h(x)
Aproximação Inicial
๏ฎ
โ€ข
Exemplo: ๐’‡ ๐’™ = ๐’™ โˆ’ ๐Ÿ“๐’†โˆ’๐’™
๐’ˆ ๐’™ = ๐’™ e ๐ก ๐ฑ = ๐Ÿ“๐’†โˆ’๐’™
g(x) > h(x)
Curvas se
interceptam no
intervalo [1,2] !
Métodos iterativos
Como definir uma aproximação suficientemente
boa?
๏ฎ
Critérios de parada
1. ๐‘ฅ๐‘˜ está suficientemente próximo da raíz exata?
1. ๐‘ฅ โ€ฒ โˆ’ ๐œ€ < ๐ธ ou
2. Número máximo de iterações
๏ฎ
Como efetuar o teste 1.1 se não conhecemos ๐œบ ?
โ€“ Testamos se encontramos um intervalo [a,b] tal que:
โ€“
๐œ€ ๐œ– ๐‘Ž, ๐‘ ๐‘’
๐‘โˆ’๐‘Ž <๐ธ
โ€“ Então, โˆ€๐‘ฅ โˆˆ ๐‘Ž, ๐‘ , ๐‘ฅ โˆ’ ๐œ€ < ๐ธ
Como definir uma aproximação suficientemente
boa?
MÉTODO DA BISSEÇÃO
Método da bissecção
๏ฎ
๏ฎ
Seja a função f(x) contínua no intervalo [a,b] e tal
que f(a).f(b) < 0, supondo que este intervalo
contenha apenas uma única raíz
O objetivo deste método é reduzir a amplitude do
intervalo que contém a raíz até se atingir a
precisão requerida: ๐‘ โˆ’ ๐‘Ž < ๐ธ
โ€“ Utiliza a sucessiva divisão de [a,b] ao meio
Método da bissecção
1.
2.
3.
4.
Escolha a,b (extremos do intervalo de separação) e l (amplitude
final de [a;b])
Faça
1. ๐‘ = ๐‘ โˆ’ ๐‘Ž
2. ๐‘ฅ0 = (๐‘Ž + ๐‘)/2
Fim;
Enquanto ๐‘ > ๐‘™ ou ๐‘“(๐‘ฅ0 ) โ‰  0
1. Se ๐‘“ ๐‘Ž × ๐‘“ ๐‘ฅ0 < 0
1.
๐‘ = ๐‘ฅ0
2. Se ๐‘“ ๐‘ × ๐‘“ ๐‘ฅ0 < 0
1.
๐‘Ž = ๐‘ฅ0
3. ๐‘ = ๐‘ โˆ’ ๐‘Ž
4. ๐‘ฅ0 = (๐‘Ž + ๐‘)/2
5.
6.
Fim.
(๐‘ฅ0 contem o valor aproximado da raíz procurada)
Método da bissecção
Exercícios
๏ฎ
Calcule a raiz real da equação ๐’™๐Ÿ + ๐’๐’ ๐’™ = ๐ŸŽ com
tolerância máxima de ๐‘ฌ < ๐Ÿ๐ŸŽโˆ’๐Ÿ usando o método
da bisseção
โ€“ *Necessário calcular aproximação inicial
๏ฎ
Calcule a raíz real da equação xlog(x)-1 que possui
zero em [2,3] para um erro menor que 0,001
Precisão preestabelecida
๏ฎ
๏ฎ
Podemos prever exatamente qual o menor numero
de iterações para que a precisão estabelecida seja
alcançada
Dada uma função ๐’‡, um intervalo de separação ๐‘ฐ =
[๐’‚; ๐’ƒ], e a amplitude final ๐’:
โ€“ ๐‘๐‘˜ โˆ’ ๐‘Ž๐‘˜ =
โ€“
๐‘0 โˆ’๐‘Ž0
2๐‘˜
โ€“๐‘˜โ‰ฅ
๏ฎ
๐‘0 โˆ’๐‘Ž0
, โˆ€๐‘˜
2๐‘˜
โˆˆ ๐‘, logo
โ‰ค ๐‘™ โ‡’ 2๐‘˜ โ‰ฅ (๐‘0 โˆ’ ๐‘Ž0 )/๐‘™, assim
(ln ๐‘0 โˆ’๐‘Ž0 โˆ’ln(๐‘™))
ln 2
=๐‘ก
K é o menor inteiro maior ou igual a t
MÉTODO DA FALSA POSIÇÃO
Falsa posição
๏ฎ
Assim como o método da bisseção, também é um
método de quebra
โ€“ Quebra é realizada no ponto de interseção da reta
definida pelos pontos ๐‘Ž, ๐‘“ ๐‘Ž e (๐‘, ๐‘“(๐‘)) com o eixo x
โ€“ Substituímos a função f por uma reta
Falsa posição
๏ฎ
Idéia geométrica
Falsa posição
๏ฎ
Idéia geométrica
โ€“ ๐‘ฆ=๐‘“ ๐‘Ž +
๐‘“ ๐‘ โˆ’๐‘“(๐‘Ž)
๐‘โˆ’๐‘Ž
× (๐‘ฅ โˆ’ ๐‘Ž)
โ€“ Podemos definir a interseção com o eixo x (fazendo
y=0)
โ€“๐‘ฅ=
๐‘Ž๐‘“ ๐‘ โˆ’๐‘๐‘“(๐‘Ž)
๐‘“ ๐‘ โˆ’๐‘“(๐‘Ž)
Falsa posição
๏ฎ
Podemos obter um algoritmo para o método da
falsa posição similar ao algoritmo para o método
da bisseção
โ€“ Mudamos apenas a โ€œmáquina geradoraโ€
โ€“ Ao invés de calcular o valor de x como
โ€“ ๐‘ฅ = (๐‘Ž + ๐‘)/2, teremos
โ€“ ๐‘ฅ=
๐‘Ž๐‘“ ๐‘ โˆ’๐‘๐‘“(๐‘Ž)
๐‘“ ๐‘ โˆ’๐‘“(๐‘Ž)
Exercícios
1.
2.
3.
Aplique o método da falsa posição na função
๐’‡ ๐’™ = ๐’™๐Ÿ‘ โˆ’ ๐Ÿ—๐’™ + ๐Ÿ‘ no intervalo [0,1] considerando
๐ฅ = ๐Ÿ“. ๐Ÿ๐ŸŽโˆ’๐Ÿ’
Aplique o método da falsa posição na função
๐’‡ ๐’™ = ๐’„๐’๐’”(๐’™) + ๐’™ considerando ๐’ = ๐ŸŽ. ๐ŸŽ๐ŸŽ๐Ÿ
Aplique o método da falsa posição na função
๐’‡ ๐’™ = ๐’†๐’™ + ๐’™ considerando ๐’ = ๐ŸŽ. ๐ŸŽ๐Ÿ
Download

Aula 4 - Aproximaรงรฃo Inicial, Bolzano, Bissecao, falsas cordas