Zero de função Problema O cálculo de raízes de funções encontra um grande emprego na obtenção da solução de uma vasta gama de problemas de engenharia. Em geral, trata-se de determinar o(s) valores de x tal que f(x)=0, onde f é a função cujo raízes são a determinar. Métodos matemáticos A matemática fornece métodos formais que permite a determinação exata das raízes em diversos casos. Os métodos mais conhecido permitem a determinação de raízes de polinômios ate grau 3, ou grau maior mais em certas condições. Em muitas situações, a resolução matemática necessita de intuição para que elas sejam transformadas em casos resolvíveis. Exemplos Polinômios do primeiro e segundo grau ou transformáveis em polinômios do primeiro ou segundo grau: 2 x 3 0; x 2 3x 5 0 2sin x 3 0;sin 2 x 3sin x 5 0 Funções cuja a recíproca é conhecida: log10 x 5 0 Determinação gráfica A representação gráfica de uma função é uma fonte de informações úteis sobre o comportamento da função, particularmente para a determinação das raízes. Além disso, o grafo permite de compreender o funcionamento dos métodos numéricos para determinar as raízes. Raízes com gráfico Raízes são dadas pelos pontos de interseção do grafo com o eixo dos x. Métodos numéricos Mesmo com um método formal, o(s) valor(es) calculado(s) pelo computador é aproximado, a não seja usar um CAS. Existem métodos numéricos que permite aproximar as raízes em casos gerais, inclusivo casos que a matemática não resolva de formalmente. Métodos numéricos Vamos estudar três métodos de determinação de raízes: Bisseção Secante Newton-Raphson Bisseção Th: Se y=f(x) é uma função contínua e muda de sinal no intervalo [a,b] (isto é se f(a).f(b)<0), então existe pelo menos um ponto x0 [a,b] tal que f(x0)=0. Além disso, se f’(x) não muda de sinal em [a,b], x0 é a única raiz de f(x) nesse intervalo. Bisseção Para se aproximar de uma raiz, o princípio da bisseção consista em reduzir o intervalo inicial testando o sinal de f(x) para o ponto médio do intervalo. Considerando o intervalo [a,b] ab Se f (a). f ( ) 0 , o novo intervalo e [a,(a+b)/2] 2 Se f (b). f ( ab ) 0, o novo intervalo e [(a+b)/2,b] 2 Algoritmo Raiz(f,a,b,tol) Enquanto (|a-b|>tol) x=(a+b)/2 Se f(x).f(a)<0 b=x Senão a=x Resultado=(a+b)/2 Bisseção Esse método, com um bom escolhe do intervalo inicial, é adaptado com a representação dos números do computador: a divisão por 2 a cada passo é uma operação simples. A convergência do algoritmo é garantida, o algoritmo não saia do intervalo inicial, esse intervalo é cada vez dividido por dois, A convergência é muito lenta: para ganhar uma decimal (base 10), preciso de 3 a 4 passos. Secante O método da secante funciona sobre o mesmo princípio que a bisseção e necessita da mesma condição inicial: continuidade da função. Secante Com esse método, determinamos um ponto a partir da assimilação da curva com um segmento passando pelos pontos (XE, f(XE)) e (XD, f(YD)). O candidato para ser raiz é o ponto de interseção desse segmento com o eixo x. Secante Determinação de XN: Temos a relação: f (XD) XD XN f (XE ) XE XN De onde podemos extrair XN: f (XD)XE f (XE )XD XN f (XD) f (XE ) Secante O segmento (XN,f(XN)); (XD,f(XD)) é usado para determinar o valor do passo seguinte. Algoritmo Raiz(f,a,b,iter) Repete iter vezes b=(b.f(a)-a.f(b))/(f(a)-f(b)) Resultado=b Falsa posição O método da falsa posição aparece como uma combinação entre o método da secante e a bisseção. As condições iniciais são as mesma que no caso da bisseção (intervalo onde a função troca de sinal). Falsa posição Como no caso da secante, determinamos um ponto a partir da assimilação da curva com um segmento passando pelos pontos (XE, f(XE)) e (XD, f(YD)). Temos: f (XD)XE f (XE )XD XN f (XD) f (XE ) Falsa posição No caso da falsa posição, o novo segmento é determinado em função dos sinais de f(XN)f(XD) e f(XN)f(XE). Se f troca de sinal entre XE e XN, o novo intervalo é [XE, XN], senão o novo intervalo é [XN, XE]. Algoritmo Raiz(f,a,b,iter) Repete iter vezes x=(b.f(a)-a.f(b))/(f(a)-f(b)) Se f(x).f(a)<0, b=x Senão a=x Resultado=x Newton-Raphson O método de NewtonRaphson não precisa de um intervalo inicial. Ela considera que a curva no ponto inicial pode ser aproximada com a reta tangente à curva nesse ponto. Newton-Raphson De forma equivalente, consista também a considerar a função como aproximada nesse ponto pela série de Taylor de 1° grau: f(x1)=f(x0)+(x1-x0).f’(x0) Determinação de XN: XN=XD-f(XD)/f’(XD) Newton-Raphson Por um processo iterativo, a raiz pode ser aproximada: xi+1=xi-f(xi)/f’(xi) Algoritmo Raiz(f,x0,iter) X=x0 Repete iter vezes X=X-f(X)/f’(X) Resultado=X Newton-Raphson e Secante Os dois métodos de secante e NewtonRaphson são próximos. O método da secante é o método de Newton-Raphson aonde a derivada no ponto inicial é substituída pela diferencia finita. A vantagem da secante é que não é necessário conhecer a função derivada. Convergência A convergência desses métodos é em geral mais rápida que no caso da bisseção. O método da bisseção usa sempre o mesmo algoritmo para qualquer função enquanto os outros métodos usam o comportamento da curva (diferencia finita ou derivada) para se aproximar da raiz. Convergência Se Newton-Raphson e Secante podem ser mais eficiente, elas podem ser também com dificuldade de convergência se a função tem variação do sinal da derivada próxima da raiz procurada. Convergência Vários critérios podem ser usados para decidir de para a aplicação do algoritmo: um número dado de iterações, quando a diferencia entre dois passo de uma iteração é menos que um erro |xi+1-xi|< e, quando o valor da função em xi é perto de 0 |f(xi)|<e, quando os dois últimos critérios não para o algoritmo, ele pode ser parado porque considerado como não convergente.