Mathematical Problems in Engineering
Volume 2017 (2017), Article ID 1805091, 14 pages
https://doi.org/10.1155/2017/1805091

Research Article

Applications of an Improved Time-Frequency Filtering Algorithm to Signal Reconstruction

Junbo Long,¹ Haibin Wang,² Daifeng Zha,³ Hongshe Fan,¹ Zefeng Lao,¹ and Huajie Wu¹

¹College of Electronic and Engineering, Jiujiang University, Jiujiang, China
²College of Information Science and Engineering Technology, Jiujiang University, Jiujiang, China
³College of Science, Jiujiang University, Jiujiang, China

Correspondence should be addressed to Junbo Long; 18488870@qq.com

Received 11 January 2017; Revised 28 March 2017; Accepted 6 April 2017; Published 8 May 2017

Academic Editor: Denis Benasciutti

Copyright © 2017 Junbo Long et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The short time Fourier transform time-frequency representation (STFT-TFR) method degenerates, and the corresponding short time Fourier transform time-frequency filtering (STFT-TFF) method fails under stable distribution noise environment. A fractional low order short time Fourier transform (FLOSTFT) which takes advantage of fractional order moment is proposed for stable distribution noise environment, and the corresponding FLOSTFT time-frequency representation (FLOSTFT-TFR) algorithm is presented in this paper. We study vector formulation of the FLOSTFT and inverse FLOSTFT (IFLOSTFT) methods and propose a FLOSTFT time-frequency filtering (FLOSTFT-TFF) method which takes advantage of time-frequency localized spectra of the signal in time-frequency domain. The simulation results show that, employing the FLOSTFT-TFR method and the FLOSTFT-TFF method with an adaptive weight function, time-frequency distribution of the signals can be better gotten and time-frequency localized region of the signal can be effectively extracted from stable distribution noise, and also the original signal can be restored employing the IFLOSTFT method. Their performances are better than the STFT-TFR and STFT-TFF methods, and MSEs are smaller in different and GSNR cases. Finally, we apply the FLOSTFT-TFR and FLOSTFT-TFF methods to extract fault features of the bearing outer race fault signal and restore the original fault signal from stable distribution noise; the experimental results illustrate their performances.

1. Introduction

The time-frequency analysis method is an effective tool for nonstationary signal processing, which has been broadly used in communication, acoustics, mechanical signal processing, medical signal processing, and others [1–5]. Short time Fourier transform time-frequency representation (STFT-TFR) has a rich structure with dependencies along both the time and frequency axes, which is widely applied in speech [6, 7], communication [8, 9], and others [10, 11]. A speech enhancement algorithm employing the STFT-TFR method is presented in [6], which takes advantage of the time and frequency dependencies of speech signals. An improved method which can reconstruct the spectral phase of the speech signal is proposed in [7]; the method only depends on the fundamental frequency and the noisy observation from the STFT domain. An interference excision algorithm is proposed for direct-sequence spread spectrum communication systems [8, 9], which analyze the interference in STFT time-frequency domain and detect and estimate the high-power jamming signal employing inverse STFT (ISTFT). Abdullah et al. apply the STFT-TFR method to analyze the switches faults of three-phase voltage source inverter [10]. Zhang et al. propose a time-frequency MUSIC algorithm employing the STFT-TFR method in [11], which is applied to estimate the DOAs of the closely spaced sources.

The time-frequency filtering technology which takes advantage of time-frequency localized spectra of the data provides an adaptive filtering method for the nonstationary signals. The time-frequency filtering method applies an adaptive weighting function to separate out the signals from the noise. The higher weighting parts localize the regions which are expected to be the signal components, and the lower weighting parts attenuate the noise in the time-frequency domain. The inverse transform of time-frequency representation is used to reconstruct the original signals. Recently, an adaptive time-frequency filtering method employing the Stockwell transform time-frequency representation (ST-TFR) method and its inverse transform are used to analyze the earthquake data [12, 13]; WANG et al. propose a new time-frequency filtering method employing normalized window S transform and TT transform in [14], which is applied to attenuate high frequency noise and random noise for the radar echo signal. An adaptive filtering method employing generalized S transform is presented in [15], which is constructed by a new adaptive time-frequency filtering factor, and the method is applied to restore the LFM signals from the noise. A generalized S transform is presented in [16], which is used to extract the special signal components and filter out the noise. CHI et al. propose an adaptive time-frequency filtering method based on the STFT-TFR method and apply it to filter out the noise in FM signal; the results demonstrate that the time-frequency filtering method is effective for FM signal filtering [17]. A FM interference suppression method employing the STFT-TFR and time-frequency filtering is presented in [18], which is applied to suppress the interference in the FM signals and improve the correlation output. Yang proposed an inverse STFT method and studied its vector formulation in [19].

The Gaussian noise is assumed in the above time-frequency representation methods and time-frequency filtering methods, and second-order statistics is applied to analyze the signals. In general cases, Gaussian hypothesis is reasonable, and the above methods are effective. However, in some actual signals such as the mechanical bearing fault signals, electroencephalogram, and radar signal, the noise has obvious pulsing characteristics, and the noise is non-Gaussian and nonstationary stable distribution [20–26]. When , the noise belongs to Gaussian distribution, and when , it is low order stable distribution. In low-order stable distribution environment, the mentioned methods in [12–19] are based on second-order statistics, and their performance will degenerate when Gaussian model is employed to analyze the non-Gaussian signals. Recently, stable distribution model is used for the bearing fault signals diagnosis in [27–30]. Fractional lower order statistic is applied to extract the evoked potential in [31, 32]. In stable distribution environment, the STFT-TFR method in the literature [8–11] degrades and even fails. Therefore, the corresponding time-frequency filtering method employing the STFT-TFF method in [17–19] degenerates. Hence, we propose fractional lower order short time Fourier transform and its inverse transform and propose a fractional lower order short time Fourier transform time-frequency representation (FLOSTFT-TFR) method for stable distribution noise environment. Based on the proposed FLOSTFT-TFR and IFLOST methods, we present a novel fractional lower order short time Fourier transform time-frequency filtering (FLOSTFT-TFF) algorithm which takes advantage of time-frequency localized spectra of the signal in time-frequency domain in this paper. The simulations show that the FLOSTFT-TFR method can better demonstrate time-frequency distribution of the signals, the FLOSTFT-TFF method can effectively restore the original signal under stable distribution noise environment, their performances are better than the STFT-TFR and STFT-TFF methods, and their MSEs are smaller under different and GSNR. We apply the FLOSTFT-TFR and FLOSTFT-TFF methods to extract fault features of the bearing outer race fault signal and restore the original fault signal from stable distribution noise, which illustrate their performances. Hence, the FLOSTFT-TFR method provides an approach for the special time-frequency analysis cases. IFLOSTFT is inverse transform of FLOSTFT, which provides a computational efficient way to restore the original signal from its time-frequency distribution when the undesired parts are removed.

In this paper, an improved time-frequency representation method employing fractional low order moment is proposed for stable distribution noise environment, and a novel time-frequency filtering method is presented. The paper is structured in the following manner. stable distribution is introduced in Section 2. In Section 3 we demonstrate a fractional low order short time Fourier transform and its inverse transform, the vector formulation of FLOSTFT and IFLOSTFT are derived, and the improved time-frequency representation is studied. The improved fractional lower order time-frequency filtering method is demonstrated in Section 4. The simulations are performed to demonstrate justifiability of the proposed time-frequency representation and time-frequency filtering methods based on FLOSTFT and IFLOSTFT, and the simulations of the outer race fault signals diagnosis are demonstrated in Section 5. Finally, the discussions and conclusions are given in Sections 6 and 7, respectively.

2. Stable Distribution

The characteristic function of stable distribution can be defined as [20, 21]

where is the characteristic index; it determines the pulse degree of stable distribution. If , is Gaussian distribution; if , is non-Gaussian fractional lower order distribution. Parameter is the dispersion coefficient, parameter is the location parameter, and parameter is the symmetry coefficient. If the parameter , is called symmetric stable distribution (), its characteristic function can be expressed as (2), and when , which is the same as the characteristic function of Gaussian distribution,

Probability density functions (PDFs) of SαS are given under , , , and in Figure 1.

Figure 1: PDFs of SαS stable distribution in different .

3. FLOSTFT and Its Inverse Transform

3.1. FLOSTFT

Fractional order moment of the signal is defined as

where is order moment operation of . When is a real signal, , and when is a complex signal, . is a real parameter, and [21, 22], demonstrate conjugate. Both sides of (3) are multiplied by a real window function centered at time , and localized spectrum of at is obtained, as shown in

After (4) has taken Fourier transform to time , we obtain

Equation (6) is called fractional low order short time Fourier transform (FLOSTFT); the corresponding discrete FLOSTFT can be written as

3.2. Computation of FLOSTFT

The discrete form of can be expressed as , ; then its discrete fractional order moment in (3) is written as

We divide into overlapping parts employing a real weighted function ; then the length of each part is , and the shift length from part to part is . The th part of can be written as

where . We define matrices   and as

where is identity matrix. is a diagonal matrix, which is composed of elements; the shift distance between and is . From (8) to (10), we can get the relationship between and as

We compute () points discrete Fourier transform to (9) and obtain

When all perform discrete Fourier transform, Fourier transform vector of is gotten:

We define Fourier transform matrix as

And let

is ) matrix. The weight function can be written as diagonal matrix, as shown in

By substituting (14)–(16) into (12), we can obtain

Similarly, (13) can be written as

where is identity matrix, denotes Kronecker tensor product, is matrix, and is matrix. , and is matrix. We call (18) fractional low order short time Fourier transform (FLOSTFT).

3.3. Computation of Inverse FLOSTFT

In this section, we introduce how to solve fractional order moments of from FLOSTFT. Both sides of (18) are multiplied matrix  ; we can get

Then

where is the estimation of . is matrix, and is matrix. We call (20) inverse fractional low order short time Fourier transform (IFLOSTFT).

3.4. FLOSTFT-TFR Method

The squared magnitude of in (7), denoted by , is

And the corresponding discrete form is written as

Equations (21)-(22) are called fractional low order short time Fourier transform time-frequency representation (FLOSTFT-TFR). By solving (22), we can get time-frequency representation of .

4. FLOSTFT Time-Frequency Filtering Algorithm

Each signal has certain regions and feature in time-frequency domain. The problem is illustrated in Figure 2, which shows that the FM signals gather in the time regions from to and the frequency regions from to . Then, we can separate the signal from the noise or eliminate undesired parts of the signal in time-frequency domain. In this section, we will study the time-frequency filtering algorithm based on the ISTFT and IFLOSTFT method under stable distribution noise.

Figure 2: Time-frequency representation of FM signal employing the FLOSTFT-TFR method under stable distribution noise.

4.1. STFT-TFF Algorithm

Short time Fourier transform (STFT) of is defined as

And the corresponding short time Fourier transform time-frequency representation (STFT-TFR) is expressed as

When both sides of (24) are multiplied by a weight function , we have

is designed according to time and frequency characteristics of the signal in time domain. From (25), we calculate inverse transform, yielding

where is the estimation of , which is filtered by the filtering function in time-frequency domain. When all is solved, we obtain . Following, we can calculate its inverse STFT employing the method in [19], and the estimated is gotten. The time-frequency filtering method is called short time Fourier transform time-frequency filtering (STFT-TFF) algorithm.

4.2. FLOSTFT-TFF Algorithm

The STFT-TFR method fails under stable distribution noise environment, and the STFT-TFF time-frequency filtering method based on the STFT-TFR method degenerates; hence, the FLOSTFT-TFR method has been demonstrated in Section 3. In this section, we will study how to restore the signal from its time-frequency representation employing the FLOSTFT-TFR method.

We multiply a weight function on both sides of (21) and then obtain

is a filter function, which is selected according to time-frequency characteristics of the signal. We define the estimation of as

where is th time point and th frequency point. When all points are calculated, we can obtain . Then, we can solve employing the IFLOSTFT method in Section 3.3. For the real signal , we have

From (29), we obtain

in (30) is the estimation of . We call the filtering method fractional low order short time Fourier transform time-frequency filtering (FLOSTFT-TFF) algorithm.

4.3. FLOSTFT-TFF Algorithm Steps

Step 1. Calculate FLOSTFT with (6).

Step 2. Compute FLOSTFT-TFR employing (21).

Step 3. Select appropriate filter function and extract the time-frequency regions of the signal employing (27).

Step 4. Calculate inverse FLOSTFT-TFR with (28), and compute IFLOSTFT employing (19)-(20).

Step 5. Perform inverse operation of with (30).

5. Simulations

The real signal contaminated by the noise may be written as

And two FM signals are expressed as

where , , , , , and is Gaussian noise or stable distribution noise. We let , . When is Gaussian noise, signal to noise ratio (SNR) can be used, however. If is stable distribution noise, SNR is unusable. Hence, we apply generalized signal noise ratio (GSNR) to substitute SNR, and . To measure the performance of the proposed methods, we define mixed mean square error (MSE) as

where is the number of Monte-Carlo experiments and and denote the signals estimated by the STFT-TFF or FLOSTFT-TFF methods.

5.1. Study of the FLOSTFT-TFR and STFT-TFR Methods

In this simulation, we apply the FLOSTFT-TFR and STFT-TFR methods to demonstrate the time-frequency distribution of FM signals under Gaussian noise and SαS stable distribution noise environment, respectively. If is Gaussian noise (), we set , . And if is SαS distribution noise, we set the characteristic index of SαS distribution noise , fractional low order moment parameter , , and length of the weighted function . The experimental results are shown in Figures 3(a)–3(d).

(a)

(b)

(c)

(d)

Figure 3: Time-frequency distributions of employing the STFT-TFR and FLOSTFT-TFR methods. (a) Time-frequency distribution employing the STFT-TFR method under Gaussian noise environment. (b) Time-frequency distribution employing the FLOSTFT-TFR method under Gaussian noise environment. (c) Time-frequency distribution employing the STFT-TFR method under SαS distribution noise environment. (d) Time-frequency distribution employing the FLOSTFT-TFR method under SαS distribution noise environment.

Figures 3(a) and 3(b) are the time-frequency distributions of employing the STFT-TFR and FLOSTFT-TFR methods under Gaussian noise environment, respectively. The figures show that the methods can also obtain time-frequency distribution of . Time-frequency distributions of employing the STFT-TFR and FLOSTFT-TFR methods under stable distribution noise environment are shown in Figures 3(c) and 3(d), respectively. The results show that the STFT-TFR method fails, and the FLOSTFT-TFR method has good performance. Hence, the FLOSTFT-TFR method has better performance in demonstrating the time-frequency distribution of the signals than the STFT-TFR method under stable distribution noise environment and is robust.

5.2. Study of the FLOSTFT-TFF and STFT-TFF Algorithm

In this simulation, we set as SαS distribution noise, , , , and . We apply the FLOSTFT-TFF and STFT-TFF methods to reconstruct the original FM signals from SαS distribution noise; the simulations are shown in Figure 4.

(a)

(b)

Figure 4: The filtered time-frequency representation employing the STFF-TFF and FLOSTFT-TFF method under SαS distribution noise environment. (a) Time-frequency representation of filtered by the STFF-TFF method. (b) Time-frequency representation of filtered by the FLOSTFF-TFF method.

Figures 4(a) and 4(b) are the time-frequency representations of filtered by the STFF-TFF and FLOSTFF-TFF methods under SαS distribution noise environment, respectively. The results show that the time-frequency representation filtered by the STFF-TFF method is incorrect, but the FLOSTFF-TFF method can better extract the time-frequency distribution of . Hence, the FLOSTTF-TFF method has better performance in extracting time-frequency distribution of the signal under SαS distribution noise environment.

In Figure 5 we show from (a) to (d) the original FM signal, the original FM signal contaminated by SαS distribution noise, the estimated FM signal employing the STFT-TFF method, and the FM signal estimated by the FLOSTFT-TFF method. The experimental results show that the contaminated original signal is overwhelmed by SαS distribution noise. The restored signal employing the STFT-TFF method has larger deviation, but the signal reconstructed by the FLOSTFT-TFF method is similar to the original signal. Hence, the FLOSTFT-TFF method has a better quality of signal reconstruction under stable distribution noise environment.

(a)

(b)

(c)

(d)

Figure 5: The real waveforms of two FM signals in time domain. (a) The original signal. (b) The original signal contaminated by SαS distribution noise. (c) The estimated signal employing the STFT-TFF method. (d) The estimated signal employing the FLOSTFT-TFF method.

5.3. MSE Analysis under Different and GSNR

In this simulation, is set as SαS distribution noise, , , , and . We apply the FLOSTFT-TFF and STFT-TFF methods to restore the original FM signals and compare MSEs of the methods under different GSNR and ; the simulations are shown in Figures 6 and 7.

Figure 6: The mixed MSEs employing the STFT-TFF and FLOSTFT-TFF methods under different (14 dB–26 dB).

Figure 7: The mixed MSEs employing the STFT-TFF and FLOSTFT-TFF methods under different (0–2).

Figure 6 is MSE analysis under different GSNR; the experimental results show that MSEs employing the STFT-TFF method vary from 8 dB to 80 dB when GSNR change from 26 dB to 14 dB, but MSEs employing the FLOSTFT-TFF method are in the range of −20 dB to −40 dB. Hence, the FLOSTFT-TFF method has obvious advantage to restore the original signals under different GSNR. In particular, the advantage of the FLOSTFT-TFF method is more obvious when .

The MSE analysis of the STFT-TFF and FLOSTFT-TFF methods under different is shown in Figure 7. We see that MSEs of the STFT-TFF method have a large variation, which vary from −20 dB to 220 dB, but MSE of the FLOSTFT-TFF method is stable about −30 dB. Hence, the FLOSTFT-TFF method has better stability in reconstructing the original signals.

5.4. Applications of FLOSTFT-TFF Method to Machine Fault Signals

In this simulation, the experiment signal is adopted from the Case Western Reserve University (CWRU) bearing data center [33]. The faulted bearings with fault diameters of 0.007 inches are reinstalled into the test motor, and the motor speed is 1797 RPM (revolutions per minute). The digital data is collected at 12,000 samples per second, and 0.2-second data is selected as test signal from the outer race bearing fault in BA, DE, and FE; then . SαS distribution noise (, ) is added in the fault signals, which is assumed as the actual working environment background noise. The STFT-TFR and FLOSTFT-TFR methods are used to demonstrate time-frequency distribution of the fault signals for fault analysis. Similarly, the STFT-TFF and FLOSTFT-TFF methods are also applied to restore the original bearing outer race fault signals from SαS distribution noise. The experimental simulations are shown in Figures 8–13.

(a)

(b)

(c)

(d)

Figure 8: Time-frequency representations of the machine fault signals in BA under SαS distribution noise environment. (a) Time-frequency representation of the machine fault signals in BA employing the STFT-TFR. (b) Time-frequency representation of the machine fault signals in BA employing the FLOSTFT-TFR. (c) The machine fault signals in BA restored by the STFT-TFF method. (d) The fault signals in BA reconstructed by the FLOSTFT-TFF method.

(a)

(b)

(c)

(d)

Figure 9: Time-frequency representations of the machine fault signals in DE under SαS distribution noise environment. (a) Time-frequency representation of the machine fault signals in DE employing the STFT-TFR method. (b) Time-frequency representation of the machine fault signals in DE employing the FLOSTFT-TFR method. (c) The machine fault signals in DE restored by the STFT-TFF method. (d) The fault signals in DE reconstructed by the FLOSTFT-TFF method.

(a)

(b)

(c)

(d)

Figure 10: Time-frequency representations of the machine fault signals in FE under SαS distribution noise environment. (a) Time-frequency representation of the machine fault signals in FE employing the STFT-TFR method. (b) Time-frequency representation of the machine fault signals in FE employing the FLOSTFT-TFR method. (c) The machine fault signals in FE restored by the STFT-TFF method. (d) The fault signals in FE reconstructed by the FLOSTFT-TFF method.

(a)

(b)

(c)

(d)

Figure 11: The waveforms of the machine fault signals in time domain. (a) The original bearing out race fault signals in BA. (b) The original bearing outer race fault signals in BA contaminated by SαS distribution noise. (c) The restored bearing outer race fault signals in BA employing the STFT-TFF method. (d) The bearing outer race fault signals in BA resconstructed by the FLOSTFT-TFF method.

(a)

(b)

(c)

(d)

Figure 12: The waveforms of the machine fault signals in time domain. (a) The original bearing outer race fault signals in DE. (b) The original bearing outer race fault signals in DE contaminated by SαS distribution noise. (c) The restored bearing outer race fault signals in DE employing the STFT-TFF method. (d) The bearing outer race fault signals in DE reconstructed by the FLOSTFT-TFF method.

(a)

(b)

(c)

(d)

Figure 13: The waveforms of the machine fault signals in time domain. (a) The original bearing outer race fault signals in FE. (b) The original bearing outer race fault signals in FE contaminated by SαS distribution noise. (c) The restored bearing outer race fault signals in FE employing the STFT-TFF method. (d) The bearing outer race fault signals in FE reconstructed by the FLOSTFT-TFF method.

Figures 8(a) and 8(b) are time-frequency distributions of the bearing outer race fault signal in BA employing the STFT-TFR and FLOSTFT-TFR methods, respectively. The STFT-TFR and FLOSTFT-TFR time-frequency distributions of the bearing outer race fault signals in DE and FE are shown in Figures 9(a), 10(a), 9(b), and 10(b), respectively. The figures show that the STFT-TFR method fails, but the FLOSTFT-TFR method can better demonstrate the time-frequency representations of the bearing outer race fault signals, which has good performance.

Figures 8(c)-10(c) and Figures 8(d)-10(d) are the time-frequency distributions filtered by the STFT-TFF and FLOSTFT-TFF methods, respectively. In the figures, we see that the extracted time-frequency representations employing the STFF-TFF method are incorrect, as shown in Figures 8(c)-10(c). But the FLOSTFT-TFF method in Figures 8(d)-10(d) can better fetch the time-frequency distribution of the fault signals. In Figures 8(d)-10(d) we may know that the gap regularly changes between the impacts, the fault vibration interval is approximately 30 ms, and the interval corresponds to the characteristic frequency of outer race which is 33.33 Hz. In Figure 8(d) the shock pulse is mainly distributed in low-frequency band from 500 Hz to 4500 Hz, and the transient harmonic vibration components of about 600 Hz, 2800 Hz, and 4200 Hz dominate frequency-domain. In Figure 9(d) the transient harmonic vibration components are about 600 Hz, 2800 Hz, and 4200 Hz, and the transient harmonic vibration components of about 600 Hz, 3400 Hz, 4500 Hz, and 5200 Hz are shown in Figure 10(d).

In Figures 11–13 we show from (a) to (d) the waveforms of the original bearing outer race fault signals in BA, DE, and FE, the original bearing outer race fault signals contaminated by SαS distribution noise, the bearing outer race fault signals reconstructed by the STFT-TFF method, and the reconstructed bearing outer race fault signals employing the FLOSTFT-TFF method, respectively. The simulations show that the contaminated bearing outer race fault signals in BA, DE, and FE are overwhelmed by SαS distribution noise in Figures 11(b)-13(b), the restored outer race fault signals employing the STFT-TFF method have larger deviation in Figures 11(c)-13(c), and the reconstructed outer race fault signals employing the FLOSTFT-TFF method are similar to the original signal in Figures 11(d)-13(d), respectively. Hence, the FLOSTFT-TFF method has a better quality of signal recovery for the bearing outer race fault signals.

Combined with the results of upper experiments, we see that the FLOSTFT-TFF method can better demonstrate time-frequency distribution of the machine fault signals and extract their fault features than the existing STFT-TFF method under SαS distribution noise environment.

6. Discussion

(i)Fractional order moment of the nonstationary signal is stationary and integrable within the time window ; however, the traditional STFT method is still no-stationary and integrable because when , .(ii)The shorter the length of weight function , the higher time resolution, but the lower the frequency resolution. In a practical application, we should set a reasonable balance value between time resolution and frequency resolution and select the proper length of window function.(iii)We can set the parameter according to characteristic index of stable distribution noise ; the smaller , the smaller . When , the FLOSTFT method degenerates into the STFT method, and the corresponding IFLOSTFT evolves into the ISTFT method,;hence, the improved FLOSTFT is a generalized algorithm.(iv)The FLOSTFT and IFLOSTFT methods are a linear transformation, but the FLOSTFT-TFR algorithm is a nonlinear transformation.

7. Conclusions

In this paper, we consider stable distribution noise as the background noise. Hence, the improved FLOSTFT and IFLOSTFT methods employing fractional p order moment are proposed, and FLOSTFT-TFR is presented for time-frequency analysis. We take advantage of time-frequency localized spectra of the signal and propose a novel FLOSTFT-TFF method, which can restore the original signal from stable distribution noise. We apply the FLOSTFT-TFR method to get time-frequency distribution of two FM signals and use the FLOSTFT-TFF method to extract the time-frequency distribution of two FM signals from SαS distribution noise; finally, the IFLOSTFT algorithm is used to restore the original FM signals; the results demonstrate their effectiveness. We compare the performance of the STFT-TFF and FLOSTFT-TFF methods; the results show that the performance of the FLOSTFT-TFF method is better than the STFT-TFF method, which can effectively suppress SαS distribution noise and work in low GSNR. In practical applications, we analyze the fault features from time-frequency representation of the mechanical bearing fault signals employing FLOSTFT-TFR method, extract time-frequency region of the fault signals from stable distribution noise employing FLOSTFT-TFF method, and reconstruct the original fault signal employing IFLOSTFT.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work is financially supported by Natural Science Foundation of China (61261046 and 61362038), the Natural Science Foundation of Jiangxi Province, China (20142BAB207006 and 20151BAB207013), the Research Foundation of Health Department of Jiangxi Province, China (20175561), and Science and Technology Project of Jiujiang University, China (2013KJ01, 2016KJ001, 2016KJ002, and 2016KJ003).

References

T. Higuchi, N. Ito, T. Yoshioka, and T. Nakatani, “Robust MVDR beamforming using time-frequency masks for online/offline ASR in noise,” in Proceeding of the 41st IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '16), pp. 5210–5214, China, March 2016.

X. Wang, Y. Zheng, Z. Zhao, and J. Wang, “Bearing fault diagnosis based on statistical locally linear embedding,” Sensors, vol. 15, no. 7, pp. 16225–16247, 2015.

A. Katunin, “Application of time-frequency distributions in diagnostic signal processing problems: a case study,” Diagnostyka, vol. 17, no. 2, pp. 95–103, 2016.

A. S. Alfahoum and A. Alfraihat, “Methods of EEG Signal Features Extraction Using Linear Analysis in Frequency and Time-Frequency Domains,” Isrn Neuroscience, vol. 2014, no. 3, Article ID 730218, 2014.

N. H. Chandra, A. S. Sekhar, N. H. Chandra et al., “Fault detection in rotor bearing systems using time frequency techniques,” Mechanical Systems & Signal Processing, vol. 72-73, pp. 105–133, 2015.

Y. Andrianakis and P. R. White, “A speech enhancement algorithm based on a chi mrf model of the speech STFT amplitudes,” IEEE Transactions on Audio Speech & Language Processing, vol. 17, no. 8, pp. 1508–1517, 2009.

M. Krawczyk and T. Gerkmann, “STFT phase reconstruction in voiced speech for an improved single-channel speech enhancement,” IEEE/ACM Transactions on Audio Speech & Language Processing, vol. 22, no. 12, pp. 1931–1940, 2014.

S. Aldirmaz and L. Durak, “Broadband interference excision in spread spectrum communication systems based on short-time fourier transformation,” Pier B, vol. 7, pp. 309–320, 2008.

S. Aldırmaz and L. Durak, “Performance analysis of an STFT-based broadband interference excision algorithm in DS-SS systems,” in Proceeding of Siu IEEE 15th Signal Processing & Communications Applications, pp. 1–4, June 2007.

A. R. Abdullah, N. S. Ahmad, N. H. Bahari et al., “Switched Faults Analysis of Voltage Source Inverter (VSI) using Short Time Fourier Transform (STFT),” International Review on Modelling & Simulations, vol. 7, no. 3, pp. 409–415, 2014.

H. Zhang, G. Bi, Y. Cai, S. G. Razul, and C. M. S. See, “DOA estimation of closely-spaced and spectrally-overlapped sources using a STFT-based MUSIC algorithm,” Digital Signal Processing: A Review Journal, vol. 52, pp. 25–34, 2016.

M. Schimmel and J. Gallart, “Authors’ Reply to Comments on 'The Inverse S-Transform in Filters With Time-Frequency Localization',” IEEE Transactions on Signal Processing, vol. 55, no. 10, pp. 5120–5121, 2007.

M. Schimmel and J. Gallart, “The inverse S-transform in filters with time-frequency localization,” IEEE Transactions on Signal Processing, vol. 53, no. 11, pp. 4417–4422, 2005.

D. Wang, J. Wang, Y. Liu et al., “An adaptive time-frequency filtering algorithm for multi-component LFM signals based on generalized S-transform,” in Proceeding of the International Conference on Automation & Computing, Glasgow UK, September 2015.

J. Wang, Y. G. He, G. F. Fang et al., “NWST and TT transform joint time-frequency filtering,” Application Research of Computers, vol. 29, no. 11, pp. 4254–4256, 2012.

X. H. Chen, H. E. Zhen, and D. J. Huang, “Generalized S Transform and its time-frequency filtering,” Signal Processing, vol. 24, no. 1, pp. 28–31, 2008.

H. S. Chi, H. X. Wang, Q. Guo et al., “Application of STFT in time frequency filtering of lfm signals,” Telecommunication Engineering, vol. 52, no. 2, pp. 155–159, 2012.

Z. Zha and X. Shi, “FM interference suppression based on STFT and time varying filtering,” Journal of Detection & Control, vol. 31, no. 3, pp. 46–50, 2009.

B. Yang, “A study of inverse short-time fourier transform,” in Proceeding of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '08), pp. 3541–3544, Las Vegas, NV, USA, April 2008.

X. Ma and C. L. Nikias, “Joint estimation of time delay and frequency delay in impulsive noise using fractional lower order statistics,” IEEE Transactions on Signal Processing, vol. 44, no. 11, pp. 2669–2687, 1996.

M. Shao and C. L. Nikias, “Signal processing with fractional lower order moments: stable processes and their applications,” Proceedings of the IEEE, vol. 81, no. 7, pp. 986–1010, 1993.

X. Ma and C. L. Nikias, “Parameter estimation and blind channel identification in impulsive signal environments,” IEEE Transactions on Signal Processing, vol. 43, no. 12, pp. 2884–2897, 1995.

T.-H. Liu and J. M. Mendel, “A subspace-based direction finding algorithm using fractional lower order statistics,” IEEE Transactions on Signal Processing, vol. 49, no. 8, pp. 1605–1613, 2001.

L. R. Padovese, “Hybrid time-frequency methods for non-stationary mechanical signal analysis,” Mechanical Systems and Signal Processing, vol. 18, no. 5, pp. 1047–1064, 2004.

A. G. Poulimenos and S. D. Fassois, “Parametric time-domain methods for non-stationary random vibration modelling and analysis—a critical survey and comparison,” Mechanical Systems and Signal Processing, vol. 20, no. 4, pp. 763–816, 2006.

J. B. Long and H. B. Wang, “Parameter estimation and time-frequency distribution of fractional lower order time-frequency auto-regressive moving average model algorithm based on SαS process,” Journal of Electronics & Information Technology, vol. 38, no. 7, pp. 1710–1716, 2016.

C. Li and G. Yu, “A new statistical model for rolling element bearing fault signals based on alpha-stable distribution,” in Proceeding of the International Conference on Computer Modeling and Simulation (ICCMS '10), pp. 386–390, China, January 2010.

X. M. Yu and T. Shu, “Fault diagnosis method for gearbox based on α-stable distribution parameters and support vector machines,” Measurement and Control Technology, vol. 31, no. 8, pp. 21–30, 2012.

G. Yu, C. N. Li, and J. F. Zhang, “A new statistical modeling and detection method for rolling element bearing faults based on alpha-stable distribution,” Mechanical Systems and Signal Processing, vol. 41, no. 1-2, pp. 155–175, 2013.

J. B. Long, H. B. Wang, P. Li, and H. S. Fan, “Applications of fractional lower order time-frequency representation to machine bearing fault diagnosis,” IEEE/CAA Journal of Automatica Sinica, vol. 99, pp. 1–17, 2016.

J. B. Long, H. B. Wang, and D. F. Zha, “Evoked potential blind extraction based on fractional lower order spatial time-frequency matrix,” Journal of Biomedical Engineering, vol. 32, no. 2, pp. 269–274, 2015.

Z. J. Lin, D. F. Zha, and J. Sheng, “Blind estimation of evoked potentials based on covariations in non-gaussian noise,” Journal of Biomedical Engineering, vol. 27, no. 4, pp. 727–730, 2010.

CWRU bearing data center: http://csegroups.case.edu/bearingdatacenter/pages/download-data-file.