My title - Departamento de Matemática da Universidade do Minho
My title - Departamento de Matemática da Universidade do Minho
My title - Departamento de Matemática da Universidade do Minho
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
5 ÓRBITAS REGULARES E PERTURBAÇÕES 37<br />
A família quadrática. A família quadrática, ou logística (<strong>do</strong> francês “logement” 9 ), é a família<br />
<strong>de</strong> transformações f λ : [0, 1] → [0, 1] <strong>de</strong>fini<strong>da</strong>s por<br />
x ↦→ λx(1 − x)<br />
on<strong>de</strong> o parâmetro λ tem valores no intervalo [0, 4].<br />
Os pontos fixos <strong>de</strong> f λ são 0, que é atrativo quan<strong>do</strong> 0 ≤ λ < 1, e p λ = λ−1<br />
λ<br />
, que é atrativo<br />
quan<strong>do</strong> 1 < λ < 3.<br />
Se λ ∈ [0, 1] então to<strong>da</strong> trajetória (fλ n (x)) converge para 0. De fato, to<strong>da</strong> trajetória é uma<br />
sucessão <strong>de</strong>crescente e limita<strong>da</strong>, logo convergente, e o limite é o único ponto fixo 0.<br />
Se λ ∈ ]1, 3] então to<strong>da</strong> trajetória (fλ n (x)) converge para p λ. De fato, se 1 < λ < 3, existe uma<br />
vizinhança V <strong>de</strong> p λ tal que f λ | V é uma contração e tal que para to<strong>do</strong> x ∈ [0, 1] existe um tempo<br />
n ≥ 0 tal que fλ n (x) ∈ V . O caso em que λ = 3 não é muito diferente ...<br />
5.6 Convergência no méto<strong>do</strong> <strong>de</strong> Newton<br />
Newton method. Seja p ∈ R[x] um polinómio com coeficientes reais. O méto<strong>do</strong> <strong>de</strong> Newton<br />
para <strong>de</strong>terminar as raizes <strong>de</strong> p(x), ou seja, resolver a equação p(x) = 0, consiste em escolher uma<br />
primeira aproximação x 0 , e iterar<br />
x n+1 = x n − p(x n)<br />
p ′ (x n ) .<br />
Ou seja, se x 0 é uma primeira conjectura, uma conjectura melhor é o zero <strong>da</strong> aproximação linear<br />
p(x) ≃ p(x 0 ) + p ′ (x 0 ) · (x − x 0 ).<br />
Se a sucessão (x n ) converge, i.e. x n → r, e se p ′ (r) ≠ 0, então o limite r é um zero <strong>do</strong><br />
polinómio p. Vice-versa, se r é uma raiz <strong>do</strong> polinómio p, e se p ′ (r) ≠ 0, então r é um ponto fixo<br />
<strong>da</strong> transformação<br />
x ↦→ f(x) = x − p(x)<br />
p ′ (x)<br />
A <strong>de</strong>riva<strong>da</strong> <strong>de</strong> f em r é<br />
f ′ (r) = 1 − (p′ (r)) 2 − p(r)p ′′ (r)<br />
(p ′ (r)) 2 = 0<br />
Portanto, r é um ponto fixo atrativo <strong>de</strong> f: a trajetória <strong>de</strong> to<strong>da</strong> conjectura inicial x 0 suficientemente<br />
próxima <strong>da</strong> raiz r converge para r.<br />
Quadratic convergence. De facto, sen<strong>do</strong> a <strong>de</strong>riva<strong>da</strong> f ′ (r) = 0, a raiz é um ponto fixo superatractivo<br />
<strong>de</strong> f, e a convergência é, em geral, mais <strong>do</strong> que exponencial.<br />
Theorem. Let r be a non-critical root of the polynomial p ∈ R[x], i.e. a root where p ′ (r) ≠ 0.<br />
Then Newton’s iterations starting with some x 0 sufficiently near the root r converge to this root,<br />
and the convergence is “quadratic”, i.e. the error ε n = |x n − r| <strong>de</strong>creases as<br />
for some K > 0.<br />
ε n+1 ≤ K · ε 2 n<br />
9 A família quadrática nasceu como mo<strong>de</strong>lo <strong>de</strong> dinâmica <strong>de</strong> populações. Num meio ambiente ilimita<strong>do</strong>, uma<br />
população tem crescimento exponencial, e um mo<strong>de</strong>lo razoável a tempo discreto é z n+1 = λz n, on<strong>de</strong> z n é o<br />
tamanho <strong>da</strong> população no tempo n e λ > 0 é um parâmetro que caracteriza a ”fertili<strong>da</strong><strong>de</strong>” <strong>da</strong> especie. Se o meio<br />
ambiente (a disponibili<strong>da</strong><strong>de</strong> <strong>de</strong> espaço e comi<strong>da</strong>) é limita<strong>do</strong>, parece razoável acrescentar um termo negativo −βzn,<br />
2<br />
on<strong>de</strong> β > 0, que toma conta <strong>da</strong> mortali<strong>da</strong><strong>de</strong> <strong>de</strong>vi<strong>da</strong> à falta <strong>de</strong> recursos (a probabili<strong>da</strong><strong>de</strong> <strong>de</strong> <strong>do</strong>is individuos estar<br />
num mesmo sitio, e portanto disputar a comi<strong>da</strong>, é proporcional a zn). 2 O mo<strong>de</strong>lo é portanto z n+1 = λz n − βzn. 2 A<br />
maior população suporta<strong>da</strong> pelo meio ambiente resulta ser λ/β e, se chamamos x n = z nβ/α a população relativa,<br />
obtemos o mo<strong>de</strong>lo x n+1 = λx n (1 − x n).