La transformada S como distribución tiempo-frecuencia se desarrolló en 1994 para analizar datos geofísicos. [ 1 ] [ 2 ] De esta manera, la transformada S es una generalización de la transformada de Fourier de tiempo corto (STFT), extendiendo la transformada wavelet continua y superando algunas de sus desventajas. Por un lado, las sinusoides de modulación están fijas con respecto al eje del tiempo; esto localiza las dilataciones y traslaciones de la ventana gaussiana escalable en la transformada S. Además, la transformada S no tiene un problema de término cruzado y produce una mejor claridad de señal que la transformada de Gabor . Sin embargo, la transformada S tiene sus propias desventajas: la claridad es peor que la función de distribución de Wigner y la función de distribución de clase de Cohen .
En 2010 se inventó un algoritmo de transformada S rápida. [ 3 ] [ 4 ] Este reduce la complejidad computacional de O[N 2 ·log(N)] a O[N·log(N)] y hace que la transformada sea biyectiva, donde la transformada tiene el mismo número de puntos que la señal o imagen de origen, en comparación con la complejidad de almacenamiento de N 2 para la formulación original. [ 4 ] [ 5 ] Una implementación está disponible para la comunidad de investigación bajo una licencia de código abierto . [ 6 ] [ 7 ]
Una formulación general de la transformada S [ 4 ] aclara la relación con otras transformadas tiempo-frecuencia como la transformada de Fourier, la transformada de Fourier de tiempo corto y la transformada wavelet. [ 4 ]
Definición
Hay varias formas de representar la idea de la transformada S. Aquí, la transformada S se deriva como la corrección de fase de la transformada wavelet continua con ventana siendo la función gaussiana .
- Transformada S
- Transformada S inversa
Forma modificada
- Forma del espectro
La definición anterior implica que la función de transformación s puede expresarse como la convolución dey. Aplicando la transformada de Fourier a ambosyda
- .
- Transformada S de tiempo discreto
A partir de la forma espectral de la transformada S, podemos derivar la transformada S de tiempo discreto. Sea, dóndees el intervalo de muestreo yes la frecuencia de muestreo. La transformada S de tiempo discreto se puede expresar entonces como:
Implementación de la transformada S de tiempo discreto
A continuación se muestra el pseudocódigo de la implementación.
Paso 1. Calcular bucle sobre m (voces) Paso 2. Calcularpara Paso 3. Movera Paso 4. Multiplica el Paso 2 y el Paso 3. Paso 5.IDFT(). Repetir.}
Comparación con otras herramientas de análisis tiempo-frecuencia
Comparación con la transformada de Gabor
La única diferencia entre la transformada de Gabor (GT) y la transformada S es el tamaño de la ventana. Para GT, el tamaño de la ventana es una función gaussiana.Mientras tanto, la función de ventana para la transformada S es una función de f. Con una función de ventana proporcional a la frecuencia, la transformada S funciona bien en el análisis del dominio de la frecuencia cuando la frecuencia de entrada es baja. Cuando la frecuencia de entrada es alta, la transformada S ofrece mayor claridad en el dominio del tiempo . Véase la tabla a continuación.
Este tipo de propiedad convierte a la transformada S en una herramienta poderosa para analizar el sonido, ya que el ser humano es sensible a la parte de baja frecuencia de una señal sonora.
Comparación con la transformada de Wigner
El principal problema de la transformada de Wigner reside en el término cruzado, que proviene de la función de autocorrelación. Este término cruzado puede generar ruido y distorsiones en los análisis de señales. Los análisis mediante la transformada S evitan este problema.
Comparación con la transformada de Fourier de corto tiempo
Podemos comparar la transformada S y la transformada de Fourier de tiempo corto (STFT). [ 2 ] [ 8 ] Primero, se utilizan una señal de alta frecuencia, una señal de baja frecuencia y una señal de ráfaga de alta frecuencia en el experimento para comparar el rendimiento. La característica de resolución dependiente de la frecuencia de la transformada S permite la detección de la ráfaga de alta frecuencia. Por otro lado, como la STFT consiste en un ancho de ventana constante, conduce a un resultado con una definición más pobre. En el segundo experimento, se agregan dos ráfagas más de alta frecuencia a chirps cruzados. En el resultado, las cuatro frecuencias fueron detectadas por la transformada S. Por otro lado, las dos ráfagas de alta frecuencia no son detectadas por la STFT. El término cruzado de las ráfagas de alta frecuencia hizo que la STFT tuviera una sola frecuencia en la frecuencia más baja.
Aplicaciones
- Filtrados de señal [ 9 ]
- Imágenes por resonancia magnética (IRM) [ 10 ]
- reconocimiento de perturbaciones en el sistema eléctrico
- Se ha demostrado que la transformada S puede identificar algunos tipos de perturbaciones, como caídas de tensión, aumentos de tensión, interrupciones momentáneas y transitorios oscilatorios. [ 11 ]
- La transformada S también se puede aplicar a otros tipos de perturbaciones, como muescas, armónicos con hundimientos y abultamientos, etc.
- La transformada S genera contornos adecuados para una simple inspección visual. Sin embargo, la transformada wavelet requiere herramientas específicas como el análisis multirresolución estándar .
- Análisis de señales geofísicas
Véase también
Referencias
- ↑ Stockwell, RG; Mansinha, L; Lowe, RP (1996). "Localización del espectro complejo: la transformada S". IEEE Transactions on Signal Processing . 44 (4): 998– 1001. Bibcode : 1996ITSP...44..998S . CiteSeerX 10.1.1.462.1500 . doi : 10.1109/78.492555 . S2CID 30202517 .
- 1 2 Stockwell, RG (1999). Análisis de la transformada S de la actividad de ondas gravitatorias a partir de una red a pequeña escala de generadores de imágenes de luminiscencia atmosférica. Tesis doctoral, Universidad de Western Ontario, London, Ontario, Canadá.
- ↑ Brown, RA; Frayne, R (2008). "Una transformada S discreta rápida para el procesamiento de señales biomédicas". 30.ª Conferencia Internacional Anual de la Sociedad de Ingeniería en Medicina y Biología del IEEE , 2008. Vol. 2008. págs. 2586–2589 . doi : 10.1109/IEMBS.2008.4649729 . ISBN 978-1-4244-1814-5. PMID 19163232 . S2CID 29974786 .
- 1 2 3 4 Brown, Robert A.; Lauzon, M. Louis; Frayne, Richard (enero de 2010). "Una descripción general de las transformadas tiempo-frecuencia lineales y formulación de una transformada rápida e invertible que muestrea el espectro de la transformada S continua de forma no redundante". IEEE Transactions on Signal Processing . 58 (1): 281– 290. Bibcode : 2010ITSP...58..281B . doi : 10.1109/tsp.2009.2028972 . ISSN 1053-587X . S2CID 16074001 .
- ↑ Kelly Sansom, "Transformación S rápida", Universidad de Calgary, https://www.ucalgary.ca/news/utoday/may31-2011/computing
- ↑ "Descarga rápida de S-Transform | SourceForge.net" . 13 de agosto de 2018.
- ↑ Brown, Robert (1 de diciembre de 2024), robb-brown/GFT , consultado el 10 de diciembre de 2024
- ↑ E. Sejdić, I. Djurović, J. Jiang, "Representación de características tiempo-frecuencia mediante concentración de energía: una visión general de los avances recientes", Procesamiento de señales digitales , vol. 19, n.º 1, págs. 153-183, enero de 2009.
- ↑ Ditommaso, Rocco; Mucciarelli, Marco; Ponzo, Felice Carlo (2012). "Análisis de sistemas estructurales no estacionarios mediante un filtro de banda variable" (PDF) . Boletín de Ingeniería Sísmica . 10 (3): 895– 911. Bibcode : 2012BuEE...10..895D . doi : 10.1007/s10518-012-9338-y .Véase también el archivo MATLAB .
- ↑ Hongmei Zhu y J. Ross Mitchell, "La transformada S en imágenes médicas", Centro de Investigación de RM de la Familia Seaman de la Universidad de Calgary, Centro Médico Foothills, Canadá.
- ↑ Ray, Prakash K.; Mohanty, Soumya R.; Kishor, Nand; Dubey, Harish C. (2010). "Determinación de la coherencia en un sistema híbrido basado en generación distribuida conectado a la red bajo escenarios de funcionamiento en isla". 2010 IEEE International Conference on Power and Energy . pp. 85–88 . doi : 10.1109/PECON.2010.5697562 . ISBN 978-1-4244-8947-3.
Lecturas adicionales
- Rocco Ditommaso, Felice Carlo Ponzo, Gianluca Auletta (2015). Detección de daños en estructuras reticuladas: evaluación de la curvatura modal mediante la transformada de Stockwell bajo excitación sísmica. Ingeniería sísmica y vibraciones en ingeniería. Junio de 2015, volumen 14, número 2, págs. 265-274.
- Rocco Ditommaso, Marco Mucciarelli, Felice C. Ponzo (2010). Filtro basado en la transformada S aplicado al análisis del comportamiento dinámico no lineal del suelo y los edificios. XIV Conferencia Europea de Ingeniería Sísmica. Actas. Ohrid, República de Macedonia. 30 de agosto - 3 de septiembre de 2010. (Disponible para descargar en http://roccoditommaso.xoom.it )
- M. Mucciarelli, M. Bianca, R. Ditommaso, MR Gallipoli, A. Masi, C Milkereit, S. Parolai, M. Picozzi, M. Vona (2011). DAÑOS EN CAMPO LEJANO EN EDIFICIOS RC: EL ESTUDIO DE CASO DE NAVELLI DURANTE LA SECUENCIA SÍSMICA DE L'AQUILA (ITALIA), 2009. Boletín de Ingeniería Sísmica . doi : 10.1007/s10518-010-9201-y .
- JJ Ding, "Apuntes del curso sobre análisis tiempo-frecuencia y transformada wavelet", Departamento de Ingeniería Eléctrica, Universidad Nacional de Taiwán (NTU), Taipéi, Taiwán, 2007.
- Jaya Bharata Reddy, Dusmanta Kumar Mohanta y BM Karan, "Reconocimiento de perturbaciones en sistemas de potencia mediante técnicas de ondículas y transformada S", Instituto Tecnológico Birla, Mesra, Ranchi-835215, 2004.
- B. Boashash, "Notas sobre el uso de la distribución de Wigner para el análisis de señales tiempo-frecuencia", IEEE Trans. on Acoust. Speech. and Signal Processing, vol. 26, n.º 9, 1987.
- RN Bracewell, La transformada de Fourier y sus aplicaciones, McGraw Hill Book Company, Nueva York, 1978
- EO Brigham, La transformada rápida de Fourier , Prentice-Hall Inc., Englewood Cliffs, Nueva Jersey, 1974
- Cohen, L. (1989). "Distribuciones tiempo-frecuencia: una revisión". Proc. IEEE . 77 (7): 941– 981. CiteSeerX 10.1.1.1026.2853 . doi : 10.1109/5.30749 .
- I. Daubechies, "La transformada wavelet, la localización tiempo-frecuencia y el análisis de señales", IEEE Trans. on Information Theory , vol. 36, n.º 5, septiembre de 1990.
- Farge, M. (1992). "Transformadas wavelet y su aplicación a la turbulencia" . Annual Review of Fluid Mechanics . 24 : 395–457 . doi : 10.1146/annurev.fluid.24.1.395 .
- D. Gabor, "Teoría de la comunicación", J. Inst. Elect. Eng., vol. 93, n.º 3, págs. 429–457, 1946
- Goupillaud, P.; Grossmann, A.; Morlet, J. (1984). "Transformadas ciclo-octava y relacionadas en el análisis sísmico". Geoexploration . 23 : 85–102 . doi : 10.1016/0016-7142(84)90025-5 .
- F. Hlawatsch y GF Boudreuax-Bartels, 1992 "Representaciones lineales y cuadráticas de señales tiempo-frecuencia", IEEE Signal Processing Magazine, págs. 21–67
- Rioul, O.; Vetterli, M. (1991). "Ondículas y procesamiento de señales" (PDF) . IEEE Signal Processing Magazine . 8 (4): 14– 38. Bibcode : 1991ISPM....8...14R . doi : 10.1109/79.91217 . S2CID 13266737 .
- RK Young, Teoría de las ondículas y sus aplicaciones, Kluwer Academic Publishers, Dordrecht, 1993
- Transformaciones integrales
- Análisis de Fourier
- Análisis tiempo-frecuencia