Nonlinear Vibrations, Stability Analysis and ControlView this Special Issue
Research Article | Open Access
Shannon Wavelets for the Solution of Integrodifferential Equations
Shannon wavelets are used to define a method for the solution of integrodifferential equations. This method is based on (1) the Galerking method, (2) the Shannon wavelet representation, (3) the decorrelation of the generalized Shannon sampling theorem, and (4) the definition of connection coefficients. The Shannon sampling theorem is considered in a more general approach suitable for analysing functions ranging in multifrequency bands. This generalization coincides with the Shannon wavelet reconstruction of functions. Shannon wavelets are -functions and their any order derivatives can be analytically defined by some kind of a finite hypergeometric series (connection coefficients).
In recent years wavelets have been successfully applied to the wavelet representation of integro-differential operators, thus giving rise to the so-called wavelet solutions of PDE and integral equations. While wavelet solutions of PDEs can be easily find in a large specific literature, the wavelet representation of integro-differential operators cannot be considered completely achieved and only few papers discuss in depth this question with particular regards to methods for the integral equations. Some of them refer to the Haar wavelets [1–3] to the harmonic wavelets [4–9] and to the spline-Shannon wavelets [10–13]. These methods are mainly based on the Petrov-Galerkin method with a suitable choice of the collocation points . Alternatively to the collocation method, there has been also proposed, for the solution of PDEs, the evaluation of the differential operators on the wavelet basis, thus defining the so-called connection coefficients [6, 15–21].
Wavelets  are localized functions which are a useful tool in many different applications: signal analysis, data compression, operator analysis, PDE solving (see, e.g., [15, 23] and references therein), vibration analysis, and solid mechanics . Very often wavelets have been used only as any other kind of orthogonal functions, without taking into consideration their fundamental properties. The main feature of wavelets is, in fact, their possibility to split objects into different scale components [22, 23] according to the multiscale resolution analysis. For the functions, that is, functions with decay to infinity, wavelets give the best approximation. When the function is localized in space, that is, the bottom length of the function is within a short interval (function with a compact support), such as pulses, any other reconstruction, but wavelets, leads towards undesirable problems such as the Gibbs phenomenon when the approximation is made in the Fourier basis. Wavelets are the most expedient basis for the analysis of impulse functions (pulses) [24, 25].
Among the many families of wavelets, Shannon wavelets  offer some more specific advantages, which are often missing in the others. In fact, Shannon wavelets(1)are analytically defined; (2)are infinitely differentiable;(3)are sharply bounded in the frequency domain, thus allowing a decomposition of frequencies in narrow bands;(4)enjoy a generalization of the Shannon sampling theorem, which extend to all range of frequencies (5)give rise to the connection coefficients which can be analytically defined [15–17] for any order derivatives, while for the other wavelet families they were computed only numerically and only for the lower order derivatives [18, 19, 21].
The (Shannon wavelet) connection coefficients are obtained in  as a finite series (for any order derivatives). In Latto's method [18, 20, 21], instead, these coefficients were obtained only (for the Daubechies wavelets) by using the inclusion axiom but in approximated form and only for the first two-order derivatives. The knowledge of the derivatives of the basis enables us to approximate a function and its derivatives and it is an expedient tool for the projection of differential operators in the numerical computation of the solution of both partial and ordinary differential equations [6, 15, 23, 26].
The wavelet reconstruction by using Shannon wavelets is also a fundamental step in the analysis of functions-operators. In fact, due to their definition Shannon wavelets are box functions in the frequency domain, thus allowing a sharp decorrelation of frequencies, which is an important feature in many physical-engineering applications. In fact, the reconstruction by Shannon wavelets ranges in multifrequency bands. Comparing with the Shannon sampling theorem where the frequency band is only one, the reconstruction by Shannon wavelets can be done for functions ranging in all frequency bands (see, e.g., ). The Shannon sampling theorem , which plays a fundamental role in signal analysis and applications, will be generalized, so that under suitable hypotheses a few set of values (samples) and a preliminary chosen Shannon wavelet basis enable us to completely represent, by the wavelet coefficients, the continuous signal and its frequencies.
The Shannon wavelet solution of an integrodifferential equation (with functions localized in space and slow decay in frequency) will be computed by using the Petrov-Galerkin method and the connection coefficients. The wavelet coefficients enable to represent the solution in the frequency domain singling out the contribution to different frequencies.
This paper is organized as follows. Section 2 deals with some preliminary remarks and properties of Shannon wavelets also in frequency domain; the reconstruction of a function is given in Section 3 together with the generalization of the Shannon sampling theorem; the error of the wavelet approximation is computed. The wavelet reconstruction of the derivatives of the basis and the connection coefficients are given in Section 4. Section 5 deals with the Shannon wavelet solution of an integrodifferential equation and an example is given at last in Section 6.
2. Shannon Wavelets
From these functions a multiscale analysis  can be derived. The dilated and translated instances, depending on the scaling parameter and space shift , are
2.1. Properties of the Shannon Scaling and Wavelet Functions
By a direct computation it can be easily seen that
with Kroneker symbol, so that
It is also
where either the scaling functions or the wavelet vanishes, but it is important to notice that when the scaling function is zero, the wavelet is not and viceversa. As we shall see later, this property will simplify the numerical methods based on collocation point.
it can be also easily checked that for a fixed
Analogously we have
The maximum and minimum values of these functions can be easily computed. The maximum value of the scaling function can be found in correspondence of
The min value of can be computed only numerically and it is
The minimum of the wavelet can be found in correspondence of the middle point of the zeroes (2.7) so that
and the max values of are
2.2. Shannon Wavelets Theory in the Fourier Domain
be the Fourier transform of the function , and
its inverse transform.
Analogously for the dilated and translated instances of scaling/wavelet function, in the frequency domain, it is
It can be seen that
so that by using the function and there is a decorrelation into different non-overlapping frequency bands.
For each and , the inner product is defined as
which, according to the Parseval equality, can be expressed also as
where the bar stands for the complex conjugate.
Theorem 2.1. Shannon wavelets are orthonormal functions, in the sense that With being the Kroenecker symbols.
Theorem 2.2. The translated instances of the Shannon scaling functions , at the level , are orthogonal, in the sense that being .
See the proof in .
Theorem 2.3. The translated instances of the Shannon scaling functions , at the level , are orthogonal to the Shannon wavelets, in the sense that being .
Proof is in .
3. Reconstruction of a Function by Shannon Wavelets
Let be a function such that for any value of the parameters , it is
and the Paley-Wiener space, that is, the space of band limited functions, that is,
According to the sampling theorem (see, e.g.,  and references therein) we have the following.
Theorem 3.1 (Shannon). If and the series uniformly converges to , and
Proof. In order to compute the values of the coefficients we have to evaluate the series in correspondence of the integer:
having taken into account (2.5).
The convergence follows from the hypotheses on . In particular, the importance of the band limited frequency can be easily seen by applying the Fourier transform to (3.3): so that In other words, if the function is band limited (i.e., with compact support in the frequency domain), it can be completely reconstructed by a discrete Fourier series. The Fourier coefficients are the values of the function sampled at the integers.
As a generalization of the Paley-Wiener space, and in order to generalize the Shannon theorem to unbounded intervals, we define the space of functions such that the integrals
For the unbounded interval, let us prove the following.
Proof. The representation (3.11) follows from the orthogonality of the scaling and Shannon wavelets (Theorems 2.1, 2.2, and 2.3). The coefficients, which exist and are finite, are given by (3.8). The convergence of the series is a consequence of the wavelet axioms.
It should be noticed that
so that for a band limited frequency signal, that is, for a signal whose frequency belongs to the band , this theorem reduces to the Shannon sampling theorem. More in general, the representation (3.11) takes into account more frequencies ranging in different bands. In this case we have some nontrivial contributions to the series coefficients from all bands, ranging from :
In the frequency domain, (3.11) gives
3.1. Error of the Shannon Wavelet Approximation
Let us fix an upper bound for the series of (3.11) in a such way that we can only have the approximation
This approximation can be estimated by the following
Theorem 3.3 (Error of the Shannon wavelet approximation). The error of the approximation (3.19) is given by
4. Reconstruction of the Derivatives
Let and let be a differentiable function with sufficiently high. The reconstruction of a function given by (3.11) enables us to compute also its derivatives in terms of the wavelet decomposition:
so that, according to (3.11), the derivatives of are known when the derivatives
Indeed, in order to represent differential operators in wavelet bases, we have to compute the wavelet decomposition of the derivatives:
Their computation can be easily performed in the Fourier domain, thanks to the equality (2.27). In fact, in the Fourier domain the -order derivative of the (scaling) wavelet functions is
and according to (2.24),
Taking into account (2.27), we can easily compute the connection coefficients in the frequency domain
with the derivatives given by (4.6).
Theorem 4.1. The any order connection coefficients (4.4)1 of the Shannon scaling functions are
For the proof see .
Analogously for the connection coefficients (4.4)2 we have the following.
Theorem 4.2. The any order connection coefficients (4.7)2 of the Shannon scaling wavelets are respectively, for , and .
For the proof see .
Theorem 4.3. The connection coefficients are recursively given by the matrix at the lowest scale level:
Moreover it is
If we consider a dyadic discretisation of the -axis such that
For instance, in
Analogously it is
from where, in , it is
5. Wavelet Solution of the Integrodifferential Equation
Let us consider the following linear integrodifferential equation:
which includes as special cases the integral equation () and the differential equation (). When there is the trivial solution .
It is assumed that the kernel is in the form:
and the given functions , , , so that, according to (3.11)
with the wavelet coefficients given by (3.8).
The analytical solution of (5.1) can be obtained as follows.
Although the existence of solution is proven, the computation of the Fourier transform could not be easily performed. Therefore the numerical computation is searched in the wavelet approximation.
The wavelet solution of (5.1) can be obtained as follows: it is assumed that the unknown function and its derivative can be written as
and the integral can be written as
There follows the system
and, according to the definition of the connection coefficients,
Analogously, it is
or, according to (4.11), and rearranging the indices
and up to a fixed scale of approximation :
Let us consider the following equation:
with the condition
The analytical solution, as can be directly checked, is
At the level of approximation , from (5.3) we have
System (5.18) becomes
and, since and , according to (6.6) we have
whose solution is
In this paper the theory of Shannon wavelets combined with the connection coefficients methods and the Petrov-Galerkin method has been used to find the wavelet approximation of integrodifferential equations. Among the main advantages there is the decorrelation of frequencies, in the sense that the differential operator is splitted into its different frequency bands.
- H.-T. Shim and C.-H. Park , “An approximate solution of an integral equation by wavelets,” Journal of Applied Mathematics and Computing, vol. 17, no. 1-2-3, pp. 709–717, 2005.
- U. Lepik, “Numerical solution of evolution equations by the Haar wavelet method,” Applied Mathematics and Computation, vol. 185, no. 1, pp. 695–704, 2007.
- U. Lepik, “Solving fractional integral equations by the Haar wavelet method,” Applied Mathematics and Computation, vol. 214, no. 2, pp. 468–478, 2009.
- C. Cattani and A. Kudreyko, “Application of periodized harmonic wavelets towards solution of egenvalue problems for integral equations,” Mathematical Problems in Engineering, vol. 2010, Article ID 570136, 8 pages, 2010.
- C. Cattani and A. Kudreyko, “Harmonic wavelet method towards solution of the Fredholm type integral equations of the second kind,” Applied Mathematics and Computation, vol. 215, no. 12, pp. 4164–4171, 2010.
- S. V. Muniandy and I. M. Moroz, “Galerkin modelling of the Burgers equation using harmonic wavelets,” Physics Letters A, vol. 235, no. 4, pp. 352–356, 1997.
- D. E. Newland, “Harmonic wavelet analysis,” Proceedings of the Royal Society of London A, vol. 443, pp. 203–222, 1993.
- J.-Y. Xiao, L.-H. Wen, and D. Zhang, “Solving second kind Fredholm integral equations by periodic wavelet Galerkin method,” Applied Mathematics and Computation, vol. 175, no. 1, pp. 508–518, 2006.
- S. Yousefi and A. Banifatemi, “Numerical solution of Fredholm integral equations by using CAS wavelets,” Applied Mathematics and Computation, vol. 183, no. 1, pp. 458–463, 2006.
- Y. Mahmoudi, “Wavelet Galerkin method for numerical solution of nonlinear integral equation,” Applied Mathematics and Computation, vol. 167, no. 2, pp. 1119–1129, 2005.
- K. Maleknejad and T. Lotfi, “Expansion method for linear integral equations by cardinal B-spline wavelet and Shannon wavelet as bases for obtain Galerkin system,” Applied Mathematics and Computation, vol. 175, no. 1, pp. 347–355, 2006.
- K. Maleknejad, M. Rabbani, N. Aghazadeh, and M. Karami, “A wavelet Petrov-Galerkin method for solving integro-differential equations,” International Journal of Computer Mathematics, vol. 86, no. 9, pp. 1572–1590, 2009.
- A. Mohsen and M. El-Gamel, “A sinc-collocation method for the linear Fredholm integro-differential equations,” ZAMP, vol. 58, no. 3, pp. 380–390, 2007.
- N. Bellomo, B. Lods, R. Revelli, and L. Ridolfi, Generalized Collocation Methods: Solutions to Nonlinear Problem, Modeling and Simulation in Science, Engineering and Technology, Birkhäuser, Boston, Mass, USA, 2008.
- C. Cattani, “Harmonic wavelets towards the solution of nonlinear PDE,” Computers & Mathematics with Applications, vol. 50, no. 8-9, pp. 1191–1210, 2005.
- C. Cattani, “Connection coefficients of Shannon wavelets,” Mathematical Modelling and Analysis, vol. 11, no. 2, pp. 117–132, 2006.
- C. Cattani, “Shannon wavelets theory,” Mathematical Problems in Engineering, vol. 2008, Article ID 164808, 24 pages, 2008.
- A. Latto, H. L. Resnikoff, and E. Tenenbaum, “The evaluation of connection coefficients of compactly supported wavelets,” in Proceedings of the French-USA Workshop on Wavelets and Turbulence, Y. Maday, Ed., pp. 76–89, Springer, New York, NY, USA, June 1992.
- E. B. Lin and X. Zhou, “Connection coefficients on an interval and wavelet solutions of Burgers equation,” Journal of Computational and Applied Mathematics, vol. 135, no. 1, pp. 63–78, 2001.
- J. Restrepo and G. K. Leaf, “Wavelet-Galerkin discretization of hyperbolic equations,” Journal of Computational Physics, vol. 122, no. 1, pp. 118–128, 1995.
- C. H. Romine and B. W. Peyton, “Computing connection coefficients of compactly supported wavelets on bounded intervals,” Tech. Rep. ORNL/TM-13413, Oak Ridge, Computer Science and Mathematical Division, Mathematical Sciences Section, Oak Ridge National Laboratory, 1997, http://citeseer.ist.psu.edu/romine97computing.html.
- I. Daubechies, Ten Lectures on Wavelets, vol. 61 of CBMS-NSF Regional Conference Series in Applied Mathematics, SIAM, Philadelphia, Pa, USA, 1992.
- C. Cattani and J. Rushchitsky, Wavelet and Wave Analysis as Applied to Materials with Micro or Nanostructure, vol. 74 of Series on Advances in Mathematics for Applied Sciences, World Scientific, Singapore, 2007.
- E. Bakhoum and C. Toma, “Mathematical transform of travelling-wave equations and phase aspects of quantum interaction,” Mathematical Problems in Engineering, vol. 2010, Article ID 695208, 15 pages, 2010.
- G. Toma, “Specific differential equations for generating pulse sequences,” Mathematical Problems in Engineering, vol. 2010, Article ID 324818, 11 pages, 2010.
- C. Cattani, “Harmonic wavelet solutions of the schrödinger equation,” International Journal of Fluid Mechanics Research, vol. 5, pp. 1–10, 2003.
- S. Unser, “Sampling-50 years after Shannon,” Proceedings of the IEEE, vol. 88, no. 4, pp. 569–587, 2000.
- C. Cattani, “Shannon wavelet analysis,” in Proceedings of the International Conference on Computational Science (ICCS '07), Y. Shi et al., Ed., vol. 4488 of Lecture Notes in Computer Science Part II, pp. 982–989, Springer, Beijing, China, May 2007.
- E. Deriaz, “Shannon wavelet approximation of linear differential operators,” Institute of Mathematics of the Polish Academy of Sciences, no. 676, 2007.
Copyright © 2010 Carlo Cattani. 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.