Journal of Applied Mathematics

Volume 2014 (2014), Article ID 526940, 6 pages

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

Research Article

Expansions of Functions Based on Rational Orthogonal Basis with Nonnegative Instantaneous Frequencies

Xiaona Cui and Suxia Yao

College of Mathematics and Information Science, Henan Normal University, Xinxiang 453007, China

Received 16 March 2014; Revised 22 July 2014; Accepted 5 August 2014; Published 21 August 2014

Academic Editor: Didier Felbacq

Copyright © 2014 Xiaona Cui and Suxia Yao. 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.

We consider in this paper expansions of functions based on the rational orthogonal basis for the space of square integrable functions. The basis functions have nonnegative instantaneous frequencies so that the expansions make physical sense. We discuss the almost everywhere convergence of the expansions and develop a fast algorithm for computing the coefficients arising in the expansions by combining the characterization of the coefficients with the fast Fourier transform.

1. Introduction

A common approach for understanding a function is to expand it as a sum of basic functions. These basic functions should be chosen for different practical purpose. In the time-frequency analysis, the instantaneous frequency is one of the most important information for understanding a given signal. Then the basic functions used to express a complicated signal are expected to have nonnegative instantaneous frequencies so that the expansion of a signal makes physical sense. In this point of view, recent works have contributed to characterize and construct the basic functions [1–7] and establish numerical algorithms for decomposing a signal based on the basic functions [8–10].

Motivated by the analytic signal approach [11, 12], the authors in [7] constructed rational orthogonal basis for , the space of square integrable functions on , with each element of the basis enjoying physically meaningful instantaneous frequency. Specifically, let denote the set of all positive integers. The sequence , , forms an orthogonal basis for , where for each

and for each , . As pointed out in [7], the basis functions all have nonnegative instantaneous frequencies in the sense of the analytic signal. Based upon this basis, we focus in this paper on decomposing a function in into a sum of the basis functions , . For simplicity, we first transform the basis into an orthonormal one by normalization. A direct computation shows that there holds, for each ,

Hence, we obtain the resulting orthonormal basis for , where

In fact, due to the completeness of in , each can be represented as the infinite series

where the coefficients are defined by

However, for practical purpose, one often needs to decompose a function into a finite sum of basic functions. Thus, we have to approximate a function by the partial sum of the series (4) defined by

From this point of view, we have to deal with two issues. The first one is to study the convergence of the partial sum (6) as tends to infinity. Since it is known that the sequence of the partial sums converges to in the sense of norm, we will aim in this paper at discussing the almost everywhere convergence of them. The second one is to develop a fast algorithm for computing the partial sum so that this type of decomposition can be carried out in practical applications.

The organization of this paper is as follows. In Section 2, we establish the relation between the coefficients in (5) and the classical Fourier coefficients and then discuss the almost everywhere convergence of the sequence of the partial sums in (6). We develop in Section 3 an algorithm for computing the coefficients in the expansion (6) by employing the relation between the coefficients and the classical Fourier coefficients and the fast Fourier transform. To show that the proposed algorithm is fast, we also give an estimate of its computational complexity.

2. Almost Everywhere Convergence

We discuss in this section the almost everywhere convergence of the sequence , . This will be done by making use of some classical results about the Fourier series. To this end, with respect to any , , a periodic function is introduced as follows.

Let be defined by

Then , means is a one-to-one mapping from onto . Then for any , , we set

where . The following lemma gives a desired properties of the functions and .

Lemma 1. Let . If satisfies

then .

Proof. Taking the change of variables , we have

Similarly, for the function , we get

Together with the assumption (9), the above equalities lead to the desired conclusions.

We need one more lemma to represent the coefficients , , in terms of the classical Fourier coefficients of . For any , the Fourier coefficients of are defined by

Lemma 2. Suppose that satisfies

and and are defined in (8). Then the coefficients , defined in (5) are determined by

or

Proof. For simplicity, we only give the proof for the cases of and , . The proof for the other cases can be similarly handled. By definition we have for each

Taking the change of variables , we get

where we used the fact that . Combining the assumptions with Lemma 1, we point out that . Then the above equation leads to

In the same method, for the case of , , and by the fact that , we obtain

According to the relationship between and , we have another form of , which completes the proof.

It is known that for any , its Fourier coefficients tend to zero as . This is the well-known Riemann-Lebesgue lemma [13]. As a direct consequence of Lemma 2, we have the same decay property of the coefficients .

Corollary 3. If satisfies (13) then there holds

Proof. It follows from Lemma 1 that . Thus, there holds , . The relation between , and , , which is given in Lemma 2, yields the desired result.

We now turn to the main theorem in this section about the almost everywhere convergence of the sequence , . The almost everywhere convergence of the Fourier series is one of the most important results in the theory of Fourier series [14, 15]. It is said that for any , , the partial sum of its Fourier series

converges to almost everywhere on as . We will show that the sequence of , , enjoys the same convergence property if satisfies condition (13).

Theorem 4. If satisfies (13), then the sequence , , converges to almost everywhere on .

Proof. Associated with any satisfying (13), due to Lemma 1, we know that ; it suffices to prove that as there holds almost everywhere on

Consequently, the partial sum of its Fourier series tends to almost everywhere on as . Upon this observation, (22) can be obtained by proving that as ,

To get the difference between and , we need to express the two items, respectively. By definition and change of variables , we have

Applying Lemma 2 to the above equation, we can reexpress through a direct computation. Specifically, upon introducing four summations

we rewrite as follows:

To be compared with , the item should be represented in terms of the Fourier coefficients of . For this purpose, we rewrite , , as follows:

Substituting the above formula into the partial sum of the Fourier series of , we have

In a similar manner, we split the above equation into five items as follows:

where

It is clear that , . For , there holds

Then by subtracting (29) from (26), we obtain

Together with the fact that , Riemann-Lebesgue lemma yields for any

3. Numerical Algorithm

Since the expansion (6) can serve as a representation of with both mathematical and physical meanings, developing a fast approach for computing the coefficients , is practically important. Lemma 2 states that evaluating the expansion coefficients of can be carried by calculating the Fourier coefficients of , and then the well-known fast Fourier transform [16, 17] can be applied to compute the latter ones.

We first recall the fast Fourier transform. The algorithm is developed to reduce the computational cost of the discrete Fourier transform which is defined as follows:

where . It requires number of multiplications to compute , directly by formula (34). When is large, the computational cost is huge. The key idea of the fast Fourier transform is to compute by splitting the sum into two parts according to the even and odd indices. More precisely, assume that is an positive even integer, that is, . We rewrite formula (34) as follows:

where

Note that there holds the relations and , for . Thus, one can obtain the same result for half of the cost by the following three steps.

Step 1. For , compute and .

Step 2. Compute , by .

Step 3. Compute , by .

Assume that is a power of 2. It is clear that and are still two discrete Fourier transform of order . Then the above process can be iterated until we arrive at the discrete Fourier transform of order 2. Consequently, we get the celebrated fast Fourier transform, which brings the computational cost for computing , down to .

Applying the fast Fourier transform to compute the Fourier coefficients which are used to represent the coefficients , we have the following fast algorithm for computing , .

Algorithm 5. Let be a power of 2.

Step 1. Compute , by applying the fast Fourier transform to , . Set , and .

Step 2. Compute , according to formula .

Step 3. Compute , according to formula .

Step 4. Compute , according to formula .

The output , of Algorithm 5 is regarded as an approximation of the coefficients . The number of multiplications used in Algorithm 5 is estimated in the next theorem.

Theorem 6. The total number of multiplications needed for Algorithm 5 is .

Proof. The number of multiplications of Step 1, used to compute the discrete Fourier transform of order by the fast Fourier transform, is . In the other steps, the number of multiplications for evaluating a family of can be obtained as . We use to denote the number of multiplications used in Algorithm 5. Then we have

which gives the desired result.

It is noticed that the change of variable could impact the regularity of the related function in (8), thus impacting the accuracy of the numerical scheme. Therefore we require the conditions of Lemma 2 to keep the boundedness of the coefficients , . Other kinds of change of variable should be tried in the future study.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This research is supported by the Natural Science Foundation of China under Grants 11171092 and 41374158.

References

N. E. Huang, Z. Shen, S. R. Long et al., “The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis,” The Royal Society of London. Proceedings A: Mathematical, Physical and Engineering Sciences, vol. 454, no. 1971, pp. 903–995, 1998.

Y. Liu and Y. Xu, “Piecewise linear spectral sequences,” Proceedings of the American Mathematical Society, vol. 133, no. 8, pp. 2297–2308, 2005.

T. Qian, “Characterization of boundary values of functions in Hardy spaces with applications in signal analysis,” Journal of Integral Equations and Applications, vol. 17, no. 2, pp. 159–198, 2005.

T. Qian, “Mono-components for decomposition of signals,” Mathematical Methods in the Applied Sciences, vol. 29, no. 10, pp. 1187–1198, 2006.

T. Qian, R. Wang, Y. Xu, and H. Zhang, “Orthonormal bases with nonlinear phases,” Advances in Computational Mathematics, vol. 33, no. 1, pp. 75–95, 2010.

R. C. Sharpley and V. Vatchev, “Analysis of the intrinsic mode functions,” Constructive Approximation, vol. 24, no. 1, pp. 17–47, 2006.

L. Tan, L. Shen, and L. Yang, “Rational orthogonal bases satisfying the Bedrosian identity,” Advances in Computational Mathematics, vol. 33, no. 3, pp. 285–303, 2010.

T. Qian, L. Zhang, and Z. Li, “Algorithm of adaptive Fourier decomposition,” IEEE Transactions on Signal Processing, vol. 59, no. 12, pp. 5899–5906, 2011.

T. Qian and Y. Wang, “Adaptive Fourier series—a variation of greedy algorithm,” Advances in Computational Mathematics, vol. 34, no. 3, pp. 279–293, 2011.

R. Wang, Y. Xu, and H. Zhang, “Fast nonlinear fourier expansions,” Advances in Adaptive Data Analysis: Theory and Applications, vol. 1, no. 3, pp. 373–405, 2009.

L. Cohen, Time-Frequency Analysis, Prentice-Hall, Englewood Cliffs, NJ, USA, 1995.

D. Gabor, “Theory of communication,” Journal of Electrical Engineering, vol. 93, pp. 426–457, 1946.

A. Zygmund, Trigonometric Series, Cambridge University Press, Cambridge, Mass, USA, 1959.

L. Carleson, “On convergence and growth of partial sums of Fourier series,” Acta Mathematica, vol. 116, pp. 135–157, 1966.

R. A. Hunt, “On the convergence of fourier series,” in Proceedings of the Southern Illinois University Conference, Southern Illinois University, Edwardsville, Ill, USA, 1967.

C. Gasquet and P. Witomski, Fourier Analysis and Applications, Springer, New York, NY. USA, 1999.

D. Kincaid and W. Cheney, Numerical Analysis, Brooks/Cole, Pacific Grove, Calif, USA, 2nd edition, 1996.