Articulo de referencia

Teorema de factorización polar

En el transporte óptimo , una rama de las matemáticas, la factorización polar de campos vectoriales es un resultado básico debido a Brenier (1987), [ 1 ] con antecedentes de Kno...

En el transporte óptimo , una rama de las matemáticas, la factorización polar de campos vectoriales es un resultado básico debido a Brenier (1987), [ 1 ] con antecedentes de Knott-Smith (1984) [ 2 ] y Rachev (1985), [ 3 ] que generaliza muchos resultados existentes entre los que se encuentran la descomposición polar de matrices reales y el reordenamiento de funciones de valor real.

El teorema

Notación. Denotarξ#μ{\displaystyle \xi _{\#}\mu }la medida de la imagen deμ{\displaystyle \mu }a través del mapaξ{\displaystyle \xi }.

Definición : Mapa que conserva la medida.(incógnita,μ){\displaystyle (X,\mu )}y(Y,ν){\displaystyle (Y,\nu )}sean algunos espacios de probabilidad yσ:incógnitaY{\displaystyle \sigma :X\rightarrow Y}un mapa medible. Luego,σ{\displaystyle \sigma }Se dice que es una medida que preserva si y solo siσ#μ=ν{\displaystyle \sigma _{\#}\mu =\nu }, dónde#{\displaystyle \#}es la medida de impulso . En otras palabras: por cadaν{\displaystyle \nu }-subconjunto medibleΩ{\displaystyle \Omega }deY{\displaystyle Y},σ1(Ω){\displaystyle \sigma ^{-1}(\Omega )}esμ{\displaystyle \mu }-medible yμ(σ1(Ω))=ν(Ω){\displaystyle \mu (\sigma ^{-1}(\Omega ))=\nu (\Omega )}Esto último es equivalente a:

incógnita(Fσ)(incógnita)μ(dincógnita)=incógnita(σF)(incógnita)μ(dincógnita)=YF(y)(σ#μ)(dy)=YF(y)ν(dy){\displaystyle \int _{X}(f\circ \sigma )(x)\mu (dx)=\int _{X}(\sigma ^{*}f)(x)\mu (dx)=\int _{Y}f(y)(\sigma _{\#}\mu )(dy)=\int _{Y}f(y)\nu (dy)}

dóndeF{\displaystyle f}esν{\displaystyle \nu }-integrable yFσ{\displaystyle f\circ \sigma }esμ{\displaystyle \mu }-integrable.

Teorema. Consideremos un mapaξ:ΩRd{\displaystyle \xi :\Omega \rightarrow R^{d}} dondeΩ{\displaystyle \Omega }es un subconjunto convexo deRd{\displaystyle R^{d}}, yμ{\displaystyle \mu }una medida sobreΩ{\displaystyle \Omega }que es absolutamente continua. Supongamos queξ#μ{\displaystyle \xi _{\#}\mu }es absolutamente continua. Entonces hay una función convexa.φ:ΩR{\displaystyle \varphi :\Omega \rightarrow R} y un mapaσ:ΩΩ{\displaystyle \sigma :\Omega \rightarrow \Omega } preservandoμ{\displaystyle \mu }de tal manera que

ξ=(φ)σ{\displaystyle \xi =\left(\nabla \varphi \right)\circ \sigma }

Además,φ{\displaystyle \nabla \varphi}yσ{\displaystyle \sigma }están definidos de forma única en casi todas partes. [ 1 ] [ 4 ]

Aplicaciones y conexiones

Dimensión 1

En la dimensión 1, y cuandoμ{\displaystyle \mu }es la medida de Lebesgue sobre el intervalo unitario, el resultado se especializa en el teorema de Ryff. [ 5 ] Cuandod=1{\displaystyle d=1}yμ{\displaystyle \mu }es la distribución uniforme sobre[0,1]{\displaystyle \left[0,1\right]}, la descomposición polar se reduce a

ξ(t)=Fincógnita1(σ(t)){\displaystyle \xi \left(t\right)=F_{X}^{-1}\left(\sigma \left(t\right)\right)}

dóndeFincógnita{\displaystyle F_{X}}es la función de distribución acumulativa de la variable aleatoriaξ(U){\displaystyle \xi \left(U\right)}yU{\displaystyle U}tiene una distribución uniforme sobre[0,1]{\displaystyle \left[0,1\right]}.Fincógnita{\displaystyle F_{X}}se supone que es continuo, yσ(t)=Fincógnita(ξ(t)){\displaystyle \sigma \left(t\right)=F_{X}\left(\xi \left(t\right)\right)}conserva la medida de Lebesgue sobre[0,1]{\displaystyle \left[0,1\right]}.

Descomposición polar de matrices

Cuandoξ{\displaystyle \xi }es un mapa lineal yμ{\displaystyle \mu }es la distribución normal gaussiana , el resultado coincide con la descomposición polar de matrices . Suponiendoξ(incógnita)=METROincógnita{\displaystyle \xi \left(x\right)=Mx}dóndeMETRO{\displaystyle M}es un invertibled×d{\displaystyle d\times d}matriz y considerandoμ{\displaystyle \mu }elnorte(0,Id){\displaystyle {\mathcal {N}}\left(0,I_{d}\right)}medida de probabilidad, la descomposición polar se reduce a

METRO=SO{\displaystyle M=SO}

dóndeS{\displaystyle S}es una matriz simétrica definida positiva , yO{\displaystyle O}una matriz ortogonal . La conexión con la factorización polar esφ(incógnita)=incógnitaSincógnita/2{\displaystyle \varphi \left(x\right)=x^{\top }Sx/2}que es convexa, yσ(incógnita)=Oincógnita{\displaystyle \sigma \left(x\right)=Ox}que conserva elnorte(0,Id){\displaystyle {\mathcal {N}}\left(0,I_{d}\right)}medida.

descomposición de Helmholtz

Los resultados también permiten recuperar la descomposición de Helmholtz . DejandoincógnitaV(incógnita){\displaystyle x\rightarrow V\left(x\right)}Si es un campo vectorial suave, entonces se puede escribir de una manera única como

V=w+pag{\displaystyle V=w+\nabla p}

dóndepag{\displaystyle p}es una función real suave definida enΩ{\displaystyle \Omega }, único salvo una constante aditiva, yw{\displaystyle w}es un campo vectorial suave sin divergencia, paralelo al límite deΩ{\displaystyle \Omega }.

La conexión se puede ver asumiendoμ{\displaystyle \mu }es la medida de Lebesgue en un conjunto compactoΩRnorte{\displaystyle \Omega \subset R^{n}}y escribiendoξ{\displaystyle \xi }como una perturbación del mapa identidad

ξϵ(incógnita)=incógnita+ϵV(incógnita){\displaystyle \xi _{\epsilon }(x)=x+\epsilon V(x)}

dóndeϵ{\displaystyle \epsilon }es pequeño. La descomposición polar deξϵ{\displaystyle \xi _{\epsilon }}es dado porξϵ=(φϵ)σϵ{\displaystyle \xi _{\epsilon }=(\nabla \varphi _{\epsilon })\circ \sigma _{\epsilon }}. Luego, para cualquier función de pruebaF:RnorteR{\displaystyle f:R^{n}\rightarrow R}Se cumple lo siguiente:

ΩF(incógnita+ϵV(incógnita))dincógnita=ΩF((φϵ)σϵ(incógnita))dincógnita=ΩF(φϵ(incógnita))dincógnita{\displaystyle \int _{\Omega }f(x+\epsilon V(x))dx=\int _{\Omega }f((\nabla \varphi _{\epsilon })\circ \sigma _{\epsilon }\left(x\right))dx=\int _{\Omega }f(\nabla \varphi _{\epsilon }\left(x\right))dx}

donde el hecho de queσϵ{\displaystyle \sigma _{\epsilon }}Se estaba conservando la medida de Lebesgue que se utilizó en la segunda igualdad.

De hecho, comoφ0(incógnita)=12incógnita2{\displaystyle \textstyle \varphi _{0}(x)={\frac {1}{2}}\Vert x\Vert ^{2}}, uno puede expandirφϵ(incógnita)=12incógnita2+ϵpag(incógnita)+O(ϵ2){\displaystyle \textstyle \varphi _{\epsilon }(x)={\frac {1}{2}}\Vert x\Vert ^{2}+\epsilon p(x)+O(\epsilon ^{2})}y por lo tantoφϵ(incógnita)=incógnita+ϵpag(incógnita)+O(ϵ2){\displaystyle \textstyle \nabla \varphi _{\epsilon }\left(x\right)=x+\epsilon \nabla p(x)+O(\epsilon ^{2})}. Como resultado,Ω(V(incógnita)pag(incógnita))F(incógnita))dincógnita{\displaystyle \textstyle \int _{\Omega }\left(V(x)-\nabla p(x)\right)\nabla f(x))dx}para cualquier funcionamiento suaveF{\displaystyle f}, lo que implica quew(incógnita)=V(incógnita)pag(incógnita){\displaystyle w\left(x\right)=V(x)-\nabla p(x)}es libre de divergencia. [ 1 ] [ 6 ]

Véase también

Referencias

  1. 1 2 3 Brenier, Yann (1991). "Factorización polar y reordenamiento monótono de funciones con valores vectoriales" (PDF) . Communications on Pure and Applied Mathematics . 44 (4): 375– 417. doi : 10.1002/cpa.3160440402 . Recuperado el 16 de abril de 2021 .
  2. Knott, M.; Smith, CS (1984). "Sobre la asignación óptima de distribuciones" . Journal of Optimization Theory and Applications . 43 : 39–49 . doi : 10.1007/BF00934745 . S2CID 120208956. Consultado el 16 de abril de 2021 . 
  3. Rachev, Svetlozar T. (1985). "El problema de transferencia de masa de Monge-Kantorovich y sus aplicaciones estocásticas" (PDF) . Theory of Probability & Its Applications . 29 (4): 647– 676. doi : 10.1137/1129093 . Recuperado el 16 de abril de 2021 .
  4. Santambrogio, Filippo (2015). Transporte óptimo para matemáticos aplicados . Nueva York: Birkäuser. CiteSeerX 10.1.1.726.35 . 
  5. Ryff, John V. (1965). "Órbitas de funciones L1 bajo transformación doblemente estocástica" . Transactions of the American Mathematical Society . 117 : 92–100 . doi : 10.2307/1994198 . JSTOR 1994198. Consultado el 16 de abril de 2021 . 
  6. Villani, Cédric (2003). Temas en transporte óptimo . Sociedad Matemática Americana.