3. Clasificación de puntos críticos

De manera análoga al caso de una variable, usamos el polinomio de Taylor de segundo orden (como una primera opción) para clasificar puntos críticos.

La fórmula de Taylor de segundo orden en dos variables, alrededor de (a,b) se puede definir en un entorno U abierto alrededor de (a,b), si f : U 2, es de clase C2. En ese caso,

T2(x,y) = f(a,b) + fx(a,b)(x a) + fy(a,b)(y b) + 1 2fxx(a,b)(x a)2 + f xy(a,b)(x a)(y b)

Haciendo x = (x,y),p = (a,b), (h,k) = (x a,y b), A = fxx(a,b),B = fxy(a,b) y C = fyy(a,b) se tiene

f(x) = f(p) + fx(p)h + fy(p)k + 1 2 [fxx(p)h2 + 2f xy(p)hk + fyy(p)k2] = f(p) + f(p) (h,k) + 1 2 [Ah2 + 2Bhk + Ck2] + r(x,p)

donde el resto r(x,p) disminuye más rápido que ||x p||2, es decir, r(x,p) ||x p||2 0 si x p

Ahora, si p es un punto crítico de f, entonces

f(x)f(p) = 1 2 [Ah2 + 2Bhk + Ck2]+r(x,p) con  { A = fxx(p) B = fxy(p)

Para determinar si en el el punto crítico p, la función f alcanza un máximo o mínimo local, debemos determinar el signo de f(x) f(p) en un entorno de p, para saber si f(x) f(p) o f(x) f(p). Si en cualquier entorno de p hay puntos donde f cambia de signo, entonces tenemos un punto de silla.

La teoría es similar al caso de una variable: En presencia de extremos locales, el signo de f(x) f(p) es el signo de Ah2 + 2Bhk + Ck2 en un entorno suficientemente pequeño de cada punto crítico donde f alcanza máximos o mínimos locales1 . Si el determinante D2 = AC B2 > 0, entonces la forma cuadrática Ah2 + 2Bhk + Ck2, es siempre positiva o siempre negativa. Si D2 = AC B < 0 la forma cuadrática cambia de signo.

Podemos analizar el signo de Ah2 + 2Bhk + Ck2 completando cuadrados. Sea D2 = AC B2, entonces

A(Ah2 + 2Bhk + Ck2) = (Ah + Bk)2 + (AC B2)k2A(Ah2 + 2Bhk + Ck2) 0 si D 2 > 0

Entonces

Ah2 + 2Bhk + Ck2 0  si  D 2 > 0 y A > 0 Ah2 + 2Bhk + Ck2 0  si  D 2 > 0 y A < 0

con esto se puede establecer que

1.
f(x) f(p) 0 en un entorno de p si D2 > 0 y si A > 0
2.
f(x) f(p) 0 en un entorno de p si D2 > 0 y si A < 0

Si D2 < 0 entonces, podemos razonar con varios casos. Por simplicidad solo consideremos dos casos con A0 y C0. Si B = 0 entonces A y C tienen signos contrarios, por lo que la forma cuadrática cambia de signo sobre las rectas x = 0 y y = 0. Si B0, entonces la forma cuadrática cambia de signo sobre las rectas y = 0 y By = Ax. Entonces se podría establecer, comparando f(x) f(p) 0 con la forma cuadrática, que p es un punto de silla.

Si D2 = 0 entonces en general, si A0 y B0, hay rectas que pasan por el origen, sobre las que el término (Ah + Bk)2 se anula, entonces f(x) f(p) = r(x,p) sobre estas rectas, es decir, no podemos decir algo acerca del signo de f(x) f(p).

PIC

Figura 5.6: Forma cuadrática Ah2 + 2Bhk + Ck2 para distintos valores de D2 y A

Específicamente: Si f tiene segundas derivadas continuas y si el determinante D2 = AC B2 es positivo, entonces el signo de f(x) f(p) es el signo de Ah2 + 2Bhk + Ck2, en un entorno suficientemente pequeño de p. Si el determinante H es negativo, entonces f(x) f(p) cambia de signo con la forma cuadrática, en trayectorias contenidas en un entorno de p .

En vista de la anterior afirmación, la clasificación de los puntos críticos depende del signo de la forma cuadrática (si D20). El signo de la forma cuadrática es fácil de establecer usando el discriminante D2. La idea geométrica es la que se muestra a continuación.

Teorema 13 — (Condición suficiente)..

Sea f : U 2 de clase C2 en un subconjunto abierto U de 2. Consideremos el “discriminante” D2(x,y) = fxx(x,y) fyy(x,y) [fxy(x,y)]2. Si (x0,y0) U es punto crítico de f, entonces

1.
si D2(x0,y0) > 0 y fxx(x0,y0) > 0, entonces f alcanza un mínimo local en (x0,y0).
2.
si D2(x0,y0) > 0 y fxx(x0,y0) < 0, entonces f alcanza un máximo local en (x0,y0).
3.
Si D2(x0,y0) < 0, entonces (x0,y0,f(x0,y0)) es un punto de silla.

PIC

Figura 5.7: Polinomio de Taylor P2(x,y) y la forma cuadrática, sobre la superficie S

El teorema solo da condiciones suficientes: No nos dice algo si D2(x0,y0) = 0 pues, como ya vimos en este caso, en un entorno de p,

f(x) f(p) = r(x,p)

por lo que no podemos determinar el signo de f(x) f(p) (ver el ejercicio ). En este caso se podría usar otros métodos para clasificar.

En el teorema se puede usar fyy en vez de fxx pues si D2(x0,y0) > 0, ambas tienen el mismo signo.

Ejemplo116(Polinomio de Taylor y clasisificación de puntos críticos).

Considere la superficie S : z = xex2y2+1 + 1.

Puntos críticos.
{ fx = 0 ex2y2+1 (1 2x2) = 0x = ± 1 2(E1)

Una solución (x0,y0) del sistema requiere que se anule la ecuación E1 y la ecuación E2. La ecuación E1 solo tiene como soluciones posibles x = ± 1 2. Por lo tanto, en la E2 tenemos que descartar x = 0. La solución dels sistema es

( 1 2,0) ( 1 2,0)

El polinomio de Taylor en (a,b) es

T2(x,y) = f(a,b) + fx(a,b)(x a) + fy(a,b)(y b) + 1 2fxx(a,b)(x a)2 + f xy(a,b)(x a)(y b)

El polinomio de Taylor en un punto crítico Pc = (a,b) es el que sigue, la expresión en magenta es la forma cuadrática que determina la clasificación del punto crítico.

T2(x,y) = f(a,b) + 1 2 [fxx(a,b)(x a)2 + 2f xy(a,b)(x a)(y b) + fyy(a,b)(y b)2]

donde { fxx(x,y) = ex2y2+1 (4x3 6x) fxy(x,y) = ex2y2+1 (4x2y 2y) fyy(x,y) = ex2y2+1 (4xy2 2x)

En el código que sigue, en Mathematica, f[x,y] es la función y T[x,y] es cada polinomio de Taylor. Se realiza la representación gráfica en un entorno adecuado de cada punto crítico.

dx=0.5; dy=0.5; 
Plot3D[f[x,y],{x,-2,2},{y, -1,1}, Mesh->None], 
 
T[x_,y_]=1+Sqrt[E/2]+1/2 (-2 Sqrt[2 E] (-(1/Sqrt[2])+x)^2-Sqrt[2 E] y^2); 
Plot3D[T[x,y],{x,0.7-dx,0.70+dx},{y, 0-dy,0+dy}, Mesh->None] 
 
T[x_,y_]=1-Sqrt[E/2]+1/2 (2 Sqrt[2 E] (1/Sqrt[2]+x)^2+Sqrt[2 E] y^2); 
Plot3D[T[x,y],{x,0.7-dx,0.70+dx},{y, 0-dy,0+dy}, Mesh->None]

En el código que sigue, en Mathematica, F_i[x,y] son las formas cuadráticas:

F_1[x_,y_]=1/2 (-2 Sqrt[2 E] (-(1/Sqrt[2])+x)^2-Sqrt[2 E] y^2); 
 
F_2[x_,y_]=1/2 (2 Sqrt[2 E] (1/Sqrt[2]+x)^2+Sqrt[2 E] y^2);

PIC

PIC

Figura 5.8: Polinomio de Taylor y las formas cuadráticas, con un solo signo!
Clasificación de los puntos críticos usando el teorema

D1 = fxx = ex2y2+1 (4x3 6x)

D2 = e2x22y2+2 | 4x3 6x00 4x2y 2y 4x2y 2y00 4xy2 2x | = e2x22y2+2 (12x28x2y24y28x4)

Pc = ( 1 2,0) D2(Pc) = 4e > 0 fxx(Pc) = 22e < 0 ( 1 2,0,1 + e 2 )  máximo local  Pc = ( 1 2,0) D2Pc(Pc) = 4e > 0 fxx(Pc) = 22e > 0 ( 1 2,0,1 e 2 )  mínimo local 

Ejemplo117

Calcule y clasifique los puntos críticos de la función f(x,y) = x3 + 3y y3 3x.

Solución. :

Puntos críticos: { ∂f ∂x = 0 ∂f ∂y = 0 { 3x2 3 = 0x = ±1 3 3y2 = 0y = ±1

En este caso, cualquier combinación de signos anula el sistema, por eso los puntos críticos son (1,1),(1,1),(1,1) y (1,1).

Clasificación. D2(x,y) = (6x)(6y) (0)2

(x0,y0) D1 = fxx(x0,y0) H = D2(x0,y0)  Clasificación 
(1,1) -16 36 (1,1,0) es punto de silla 
(1,1) -16 -136 (1,1,1.2) es mínimo local. 
(1,1) 6 -1 36 (1,1,1.2) es máximo local. 
(1,1) 6 36 (1,1,0) es punto de silla 

La representación gráfica de f se muestra en al figura.

PIC

Figura 5.9: f(x,y) = x3 + 3y y3 3x
Ejemplo118

Calcule y clasifique los puntos críticos de la función f(x,y) = x4 + y4 2x2 + 4xy 2y2.

Solución. :

Puntos críticos: { ∂f ∂x = 0 ∂f ∂y = 0 { 4x3 + 4y 4x = 0 4y3 4y + 4x = 0 (E2) 

Sumando miembro a miembro obtenemos x3 + y3 = 0x = y.

Ahora sustituimos en la ecuación (E2), queda 4x3 4x 4x = 0x(x2 2) = 0; con lo cual obtenemos los puntos críticos

(0,0),(2, 2),(2,2)


Clasificación. D2(x,y) = (12x2 4)(12y2 4) (4)2

(x0,y0) D1 = fxx(x0,y0) H = D2(x0,y0)  Clasificación 
(0,0) 4 0  Criterio no decide. 
(2, 2) 20 384 (2, 2, 8) es mínimo local. 
(2,2) 20 384 (2,2, 8) es mínimo local. 

La representación gráfica de f se muestra en al figura. Aunque D2(0,0) = 0 y el criterio no proporciona información, la gráfica a la derecha nos indica que se trata de un punto de silla.

PIC

Figura 5.10: f(x,y) = x4 + y4 2x2 + 4xy 2y2
Ejemplo119

Sea f(x,y) = 6xy 2x2y 3xy2. Calcule y clasifique los puntos críticos de f.

Solución. : Los puntos críticos se obtienen resolviendo el sistema f = (0,0),

{ ∂f ∂x = 0 ∂f ∂y = 0 { 6y 4xy 3y2 = 0 6x 2x2 6xy = 0 { y(6 4x 3y) = 0y = 0 o y = 6 4x 3 6x 2x2 6xy = 0 (E2) 

Si y = 0, al sustituir en la ecuación (E2) obtenemos los puntos críticos (0,0),(3,0).

Si y = 6 4x 3 , al sustituir en la ecuación (E2) obtenemos los puntos críticos (0,2), (1, 2 3 ).

Clasificación. D2(x,y) = (4y)(6x) [6 4x 6y]2

(x0,y0) D1 = fxx(x0,y0) H = D2(x0,y0)  Clasificación 
(0,0) 0 36 (0,0,0) es punto de silla 
(3,0) 0 36 (3,0,0) es punto de silla 
(0,2) 8 36 (0,2,0) es punto de silla 
(1,23) 83 12 (1,23,43) es máximo local. 

PIC

Figura 5.11: f(x,y) = 6xy 2x2y 3xy2

Ejemplo120

Sea z = yxex + y2. Calcule y clasifique los puntos críticos de z.

Solución. : Los puntos críticos se obtienen resolviendo el sistema z = (0,0),

{ ∂z ∂x = 0 ∂z ∂y = 0 { exy exxy = 0 exx + 2y = 0 { y(1 x) = 0y = 0 o x = 1 exx + 2y = 0 (E2) 

I caso. Si y = 0, sustituimos en (E2) y obtenemos x = 0.

II caso Si x = 1 sustituimos en (E2) y obtenemos y = 1 2e.

Clasificación. D2(x,y) = 2yex (x 2) (ex exx)2

(x0,y0) D1 = fxx(x0,y0) H = D2(x0,y0)  Clasificación 
(0,0) 0 1 < 0 (0,0,0) es punto de silla 
(1, 1 2e) 1 2e2 > 0 1 e2 > 0 (1,12e, 14e2) es mínimo local 
Ejemplo121

Sea z = x2y2 x y. Calcule y clasifique los puntos críticos de z.

Solución. : Los puntos críticos se obtienen resolviendo el sistema z = (0,0),

{ ∂z ∂x = 0 ∂z ∂y = 0 { 2xy2 1 = 0 (E1)  2x2y 1 = 0 (E2) 

Como y = 0 no es solución, podemos despejar x = 1 2y2 de (E1). Ahora sustituimos en (E2) y obtenemos y = 1 23. Entonces tenemos el punto crítico 2 (1 23,1 23).

Clasificación. D2(x,y) = 2y2 2x2 (4xy)2

(x0,y0) D1 = fxx(x0,y0) H = D2(x0,y0)  Clasificación 
(1 23,1 23) 1 23 343 < 0 (1 23,1 23, 3 223) es punto de silla 
Ejemplo122
Calcule el volumen de la caja rectangular más grande que esté en el primer octante con tres de sus caras en los planos coordenados y un vértice en el plano x + 2y + 3z = 6.
Solución. : Debemos maximizar V = xyz. Como z = 2 x3 2y3, el volumen de la caja es V = xy (2 x3 2y3).

PIC

Puntos críticos. Nos interesa solo x > 0 y y > 0. Entonces,

{ V x = 0 V y = 0 { 2y 3 (3 + x + y) = 0 x 3(6 + x + 4y) = 0 { 3 + x + y = 0 6 + x + 4y = 0 x = 2,y = 1.

Clasificación. D2(x,y) = V xxV yy V xy2 = 2y 3 4x 3 [2 3(x + 2y 3)]2. Así D2(2,1) = 43 > 0 y D1 = V xx(2,1) = 23 < 0. Esto nos dice que el volumen es máximo cuando las dimensiones de la caja son x = 2,y = 1 y z = 2 3. Por otro lado, el volumen máximo es 4 3ul3.

Ejemplo123
Calcule y clasifique los puntos críticos de la función f(x,y) = x2.
Solución. : Primero calculamos los puntos críticos
{ fx = 2x = 0 fy = 0 = 0
PIC
El sistema tiene infinitas soluciones de la forma (0,y). Así que tenemos un número infinito de puntos críticos. D2(x,y) = (2)(0) (0)2 = 0 así que este criterio no da información aunque, de acuerdo a la gráfica , se trata de puntos donde f alcanza mínimos locales.