Deducción paso a paso de la Transformada
de Stockwell (ST)
Introducción
La transformada de Stockwell (ST) se desarrolló para analizar señales no estacionarias,
proporcionando información tanto de frecuencia como de tiempo, y conservando además
la fase, algo que no ocurre en otras transformadas como la wavelet. A continuación se
describe de manera formal cómo se construye su fórmula.
1. Transformada de Fourier (FT)
La transformada de Fourier permite obtener el contenido en frecuencia de una señal x(t):
Z +∞
X(f ) = x(t) e−j2πf t dt
−∞
Sin embargo, se pierde toda información temporal sobre cuándo ocurren los cambios en
la señal.
2. Transformada de Fourier de Tiempo Corto (STFT)
Para añadir localización temporal, se introduce una ventana w(τ −t) centrada en el tiempo
τ , obteniendo: Z+∞
SSTFT (τ, f ) = x(t) w(τ − t) e−j2πf t dt
−∞
El problema es que la ventana es siempre la misma y no se adapta a diferentes frecuencias,
lo que limita la resolución tiempo-frecuencia.
3. Idea clave de la ST
La transformada de Stockwell propone que la ventana dependa de la frecuencia:
w(τ − t, f )
1
Ası́, la ventana se ajusta automáticamente:
• A frecuencias altas, la ventana es más estrecha (mejor resolución temporal).
• A frecuencias bajas, la ventana es más ancha (mejor resolución frecuencial).
4. Definición formal de la ST
La transformada de Stockwell se define como:
Z +∞
S(τ, f ) = x(t) w(τ − t, f ) e−j2πf t dt
−∞
Esta expresión permite saber qué frecuencia hay (f ) y cuándo aparece (τ ), manteniendo
además la información de fase.
5. Ventana gaussiana dependiente de la frecuencia
Para realizar el análisis multirresolución, se emplea tı́picamente una ventana gaussiana
definida por:
|f | t2 f 2
w(t, f ) = √ e− 2
2π
Gracias a esta dependencia en f , la ventana se adapta de manera automática a cada
frecuencia.
6. Forma equivalente en el dominio de Fourier
La transformada de Stockwell puede expresarse también en términos de la transformada
de Fourier X(·) de la señal x(t):
Z +∞
2 α2 /f 2
S(τ, f ) = X(α + f ) e−2π ej2πατ dα
−∞
Esta forma es útil para implementaciones computacionales.
7. Versión discreta
Considerando una señal muestreada x[kT ] con N muestras y periodo de muestreo T , la
transformada discreta de Stockwell se expresa como:
N −1 2πmp
2 m2 /n2
S[pT, NnT ] = e−2π ej
X
m+n
X NT
N
m=0
donde p, n, m = 0, 1, . . . , N − 1.
2
8. Transformada inversa
La señal original puede recuperarse a partir de la transformada de Stockwell mediante:
N
1 NX
−1 −1 h i
−j
2πnk
x[kT ] =
X
n
N
S pT, N T
e N
n=0 p=0
Conclusión
La transformada de Stockwell parte de la idea de la STFT pero soluciona su limitación
mediante una ventana gaussiana dependiente de la frecuencia, lo que permite un análisis
multirresolución. Su fórmula:
Z +∞
S(τ, f ) = x(t) w(τ − t, f ) e−j2πf t dt
−∞
proporciona simultáneamente información sobre frecuencia, momento en que ocurre y
mantiene la fase de la señal original.