Intuición detrás de la conmutatividad de convolución en sistemas LTI

7

¿Por qué la convolución es conmutativa, ya que parece tratar dos señales de manera diferente en un sistema LTI?

Si imagina con como señal de entrada como respuesta al impulso de un sistema LTI A, ¿cómo tiene sentido ese sistema LTI? B con la entrada y la respuesta al impulso genera exactamente la misma salida ?y[norte]=X[norte]h[norte]X[norte]h[norte]h[norte]X[norte]y[norte]

NocheLluvia23
fuente
1
y[norte]=X[norte]h[norte] no es una notación sensible. Debería ser solo . y=Xh
Leftaroundabout
3
En realidad, he visto la notación utilizada con frecuencia. No veo lo que no es sensato al respecto; se lee como "la secuencia se encuentra convolucionando la secuencia con la secuencia ". y[norte]=X[norte]h[norte]y[norte]X[norte]h[norte]
Jason R

Respuestas:

3

En un sistema de tiempo discreto como el que tiene, el número (aquí es un entero fijo) es una suma de la forma que se puede reorganizar mediante un cambio de variables (reemplace por ) aEntonces, la conmutatividad de la convolución es trivial. El problema es la interpretación que le pones. Como señala la respuesta de Laurent Duval / s, los sistemas A y B no son equivalentes en ningún sentido del término. Si la señaly[n0]n0

k=h[k]x[n0k]
kn0
=h[n0]x[].
xfueron reemplazados por una señal diferente , entonces el sistema A tendría salida , pero no obtendría la misma salida si el sistema B estaban emocionados por ; la respuesta al impulso del sistema B continúa siendo , y el sistema B tiene por lo tanto salida .x^y^=hx^y^hxxx^=x^xhx^

Dilip Sarwate
fuente
Gracias por dar (algún) sentido a mi respuesta, que ahora parece poco clara.
Laurent Duval
2

Imagine un sistema que acepta un solo número como su entrada, y multiplica ese número con otro número . ¿Le sorprendería que otro sistema que multiplica su entrada con el número da la misma salida que el primer sistema cuando se alimenta con el número como entrada? Si no es así, tampoco debería sorprendernos que la salida de un sistema LTI con respuesta de impulso y entrada proporcione la misma salida que otro sistema LTI con respuesta de impulso y entrada .XhXhh[norte]X[norte]X[norte]h[norte]

O, en lenguaje matemático, para el caso de tiempo discreto:

(xh)[n]=kx[k]h[nk]|m=nk=mx[nm]h[m]=(hx)[n]

Matt L.
fuente
1

Piensa en la convolución xhcomo tomar algunas copias demoradas dex resumidos juntos, cada uno con una amplitud leída de la entrada de hcon ese retraso Imaginémoslo con señales de impulso disperso: (ignore el desplazamiento vertical, eso es solo para desordenar la trama)Visualización de cómo se construye una convolución

Ahora, si cambia la circunvolución, todo lo que cambia es la noción de qué escala de tiempo es "retraso de la señal" y cuál simplemente "tiempo de la señal original". Terminas con el mismo resultado.

Visualizando la convolución al revés.

Código fuente (Haskell con diagrama dinámico ):

import Graphics.Dynamic.Plot.R2
import Control.Monad
import Data.List
import System.Random

main :: IO ()
main = do
   -- Times of the impulses in signal 𝑥
   txs <- scanl (+) 0 <$> replicateM 16 (randomRIO (0,0.1))
   -- Amplitudes of the impulses in 𝑥
   xs <- replicateM 8 (randomRIO (0,1))
   -- Times of the echoes in convolution kernel ℎ
   ths <- scanl (+) 0 <$> replicateM 16 (randomRIO (0,0.4))
   -- Amplitudes of the echoes in ℎ
   hs <- (1:) <$> replicateM 9 (randomRIO (0,1.0))

   plotWindow [plotLatest
     [ plotDelay 0.5 $ plotMultiple
    [ legendName name $ lineSegPlot
                             [(t,y+y0) | (t,x) <- sig, y<-[0,x,0]]
        | (    sig,                     name,   y0 ) <-
           [ ( zip txs xs,              "𝑥" ,   0  )
           , ( zip ths hs',             "ℎ" ,   1  )
           , ( [ (tx+th,x*h)
               | (tx,x)<-zip txs xs
               , (th,h)<-zip ths hs' ], "𝑥⋆ℎ", 0.5 )
           ]
        ]
      | hs' <- cycle . tail $ inits hs]
    , xInterval (0,2) ]
   return ()
a la izquierda
fuente
1

Recomiendo encarecidamente 1 definir la convolución de una manera ligeramente diferente, a saber, como:

(Xh)[norte]=yo+j=norteX[yo]h[j]

Puede ver que esta es solo la definición estándar con un cambio de variables, es decir, elegir j=norte-yo. (Necesitas ser un poco más claro sobre el rangoyo y j están atropellando, es decir, hace una diferencia si yo,jnorte versus yo,jZ.)

Esta definición hace obvia la conmutatividad. Sin embargo, otra perspectiva sobre esto es que si tienes un polinomioX(z)=yoX[yo]zyo y un polinomio H(z)=yoh[yo]zyo y los multiplicas como polinomios, X(z)H(z), usted obtiene

X(z)H(z)=yo(Xh)[yo]zyo
(Esto se generaliza a las series de poder formales). Simplemente escriba la multiplicación usted mismo en algunos ejemplos para verificar esto y tal vez probarlo inductivamente. Esta perspectiva nuevamente hace que la conmutatividad sea completamente obvia. Esta perspectiva también hace que el teorema de convolución para la transformación Z también sea bastante obvio.

1 Escribí una publicación de blog completa sobre esto.

Derek Elkins dejó SE
fuente
bueno, a alguien no le gustó tu respuesta, pero a mí sí (y cancelé su "-1"). Nunca he visto esta notación antes:
yo+j=norteX[yo]h[j]
Y creo que es algo lindo. tal vez, en el futuro, podría pensar que es desagradable, pero ahora creo que es lindo.
Robert Bristow-Johnson
Hay algo interesante en este argumento "diagonal". Intentaré encontrar una versión gráfica
Laurent Duval,
1

Como complemento de las respuestas anteriores, apenas mirar la salida no le permite identificar lo que tenía como entrada / función del sistema. Tomemos, por ejemplo, los siguientes sistemas (interconexión de) que producen la misma salida.


ingrese la descripción de la imagen aquí

Para el sistema anterior tenemos:

y[norte]=(Xh)[norte]
Y(z)=X(z).H(z)=H(z).X(z)


ingrese la descripción de la imagen aquí

Para el sistema anterior tenemos:

y[norte]=((δX)h)[norte]
Y(z)=1.X(z).H(z)=X(z).H(z)=H(z).X(z)


ingrese la descripción de la imagen aquí Para el sistema anterior tenemos:

y[norte]=((δαX)1αh)[norte]
Y(z)=1)αX(z).1αH(z)=X(z).H(z)=H(z).X(z)


Por supuesto, son posibles muchas otras combinaciones que dan el mismo resultado, incluido el intercambio de roles de h[norte] y X[norte].

Tenga en cuenta que esta observación es la principal preocupación del campo llamado Identificación del sistema ciego : intente identificar la entrada y / o la función del sistema utilizando solo los datos de salida.

Como no podemos resolver el problema como se indicó, se necesita más información secundaria (por ejemplo, estadísticas de salida), canales más diversos también (para compensar los ceros, es decir, frecuencias por las que no puede pasar ninguna señal) y una señal de entrada suficientemente rica que puede Ser útil para recuperar la función del sistema.

Tener una intuición para su pregunta sobre

¿Por qué la convolución es conmutativa, ya que parece tratar dos señales de manera diferente en un sistema LTI?

mirar Y(z)=X(z)H(z)como la salida de dos sistemas LTI en serie. Puedes intercambiarlos, aplicando el segundo y luego el primero; Eso es conmutatividad. También siX(z) y H(z) ambos tienen la misma (o infinita) banda de frecuencia, entonces podemos recuperarlos (hasta una multiplicación escalar), de lo contrario, definitivamente perderemos algunas frecuencias y la recuperación total no es posible.

Aprender_y_ Compartir
fuente
1

Existen diferencias fundamentales de concepto entre señales y sistemas . Explicaré esto a través de la idea de consistencia de la unidad (ver por ejemplo). Sin embargo, para los sistemas LTI, las señales y los sistemas se vuelven duales por convolución, ya que este último es conmutativo. Dos digresiones primero, debido a la mención en la respuesta de @Dilip Sarwate .

  • Digresión 1: los sistemas LTI pueden tener la misma salida para diferentes señales

Si dos sistemas diferentes proporcionan las mismas salidas para algunas señales de entrada, esto significa que comparten algunas propiedades. Pero si sus salidas son iguales para todas las entradas, entonces esencialmente tienen la misma respuesta de impulso, y son prácticamente los mismos sistemas.

Por ejemplo, imagine que tiene un seno de entrada a frecuencia F. Si ambos sistemas cortan la frecuencia arribaF-ϵ, ambos tienen el mismo comportamiento para esa señal, pero pueden ser dos sistemas de paso bajo diferentes, se necesitan más señales para distinguirlos.

  • Digresión 2: dos señales de entrada diferentes pueden tener la misma salida a través de un sistema LTI dado

Por ejemplo, una señal constante igual a uno, o una señal de 2 periódicos con {2,0 0} valores producen la misma salida para 2norte-procesadores de promedio.

De vuelta a tu pregunta . Un sistemaS convierte entradas X en salidas Y, respectivamente con unidades físicas tuX y tuY. Entonces, un sistema puede verse como un convertidor de unidades, formalmente con unidad internatuY/ /tuX. En general, el sistema es "fijo", mientras que las entradas pueden variar. Entonces, no hay razón por la cualS y X debería jugar el mismo papel.

Sin embargo, cuando uno considera los sistemas LTI, de repente las propiedades del sistema pueden transferirse de alguna manera a las señales, y viceversa (siempre que la convolución esté bien definida). Esto está relacionado con el hecho de que la convolución conmuta con los cambios. Para simplificar, imagine un sistema de "tres toques", conz-transformar la respuesta hlz-l+hmetroz-metro+hnortez-norte. Puede convertir esto directamente en un banco de filtros de tres bandas, con una sola entrada y respuestas respectivashlz-l, hmetroz-metro y hnortez-norte. Cada rama solo proporciona, para cada entrada, un factor de escala y un retraso.

Pero lo mismo sucede con las señales: cada entrada X={...,Xl,...,Xmetro,...,Xnorte,...} se puede dividir en componentes escalares:

X=...+Xlδl+...+Xmetroδmetro+Xnorteδnorte+...
dónde δdenota el símbolo de Kronecker. Debido a la linealidad, cada componente podría alimentarse a través del sistema lineal. Cuando todo (señal y sistema) se divide de esta manera, los cálculos son solo un montón deXkδk pasando por un par de hyoz-yo, que son fundamentalmente las mismas operaciones: un factor / una amplitud y una muestra retardada / un operador de retardo. En otras palabras,Xkδk ir a través hyoz-yo produce el mismo resultado que hkδk ir a través Xyoz-yo, porque el producto hkXyo es conmutativo (y conserva la consistencia de la unidad), y también retrasa los desplazamientos.

En otras palabras, un LTI solo produce una suma ponderada con pesos h en muestras de entrada de X: hyoXk-yo, que puede leerse así como una suma ponderada con pesas X en muestras de entrada de h: Xyohk-yo. Sin embargo, para la consistencia de la unidad, uno debe cambiar las unidades deX y h.

Esta intercambiabilidad entre señales y sistemas en el LTI parece estar en juego (a primera vista) en la expresión de polifase / modulación de los bancos de filtros , o en el filtrado coincidente.

Laurent Duval
fuente
2
Esto parece estar respondiendo una pregunta diferente.
user253751
Realmente creo que estaba respondiendo una parte de la pregunta en su formulación inicial. Ya no, déjame modificar eso
Laurent Duval
@LaurentDuval ¿por qué esto no responde la pregunta (más)?
AlexTP
0

Tienes razón. Es completamente absurdo pensar que la respuesta al impulso de un sistema LTI puede ser reemplazada por la señal de entrada y viceversa y, sin embargo, producen el mismo resultado.

Como ejemplo, considere un filtro de paso bajo con respuesta de impulso IIR h[norte] que es alimentado por las muestras de forma de onda del habla X[norte]para producir una versión filtrada de paso bajo del discurso. Sin embargo, intercambiando los roles de entrada de voz y respuesta del impulso del sistema LTIh[norte] se convierte en un absurdo en un entorno práctico.

Sin embargo , matemáticamente ese es el caso. E incluso puede encontrar aplicaciones de ejemplo que pueden beneficiarse de dicho intercambio. Se da una explicación matemática en la respuesta de Matt.

Grasa32
fuente