¿Cómo cambia la similitud de coseno después de una transformación lineal?


9

¿Existe una relación matemática entre:

  • la similitud del coseno de dos vectores y , ysim(A,B)AB
  • ¿La similitud del coseno de y , escalada de manera no uniforme a través de una matriz dada ? Aquí es una matriz diagonal dada con elementos desiguales en la diagonal.sim(MA,MB)ABMM

Intenté repasar los cálculos, pero no pude alcanzar un enlace simple / interesante (expresión). Me pregunto si hay uno.


Por ejemplo, los ángulos no se conservan en una escala no uniforme, pero ¿cuál es la relación entre los ángulos originales y los posteriores a la escala no uniforme? ¿Qué se puede decir sobre el enlace entre un conjunto de vectores S1 y otro conjunto de vectores S2, donde S2 se obtiene escalando de manera no uniforme S1?


@whuber, gracias! Sí, M es una matriz dada (una matriz de escala, por lo tanto, una matriz diagonal, sin otras restricciones). En cierto sentido, quería saber qué sucede (en términos de la similitud del coseno para cualquier par de vectores) con un espacio vectorial que sufre una escala no lineal.
turdus-merula

2
Vale la pena señalar que si todos los factores de escala no son negativos (como uno supondría naturalmente), entonces todas las matrices simétricas positivas definidas pueden considerarse matrices de "escala". La relación que busca se usa ampliamente, entre otras cosas , en el estudio y la descripción de la distorsión en las proyecciones de mapas. Allí, el interés se centra en los ángulos máximo y mínimo en la superficie de la tierra que se asociarían con dos direcciones perpendiculares en el mapa. Existe una relación directa entre estos ángulos y las proporciones de los dos factores de escala.
whuber

Respuestas:


8

Debido a que es bastante general, y el cambio en la similitud del coseno depende de y particulares y su relación con , no es posible una fórmula definida. Sin embargo, existen límites prácticamente calculables de cuánto puede cambiar la similitud del coseno . Se pueden encontrar extremizando el ángulo entre y dado que la similitud del coseno entre y es un valor especificado, digamos (donde es el ángulo entre y ). La respuesta nos dice cuánto cualquier ánguloA B M M A M B A B cos ( 2 ϕ ) 2 ϕ A B 2 ϕ MMABMMAMBABcos(2ϕ)2ϕAB2ϕposiblemente se puede doblar por la transformación .M

Los cálculos amenazan con ser desordenados. Algunas opciones inteligentes de notación, junto con algunas simplificaciones preliminares, reducen el esfuerzo. Resulta que la solución en dos dimensiones revela todo lo que necesitamos saber. Este es un problema manejable, que depende solo de una variable real , que se resuelve fácilmente utilizando técnicas de cálculo. Un argumento geométrico simple extiende esta solución a cualquier cantidad de dimensiones .nθn

Preliminares matemáticos

Por definición, el coseno del ángulo entre dos vectores y se obtiene normalizándolos a la unidad de longitud y tomando su producto. Así,BAB

AB(AA)(BB)=cos(2ϕ)

y, escribiendo , el coseno del ángulo entre las imágenes de y bajo la transformación esA B MΣ=MMABM

(1)(MA)(MB)((MA)(MA))((MB)(MB))=AΣB(AΣA)(BΣB).

Observe que solo importa en el análisis,Σ no sí. Por lo tanto, podemos explotar la Descomposición de valor singular (SVD) de para simplificar el problema. Recuerde que esto expresa como un producto (de derecha a izquierda) de una matriz ortogonal , una matriz diagonal y otra matriz ortogonal :M M V D UMMMVDU

M=UDV.

En otras palabras, existe una base de vectores privilegiados (las columnas de ) sobre los cuales actúa reescalando cada por separado mediante la entrada diagonal de (que llamaré ) y luego aplicando una rotación (o antirrotación) al resultado. Esa rotación final no cambiará ninguna longitud o ángulo y, por lo tanto, no debería afectar a . Puedes ver esto formalmente con el cálculo V M e i i th D d i U Σe1,,enVMeiithDdiUΣ

Σ=MM=(UDV)(UDV)=VD(UU)DV=VD2V.

En consecuencia, para estudiar podemos reemplazar libremente por cualquier otra matriz que produzca los mismos valores en . Al ordenar para que disminuya de tamaño (y suponiendo que no sea idénticamente cero), una buena elección de esM ( 1 ) e i d i M MΣM(1)eidiMM

M=1d1DV.

Los elementos diagonales de son(1/d1)D

1=d1/d1λ2=d2/d1λ3=d3/d1λn=dn/d10.

Específicamente, el efecto de (ya sea en su forma original o modificada) en todos los ángulos está completamente determinado por el hecho de queM

Mei=λiei.

Análisis de un caso especial.

Deje . Debido a que cambiar las longitudes de los vectores no cambia el ángulo entre ellos, podemos suponer que y son vectores unitarios. En el plano, todos estos vectores pueden designarse por el ángulo que forman con , lo que nos permite escribirA B e 1n=2ABe1

A=cos(θϕ)e1+sin(θϕ)e2.

Por lo tanto

B=cos(θ+ϕ)e1+sin(θ+ϕ)e2.

(Vea la figura a continuación).

Aplicar es simple: corrige las primeras coordenadas de y y multiplica sus segundas coordenadas por . Por lo tanto, el ángulo de a esA B λ 2 M A M BMABλ2MAMB

f(θ)=arctan(λ2tan(θ+ϕ))arctan(λ2tan(θϕ)).

Como es una función continua, esta diferencia de ángulos es una función continua de . De hecho, es diferenciable. Esto nos permite encontrar los ángulos extremos inspeccionando los ceros de la derivada . Esa derivada es fácil de calcular: es una relación de funciones trigonométricas. Los ceros solo pueden aparecer entre los ceros de su numerador, así que no nos molestemos en calcular el denominador. Obtenemosθ f ( θ )Mθf(θ)

f(θ)=λ2(1λ2)(λ2+1)sin(2θ)sin(2ϕ).

Los casos especiales de , y se entienden fácilmente: corresponden a las situaciones en las que es de rango reducido (y así aplasta todos los vectores en una línea); donde es un múltiplo de la matriz de identidad; y donde y son paralelos (de donde el ángulo entre ellos no puede cambiar, independientemente de ). El caso está excluido por la condición .λ 2 = 1λ2=0λ2=1M M A B θ λ 2 = - 1 λ 20ϕ=0MMABθλ2=1λ20

Además de estos casos especiales, los ceros ocurren solo donde : es decir, o . Esto significa que la línea determinada por divide el ángulo . Ahora sabemos que los valores extremos del ángulo entre y deben estar entre los valores de , así que calculemos:sin(2θ)=0θ = π / 2 e 1 A B M A M B f ( θ )θ=0θ=π/2e1ABMAMBf(θ)

f(0)=arctan(λ2tan(ϕ))arctan(λ2tan(ϕ))=2arctan(λ2tan(ϕ));f(π/2)=arctan(λ2tan(π/2+ϕ))arctan(λ2tan(π/2ϕ))=2arctan(λ2cot(ϕ)).

Los cosenos correspondientes son

(2)cos(f(0))=1λ22tan(ϕ)21+λ22tan(ϕ)2

y

(3)cos(f(π/2))=1λ22cot(ϕ)21+λ22cot(ϕ)2=tan(ϕ)2λ22tan(ϕ)2+λ22.

A menudo es suficiente entender cómo distorsiona los ángulos rectos. En este caso, , lo que lleva a , que puede insertar en las fórmulas anteriores.2 ϕ = π / 2 tan ( ϕ ) = cot ( ϕ ) = 1M2ϕ=π/2tan(ϕ)=cot(ϕ)=1

Tenga en cuenta que cuanto menor sea , más extremos serán estos ángulos y mayor será la distorsión.λ2

Figura que muestra cuatro configuraciones

Esta figura muestra cuatro configuraciones de los vectores y separadas por un ángulo de . El círculo unitario y su imagen elíptica debajo de están sombreados como referencia (con la acción de reescalada uniformemente para hacer ). Las figura encabezamientos indican el valor de , el punto medio de y . Lo más parecido que puede llegar a ser y cuando se transforma por es una configuración como la de la izquierda conB 2 ϕ = π / 3 M M λ 1 = 1 θAB2ϕ=π/3MMλ1=1θB A B M θ = 0 θ = π / 2ABABMθ=0. Lo más alejado que pueden estar es una configuración como la de la derecha con . Se muestran dos posibilidades intermedias.θ=π/2

Solución para todas las dimensiones.

Hemos visto cómo actúa expandiendo cada dimensión por un factor . Esto distorsionará la unidad de esfera en un elipsoide. El determina sus ejes principales. Los son las distancias desde el origen, a lo largo de estos ejes, al elipsoide. En consecuencia, el más pequeño, , es la distancia más corta (en cualquier dirección) desde el origen al elipsoide y el más grande, , es la distancia más lejana (en cualquier dirección) desde el origen al elipsoide.i λ i {Miλie i λ i λ n λ 1{A|AA=1}eiλiλnλ1

En dimensiones superiores , y son parte de un subespacio bidimensional. asigna el círculo unitario en este subespacio en la intersección del elipsoide con un plano que contiene y . Esta intersección, al ser una distorsión lineal de un círculo, es una elipse. Obviamente, la distancia más lejana a esta elipse no es más que y la distancia más corta no es menos que .n>2B M M A M B λ 1 = 1 λ nABMMAMBλ1=1λn

Como observamos al final de la sección anterior, la posibilidad más extrema es cuando y están situados en un plano que contiene dos de para los cuales la relación de correspondiente es lo más pequeña posible. Esto sucederá en el plano . Ya tenemos la solución para ese caso.B e i λ i e 1 , e nABeiλie1,en

Conclusiones

Los extremos de similitud de coseno alcanzables aplicando a dos vectores que tienen similitud de coseno están dados por y . Se logran situando y en ángulos iguales a una dirección en la que alarga al máximo cualquier vector (como la dirección ) y separándolos en una dirección en la que alarga mínimamente cualquier vector ( como la dirección ).cos ( 2 ϕM( 2 ) ( 3 )cos(2ϕ)(2)(3)B Σ = M M e 1 Σ e nABΣ=MMe1Σen

Estos extremos se pueden calcular en términos de la SVD de .M


¡Esta es una respuesta fantástica! ¡Muchas gracias por esta discusión detallada! Creo que tiene un error de signo en la ecuación (3) donde debería tener un signo menos general.
LFH

Estoy interesado en el caso en que el ángulo aproxima a cero y me gustaría obtener una desigualdad entre y . ¿Es cierto que, en función de su cálculo, solo necesito encontrar el más extremo (que es el más pequeño) y en este caso, la desigualdad asintótica viene dada por como ? 2 ϕ f λ n 2 λ n ϕ f 2 λ2ϕ2ϕfλnϕ2λnϕf2λn1ϕϕ0
LFH

6

Probablemente le interese:

(MA,MB)=AT(MTM)B,

Puede diagonalizar (o como lo llaman, PCA), lo que le dice que la similitud de bajo la transformación comporta proyectando en sus componentes principales, y posteriormente calculando similitud en este nuevo espacio. Para desarrollar esto un poco más, dejemos que los componentes principales sean con valores propios . Entonces A , B M A , B u i λ iMTM=UΣUTA,BMA,Buiλi

UB=i(ui,bi)ui, UA=i(ui,ai)ui,

que te da:

(MA,MB)=i=1n(ui,ai)(ui,bi)λi.

Observe que hay una escala aquí: los se . Cuando son vectores unitarios y si cada , entonces corresponde a una rotación, y obtienes: , que es equivalente a decir que los productos internos son invariantes bajo rotaciones. En general, el ángulo permanece igual cuando es una transformación conforme, que en este caso requiere que sea ​​invertible y que la descomposición polar de satisfaga con , es decir, . A , B λ i = 1 M sim ( M A , M B ) = sim ( A , B ) M M M M = O P P = a I M T M = a 2 IλiA,Bλi=1Msim(MA,MB)=sim(A,B)MMMM=OPP=aIMTM=a2I


1
Su declaración inicial del problema descuida la normalización de los vectores , , y necesarios para calcular la similitud del coseno. Tampoco parece que el análisis posterior aborde esta normalización. Tenga en cuenta, en particular, que las similitudes del coseno se conservan incluso cuando todos los valores propios son iguales a algún valor (positivo) que difiere de . Eso demuestra, incluso en este caso simple, que se puede decir mucho más. AM A M B 1BMAMB1
whuber

@whuber: la similitud del coseno se conserva exactamente cuando es una transformación conforme, que en este caso es equivalente a requerir que sea ​​invertible y , un múltiplo de la identidad. Dicho de otra manera, la descomposición polar de satisface , donde . Tienes razón en la normalización, pero, parece absurdo hablar de coseno similitud con vectores no normalizados . M M T M = a 2 I M M =MMMTM=a2IMP = a I A , BM=OPP=aIA,B
Alex R.

2
No es tonto en absoluto! Dado que esta "similitud" viene dada por el coseno del ángulo entre los vectores, tiene sentido para cualquiera de los dos vectores distintos de cero. Lo que quería decir con "mucho más se puede decir" es que los límites efectivos en el ángulo entre las imágenes de y se puede solicitar en términos del ángulo entre y y los valores propios de . B A B MABABM
whuber
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.