Efficient Time-Frequency Localization of a Signal

A new representation of the Fourier transform in terms of time and scale localization is discussed that uses a newly coinedA-wavelet transform (Grigoryan 2005). The A-wavelet transform uses cosineand sine-wavelet type functions, which employ, respectively, cosine and sine signals of length 2π. For a given frequencyω, the cosineand sine-wavelet type functions are evaluated at time points separated by 2π/ω on the time-axis. This is a two-parameter representation of a signal in terms of time and scale (frequency), and can find out frequency contents present in the signal at any time point using less computation. In this paper, we extend this work to provide further signal information in a better way and name it as A-wavelet transform. In our proposed work, we use cosine and sine signals defined over the time intervals, each of length 2πm/(2ω), m ≤ 2, m and n are nonnegative integers, to develop cosineand sine-type wavelets. Using smaller time intervals provides sharper frequency localization in the time-frequency plane as the frequency is inversely proportional to the time. It further reduces the computation for evaluating the Fourier transform at a given frequency. The A-wavelet transform can be derived as a special case of the A-wavelet transform.


Introduction
There are two domains for representing a signal: time and frequency domains.Depending upon the information required, either representation can be used.Fourier analysis has been the main technique for transforming a signal from one representation into another.In spite of the fact that the Fourier analysis is an ideal solution for deterministic and stationary signals, it is hardly of any use for time-varying signals or nonstationary signals, because analysis of these types of signals compromises between their transition and long term behaviors.For these types of signals, a transform is desired to represent the signal in a two-parameter form.The very first such transform in literature is the short time Fourier transform (STFT) [1].This transform uses a time-window function to decompose the signal into segments and then the Fourier analysis is carried out on individual segments.The STFT provides local features that are present in the signal in a limited form because it uses the time-window function of fixed width for all frequency contents and thus it is unable to extract the required information in any given signal.So, a transform that can represent the signal in two-parameter form and uses the time-window function of different lengths is needed.Application of the short time Fourier transform multiple times by using different time-window functions can be a possible solution, but it is a cumbersome process and hence may not be practically feasible.So, a transform that can represent the signal in two-parameter form and uses a time-window function of varying length is required.Such functions/transforms exist in literature since long time.Prior to a decade and so, their applications had hardly been explored in signal processing.These transforms are called wavelet transforms [2][3][4][5][6][7][8].In the last two decades and so, large numbers of articles have appeared on the wavelet transforms in literature [9][10][11][12][13].A wavelet transform employs a set of variable length time-window functions derived from a single time-window function, called the mother wavelet.The advantage of this type of representation of a signal is that the signal information can be obtained in different frequency bands.
There have been some studies to represent the Fourier transform in such a way that one does not need to know the function on the entire time-axis in order to find its Fourier transform in a given frequency range, possibly at a given frequency.One such study has been reported in [14].In that study, the cosine and sine signals truncated over a period of International Journal of Digital Multimedia Broadcasting 2 have been considered to develop the cosine-and sinetype wavelets, which have been used to represent the Fourier transform in two-parameter form unlike the traditional form.This representation is important because one does not need to compute the values of cosine-and sine-type wavelets at all time points.Their values are computed at certain time points, which are multiple of 2/.In other words, this provides the signal information along the curves that are separated out by 2/ in the time-frequency plane, thus requiring less computation.Besides, it provides multiresolution signal processing.This work has motivated to develop a new representation of the Fourier transform that can provide the signal information in time-frequency plane along the curves separated out by less than 2/ distance.The proposed representation, named as  * -wavelet transform, further reduces the computations and provides sharper frequency localization.In this paper, the cosine-and sine-type wavelets are derived from the cosine and sine signals.The length of the time intervals over which the cosine and sine signals are defined can be adapted as per the application requirement.Thus, the Fourier transform represented in terms of cosine-and sine-type wavelet transform can provide the required signal information by adjusting the time interval length of the cosine and sine signals.For example, the wavelets derived from the (co)sinusoidal signals with smaller period give sharper frequency localization and vice versa.
The rest of the paper is organized as follows.Section 2 reviews A-wavelet transform.In Section 3, the  * -wavelet transform is proposed.Section 4 provides results and discussions and Section 5 elucidates the theory by using examples.Lastly, the paper is concluded in Section 6.

A-Wavelets
The Fourier analysis describes the global nature of the signals.
It is beneficial if the signal does not have time-varying nature; that is, it is stationary.There are many time-varying signals, which have practical applications such as speech signals, audio signals, and video signals.For such applications, the Fourier analysis is less beneficial as these signals contain local characteristics different from global ones, which are required for further analysis.The short-time Fourier transform (STFT) provides solutions for such types of applications by using a real and symmetric time-window function having unit norm in  2 () and nonzero only in the region of interest.Here  2 () refers to the space of square-integrable functions.The STFT of a signal (), denoted by (, ), using the window function () and reconstruction of the original function () from its STFT are given by The time-window function () is of fixed size, which limits the STFT to determine the locations of occurrence of many frequency contents.This problem has been alleviated by applications of the wavelet transforms.A wavelet transform employs a varying size time-window function, called wavelet function, by dilating and translating in its time parameter.Mathematically, the wavelet transform, denoted by (, ), of an arbitrary function () using the wavelet function () is given by where  and  are called dilation and translation parameters, respectively.
In [14], A-wavelet transform has been discussed that represents the Fourier transform in terms of cosine-and sinetype wavelet transforms.The cosine-and sine-type wavelet transforms are based on the truncated cosinusoidal and sinusoidal signals, denoted by () and (), respectively, and they are defined over a time interval of length 2; that is, The functions   (,   ) and   (,   ) are called cosinewavelet (C-wavelet) and sine-wavelet (S-wavelet), respectively.The A-wavelet transform representation of the Fourier transform () is given by where  is an imaginary quantity. ( This representation (i.e., C-wavelet and S-wavelet) reduces the computations because, for a given frequency , the functions   (,   ) and   (,   ) are evaluated at the centers of the prespecified time intervals of length 2/ in time parameter.The frequency and time parameters are related by the relation  = 2.The number of time intervals on time-axis is given by 2() + 1 (i.e., −(), −() + 1, . . ., −1, 0, 1, . . ., (), () + 1), where () = 0 for  ∈ [0, 1/3], and for  > 1/3, () = ⌊(3 − 1)/2⌋.For details on () refer to [14].The A-wavelet transform provides the signal information in the time-frequency plane along the curves  = 2,  = 0, ±1, ±2, . ... In the next section, we discuss  * -wavelet transform, an extension of A-wavelet transform, that provides the signal information along the curves separated by less 2 in the time-frequency plane.

𝐴 * -Wavelet Transform
The Fourier transform of a function in terms of A-wavelet transform is expressed in terms of C-and S-type wavelet transforms and these transforms have been derived from the cosine and sine functions over a time interval of length 2 (refer to (3a) and (3b)).The A-wavelet transform provides the signal information along the curves  = 2, where  is an integer, that is, all such curves that are separated by 2 in the time-frequency plane.We call the curves  = 2,  = 0, ±1, ±2, . .., as primary curves and the curves between the primary curves as the secondary curves.The signal information along the primary curves can be obtained using the A-wavelet transform, but the A-wavelet transform cannot provide the signal information along the secondary curves.We want to get the signal information along the secondary curves in the current paper.For doing this, we develop C-and S-type wavelets using the cosine and sine signals defined over a time interval of length less than 2.
The interval length is parameterized by introducing two new integer variables  and .Define the functions  , () and  , () as follows: We look at different values of  and .For  = 0,  , () is the Dirac Delta function and  , () is identically zero and the outcome is simply a sampled signal.For  = 2  , the functions  , () and  , () are exactly the same as the () and () functions, respectively, which are given in (3a) and (3b).For  > 2  , the functions  , () and  , () have nonzero values over the time duration, that is, more than one time interval.We are interested to have these signals as nonzero over a single time interval, which gives a condition  ≤ 2  .Without loss of generality, we may assume that  and  assume nonnegative integer values.For negative , the same discussion will hold as that for positive , but on the negative time-axis.Negative  will increase the number of time intervals over which the functions  , () and  , () have nonzero values.The condition  < 2  makes the time interval length smaller over which the functions  , () and  , () are nonzero and this is the requirement to develop the current work.By doing this, we intend to get the signal information along the secondary curves.We construct a family of functions { ,,, (),  ,,, ()} from  , () and  , () by using the time-scaling and shift transformation; that is,  →  and  → ( −   ); that is, where the frequency variable  assumes real values and   = 2/(2  ),  = 0, ±1, ±2, . . ., ±⌊(3−1)2  /2⌋.In general,  is infinite; however, for finite support signals, its value is finite.
Consider a function () that satisfies all the conditions for existence of Fourier transform.Its Fourier transform () is given by For  = 0, the Fourier transform is a real quantity that is given by integrating the function.This is a trivial case and needs no further discussion.Our further discussion is meant for  ̸ = 0 frequencies.We write () that is defined in (8) as follows: We try to write the integrals in (9) in terms of  ,,, () and  ,,, () functions that are defined in (7).For this purpose, we uniformly divide the time-axis into disjoint intervals of length 2/(2  ) that are centered at 2/(2  ),  = 0, ±1, ±2, . ... Denote th interval by  ,, , which can be defined as follows: In (10),  assumes all integer values in general; however, its maximum value is determined by the frequencies, denoting them by Θ , , that need be analyzed.For example,  = 1 and  = 3 gives the frequency range Θ 1,3 = [/(3 * 2  ), 3/(3 * 2  )).The secondary curves are given by the equation  = 2/2  , ( ̸ = 2  ), where  is a positive integer and  is nonnegative integer.We can write the Fourier transform () defined in (9) over the disjoint intervals  ,, as follows: Changing the limits of integration, we have International Journal of Digital Multimedia Broadcasting On simplifying, it gives which can further be written as follows: In ( 14) the summation contains the term  −2/2  .We can write 2  / =  + , where  is a nonnegative integer and  is an integer such that 0 ≤  < .Thus, we have  −2/2  =  −2/(+) ,  = 0, 1, 2, . . .,  − 1,  = 0, ±1, ±2, . ... Since  assumes all integer values including zero, the expression  −2/(+) will have ( + ) distinct values for all values of .For positive and negative values of , we will have same results.For positive , we will have intervals along positive axis and for negative , we have intervals on negative axis; henceforth we will confine  as nonnegative integer.For fixed , increasing  will increase the value of , but  will have only  distinct values.Increasing  will help in analyzing the high frequencies.Increasing  will increase the values of  ( ̸ = 0) and that in turn will increase the secondary curves.It may be noticed that the curves for  = 0 correspond to the primary curves.Increasing  and  will increase both primary and secondary curves in the time-frequency plane.The primary curves depend upon the time interval in which we want to analyze the desired frequency contents and the secondary curves provide further details.In limiting case, the distance between two adjacent curves can be made arbitrarily small by increasing .If the signal is periodic, then the primary curves are sufficient enough to provide the desired information present in the signal; otherwise the secondary curves are also needed.The information corresponding to the primary curves is the same as that provided by the Fourier series for the periodic signals.The number of secondary curves in a given time interval is decided by the frequencies to be analyzed in the given signal.

Results and Discussions
The Fourier transform in (16) has been decomposed into different components and each component that has been represented in terms of time and scale localization can be analyzed explicitly and independently.The work [14] contains only one component; whereas the current work contains multiple components.Thus, it provides much more information about the desired signal.In fact, the signal information obtained in [14] can be given by the first component  0 () in a better localized form using the current work.The signals   () ( ̸ = 0) give the information about the signal along the curves  , = 2/( + 1),  = 1, 2, 3, . . ., ( +  − 1) in the time-frequency plane.For  = 0, it gives information along the primary curves.The equation  , = 2/ indeed represents a family of hyperbolic curves of second kind, one for each value of , in the time-frequency plane.Maximum value of  or the number of hyperbolic curves is determined by the values of  and .The set of functions { ,,, (),  ,,, ()} constitutes a basis to represent the Fourier transform of a given signal in terms of scale and time localization.The Fourier transform is given by the functions  ,,, (,  (+)+ ) and  ,,, (,  (+)+ ) (these functions include the signal whose Fourier transform is determined).For a given frequency, these functions are evaluated at various points, each lying on a different curve.In this work, we have discretized the frequency-time plane in different curves.The value of  used in (6a) and (6b) determines the number of curves in the time-frequency plane.The maximum distance between two consecutive curves can be 2, the coarsest discretization that has been considered in [14].The minimum distance (in limiting case) can be zero, which corresponds to the continuous domain and in that case the number of curves is infinite.In conventional Fourier transform, we implicitly consider infinite number of curves and, in the Fourier series, finite number of curves, each separated by 2.That is why infinite frequencies (or time) are needed to represent the Fourier transform at a given point in time (or frequency).
The results derived in this paper help finding the values of  ,,, (,  (+)+ ) and  ,,, (,  (+)+ ) over the time interval [−/2  , /2  ) = [−/( + ), /( + )).This can be analogous to a stack of rectangular blocks, each having time and scale (frequency) as its two sides.The zeroth block ( = 0) has the time parameter length equal to the signal support and the frequency length is decided by the frequencies to be analyzed.If the desired frequencies are not resolved in the zeroth block, then the first block ( = 1) is considered whose time parameter length is half of that of the zeroth block.If the first block fails to resolve the desired frequencies, then the second block ( = 2) is considered whose time parameter length is half of that of the first block or one-fourth of the signal support, and so on.It may be noted that the segment in time parameter of a block is appropriately chosen from the support of the signal maintaining the length of th block as /2  , where  is the signal support.The zeroth block is used to analyze the smallest frequencies starting from zero.The parameter  specifies distinct (type) curves in a block and  specifies a particular (type) curve in that block.A block may contain several curves of the same type.Theoretically, a signal contains infinite frequencies, which means there are infinite numbers of blocks.For analyzing high frequencies, we need higher indexed blocks.For taking  sufficiently large, we can decrease the distance between two curves as small as we please.Practical signals are generally both band-limited and time-limited.Thus, most of the time, we have both the number of blocks and number of time samples (curves) as finite.The work [14] considers only primary curves for representing the signal information and the work [16] considers both primary and secondary curves for representing the signal information, but only one block that corresponds to the first block of our proposed work that consists of multiple blocks, each containing both primary as well as secondary curves.Figures 1(a)-1(d) show frequency contents in different blocks for a signal with support for  ∈ (4,36).We want to determine the frequency contents in the vicinity of the time point 4.5.Figure 1(d) shows that the frequency contents up to three decimal points can be   resolved.If we want to resolve further closer frequencies, we need to increase the number of blocks.
(23) For  = 0 and  = 1, (22) gives It is not difficult to show that the signal information given by  0 () is the same as that obtained in [14].For  = 4 and  = 4, the plot of the equation  , = 2/2  is shown in Figure 2. (Two consecutive curves are /4 apart).Figure 2 helps determining the number of points in set R defined in (23) in the time-frequency plane for a given frequency range.For example, for all frequencies less than 1/12, only one point is needed; for frequencies less than 1/4, two points are needed, and so on.In this figure, the curves drawn correspond to the loci of points in Y e s the time-frequency plane for  ∈ (0, 5),   ∈ (−5, 5), where   signifies the shift in time and  = 4.In this figure, there are four types of curves, each corresponding to   (),  = 0, 1, 2, 3.The dark black curves are the primary curves that correspond to first component  0 () of ( 22).These are the curves that can also be identified using the A-wavelet transform [14].The important point here is that though the curves are the same, yet the time intervals considered in the current paper are of smaller size and as a result this gives better frequency localization.The dashed, dash-dotted, and dotted curves correspond to second component  1 (), third component  2 (), and fourth component  3 (), respectively.These are the secondary curves which cannot be obtained using the A-wavelet transform.

Illustration
We will elucidate the above discussed theory with an example and  ,,, (,  (+)+ ) are to be evaluated, for  = 1.75, is 21, that is, 0, ±1, ±2, ±3, . . ., ±10.Since  = 4, the Fourier transform contains four components, that is,  0 (),  1 (),  2 (), and  3 (), and since  = 1, we need to consider the curves corresponding to the first block. 0 () needs be computed at the center of every fourth time interval in both directions ( = 0, ±4, ±8) starting from the first one (assuming the index of the first interval as 0), that is, (,  0 ), (, ± 4 ), and (, ± 8 ), where   represents the time (location) parameter;  1 () is computed at the center of every fourth time interval in both directions starting from the second one, that is, (, ± 1 ), (, ± 5 ), and (, ± 9 );  2 () is computed at the center of every fourth time interval in both directions starting from the third one, that is, (, ± 2 ), (, ± 6 ), and (, ± 10 );  3 () is computed at the center of every fifth time interval in both directions starting from the fourth one, that is, (, ± 3 ) and (, ± 7 ).Apart from these computations, no values of any other function or translation of any function is needed.The above discussion has been summarized in Table 1.
We now illustrate the computation of the Fourier transform and compare the results with those of the Awavelet transform [14] and SC-and SS-wavelet transforms [16].We take the signal () = cos(1.3)+ cos(2.5)+ cos(3.6)+ cos(4.8)and evaluate  ,,, (,  (+)+ ) and  ,,, (,  (+)+ ),  = 0,  = 1,  = 4, and  = 0, 1, . . ., ( +  − 1) using the current work.Substituting the values of () and  ,,, () in (20a), we have where The corresponding plots in 3D and 2D are shown in Figures 3(a) and 3(d), respectively, for the first block.In fact, these figures can also be obtained by using the work [16].However, the plots shown in Figure 3(b) can only be obtained by using the current work.We have omitted the 2D graphs for the current work because Figure 3(b) is good enough to show the advantage over the works [14,16].The cosine wavelet transform corresponding to the A-wavelet transform [14] for () is given by The corresponding plots in 3D and 2D are shown in Figures 3(c) and 3(e), respectively.As discussed in the beginning of this section, the first component  0 () provides the same information as that obtained in [14] but along the primary curves and also the work [16] but in the first block along both primary and secondary curves.So, we have computed only the first component  0 ().The graphs in Figure 3(a) have sharper peaks in the time-frequency plane than that of Figure 3(c).Comparing 2D plots, the frequencies are better resolved in Figure 3(d) than that in Figure 3(e).In Figure 3(e), the frequencies are either not resolved at all or are poorly resolved.Figure 3(b) has sharper peaks in the time-frequency plane than that of Figure 3(a) and better resolves the frequency contents.
We compute the values of  ,,, (,  (+)+ ), sine wavelet transform, developed in the current paper for the International Journal of Digital Multimedia Broadcasting localization is better than both Figures 4(a) and 4(c).Using the proposed work, we can resolve the very close frequencies by taking appropriate block.Important unitary transforms have been discussed in [15].The results obtained in this paper can be applied to them too.

Conclusion
In this paper, we have proposed a new representation of the Fourier transform,  * -wavelet transform, which provides better frequency localization than that of A-wavelet transform.The A-wavelet transform is a particular case of the  *wavelet transform that provides the signal information along the primary curves, which are separated out by 2 in the timefrequency plane.The proposed work can provide the signal information along the secondary curves and the separation between two curves can be made arbitrarily small, which is not possible in the A-wavelet transform.In this work, we can theoretically resolve the frequency contents as small as we please.This study can be useful in many areas such as image processing and computer vision.

Figure 2 :
Figure 2: Loci of points in time-frequency plane of S-wavelet transform for  = 4.

Table 1 :
Fourier transform computation for the frequency  = 1.75 of a function defined in time duration (−10, 10).