Filtro Gaussiano Medio Móvil

El promedio móvil como filtro El promedio móvil se utiliza a menudo para suavizar los datos en presencia de ruido. La media móvil simple no siempre se reconoce como el filtro de respuesta de impulso finito (FIR) que es, mientras que es realmente uno de los filtros más comunes en el procesamiento de señales. Tratarlo como un filtro permite compararlo con, por ejemplo, filtros de sinc de ventana (véanse los artículos sobre filtros de paso bajo, paso alto, paso de banda y rechazo de banda para ejemplos de los mismos). La principal diferencia con estos filtros es que el promedio móvil es adecuado para señales para las cuales la información útil está contenida en el dominio del tiempo. De las cuales las mediciones de suavizado por promediado son un excelente ejemplo. Sin embargo, los filtros windowed-sinc son fuertes en el dominio de la frecuencia. Con la ecualización en el procesamiento de audio como un ejemplo típico. Hay una comparación más detallada de ambos tipos de filtros en el dominio del tiempo frente al rendimiento de los dominios de frecuencia de los filtros. Si tiene datos para los que tanto el tiempo como el dominio de frecuencia son importantes, entonces puede que desee echar un vistazo a Variaciones en el promedio móvil. Que presenta una serie de versiones ponderadas de la media móvil que son mejores en eso. El promedio móvil de longitud (N) puede definirse como escrito tal como se implementa típicamente, con la muestra de salida actual como el promedio de las muestras (N) anteriores. Visto como un filtro, el promedio móvil realiza una convolución de la secuencia de entrada (xn) con un pulso rectangular de longitud (N) y altura (1 / N) (para hacer el área del pulso y, por tanto, la ganancia de El filtro, uno). En la práctica, es mejor tomar (N) impar. Aunque un promedio móvil también puede calcularse usando un número par de muestras, usar un valor impar para (N) tiene la ventaja de que el retardo del filtro será un número entero de muestras, ya que el retardo de un filtro con (N) Muestras es exactamente ((N-1) / 2). El promedio móvil puede entonces alinearse exactamente con los datos originales desplazándolo por un número entero de muestras. Dominio de tiempo Dado que el promedio móvil es una convolución con un pulso rectangular, su respuesta de frecuencia es una función sinc. Esto hace que sea algo así como el dual del filtro windowed-sinc, ya que es una convolución con un pulso sinc que da como resultado una respuesta de frecuencia rectangular. Es esta respuesta de frecuencia de sinc que hace que el promedio móvil sea un pobre intérprete en el dominio de la frecuencia. Sin embargo, funciona muy bien en el dominio del tiempo. Por lo tanto, es perfecto para suavizar los datos para eliminar el ruido, mientras que al mismo tiempo sigue manteniendo una respuesta de paso rápido (Figura 1). Para el ruido gaussiano blanco aditivo típico (AWGN) que se asume a menudo, las muestras del promedio (N) tienen el efecto de aumentar el SNR por un factor de (sqrt N). Dado que el ruido para las muestras individuales no está correlacionado, no hay razón para tratar cada muestra de manera diferente. Por lo tanto, el promedio móvil, que da a cada muestra el mismo peso, eliminará la cantidad máxima de ruido para una nitidez de respuesta dada. Implementación Debido a que es un filtro FIR, el promedio móvil puede implementarse a través de la convolución. Entonces tendrá la misma eficiencia (o falta de ella) como cualquier otro filtro FIR. Sin embargo, también se puede implementar recursivamente, de una manera muy eficiente. Se deduce directamente de la definición que esta fórmula es el resultado de las expresiones para (yn) y (yn1), es decir, donde observamos que el cambio entre (yn1) y (yn) es que un término extra (xn1 / N) Aparece al final, mientras que el término (xn-N1 / N) se elimina desde el principio. En aplicaciones prácticas, a menudo es posible omitir la división por (N) para cada término, compensando la ganancia resultante de (N) en otro lugar. Esta implementación recursiva será mucho más rápida que la convolución. Cada nuevo valor de (y) se puede calcular con sólo dos adiciones, en lugar de las (N) adiciones que serían necesarias para una implementación directa de la definición. Una cosa a tener en cuenta con una implementación recursiva es que se acumularán errores de redondeo. Esto puede o no ser un problema para su aplicación, pero también implica que esta implementación recursiva funcionará mejor con una implementación entera que con números de coma flotante. Esto es bastante inusual, ya que una implementación en coma flotante suele ser más simple. La conclusión de todo esto debe ser que usted nunca debe subestimar la utilidad del filtro de media móvil simple en aplicaciones de procesamiento de señales. Herramienta de diseño de filtros Este artículo se complementa con una herramienta de diseño de filtros. Experimente con diferentes valores para (N) y visualice los filtros resultantes. Pruébelo ahora El científico y los ingenieros dirigen el procesamiento de señales digitales Por Steven W. Smith, Ph. D. En un mundo perfecto, los diseñadores de filtros sólo tendrían que ocuparse de la información codificada en el dominio del tiempo o en el dominio de la frecuencia, pero nunca una mezcla de los dos en la misma señal. Desafortunadamente, hay algunas aplicaciones donde ambos dominios son simultáneamente importantes. Por ejemplo, las señales de televisión caen en esta categoría desagradable. La información de vídeo se codifica en el dominio de tiempo, es decir, la forma de la forma de onda corresponde a los patrones de brillo en la imagen. Sin embargo, durante la transmisión, la señal de vídeo se trata según su composición de frecuencia, tal como su anchura de banda total, cómo se añaden las ondas portadoras para el sonido y el color del amplificador, restauración de amplificación de eliminación del componente de CC, etc. Se entiende mejor en el dominio de la frecuencia, incluso si la información de las señales está codificada en el dominio del tiempo. Por ejemplo, el monitor de temperatura en un experimento científico podría estar contaminado con 60 hercios de las líneas eléctricas, 30 kHz de una fuente de alimentación de conmutación, o 1320 kHz de una emisora ​​local de radio AM. Los familiares del filtro de media móvil tienen un mejor rendimiento en el dominio de la frecuencia, y pueden ser útiles en estas aplicaciones de dominio mixto. Los filtros de media móvil de paso múltiple implican pasar la señal de entrada a través de un filtro de media móvil dos o más veces. La figura 15-3a muestra el núcleo del filtro resultante de una, dos y cuatro pasadas. Dos pasadas son equivalentes a usar un núcleo de filtro triangular (un núcleo de filtro rectangular convolucionado con sí mismo). Después de cuatro o más pases, el kernel de filtro equivalente parece un Gaussiano (recuerde el Teorema del Límite Central). Como se muestra en (b), múltiples pasadas producen una respuesta de paso en forma de s, en comparación con la línea recta de la única pasada. Las respuestas de frecuencia en (c) y (d) están dadas por la Ec. 15-2 multiplicado por sí mismo para cada pase. Es decir, cada vez que la convolución del dominio da como resultado una multiplicación de los espectros de frecuencia. La figura 15-4 muestra la respuesta en frecuencia de otros dos familiares del filtro de media móvil. Cuando un Gaussiano puro es usado como un núcleo de filtro, la respuesta de frecuencia es también Gaussiana, como se discutió en el Capítulo 11. El Gaussiano es importante porque es la respuesta de impulso de muchos sistemas naturales y artificiales. Por ejemplo, un breve impulso de luz que entra en una línea de transmisión de fibra óptica larga saldrá como un pulso gaussiano, debido a las diferentes trayectorias tomadas por los fotones dentro de la fibra. El kernel de filtro gaussiano también se utiliza ampliamente en el procesamiento de imágenes porque tiene propiedades únicas que permiten convoluciones bidimensionales rápidas (véase el Capítulo 24). La segunda respuesta de frecuencia en la Fig. 15-4 corresponde a usar una ventana de Blackman como un núcleo de filtro. (El término ventana no tiene significado aquí es simplemente parte del nombre aceptado de esta curva). La forma exacta de la ventana de Blackman se da en el Capítulo 16 (Ec. 16-2, Fig. 16-2) sin embargo, se parece mucho a un Gaussiano. ¿Cómo son estos parientes del filtro de media móvil mejor que el filtro de media móvil en sí Tres maneras: En primer lugar, y lo más importante, estos filtros tienen mejor atenuación de banda de detención que el filtro de media móvil. En segundo lugar, los granos de filtro se estrechan hasta una amplitud más pequeña cerca de los extremos. Recuerde que cada punto en la señal de salida es una suma ponderada de un grupo de muestras de la entrada. Si el núcleo del filtro se estrecha, las muestras en la señal de entrada que están más alejadas reciben menos peso que las cercanas. En tercer lugar, las respuestas de paso son curvas suaves, en lugar de la línea recta brusca de la media móvil. Estos últimos dos son generalmente de beneficio limitado, aunque usted puede ser que encuentre aplicaciones donde son ventajas genuinas. El filtro de media móvil y sus familiares son todos aproximadamente iguales en la reducción del ruido aleatorio mientras que mantiene una respuesta aguda del paso. La ambigüedad radica en cómo se mide el tiempo de subida de la respuesta escalonada. Si el tiempo de subida se mide de 0 a 100 del paso, el filtro de media móvil es lo mejor que puede hacer, como se mostró anteriormente. En comparación, medir el tiempo de subida de 10 a 90 hace que la ventana de Blackman sea mejor que el filtro de media móvil. El punto es, esto es sólo disputas teóricas considerar estos filtros iguales en este parámetro. La mayor diferencia en estos filtros es la velocidad de ejecución. Utilizando un algoritmo recursivo (descrito a continuación), el filtro de media móvil funcionará como un rayo en su computadora. De hecho, es el filtro digital más rápido disponible. Múltiples pases del promedio móvil serán correspondientemente más lentos, pero aún así muy rápidos. En comparación, los filtros Gaussiano y Blackman son extremadamente lentos, porque deben usar convolución. Piense un factor de diez veces el número de puntos en el núcleo del filtro (basado en la multiplicación es aproximadamente 10 veces más lento que la adición). Por ejemplo, espere que un Gaussiano de 100 puntos sea 1000 veces más lento que un promedio móvil usando recursión. Considere tener una señal en el dominio del tiempo, y desea suavizar la señal. Media móvil y filtros gaussianos que se utilizan. ¿Cómo elegir qué se utiliza para qué? ¿Cuáles son las condiciones bajo las cuales Gaussian es mejor y las condiciones en las que media móvil es mejor Lo que estoy tratando de hacer con esta señal es, detección de pico inicialmente, luego aplicar pequeñas ventanas en cada parte y figura Los cambios de frecuencia (desplazamientos Doppler) para cada parte para determinar la dirección del movimiento a partir del cambio de frecuencia. Quiero suavizar la señal en el dominio del tiempo sin pérdida de información en el dominio de la frecuencia. Pensé para la parte de calcular hacia fuera los cambios de Doppler, usando STFT sería una idea buena. Si la referencia pudo ser dada a un cierto papel, eso también sería realmente provechoso. Preguntó Sep 29 a las 11:12