Articulo de referencia

Diagonalización de aproximación conjunta de matrices propias

La Diagonalización Conjunta de Matrices Propias (JADE) es un algoritmo para el análisis de componentes independientes que separa las señales mixtas observadas en señales fuente ...

La Diagonalización Conjunta de Matrices Propias (JADE) es un algoritmo para el análisis de componentes independientes que separa las señales mixtas observadas en señales fuente latentes mediante la explotación de momentos de cuarto orden . [ 1 ] Los momentos de cuarto orden son una medida de no gaussianidad, que se utiliza como un indicador para definir la independencia entre las señales fuente . La motivación para esta medida es que las distribuciones gaussianas poseen curtosis de exceso cero , y dado que la no gaussianidad es una suposición canónica del ICA, JADE busca una rotación ortogonal de los vectores mixtos observados para estimar vectores fuente que poseen altos valores de curtosis de exceso.

Algoritmo

Dejarincógnita=(incógnitaij)Rmetro×norte{\displaystyle \mathbf {X} =(x_{ij})\in \mathbb {R} ^{m\times n}}denota una matriz de datos observados cuyanorte{\displaystyle n}Las columnas corresponden a observaciones demetro{\displaystyle m}Vectores mixtos de -varianzas. Se supone queincógnita{\displaystyle \mathbf {X} }está preblanqueada , es decir, sus filas tienen una media muestral igual a cero y una covarianza muestral es lametro×metro{\displaystyle m\times m}matriz identidad dimensional, es decir,

1nortej=1norteincógnitaij=0y1norteincógnitaincógnita=Imetro{\displaystyle {\frac {1}{n}}\sum _{j=1}^{n}x_{ij}=0\quad {\text{y}}\quad {\frac {1}{n}}\mathbf {X} {\mathbf {X} }^{\prime }=\mathbf {I} _{m}}.

Aplicación de JADE aincógnita{\displaystyle \mathbf {X} }implica

  1. cálculo de cumulantes de cuarto orden deincógnita{\displaystyle \mathbf {X} }y luego
  2. optimizar una función de contraste para obtener unametro×metro{\displaystyle m\times m}matriz de rotaciónO{\displaystyle O}

para estimar los componentes de la fuente dados por las filas de lametro×norte{\displaystyle m\times n}matriz dimensionalZ:=O1incógnita{\displaystyle \mathbf {Z} :=\mathbf {O} ^{-1}\mathbf {X} } . [ 2 ]

Referencias

  1. Cardoso, Jean-François; Souloumiac, Antoine (1993). "Formación de haces ciega para señales no gaussianas". IEE Proceedings F - Radar and Signal Processing . 140 (6): 362– 370. CiteSeerX 10.1.1.8.5684 . doi : 10.1049/ip-f-2.1993.0054 . 
  2. Cardoso, Jean-François (enero de 1999). "Contrastes de alto orden para el análisis de componentes independientes". Neural Computation . 11 (1): 157– 192. CiteSeerX 10.1.1.308.8611 . doi : 10.1162/089976699300016863 .