Exercicis de models amb equacions diferencials

Pas 1 — el balanç. El capital creix per dos motius: l'interès, proporcional al que ja hi ha, i l'aportació, constant.

$K' = \frac{3}{100}K$

Pas 2 — resoldre la lineal. Factor integrant $e^{-\frac{3}{100}t}$, i surt

$K(t) = 5000e^{\frac{3}{100}t}$

Pas 3 — avaluar a $T = 5$.

$K(5) \approx 5809.1712\ \text{€}$

És un model idealitzat amb taxa fixa i temps continu: serveix per entendre l'estructura de l'interès compost, no per predir rendiments reals.

Pas 1 — discriminant.

$c^{2} - 4mk = 6^{2} - 4\cdot1\cdot9 = 0$

Pas 2 — llegir-lo. Negatiu: arrels complexes, el sistema oscil·la amb amplitud decreixent. Zero: arrel doble, torna a l'equilibri sense oscil·lar i el més ràpid possible. Positiu: dues arrels reals, torna sense oscil·lar però més lentament.

Aquí val $0$, o sigui Crític.

El detall contraintuïtiu: passat el punt crític, posar més amortiment fa que el sistema torni a l'equilibri més lentament, no més ràpid. Per això els amortidors es dissenyen a prop del cas crític.

Pas 1 — el balanç. Variació = entrades menys sortides. No entra sal, i la que surt és el cabal per la concentració del dipòsit (mescla homogènia):

$M' = 0 - 4\cdot\frac{M}{100} = -\frac{1}{25}M$

Comprova les unitats: kg/min als dos costats. El volum és constant perquè entra i surt el mateix cabal.

Pas 2 — resoldre. És l'exponencial de sempre:

$M(t) = 20e^{-\frac{1}{25}t}$

Pas 3 — avaluar a $t = 30$.

$M(30) = 20e^{-\frac{6}{5}} \approx 6.0239\ \text{kg}$

El temps de mitja vida és $\frac{\ln 2}{\frac{1}{25}} \approx 17.33$ minuts, i no depèn de la sal inicial.

Truc que estalvia tot el càlcul. A llarg termini el dipòsit arriba a l'equilibri, o sigui $M' = 0$. Amb cabal $q$,

$0 = q\cdot1 - q\cdot\frac{M}{300} \quad\Longrightarrow\quad M = 300\cdot1 = 300.0000\ \text{kg}$

El cabal se simplifica: no afecta el valor final, només la velocitat amb què s'hi arriba.

Interpretació. Té tot el sentit: al final el dipòsit acaba tenint exactament la concentració del que hi entra, i la sal total és el volum per aquesta concentració. Fer aquesta comprovació abans o després de resoldre és el millor control de qualitat d'un problema de mescles.

Pas 1 — l'equació. La llei de Newton diu que la velocitat de refredament és proporcional a la diferència amb l'ambient:

$T' = -\frac{1}{20}(T - 20)$

Pas 2 — el canvi que la fa fàcil. Si crides $u = T - 20$, aleshores $u' = -\frac{1}{20}u$, l'exponencial de sempre. Desfent el canvi:

$T(t) = 20 + (80 - 20)e^{-\frac{1}{20}t}$

Pas 3 — avaluar.

$T(20) = 20 + 60\,e^{-1} \approx 42.0728\ \text{°C}$

Comprovació de sentit comú: el resultat ha d'estar entre $20$ i $80$, i acostar-se a l'ambient sense arribar-hi mai. Si et surt per fora, has perdut un signe.

Pas 1 — la constant. De la condició inicial surt

$A = \frac{K}{N_0} - 1 = \frac{500}{50} - 1 = 9$

Pas 2 — la fórmula.

$N(t) = \frac{500}{1 + 9e^{-\frac{1}{2}t}}$

Pas 3 — avaluar a $t = 2$.

$N(2) \approx 115.9847$

Lectura del model. Els equilibris són $0$ (inestable) i $500$ (atractor): tota població positiva acaba tendint a $500$. I el creixement absolut és màxim a $N = K/2 = 250$, que no s'ha de confondre amb el creixement per capita, que decreix sempre.

Pas 1 — el balanç. El capital creix per dos motius: l'interès, proporcional al que ja hi ha, i l'aportació, constant.

$K' = \frac{3}{100}K$

Pas 2 — resoldre la lineal. Factor integrant $e^{-\frac{3}{100}t}$, i surt

$K(t) = 5000e^{\frac{3}{100}t}$

Pas 3 — avaluar a $T = 10$.

$K(10) \approx 6749.2940\ \text{€}$

És un model idealitzat amb taxa fixa i temps continu: serveix per entendre l'estructura de l'interès compost, no per predir rendiments reals.

Pas 1 — discriminant.

$c^{2} - 4mk = 3^{2} - 4\cdot1\cdot9 = -27$

Pas 2 — llegir-lo. Negatiu: arrels complexes, el sistema oscil·la amb amplitud decreixent. Zero: arrel doble, torna a l'equilibri sense oscil·lar i el més ràpid possible. Positiu: dues arrels reals, torna sense oscil·lar però més lentament.

Aquí val $-27$, o sigui Subamortit (oscil·la).

El detall contraintuïtiu: passat el punt crític, posar més amortiment fa que el sistema torni a l'equilibri més lentament, no més ràpid. Per això els amortidors es dissenyen a prop del cas crític.

Pas 1 — el balanç. Variació = entrades menys sortides. No entra sal, i la que surt és el cabal per la concentració del dipòsit (mescla homogènia):

$M' = 0 - 5\cdot\frac{M}{100} = -\frac{1}{20}M$

Comprova les unitats: kg/min als dos costats. El volum és constant perquè entra i surt el mateix cabal.

Pas 2 — resoldre. És l'exponencial de sempre:

$M(t) = 10e^{-\frac{1}{20}t}$

Pas 3 — avaluar a $t = 10$.

$M(10) = 10e^{-\frac{1}{2}} \approx 6.0653\ \text{kg}$

El temps de mitja vida és $\frac{\ln 2}{\frac{1}{20}} \approx 13.86$ minuts, i no depèn de la sal inicial.

Truc que estalvia tot el càlcul. A llarg termini el dipòsit arriba a l'equilibri, o sigui $M' = 0$. Amb cabal $q$,

$0 = q\cdot2 - q\cdot\frac{M}{300} \quad\Longrightarrow\quad M = 300\cdot2 = 600.0000\ \text{kg}$

El cabal se simplifica: no afecta el valor final, només la velocitat amb què s'hi arriba.

Interpretació. Té tot el sentit: al final el dipòsit acaba tenint exactament la concentració del que hi entra, i la sal total és el volum per aquesta concentració. Fer aquesta comprovació abans o després de resoldre és el millor control de qualitat d'un problema de mescles.