Modelo de Malthus ECUUU

Modelo de Malthus Fue al parecer Euler quien desarrolló los primeros modelos de población, pero comúnmente se atribuye a

Views 63 Downloads 0 File size 158KB

Report DMCA / Copyright

DOWNLOAD FILE

Recommend stories

Citation preview

Modelo de Malthus Fue al parecer Euler quien desarrolló los primeros modelos de población, pero comúnmente se atribuye a Malthus 1 el desarrollo y análisis del primer modelo de evolución de 𝑝! (𝑡), según el cual Es decir, en cada instante la rapidez de cambio de la población es proporcional al total de la población presente. Por ejemplo, si 𝑃(𝑡) > 0 𝑦 𝑃(𝑡) creciente, esto implica que k > 0. Resolvemos la ecuación diferencial (3.5): 𝑑𝑃 𝑑𝑡

= 𝑘𝑃 →

𝑑𝑃 𝑃

= 𝑘𝑑𝑡

Integrando se tiene ∫

𝑑𝑃 𝑑𝑡

= ∫ 𝑘𝑑𝑡 → 𝑙𝑛𝑃 = 𝑘𝑡 + 𝑐1 → 𝑃 = 𝑒 𝑘𝑡+𝑐1 = 𝑒 𝑘𝑡 𝑒 𝑐1 = 𝑒 𝑘𝑡 𝑐

𝑃(𝑡) = 𝑐𝑒 𝑘𝑡 Ésta es la solución general de la ecuación diferencial (3.5). Es común conocer la población inicial, 𝑃(0) = 𝑝0 . Con esto podemos calcular la constante C: 𝐏(𝟎) = 𝐩𝟎 = 𝐜𝐞𝐤.𝟎 = 𝐂 → 𝐂 = 𝐩𝟎 → 𝐏(𝐭) = 𝐩𝟎 𝐞𝐤𝐭 Para calcular k es necesario conocer la cantidad de población existente en un tiempo 𝑡1 > 𝑡0 , digamos 𝑝(𝑡1 ) = 𝑝1 𝐩

𝐩

𝐩(𝐭 𝟏 ) = 𝐩𝟏 = 𝐩𝟎 𝐞𝐤𝐭 𝟏 → 𝐩𝟏 = 𝐞𝐤𝐭 𝟏 → 𝐥𝐧 (𝐩𝟏 ) = 𝐤𝐭 𝟏 𝟎

→𝐤=

𝟎

𝐥𝐧𝐩𝟏 −𝐥𝐧𝐩𝟎 𝐭𝟏

Observaciones: 1. 𝑡𝑑 Si es el tiempo en el que la población se duplica, 𝑃(𝑡𝑑 ) = 2𝑝0; entonces tenemos, de acuerdo con la ecuación previa: 𝐩(𝐭 𝐝 ) = 𝟐𝐩𝟎 = 𝐩𝟎 𝐞𝐤𝐭 𝐝 → 𝟐 = 𝐞𝐤𝐭 𝐝 → 𝐥𝐧 𝟐 = 𝐤𝐭 𝐝 Y ahora: 𝐤=

𝐥𝐧𝟐 𝐭𝐝

𝐲 𝐭𝐝 =

𝐥𝐧𝟐 𝐤

Lo anterior indica que el tiempo para que una población se duplique no depende de la cantidad inicial de la misma. 2. Si se proporcionan𝑝(𝑡1 ) = 𝑝1 𝑦 𝑝(𝑡2 ) = 𝑝2 para dos tiempos t1 < t2, obtenemos los siguientes resultados: 𝐏(𝐭 𝟏 ) = 𝐩𝟏 = 𝐜𝐞𝐤𝐭 𝟏 𝐩(𝐭 𝟐 ) = 𝐩𝟐 = 𝐜𝐞𝐤𝐭 𝟐

Para resolver este sistema de dos ecuaciones con dos incógnitas, C y k, dividimos la segunda ecuación entre la primera y obtenemos: 𝐩𝟐 𝐩𝟏

𝐜𝐞𝐤𝐭 𝟐

=

𝐜𝐞𝐤𝐭 𝟏

=

𝐞𝐤𝐭𝟐 𝐞𝐤𝐭𝟏

= 𝐞𝐤(𝐭 𝟐 −𝐭 𝟏 )

Entonces: 𝒑

𝐥𝐧 (𝒑𝟐 ) = 𝒌(𝒕𝟐 − 𝒕𝟏 ) 𝟏

De donde: 𝐤=

𝐥𝐧𝐩𝟐 −𝐥𝐧𝐩𝟏 𝐭 𝟐 −𝐭 𝟏

Tenemos también, usando (3.7) que 𝐜 = 𝐩𝟏 𝐞−𝐤𝐭 𝟏 como 𝐩(𝐭) = 𝐜𝐞𝐤𝐭 𝟏 𝒑(𝒕) = 𝒑𝟏 𝒄𝒆𝒌(𝒕−𝒕𝟏 ) Con k dado por (3.8). 3.

Hemos mencionado que la derivada

𝑑𝑃 𝑑𝑡

es la rapidez de cambio de la población

𝑃(𝑡) A esta derivada también se le denomina tasa de cambio de la población. De aquí surge una expresión frecuentemente usada en los problemas de población: tasa de crecimiento. La tasa de crecimiento de una población en cierto tiempo T se define como la razón 𝑃!(𝑡) 𝑃(𝑡)

y se da comúnmente en términos porcentuales anuales. Es decir, la tasa de

crecimiento de una población es precisamente la constante de proporcionalidad 𝒌=

𝑷!(𝒕) 𝑷(𝒕)

: Recuérdese que 𝑷! (𝒕) = 𝒌𝒑(𝒕)

4. Así como 𝑝(𝑡) = 𝑝0 𝑒 𝒌𝒕 nos sirve para calcular la población creciente de una comunidad, es posible utilizarla para calcular una población que disminuye al paso del tiempo. Sólo debemos tener presente que: 𝒅

a. Si la población 𝑷(𝒕) aumenta, entonces 𝒅𝒕 𝑷(𝒕) = 𝒌𝑷(𝒕) 𝒄𝒐𝒏 𝒌 > 𝟎. 𝒅

b. Si la población 𝑷(𝒕) disminuye, entonces 𝒅𝒕 𝑷(𝒕) = 𝒌𝑷(𝒕) 𝒄𝒐𝒏 𝒌 < 𝟎.

Ejemplo En un cultivo de bacterias, se estimó que inicialmente había 150 bacterias y 200 después de una hora (h). Suponiendo una rapidez de crecimiento proporcional a la cantidad de bacterias presente, determinar: 1. 2. 3.

La cantidad de bacterias después de t horas. La cantidad de bacterias después de 2 h. El tiempo que debe transcurrir para que la población se triplique.

1. Si 𝑃(𝑡) es la cantidad de bacterias presentes después de 𝑡 horas, entonces 𝑃(0) = 𝑃0 = 150 𝑌 𝑃(1) = 𝑃1 = 200. Luego, P (t) está dada por la solución del PVI:

𝑃! (𝑡) = 𝑘𝑃(𝑡).

𝑐𝑜𝑛 𝑃(0) = 150 𝑌 𝑎𝑑𝑒𝑚𝑎𝑠 𝑃(1) = 200

Puesto que𝑃(𝑡) = 𝑐𝑒 𝑘𝑡 , 𝑠𝑒 𝑡𝑖𝑒𝑛𝑒 y remplazo 𝑃(0) = 150 → 𝑐𝑒 0 = 150 → 𝐶 = 150 → 𝑃(𝑡) = 150𝑒 𝑘𝑡 𝑃(1) = 200 → 150𝑒 𝑘 = 200 → 𝑒 𝑘 =

200 4 4 = → 𝑘 = ln ( ) = 0.2877 150 3 3

→ 𝑃(𝑡) = 150𝑒 0.2877𝑡 Que es la solución del PVI y es la cantidad de bacterias después de 𝑡 horas 2. La cantidad de bacterias después de 2 h es 𝑷(𝒕) = 𝟏𝟓𝟎𝒆𝟎.𝟐𝟖𝟕𝟕𝒕 Entonces diremos que el tiempo igual 2h 𝑷(𝟐) = 𝟏𝟓𝟎𝒆𝟎.𝟐𝟖𝟕𝟕(𝟐) =266.6666 𝑷(𝟐) = 𝟐𝟔𝟔 𝒃𝒂𝒄𝒕𝒆𝒓𝒊𝒂𝒔 3. Para que la población se triplique: 𝑷(𝒕) = 𝟑𝑃0 → 𝟏𝟓𝟎𝒆𝟎.𝟐𝟖𝟕𝟕𝒕 = 𝟑(𝟏𝟓𝟎) → 𝒆𝟎.𝟐𝟖𝟕𝟕𝒕 = 𝟑 → 𝟎. 𝟐𝟖𝟕𝟕𝒕 = 𝒍𝒏(𝟑) →𝑡=

ln(3) = 3.8186ℎ → 𝑡 = 3ℎ𝑜𝑟𝑎𝑠, 49 𝑚𝑖𝑛𝑢𝑡𝑜𝑠, 7𝑠𝑒𝑔𝑢𝑛𝑑𝑜𝑠 0.2877