19.01.2015 Views

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

SHOW MORE
SHOW LESS

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).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!