Módulo II (a). Función Derivada y Diferenciación
- Comisión Ecología: Dr. Leandro Giordano / Prof. Silvana Ávila
- Hacer uso de las herramientas de cálculo de una variable para la formulación de modelos matemáticos en la estimación de variables biofísicas, mediante las técnicas de derivación e integración (i.e. cálculo de tasas/valores extremos y de valores acumulados/valor medio, estimación de áreas y volúmenes), incorporando las nociones de diferencial de una variable y una función, por un lado, y de integración definida, por otro.
-
Problemas de Extremos (máximos/mínimos, locales/globales, puntos de inflexión)
- Ejemplos: Estimar rendimientos o puntos de localización óptimos. Estimar el pico y el tiempo al pico de distintas señales (e.g. concentración de un contaminante, caudal líquido en una sección fluvial). Estimar puntos de inflexión (función curvatura).
-
Estimación del área bajo una curva y del valor medio. Estimación de volúmenes a partir de aplicación de la técnica de integración sobre curvas Área-Elevación.
- Ejemplos: Estimar el stock o resevas acumuladas o su valor medio en un intervalo de tiempo, a partir del conocimiento de sus tasas de varaición. Extraer Curvas Área-Volumen a partir de Información Topográfica. Resampleo de señales a distinta resolución temporal.
En el desarrollo de esta sección se asumirá que se está familiarizado con los conceptos básicos de función
$f(x):R \to R$ de una variable real$x$ (e.g. definición, univocacidad/biyectividad, teorema de continuidad de Bolzano, funciones polinómicas y funciones trascendentales - exponenciales, logarítmicas, trigonométricas -, límites de una función, funciones invertibles y composición de funciones), propios del desarrollo de una Introducción a las Matemáticas Superiores o de las Matemáticas Generales. En todo caso, se recomienda realizar un breve repaso por los ítems mencionados.
Es sabido que muchos procesos ecológicos varían continuamente en el espacio o en el tiempo, alterando las magnitudes de los factores involucrados. Por ejemplo, se sabe que la cantidad de biomasa en un organismo irá variando de acuerdo a su estado ontogenético, el cual depende del tiempo. Luego, si el proceso puede describirse a partir de una relación, al menos unívoca, entre el tiempo y la biomasa, y si llamamos
En este ejemplo, además debiera considerarse que
$t$ estará definido para un subconjunto en el intervalo$[0,T]$ (en donde$T$ es la extensión esperada del ciclo de vida del organismo), tanto como$M$ estará definido en el intervalo$[0,M_0]$ (en donde$M_0$ será la biosama máxima que pueda alcanzar dicho organismo). Así, sería más apropiado formular$M(t):[0,T] \to [0,M_0]$ .
Particularmente, el interés en la formulación de este tipo de modelos puede residir en responder a preguntas tales cómo:
¿Cuál es la tasa de crecimiento de la función
$M(t)$ ? ¿Es caracterizable por otra función$f(t)$ ?¿Hay algún instante
$t$ para el cual esta tasa sea máxima, nula o mínima?¿Para qué instantes
$t_1,t_2,...,t_n$ la función$M(t)$ presenta máximos o mínimos locales/globales?¿Cambia la curvatura de la función
$M(t)$ en algún instante?
Por lo general, una función
-
(a) Como una tabla con dos columnas, una para la variable
$x$ y otra para$f(x)$ (i.e. una matriz con 2 vectores columnas$\mathbf{X}=[x_1,...,x_n]$ y$\mathbf{Y}=[f(x_1),...,f(x_n)]$ ) -
(b) Como gráfica en un sistema de coordenadas rectangulares
$(x,y)$ con$y=f(x)$ -
(c) Como expresión análitica o su forma explícita. Por ejemplo, la función lineal
$f(x)=a_1 x+a_0$
Cualquiera de estas 3 representaciones es una función matemática, siempre y cuando se satisfaga el principio básico que establece que a cada elemento del conjunto dominio
La Fig.1 muestra la representación en un sistema de coordenadas rectangulares de una función
Fig.1 - Representación gráfica de una función f(x) en un sistemas de coordenadas rectangular. Puede notarse que un incremento finito de la variable x se asocia a un incremento finito de la función f(x). El cociente de estos incrementos representa la tasa promedio de cambio de la función f(x) entre los dos valores x considerados (tasa promedio en intervalo en x)
En principio, introduzcamos la noción de incremento o diferencia finita de la variable
A la vez, introduxcamos el concepto de incremento o diferencia finita de la función
Resulta evidente que el cociente de estos incrementos finitos, definido mediante:
es la pendiente de la recta secante que une los puntos
Esto es, la tasa de incremento medio de
$f(x)$ en el intervalo con límites$x_1$ y$x_2$ .
Luego, se introduce el concepto de función derivada
La interpretación geométrica es la siguiente: si
$\Delta x \to 0$ , luego también$\Delta f(x) \to 0$ . Consecuentemente, en el límite, la recta secante de la Fig. 1 se transformará en una recta tangente, aquella que pasa solamente por el punto con coordenadas$x_1$ y$f(x_1)$ . Luego, el límite del cociente de estos incrementos finitos representa la tasa de cambio local o instantánea de$f(x)$ en$x_1$ , cuyo valor será igual a la pendiente de esta recta tangente a$f(x)$ en$x_1,f(x_1)$ . En la Fig. 2 se presenta una animación que da cuenta de esto y lo pone en relación con la noción d1. e diferencial.
Asimismo, se introduce la noción de diferencial de una variable
En otras palabras, el diferencial de una variable es una diferencia infinitesimal. Coloquialmente, 'tan pequeña como podamos imaginar y aun más, hasta dejar de imaginarlo'. Estrictamente, una diferencia infinitamente pequeña.
Del mismo modo, se introduce la noción de diferencial de una función
Y así puede observarse también que se sostiene la siguiente igualdad:
que establece que la derivada
Fig.2 - Derivada como límite del cociente de incrementales de la función f(x) y la variable independiente x, relación con el cociente de diferenciales
Asimismo, puede observase que se deduce
$f(x)'dx=df(x)$ . Esto es, se establece que el incremento infinitesimal de una función es igual al producto entre la derivada de la función$f(x)$ y el diferencial de la variable$dx$ . Esto será fundamental para entender en mayor profundidad, posteriormente, el cálculo integral (puesto que$f(x)'dx$ será el integrando).
Luego, la pregunta en términos de aplicación será por qué es conveniente la adopción de estas nociones, de la definición de la derivada en términos de un cociente de diferenciales. La respuesta directa es, la diferenciación es una operación y luego
Así también, como la derivada
- Toda función
$f(x)'$ puede derivarse recursivamente al menos hasta que el valor de la función derivada obtenida sea una constante$k$ , como ya veremos (la derivada de una constante es igual a 0, pues la pendiente es 0).
Así, la
$n$ -énesima derivada de$f(x)$ será aquella función$g(x)$ que se obtiene al iterar(repetir)$n$ -veces la operación de derivación sobre$f(x)$ (una vez que se obtiene una derivada, se vuelve a realizar al operación hasta cumplir con$n$ operaciones).
Una aplicación de importancia de esto es el estudio de la segunda derivada de
Así como primer derivada se denota mediante
$f(x)'$ , la segunda se denota mediante$f(x)''$ y así en adelante. Asimismo, la notación diferencial es más conveniente, pues permite una mayor generalización. En esta se define a la derivada de orden$n$ como aquella función definida mediante la operación:$$\dfrac{d^n f(x)}{dx^n}$$ .
Antes de las propiedades de la diferenciación y su aplicación para el cómputo de derivadas, es muy posible que para esto haya tenido que desarrollar el cálculo del límite de
Por ejemplo, si consideramos:
Para esto se procedía de la siguiente forma. Se sabe que:
Luego si reemplazamos
Desarrollando el binomio y el producto
Aplicando propiedad distributiva y eliminando los factores con igual valor se obtiene:
Luego, finalmente puede establecerse la igualdad:
Los pasos del método pueden repertirse para otras funciones. En lenguaje coloquial, la estrategia consiste en eliminar a
Ciertamente, la aplicación de este método para polinomios de bajo orden puede no ser muy dificultosa para personas con ganas y tiempo (si quiere probar con un ejemplo inicial le recomendamos el cálculo para la función lineal. Intuitivamente podemos pensar que es una constante e igual a su pendiente, pues la recta tangente a los puntos de una recta, es la misma recta).
Como siempre hay entusiastas de la ciencia, muchas personas se dedicaron a estudiar estos límites paralas funciones polinómicas, trignométricas y exponenciales. Gracias a ello, siempre se pudo disponer de tablas de derivadas. Aun así, la introducción del uso del operador diferencial y sus reglas (i.e. las mismas de la derivación), nos permitirá fácilmente deducir derivadas, puesto el carácter directo de la operación. Una vez calculado el diferencial de una función
Al inicio de este módulo se presentó la siguiente igualdad, introduciendo el concepto de diferencial de una variable
Que como se ha visto, también puede expresarse como:
Quedando en claro que la diferenciación es una operación y, por tanto,
Ya se mencionó que esta expresión tiene importantes consecuencias operativas en el cálculo (el diferencial de una función es el integrando en el cálculo integral). Específicamente se desarrollarán a continuación las vinculadas al cómputo de la derivada. Aún así, no quisiera dejar de notarse que un aspecto bello de esta igualdad es que zanja algunos aspectos en la guerra del cálculo, librada por Newton y Leibniz (ver aquí!), puesto que pone en relación la notación operativa de ambos autores. Luego de terminar la lectura sobre diferenciales, bien puede entretenerse (aquí!)
En primer lugar:
- El diferencial de una constante escalar
$k$ es siempre:
Esto guarda sentido, puesto que si el valor es constante, cualquier diferencia finita es igual a 0 y, de ahí, que el diferencial también lo sea.
Además:
- El diferencial del producto de un escalar
$k$ y una función$f(x)$ es igual al producto de$k$ por el diferencial de la función$f(x)$ :
- La diferenciación de funciones es distributiva y conmutativa en la suma (resta). Luego:
- El diferencial del producto de 2 funciones queda dado mediante (regla del producto):
- Finalmente, se introduce la regla del cociente:
Estas propiedades serán utilizadas para simplificar el cómputo de la derivada.
A continuación, se introducen las definiciones de la operación diferenciación en funciones polinómicas, trigonométricas, potenciación y logaritmación.
Diferenciación de Polinomios
Para el caso de los polinomios, se introduce el caso general
Nótese que el caso más simple será
$f(x)=x$ (función identidad de la variable independiente), de manera tal que$g(x)=f(x)^n=x^n$ , en dicho caso. En esta sección utilziamos$f(x)$ en lugar de$x$ pues se está trabajando con el operador 'diferencial de una función'. Asimismo, si$f(x)=x$ , luego queda claro que$df(x)=dx$ .
Luego, la diferenciación de polinomios queda dada mediante:
Aquí surgen 2 aspectos interesantes. El primero, es que si recuerda la tabla de derivadas de los polinomios, la forma es exactamente la misma (ahí la relación entre operaciones). El segundo, es que siempre al diferenciar una función, el resultado es el producto de otra función y el diferencial de la primera. Esto último tiene consecuencias importantes en el cómputo de la derivada, al eliminar el denominador del cociente.
Veamos un ejemplo, para el caso más simple
$f(x)=x$ , entonces$dg(x)=d(f(x)^n)=nx^{n-1}dx$ .
Y de ahí,
$\dfrac{dg(x)}{dx}=\dfrac{nx^{n-1}dx}{dx}=nx^{n-1}$ , tal cual puede observarse en una simple tabla de derivadas.
Diferenciación de funciones trigonométricas
Si es necesario repasar sobre la definición de funciones trigonométricas, siempre es conveniente referirse a la circunferencia de radio 1 (ver aquí!). En caso que quieran consultarse las identidades trigonométricas (relaciones entre operaciones), útiles para la simplificación de problemas, pueden consultar en páginas como esta.
Para la forma general
Así, en el caso maś simple si
$f(x)=x$ , luego$d(sen(x))=cos(x)dx$ . De forma tal que:$$\dfrac{d(sen(x))}{dx}=cos(x)$$ Tal cual se puede apreciar en una simple tabla de derivadas. Esto también podrá verificarlo de manera análoga en todas las definiciones posteriores.
Asimismo para
Finalmente para
Diferenciación de funciones potenciales/exponenciales
Sea
Así, para el caso más simple
$f(x)=x$ , luego$dg(x)=d(a^{f(x)})=ln(a)a^xdx$ y, de ahi que:$$\dfrac{d(a^x)}{dx}=ln(a)a^{x}$$ como se puede apreciar en una simple tabla de derivadas.
Queda claro que si el exponente
Así, para el caso más simple
$f(x)=x$ , puede comprobrase que (verifique):$$\dfrac{d(e^x)}{dx}=e^x$$ .
Finalmente para la función inversa de la exponenciación, la logaritmación
Así, puede demostrarse que para el caso más simple
$f(x)=x$ , será (verifique):$$\dfrac{dg(x)}{dx}=\dfrac{d(ln(x))}{dx}=\dfrac{1}{x}$$
Pues bien, con la definición de estas operaciones y de las propiedades fundamentales de la diferenciación de funciones, se ha visto que para el caso más simple
Regla de la cadena
Sea
Que en términos operativos equivale a:
Que también puede expresarse como:
En pocas palabras, en dicha situación la derivada
Ejemplo de aplicación de la regla de la cadena
Sea
$g(x)=cos(u)$ con$u=f(x)$ y$f(x)=ax^n+b$ . Luego, puede plantearse de acuerdo a lo formulado:$$\dfrac{dg(x)}{dx}=\dfrac{dg(u)}{du}\dfrac{du}{dx}=\dfrac{dg(u)}{du}\dfrac{df(x)}{dx}$$ Primeramente puede resolverse:$$\dfrac{dg(u)}{du}=\dfrac{-sen(u)du}{du}=-sen(u)$$ Luego, mediante la aplicación de las propiedades y definiciones fundamentales de la diferenciación de funciones, también se procede a resolver:$$\dfrac{du}{dx}=\dfrac{df(x)}{dx}=\dfrac{d(ax^n+b)}{dx}$$ $$\dfrac{du}{dx}=\dfrac{d(ax^n+b)}{dx}=\dfrac{d(ax^n)}{dx}+\dfrac{d(b)}{dx} $$ $$\dfrac{du}{dx}=anx^{n-1}\dfrac{dx}{dx}+0$$ $$\dfrac{du}{dx}=anx^{n-1}$$ Finalmente, reemplazando ambos resultados en la igualdad formulada originalmente, y sustituyendo mediante la definición explícita de$u=f(x)$ se obtiene:
$$\begin{align} \dfrac{dg(x)}{dx}&=-ansen(u)x^{n-1} \\ \dfrac{dg(x)}{dx}&=-ansen(ax^n+b)x^{n-1} \\ \end{align}$$ O lo que es lo mismo: $$g(x)'=-ansen(ax^n+b)x^{n-1}$$ En los slides de los talleres precedentes (NotasdeAula3.pdf) podrán obserar un desarrollo alternativo de esta regla, aplicando directamente las propiedades de diferenciación de forma recursiva, primeramente al diferenciar
$g(u)$ y al reemplazar$u=f(x)$ y continuar diferenciado hasta reducir la diferenciación a$dx$ , el diferencial de la variable independiente de la función (composición) explícita$g(x)$ (debe tenerse en cuenta que hay una errata en la primer linea del último paso del desarrollo - una$n$ de más, que luego se corrige -).
Ejemplo 1. Determinación del tiempo al pico y su relación con el tiempo de residencia para un modelo de tránsito de caudales (función impulso-respuesta)
Para la simulación del tránsito de caudales o de contaminantes sobre una red fluvial, más específicamente para la elaboración de funciones de transferencia en sistemas lineales, es común asumir una función de impulso-respuesta, que muestra cómo se distribuiría temporalmente un impulso de escorrentía generado en una cuenca sobre el punto de cierre de esta, el traslado de un impulso de caudal registrado aguas arriba hacia un punto situado aguas abajo o el vuelco instantáneo de un contaminante, desde un punto situado aguas arriba en su arribo a un sitio aguas abajo, entre otras situaciones (e.g. también puede simular la distribución temporal en respuesta a la ingesta instantánea de una sustancia).
Para el tránsito de caudales o de sustancia contaminantes, una aproximación para obtener dicha función consiste en conceptualizar la red fluvial, o el tramo considerado, como un sistema de
Fig.3 - Modelo de tránsito de caudal o sustancias contaminantes (función impulso-respuesta) sobre la base de una cascada de reservorios lineales en serie (ver Nash, J. E., The form of the instantaneous unit hydrograph, 1957). Se asume que un impulso instantáneo de información o masa ingresa al sistema (un pulso de duración infinitesimal) y el tránsito por la cascada lo propaga en el tiempo, dando por resultado una función de distribución temporal. Nótese que a medida que se incrementa la cantidad
Como veremos luego, cuando se desarrollen las Ecuaciones Diferenciales Ordinarias Lineales de Primer Grado, esta asunción conduce al siguiente modelo matemático (por ahora, tómelo por hecho!):
En donde
Así por ejemplo, la Fig. 4 muestra como un ingreso instantáneo de caudal de 500
Fig.4 - Modelo de tránsito de un caudal instantáneo de 500 m³/s sobre un tramo compuesto por 3 reservorios lineales, cada uno con un tiempo de residencia promedio de 2.5 días. El modelo representa como este ingreso instantáneo se distribuiría en su arribo a la salida del sistema (la descarga del 3 reservorio). Nótese el claro desarrollo de un pico, luego de la fase de ascenso y la posterior recesión (extinción exponencial), luego de un punto de inflexión o cambio de curvatura, durante la fase de descenso. Esto es común al apreciar 'tránsitos de avenidas' y, por tanto, el modelo puede utilizarse en varias situaciones, sobre todo, en aquellas que pueda asumirse una respuesta lineal (no en el caso de propagación rápida por rotura de una presa!).
Para el caso de este modelo, bien podríamos indagar si bien pudiera establecerse alguna relación directa entre el tiempo al pico
En efecto, si el tiempo al pico
$t_p$ es el instante para el cual la función$Q(t)$ se maximiza luego debe satisfacerse:
Por lo que primero se procede a calcular la primer derivada. Aplicando las propiedades de diferenciación podemos desarrollar esto mediante:
Aplicando la regla del producto y la definición de diferenciación de un polinomio y de la función exponencial, se obtiene:
De manera tal que igualando a 0 la derivada y utilizando el factor común
Luego, por un lado se sabe que:
si y sólo si
Concluyendo que
Como puede observarse, en este tipo de función este valor será un máximo absoluto, por tanto bien puede corrobarse que para este modelo, siempre será:
Así, para el caso del ejemplo de la Fig. 4 si
En principio diremos que un espacio isotrópico es aquel para el cual la variación de una magnitud
Fig.5 - Modelo de localización de la Teoría del Lugar Central. La renta urbana depende únicamente de la distancia al centro de servicios (CBD), sin importar la dirección considerada.
En consecuencia, podrá decirse que el modelo se simplifica a una relación
- El costo de de desplazarse al centro desde el lugar de residencia. El cual se incrementa a medida que la distancia
$x$ también crece. En principio consideraremos que este costo se incrementa linealmente a una tasa$a$ , esto es:
En donde
- El costo de residencia o de establecerse en un sitio particular. Este costo dependerá de varios factores. En una primera y simple aproximación supondremos que este decrece exponencialmente desde el centro a la periferia (se asume que la competencia por el espacio es mayor), siguiendo la ley:
En donde
El costo total
Luego, si quiéramos evaluar en qué sitio se minimizaría la renta y, por tanto, el costo de situarse en el mismo, se encontraría un sitio óptimo. Nuevamente, para esto se procede evaluando:
Desarrollando, mediante la propiedad distributiva de la diferenciación y mediante las definiciones para polinomios (una recta lo es, de grado 1) y exponenciales, se obtiene:
de modo tal que operando se obtiene, que el candidato a mínimo debe satisfacer:
La Fig. 6 muestra un modelo de este tipo, con parámetros a=1.5 ARG/km k=0.5 km y C0=100 ARG (En donde ARG representa pesos, sí... un poco desactualizado... pero en 1925 debía ser bastante!)
Fig.6 - Ejemplo del modelo desarrollado.