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 ๐ = ๐. ๐๐