Modelagem de Populações por Equações Diferenciais e por Mapas Discretos

1. Modelagem de Populações por Equações Diferenciais e por Mapas Discretos

2. Conteúdo • Modelos discretos para uma única população • Modelos contínuos para uma única população • Modelos para populações interagentes • Estudo da Estabilidade dos Regimes de Equilíbrio. (opcional) • Exemplo de sistema com três espécies e caos.

3. 1. Modelos discretos para uma única população X representa a quantidade de indivíduos em uma população. Se a população for grande podemos tratar X como uma variável contínua. Se a escala de tempo de interesse for da ordem de algumas gerações, e se assumirmos que a cada geração a população é completamente substituída pelos filhos, podemos escrever onde m é o número de descendentes médio por indivíduo. A solução dessa equação é: Assim, se m > 1 a população cresce indefinidamente e se m < 1 ela desaparece.

4. Esse modelo Mauthusiano (Thomas Mauthus – 1798) é claramente muito simplista e irreal. Podemos melhorá-lo acrescentando uma capacidade de suporte K, que representa a população máxima que pode ser sustentada com os recursos disponíveis: A nova taxa de natalidade é Essa equação pode ainda ser simplificada definido-se , i.e., Medindo-se a população em unidades da população máxima. A equação resultante é conhecida como Mapa Logístico e será estudada em detalhes adiante.

5. Equações do tipo são chamadas de mapas discretos. Uma das informações básicas que podemos obter de um mapa diz respeito aos pontos de equilíbrio da dinâmica, i.e., os valores de x que são estacionários. Esse valores são dados pela equação Para o caso Maltusiano isso implica em x=0. Para o mapa Logístico obtemos dois valores: x0=0 e xm=(m-1)/m. O primeiro representa a extinção da espécie, enquanto que o segundo representa uma população cujo valor está entre 0 (para m=1) e 1 (para m muito grande). Note que os valores de m relevantes são restritos a m > 1, senão as populações ficam negativas. A figura na próxima página mostra um exemplo numérico de xn versus n com m = 2.9 e população inicial x0=0.5 (ou X0=K/2). Veja que a população converge para o valor estacionário xm=(2.9-1)/2.9 ~ 0.655. Porque a população não convergiu para x0=0, que também é uma solução de equilíbrio do sistema? A resposta é que x0 é uma solução instável, enquanto que xm é estável! Portanto, não basta achar os valores de equilíbrio, mas temos também que estudar sua estabilidade.

7. 2. Modelos contínuos para uma única população Novamente chamaremos de X a quantidade de indivíduos em uma população. Se a escala de tempo de interesse for da ordem de muitas gerações, como no caso das lebres no Canadá, e se houver superposição de gerações, podemos tratar o tempo como uma variável contínua. Nesse caso escrevemos, por exemplo, pois em um pequeno intervalo de tempo a população não deve mudar muito.

8. Nesse caso é mais conveniente falar da taxa de variação da população com o tempo. Escrevemos que é uma equação diferencial de primeira ordem (só aparece a primeira derivada de X em relação ao tempo. O famoso modelo de Lotka e Volterra, e muitos outros modelos criados depois, é descrito em termos de equações diferenciais. O sinal de aproximação é para lembrar que a derivada só aparece de fato no limite em que e vai a zero.

9. Exemplo: a população na Terra é da ordem de 6 bilhões de pessoas. Estima-se que a cada hora cerca de 12.000 novos indivíduos sejam acrescentados ao planeta, o que da aproximadamente 3 pessoas por segundo. Então a variação da população por unidade de tempo é Escrevendo, como no modelo logístico, , vemos que a taxa R fica dada em termos de pessoas/segundo por pessoa: R indica que cada 1 milhão de pessoas contribui em média com 2 novos habitantes a cada hora.

10. Trecho do livro do Lotka

11. Resolver equações diferencias não é tarefa fácil. Obter uma solução X(t) significa escrever a evolução temporal em termos de funções simples, como o seno, o cosseno e a exponencial. É então importante lembrar do conceito de derivada e rever alguns resultados do Cálculo. Derivada de uma função: f(x) f(x0+dx) f(x0) A interpretação geométrica da derivada é a tangente do ângulo formado entre a reta vermelha (que tangencia f(x) em x=x0) e a horizontal. x0 x0+dx x

12. Formula Importante: Série de Taylor Se dx é pequeno Resolvendo essa equação para f(x0+dx) obtemos a famosa expansão de Taylor da função f(x) em torno de x0:

13. Derivadas de algumas funções simples: Notação: em geral vamos escrever

14. Exemplos: 1 – 2 – 3 – A solução é pois A solução é pois A solução é Mostre que x(t) satisfaz a equação diferencial e que x(0)=x0 . Mostre ainda que no limite em que t vai a infinito, x(t) vai a 1.

15. 3. Modelos contínuos para várias populações O modelo de Lotka-Volterra para interações entre predador x e presa y Se não há interações, b=d=0 o predador P morre de fome, o que é representado pelo coeficiente negativo – c. A presa N por sua vez, cresce indefinidamente se não há predadores, por isso o termo positivo +a. O encontro entre predadores e presas é proporcional ao número de indivíduos de cada tipo, por isso a interação aparece como o produto N P. O coeficiente b mede a taxa com que as presas são comidas e o coeficiente d a taxa com que os alimentos são efetivamente transformados em novos indivíduos predadores.

16. Essas equações podem ser simplificadas com algumas mudanças de escalas: Com isso obtemos

17. Os pontos de equilíbrio são obtidos zerando as taxas de variação de x e de y: • o que resulta em dois valores: • u = v = 0 (extinção) • u = v = 1

18. 4. Estudo da Estabilidade das Soluções de Equilíbrio 4.1 – Mapas discretos 4.2 – Modelos contínuos para uma espécie 4.3 – Modelos contínuos para duas espécies Obs. Essa parte da aula requer um pouco mais de familiaridade com cálculo e pode ser pulada sem problemas. Uma discussão mais qualitativa sobre estabilidade será feita na próxima aula sobre dinâmica caótica.

19. 4.1 - Estabilidade de soluções estacionárias em mapas discretos Para sabermos se uma solução estacionaria qualquer xe é estável ou não, temos que estudar o que acontece com soluções próximas de xe: se elas se aproximarem de xe dizemos que xe é estável. Se elas se distanciarem, xe será instável. Seja então x0 = xe + d x0 uma condição inicial próxima de xe. Então o próximo valor de x, x1, será também próximo de xe, ou seja, podemos escrever x1 = xe + d x1: A condição para que x1 esteja mais próximo de xe do que o ponto inicial x0 é que |d x1| < |d x0|, ou onde f ‘(xe) indica a derivada de f(x) calculada no ponto xe.

20. Para o caso do mapa Logístico temos f ’(x)= m(1-2x). Para o ponto x0=0 obtemos f ‘(0) = m. O ponto x0=0 será instável para m > 1. Para o ponto xe = (m-1)/m obtemos f ‘(xe) = 2 - m. Quando m=1 o valor da derivada é 1. Para 1 < m < 3 a derivada fica entre 1 e -1 e, portanto, seu módulo é menor do que 1. Concluímos que xe é estável nesse intervalo. Para m > 3 obtemos f ‘(xe) < -1 e | f ‘(xe) | > 1 fazendo com que xe fique instável. A figura abaixo mostra xn versus n para m=3.1. A solução não converge mais para xe~0.677. xe

21. 4.2 - Soluções estacionárias e estabilidade para modelos contínuos de uma única espécie. As soluções estacionárias são dadas por Seja xe uma solução estacionária, i.e., f(xe)=0. Como estudar a estabilidade de xe? Novamente olhamos para o comportamento de soluções vizinhas ao ponto de equilíbrio: fazemos x0 = xe + d x0 e calculamos a trajetória x(t) escrevendo x(t) = xe + d x. Se x(t) se aproximar de xe, então xe é estável. Senão xe é instável: Então cuja solução nos já conhecemos: Se f ‘(xe) > 0 d x(t) cresce e x(t) se afasta de xe e o ponto é instável Se f ‘(xe) < 0 d x(t) decresce e x(t) se aproxima de xe e o ponto é estável

22. Exemplo: Nesse caso tem duas soluções: x0 = 0 e x1 = 1. A derivada de f(x) é Então f ‘ (0) = m > 0 i.e., x0 = 0 é instável Então f ‘ (1) = - m < 0 i.e., x1 = 0 é estável

23. 4.3 - Soluções estacionárias e estabilidade para modelos contínuos de duas espécies. Solução de equilíbrio (X1e , X2e) soluções vizinhas ao equilíbrio X1 = X1e + d X1, X2 = X2e + d X2

24. ou onde: Essa equação é muito parecida com a equação para uma única população, mas Z é um vetor e A uma matriz. De qualquer forma, a equação é linear e podemos resolve-la procurando uma solução da forma exponencial. Tentamos Z(t) = Z0 el t. Como a derivada em relação ao tempo é lZ(t) obtemos a relação

25. Se a matriz (A - l 1) puder ser invertida, podemos multiplicar os dois lados pela sua inversa e vamos obter Z(t)=0. Essa solução não interessa, pois refere-se ao ponto de equilíbrio apenas, e não a sua vizinhança. A condição para que outras soluções existam é que não exista a inversa, ou seja, que o determinante de (A - l 1) seja nulo: Os valores do expoente podem ser reais (positivos ou negativos) ou mesmo complexos:

26. reais, ambos positivos; (b) reais, ambos negativos; (c) reais, um positivo e um negativo • (d) complexos parte real positiva; (e) complexos parte real negativa;

27. 5 - Um sistema com três espécies.