Página web del curso

Contenidos

Ejercicio

Mínimo local y mínimo global. Cálculo simbólico

Dado un conjunto de puntos $(x_i,y_i)$ con $i = 1,2,\ldots,n$ encontrar el punto central $(m_x,m_y)$ de forma que la suma de las distancias cuadráticas

$$D = \sum\limits_{i=1}^n\left[(x_i-m_x)^2+(y_i-m_y)^2\right]$$

a dicho punto sea mínima. El conjunto de puntos es

$$ \begin{array}{|c|rrrrr|} \hline x & 1 & 3 & 4 & 5\\ \hline y & 2 & -1 & 3 & 1\\ \hline \end{array} $$

La función a minimizar es

$$ \begin{array}{lll} D(m_x,m_y)& = & d_1^2 +d_2^2 +d_3^2 +d_4^2 \\ D(m_x,m_y)& = &[(1-m_x)^2+(2-m_y)^2]+[(3-m_x)^2+(-1-m_y)^2]+\\ & &[(4-m_x)^2+(3-m_y)^2]+[(5-m_x)^2+(1-m_y)^2] \end{array} $$

La condición necesaria de óptimo es que las derivadas parciales sean cero

$$ \begin{array}{lll}\\ \dfrac{\partial D}{\partial m_x}=0 \\ \dfrac{\partial D}{\partial m_y}=0 \end{array} $$

Derivando $$ \begin{array}{rcl}\\ -2(1-m_x)-2(3-m_x)-2(4-m_x)-2(5-m_x)& = & 0\\ -2(2-m_y)-2(-1-m_y)-2(3-m_y)-2(1-m_y)& = & 0 \end{array} $$ Quitamos el factor $-2$ $$ \begin{array}{rcl}\\ (1-m_x)+(3-m_x)+(4-m_x)+(5-m_x)& = & 0\\ (2-m_y)+(-1-m_y)+(3-m_y)+(1-m_y)& = & 0 \end{array} $$

Reorganizamos las ecuaciones $$ \begin{array}{lll}\\ (1+3+4+5)& = & 4 m_x\\ (2-1+3+1)& = & 4 m_y \end{array} $$ La solución es $$ \begin{array}{lll}\\ m_x = \dfrac{1+3+4+5}{4}=3.25\\ m_y = \dfrac{2-1+3+1}{4}=1.25 \end{array} $$

La condición suficiente es que la matriz Hessiana en el punto $\left(3.25,1.25\right)$ sea definida positiva. La matriz Hessiana para la función $D$ en cualquier punto $\left(m_{x},m_{y}\right)$ es

$$ H(m_{x},m_{y})=\left(\begin{array}{cc} \dfrac{\partial^2 D}{\partial m_x^2} & \dfrac{\partial^2 D}{\partial m_x\partial m_y}\\ \dfrac{\partial^2 D}{\partial m_y\partial m_x} & \dfrac{\partial^2 D}{\partial m_y^2} \end{array}\right)= \left(\begin{array}{cc} 2+2+2+2 & 0 \\ 0 & 2+2+2+2 \end{array}\right) = \left(\begin{array}{cc} 8 & 0 \\ 0 & 8 \end{array}\right) $$

Y la matriz Hessiana es constante y la misma para todos los puntos del plano. En particular, es también el valor para el puntos $(3.25,1.25).$

Una condición para que una matriz sea definida positiva es que los determinantes de los menores principales han de ser estrictamente positivos. Es decir

$$ \left|\begin{array}{cc} 8 & 0 \\ 0 & 8 \end{array}\right|=64\gt0\quad\mathrm{y} \qquad \left|8\right|=8\gt0. $$

Y se cumple la condición suficiente de mínimo.

En el dibujo siguiente vemos los valores que la función $D$ toma para diferentes puntos y que, efectivamente, el punto $(3.25,1.25)$ (en rojo) es el punto donde la suma de distancias a los otro cuatro puntos es mínima.