La PizarraMatemática.

Máster y doctorado

Optimización no diferenciable

Una esquina puede ser un mínimo aunque no tenga derivada. En lugar de forzar una pendiente única, buscaremos todas las rectas que tocan la gráfica y quedan por debajo. Esa idea permitirá certificar mínimos, incorporar restricciones y entender dos familias de algoritmos.

1. Una pendiente de soporte compara con todos los puntos

Para f propia convexa y finita en x, g es subgradiente si f(y)\ge f(x)+\langle g,y-x\rangle\quad\text{para todo }y. El lado derecho es una recta o plano que toca en x y queda debajo globalmente. El conjunto de tales g se llama subdiferencial ∂f(x).

El conjunto ∂f(x) es cerrado y convexo; puede contener varios vectores o estar vacío en la frontera del dominio. Una función convexa finita en un entorno de x sí tiene subgradientes allí. No basta que una recta parezca quedar debajo en dos puntos de una tabla.

Un soporte que se comprueba con un cuadrado

f(y)=y² y punto x=2; proponemos g=4.

  1. Recta

    f(2)+4(y−2)=4y−4.
  2. Diferencia

    y²−(4y−4)=(y−2)².
  3. Global

    El cuadrado es no negativo para cualquier y.

Resultado

4∈∂f(2), con soporte que toca en (2,4).

Lo comprobamos

La comprobación cubre toda la recta, no solo valores próximos.

Un borde sin subgradientes finitos

f(y)=−√y en y≥0, y +∞ fuera; estudiamos x=0.

  1. Exigir

    Un soporte g tendría −√y≥gy para todo y>0.
  2. Dividir

    g≤−1/√y para todo y>0.
  3. Acercarse

    El lado derecho baja sin límite cuando y→0⁺. Ningún g real cumple todas las desigualdades.

Resultado

∂f(0) está vacío aunque f sea propia convexa y semicontinua inferior.

Lo comprobamos

El punto es de frontera, no un punto con entorno de valores finitos.

2. El valor absoluto enseña las tres regiones de una esquina

Si f es convexa y diferenciable en un punto interior, ∂f(x)={∇f(x)}.
Para |x|, en x>0 la pendiente es uno y en x<0 es menos uno.

En cero no hay derivada ordinaria, pero sí un intervalo de soportes.

En una esquina caben muchas rectas por debajo de la función
f(x)=|x|Todas las pendientes g∈[−1,1]|y|≥gy para cualquier y.Una esquina admite varios soportes.
Valor absoluto y cinco rectas representativas de las infinitas pendientes entre menos uno y uno que soportan la gráfica en cero. La desigualdad valor absoluto de y mayor o igual que g por y se cumple globalmente para todas esas pendientes. Las rectas de pendientes uno y menos uno coinciden con una rama, y las interiores solo tocan en cero. Se usan escalas iguales en ambos ejes; no se está afirmando que el valor absoluto tenga derivada ordinaria en cero.
Para hallar ese intervalo pedimos |y|≥gy.
Con y>0 da g≤1; con y<0 da g≥−1.

Cada g entre ambos extremos funciona en las dos regiones. Por eso ∂|·|(0)=[−1,1], sin elegir una pendiente única arbitraria.

Probar una pendiente interior

g=−1/2 en x=0 para el valor absoluto.

  1. y≥0

    |y|−gy=y+y/2=3y/2≥0.
  2. y<0

    |y|−gy=−y+y/2=−y/2>0.
  3. Tocar

    En y=0 ambos lados son cero.
Resultado

−1/2 es un subgradiente válido.

Lo comprobamos

También son válidas pendientes distintas; eso no crea una derivada ordinaria.

Excluir una pendiente demasiado grande

Proponemos g=3/2 en x=0.

  1. Probar

    Basta encontrar un y que viole la desigualdad global.

  2. y=2

    |2|=2 pero gy=3.
  3. Fallo

    Dos no es mayor o igual que tres.

Resultado

3/2 no pertenece al subdiferencial.

Lo comprobamos

Con pendientes menores que −1 se encuentra el fallo tomando y negativo.

3. La pendiente cero certifica un mínimo global convexo

La definición da 0∈∂f(x) si y solo si f(y)≥f(x) para todo y.

Es la condición de Fermat convexa: certifica un mínimo global, incluidas esquinas y restricciones incorporadas al dominio.

Existencia y unicidad siguen siendo preguntas distintas. Subniveles compactos y semicontinuidad inferior, con función propia y valor ínfimo finito, dan existencia. Convexidad estricta da a lo sumo un minimizador, pero no garantiza que se alcance el ínfimo.

Un mínimo en una esquina

f(x)=|x|.

  1. Subdiferencial

    En cero contiene [−1,1].

  2. Cero

    La pendiente cero está incluida.

  3. Certificado

    Para todo y, |y|≥0=f(0).
Resultado

x=0 es mínimo global; además es único porque |y|=0 exige y=0.

Lo comprobamos

La unicidad se ha comprobado directamente, sin afirmar que |x| sea estrictamente convexa.

Estricta convexidad sin minimizador

f(x)=eˣ en toda ℝ.

  1. Valores

    Son siempre positivos.

  2. Ínfimo

    Al tomar x→−∞ se acercan a cero sin alcanzarlo.
  3. Fermat

    La derivada eˣ nunca es cero, coherente con no tener un mínimo.

Resultado

La forma convexa no sustituye una condición de alcanzabilidad.

Lo comprobamos

La segunda derivada positiva no añade un punto en el infinito.

4. Las reglas de suma y máximo necesitan su dominio

Para funciones convexas finitas en todo ℝⁿ, ∂(f+h)=∂f+∂h.

Si h es diferenciable, sumamos ∇h. Con dominios restringidos se requieren cualificaciones, como intersección de interiores relativos; no se presupone la igualdad en cualquier frontera.

Si f=maxᵢfᵢ con un número finito de funciones convexas diferenciables, su subdiferencial es la envolvente convexa de los gradientes activos, aquellos cuyo valor alcanza el máximo en ese punto. Un gradiente de una rama que queda por debajo no se añade automáticamente.

Sumar una esquina y una parte suave

f(x)=|x|+x² en cero.

  1. Primera parte

    ∂|·|(0)=[−1,1].
  2. Suave

    La derivada de x² en cero es cero.

  3. Suma

    ∂f(0)=[−1,1]+{0}=[−1,1].
Resultado

Cero pertenece y certifica el mínimo global.

Lo comprobamos

Las dos funciones son finitas y convexas en toda ℝ.

Identificar ramas activas

f(x)=max{x,−x,−2}.

  1. En x=0

    Las dos primeras ramas valen cero y la tercera −2: solo las dos primeras son activas.

  2. Pendientes

    Son uno y menos uno.

  3. Envolvente

    ∂f(0)=[−1,1].
    En x=2 solo gana la rama x, así ∂f(2)={1}.
Resultado

Las pendientes se eligen por actividad en el punto.

Lo comprobamos

La tercera rama es inactiva; que su pendiente cero caiga dentro del intervalo no la convierte en activa.

5. Una restricción aporta un cono normal

Para C convexo cerrado no vacío, su indicador I_{C} vale cero en C e infinito fuera.

En x∈C, ∂I_{C}(x)=N_{C}(x)={n:⟨n,y−x⟩≤0 para todo y∈C}.

Son direcciones exteriores que la restricción puede equilibrar.

Bajo una regla de suma válida, minimizar f sobre C equivale a 0∈∂f(x)+N_{C}(x). Si f es convexa diferenciable en un entorno de C, equivale a ⟨∇f(x),y−x⟩≥0 para todo y∈C. En una frontera el gradiente puede ser distinto de cero.

Una pendiente positiva en un mínimo restringido

Minimizar f(x)=x en C=[0,∞).

  1. Cono en cero

    N_{C}(0)=(−∞,0], porque ny≤0 para todo y≥0.
  2. Objetivo

    ∇f=1.
  3. Equilibrar

    Elegimos n=−1; 1+(−1)=0.
    El mínimo es x=0.
Resultado

La restricción impide desplazarse hacia valores negativos.

Lo comprobamos

Exigir derivada cero habría perdido el mínimo de frontera.

Proyectar sobre un intervalo

Minimizar (u−3)²/2 sobre C=[0,1].

  1. Candidato

    u=1 es el extremo más próximo a tres.
  2. Gradiente

    u−3=−2.
  3. Cono

    N_{C}(1)=[0,∞) contiene dos, que equilibra −2.
Resultado

El óptimo restringido es uno.

Lo comprobamos

El punto crítico libre tres no es admisible.

6. El umbral suave reúne los candidatos de las tres regiones

Minimizamos Q(x)=λ|x|+(x−a)²/2 con λ>0.
En x>0, 0=λ+x−a da x=a−λ y exige a>λ.
En x<0, 0=−λ+x−a da x=a+λ y exige a<−λ.
En x=0, 0∈λ[−1,1]−a equivale a |a|≤λ.
Así x^*=\operatorname{sign}(a)\max\{|a|-\lambda,0\}.

La parte cuadrática da estricta convexidad y coercividad; el mínimo existe y es único, incluidos los bordes a=±λ.

Una señal negativa fuera del umbral

a=−5 y λ=2.

  1. Región

    a<−λ, así esperamos x<0.
  2. Resolver

    x=a+λ=−3.
  3. Verificar

    −2+(−3)−(−5)=0 y el candidato es negativo.
Resultado

El módulo baja de cinco a tres y el signo se conserva.

Lo comprobamos

No se resta dos al número negativo, porque esa sería la ecuación de otra región.

El punto exacto del umbral

a=2 y λ=2.

  1. Cero

    a pertenece a λ[−1,1]=[−2,2].
  2. Optimalidad

    En x=0 elegimos subgradiente s=1: λs−a=0.
  3. Unicidad

    La estricta convexidad del objetivo impide otro minimizador.

Resultado

El resultado es cero también en el borde.

Lo comprobamos

No se exige |a|<λ estrictamente.

7. Un paso de subgradiente no garantiza descenso

La iteración xₖ₊₁=xₖ−αₖgₖ, con gₖ∈∂f(xₖ) y αₖ>0, usa una pendiente de soporte. A diferencia del gradiente suave con paso adecuado, la dirección −gₖ no asegura que f baje para toda elección admisible.

Un paso de subgradiente puede salir de un mínimo
x₀=0; g=1; α=3/4−3/40f(x₀)=0 → f(x₁)=3/4El mejor valor visitado sigue en cero.Subgradiente válido ≠ descenso seguro.
Función valor absoluto, punto inicial cero, subgradiente elegido g uno y paso alpha tres cuartos. La iteración da menos tres cuartos, cuyo valor es tres cuartos: aumenta desde el mínimo cero. El tramo horizontal rojo muestra el cambio de coordenada y el vertical punteado lee el valor de la función; no representan el recorrido de una dinámica continua sobre la gráfica. Elegir g cero habría dejado el mínimo fijo, de modo que este ejemplo muestra ausencia de garantía de descenso para cualquier elección admisible, no que todo paso de subgradiente aumente.

Si se imponen restricciones convexas cerradas, proyectamos el nuevo punto sobre ellas. Para garantías hacen falta mínimo existente, un control de los subgradientes y una regla de pasos apropiada. Se controla el mejor valor visitado o ciertos promedios; no se promete una bajada monótona de cada valor.

Moverse desde un mínimo en una dirección permitida

f=|x|, x₀=0, g₀=1 y α₀=3/4.

  1. Validez

    Uno pertenece al intervalo subdiferencial [−1,1].

  2. Paso

    x₁=0−3/4=−3/4.
  3. Valor

    f(x₁)=3/4>f(x₀)=0.
Resultado

Una pendiente válida puede aumentar el objetivo.

Lo comprobamos

Elegir g₀=0 habría mantenido el mínimo; el fallo es la garantía para cualquier elección.

Un paso constante puede oscilar

f=|x|, x₀=2/5 y α=1, usando el signo en puntos no nulos.

  1. Primero

    x₁=2/5−1=−3/5.
  2. Segundo

    x₂=−3/5+1=2/5.
  3. Repetir

    La órbita alterna y los valores son 2/5 y 3/5, sin alcanzar cero.

Resultado

La elección de pasos importa incluso en una función muy sencilla.

Lo comprobamos

No se infiere convergencia de que cada g sea un subgradiente.

8. La cota de convergencia sale de expandir una distancia

Sea x* un mínimo, ||x₀−x*||≤R y ||gₖ||≤G.

Expandimos ||xₖ−αₖgₖ−x*||².

El soporte da ⟨gₖ,xₖ−x*⟩≥f(xₖ)−f*.

Si hay proyección sobre un convexo que contiene x*, su no expansión no empeora la cota de distancia.

Sumar y descartar la distancia final no negativa produce \min_{0\le k<N}f(x_k)-f^*\le\frac{R^2+G^2\sum_{k=0}^{N-1}\alpha_k^2}{2\sum_{k=0}^{N-1}\alpha_k}.

Si Σαₖ=∞ y Σαₖ²<∞, el lado derecho tiende a cero.

La garantía se refiere al mejor valor, no al descenso de cada iteración.

Elegir un paso para un horizonte fijo

N pasos, R>0,G>0 y α=R/(G√N).

  1. Suma

    Σα=Nα=R√N/G.
  2. Cuadrados

    G²Σα²=R².
  3. Cota

    (R²+R²)/(2R√N/G)=RG/√N.
Resultado

Es una garantía de valor con ese horizonte y esas constantes.

Lo comprobamos

Requiere conocer o acotar R y G; no prueba una tasa sin hipótesis.

Una regla que satisface las dos series

αₖ=1/(k+1).

  1. Primera

    Σ1/(k+1) diverge.

  2. Segunda

    Σ1/(k+1)² converge.

  3. Resultado

    Con las demás hipótesis, la cota del mejor valor tiende a cero.

Resultado

Una regla decreciente puede compensar las oscilaciones.

Lo comprobamos

Que α tienda a cero por sí solo no basta; pasos con suma finita pueden detenerse antes.

9. El proximal resuelve un problema en vez de extrapolar una pendiente

Para f propia convexa semicontinua inferior en ℝⁿ y λ>0, \operatorname{prox}_{\lambda f}(a)=\arg\min_u\{f(u)+\|u-a\|^2/(2\lambda)\} existe y es único. El término cuadrático y un minorante afín de f dan coercividad; la estricta convexidad da unicidad.

Su condición es 0∈∂f(u)+(u−a)/λ.

Con f=|·| da el umbral suave; con f=I_{C} da proyección euclídea sobre el convexo cerrado no vacío C. El parámetro λ multiplica f si escribimos de forma equivalente λf(u)+||u−a||²/2.

Una proyección como proximal

f=I_[0,2], a=3 y cualquier λ>0.

  1. Dominio

    Fuera de [0,2] el coste es infinito.

  2. Dentro

    Solo minimizamos la distancia cuadrada a tres.

  3. Óptimo

    El punto más próximo es dos, independiente de λ en este indicador.

Resultado

prox de una restricción es su proyección.

Lo comprobamos

Se necesita un conjunto no vacío cerrado convexo para esta afirmación de existencia y unicidad.

Ver cómo cambia el parámetro en una cuadrática

f(u)=u²/2 y a=6.

  1. Condición

    u+(u−6)/λ=0.
  2. Despejar

    u=6/(1+λ).
  3. Casos

    λ=1 da tres y λ=2 da dos.
Resultado

La penalización se hace más fuerte al aumentar λ en esta convención.

Lo comprobamos

No se confunde λf con la función f evaluada en λu.

10. El gradiente proximal conserva la estructura del objetivo

Para F=h+f, con h convexa diferenciable y gradiente L-Lipschitz, f propia convexa semicontinua inferior y mínimo existente, el paso es x_{k+1}=\operatorname{prox}_{\alpha f}(x_k-\alpha\nabla h(x_k)). Una elección 0<α≤1/L, con L>0, da el marco usual de descenso y convergencia de valores; no se traslada sin revisar hipótesis.

En problemas no convexos se usan otras subdiferenciales, como Clarke, y distintas condiciones estacionarias. La presencia del cero ya no certifica por sí sola un mínimo global. El certificado global que hemos usado viene de convexidad y soporte, no de la palabra subgradiente.

Un paso proximal que resuelve este ejemplo

h(x)=(x−3)²/2 y f(x)=|x|; L=1 y α=1.

  1. Gradiente

    h′(x)=x−3.
  2. Punto intermedio

    x−h′(x)=3, para cualquier x inicial.
  3. Proximal

    soft(3,1)=2, el minimizador de F.
Resultado

Aquí se alcanza el óptimo en un paso por la estructura especial de la cuadrática.

Lo comprobamos

No se promete un paso para una suma suave general.

Un cero estacionario no convexo puede ser un máximo

F(x)=−x², suave pero no convexa.

  1. Derivada

    F′(0)=0.
  2. Comparación

    Para x≠0, F(x)<F(0)=0.
  3. Clasificar

    Cero es un máximo global y no un mínimo.

Resultado

Una generalización no convexa requiere otra conclusión.

Lo comprobamos

La desigualdad de soporte convexo falla para esta función.

Ahora lo resolvemos de principio a fin.

EJEMPLO RESUELTO 1

Subgradiente en una esquina

Demuestra que g=1/2 es subgradiente de |x| en 0.

  1. Definir

    En x=0 pedimos |y|≥(1/2)y para todo y.
  2. y≥0

    La diferencia es y−y/2=y/2≥0.
  3. y<0

    La diferencia es −y−y/2=−3y/2>0.
  4. Contacto

    En y=0 ambas expresiones valen cero.
  5. Concluir

    1/2∈∂|·|(0).
  6. Comprobar

    La pendiente queda dentro del intervalo completo [−1,1].

Resultado

1/2∈∂|·|(0).

Lo comprobamos

La diferencia|y|−y/2 vale y/2 si y≥0 y−3y/2 si y<0, siempre no negativa.

EJEMPLO RESUELTO 2

Mínimo con penalización

Minimiza |x|+(x−3)²/2.

  1. Objetivo

    Q(x)=|x|+(x−3)²/2 es estrictamente convexo y coercivo.
  2. Región positiva

    La derivada allí es 1+x−3=x−2, con candidato x=2>0.
  3. Región negativa

    La derivada x−4=0 daría cuatro, que no es negativo.
  4. En cero

    El subdiferencial es [−1,1]−3=[−4,−2], sin cero.
  5. Certificado

    En x=2 la condición cero del subdiferencial sí se cumple: es mínimo global.

  6. Unicidad

    La estricta convexidad lo hace único; Q(2)=1/2+2=5/2.
Resultado

El mínimo está en x=2.

Lo comprobamos

Las tres regiones se han revisado y la condición de Fermat certifica el candidato; la convexidad estricta lo hace único.

Los símbolos que hemos usado

∂f(x)
Conjunto de subgradientes convexos de f en un punto de valor finito.
g
Vector que da una pendiente de soporte global.
I_C
Indicador convexo: cero en C y+∞ fuera.
N_C(x)
Cono normal: vectores v con⟨v,y−x⟩≤0 para todo y∈C.
αₖ
Paso positivo de una iteración de subgradiente.
prox_λf(a)
Minimizador de f(u)+‖u−a‖²/(2λ), conλ>0.
sign a
Signo de a; vale−1,0 o1.

Ahora te toca a ti.

Tienes cinco ejercicios. Inténtalos; las pistas te ayudan si te atascas.

Practicar este tema →

Para ampliar