En el procesamiento de señales , la salida del filtro adaptado se obtiene correlacionando una señal retardada conocida , o plantilla , con una señal desconocida para detectar la presencia de la plantilla en la señal desconocida. [1] [2] Esto es equivalente a convolucionar la señal desconocida con una versión conjugada en el tiempo invertido de la plantilla. El filtro adaptado es el filtro lineal óptimo para maximizar la relación señal-ruido (SNR) en presencia de ruido estocástico aditivo .
Los filtros adaptados se utilizan habitualmente en radares , en los que se envía una señal conocida y se examina la señal reflejada en busca de elementos comunes de la señal saliente. La compresión de pulsos es un ejemplo de filtrado adaptado. Se denomina así porque la respuesta al impulso se adapta a las señales de pulso de entrada. Los filtros adaptados bidimensionales se utilizan habitualmente en el procesamiento de imágenes , por ejemplo, para mejorar la relación señal-ruido de las observaciones de rayos X. Otras aplicaciones destacables son la sismología y la astronomía de ondas gravitacionales .
El filtrado adaptado es una técnica de demodulación con filtros LTI (invariantes en el tiempo lineales) para maximizar la relación señal-ruido. [3] Originalmente también se conocía como filtro Norte . [4]
Derivación
Derivación mediante álgebra matricial
En la siguiente sección se obtiene el filtro adaptado para un sistema de tiempo discreto . La obtención para un sistema de tiempo continuo es similar, pero las sumas se reemplazan por integrales.
El filtro adaptado es el filtro lineal, , que maximiza la relación señal-ruido de salida .
donde es la entrada en función de la variable independiente , y es la salida filtrada. Aunque la mayoría de las veces expresamos los filtros como la respuesta al impulso de los sistemas de convolución, como se indicó anteriormente (consulte la teoría de sistemas LTI ), es más fácil pensar en el filtro adaptado en el contexto del producto interno , que veremos en breve.
Podemos derivar el filtro lineal que maximiza la relación señal-ruido de salida invocando un argumento geométrico. La intuición detrás del filtro adaptado se basa en correlacionar la señal recibida (un vector) con un filtro (otro vector) que es paralelo a la señal, maximizando el producto interno. Esto mejora la señal. Cuando consideramos el ruido estocástico aditivo, tenemos el desafío adicional de minimizar la salida debido al ruido al elegir un filtro que sea ortogonal al ruido.
Definamos formalmente el problema. Buscamos un filtro, , tal que maximicemos la relación señal-ruido de salida, donde la salida es el producto interno del filtro y la señal observada .
Nuestra señal observada consta de la señal deseada y ruido aditivo :
Definamos la matriz de autocorrelación del ruido, recordando que esta matriz tiene simetría hermítica , propiedad que será útil en la derivación:
donde denota la transpuesta conjugada de , y denota la expectativa (nótese que en caso de que el ruido tenga media cero, su matriz de autocorrelación es igual a su matriz de covarianza ).
Llamemos a nuestra salida, , el producto interno de nuestro filtro y la señal observada tal que
Ahora definimos la relación señal-ruido, que es nuestra función objetivo, como la relación entre la potencia de la salida debida a la señal deseada y la potencia de la salida debida al ruido:
Reescribimos lo anterior:
Deseamos maximizar esta cantidad eligiendo . Desarrollando el denominador de nuestra función objetivo, tenemos
Ahora, nuestro se convierte
Reescribiremos esta expresión con alguna manipulación de matrices. La razón de esta medida aparentemente contraproducente se hará evidente en breve. Aprovechando la simetría hermítica de la matriz de autocorrelación , podemos escribir
Nos gustaría encontrar un límite superior para esta expresión. Para ello, primero reconocemos una forma de la desigualdad de Cauchy-Schwarz :
es decir, que el cuadrado del producto interno de dos vectores solo puede ser tan grande como el producto de los productos internos individuales de los vectores. Este concepto vuelve a la intuición detrás del filtro adaptado: este límite superior se logra cuando los dos vectores y son paralelos. Reanudamos nuestra derivación expresando el límite superior de nuestro a la luz de la desigualdad geométrica anterior:
Nuestra valiente manipulación de matrices ha dado sus frutos. Vemos que la expresión para nuestro límite superior se puede simplificar enormemente:
Podemos alcanzar este límite superior si lo elegimos,
donde es un número real arbitrario. Para verificarlo, introducimos en nuestra expresión la salida :
Por lo tanto, nuestro filtro adaptado óptimo es
A menudo optamos por normalizar el valor esperado de la potencia de la salida del filtro debido al ruido a la unidad. Es decir, limitamos
Esta restricción implica un valor de , para el cual podemos resolver:
flexible
dándonos nuestro filtro normalizado,
Si queremos escribir la respuesta al impulso del filtro para el sistema de convolución, es simplemente la inversión temporal conjugada compleja de la entrada .
Aunque hemos derivado el filtro adaptado en tiempo discreto, podemos extender el concepto a sistemas de tiempo continuo si lo reemplazamos con la función de autocorrelación de tiempo continuo del ruido, asumiendo una señal continua , un ruido continuo y un filtro continuo .
Derivación por Lagrangiano
Como alternativa, podemos resolver el problema de maximización del filtro adaptado con un lagrangiano. Nuevamente, el filtro adaptado intenta maximizar la relación señal-ruido de salida ( ) de una señal determinista filtrada en ruido aditivo estocástico. La secuencia observada, nuevamente, es
con la matriz de autocorrelación de ruido,
La relación señal-ruido es
donde y .
Evaluando la expresión en el numerador, tenemos
y en el denominador,
La relación señal-ruido se convierte en
Si ahora restringimos el denominador a 1, el problema de maximizar se reduce a maximizar el numerador. Podemos entonces formular el problema utilizando un multiplicador de Lagrange :
que reconocemos como un problema de valor propio generalizado
Como es de rango unitario, solo tiene un valor propio distinto de cero. Se puede demostrar que este valor propio es igual a
obteniendo el siguiente filtro adaptado óptimo
Este es el mismo resultado encontrado en la subsección anterior.
Interpretación como estimador de mínimos cuadrados
Derivación
El filtrado emparejado también se puede interpretar como un estimador de mínimos cuadrados para la ubicación y escala óptimas de un modelo o plantilla determinados. Una vez más, definamos la secuencia observada como
donde es ruido de media cero no correlacionado. Se supone que la señal es una versión escalada y desplazada de una secuencia de modelo conocida :
Queremos encontrar estimaciones óptimas para el cambio y escala desconocidos minimizando los residuos de mínimos cuadrados entre la secuencia observada y una "secuencia de sondeo" :
El filtro apropiado se convertirá más tarde en el filtro adecuado, pero aún no se ha especificado. Al expandir y elevar al cuadrado la suma se obtiene
El primer término entre paréntesis es una constante (ya que se da la señal observada) y no tiene influencia en la solución óptima. El último término tiene un valor esperado constante porque el ruido no está correlacionado y tiene media cero. Por lo tanto, podemos eliminar ambos términos de la optimización. Después de invertir el signo, obtenemos el problema de optimización equivalente.
Al establecer la derivada con respecto a cero se obtiene una solución analítica para :
Insertar esto en nuestra función objetivo produce un problema de maximización reducido para solo :
El numerador se puede acotar superiormente mediante la desigualdad de Cauchy-Schwarz :
El problema de optimización supone su máximo cuando se cumple la igualdad en esta expresión. Según las propiedades de la desigualdad de Cauchy-Schwarz, esto sólo es posible cuando
para constantes arbitrarias distintas de cero o , y la solución óptima se obtiene en como se desea. Por lo tanto, nuestra "secuencia de sondeo" debe ser proporcional al modelo de señal , y la elección conveniente produce el filtro adaptado
Tenga en cuenta que el filtro es el modelo de señal reflejado. Esto garantiza que la operación que se debe aplicar para encontrar el óptimo es, de hecho, la convolución entre la secuencia observada y el filtro adaptado . La secuencia filtrada asume su máximo en la posición donde la secuencia observada coincide mejor (en un sentido de mínimos cuadrados) con el modelo de señal .
Trascendencia
El filtro adaptado se puede derivar de diversas maneras, [2] pero como un caso especial de un procedimiento de mínimos cuadrados también se puede interpretar como un método de máxima verosimilitud en el contexto de un modelo de ruido gaussiano (coloreado) y la verosimilitud de Whittle asociada . [5] Si la señal transmitida no poseyera parámetros desconocidos (como tiempo de llegada, amplitud,...), entonces el filtro adaptado minimizaría, según el lema de Neyman-Pearson , la probabilidad de error. Sin embargo, dado que la señal exacta generalmente está determinada por parámetros desconocidos que efectivamente se estiman (o ajustan ) en el proceso de filtrado, el filtro adaptado constituye una estadística (de prueba) de máxima verosimilitud generalizada . [6] La serie temporal filtrada se puede interpretar entonces como (proporcional a) la verosimilitud del perfil , la verosimilitud condicional maximizada como una función del parámetro de tiempo. [7] Esto implica en particular que la probabilidad de error (en el sentido de Neyman y Pearson, es decir, en relación con la maximización de la probabilidad de detección para una probabilidad de falsa alarma dada [8] ) no es necesariamente óptima. Lo que comúnmente se conoce como la relación señal-ruido (SNR) , que se supone que se maximiza mediante un filtro adaptado, en este contexto corresponde a , donde es la razón de verosimilitud maximizada (condicionalmente). [7] [nb 1]
La construcción del filtro adaptado se basa en un espectro de ruido conocido . Sin embargo, en la realidad, el espectro de ruido suele estimarse a partir de datos y, por lo tanto, solo se conoce con una precisión limitada. En el caso de un espectro incierto, el filtro adaptado se puede generalizar a un procedimiento iterativo más robusto con propiedades favorables también en ruido no gaussiano. [7]
Interpretación en el dominio de la frecuencia
Cuando se observa en el dominio de frecuencia, es evidente que el filtro adaptado aplica la mayor ponderación a los componentes espectrales que exhiben la mayor relación señal-ruido (es decir, gran peso donde el ruido es relativamente bajo, y viceversa). En general, esto requiere una respuesta de frecuencia no plana, pero la "distorsión" asociada no es motivo de preocupación en situaciones como el radar y las comunicaciones digitales , donde se conoce la forma de onda original y el objetivo es la detección de esta señal contra el ruido de fondo. Desde el punto de vista técnico, el filtro adaptado es un método de mínimos cuadrados ponderado basado en los datos del dominio de frecuencia ( heteroscedásticos ) (donde los "pesos" se determinan a través del espectro de ruido, consulte también la sección anterior), o equivalentemente, un método de mínimos cuadrados aplicado a los datos blanqueados .
Ejemplos
Radar y sonar
Los filtros adaptados se utilizan a menudo en la detección de señales . [1] Como ejemplo, supongamos que deseamos calcular la distancia de un objeto reflejando una señal de él. Podemos optar por transmitir una sinusoide de tono puro a 1 Hz. Suponemos que nuestra señal recibida es una forma atenuada y desfasada de la señal transmitida con ruido añadido.
Para juzgar la distancia del objeto, correlacionamos la señal recibida con un filtro adaptado, que, en el caso del ruido blanco (no correlacionado) , es otra sinusoide de 1 Hz de tono puro. Cuando la salida del sistema de filtro adaptado supera un cierto umbral, concluimos con alta probabilidad que la señal recibida se ha reflejado en el objeto. Usando la velocidad de propagación y el momento en que observamos por primera vez la señal reflejada, podemos estimar la distancia del objeto. Si cambiamos la forma del pulso de una manera especialmente diseñada, la relación señal-ruido y la resolución de la distancia pueden incluso mejorarse después del filtrado adaptado: esta es una técnica conocida como compresión de pulsos .
Además, los filtros adaptados se pueden utilizar en problemas de estimación de parámetros (véase teoría de estimación ). Para volver a nuestro ejemplo anterior, podemos desear estimar la velocidad del objeto, además de su posición. Para explotar el efecto Doppler , nos gustaría estimar la frecuencia de la señal recibida. Para ello, podemos correlacionar la señal recibida con varios filtros adaptados de sinusoides a frecuencias variables. El filtro adaptado con la salida más alta revelará, con alta probabilidad, la frecuencia de la señal reflejada y nos ayudará a determinar la velocidad radial del objeto, es decir, la velocidad relativa ya sea directamente hacia o lejos del observador. Este método es, de hecho, una versión simple de la transformada discreta de Fourier (DFT) . La DFT toma una entrada compleja de valor y la correlaciona con filtros adaptados, correspondientes a exponenciales complejos a diferentes frecuencias, para producir números de valor complejo correspondientes a las amplitudes y fases relativas de los componentes sinusoidales (véase Indicación de objetivo móvil ).
Comunicaciones digitales
El filtro adaptado también se utiliza en comunicaciones. En el contexto de un sistema de comunicación que envía mensajes binarios desde el transmisor al receptor a través de un canal ruidoso, se puede utilizar un filtro adaptado para detectar los pulsos transmitidos en la señal recibida ruidosa.
Imaginemos que queremos enviar la secuencia “0101100100” codificada en no polar sin retorno a cero (NRZ) a través de un canal determinado.
Matemáticamente, una secuencia en código NRZ se puede describir como una secuencia de pulsos unitarios o funciones rect desplazadas , cada pulso se pondera en +1 si el bit es "1" y en -1 si el bit es "0". Formalmente, el factor de escala para el bit es,
Podemos representar nuestro mensaje, , como la suma de pulsos unitarios desplazados:
donde es la longitud de tiempo de un bit y es la función rectangular .
Por lo tanto, la señal que debe enviar el transmisor es
Si modelamos nuestro canal ruidoso como un canal AWGN , se agrega ruido gaussiano blanco a la señal. En el extremo del receptor, para una relación señal-ruido de 3 dB, esto podría verse así:
A primera vista no se puede ver la secuencia original transmitida. Hay una alta potencia de ruido en relación con la potencia de la señal deseada (es decir, hay una baja relación señal-ruido ). Si el receptor muestreara esta señal en los momentos correctos, el mensaje binario resultante podría ser incorrecto.
Para aumentar nuestra relación señal-ruido, pasamos la señal recibida a través de un filtro adaptado. En este caso, el filtro debe adaptarse a un pulso NRZ (equivalente a un "1" codificado en código NRZ). Precisamente, la respuesta al impulso del filtro adaptado ideal, suponiendo que hay ruido blanco (no correlacionado), debe ser una versión escalada compleja conjugada en el tiempo invertida de la señal que estamos buscando. Elegimos
En este caso, debido a la simetría, el conjugado complejo invertido en el tiempo de es de hecho , lo que nos permite llamar a la respuesta al impulso de nuestro sistema de convolución de filtro adaptado.
Después de convolucionar con el filtro adaptado correctamente, la señal resultante es,
donde denota convolución.
Que ahora el receptor puede muestrear de forma segura en los instantes de muestreo correctos y compararlo con un umbral apropiado, dando como resultado una interpretación correcta del mensaje binario.
Astronomía de ondas gravitacionales
Los filtros adaptados desempeñan un papel central en la astronomía de ondas gravitacionales . [9] La primera observación de ondas gravitacionales se basó en el filtrado a gran escala de la salida de cada detector para señales que se asemejaran a la forma esperada, seguido de una detección posterior de disparadores coincidentes y coherentes entre ambos instrumentos. [10] Las tasas de falsas alarmas y, con ello, la significación estadística de la detección se evaluaron utilizando métodos de remuestreo . [11] [12] La inferencia sobre los parámetros de la fuente astrofísica se completó utilizando métodos bayesianos basados en modelos teóricos parametrizados para la forma de onda de la señal y (de nuevo) en la probabilidad de Whittle . [13] [14]
Sismología
Los filtros adaptados se utilizan en sismología para detectar terremotos similares u otras señales sísmicas, a menudo utilizando plantillas determinadas empíricamente de múltiples componentes y/o múltiples canales. [15] Las aplicaciones de filtrado adaptado en sismología incluyen la generación de grandes catálogos de eventos para estudiar la sismicidad de los terremotos [16] y la actividad volcánica, [17] [18] y en la detección global de explosiones nucleares. [19]
Biología
Los animales que viven en entornos relativamente estáticos tendrían que percibir características relativamente fijas del entorno, lo que permite la evolución de filtros que hacen coincidir la señal esperada con la mayor relación señal-ruido, el filtro adaptado. [20] Los sensores que perciben el mundo "a través de un 'filtro adaptado' de este tipo limitan severamente la cantidad de información que el cerebro puede captar del mundo exterior, pero lo liberan de la necesidad de realizar cálculos más complejos para extraer la información que finalmente se necesita para cumplir una tarea en particular". [21]
Véase también
- Periodograma
- Retroproyección filtrada (transformación de Radon)
- Filtro digital
- Procesamiento estadístico de señales
- Reducir la probabilidad
- Probabilidad de perfil
- Teoría de la detección
- Problema de comparaciones múltiples
- Capacidad del canal
- Teorema de codificación de canal ruidoso
- Estimación de la densidad espectral
- Filtro de mínimos cuadrados medios (LMS)
- Filtro de Viena
- Clasificación de múltiples señales (MUSIC), un método de superresolución paramétrica popular
- SAMV
Notas
- ^ La referencia común a la relación señal/ruido ha sido criticada por ser un tanto engañosa: " La característica interesante de este enfoque es que se alcanza la perfección teórica sin apuntar conscientemente a una relación señal/ruido máxima. Como asunto de interés bastante incidental, sucede que la operación [...] maximiza la relación señal/ruido pico, pero este hecho no juega ningún papel en la teoría actual. La relación señal/ruido no es una medida de información [...]. " ( Woodward , 1953; [1] Sec.5.1).
Referencias
- ^ abc Woodward, PM (1953). Probabilidad y teoría de la información con aplicaciones al radar . Londres: Pergamon Press .
- ^ ab Turin, GL (1960). "Introducción a los filtros adaptados". IRE Transactions on Information Theory . 6 (3): 311–329. doi :10.1109/TIT.1960.1057571. S2CID 5128742.
- ^ "Demodulación". OpenStax CNX . Consultado el 18 de abril de 2017 .
- ^ En honor a DO North, que fue uno de los primeros en introducir el concepto: North, DO (1943). "An analysis of the factors which determine signal/noise discrimination in pulsed carrier systems" (Análisis de los factores que determinan la discriminación señal/ruido en sistemas de portadora pulsada). Informe PPR-6C, RCA Laboratories, Princeton, NJ .
Reimpresión: North, DO (1963). "Análisis de los factores que determinan la discriminación señal/ruido en sistemas de portadora pulsada". Actas del IEEE . 51 (7): 1016–1027. doi :10.1109/PROC.1963.2383.
Véase también: Jaynes, ET (2003). "14.6.1 El filtro clásico emparejado ". Teoría de la probabilidad: la lógica de la ciencia . Cambridge: Cambridge University Press . - ^ Choudhuri, N.; Ghosal, S.; Roy, A. (2004). "Contigüidad de la medida de Whittle para una serie temporal gaussiana". Biometrika . 91 (4): 211–218. doi : 10.1093/biomet/91.1.211 .
- ^ Mood, AM; Graybill, FA; Boes, DC (1974). "IX. Pruebas de hipótesis ". Introducción a la teoría de la estadística (3.ª ed.). Nueva York: McGraw-Hill.
- ^ abc Röver, C. (2011). "Filtro basado en Student-t para detección robusta de señales". Physical Review D . 84 (12): 122004. arXiv : 1109.0442 . Código Bibliográfico :2011PhRvD..84l2004R. doi :10.1103/PhysRevD.84.122004.
- ^ Neyman, J.; Pearson, ES (1933). "Sobre el problema de las pruebas más eficientes de hipótesis estadísticas". Philosophical Transactions of the Royal Society of London A . 231 (694–706): 289–337. Bibcode :1933RSPTA.231..289N. doi : 10.1098/rsta.1933.0009 .
- ^ Schutz, BF (1999). "Astronomía de ondas gravitacionales". Gravedad clásica y cuántica . 16 (12A): A131–A156. arXiv : gr-qc/9911034 . Código Bibliográfico :1999CQGra..16A.131S. doi :10.1088/0264-9381/16/12A/307. S2CID 19021009.
- ^ "LIGO: Cómo buscamos agujeros negros fusionados y encontramos GW150914".
Se utiliza una técnica conocida como filtrado adaptado para ver si hay señales contenidas en nuestros datos. El objetivo del filtrado adaptado es ver si los datos contienen señales similares a un miembro del banco de plantillas. Dado que nuestras plantillas deben describir las formas de onda gravitacionales para el rango de diferentes sistemas fusionados que esperamos poder ver, cualquier señal lo suficientemente fuerte debería encontrarse mediante este método.
- ^ Usman, Samantha A. (2016). "La búsqueda de ondas gravitacionales de la coalescencia binaria compacta en PyCBC". Clase. Quantum Grav . 33 (21): 215004. arXiv : 1508.02357 . Código Bibliográfico :2016CQGra..33u5004U. doi :10.1088/0264-9381/33/21/215004. S2CID: 53979477.
- ^ Abbott, BP; et al. (La Colaboración Científica LIGO, la Colaboración Virgo) (2016). "GW150914: Primeros resultados de la búsqueda de coalescencia de agujeros negros binarios con Advanced LIGO". Physical Review D . 93 (12): 122003. arXiv : 1602.03839 . Bibcode :2016PhRvD..93l2003A. doi :10.1103/PhysRevD.93.122003. PMC 7430253 . PMID 32818163.
- ^ Abbott, BP; et al. (La Colaboración Científica LIGO, la Colaboración Virgo) (2016). "Propiedades de la fusión de agujeros negros binarios GW150914". Physical Review Letters . 116 (24): 241102. arXiv : 1602.03840 . Código Bibliográfico :2016PhRvL.116x1102A. doi :10.1103/PhysRevLett.116.241102. PMID 27367378. S2CID 217406416.
- ^ Meyer, R.; Christensen, N. (2016). "Ondas gravitacionales: una autopsia estadística de la fusión de un agujero negro". Significance . 13 (2): 20–25. doi : 10.1111/j.1740-9713.2016.00896.x .
- ^ Senobari, N.; Funning, G.; Keogh, E.; Zhu, Y.; Yah, C.; Zimmerman, Z.; Mueen, A. (2018). "Correlación cruzada supereficiente (SEC-C): un código de filtrado rápido y adaptado adecuado para computadoras de escritorio". Seismological Research Letters . 90 (1).
- ^ Shelly, D. (2017). "Un catálogo de 15 años de más de 1 millón de terremotos de baja frecuencia: seguimiento de temblores y deslizamientos a lo largo de la falla profunda de San Andrés". Revista de investigación geofísica . 122 (5).
- ^ Shelly, D.; Thelen, W. (2019). "Anatomía del colapso de una caldera: secuencia de sismicidad en la cumbre del Kilauea 2018 en alta resolución". Geophysical Research Letters . 46 (24).
- ^ Knox, H.; Chaput, J.; Aster, R.; Kyle, P. (2018). "Cambios plurianuales en conductos poco profundos observados con sismogramas de erupción de lagos de lava en el volcán Erebus, Antártida". Journal of Geophysical Research: Solid Earth . 123 .
- ^ Robinson, E. (1963). "Desarrollo matemático de filtros discretos para la detección de explosiones nucleares". Journal of Geophysical Research . 68 (19).
- ^ Warrant, Eric J. (octubre de 2016). "Filtros sensoriales adaptados". Current Biology . 26 (20): R976–R980. doi : 10.1016/j.cub.2016.05.042 . ISSN 0960-9822. PMID 27780072.
- ^ Wehner, Rüdiger (1987). «'Filtros emparejados': modelos neuronales del mundo externo». Journal of Comparative Physiology A . 161 (4): 511–531. doi :10.1007/bf00603659. ISSN 0340-7594. S2CID 32779686.
Lectura adicional
- Turín, GL (1960). "Introducción a los filtros adaptados". IRE Transactions on Information Theory . 6 (3): 311–329. doi :10.1109/TIT.1960.1057571. S2CID 5128742.
- Wainstein, LA; Zubakov, VD (1962). Extracción de señales a partir del ruido . Englewood Cliffs, NJ: Prentice-Hall .
- Melvin, WL (2004). "Una visión general de STAP". Revista IEEE Aerospace and Electronic Systems . 19 (1): 19–35. doi :10.1109/MAES.2004.1263229. S2CID 31133715.
- Röver, C. (2011). "Filtro basado en Student-t para detección robusta de señales". Physical Review D . 84 (12): 122004. arXiv : 1109.0442 . Código Bibliográfico :2011PhRvD..84l2004R. doi :10.1103/PhysRevD.84.122004.
- Fish, A.; Gurevich, S.; Hadani, R.; Sayeed, A.; Schwartz, O. (diciembre de 2011). "Cálculo del filtro adaptado en tiempo lineal". arXiv : 1112.4883 [cs.IT].