Rectangular mask short-time Fourier transform

From HandWiki

{{multiple issues|

In mathematics and Fourier analysis, a rectangular mask short-time Fourier transform (rec-STFT) has the simple form of short-time Fourier transform. Other types of the STFT may require more computation time than the rec-STFT.

The rectangular mask function can be defined for some bound (B) over time (t) as

[math]\displaystyle{ w(t) =\begin{cases} \ 1; & |t|\leq B \\ \ 0; & |t|\gt B \end{cases} }[/math]
B = 50, x-axis (sec)

We can change B for different tradeoffs between desired time resolution and frequency resolution.

Rec-STFT

[math]\displaystyle{ X(t,f)=\int_{t-B}^{t+B} x(\tau) e^{-j2\pi f\tau} \, d\tau }[/math]

Inverse form

[math]\displaystyle{ x(t)=\int_{-\infty}^\infty X(t_1,f)e^{j2\pi ft} \, df\text{ where } t-B\lt t_1\lt t+B }[/math]

Property

Rec-STFT has similar properties with Fourier transform

  • Integration

(a)

[math]\displaystyle{ \int_{-\infty}^\infty X(t, f)\, df = \int_{t-B}^{t+B} x(\tau)\int_{-\infty}^\infty e^{-j 2 \pi f \tau}\, df \, d\tau = \int_{t-B}^{t+B} x(\tau)\delta(\tau) \, d\tau=\begin{cases} \ x(0); & |t|\lt B \\ \ 0; & \text{otherwise} \end{cases} }[/math]

(b)

[math]\displaystyle{ \int_{-\infty}^\infty X(t, f)e^{-j 2 \pi f v} \,df =\begin{cases} \ x(v); & v-B\lt t\lt v+B \\ \ 0; & \text{otherwise} \end{cases} }[/math]
  • Shifting property (shift along x-axis)
[math]\displaystyle{ \int_{t-B}^{t+B} x(\tau+\tau_0) e^{-j 2 \pi f \tau}\, d\tau = X(t+\tau_0,f)e^{j 2 \pi f \tau_0} }[/math]
  • Modulation property (shift along y-axis)
[math]\displaystyle{ \int_{t-B}^{t+B} [x(\tau) e^{j 2 \pi f_0 \tau}] d\tau = X(t,f-f_0) }[/math]
  • special input
  1. When [math]\displaystyle{ x(t)=\delta(t), X(t,f)=\begin{cases} \ 1; & |t|\lt B \\ \ 0; & \text{otherwise} \end{cases} }[/math]
  2. When [math]\displaystyle{ x(t)=1,X(t,f)=2B\operatorname{sinc}(2Bf)e^{j 2 \pi f t} }[/math]
  • Linearity property

If [math]\displaystyle{ h(t)=\alpha x(t)+\beta y(t) \, }[/math],[math]\displaystyle{ H(t,f), X(t,f), }[/math]and [math]\displaystyle{ Y(t,f) \, }[/math]are their rec-STFTs, then

[math]\displaystyle{ H(t,f)=\alpha X(t,f)+\beta Y(t,f) . }[/math]
  • Power integration property
[math]\displaystyle{ \int_{-\infty}^\infty |X(t, f)|^2\, df = \int_{t-B}^{t+B} |x(\tau)|^2\,d\tau }[/math]
[math]\displaystyle{ \int_{-\infty}^\infty \int_{-\infty}^\infty |X(t, f)|^2\,df\,dt = 2B \int_{-\infty}^\infty |x(\tau)|^2\,d\tau }[/math]
[math]\displaystyle{ \int_{-\infty}^\infty X(t,f)Y^*(t,f)\,df = \int_{t-B}^{t+B} x(\tau)y^*(\tau)\,d\tau }[/math]
[math]\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 }[/math]

Example of tradeoff with different B

Spectrograms produced from applying a rec-STFT on a function consisting of 3 consecutive cosine waves. (top spectrogram uses smaller B of 0.5, middle uses B of 1, and bottom uses larger B of 2.)

From the image, when B is smaller, the time resolution is better. Otherwise, when B is larger, the frequency resolution is better.

Advantage and disadvantage

Compared with the Fourier transform:

  • Advantage: The instantaneous frequency can be observed.
  • Disadvantage: Higher complexity of computation.

Compared with other types of time-frequency analysis:

  • Advantage: Least computation time for digital implementation.
  • Disadvantage: Quality is worse than other types of time-frequency analysis. The jump discontinuity of the edges of the rectangular mask results in Gibbs ringing artifacts in the frequency domain, which can be alleviated with smoother windows.

See also

References

  1. Jian-Jiun Ding (2014) Time-frequency analysis and wavelet transform