Mundo exterior - Matemáticas: Tema II (derivadas segundas)

Buenas!

Llevaba algún tiempo sin escribir... y es que vuelvo a andar liado y hasta arriba!

La parte positiva es que me siento inusualmente activo... jejeje. Por otra parte, queda poco para acabar este tema matemático y volver a saco con COBOL.

En mis ratos libres hago mis pequeñas investigaciones con COBOL para adaptar los requerimientos de la app... también he tenido (recientemente) una revelación: ¿se podría, de alguna manera, tratar un sudoku como una matriz - o conjunto de ellas? en ello estoy, también.

En fin, dejémonos de divagaciones y vamos a las matemáticas.

Esta entrada tratará de segundas derivadas.

Definimos las derivadas segundas como una función lineal que transforma el espacio de R2 en R.

Podríamos aventurarnos a lanzar una primera ecuación correspondiente a una segunda derivada en una función de varias variables, a saber:

f(x,y) ≈ f(a,b) + f_x(a,b)(x-a) + f_y(a,b)(y-b) + ½ [f_xx(a,b)(x-a)² + 2f_xy(a,b)(y-b) + f_yy(a,b)(y-b)²]

f(x,y) ≈ f(a,b) + grad(f(a,b)) * R1 + ½ * R2 * H * R1

donde (por ejemplo) fxx es la derivada parcial de "x" sobre la derivada parcial de "x"; es decir, se aplican varias derivadas parciales sucesivas, tomando como variable la variable indicada.

Por otra parte, R1, R2, y H son matrices, tales que:

R1 = x-a

         y-b

R2 = (x-a   y-b)

H = f_xx    f_xy

       f_yx   f_yy

A "H" se le conoce como matriz Hessiana, y se suele denotar como "H(a,b)".

En el caso que nos atañe (funciones de dos variables), "H" es una matriz 2x2, y se dice que tiene dos valores propios (vaps): fxx y fyy (ya que fxy = fyx).

Veamos qué pasa con estos vaps:

• si det(H) = 0: alguno de los vaps. es 0
• si det(H) > 0: los dos vaps son mayores que 0, por tanto la función está por encima del plano "z"
• si (a+b) > 0: H(a,b) es simétrica, cuadrada, y se puede diagonalizar.
• si det(H) < 0: los dos vaps son menores que 0, por tanto la función está por debajo del plano "z"
• si (a+b) > 0: H(a,b) es simétrica, cuadrada y se puede diagonalizar.

En caso que sea 0, no nos reporta ninguna información útil (al parecer).

Veamos qué son los vaps. de H:

vaps(H(a,b)) = l₁, l₂

l_{1 *} l₂  = det(H) = f_xx*f_yy-f_xy2

l_{1 +} l₂  = f_{xx +} f_yy

veamos un ejemplo en el que se aplique todos estos conceptos:

f(x,y) = x² + y/3 + xy²

Derivada segunda en el punto (7,2)

f(7,2) = 49 + 2/3 + 7*4 = 233/3

f_x = 2x + y² => en (7,2) = 18

f_y = 1/3 + 2xy => en (7,2) = 85/3

f_xx = 2

f_xy = f_yx = 2x

f_yy = -2y

aplicamos fórmulas:

z = 233/3 + 18(x-7) + 85/3(y-2)

H_(a,b) = 2    2y

            2y  2x

H_(7,2) = 2   4

            4  14

det(H_(7,2)) = 12

2 + 14 = 16

_____________

vaps > 0

Plano tangente:

Z = z + ½(2(x-7)² + 2*4(x-7)(y-2) + 14(y-2)²)

Visto el ejemplo, pasemos a un inciso en las derivadas segundas: extremos relativos.

Para obtener los extremos relativos de una función de dos variables debe cumplirse lo siguiente:

fx = 0

fy = 0

Entonces, usando las fórmulas anteriores, obtendremos varios "juegos" de puntos. Basándose en los vaps de H(a,b), podemos saber qué tipo de extremo relativo estamos tratando:

• si los dos vaps > 0: mínimo
• si los dos vaps < 0: máximo
• si un vap > 0 y otro vap < 0: intersección (hacia un eje es mínimo, hacia otro es máximo)
• si algún vap = 0: no aporta información

Para obtener esto puntos, debemos ejecutar lo siguiente:

 - Obtenemos "fx" y "fy"

 - Obtenemos "fxx", "fxy", "fyy"

 - Igualamos "fx" y "fy" a 0, creando un sistema de ecuaciones

 - Resolvemos el sistema de ecuaciones

 - Para cada juego de soluciones, aplicamos H(a,b)

Y bien, con esto acabamos la cuarta entrega del segundo tema... las dos próximas serán cortas (ejemplos y utilidades).

Os dejo, y como siempre, ¡Hasta la próxima!