Abstract

The short-time linear canonical transform (STLCT) is a novel time-frequency analysis tool, which has attracted some attention recently. However, its applications in signal processing are limited because the time-frequency properties of the STLCT are still little known. Most existing studies focus on mathematical properties rather than time-frequency properties in signal processing. To handle this problem, first, we investigate some basic time-frequency properties such as 2-D resolution of the time-frequency plane, the STLCT domain support, computation of the STLCT, etc by generalizing the short-time Fourier transform to the STLCT. Then, based on these derived properties, we find that the Gaussian window is optimal window of the STLCT. Signal separation verified the results.

1. Introduction

The short-time Fourier transform (STFT) has been the widely used as time-frequency analysis tool for many years. It can extract the characteristic behavior of a signal in joint time-frequency domains. Recently, the linear canonical transform (LCT), a generalization of the traditional Fourier transform (FT) and well-known fractional Fourier transform (FRFT), has received much attention in signal processing [18]. The LCT is more flexible because it has a total of four free parameters and one constraint. Many theoretical issues of the FT and FRFT have been extended to the LCT, which has found many research results in signal processing [922]. But the LCT fails in locating frequency-domain contents due to its global kernel function. For some practical applications, we want to know not only the LCT frequency contents but also how they change by time. For this reason, Kou and Xu [23] first systematically introduced the short-time linear canonical transform (STLCT), in which two mathematical theorems were developed.

The motivation of this paper is to study the time-frequency properties of the STLCT, especially the optimal window selection. Since the STLCT is a generalized form of classical transformation, such as the STFT, its possible applications theoretically include nonstationary signal filtering, radar signal analysis, and so on. The premise of these applications is the property of the STLCT, such as sampling theorem, calculation algorithm, and so on [18]. To our knowledge, few achievements on time-frequency properties associated with the STLCT have been reported. Thus, we account for time-frequency properties of the STLCT from the point of view of signal processing. Serving as a novel time-frequency analysis tool, window selection and applications of the STLCT are still little known because of the following two reasons.

The first reason is most related work focused more on its mathematical theories than time-frequency properties. For instance, in [23], Kou and Xu proposed the STLCT and developed many useful theorems including covariance property, orthogonality property, and limit theorem in function space. Then, in [24], they also proved the Paley–Wiener theorems and the uncertainty principles for the (inverse) windowed linear canonical transform. Subsequently, the STLCT has also received attention by scholars. In [25], Zhang addressed a sampling theorem for the STLCT by means of generalized Zak transforms associated with the LCT and investigated biorthogonal condition of series expansion of the STLCT. Besides, Guanlei et al. discussed the uncertainty inequation relation in two STLCT domains [26]. Recently, in [27], we studied the uncertainty principle in one STLCT domain based on the definition of time width and frequency bandwidth of its time-frequency plane. Most studies paid more attention to its mathematic properties. We study time-frequency performance of the STLCT in this paper. In fact, the LCT has been introduced to classic time-frequency analysis tool, such as Wigner distribution (WD) and STFT [2830]. For instance, in [29], Zhang and Luo considered the time-frequency tool and its application using the Gaussian windowed STFT and LCT. The performance of frequency location depends on sliding window of the STFT.

The second issue is that there are few studies on the time-frequency properties of the STLCT. As a novel time-frequency analysis tool, these basic properties are prerequisites of applications. For example, the principle of selecting the optimal window function was studied in [31]. So, we want to know the upper bound of 2-D resolution of the time-frequency plane and the minimum STLCT domain support. That is to say, how to select an appropriate window will have an impact on performance of the STLCT. For practical applications, digital computation of the STLCT is also needed to study. Besides, some other time-frequency features such as inverse STLCT, time delay, and frequency shift properties are necessary supplementary principles to its applications.

The contribution of this paper is twofold. First, we develop the time-frequency properties of the STLCT by generalizing the short-time Fourier transform (STFT) to the STLCT by substituting the FT kernel function with the LCT kernel function, including linearity, time delay and frequency shift, inverse STLCT, etc. Second, since performance of the STLCT depends largely on the window function, we prove that the Gaussian window function is optimal. The Gaussian-windowed STLCT has both maximum 2D resolution and minimum time-frequency domain support. Furthermore, performance of the STLCT affected by different windows is studied by signal separation.

Our paper is organized as follows. In Section 2, the LCT and the STLCT as well as the previous works about time-frequency analysis are reviewed. In Section 3, framework of the STLCT is established, including the idea of the STLCT, inverse STLCT, 2-D resolution of the time-frequency plane, digital computation algorithm, and STLCT domain support. In Section 4, an optimal STLCT is developed and some potential applications are discussed. In Section 5, we summarize the paper and make a future direction.

2. Preliminaries

2.1. Linear Canonical Transform

The linear canonical transform (LCT) of a signal with parameter A is defined as [3]where a, b, c, and d are real numbers satisfying . The kernel function of the LCT is defined as

The LCT satisfies the additivity property as follows:where

Thus, inverse LCT can be deduced through the additivity property, particularly,where is the inverse of matrix. Other properties of the LCT, such as sampling, eigenfunction, uncertainty principles, fast computation, and relationships with other transforms, can be found in detail in [28].

The LCT is a generalized linear transform. When , it reduces to the FRFT. When , it becomes the classic FT. When the parameter , the LCT is a scaling transform operation multiplying a linear frequency modulation signal, and we only consider in our paper.

2.2. Short-Time Linear Canonical Transform

The LCT has the same limitations as the FT. It means that the LCT fails in locating the linear canonical domain frequency contents due to its global kernel function. So, in [23], Kit proposed the short-time linear canonical transform (STLCT) to handle this problem. Actually, a STLCT was mentioned in uncertainty principles study even earlier [17]. The definition of the STLCT iswhere is the selected window function with a short time support and is the LCT kernel function. By definition, the result of the integral transform contains information about time and frequency . So, we can know not only the STLCT frequency-domain contents but also how they change with time. The basis of the STLCT is defined as

3. Theoretical Framework of the STLCT

In this section, we will establish the time-frequency analysis theories of the STLCT. Firstly, basic time-frequency properties of the STLCT including inverse STLCT, time delay, and frequency shift properties are deduced. Then, in order to achieve a better implementation of the STLCT, an optimal window selection for STLCT is developed. According to the classic STFT, there are several kinds of window function such as rectangular window, Hamming window, Gaussian window, and so on. But only the Gaussian window can obtain the lower bound of the uncertainty principle of STFT and achieve minimum STFT domain support, which implies that Gaussian windowed FT is optimal STFT [31]. Secondly, as mentioned above, the STFT is a special case of the STLCT. So, we want to find out if these properties of the STFT can be generalized to the STLCT, and we seek the results via both mathematical proof and simulation verification.

3.1. Time-Frequency Properties of the STLCT
3.1.1. Property of Linearity

This property was mentioned in [23], and the STLCT is linear:where and are arbitrary real constants. This property is convenient to process multicomponent signals and will be used in application part in this paper.

3.1.2. Time Delay and Frequency Shift Properties

STLCT is windowed LCT and can be viewed as a result of finitely repeated LCT. Thus, the time delay and frequency shift properties of the STLCT remain unchanged with the LCT, respectively. They are shown as follows:

These properties can be used to deduce the case that both time delay and frequency shift occur to signal. In this case, the result is shown below.where

3.1.3. Inverse STLCT

1-D Inverse STLCT. We first assume that the reconstruction formula is

Substituting (6) into (12), we obtain

Then, simplifying (13) yields

In order to obtain the perfect reconstruction, let ; therefore, the constant coefficient should satisfy the following requirement:

So, the 1-D inverse STLCT is

2-D Inverse STLCT. We first assume that the reconstruction formula is

Substituting (6) into (17) yields

Then, by simplifying (18), we have

A simple form can be obtained through variable substitution:

Similarly, to ensure , the following equation should be satisfied:

Therefore, we usually choose an appropriate window function to meet the condition . The 2-D inverse STLCT is

3.1.4. Uncertainty Principle

In time-frequency analysis field, Heisenberg’s principle is a fundamental result in signal processing. Different time-frequency analysis tools such as FT and FRFT have different lower bounds. Let us s take the FT for example; a signal cannot be both narrowband and short duration, as the variances of FT pairs cannot be arbitrarily small. So, the lower bound is largely dependent upon the selected technology. The STFT has the same lower bound with the FT (Section 3, [26]). However, the STLCT includes the STFT, and we want to know not only whether the STLCT keeps the uncertainty property of its special case unchanged but also the reality verified by the simulation. Our previous study in [20] showed that lower bound of the STLCT remains the same with the LCT. Also, in this section, in order to get actual performance of the STLCT which could guide potential engineering practice, we choose the signal , where , , and , as our test signal. We study the Gaussian-windowed STLCT with different window lengths (see Figure 1).

3.1.5. 2-D Resolution of the STLCT

The STLCT can map the time-domain signal into the joint time and frequency domain. From this 2-D transformation, we can know not only the LCT domain frequency contents but also how they modify with time via its time-frequency plane. This plane is divided by many parallelograms while the plane of the GT is divided by numerous rectangles. The two sides of parallelogram represent the time width and the LCT domain bandwidth, respectively, based on uncertainty principle analysis in [9]. In a similar manner, they are defined bywhere E denotes the energy of signal, and specially, for an energy normalized signal, we have . is the LCT of . Thus, the definition of is given as follows:

We can rewrite the above equation by substituting the kernel function of the LCT. The computation process can be found in Appendix A.where is the LCT of with .

Similar to the classic uncertainty principle in FT domain, and are the mean values of time domain and frequency domain, respectively. They are defined as

The 2-D resolution of the STLCT is defined as

From the definition of the resolution, the performance is affected by selected window function and free parameters of the STLCT. In [20], we have proved that the uncertainty principle of the STLCT is

The equality holds if and only if Gaussian window function is selected. That is to say, the 2-D resolution can be obtained as

Clearly, it has an upper bound when we chose Gaussian window. We summarize this section as a theorem as follows.

Theorem 1. The upper bound of 2-D resolution of the STLCT is , which is achieved if and only if the Gaussian window is used.

3.1.6. STLCT Domain Support

It is known that classic time-frequency plane of the STFT is divided by many cell grids (see Figure 2), which is also called time-frequency cells. The two sides of the cell represent the STFT time width and band width. Moreover, each cell grid is bounded by rectangles. The coordinates of the center of rectangle are mean time and frequency, and two sides of the rectangle are the time width and bandwidth of the STFT. The product of the two sides (time-frequency product) is often defined as a STFT domain support.

The STLCT is generalized STFT. Therefore, similar to the STFT, the time width, bandwidth, mean time, and mean bandwidth in STLCT domain are defined as

To seek the relationship between the selected window function and the STLCT domain, we have the following equation, and the computation process can be found in Appendix B.where and are the squared time and bandwidth of the given window function. Also, the STLCT domain support is defined as squared product of and .

3.2. Computation of the STLCT

The basic idea of computation of the STLCT is as follows. First, original signal is intercepted and captured by a sliding window in time domain. Second, the instantaneous STLCT spectrum in time domain is obtained by LCT of the signal portion captured by the window. Then, by moving the window along time domain, we can see how the frequency contents change by time. As the STLCT is based on LCT, to utilize the STLCT to the practical application, we choose a simple and efficient method of the LCT proposed in [29]. The discrete LCT and its inverse are as follows. Also, the computation steps of the STLCT are summarized as follows.Step 1: select an appropriate window function and initialize a set of parameters of the LCT.Step 2: choose the length of window.Step 3: add the window to the given signal and compute the LCT of the signal cut by the window.Step 4: move the window along the time till the end of the signal.

From the idea of the STLCT, we can see that computation complexity of the STLCT is the same as the LCT because the STLCT is the sum of finite LCT result of each segment signal (see Figure 3).

4. Results

4.1. Optimal STLCT

In this section, we demonstrate that the Gaussian window is the optimal window for the STLCT. The performance of the STLCT is determined by selected window and free parameters according to the above study. In reality, we usually preset the parameters and then select a window based on the situation. The time and frequency resolution of the STLCT is constrained by the uncertainty principle, which means that a high-frequency resolution leads to a coarse time resolution and vice versa. Thus, there is a trade-off between time and frequency resolution in choosing an appropriate window for the STLCT. Window with short length is more appropriate for signals of short duration, whereas a long length window with high-frequency resolution is more appropriate for narrowband signals.

Durak and Arikan put forward a criterion of optimal window selection in [31]. Optimal windowed STLCT has both maximum 2D resolution and minimum time-frequency domain support. According to the time-frequency properties (5) and (6) in Section 3, we demonstrate that the Gaussian window also has the minimum STLCT domain support, which is concluded as the following theorem.

Theorem 2. The minimum STLCT domain support exists if and only if the selected window is Gaussian, and the expression of the window function is

Proof. According to Theorem 1, the Gaussian-windowed STLCT domain has maximum 2D resolution, so we suppose an energy normalized window functionThen, squared time and bandwidth of the window can be computed by integral formula . We obtainIn addition, by substituting (30) and (33) into (31), we haveBased on basic inequality property,The equality is achieved if and only if . Furthermore, the result can be reduced to the STFT and STFRFT cases when and , respectively.
It is known that there are some commonly used windows such as rectangular window, Hamming window, and Gaussian window. We choose the signal used in Section 3 to test the performance of the STLCT with different windows under the same parameters (a = 10, b = 0.5, d = 10, and window length = 128). It demonstrates that the Gaussian window leads to fine resolution and concentration of the time-frequency plane, relative to other windows (see Figure 4).

4.2. Signal Separation Application

Signal separation is a vital setup in many real applications such as fault diagnosis, feature selection, and classification. In this section, the STLCT is applied to signal separation. Especially, the Gaussian window is verified the optimal window of the STLCT through the application. Here we use the nonstationary mixed signals (see Figure 5).

The simulation results show that the Gaussian window function has the best time-frequency plane resolution and aggregation. In this section, the selective signal separation is applied to verify the superiority of Gaussian window. Signal separation is an important task in signal processing, which is often used in practical processing tasks such as filter design and fault diagnosis. For example, in fault diagnosis domain, it is often necessary to extract features from signals. For signals with fault and abnormal features, people want to separate them from normal signals for diagnosis. In this section, we will take the mixed signal of the following two nonstationary signals as an example (as shown in Figure 5) and separate them in the time-frequency plane by using the time-frequency analysis technology based on short-time linear canonical transformation.

The main steps of simulation are as follows.Step 1: initialize the parameters of the short-term linear canonical transformation matrix and start the type and width of window function to process the original signal.Step 2: the time-frequency distribution of the two signals is obtained by short-time linear canonical transformation.Step 3: mask the signals to be removed to separate them.Step 4: do the inverse transform of short-time linear canonical transform for the reserved signal.

The experimental results are shown in Figure 6. The results show that in terms of time-frequency distribution, the resolution and aggregation degree of time-frequency distribution of traditional short-time Fourier transform with window length of 256 are worse than those of short-time linear canonical transformation with window length of 512 (as shown in Figure 1). The longer the window length is, the less likely it is to reflect the frequency location characteristics theoretically. Compared with rectangular window and Hamming window, Gaussian window has the best time-frequency resolution and aggregation for short-time linear canonical transform. Especially in the third step, mask filter has better filtering effect because of its high aggregation. Therefore, qualitative analysis based on simulation verifies the conclusion that Gaussian window is the optimal window function. When the parameters are given, the short-term linear canonical transform will degenerate into the traditional optimal short-time Fourier transform, namely, Gabor transform.

4.3. Discussion and Some Potential Applications
4.3.1. Advantages of the STLCT

From the properties deduced in Section 3, the advantage of the STLCT is that it can obtain a higher resolution compared to the STFT and also has no cross terms in the time-frequency plane. If the Gaussian window is taken, the optimal STLCT will become the generalized Gabor transform (GT) with three free parameters. It is known that the GT has found many applications in signal analysis, image processing, and pattern recognition. As its generalized form, the STLCT has more degree of freedom.

4.3.2. Weakness of the STLCT

Although the STLCT is a new time-frequency analysis tool, it has not been used in practical scenarios so far. The classic short-time Fourier transform (STFT) (special case of the STLCT) has been widely used in the field of nonstationary signal processing. In particular, when the window function is Gaussian function, the short-time Fourier transform is Gabor transform (GT). However, the GT did not attract scholars’ interest in the decades after it was proposed. The main reason is that the expansion basis of the GT composed of window function and Fourier kernel function is nonorthogonal, which makes it difficult to calculate Gabor expansion coefficients. The STLCT is the generalized form of the STFT, so the STLCT also suffers from the difficulty of calculating its coefficients.

One of the potential applications of the STLCT is nonstationary signal filtering. Image or signal often has noise and interference. If different regions of a signal or image are contaminated with noise of different distribution characteristics, we can separate the target signal from noise in time-frequency plane based on high concentration by selecting suitable parameter matrix. We can also choose different free parameters for different regions to achieve a better performance rather than working with a fixed parameters for the whole signal or image. For example, we consider the LFM signal, the most important nonstationary signals in radar and sonar systems. Based on optimization algorithm, maximization local energy concentration can be obtained by choosing an optimal parameter matrix [14, 15]. Another potential application is in radar signal analysis. Moving target detection technology in complex background is one of the key technologies and difficulties of radar target detection, and the STFRFT has excellent performance in processing multicomponent signals, especially in high-resolution sparse domain signal processing. Since classic FRFT and STFRFT are special forms of the STLCT. Therefore, in theory, the STLCT should also have their advantages in radar target detection and parameter estimation.

4.3.3. Calculation Acceleration of the STLCT

Based on the STLCT calculation algorithm in this paper, the calculation complexity of the STLCT mainly depends on the calculation of LCT. However, the large-scale signals generated in practical engineering restrict the application of the STLCT. Fortunately, in [32, 33], Zhang et al. and Liu et al. proposed a sparse discrete FRFT algorithm to reduce the computational complexity when dealing with large signal data. Inspired by Liu et al.’s work, a feasible future direction is to study sparse LCT because LCT is the generalized form of FRFT.

5. Conclusions

In this paper, we have studied the time-frequency properties of the STLCT including linearity, time delay and frequency shift, inverse STLCT, computation of the STLCT, 2-D resolution, and time-frequency plane support. Gaussian window is proved to be the optimal window for the STLCT. The performance of uncertainty principle and optimal STLCT is qualitatively verified via simulation. Some potential applications also have been discussed. STLCT-based filtering design and IFE are now under investigation.

Appendix

A. Computation of (26)

B. Proof of (31)

To simplify (31), we first compute as follows:

On the other hand, suppose the window function is real and symmetric, and we obtain

So, we have the following equation:

By variable substitution,

The Jacobian partial derivative is

Then, (38) can be rewritten as

Thus, and are obtained:

Data Availability

The selected signal and parameter setting data are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (62006093) and Natural Science Research Project of the Higher Education Institutions of Jiangsu Province, China (18KJB510006).