Computao Aplicada Geologia
DEPARTAMENTO DE GEOLOGIA
Matlab Equaes no lineares
FACULDADE DE CINCIAS
UNIVERSIDADE DE LISBOA
Problema
Resolver f(x)=0 quando no se tem uma frmula explcita para a
soluo.
oloss
l
pp
m
em
EExxe
sin( x) + e x = 0
x cos( x) = 0
x5 4 x + 3 = 0
Computao Aplicada Geologia
DEPARTAMENTO DE GEOLOGIA
FACULDADE DE CINCIAS
UNIVERSIDADE DE LISBOA
Matlab Equaes no lineares
Definio
Se f(z) = 0, ento diz-se que z uma raiz da equao f(x) = 0 ou
que z um zero da funo f.
Teorema do Valor Intermdio
Se f uma funo contnua no intervalo [a, b] e d um
nmero entre f(a) e f(b), ento existe um ponto x X[a, b]
tal que f(x) = d.
Computao Aplicada Geologia
DEPARTAMENTO DE GEOLOGIA
FACULDADE DE CINCIAS
UNIVERSIDADE DE LISBOA
Matlab Equaes no lineares
Bisseco
Considere-se o intervalo [a,b], em que f(a).f(b)<0 ou seja
onde existe pelo menos um zero da funo f.
O mtodo da bisseco consiste em subdividir,
sucessivamente, este intervalo a meio, considerando o
subintervalo onde se verifica a condio de mudana de
sinal dos extremos. Este processo pode ser repetido at
o intervalo ser to pequeno quanto o desejado.
Computao Aplicada Geologia
DEPARTAMENTO DE GEOLOGIA
FACULDADE DE CINCIAS
UNIVERSIDADE DE LISBOA
Matlab Equaes no lineares
Bisseco - Algoritmo
Dados entrada: a, b, tolerncia (tol), n max. iteraes
(maxIter)
Sada: soluo aproximada
Passo 1: Fazer k = 1:maxIter
Passo 2: c = a + (b a) / 2;
Passo 3: Se f(c) = 0 ou |b-a|/2 < tol ento
devolver c, [k]
rotina completada com sucesso,
Passo 4: Se f(a).f(c) > 0 ento a = c seno b = c
Passo 5:
devolver c, [k]
no foi encontrada uma soluo que garanta a
tolerncia especificada
Computao Aplicada Geologia
DEPARTAMENTO DE GEOLOGIA
FACULDADE DE CINCIAS
UNIVERSIDADE DE LISBOA
Matlab Equaes no lineares
Bisseco Erros e convergncia
O mtodo da bisseco converge sempre.
Ao fim de k iteraes o intervalo tem comprimento (b-a)/2k, ou seja,
ba
k
2
Assim, para uma dada tolerncia o nmero de iteraes
necessrias e suficientes para atingir o critrio de paragem o
menor inteiro superior a
log(b a ) log( )
log(2)
Computao Aplicada Geologia
DEPARTAMENTO DE GEOLOGIA
FACULDADE DE CINCIAS
UNIVERSIDADE DE LISBOA
Matlab Equaes no lineares
Newton
O mtodo de Newton, ou Newton-Raphson, um dos
mais conhecidos e poderosos mtodos para a
determinao de razes de equaes no lineares.
Neste mtodo utiliza-se o conhecimento da derivada de f
para se chegar mais rapidamente soluo.
Computao Aplicada Geologia
DEPARTAMENTO DE GEOLOGIA
FACULDADE DE CINCIAS
UNIVERSIDADE DE LISBOA
Matlab Equaes no lineares
Newton
yy
f ( xk )
xk+1
xk
f ( xk )
f ( xk )
x
=
x
ou
tan( ) = f ' ( xk ) =
k
ouseja
seja k +1
f ' ( xk )
xk +1 xk
xx
Computao Aplicada Geologia
DEPARTAMENTO DE GEOLOGIA
FACULDADE DE CINCIAS
UNIVERSIDADE DE LISBOA
Matlab Equaes no lineares
Newton - Algoritmo
Dados entrada: estimativa inicial (x0), tolerncia (tol),
n max. iteraes (maxIter)
Sada: soluo aproximada
Passo 1: Fazer k = 1:maxIter
Passo 2: xk+1 = xk f(xk) / f(xk)
Passo 3: Se |xk+1 - xk| < tol ento
devolver xk+1, [k]
rotina completada com sucesso,
Passo 4: xk = xk+1
Passo 5:
devolver xk, [k]
o mtodo no convergiu
Computao Aplicada Geologia
DEPARTAMENTO DE GEOLOGIA
FACULDADE DE CINCIAS
UNIVERSIDADE DE LISBOA
Matlab Equaes no lineares
Newton Convergncia
Condies
Condiessuficientes
suficientesde
deconvergncia
convergncia
Teorema 1
Seja z uma raiz da equao f(x) = 0. Se f(x), f(x) e f(x) forem
contnuas numa vizinhana de z, f(z) = 0 e x0 estiver suficientemente
prximo de z, ento a sucesso do mtodo de Newton converge para
z.
Computao Aplicada Geologia
DEPARTAMENTO DE GEOLOGIA
FACULDADE DE CINCIAS
UNIVERSIDADE DE LISBOA
Matlab Equaes no lineares
Newton Convergncia
Teorema 2
Seja z uma raiz da equao f(x) = 0. Se:
f(x), f(x) e f(x) forem contnuas num intervalo [a, b] que
contm z,
f(x) no se anular em [a, b],
f(a).f(b) < 0,
f(x)<0 ou f(x) > 0 em [a,b]
ento, se se partir de um ponto x0 [a, b] tal que f(x0).f(x0) 0, a
sucesso do mtodo de Newton converge monotonamente
para z.
