# Efficient time-frequency localization of a signal.

1. IntroductionThere 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-8]. In the last two decades and so, large numbers of articles have appeared on the wavelet transforms in literature [9-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 2[pi] have been considered to develop the cosine- and sine-type 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[pi]/[omega]. In other words, this provides the signal information along the curves that are separated out by 2[pi]/[omega] 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[pi]/[omega] distance. The proposed representation, named as [A.sup.*]-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 [A.sup.*]-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.

2. 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 [L.sub.2](R) and nonzero only in the region of interest. Here [L.sub.2](R) refers to the space of square-integrable functions. The STFT of a signal f(t), denoted by F(t,[omega]), using the window function g(t) and reconstruction of the original function f(t) from its STFT are given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (1)

The time-window function g(t) 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 W(a,b), of an arbitrary function f(t) using the wavelet function [psi](t) is given by

W(a,b) = 1/[square root of (a)] [[integral].sup.[infinity].sub.-[infinity]] f(t)[psi](t - b/a)dt for a > 0, b [member of] (-[infinity], [infinity]), (2)

where a and b 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 sine-type wavelet transforms. The cosine- and sine-type wavelet transforms are based on the truncated cosinusoidal and sinusoidal signals, denoted by [psi](t) and [phi](t), respectively, and they are defined over a time interval of length 2[pi]; that is,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (3a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (3b)

Denote [[psi].sub.[omega],b] (t) = [psi]([omega][t - b]) and [[phi].sub.[omega],b](t) = [phi]([omega][t - b]) for t [member of] (-[infinity], [infinity]) and define [T.sub.[psi]] ([omega], [b.sub.n]) and [T.sub.[phi]]([omega], [b.sub.n]) using [[psi].sub.[omega],b](t) and [[phi].sub.[omega],b](t); that is,

[T.sub.[psi]]([omega], b) = [[integral].sup.[infinity].sub.-[infinity]] f (t) [[psi].sub.[omega],b](t) dt, [T.sub.[phi]]([omega], b) = [[integral].sup.[infinity].sub.-[infinity]] f (t) [[phi].sub.[omega],b](t) dt. (4)

The functions T[psi]([omega],[b.sub.n]) and T[phi]([omega],[b.sub.n]) are called cosine-wavelet (C-wavelet) and sine-wavelet (S-wavelet), respectively. The A-wavelet transform representation of the Fourier transform F([omega]) is given by

F([omega])= [[infinity].summation over (n=-[infinity])] [T.sub.[psi]]([omega], [b.sub.n]) - j [[infinity].summation over (n=-[infinity])] T[phi]([omega], [b.sub.n]), (5)

where j is an imaginary quantity.

This representation (i.e., C-wavelet and S-wavelet) reduces the computations because, for a given frequency [omega], the functions [T.sub.[psi]]([omega], [b.sub.n]) and [T.sub.[phi]]([omega], [b.sub.n]) are evaluated at the centers of the prespecified time intervals of length 2[pi]/[omega] in time parameter. The frequency and time parameters are related by the relation [omega]t = 2[pi]. The number of time intervals on time-axis is given by 2N([omega]) + 1 (i.e., -N([omega]),-N([omega]) + 1, ..., -1, 0, 1, ..., N([omega]), N([omega]) + 1), where N([omega]) = 0 for [omega] [member of] [0, 1/3], and for [omega] > 1/3, N([omega]) = [??](3[omega] -1)/2[??]. For details on N([omega]) refer to [14]. The A-wavelet transform provides the signal information in the time-frequency plane along the curves [omega]t = 2n[pi], n = 0, [+ or -] 1, [+ or -] 2,.... In the next section, we discuss [A.sup.*]-wavelet transform, an extension of A-wavelet transform, that provides the signal information along the curves separated by less 2[pi] in the time-frequency plane.

3. [A.sup.*]-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[pi] (refer to (3a) and (3b)). The A-wavelet transform provides the signal information along the curves [omega]t = 2n[pi], where n is an integer, that is, all such curves that are separated by 2[pi] in the time-frequency plane. We call the curves [omega]t = 2n[pi], n = 0, [+ or -] 1, [+ or -] 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[pi]. The interval length is parameterized by introducing two new integer variables n and m. Define the functions [[psi].sub.m,n](t) and [[phi].sub.m,n](t) as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (6a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (6b)

We look at different values of m and n. For m = 0, [[psi].sub.m,n](t) is the Dirac Delta function and [[phi].sub.m,n](t) is identically zero and the outcome is simply a sampled signal. For m = [2.sup.n], the functions [[psi].sub.m,n](t) and [[phi].sub.m,n](t) are exactly the same as the [psi](t) and [phi](t) functions, respectively, which are given in (3a) and (3b). For m > [2.sup.n], the functions [[psi].sub.m,n](t) and [[phi].sub.m,n](t) 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 m [less than or equal to] [2.sup.n]. Without loss of generality, we may assume that m and n assume nonnegative integer values. For negative m, the same discussion will hold as that for positive m, but on the negative time-axis. Negative n will increase the number of time intervals over which the functions [[psi].sub.m,n](t) and [[phi].sub.m,n](t) have nonzero values. The condition m < [2.sup.n] makes the time interval length smaller over which the functions [[psi].sub.m,n](t) and [[phi].sub.m,n](t) 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 {[[psi].sub.m,n,[omega],k](t), [[phi].sub.m,n,[omega],k](t)} from [[psi].sub.m,n](t) and [[phi].sub.m,n](t) by using the time-scaling and shift transformation; that is, t [right arrow] [omega]t and t [right arrow] (t - [b.sub.k]); that is,

[[psi].sub.m,n,[omega],k](t) = [[psi].sub.m,n]([omega](t - [b.sub.k])), [[phi].sub.m,n,[omega],k](t) = [[phi].sub.m,n]([omega](t - [b.sub.k])), for t [member of] (-[infinity], [infinity]), (7)

where the frequency variable [omega] assumes real values and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. In general, k is infinite; however, for finite support signals, its value is finite.

Consider a function f(t) that satisfies all the conditions for existence of Fourier transform. Its Fourier transform F([omega]) is given by

F([omega]) = [[integral].sup.[infinity].sub.-[infinity]] f(t)[e.sup.-j[omega]t] dt. (8)

For [omega] = 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 [omega] [not equal to] 0 frequencies. We write F([omega]) that is defined in (8) as follows:

F([omega]) = [[integral].sup.[infinity].sub.-[infinity]] f(t) cos ([omega]t) dt - j [[integral].sup.[infinity].sub.- [infinity]] f(t) sin ([omega]t) dt. (9)

We try to write the integrals in (9) in terms of [[psi].sub.m,n,[omega],k](t) and [[phi].sub.m,n,[omega],k](t) functions that are defined in (7). For this purpose, we uniformly divide the time-axis into disjoint intervals of length 2[pi]m/([2.sup.n][omega]) that are centered at 2[pi]mk/([2.sup.n][omega]), k = 0, [+ or -] 1, [+ or -] 2,.... Denote kth interval by [I.sub.mnk], which can be defined as follows:

[I.sub.m,n,k] = [(2k - 1)m[pi]/[2.sup.n][omega], (2k + 1) m[pi]/[2.sup.n][omega]), k = 0, [+ or -] 1, [+ or -] 2,.... (10)

In (10), k assumes all integer values in general; however, its maximum value is determined by the frequencies, denoting them by [[THETA].sub.k,[omega]], that need be analyzed. For example, k = 1 and [omega] = 3[pi] gives the frequency range [[THETA].sub.1,3[pi]] = [m/(3*[2.sup.n]),3m/(3* [2.sup.n])). The secondary curves are given by the equation [omega]t = 2[pi]mk/[2.sup.n], (m [not equal to] [2.sup.n]), where m is a positive integer and n is nonnegative integer. We can write the Fourier transform F([omega]) defined in 9) over the disjoint intervals [I.sub.m,n,k] as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (11)

Changing the limits of integration, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (12)

On simplifying, it gives

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (13)

which can further be written as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (14)

In (14) the summation contains the term [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. We can write [2.sub.n]/m = sm + r, where s is a nonnegative integer and r is an integer such that 0 [less than or equal to] r < m. Thus, we have [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Since k assumes all integer values including zero, the expression [e.sup.-j2[pi]k/(ms+r)] will have (sm + r) distinct values for all values of k. For positive and negative values of k, we will have same results. For positive k, we will have intervals along positive axis and for negative k, we have intervals on negative axis; henceforth we will confine k as nonnegative integer. For fixed m, increasing n will increase the value of s, but r will have only m distinct values. Increasing s will help in analyzing the high frequencies. Increasing m will increase the values of r (r [not equal to] 0) and that in turn will increase the secondary curves. It may be noticed that the curves for r = 0 correspond to the primary curves. Increasing n and m 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 m. 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.

We write (14) as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (15)

F([omega]) can be written as a sum of (sm + r) number of terms; that is,

F([omega])= [sm+r-1.summation over (q=0)] [F.sub.q]([omega])[e.sup.-j(2[pi]q/(sm+r))]. (16)

We can write k = (sm + r)p + q, 0 [less than or equal to] q < (sm + r). From (15), [F.sub.q]([omega]) is given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (17)

Using the cosine and sine functions in terms of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (refer to (6a), (6b), and (7)), we can write (17) as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (18)

where

[b.sub.(sm+r)p+q] = 2[pi]((sm + r)p + q)/(sm + r)[omega]. (19)

The integrals in (18), for a fixed value of q, are defined over ((sm + r)p + q)th interval of length 2[pi]m/([2.sup.n][omega])(= 2[pi]/(sm + r)[omega]) that is centered at [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Denote these integrals as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (20a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (20b)

At the point (0,0), [T.sub.m,n,[psi],q](0,0) = [[integral].sup.[infinity].sub.-[infinity]]f(t)dt and [T.sub.m,n,[phi],q](0,0) = 0.

Writing expressions given in (18) in terms of [T.sub.m,n,[psi],q]([omega],[b.sub.(sm+r)p+q]) and [T.sub.m,n,[phi],q]([omega],[b.sub.(sm+r)p+q]), we get

[F.sub.q]([omega]) = [[infinity].summation over (p=-[infinity])] [T.sub.m,n,[psi],q]([omega],[b.sub.(sm+r)p+q]) -j[T.sub.m,n,[phi],q]([omega],[b.sub.(sm+r)p+q]). (21)

Using [F.sub.q]([omega]) from (21) in (16) gives

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (22)

The terms in (22) are the weighted sums of different cosine- and sine- wavelet transforms and it is named as [A.sup.*]-wavelet transform, an extension of A-wavelet transform. The relation (22) is the main outcome of this paper.

4. 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 [F.sub.0]([omega]) in a better localized form using the current work. The signals [F.sub.i]([omega]) (i [not equal to] 0) give the information about the signal along the curves [omega][b.sub.m,n] = 2[pi]/(i + 1), i = 1,2,3, ...,(sm + r - 1) in the time-frequency plane. For i = 0, it gives information along the primary curves. The equation [omega][b.sub.m,n] = 2[pi]/i indeed represents a family of hyperbolic curves of second kind, one for each value of i, in the time-frequency plane. Maximum value of i or the number of hyperbolic curves is determined by the values of m and n. The set of functions {[[psi].sub.m,n,[omega],k](t), [[phi].sub.m,n,[omega],k](t)} 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 [T.sub.m,n,[psi],q]([omega],[b.sub.(sm+r)p+q]) and [T.sub.m,n,[phi],q]([omega],[b.sub.(sm+r)p+q]) (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 n 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[pi], 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[pi]. 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 [T.sub.m,n,[psi],q]([omega],[b.sub.(sm+r)p+q]) and [T.sub.m,n,[phi],q]([omega],[b.sub.(sm+r)p+q]) over the time interval [-[pi]m/[2.sup.n], [pi]m/[2.sup.n]) = [-[pi]/(sm + r),[pi]/(sm + r)). This can be analogous to a stack of rectangular blocks, each having time and scale (frequency) as its two sides. The zeroth block (s = 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 (s = 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 (s = 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 maybe noted that the segment in time parameter of a block is appropriately chosen from the support of the signal maintaining the length of sth block as t/[2.sup.s], where t is the signal support. The zeroth block is used to analyze the smallest frequencies starting from zero. The parameter m specifies distinct (type) curves in a block and r 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 s 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 t [member of] (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.

The representation (22) has (sm + r) components [F.sub.i]([omega]), i = 0,1, ..., (sm + r - 1). Depending upon the domain under consideration in the time-frequency plane, the points on the curves are considered. For example, evaluating the frequencies [omega] [member of] [0,1/12) in the time range [-5, 5), only straight line passing through the origin and parallel to the frequency-axis is needed and for [omega] = 1.0, four curves (i.e., four samples) are needed. Thus, we need to have the values of [T.sub.m,n,[psi],q]([omega],[b.sub.(sm+r)p+q]) and [T.sub.m,n,[phi],q]([omega],[b.sub.(sm+r)p+q]) at the points ([omega],2[pi]((sm + r)p + q)/((sm + r)[omega])),where m is a fixed integer; r = 0,1,2, ..., (m - 1); q = 0,1,2, ..., (ms + r - 1); and p = 0, [+ or -] 1, [+ or -] 2.... The set of points at which the signals [T.sub.m,n,[psi],q]([omega],[b.sub.(sm+r)p+q]) and [T.sub.m,n,[phi],q]([omega],[b.sub.(sm+r)p+q]) are computed is given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (23)

For s = 0 and r = 1, 22) gives

F([omega]) = [F.sub.0]([omega]). (24)

It is not difficult to show that the signal information given by [F.sub.0]([omega]) is the same as that obtained in [14]. For m = 4 and n = 4, the plot of the equation [omega][b.sub.m,n] = 2[pi]m/[2.sup.n] is shown in Figure 2. (Two consecutive curves are [pi]/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 the time-frequency plane for [omega] [member of] (0,5), [b.sub.n] [member of] (-5, 5), where [b.sub.n] signifies the shift in time and m = 4. In this figure, there are four types of curves, each corresponding to [F.sub.i]([omega]), i = 0,1,2,3. The dark black curves are the primary curves that correspond to first component [F.sub.0]([omega]) 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 [F.sub.1]([omega]), third component [F.sub.2]([omega]), and fourth component [F.sub.3]([omega]), respectively. These are the secondary curves which cannot be obtained using the A-wavelet transform.

5. Illustration

We will elucidate the above discussed theorywith an example. The function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] determines the number of points at which the values of the signals [T.sub.m,n,[psi],q]([omega],[b.sub.(sm+r)p+q]) and [T.sub.m,n,[phi],q]([omega],[b.sub.(sm+r)p+q]) for q = 0,1,2, ..., (sm + r - 1), are needed for the given frequency [omega]. Consider a time-limited signal defined in the time interval (-10,10) whose Fourier transform for the frequency [omega] = 1.75 is to be evaluated. For m = 4, s = 1, and r = 0, the maximum value of p, denoted by [p.sub.max], is given by [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Thus, the number of points at which the signals [T.sub.m,n,[psi],q]([omega],[b.sub.(sm+r)p+q]) and [T.sub.m,n,[phi],q]([omega],[b.sub.(sm+r)p+q]) are to be evaluated, for [omega] = 1.75, is 21, that is, 0, [+ or -] 1, [+ or -] 2, [+ or -] 3, ..., [+ or -] 10. Since m = 4, the Fourier transform contains four components, that is, [F.sup.0]([omega]), [F.sub.1]([omega]), [F.sub.2]([omega]), and [F.sub.3]([omega]), and since s = 1, we need to consider the curves corresponding to the first block. [F.sub.0]([omega]) needs be computed at the center of every fourth time interval in both directions (p = 0, [+ or -] 4, [+ or -] 8) starting from the first one (assuming the index of the first interval as 0), that is, ([omega], [b.sub.0]), ([omega], [+ or -] [b.sub.4]), and ([omega], [+ or -] [b.sub.8]), where [b.sub.i] represents the time (location) parameter; [F.sub.1]([omega]) is computed at the center of every fourth time interval in both directions starting from the second one, that is, ([omega], [+ or -] [b.sub.1]), ([omega], [+ or -] [b.sub.5]), and ([omega], [+ or -] [b.sub.9]); [F.sub.2]([omega]) is computed at the center of every fourth time interval in both directions starting from the third one, that is, ([omega], [+ or -] [b.sub.2]), ([omega], [+ or -] [b.sub.6]), and ([omega], [+ or -] [b.sub.10]); [F.sub.3]([omega]) is computed at the center of every fifth time interval in both directions starting from the fourth one, that is, ([omega], [+ or -] [b.sub.3]) and ([omega], [+ or -] [b.sub.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 A-wavelet transform [14] and SC- and SS-wavelet transforms [16]. We take the signal f(t) = cos(1.3t) + cos(2.5t) + cos(3.6t) + cos(4.8t) and evaluate [T.sub.m,n,[psi],q]([omega],[b.sub.(sm+r)p+q]) and [T.sub.m,n,[phi],q]([omega],[b.sub.(sm+r)p+q]), r = 0, s = 1 m = 4, and q = 0,1, ..., (sm + r - 1) using the current work. Substituting the values of f(t) and [[psi].sub.m,n,[omega],k](t) in (20a), we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (25)

where

[b.sub.(sm+r)p+q] = 2[pi]((sm+r)p+q)/(sm + r)[omega]. (26)

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 f(t) is given by

[T.sub.[psi]] ([omega],b) = [[integral].sup.[pi]/[omega].sub.-[pi]/[omega]] (cos(1.3t) + cos(2.5t) + cos(3.6t) + cos(4.8t)) x cos([omega](t - b)) dt. (27)

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 [F.sub.0]([omega]) 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 [F.sub.0]([omega]). 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 [T.sub.m,n,[phi],q]([omega],[b.sub.(sm+r)p+q]), sine wavelet transform, developed in the current paper for the same function f(t) = cos(1.3t) + cos(2.5t) + cos(3.6t) + cos(4.8t). Substituting the values of f(t) and [[phi].sub.m,n,[omega],k](t) in (20b), we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (28)

where

[b.sub.(sm+r)p+q] = 2[pi]((sm + r) p + q)/(sm + r)[omega]. (29)

The sine wavelet transform corresponding to the A-wavelet transform of f(t) is computed by

[T.sub.[phi]]([omega],b) = [[integral].sup.[pi]/[omega].sub.-[pi]/[omega]] (cos(1.3t) + cos(2.5t) + cos(3.6t) + cos(4.8t)) x sin ([omega](t - b)) dt. (30)

The corresponding 3D and 2D plots are shown in Figures 4(a) and 4(d), respectively. Here, also, the frequency localization is better than that of the A-wavelet transform as evident from Figures 4(a) and 4(c). In 2D plots, the frequencies are better resolved than that of the A-wavelet transform as shown in Figures 4(d) and 4(e). In Figure 4(b), the frequency 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.

6. Conclusion

In this paper, we have proposed a new representation of the Fourier transform, [A.sup.*]-wavelet transform, which provides better frequency localization than that of A-wavelet transform. The A-wavelet transform is a particular case of the [A.sup.*]-wavelet transform that provides the signal information along the primary curves, which are separated out by 2[pi] in the time-frequency 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.

http://dx.doi.org/10.1155/2014/529852

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

References

[1] D. Gabor, "Theory of communication," Journal of IEE, vol. 93, pp. 429-457, 1946.

[2] T. Edward, Discrete Wavelet Transforms: Theory and Implementation Stanford, Stanford University, Palo Alto, Calif, USA, 1991.

[3] A. K. Louis, D. Maass, and A. Rieder, Wavelets: Theory and Applications, John Wiley & Sons, New York, NY, USA, 1997

[4] C. S. Burrus, R. A. Gopinath, and H. Guo, Introduction to Wavelets and Wavelet Transforms: A Primer, Prentice-Hall, Upper Saddle River, NJ, USA, 1998.

[5] J. C. Goswami and A. K. Chan, Fundamentals of Wavelets: Theory, Algorithms, and Applications, John Wiley & Sons, New York, NY, USA, 1999.

[6] Y. Meyer, Wavelets: Algorithms and Applications, Society for Industrial and Applied Mathematics, Philadelphia, Pa, USA, 1993.

[7] I. Daubechies, Ten Lectures on Wavelets, Society for Industrial and Applied Mathematics, Philadelphia, Pa, USA, 1992.

[8] S. Mallat, A Wavelet Tour of Signal Processing, Academic Press, San Diego, Calif, USA, 2nd edition, 1998.

[9] M. Vetterli and C. Herley, "Wavelets and filter banks: Theory and design," IEEE Transactions on Signal Processing, vol. 40, no. 9, pp. 2207-2232, 1992.

[10] G. Strang, "Wavelets and dilation equations: a brief introduction," SIAM Review, vol. 31, no. 4, pp. 614-627, 1989.

[11] O. Rioul and M. Vetterli, "Wavelets and signal processing," IEEE Signal Processing Magazine, vol. 8, no. 4, pp. 14-38, 1991.

[12] I. W. Selesnick, "A higher density discrete wavelet transform," IEEE Transactions on Signal Processing, vol. 54, no. 8, pp. 3039-3048, 2006.

[13] C. K. Chui and W. He, "Compactly supported tight frames associated with refinable functions," Applied and Computational Harmonic Analysis, vol. 8, no. 3, pp. 293-319, 2000.

[14] A. M. Grigoryan, "Fourier transform representation by frequency-time wavelets," IEEE Transactions on Signal Processing, vol. 53, no. 7, pp. 2489-2497, 2005.

[15] A. M. Grigoryan and S. S. Agaian, Multidimensional Discrete Unitary Transforms: Representation, Partitioning and Algorithms, Marcel Dekker, New York, NY, USA, 2003.

[16] S. Chand, "SC- and SS-wavelet transforms," Signal Processing, vol. 89, no. 3, pp. 305-313, 2009.

Satish Chand

Division of Computer Engineering, Netaji Subhas Institute of Technology, Sector-3, Dwarka, New Delhi 110 078, India

Correspondence should be addressed to Satish Chand; schand86@hotmail.com

Received 12 February 2014; Revised 15 July 2014; Accepted 23 July 2014; Published 25 August 2014

Academic Editor: Ekram Khan

TABLE 1: Fourier transform computation for the frequency [omega] = 1.75 of a function defined in time duration (-10, 10). p/Component [F.sub.0] [F.sub.1] [F.sub.2] [F.sub.3] ([omega]) ([omega]) ([omega]) ([omega]) 0, [+ or -] 4, Yes [+ or -] 8 [+ or -] 1, Yes [+ or -] 5, [+ or -] 9 [+ or -] 2, Yes [+ or -] 6, [+ or -] 10 [+ or -] 3, Yes [+ or -] 7

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Chand, Satish |

Publication: | International Journal of Digital Multimedia Broadcasting |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 6065 |

Previous Article: | A model for video quality assessment considering packet loss for broadcast digital television coded in H.264. |

Next Article: | A new QIM-based watermarking method robust to gain attack. |

Topics: |