S transform

Definition

There are several ways to represent the idea of the S transform. In here, S transform is derived as the phase correction of the continuous wavelet transform with window being the Gaussian function.

• S-Transform : S_x(t,f) = \int_{-\infty}^\infty x(\tau)|f|e^{- \pi (t- \tau)^2 f^2} e^{-j2 \pi f \tau} , d \tau
• Inverse S-Transform :x(\tau) = \int_{-\infty}^\infty \left[\int_{-\infty}^{\infty}S_x(t,f), dt\right],e^{j2\pi f\tau}, df

Modified form

• Spectrum Form The above definition implies that the s-transform function can be expressed as the convolution of ( x(\tau) e^{-j2 \pi f \tau} ) and ( |f|e^{- \pi t^2 f^2} ).

Applying the Fourier transform to both ( x(\tau) e^{-j2 \pi f \tau} ) and ( |f|e^{- \pi t^2 f^2} ) gives : S_x(t,f) = \int_{-\infty}^\infty X(f+\alpha),e^{-\pi\alpha^2 /f^2},e^{j2\pi\alpha t}, d\alpha .

• Discrete-time S-transform From the spectrum form of S-transform, we can derive the discrete-time S-transform.

Let t = n\Delta_T,, f = m\Delta_F,, \alpha = p\Delta_F, where \Delta_T is the sampling interval and \Delta_F is the sampling frequency.

The Discrete time S-transform can then be expressed as:

:S_x(n\Delta_T, ,m\Delta_F) = \sum_{p=0}^{N-1} X[(p+m),\Delta_F],e^{-\pi\frac{p^2}{m^2}},e^{\frac{j2pn}{N}}

Implementation of discrete-time S-transform

Below is the Pseudo code of the implementation.

Step1.Compute X[p\Delta_{F}],

loop over m (voices) Step2.Compute e^{-\pi \frac{p^2}{m^2}}for f=m\Delta_{F}

Step3.Move X[p\Delta_{F}] to X[(p+m)\Delta_{F}]

Step4.Multiply Step2 and Step3 B[m,p] = X[(p+m)\Delta_{F}]\cdot e^{-\pi \frac{p^2}{m^2}}

Step5.IDFT(B[m,p]). Repeat.}

Comparison with other time–frequency analysis tools

Comparison with Gabor transform

The only difference between the Gabor transform (GT) and the S transform is the window size. For GT, the windows size is a Gaussian function ( e^{-\pi (t-\tau)^2} ), meanwhile, the window function for S-Transform is a function of f. With a window function proportional to frequency, S Transform performs well in frequency domain analysis when the input frequency is low. When the input frequency is high, S-Transform has a better clarity in the time domain. As table below.

This kind of property makes S-Transform a powerful tool to analyze sound because human is sensitive to low frequency part in a sound signal.

Comparison with Wigner transform

The main problem with the Wigner Transform is the cross term, which stems from the auto-correlation function in the Wigner Transform function. This cross term may cause noise and distortions in signal analyses. S-transform analyses avoid this issue.

Comparison with the short-time Fourier transform

We can compare the S transform and short-time Fourier transform (STFT). First, a high frequency signal, a low frequency signal, and a high frequency burst signal are used in the experiment to compare the performance. The S transform characteristic of frequency dependent resolution allows the detection of the high frequency burst. On the other hand, as the STFT consists of a constant window width, it leads to the result having poorer definition. In the second experiment, two more high frequency bursts are added to crossed chirps. In the result, all four frequencies were detected by the S transform. On the other hand, the two high frequencies bursts are not detected by STFT. The high frequencies bursts cross term caused STFT to have a single frequency at lower frequency.

Applications

• Signal filterings
• Magnetic resonance imaging (MRI)
• Power system disturbance recognition
  • S transform has been proven to be able to identify a few types of disturbances, like voltage sag, voltage swell, momentary interruption, and oscillatory transients.
  • S transform also be applied for other types of disturbances such as notches, harmonics with sag and swells etc.
  • S transform generates contours which are suitable for simple visual inspection. However, wavelet transform requires specific tools like standard multiresolution analysis.
• Geophysical signal analysis
  • Reflection seismology
  • Global seismology

References

References

1. (1996). "Localization of the complex spectrum: the S transform". IEEE Transactions on Signal Processing.
2. Stockwell, RG (1999). ''S''-transform analysis of gravity wave activity from a small scale network of airglow imagers. PhD thesis, University of Western Ontario, London, Ontario, Canada.
3. (2008). "2008 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society".
4. (January 2010). "A General Description of Linear Time-Frequency Transforms and Formulation of a Fast, Invertible Transform That Samples the Continuous S-Transform Spectrum Nonredundantly". IEEE Transactions on Signal Processing.
5. Kelly Sansom, "Fast S Transform", University of Calgary, https://www.ucalgary.ca/news/utoday/may31-2011/computing
6. (13 August 2018). "Fast S-Transform download {{!}} SourceForge.net".
7. Brown, Robert. (2024-12-01). "robb-brown/GFT".
8. E. Sejdić, I. Djurović, J. Jiang, "Time-frequency feature representation using energy concentration: An overview of recent advances," ''Digital Signal Processing'', vol. 19, no. 1, pp. 153-183, January 2009.
9. (2012). "Analysis of non-stationary structural systems by using a band-variable filter". Bulletin of Earthquake Engineering.
10. Hongmei Zhu, and J. Ross Mitchell, "The S Transform in Medical Imaging," University of Calgary Seaman Family MR Research Centre Foothills Medical Centre, Canada.
11. (2010). "2010 IEEE International Conference on Power and Energy".