Articulo de referencia

Transformada de Fourier de corto tiempo con máscara rectangular

En matemáticas y análisis de Fourier , la transformada de Fourier de tiempo corto con máscara rectangular (rec-STFT) es una forma simplificada de la transformada de Fourier de t...

En matemáticas y análisis de Fourier , la transformada de Fourier de tiempo corto con máscara rectangular (rec-STFT) es una forma simplificada de la transformada de Fourier de tiempo corto que se utiliza para analizar cómo cambia el contenido de frecuencia de una señal a lo largo del tiempo. En la rec-STFT, se utiliza una ventana rectangular (una función simple de limitación de tiempo de encendido/apagado) para aislar segmentos cortos de la señal. Otros tipos de STFT pueden requerir más tiempo de cálculo (el tiempo que tarda un ordenador o algoritmo en realizar un cálculo o completar una tarea) que la rec-STFT.

La función de máscara rectangular se puede definir para algún límite ( B) en el tiempo ( t ) como

w(t)={ 1;|t|B 0;|t|>B{\displaystyle w(t)={\begin{cases}\ 1;&|t|\leq B\\\ 0;&|t|>B\end{cases}}}
B = 50, eje x (seg)

Podemos modificar B para obtener diferentes compensaciones entre la resolución temporal deseada y la resolución de frecuencia.

Rec-STFT

incógnita(t,F)=tBt+Bincógnita(τ)mij2πFτdτ{\displaystyle X(t,f)=\int _{tB}^{t+B}x(\tau )e^{-j2\pi f\tau }\,d\tau }

Forma inversa

incógnita(t)=incógnita(t1,F)mij2πFtdF dónde tB<t1<t+B{\displaystyle x(t)=\int _{-\infty }^{\infty }X(t_{1},f)e^{j2\pi ft}\,df{\text{ donde }}tB<t_{1}<t+B}

Propiedad

Rec-STFT tiene propiedades similares a la transformada de Fourier.

  • Integración

(a)

incógnita(t,F)dF=tBt+Bincógnita(τ)mij2πFτdFdτ=tBt+Bincógnita(τ)δ(τ)dτ={ incógnita(0);|t|<B 0;de lo contrario{\displaystyle \int _{-\infty }^{\infty }X(t,f)\,df=\int _{tB}^{t+B}x(\tau )\int _{-\infty }^{\infty }e^{-j2\pi f\tau }\,df\,d\tau =\int _{tB}^{t+B}x(\tau )\delta (\tau )\,d\tau ={\begin{cases}\ x(0);&|t|<B\\\ 0;&{\text{en otro caso}}\end{cases}}}

(b)

incógnita(t,F)mij2πFvdF={ incógnita(v);vB<t<v+B 0;de lo contrario{\displaystyle \int _{-\infty }^{\infty }X(t,f)e^{-j2\pi fv}\,df={\begin{cases}\ x(v);&v-B<t<v+B\\\ 0;&{\text{en otro caso}}\end{cases}}}
  • Propiedad de desplazamiento (desplazamiento a lo largo del eje x)
tBt+Bincógnita(τ+τ0)mij2πFτdτ=incógnita(t+τ0,F)mij2πFτ0{\displaystyle \int _{tB}^{t+B}x(\tau +\tau _{0})e^{-j2\pi f\tau }\,d\tau =X(t+\tau _{0},f)e^{j2\pi f\tau _{0}}}
  • Propiedad de modulación (desplazamiento a lo largo del eje y )
tBt+B[incógnita(τ)mij2πF0τ]dτ=incógnita(t,FF0){\displaystyle \int _{tB}^{t+B}[x(\tau )e^{j2\pi f_{0}\tau }]d\tau =X(t,f-f_{0})}
  • entrada especial
  1. Cuandoincógnita(t)=δ(t),incógnita(t,F)={ 1;|t|<B 0;de lo contrario{\displaystyle x(t)=\delta (t),X(t,f)={\begin{cases}\ 1;&|t|<B\\\ 0;&{\text{en otro caso}}\end{cases}}}
  2. Cuandoincógnita(t)=1,incógnita(t,F)=2Bdesde(2BF)mij2πFt{\displaystyle x(t)=1,X(t,f)=2B\operatorname {sinc} (2Bf)e^{j2\pi ft}}
  • Propiedad de linealidad

Sih(t)=αincógnita(t)+βy(t){\displaystyle h(t)=\alpha x(t)+\beta y(t)\,},H(t,F),incógnita(t,F),{\displaystyle H(t,f),X(t,f),}yY(t,F){\displaystyle Y(t,f)\,}son sus rec-STFT, entonces

H(t,F)=αincógnita(t,F)+βY(t,F).{\displaystyle H(t,f)=\alpha X(t,f)+\beta Y(t,f).}
  • Propiedad de integración de energía
|incógnita(t,F)|2dF=tBt+B|incógnita(τ)|2dτ{\displaystyle \int _{-\infty }^{\infty }|X(t,f)|^{2}\,df=\int _{tB}^{t+B}|x(\tau )|^{2}\,d\tau }
|incógnita(t,F)|2dFdt=2B|incógnita(τ)|2dτ{\displaystyle \int _{-\infty }^{\infty }\int _{-\infty }^{\infty }|X(t,f)|^{2}\,df\,dt=2B\int _{-\infty }^{\infty }|x(\tau )|^{2}\,d\tau }
incógnita(t,F)Y(t,F)dF=tBt+Bincógnita(τ)y(τ)dτ{\displaystyle \int _{-\infty }^{\infty }X(t,f)Y^{*}(t,f)\,df=\int _{tB}^{t+B}x(\tau )y^{*}(\tau )\,d\tau }
incógnita(t,F)Y(t,F)dFdt=2Bincógnita(τ)y(τ)dτ{\displaystyle \int _{-\infty }^{\infty }\int _{-\infty }^{\infty }X(t,f)Y^{*}(t,f)\,df\,dt=2B\int _{-\infty }^{\infty }x(\tau )y^{*}(\tau )\,d\tau }

Ejemplo de compensación con diferentes B

Espectrogramas generados al aplicar una STFT recta a una función compuesta por 3 ondas coseno consecutivas. (El espectrograma superior utiliza un valor de B menor de 0,5, el del medio un valor de B de 1 y el inferior un valor de B mayor de 2).

En la imagen, cuando B es menor, la resolución temporal es mejor. Por el contrario, cuando B es mayor, la resolución de frecuencia es mejor.

Ventajas y desventajas

En comparación con la transformada de Fourier:

  • Ventaja: Se puede observar la frecuencia instantánea.
  • Desventaja: Mayor complejidad computacional.

En comparación con otros tipos de análisis tiempo-frecuencia :

  • Ventaja: Menor tiempo de cálculo para la implementación digital.
  • Desventaja: La calidad es peor que la de otros tipos de análisis tiempo-frecuencia. La discontinuidad abrupta de los bordes de la máscara rectangular produce artefactos de resonancia de Gibbs en el dominio de la frecuencia, que pueden mitigarse con ventanas más suaves .

Véase también

Referencias

  1. Jian-Jiun Ding (2014) Análisis tiempo-frecuencia y transformada wavelet