En matemáticas y análisis de Fourier , una transformada de Fourier de tiempo corto con máscara rectangular (rec-STFT) tiene la forma simple de transformada de Fourier de tiempo corto . Otros tipos de STFT pueden requerir más tiempo de cálculo 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
el ( a ) = { 1 ; | a | ≤ B 0 ; | a | > B {\displaystyle w(t)={\begin{cases}\ 1;&|t|\leq B\\\ 0;&|t|>B\end{cases}}} B = 50, eje x (seg)Podemos cambiar B para obtener diferentes compensaciones entre la resolución de tiempo deseada y la resolución de frecuencia.
Rec-STFT
incógnita ( a , F ) = ∫ a − B a + B incógnita ( τ ) mi − yo 2 π F τ d τ {\displaystyle X(t,f)=\int _{tB}^{t+B}x(\tau )e^{-j2\pi f\tau }\,d\tau } Forma inversa
incógnita ( a ) = ∫ − ∞ ∞ incógnita ( a 1 , F ) mi yo 2 π F a d F dónde a − B < a 1 < a + 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
(a)
∫ − ∞ ∞ incógnita ( a , F ) d F = ∫ a − B a + B incógnita ( τ ) ∫ − ∞ ∞ mi − yo 2 π F τ d F d τ = ∫ a − B a + B incógnita ( τ ) del ( τ ) d τ = { incógnita ( 0 ) ; | a | < 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{de lo contrario}}\end{cases}}} (b)
∫ − ∞ ∞ incógnita ( a , F ) mi − yo 2 π F en d F = { incógnita ( en ) ; en − B < a < en + B 0 ; de lo contrario {\displaystyle \int _{-\infty }^{\infty }X(t,f)e^{-j2\pi fv}\,df={\begin{casos}\ x(v);&v-B<t<v+B\\\ 0;&{\text{en caso contrario}}\end{casos}}} Desplazamiento de propiedad (desplazamiento a lo largo del eje x) ∫ a − B a + B incógnita ( τ + τ 0 ) mi − yo 2 π F τ d τ = incógnita ( a + τ 0 , F ) mi yo 2 π 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 ) ∫ a − B a + B [ incógnita ( τ ) mi yo 2 π F 0 τ ] d τ = incógnita ( a , F − F 0 ) {\displaystyle \int _{tB}^{t+B}[x(\tau )e^{j2\pi f_{0}\tau }]d\tau =X(t,f-f_{0})} Cuando incógnita ( a ) = del ( a ) , incógnita ( a , F ) = { 1 ; | a | < B 0 ; de lo contrario {\displaystyle x(t)=\delta (t),X(t,f)={\begin{cases}\ 1;&|t|<B\\\ 0;&{\text{de lo contrario}}\end{cases}}} Cuando incógnita ( a ) = 1 , incógnita ( a , F ) = 2 B Sincronización ( 2 B F ) mi yo 2 π F a {\displaystyle x(t)=1,X(t,f)=2B\nombre del operador {sinc} (2Bf)e^{j2\pi ft}} Si , y son sus rec-STFT, entonces yo ( a ) = alfa incógnita ( a ) + β y ( a ) {\displaystyle h(t)=\alpha x(t)+\beta y(t)\,} yo ( a , F ) , incógnita ( a , F ) , {\displaystyle H(t,f),X(t,f),} Y ( a , F ) {\displaystyle Y(t,f)\,}
yo ( a , F ) = alfa incógnita ( a , F ) + β Y ( a , F ) . {\displaystyle H(t,f)=\alpha X(t,f)+\beta Y(t,f).} Propiedad de integración de potencia ∫ − ∞ ∞ | incógnita ( a , F ) | 2 d F = ∫ a − B a + B | incógnita ( τ ) | 2 d τ {\displaystyle \int _{-\infty }^{\infty }|X(t,f)|^{2}\,df=\int _{tB}^{t+B}|x(\tau )|^{2}\,d\tau } ∫ − ∞ ∞ ∫ − ∞ ∞ | incógnita ( a , F ) | 2 d F d a = 2 B ∫ − ∞ ∞ | incógnita ( τ ) | 2 d τ {\displaystyle \int _{-\infty }^{\infty }\int _{-\infty }^{\infty }|X(t,f)|^{2}\,df\,dt=2B\int _{-\infty }^{\infty }|x(\tau )|^{2}\,d\tau } ∫ − ∞ ∞ incógnita ( a , F ) Y ∗ ( a , F ) d F = ∫ a − B a + B incó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 ( a , F ) Y ∗ ( a , F ) d F d a = 2 B ∫ − ∞ ∞ incó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 producidos a partir de la aplicación de una rec-STFT en una función que consta de 3 ondas coseno consecutivas. (el espectrograma superior utiliza un B más pequeño de 0,5, el del medio utiliza un B de 1 y el inferior utiliza un B más grande de 2).En la imagen, cuando B es más pequeño, la resolución temporal es mejor. De lo contrario, cuando B es más grande, 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 de cálculo.En comparación con otros tipos de análisis de 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 de tiempo-frecuencia. La discontinuidad de salto de los bordes de la máscara rectangular produce artefactos de timbre de Gibbs en el dominio de la frecuencia, que se pueden aliviar con ventanas más suaves .
Véase también
Referencias Jian-Jiun Ding (2014) Análisis de tiempo-frecuencia y transformada wavelet