Suavizado · Artículo 03·03

La max-relativa suavizada

Con una métrica en la mano, ya podemos definir la entropía max-relativa suavizada: el mínimo de la max-relativa sobre una bola de estados cercanos. Es monótona en ε, robusta, y finita incluso donde la versión cruda se dispara.

Reunimos las dos piezas de los artículos anteriores —la idea de tolerar un error ε\varepsilon y una forma seria de medir cercanía— y las cristalizamos en una única definición. La entropía max-relativa suavizada es lo que uno esperaría: la mejor DmaxD_{\max} alcanzable si se permite mover el estado dentro de una pequeña bola.

La definición

Sea Bε(ρ)\mathcal{B}^{\varepsilon}(\rho) la bola de suavizado bola de estados a distancia purificada como mucho ε\varepsilon de ρ\rho. La entropía max-relativa suavizada entropía max-relativa suavizada se define como

Dmaxε(ρσ)=minρ~Bε(ρ)Dmax(ρ~σ).D_{\max}^{\varepsilon}(\rho\,\|\,\sigma) = \min_{\tilde\rho \,\in\, \mathcal{B}^{\varepsilon}(\rho)} D_{\max}(\tilde\rho\,\|\,\sigma).

Es un mínimo, no un promedio: entre todos los estados indistinguibles de ρ\rho dentro de la tolerancia, nos quedamos con el que da la DmaxD_{\max} más pequeña. Se admiten estados ligeramente subnormalizados (con Trρ~1\operatorname{Tr}\tilde\rho \le 1), un tecnicismo que engrasa las demostraciones y no cambia la intuición.

Monótona, robusta y finita

De la definición salen tres propiedades sin esfuerzo. Es monótona en ε\varepsilon: una bola más grande contiene más candidatos, de modo que el mínimo solo puede bajar, y por tanto εε\varepsilon' \ge \varepsilon implica DmaxεDmaxεD_{\max}^{\varepsilon'} \le D_{\max}^{\varepsilon}. En particular, siempre DmaxεDmax0=DmaxD_{\max}^{\varepsilon} \le D_{\max}^{0} = D_{\max}: suavizar nunca empeora. Es robusta, porque optimizar sobre una bola borra los saltos que tanto molestaban. Y es finita en cuanto ε\varepsilon baste para escapar de un mal soporte, aun cuando DmaxD_{\max} valiera ++\infty.

La curva de DmaxεD_{\max}^{\varepsilon} frente a ε\varepsilon cuenta toda la historia: arranca en el valor crudo y desciende, suave y sin saltos, a medida que gastamos presupuesto de error. Muévela:

Suavizar convierte el peor caso en algo robusto

Sin suavizar Dmax = 0.000 bits  ·  con ε = 0.15: Dmaxε = 0.000 bits

Permitir un error ε (mover p hacia q) baja el peor cociente. En ε = |p₀ − q₀| llega a 0.

p(0) 0.50
q(0) 0.12
ε (error) 0.15

El caso clásico, para fijar ideas

Con dos distribuciones pp y qq, suavizar consiste en desplazar pp hacia qq tanto como el presupuesto permita, porque acercarse a qq reduce el peor cociente pi/qip_i/q_i. Si pudiéramos gastar ε\varepsilon ilimitado llegaríamos a p~=q\tilde p = q y a Dmax=0D_{\max} = 0; con presupuesto finito nos quedamos a medio camino. Ese «desplazar hacia qq» es, ni más ni menos, una minimización — la idea que el próximo artículo lleva hasta sus últimas consecuencias.

Ejemplos resueltos

Ejemplo resuelto 1 · Suavizar una rebaja concreta

Problema. Con p=(0.8, 0.2)p = (0.8,\ 0.2) y q=(0.2, 0.8)q = (0.2,\ 0.8), calcula DmaxD_{\max} y DmaxεD_{\max}^{\varepsilon} para ε=0.2\varepsilon = 0.2 (distancia de traza, para simplificar).

Sin suavizar. Dmax=log2max ⁣(0.80.2, 0.20.8)=log24=2D_{\max} = \log_2 \max\!\big(\tfrac{0.8}{0.2},\ \tfrac{0.2}{0.8}\big) = \log_2 4 = 2 bits.

Con ε=0.2\varepsilon = 0.2. Bajamos p0p_0 hacia q0=0.2q_0 = 0.2 en 0.20.2: de 0.80.8 a 0.60.6. Entonces

Dmax0.2=log2max ⁣(0.60.2, 0.40.8)=log231.585 bits.D_{\max}^{0.2} = \log_2 \max\!\big(\tfrac{0.6}{0.2},\ \tfrac{0.4}{0.8}\big) = \log_2 3 \approx 1.585 \text{ bits}.

El presupuesto del 20%20\% rebaja el peor caso de 22 a 1.5851.585 bits. Para llegar a 00 haría falta ε=0.80.2=0.6\varepsilon = |0.8 - 0.2| = 0.6.

Ejemplo resuelto 2 · Un umbral para la finitud

Problema. Sea p=(0.7, 0.3)p = (0.7,\ 0.3) y q=(1, 0)q = (1,\ 0). ¿A partir de qué ε\varepsilon deja de ser infinita DmaxεD_{\max}^{\varepsilon}?

Solución. Como q1=0q_1 = 0, el cociente (1p~0)/(1q0)(1 - \tilde p_0)/(1 - q_0) solo es finito si p~0=1\tilde p_0 = 1, es decir, si vaciamos por completo la segunda coordenada. Eso cuesta desplazar una masa de 0.30.3, o sea una distancia de traza de 0.30.3. Por tanto:

  • para ε<0.3\varepsilon < 0.3, sigue siendo Dmaxε=+D_{\max}^{\varepsilon} = +\infty;
  • en ε=0.3\varepsilon = 0.3 se alcanza p~=(1,0)=q\tilde p = (1,0) = q y Dmax0.3=0D_{\max}^{0.3} = 0.

El umbral es exactamente la masa que pp tiene fuera del soporte de qq. Suavizar solo domestica el infinito cuando el presupuesto alcanza para retirar esa masa problemática.

Ejercicios

Ejercicio 1

Justifica, sin cálculo, por qué Dmaxε(ρσ)D_{\max}^{\varepsilon}(\rho\|\sigma) es una función no creciente de ε\varepsilon.

Solución

Si εε\varepsilon' \ge \varepsilon, entonces Bε(ρ)Bε(ρ)\mathcal{B}^{\varepsilon}(\rho) \subseteq \mathcal{B}^{\varepsilon'}(\rho): la bola pequeña está contenida en la grande. Minimizar sobre un conjunto mayor no puede dar un valor más alto, así que DmaxεDmaxεD_{\max}^{\varepsilon'} \le D_{\max}^{\varepsilon}. Es el argumento general de que «más candidatos, mínimo menor o igual».

Ejercicio 2

Con p=(0.6, 0.4)p = (0.6,\ 0.4) y q=(0.3, 0.7)q = (0.3,\ 0.7), calcula DmaxεD_{\max}^{\varepsilon} para ε=0.1\varepsilon = 0.1 (distancia de traza).

Solución

El peor cociente es p0/q0=0.6/0.3=2p_0/q_0 = 0.6/0.3 = 2, así que conviene bajar p0p_0 hacia q0=0.3q_0 = 0.3: de 0.60.6 a 0.50.5. Entonces Dmax0.1=log2max ⁣(0.50.3, 0.50.7)=log2530.737D_{\max}^{0.1} = \log_2 \max\!\big(\tfrac{0.5}{0.3},\ \tfrac{0.5}{0.7}\big) = \log_2 \tfrac{5}{3} \approx 0.737 bits (partía de log22=1\log_2 2 = 1 bit sin suavizar).