Complexity
Volume 2018 (2018), Article ID 9039240, 10 pages
https://doi.org/10.1155/2018/9039240

Research Article

Effective Approach to Calculate Analysis Window in Infinite Discrete Gabor Transform

Rui Li,¹ Yong Huang,¹ and Jia-Bao Liu²

¹College of Information and Network Engineering, Anhui Science and Technology University, Bengbu 233030, China
²School of Mathematics and Physics, Anhui Jianzhu University, Hefei 230601, China

Correspondence should be addressed to Jia-Bao Liu; liujiabaoad@163.com

Received 3 December 2017; Revised 28 March 2018; Accepted 10 April 2018; Published 15 May 2018

Academic Editor: Danilo Comminiello

Copyright © 2018 Rui Li 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 long-periodic/infinite discrete Gabor transform (DGT) is more effective than the periodic/finite one in many applications. In this paper, a fast and effective approach is presented to efficiently compute the Gabor analysis window for arbitrary given synthesis window in DGT of long-periodic/infinite sequences, in which the new orthogonality constraint between analysis window and synthesis window in DGT for long-periodic/infinite sequences is derived and proved to be equivalent to the completeness condition of the long-periodic/infinite DGT. By using the property of delta function, the original orthogonality can be expressed as a certain number of linear equation sets in both the critical sampling case and the oversampling case, which can be fast and efficiently calculated by fast discrete Fourier transform (FFT). The computational complexity of the proposed approach is analyzed and compared with that of the existing canonical algorithms. The numerical results indicate that the proposed approach is efficient and fast for computing Gabor analysis window in both the critical sampling case and the oversampling case in comparison to existing algorithms.

1. Introduction

The discrete Gabor transform (DGT) [1], extended from the short-time Fourier transform (STFT), is an important time-frequency analysis tool for processing and analyzing the nonstationary signal [2–4], which multiplied a signal by a time-shifted and frequency-modulated window function with the aim of representing analyzed signals in a time-frequency localized manner. As a general case of STFT, DGT, sampled version of the STFT, is an efficient tool for analyzing the energy distribution of analyzed signals in the time-frequency domain/plane and more valuable than Fourier transform (FT) in many applications which require to extract temporal frequency information and localize simultaneously. According to the type of time-frequency atom, the DGT can be classified into two categories: the complex-valued discrete Gabor transform (CDGT) and the real-valued Discrete Gabor transform (RDGT) [5–8]. The former method for computing the DGT and its inverse transform all involves complex operators while the latter one can utilize discrete Hartley’s cas transform (DHT) kernel, discrete cosine transform (DCT) kernel, or discrete sine transform (DST) kernel as the Gabor basis function which only involves real operators. More specifically, in some real-time applications, the CDGT is more difficult to be implemented in hardware or software and the RDGT for sampled speech and image has an advantage of the computationally efficient implementation.

Over the past decades, the DGT has become a very valuable and widely used mathematic tool in diverse applications [12–19] such as audio processing, speech processing, image processing, sonar and radar signal processing, seismic signal processing, and transient signal processing. The analyzed signal, analysis window, and synthesis window in the periodic DGT usually have the same periodic length, which may confront difficulties in applying to some real-time applications of long-periodic signal sequences. To bridge these gaps, the long-periodic DGT [9, 10, 19] of complex-valued kernel and real-valued kernel has been presented to utilize the short window to analyze and process the long-periodic (or infinite) signal sequences in practical applications. Because existing canonical algorithms [11, 20–27] are mainly used in periodic/finite DGT and derived from Gabor frame theory, in this paper, a fast and effective algorithm, based on the orthogonal analysis approach and FFT algorithm, is proposed to obtain the analysis window for the long-periodic/infinite DGT in both the critical sampling case and the oversampling case.

The rest of the paper is organized as follows. In Section 2, the long-periodic/infinite DGT will be reviewed. In Section 3, a fast approach, which uses the new orthogonal constraint relationship of the analysis window and the synthesis window to deduce a series of independent linear equation sets, is introduced to calculate analysis window in the long-periodic/infinite DGT under both the critical sampling case and the oversampling case. In Section 4, the detailed comparison of computational complexity of related algorithms and numerical experiments have been given to show that the proposed method can overcome the existing algorithms in computing analysis window of the long-periodic/infinite DGT. Finally, the paper is concluded with Section 5.

2. DGT for Long-Period and Infinite Sequences

Let represent a periodic and discrete-time signal with length and the length of analysis window and synthesis window . In order to use and to analyze and process the signal , the three periodic sequences , , and of length should be constructed as follows:

where The DGT for long-periodic/infinite sequences [4, 19] is defined by the following:

and the Gabor coefficients can be calculated by

where is the imaginary unit. The DGT for long-periodic/infinite sequences can be obtained by (5) corresponding to its expansion in (4). In (4) and (5), , let the positive integers , are the sampling points in time and frequency, and let and denote a modulation step in frequency and a translation interval in time. For a numerically stable reconstruction, the constrained condition by (or ) has to be satisfied. By choosing proper parameters and , the Gabor oversampling rate is a positive integer. The biorthogonal relationship [4, 19] of and , which are derived from the completeness condition between (4) and (5), can be rewritten as

where denotes the discrete form of delta function, , and .

3. FFT-Based Method for Solving the Analysis Sequences

To derive a new type biorthogonality relationship from the completeness between (4) and (5), substituting (5) into (4) yields

recalling

and putting (8) into (7) leads to

According to the completeness of the transform, the following constraint relationship has to be satisfied:

so

Equation (11) can be expressed as a equivalent form in the following:

where and . Let , , and ; then (12) can be written as

Obviously, (13) can be rewritten in matrix form as follows:

where is an -long unit vector with first element being one and others being zeros; is a -long vector composed by

and is a matrix constructed by

where , , and . Equation (14) can be split into unrelated linear equation sets:

Due to the fact that most of elements in and were zeros according to (2) and (3), let , be the positive integer constrained by : (15) and (16) are rewritten as follows:

where

, , and . Equation (17) can then be rewritten in the following form:

where is a unit vector of length . Obviously, the following formula can be easily obtained by (20):

where .

3.1. Fast Approach for Computing Analysis Window in the Critical Sampling Case

In the critical sampling case, the oversampling rate ; (21) can be rewritten in following form:

Because is a right circulant matrix by (or ), so (22) can be written as a circular convolution in following form:

where is the first column of as the following form:

By utilizing the circular convolution theorem [28], the FFT algorithm can then be used to transform (25) into component-wise multiplication by recalling (26) as follows:

where is the -point DFT. Thus, can be fast computed by

The computation of (29) is obviously much faster than that using the standard Gaussian elimination method to solve the linear equation set.

Remark 1. More similar result to (29), the author in [11] provides a method, based on Gabor frame theory [20] and discrete Zak transform [29], to obtain dual Gabor window in periodic/finite DGT under the critical sampling. However, in present paper, the proposed method uses the orthogonal analysis approach and FFT algorithm to solve the analysis window for long-periodic/infinite signal sequences.

3.2. Fast Approach for Computing Analysis Window in the Oversampling Case

In the oversampling case, the oversampling rate can be a positive integer by setting and properly. Because the rank of is less than the number of rows, (21) is an underdetermined linear equation which could have no solution or have an infinite number of solutions. Equation (21) could be rewritten in the following form in order to use convolution theory and FFT algorithm as that in the critical sampling case.

where and is a left circulant matrix composed by

where and is a vector with length constructed by

Then we rewrite (30) as

where

Let

Then, (33) also can be rewritten as a circular convolution:

where is the first column of as follows:

By using the circular convolution theorem and the FFT operator, (36) can be expressed in the following form:

so can be solved by

Due to the fact that should be close to in the sense of the least square error (LSE) as shown in [30],

has to satisfy the least of norm which make a Gaussian-type function of the localization property:

Equation (41) can be converted into the following optimum problem and easily solved by the Lagrangian approach:

where , is a vector with length , is a Lagrangian parameter, and

Taking a derivative of with respect to leads to

then the solution of can be obtained from (44):

and putting this solution into constraint condition in (42) to solve leads to

where .

Remark 2. In sparse applications, the synthesis window without considering localization property () can be used to reconstruct the original signal, in which a certain well-localized window function will be directly utilized as an analysis window.

Remark 3. The proposed algorithm can be easily generalized and applied to other DGT for both the periodic/finite sequences and the long-periodic/infinite sequences including complex-valued (DFT-based) and real-valued (DHT-based, DCT-based) and multiwindow discrete Gabor transforms (M-DGT) because those Gabor basis functions in DGT can lead to the similar form of the biorthogonal relationship equation.

4. Computational Complexity Analysis and Numerical Experiments

When calculating the multiplication times of related algorithms, we assume that the -point FFT requires an order of multiplications. Due to the fact that the original linear equation of biorthogonal relationship can be split into unrelated subequation sets in the critical sampling case and subequation sets in the oversampling case, FFT method can be utilized to compute analysis window both in the critical sampling case and in the oversampling case. The computational complexity of the single subequation set which carries -point FFT and -point inverse FFT (IFFT). The computational complexity of is at order of in the oversampling case while it is zero in the critical sampling case. Therefore, as mentioned above, the computational complexity of proposed approach is equal to that of total subequation set, which is at the order of

in the critical sampling case and

in the oversampling case. The comparison of computational complexity and numerical experiments between proposed approach and existing methods have been given in Tables 1 and 2 which clearly shows that the total number of multiplications of proposed approach is less than that of others in both the critical sampling case and the oversampling case. In Tables 1 and 2, the symbols , CS, and OS denote , the critical sampling case, and the oversampling case, respectively. In Table 2, a detailed numeric comparison experiments on the number of multiplications related to computational time between the proposed algorithm and existing methods under different Gabor sampling scheme are given by using the formula of computational complexity of each algorithm in Table 1, which clearly demonstrates the efficiency and the advantage of decreasing the computation of the proposed approach as compared to existing algorithms. Example for computing analysis window: a Gaussian synthesis window with length is given (49) and shown in Figure 1.

The analysis windows in Figure 2 are computed by the proposed method, [9, 10] in the critical sampling case, with , and the oversampling rate . The analysis windows in Figures 3 and 4 are computed by the proposed method [9, 10] in the oversampling case, with the oversampling rate (, ) and 8 (, ), respectively. The similarity between the synthesis window and its corresponding analysis window is proportional to the oversampling rate , the reason being similar to that given in [30]. The proposed algorithm for the reconstruction of the original signal is also simulated, which is virtually error-free reconstruction (MSE = ).

Table 1: Comparison of computational complexity between proposed approach and other existing methods.

References
Computational complexity
Applicability

[9]

CS, OS

[10]

CS, OS

[11]
,
CS, OS

Proposed algorithm
,
CS, OS

Table 2: Numerical comparison of multiplications between proposed approach and other existing methods.

Total number of multiplications
Applicability

[9]
[10]
[11]
Proposed

algorithm

32
32
1
794624
258209280
2560
2560
CS

64
16
4
151552
9810880
7168
2304
OS

256
8
32
35328
122415
34304
1792
OS

64
64
1
6324224
2106322432
6144
6144
CS

128
16
8
286720
17715136
24576
4608
OS

512
8
64
70144
241455
134144
3584
OS

128
128
1
50462720
17014709760
14336
14336
CS

256
32
8
2260992
146317248
59392
11264
OS

1024
16
64
544768
2029663
335872
9216
OS

Figure 1: A Gaussian synthesis window , .

(a)

(b)

(c)

Figure 2: Analysis windows in the critical case (, , ) obtained by the proposed method [9, 10].

(a)

(b)

(c)

Figure 3: Analysis windows in the oversampling case (, , ) obtained by the proposed method [9, 10].

(a)

(b)

(c)

Figure 4: Analysis windows in the oversampling case (, , ) obtained by the proposed method [9, 10].

5. Conclusion

The DGT for long-periodic/infinite signal sequences has become an important time-frequency analysis tool in many real-time applications, which can use the short window to process and analyze the long-periodic signal sequences. This paper proposed an effective approach to calculate the analysis window of DGT for long-periodic sequences, in which the original biorthogonal equation can be decomposed into a certain number of linear subequation sets that convolution operators and FFT could be utilized to obtain the analysis window for arbitrary given synthesis window function. The computational complexity analysis and comparison between proposed approach and existing methods have been shown in a comparative study described in Section 4, which clearly indicates that the proposed approach is more competitive against the existing algorithms. In addition, the numerical experiments have been given to demonstrate the efficiency and validity of proposed approach, which make the long-periodic DGT attractive for real-time signal processing and analysis.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The work of is partly supported by the National Science Foundation of China under Grant no. 11601006, China Postdoctoral Science Foundation under Grant no. 2017M621579, the Postdoctoral Science Foundation of Jiangsu Province under Grant no. 1701081B, the Natural Science Foundation for the Higher Education Institutions of Anhui Province of China under Grant no. KJ2015A331, and the Doctoral Scientific Research Foundation of Anhui Science and Technology University under no. XWYJ201801.

References

D. Gabor, “Theory of communication,” Journal of the Institute of Electrical Engineers of Japan, vol. 93, pp. 429–457, 1946.

M. J. Bastiaans, “Gabor's Expansion of a Signal into Gaussian Elementary Signals,” Proceedings of the IEEE, vol. 68, no. 4, pp. 538-539, 1980.

J. Wexler and S. Raz, “Discrete Gabor expansions,” Signal Processing, vol. 21, no. 3, pp. 207–220, 1990.

S. Qian and D. Chen, “Discrete Gabor Transform,” IEEE Transactions on Signal Processing, vol. 41, no. 7, pp. 2429–2438, 1993.

D. F. Stewart, L. C. Potter, and S. C. Ahalt, “Computationally Attractive Real Gabor Transforms,” IEEE Transactions on Signal Processing, vol. 43, no. 1, pp. 77–84, 1995.

C. Lu, S. Joshi, and J. M. Morris, “Parallel lattice structure of block time-recursive generalized Gabor transforms,” Signal Processing, vol. 57, no. 2, pp. 195–203, 1997.

L. Tao and H. K. Kwan, “Noise reduction for NMR FID signals via oversampled real-valued discrete gabor transform,” IEICE Transaction on Information and Systems, vol. E88-D, no. 7, pp. 1511–1518, 2005.

L. Tao and H. K. Kwan, “Block Time-recursive Discrete Gabor Transform Implemented by Unified Parallel Lattice Structures,” in Proceedings of the 2005 IEEE International Symposium on Circuits and Systems, pp. 844–847, Kobe, Japan.

S. Qian, Introduction to Time-Frequency and Wavelet Transforms, China Machine Press, 2005.

L. Tao, G. H. Hu, and H. K. Kwan, “Multiwindow Real-Valued Discrete Gabor Transform and Its Fast Algorithms,” IEEE Transactions on Signal Processing, vol. 63, no. 20, pp. 5513–5524, 2015.

T. Strohmer, Numerical Algorithms for Discrete Gabor Expansions, 1998.

H. Rabbani and S. Gazor, “Local probability distribution of natural signals in sparse domains,” in Proceedings of the 36th IEEE International Conference on Acoustics, Speech, and Signal Processing, ICASSP 2011, pp. 1289–1292, May 2011.

M. Temerinac-Ott and M. Temerinac, “Discrete fourier-invariant signals: Design and application,” IEEE Transactions on Signal Processing, vol. 60, no. 3, pp. 1108–1120, 2012.

L. Bianco, D. Gottas, and J. M. Wilczak, “Implementation of a gabor transform data quality-control algorithm for UHF wind profiling radars,” Journal of Atmospheric and Oceanic Technology, vol. 30, no. 12, pp. 2697–2703, 2013.

M. Matthe, L. L. Mendes, and G. Fettweis, “Generalized frequency division multiplexing in a gabor transform setting,” IEEE Communications Letters, vol. 18, no. 8, pp. 1379–1382, 2014.

N. Saini and A. Sinha, “Face and palmprint multimodal biometric systems using Gabor–Wigner transform as feature extraction,” Pattern Analysis and Applications, vol. 18, no. 4, pp. 921–932, 2015.

B. Friedlander and B. Porat, “Detection Of Transient Signals By The Gabor Representation,” IEEE Transactions on Signal Processing, vol. 37, no. 2, pp. 169–180, 1989.

Y. Lu, S. Joshi, and J. M. Morris, “Noise reduction for NMR FID signals via Gabor expansion,” IEEE Transactions on Biomedical Engineering, vol. 44, no. 6, pp. 512–528, 1997.

L. Tao and H. K. Kwan, “Real discrete Gabor expansion for finite and infinite sequences,” Proceedings - IEEE International Symposium on Circuits and Systems, vol. 4, pp. 637–640, 2000.

M. Zibulski and Y. Y. Zeevi, “Oversampling in the Gabor Scheme,” IEEE Transactions on Signal Processing, vol. 41, no. 8, pp. 2679–2687, 1993.

H. Bölcskei, “A Necessary and Sufficient Condition for Dual Weyl-Heisenberg Frames to be Compactly Supported,” Journal of Fourier Analysis and Applications, vol. 5, no. 5, pp. 409–419, 1999.

O. Christensen, “Pairs of dual Gabor frame generators with compact support and desired frequency localization,” Applied and Computational Harmonic Analysis , vol. 20, no. 3, pp. 403–410, 2006.

P. L. Søndergaard, “Gabor frames by sampling and periodization,” Advances in Computational Mathematics, vol. 27, no. 4, pp. 355–373, 2007.

R. S. Laugesen, “Gabor dual spline windows,” Applied and Computational Harmonic Analysis , vol. 27, no. 2, pp. 180–194, 2009.

P. L. Søndergaard, “Efficient algorithms for the discrete Gabor transform with a long FIR window,” Journal of Fourier Analysis and Applications, vol. 18, no. 3, pp. 456–470, 2012.

S. Li, Y. Liu, and T. Mi, “Sparse Dual Frames and Dual Gabor Functions of Minimal Time and Frequency Supports,” Journal of Fourier Analysis and Applications, vol. 19, no. 1, pp. 48–76, 2013.

J. Zhou, X. Fang, and L. Tao, “A Sparse Analysis Window for Discrete Gabor Transform,” Circuits, Systems and Signal Processing, vol. 36, no. 10, pp. 4161–4180, 2017.

R. M. Gray, “Toeplitz and circulant matrices: a review,” Foundations and Trends in Communication and Information Theory, vol. 2, no. 3, pp. 155–239, 2006.

J. Zak, “Finite translations in solid-state physics,” Physical Review Letters, vol. 19, no. 24, pp. 1385–1387, 1967.

S. Qian, K. Chen, and S. Li, “Optimal biorthogonal functions for finite discrete-time Gabor expansion,” Signal Processing, vol. 27, no. 2, pp. 177–185, 1992.