Next Article in Journal
A Smoothed Algorithm with Convergence Analysis under Generalized Maximum Correntropy Criteria in Impulsive Interference
Next Article in Special Issue
Wavelet-Based Entropy Measures to Characterize Two-Dimensional Fractional Brownian Fields
Previous Article in Journal
Quantum Games with Unawareness with Duopoly Problems in View
Previous Article in Special Issue
Classification of Heart Sounds Based on the Wavelet Fractal and Twin Support Vector Machine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Biorthogonal-Wavelet-Based Method for Numerical Solution of Volterra Integral Equations

Department of Mathematics and Statistics, College of Natural and Health Sciences, Zayed University, 144543 Abu Dhabi, UAE
Submission received: 6 October 2019 / Revised: 26 October 2019 / Accepted: 1 November 2019 / Published: 10 November 2019
(This article belongs to the Collection Wavelets, Fractals and Information Theory)

Abstract

:
Framelets theory has been well studied in many applications in image processing, data recovery and computational analysis due to the key properties of framelets such as sparse representation and accuracy in coefficients recovery in the area of numerical and computational theory. This work is devoted to shedding some light on the benefits of using such framelets in the area of numerical computations of integral equations. We introduce a new numerical method for solving Volterra integral equations. It is based on pseudo-spline quasi-affine tight framelet systems generated via the oblique extension principles. The resulting system is converted into matrix equations via these generators. We present examples of the generated pseudo-splines quasi-affine tight framelet systems. Some numerical results to validate the proposed method are presented to illustrate the efficiency and accuracy of the method.

1. Introduction

Many natural science problems are modeled by Volterra integral equations, which therefore has brought them much attention from scientists in numerical analysis. Yet, many numerical schemes used wavelet representation to numerically solve some integral equations. However, some approximations work better with redundant expansions such as the biorthogonal wavelet (or simply, framelet) expansions. The redundancy property of framelets has been used for many applications in science and engineering disciplines, for example, in the analysis of the Gibbs phenomenon and numerical solutions of various types of integral equations (see, e.g., [1,2,3,4,5,6]), in time–frequency theory for image analysis, multifilter designs in electrical engineering, the theory of nonshift and shift-invariant spaces, and many other areas (see, e.g., [7,8,9,10,11,12,13]. It is known that the approximation accuracy improved via tight framelets due to their redundancy. Note that, in orthonormal expansion analysis, the redundancy is missing. Therefore, we have more freedom in building efficient and accurate recovery.
The aim of this paper is to present a numerical method by using a specific type of framelets generated using the unitary and oblique extension principles for approximating the solution of Volterra integral equations defined by
u ( x ) = f ( x ) + λ a x K ( x , t ) u ( t ) d t .
It is difficult, in most cases, to find the solution of the Volterra integral equations analytically. The collocation-type method is well known as an accurate numerical technique for integral equations.
We use a new and accurate method that generalizes the wavelet-collocation method used in the literature. We call it the framelet-collocation method.
Our paper is organized as follows. In Section 2, we provide some preliminary background on redundant systems (tight frames), their notations, and function expansion. Section 3 provides some principles in the construction of pseudo-spline quasi-affine tight framelet systems using the oblique extension principles. We then start the presentation of matrix assembly for solving Volterra integral equations based on the collocation-type pseudo-spline-quasi-affine-tight-framelets-based method in Section 4. We further test our method on a numerical example and some graphical illustrations in Section 5. In Section 6, we conclude with some comments.

2. Preliminary Results

The expansion of a function is not limited, in general, to a specific form, and we can have a redundancy for a given representation, for example, in the expansion generated via tight frames. The idea of the frame sequence was introduced in Ref. [14], where frames were used in the mathematical construction in the analysis of non-harmonic Fourier expansions. In Ref. [15], Daubechies presented these sets of tight frames of L 2 ( R ) in some applications of signal analysis.
Definition 1.
A sequence { f k } k = 1 of generators in the space L 2 ( R ) is called a framelet for L 2 ( R ) if numbers A , B > 0 such that
A f 2 k = 1 f , f k 2 B f 2 , f L 2 ( R ) .
The constants A , B are called framelet bounds. A framelet is called tight if it is possible to have A = B as a framelet bound. In fact, framelets are extensions of orthonormal bases. The space L 2 ( R ) is the set elements g ( x ) such that
g L 2 ( R ) = R | g ( x ) | 2 1 / 2 < .
Let f L 2 ( R ) , then the dilation and translation functions, D and T, are defined by D f ( x ) = 2 f ( 2 x ) and T a f ( x ) = f ( x a ) for a R , respectively. Note that for j Z , we have T a D j = D j T 2 j a and D j T a = T 2 j a D j . Define
2 ( Z ) = h [ k ] : k = | h [ k ] | 2 < , k Z .
For f L 2 ( R ) , the Fourier transform and its inverse are defined, respectively, by
f ^ ( ξ ) = R f ( t ) e i ξ t d t , ξ R
and
f ( x ) = 1 2 π R f ^ ( ξ ) e i ξ x d ξ , x R .
Definition 2.
Let ϕ be a compactly supported function in L 2 ( R ) . Then ϕ is a refinable function if there exists a sequence h 0 [ k ] 2 ( Z ) such that the following equation is satisfied:
ϕ ( x ) = 2 k Z h 0 [ k ] ϕ ( 2 x k ) .
Here, we call h 0 the low mask filter of the refinable function ϕ .
Note that we can write Equation (1) in the frequency domain as
ϕ ^ = ( h 0 ^ ϕ ^ ) · / 2 ,
for some 2 π -periodic h 0 ^ . In addition, it known that if ϕ L 2 ( R ) such that ϕ ^ ( 0 ) = 1 , then we have, h ^ 0 ( 0 ) = 1 and h ^ 0 ( π ) = 0 , (see, for example, Refs. [7,16,17]). Hence, h ^ 0 ( ξ ) can be written as
h ^ 0 ( ξ ) = 1 + e i ξ 2 n τ ( ξ ) ,
where n N refers to the highest multiplicity of the roots of h ^ 0 ( π ) such that τ ( ξ ) is a polynomial of trigonometric functions with τ ( 0 ) = 1 . One can easily show that Equation (2) can be rewritten as
ϕ ^ ( ξ ) = j = 1 h ^ 0 ( 2 j ξ ) .
Hence, by Equation (3), we have
ϕ ^ ( ξ ) = j = 1 1 + e i ( 2 j ξ ) 2 n j = 1 τ ( 2 j ξ ) = 1 e i ξ i ξ n j = 1 τ ( 2 j ξ ) .
We say that a function ψ has a vanishing moment of order m if
x p ψ ( x ) d x = 0 ; p = 0 , , m 1 .
It is known that multiresolution analysis (MRA) is a tool to generate wavelet bases. To formulate the matrix from and the numerical solution of a given Volterra integral equation, we will study and use pseudo-spline quasi-affine tight framelets and their constructions that are derived from an MRA, and in particular, the oblique extension principle (OEP) in Ref. [17]. The interested reader should consult [18,19,20] and other related references for more details.
Define Ψ = ψ = 1 r L 2 ( R ) as
ψ = 2 k Z h [ k ] ϕ ( 2 · k ) .
The sequence h [ k ] , k Z = 1 r is called the high mask filter of Ψ . Equation (5) can be expressed in terms of its Fourier representation and is given by
ψ ^ = ( h ^ ϕ ^ ) ( · / 2 ) , = 1 , , r ,
where h ^ ( · ) = h ^ ( · + 2 π ) , = 1 , , r .
Theorem 1.
Assume that ϕ is a refinable function in L 2 ( R ) with compact support. Let h 0 be its finitely supported low mask filter. Let
h [ k ] , k Z = 1 r
be sequences with finite support. Then,
X Ψ = ψ , j , k : = 1 , , r ; j , k are integers
generates a tight framelet system for L 2 ( R ) if the following equations are satisfied such that π ξ π , where
= 0 r | h ^ ( ξ ) | 2 = 1 a n d = 0 r h ^ ( ξ ) h ^ ( ξ + π ) = 0 .
Proof. 
See Ref. [17]. □
By Theorem 1, it can be concluded that for any function f L 2 ( R ) , we have the tight framelet representation given by
f = = 1 r j Z k Z f , ψ , j , k ψ , j , k .
The expansion in Equation (8) is known as the best possible expansion of the function f, where it can be truncated by S n , such that
S n f = = 1 r j < n k Z f , ψ , j , k ψ , j , k .
We will use Equation (9) to find the numerical solution of a given Volterra integral equation using quasi-affine tight framelets generated by pseudo-spline functions.

3. Pseudo-Spline Quasi-Affine Tight Framelets

B-spline tight framelets are one of the most important framelets in the framelet family. They are interesting due to their simple structure and properties, e.g., they have compact support and are given by explicit and quite simple formulas in the time and Fourier domain. The smoothness of the B-spline increases as we increase n. It has an important role in applied numerical mathematics, geometric analysis, and many other areas (see, e.g., Refs. [21,22]).
Definition 3.
The B-spline N m + 1 is defined by
N m + 1 x = 0 1 N m x t d t ,
where N 1 ( x ) is the indicator function on the interval [ 0 , 1 ) .
Definition 4.
Let m N . Then we define the B-spline B m by the following equation:
B m x : = N m ( x + m 2 ) .
Hence, we define B m by
B m + 1 : = B m B 1 , m N ,
where B 1 ( x ) = χ 1 / 2 , 1 / 2 x .
For m = 1 , , 4 , we plot the graphs of B-splines in Figure 1.
One can easily show that the Fourier transform of the B-spline, B m , of order m is given by
B ^ m ( ξ ) = sin ( ξ / 2 ) ξ / 2 m a n d h ^ 0 ( ξ ) = cos m ( ξ / 2 ) .
Note that
N ^ m ( ξ ) = B ^ m ( ξ m 2 ) = e i m ξ 2 B ^ m ( ξ ) , and its low mask filter is e i m ξ 2 h ^ 0 ( ξ ) .
We refer the reader to [23] for more details.
The unitary extension principle (UEP) is known as a method to generate tight framelets via a refinable function. It is known that the MRA is a special case of the well known UEP. In addition, the UEP was extended to the OEP in [18,20] by finding a 2 π -periodic function Θ . For any refinable function and to construct a tight framelet system, the non-negative function Θ , which is essentially bounded and continuous at zero such that Θ ( 0 ) = 1 , shall satisfy the following conditions:
Θ ( 2 ξ ) | h ^ 0 ( ξ ) | 2 = Θ ( ξ ) = 1 r | h ^ ( ξ ) | 2 ; Θ ( 2 ξ ) h ^ 0 ( ξ ) h ^ 0 ( ξ + π ) + = 1 r h ^ ( ξ ) h ^ ( ξ + π ) = 0 .
Definition 5
([24]). Suppose that the conditions of the UEP hold for Ψ. Then, the quasi-affine system X J Ψ generated using Ψ is defined by
X J ( Ψ ) = ψ , j , k : = 1 , , r ; j , k Z
such that
ψ , j , k = 2 j / 2 ψ ( 2 j · k ) , j J 2 j ψ ( 2 j · k ) , j < J .
Here, for our proposed method, we consider the system above for the case where J = 0 .
If h ^ 0 is the low mask filter of a given refinable function ϕ , then using the OEP, it is assumed [24] that
Θ ( ξ ) Θ ( 2 ξ ) | h ^ 0 ( ξ ) | 2 Θ ( 2 ξ ) | h ^ 0 ( ξ + π ) | 2 .
This condition helps to find the high mask filters of the required framelet system. Let | h | 2 = H , where
H = Θ ( ξ ) Θ ( 2 ξ ) | h ^ 0 ( ξ ) | 2 Θ ( 2 ξ ) | h ^ 0 ( ξ + π ) | 2
and
| θ | 2 = Θ .
Here, the square root is obtained by the spectral factorization in Ref. [7]. Assume that c 2 , c 3 are two 2 π -periodic trigonometric functions/polynomials such that
| c 2 ( ξ ) | 2 + | c 3 ( ξ ) | 2 = 1 , c 2 ( ξ ) c 2 ( ξ + π ) ¯ + c 3 ( ξ ) c 3 ( ξ + π ) ¯ = 0 .
Then, we can find three high mask filters, namely,
h ^ 1 ( ξ ) = e i ξ θ ( 2 ξ ) h ^ 0 ( ξ + π ) ¯ , h ^ 2 ( ξ ) = c 2 ( ξ ) h ( ξ ) , 4 h ^ 3 ( ξ ) = c 3 ( ξ ) h ( ξ ) ,
with a standard choice of c 2 ( ξ ) = ( 1 / 2 ) and c 3 ( ξ ) = ( 1 / 2 ) e i ξ . If we consider the UEP rather than the OEP in the construction above, i.e., Θ = 1 , then we will use the assumption that
| h ^ 0 ( ξ ) | 2 1 | h ^ 0 ( ξ + π ) | 2 .
Define the high mask filters by
h ^ 1 ( ξ ) = e i ξ h ^ 0 ( π + ξ ) ¯ , h ^ 2 ( ξ ) = ( 2 / 2 ) h ( ξ ) , h ^ 3 ( ξ ) = e i ξ h ^ 2 ( ξ ) .
Note that we can reduce the number of framelets from three to two with the new fundamental function 1 H , where
h ^ 1 ( ξ ) = e i ξ θ ( 2 ξ ) h ^ 0 ( π + ξ ) ¯ , h ^ 2 ( ξ ) = h 0 ( ξ ) h ( 2 ξ ) .
However, this will usually affect the framelet system by having less symmetry of the framelets or longer filters.
Let ϕ be as Equation (1), which generates an MRA ( V j ) j , and W 2 m ( R ) be the Sobolev space. Then, X ( Ψ ) provides approximation order m if
f S n f 2 = O ( 2 n m ) , f W 2 m ( R ) .
The approximation order of the truncated function S n was studied in [18,25]. It is known that the approximation orders rely on the behavior near zero of the function Λ ϕ , where
Λ ϕ = 1 | ϕ ^ | 2 [ ϕ ^ , ϕ ^ ]
and
[ f , g ] ( ξ ) = k 2 π Z f ( ξ + k ) g ( ξ + k ) ¯ .
Note that the refinable function ϕ satisfies m-Strang–Fix condition if the following equation is satisfied, where
ϕ ^ ( 0 ) 0 , ϕ ^ ( j ) ( 2 π k ) = 0 , j 0 , , m 1 , k is   a   non - zero   integer .
It was proved in Ref. [25] that if the function ϕ provides approximation order m, then Λ ϕ has a zero of order m at the origin. This means that Λ ϕ has a zero of order m at ξ = 0 , and then ϕ ^ has a zero of order m. In addition, Jetter et al. in [26] showed that depending on the OEP construction, the truncated S n provides approximation order m iff [ ϕ ^ , ϕ ^ ] | ϕ ^ | 2 = O ( 2 2 m ) .
Daubechies, in Ref. [18], has proved that if the system X ( Ψ ) has a vanishing moment of order m 1 and ϕ has an approximation order m, then the approximation order of X ( Ψ ) is equal to the minimum of m and 2 m 1 . To have high approximation orders, we have to construct refinable functions where the Fourier transforms are very smooth at the origin. This leads to the well-known refinable functions, pseudo-splines, and their tight framelet generators.
Pseudo-splines provide us with a nice class of refinable and compactly supported functions. The first type was introduced in [18,27] to construct a special type of tight framelets and type II were introduced in [28] to construct tight framelets with specific properties of symmetry. In the frequency domain and for non-negative integers l , m such that l < m , pseudo-splines of type I (or PS-I- ( m , l ) and type II with order m and type l (or PS-II- ( m , l ) ) can be defined by
k ϕ ^ ( ξ ) = i = 1 k h ^ 0 ( 2 i ξ ) with k ϕ ^ ( 0 ) = 1 , for k = 1 , 2 ,
where the low mask filter of the pseudo-splines of type I with order ( m , l ) is defined by
| 1 h ^ 0 ( ξ ) | 2 = cos 2 m ( ξ / 2 ) k = 0 l m + l k sin 2 k ( ξ / 2 ) cos 2 ( l k ) ( ξ / 2 )
and the low mask filter of the pseudo-splines of type II with order ( m , l ) is defined by
2 h ^ 0 ( ξ ) = cos 2 m ( ξ / 2 ) k = 0 l m + l k sin 2 k ( ξ / 2 ) cos 2 ( l k ) ( ξ / 2 ) .
Note that if r = 0 , pseudo-splines of both types are B-splines. It is known that the smoothness of the pseudo-spline increases with m and decreases with l (see Ref. [28]). According to spectral factorization, or by using the Fejér–Riesz lemma (see [7]), the low mask filter of the pseudo-spline of type I is obtained by taking the square root of type II, i.e., 2 h ^ 0 ( ξ ) = | 1 h ^ 0 ( ξ ) | 2 . In general, we have the following lemma.
Lemma 1.
Assume that L ( ξ ) is a positive valued trigonometric polynomial given by
f ( ξ ) = m = 0 M a m cos ( m ξ ) , w i t h a m R .
Then a trigonometric polynomial g of order M, where
g ( ξ ) = m = 0 M b m e i m ξ , w i t h b m R ,
such that | g ( ξ ) | 2 = f ( ξ ) .
Proof. 
See Ref. [7]. □
In Mallat’s construction (see Ref. [29]), it is shown that h ^ 0 ( ξ ) h ^ 0 ( ξ + π ) ¯ + h ^ 1 ( ξ ) h ^ 1 ( ξ + π ) ¯ = 0 , where h ^ 1 ( ξ ) = e i ξ h ^ 0 ( ξ + π ) ¯ , and that
H ( ξ ) = cos 2 m ( ξ ) k = + 1 m 1 m + k cos 2 ( k ) ( ξ / 2 ) sin 2 k ( ξ / 2 ) .
If we take m = 4 and = 1 , then we will get short filters compared with the general case. This is because of the form of H ( ξ ) , where
H ( ξ ) = k = 2 3 5 k cos 10 2 k ( ξ / 2 ) sin 2 k ( ξ / 2 ) = 10 cos 4 ( ξ / 2 ) sin 4 ( ξ / 2 ) .
In fact, we have the following fact.
Proposition 1.
For non-negative integers l , m such that l < m , If = m 3 , then
H ( ξ ) = 2 m 3 m 1 cos 2 m 4 ( ξ / 2 ) sin 2 m 4 ( ξ / 2 ) .
Proof. 
As 2 m 3 m 2 = 2 m 3 m 1 , then we have
H ( ξ ) = k = m 2 m 1 2 m 3 k cos 4 m 6 2 k ( ξ / 2 ) sin 2 k ( ξ / 2 ) .
Therefore,
H ( ξ ) = 2 m 3 m 2 cos 2 m 4 ( ξ / 2 ) sin 2 m 2 ( ξ / 2 ) + 2 m 3 m 1 cos 2 m 2 ( ξ / 2 ) sin 2 m 4 ( ξ / 2 ) = 2 m 3 m 1 cos 2 m 4 ( ξ / 2 ) sin 2 m 4 ( ξ / 2 ) .
 □
Using pseudo-splines of both types, we give some examples of quasi-affine tight framelets constructed via the OEP.
Example 1 (PS-I -(4,1)).
For the order ( 4 , 1 ) , consider the pseudo-spline of type I, 1 ϕ ^ ( ξ ) . Then, its low mask filter can be given by | 1 h ^ 0 ( ξ ) | 2 = cos 8 ( ξ / 2 ) ( 1 + 4 sin 2 ( ξ / 2 ) ) . Note that by using Lemma 1, we have
1 h ^ 0 ( ξ ) = e i 5 ξ ( 1 + e i ξ ) 4 ( 2 + ( 3 + 5 ) e i ξ ) 16 ( 1 5 ) .
Define
1 h ^ 1 ( ξ ) = e i ξ 1 h ^ 0 ( ξ + π ) , 1 h ^ 2 ( ξ ) = 5 2 sin 2 ( ξ ) , and 1 h ^ 3 ( ξ ) = e i ξ 1 h ^ 2 ( ξ ) .
Let Ψ = 1 ψ ; = 1 , 2 , 3 , where 1 ψ ^ = 1 h ^ ( ξ / 2 ) 1 ϕ ^ ( ξ / 2 ) ; = 1 , 2 , 3 . Then the system X 0 ( Ψ ) forms a quasi-affine tight framelet for L 2 ( R ) .
The pseudo-spline of type I, with its corresponding quasi-affine tight framelets generated by using 1 ϕ ^ ( ξ ) of order ( 4 , 1 ) , is depicted in Figure 2.
Example 2 (PS-II -(3,1)).
For the order ( 3 , 1 ) , consider the pseudo-spline of type II, 2 ϕ ^ ( ξ ) . Its low mask filter is given by
2 h ^ 0 ( ξ ) = cos 6 ( ξ / 2 ) ( 1 + 3 sin 2 ( ξ / 2 ) ) .
Define
2 h ^ 1 ( ξ ) = e i ξ 2 h ^ 0 ( ξ + π ) , 2 h ^ 2 ( ξ ) = h ( ξ ) + e i ξ h ( ξ )
and
2 h ^ 3 ( ξ ) = e i ξ h ( ξ ) h ( ξ ) ,
where
H ( ξ ) = 0.30108642578125305 0.2014160156250013 e 2 i ξ 0.20141601562500128 e 2 i ξ + 0.05090332031249944 e 4 i ξ + 0.050903320312499445 e 4 i ξ + 0.00024414062500026867 e 6 i ξ + 0.00024414062500026867 e 6 i ξ 0.00027465820312493067 e 8 i ξ 0.00027465820312493067 e 8 i ξ and h ( ξ ) = 0.22813823298962 + 0.00139868605052 e 2 i ξ + 0.44712319189971 e 2 i ξ + 0.00123930398199 e 4 i ξ 0.2216229489426 e 4 i ξ .
The pseudo-spline of type II with its corresponding quasi-affine tight framelets generated by 2 ϕ ^ ( ξ ) of order ( 3 , 1 ) is depicted in Figure 3.

4. Matrix Assembly via Pseudo-Spline Quasi-Affine Tight Framelets

Consider the Volterra integral equation defined by
u ( x ) = f ( x ) + λ a x K ( x , t ) u ( t ) d t ,
where λ R , f, and K are given and known functions and u is an unknown function to be approximated. K is called the kernel of Equation (11). A function u ( x ) defined over [ a , b ] can be expressed by framelets as Equation (8). To find an approximate solution u n of Equation (11), we will truncate the framelet representation of u as Equation (9). Then,
u ( x ) u n ( x ) = = 1 r j < n k Z c j , k ψ , j , k ,
where c j , k = u n , ψ , j , k . Substituting Equation (12) into Equation (11) and by using the suitable collocation points to the truncated expansion, we have
= 1 r j < n i , k Z c j , k ψ , j , k ( x i ) = f ( x ) + λ = 1 r j < n i , k Z a x K ( x i , t ) c j , k ψ , j , k ( t ) d t .
Equation (13) can be simplified to a system of equations with the unknown coefficients c j , k given by
= 1 r j < n k Z c j , k M , j , k ( x i ) = f , j , k ( x i ) ,
where
M , j , k ( x i ) = ψ , j , k ( x i ) λ a b K ( x i , t ) ψ , j , k ( t ) d t .
Now the unknown coefficients are determined by solving the resulting system of equations obtained from Equation (14), and then we get the approximate solution u n . The absolute error for this formulation is defined by
e n = u n ( x ) u ( x ) 2 , x [ a , b ] .

5. Numerical Performance

To validate the accuracy of our method, in this section we present the following example of Volterra integral equations. The numerical results obtained here using Mathematica software.
Example 3.
We consider the following Volterra–Fredholm integral equation:
u ( x ) = f ( x ) + 0 x e t cos ( x ) u ( t ) d t ,
where
f ( x ) = e x 1 2 cos ( x ) ( e 2 x 1 )
and the exact solution is u ( x ) = e x .
In Table 1 and Table 2, the error e n for different values of n and the numerical values of the exact and approximated solution u n ( x ) when n = 2 are computed, respectively. Some illustrations for the graphs of the exact and approximate solutions and the error are depicted in Figure 4 and Figure 5.
To validate the proposed method, we provide Figure 6 to show the rate of convergence of Example 3 in the log–log scale plot by using both systems of pseudo-spline quasi-affine tight framelets PS-I -(4,1) and PS-II -(3,1) generated using the OEP.

6. Conclusions

A collocation-type pseudo-spline-quasi-affine-framelet-based method is developed to numerically solve a given Volterra integral equation. This is an important research direction in the filed of framelet-based numerical schemes for integral equations. With a few orders of truncated partial sums, the results show that the proposed method is effective and accurate. It turns out that increasing the order of the truncated partial sums of the framelet system and its vanishing moments will dramatically increase the approximation orders as well as the accuracy of the solution. In addition, the accuracy orders of the approximated solution using both pseudo-spline quasi-affine tight framelet systems were close to each other, with slight differences and preferences to those with a higher order of m.

Acknowledgments

This work was supported by the Zayed University Research Fund.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Mohammad, M.; Lin, E.B. Gibbs phenomenon in tight framelet expansions. Commun. Nonlinear Sci. Numer. Simul. 2018, 55, 84–92. [Google Scholar] [CrossRef]
  2. Mohammad, M.; Lin, E.B. Gibbs Effects Using Daubechies and Coiflet Tight Framelet Systems. Contemp. Math. 2018, 706, 271–282. [Google Scholar]
  3. Mohammad, M. Special B-spline Tight Framelet and It’s Applications. J. Adv. Math. Comput. Sci. 2018, 29, 1–18. [Google Scholar] [CrossRef]
  4. Mohammad, M. On the Gibbs Effect Based on the Quasi-Affine Dual Tight Framelets System Generated Using the Mixed Oblique Extension Principle. Mathematics 2019, 7, 952. [Google Scholar] [CrossRef]
  5. Mohammad, M.; Howari, F.; Acbas, G.; Nazzal, Y.; AlAydaroos, F. Wavelets Based Simulation and Visualization Approach for Unmixing of Hyperspectral Data. Int. J. Earth. Environ. Sci. 2018, 3, 152. [Google Scholar] [CrossRef] [PubMed]
  6. Mohammad, M. A Numerical Solution of Fredholm Integral Equations of the Second Kind Based on Tight Framelets Generated by the Oblique Extension Principle. Symmetry 2019, 11, 854. [Google Scholar] [CrossRef]
  7. Daubechies, I. Ten Lectures on Wavelets; SIAM: Philadelphia, PA, USA, 1992. [Google Scholar]
  8. Aldwairi, M.; Guled, M.; Cassada, M.; Pratt, M.; Stevenson, D.; Franzon, P. Switch Architecture for Optical Burst Switching Networks. Available online: https://www.researchgate.net/profile/Monther_Aldwairi/publication/229004330_Switch_architecture_for_optical_burst_switching_networks/links/0deec51b1c89fa102d000000.pdf (accessed on 1 November 2019).
  9. Al-Ayyoub, M.; Jararweh, Y.; Rabab’ah, A.; Aldwairi, M. Feature extraction and selection for Arabic tweets authorship authentication. J. Ambient. Intell. Humaniz. Comput. 2017, 8, 383–393. [Google Scholar] [CrossRef]
  10. Gröchenig, K. Foundations of Time-Frequency Analysis; Birkhäuser Boston: Cambridge, MA, USA, 2001. [Google Scholar]
  11. Mallat, S. A Wavelet Tour of Signal Processing; Elsevier: Amsterdam, The Netherlands, 1999. [Google Scholar]
  12. Meyer, Y. Wavelets and Operators; Cambridge University Press: Cambridge, UK, 1992. [Google Scholar]
  13. Meyer, Y. Oscillating Patterns in Image Processing and Nonlinear Equations; American Mathematical Society: Providence, RI, USA, 2001. [Google Scholar]
  14. Duffin, R.J.; Schaeffer, A.C. A class of nonharmonic Fourier series. Trans. Amer. Math. Soc. 1952, 72, 341–366. [Google Scholar] [CrossRef]
  15. Daubechies, I.; Grossmann, A.; Meyer, Y. Painless nonorthogonal expansions. J. Math. Phys. 1986, 27, 1271–1283. [Google Scholar] [CrossRef]
  16. Jiang, Q.; Shen, Z. On existence and weak stability of matrix refinable functions. Constr. Approx. 1999, 15, 337–353. [Google Scholar] [CrossRef]
  17. Ron, A. Factorization theorems of univariate splines on regular grids. Israel J. Math. 1990, 70, 48–68. [Google Scholar] [CrossRef]
  18. Daubechies, I.; Han, B.; Ron, A.; Shen, Z. Framelets: MRA-based constructions of wavelet frames. Appl. Comput. Harmon. Anal. 2003, 14, 1–46. [Google Scholar] [CrossRef] [Green Version]
  19. Ron, A.; Shen, Z. Affine systems in L 2 ( R d ) II: Dual systems. J. Fourier Anal. Appl. 1997, 3, 617–637. [Google Scholar] [CrossRef]
  20. Chui, C.K. An Introduction to Wavelets; Elsevier: Amsterdam, The Netherlands, 2016. [Google Scholar]
  21. Chui, C.K.; He, W.; Stöckler, J. Compactly supported tight and sibling frames with maximum vanishing moments. Appl. Comput. Harmon. Anal. 2002, 13, 224–262. [Google Scholar] [CrossRef] [Green Version]
  22. Christensen, O. An Introduction to Frames and Riesz Bases; Birkhäuser Boston: Cambridge, MA, USA, 2003. [Google Scholar]
  23. De Boor, C. A Practical Guide to Splines; Springer: New York, NY, USA, 1978. [Google Scholar]
  24. Ron, A.; Shen, Z. Affine systems in L 2 ( R d ) : The analysis of the analysis operators. J. Funct. Anal. 1997, 148, 408–447. [Google Scholar] [CrossRef]
  25. De Boor, C.; DeVore, R.; Ron, A. Approximation from shift-invariant subspaces of L 2 ( R d ) . Trans. Amer. Math. Soc. 1994, 341, 787–806. [Google Scholar]
  26. Jetter, K.; Zhou, D.X. Order of linear approximation from shift invariant spaces. Constr. Approx. 1995, 11, 423–438. [Google Scholar] [CrossRef]
  27. Selesnick, I. Smooth wavelet tight frames with zero moments. Appl. Comput. Harmon. Anal. 2001, 10, 163–181. [Google Scholar] [CrossRef]
  28. Dong, B.; Shen, Z. Pseudo-splines, wavelets and framelets. Appl. Comput. Harmon. Anal. 2007, 22, 78–104. [Google Scholar] [CrossRef] [Green Version]
  29. Mallat, S.G. Multiresolution approximations and wavelet orthonormal bases of L 2 ( R ) . Trans. Amer. Math. Soc. 1989, 315, 69–87. [Google Scholar]
Figure 1. The B-splines B m for m = 1 , , 4 , respectively.
Figure 1. The B-splines B m for m = 1 , , 4 , respectively.
Entropy 21 01098 g001
Figure 2. The type I pseudo-spline of order ( 4 , 1 ) along with its quasi-affine tight framelets, respectively.
Figure 2. The type I pseudo-spline of order ( 4 , 1 ) along with its quasi-affine tight framelets, respectively.
Entropy 21 01098 g002
Figure 3. The pseudo-spline scaling function of type II with order ( 3 , 1 ) along with its quasi-affine tight framelets, respectively.
Figure 3. The pseudo-spline scaling function of type II with order ( 3 , 1 ) along with its quasi-affine tight framelets, respectively.
Entropy 21 01098 g003
Figure 4. The graphs of u and u n for n = 2 , 3 , respectively, of Example 3 based on the PS-I -(4,1) framelet system.
Figure 4. The graphs of u and u n for n = 2 , 3 , respectively, of Example 3 based on the PS-I -(4,1) framelet system.
Entropy 21 01098 g004
Figure 5. The graphs of the error for n = 2 , 3 , respectively, of Example 3 based on the PS-I -(4,1) framelet system.
Figure 5. The graphs of the error for n = 2 , 3 , respectively, of Example 3 based on the PS-I -(4,1) framelet system.
Entropy 21 01098 g005
Figure 6. The rate of convergence of the proposed method for Example 3.
Figure 6. The rate of convergence of the proposed method for Example 3.
Entropy 21 01098 g006
Table 1. The error e n of Example 3 using different types of pseudo-spline framelet systems.
Table 1. The error e n of Example 3 using different types of pseudo-spline framelet systems.
nPS-I -(4,1)PS-II -(3,1)
2 1.08 × 10 3 6.438 × 10 4
3 2.63 × 10 4 1.67 × 10 4
4 4.32 × 10 5 3.21 × 10 5
5 2.83 × 10 6 0.57 × 10 5
6 5.58 × 10 7 6.93 × 10 6
7 1.01 × 10 8 2.44 × 10 8
8 6.44 × 10 9 1.56 × 10 9
9 9.35 × 10 11 3.45 × 10 10
10 2.35 × 10 12 6.69 × 10 11
Table 2. Numerical results of the function u n of Example 3 using different types of pseudo-spline framelet systems and for a level n = 2 .
Table 2. Numerical results of the function u n of Example 3 using different types of pseudo-spline framelet systems and for a level n = 2 .
nExactPS-I -(4,1)PS-II -(3,1)
0.1 1.10517 1.105117381 1.105117134
0.2 1.22140 1.222175185 1.222124175
0.3 1.34986 1.350694516 1.340694516
0.4 1.49182 1.491882666 1.491802512
0.5 1.64872 1.646677909 1.646672532
0.6 1.82212 1.822082188 1.822194893
0.7 2.01375 2.015056716 2.015321421
0.8 2.22554 2.226976714 2.226776701
0.9 2.45960 2.459772459 2.459610020
1.0 2.71828 2.715929278 2.715989072

Share and Cite

MDPI and ACS Style

Mohammad, M. Biorthogonal-Wavelet-Based Method for Numerical Solution of Volterra Integral Equations. Entropy 2019, 21, 1098. https://0-doi-org.brum.beds.ac.uk/10.3390/e21111098

AMA Style

Mohammad M. Biorthogonal-Wavelet-Based Method for Numerical Solution of Volterra Integral Equations. Entropy. 2019; 21(11):1098. https://0-doi-org.brum.beds.ac.uk/10.3390/e21111098

Chicago/Turabian Style

Mohammad, Mutaz. 2019. "Biorthogonal-Wavelet-Based Method for Numerical Solution of Volterra Integral Equations" Entropy 21, no. 11: 1098. https://0-doi-org.brum.beds.ac.uk/10.3390/e21111098

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop