Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2008, Article ID 629253, 14 pages
http://dx.doi.org/10.1155/2008/629253
Research Article

Construction of Interval Wavelet Based on Restricted Variational Principle and Its Application for Solving Differential Equations

College of Information and Electrical Engineering, China Agricultural University, P.O. Box 053, 17 Qinghua Donglu Road, Beijing 100083, China

Received 10 December 2007; Revised 25 February 2008; Accepted 31 March 2008

Academic Editor: Giuseppe Rega

Copyright © 2008 Shu-Li Mei et al. 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

Based on restricted variational principle, a novel method for interval wavelet construction is proposed. For the excellent local property of quasi-Shannon wavelet, its interval wavelet is constructed, and then applied to solve ordinary differential equations. Parameter choices for the interval wavelet method are discussed and its numerical performance is demonstrated.

1. Introduction

Since the definition domain of wavelet transformation is an infinite interval, the boundary effect would occur when applied for resolving the engineering problems with bounded interval, for example, ordinary differential equations (ODEs). Consequently, it will decrease the precision and computational efficiency of the solution. Nevertheless, the boundary effect can be eliminated effectively by constructing an interval wavelet using numerical methods.

There are several ways available to construct an interval wavelet. In general, the construction method is relative to the wavelet function, that is, different interval wavelet has different construction method. A simple solution is the even 2-periodical extension 𝑓 of function 𝑓[0,1], which is usually used in image analysis. Unfortunately, this extension generally produces discontinuities at the integers that are indicated by the large wavelet coefficients near the endpoints 0 and 1. Thus, the constructed wavelet cannot exactly analyze the boundary behavior of a given function. To solve this problem, the popular method is based on special boundary and interior scaling functions as well as wavelets to reduce the numerical problem at the boundaries [1, 2].

The aim of this paper is to introduce a general construction method of interval wavelet based on the restricted variational principle. As an example, the quasi-Shannon wavelet is firstly introduced in brief. Then, its corresponding interval wavelet is constructed in detail based on the restricted variational principle [3]. After that, the wavelet collocation method is applied to obtain the interval wavelet discrete formulation of ODEs. The performance of the quasi-Shannon interval wavelet is illustrated by comparing the numerical results of Shannon, quasi-Shannon, and Shannon interval wavelets on a convection equation.

2. Quasi-Shannon Wavelet

The Shannon scaling function is smooth, this means that the function and all of its derivatives exist and are continuous. But the Shannon function does not have compact support. In order to improve the localized and asymptotic behavior of the Shannon scaling function, Wei [4] introduced a regularization factor 𝑅(𝑥) as follows:𝑤(𝑥)=𝜙(𝑥)𝑅(𝑥),𝜙(𝑥)=sin(𝜋𝑥)𝜋𝑥.(2.1) The advantage of the regularized Shannon scaling function is that its Fourier transformation is continuous [4], resulting in an excellent local property. A common and important regularization factor 𝑅(𝑥) proposed by Wei [4] is the Gaussian function𝑅𝜎𝑥(𝑥)=exp22𝜎2,𝜎>0,(2.2) where 𝜎 is the width parameter (or called window size). In practice, the best results are usually obtained, when 𝜎=𝑟Δ [4] (𝑟 is a parameter chosen in computations, Δ is the size of cell in discrete mesh). Substituting (2.2) into (2.1), the Gaussian regularized orthogonal sampling scaling function can be obtained as following:𝑤(𝑥)=sin(𝜋𝑥)𝑥𝜋𝑥exp22𝜎2,𝜎>0.(2.3) The corresponding discrete formula is𝑤𝑥𝑥𝑖=sin𝜋/Δ𝑥𝑥𝑖𝜋/Δ𝑥𝑥𝑖exp𝑥𝑥𝑛22𝜎2.(2.4)

3. Construction of Interval Wavelet Based on Restricted Variational Principle

3.1. Restricted Variational Principle

Defining a scalar fonctionelle of 𝑢 as=Ω𝐹𝑢,𝜕𝑢𝜕𝑥,dΩ+Γ𝐸𝑢,𝜕𝑢𝜕𝑥,dΓ,(3.1) where 𝑢 is an unknown function, both 𝐹 and 𝐸 are specific operators, Ω is the definition domain, and Γ is the boundary of Ω. The solution of a continuous medium problem is a stagnation point of the functional equation , that is, its variation is equal to zero𝛿=0.(3.2) The method described above for solving a continuous medium problem is called restricted variational principle [3].

The unknown function 𝑢 could be approximately expressed as𝑢̃𝑢=𝑛𝑖=1𝑁𝑖𝑎𝑖=𝑁𝑎,(3.3) where 𝑎𝑖 is the parameters to be determined and 𝑁𝑖 is the trial function.

If the unknown function 𝑢 is restricted by additional condition 𝐶(𝑢)=0 in definition domain Ω, a revised fonctionelle could be constructed as follows:=+Ω𝜆𝑇𝐶(𝑢)dΩ,(3.4) where 𝜆 is called Lagrange multiplier, which is a function vector with independent coordinates in Ω. The functions 𝑢 and 𝜆 are two unknown variables within the revised fonctionelle. Similarly as (3.3), 𝜆 could be approximately expressed by trial function 𝑁𝑖 as̃𝑁𝜆=𝑖𝑏𝑖=𝐍𝐛.(3.5) Let the variation of the revised fonctionelle equal to zero, and a system of equations could be obtained as follows:𝛿𝛿𝛿𝑐={𝛿𝛿𝑎𝑎𝑏𝛿𝑏}=0,𝑐={}.(3.6) The values of 𝑎 and 𝑏 could be obtained by solving the above system of equations.

Let the Euler equation of the fonctionelle is𝐴(𝑢)=0,(3.7) and the additional condition is𝐶(𝑢)=𝐿1(𝑢)+𝐶1=0.(3.8) Substituting (3.5), (3.7), and (3.8) into (3.4), the following function could be obtained:𝛿=𝛿𝑎𝑇Ω𝑁𝑇𝐴(̃𝑢)dΩ+𝛿𝑏𝑇Ω𝐍𝐿1̃𝑢+𝐶1dΩ+𝛿𝑎𝑇Ω𝐿𝑇1̃(𝑁)𝜆dΩ=0.(3.9) Since (3.9) is satisfied for all variations of 𝛿𝑎 and 𝛿𝑏, the following two equations can be derived:Ω𝑁𝑇𝐴(̃𝑢)dΩ+Ω𝐿𝑇1̃(𝑁)𝜆dΩ=0,(3.10)Ω𝐍𝑇𝐿1(̃𝑢)+𝐶1dΩ=0.(3.11)𝐾𝑎=𝑃.(3.12) The first item in (3.10) is the approximate equation of the natural variation of the linear system of equations:𝐾0𝐺𝐶+𝑅=𝐾𝐺𝑇0{𝑎𝑏𝑃𝑄}{}=0,(3.13) Thus, (3.10) and (3.11) could be rewritten in the following form: 𝐾=Ω𝑁𝑇𝐺𝐴(𝑁)dΩ,T=Ω𝐍𝑇𝐿1(𝑁)dΩ,𝑄=Ω𝐍𝑇𝐶1dΩ.(3.14) where𝑢(𝑥)

3.2. Construction of Quasi-Shannon Interval Wavelet

Considering a one-dimensional problem Ω, whose domain of definition [𝑚,𝑛] is 2𝑗+1,𝑗𝑍. Assuming that the number of discrete point is 𝑖, the 𝑥th discrete point of variable 𝑥𝑖=𝑚+𝑛𝑚2𝑗𝑖,𝑖𝑍.(3.15) could be written as𝑁𝑗 If the quasi-Shannon wavelet with interpolation property is employed as the trial function 𝑁𝑗=𝑤𝑗𝑥𝑥𝑖=2sin𝑗𝜋/𝑛𝑚𝑥𝑥𝑖2𝑗𝜋/𝑛𝑚𝑥𝑥𝑖2exp2𝑗1𝑥𝑥𝑖2𝑟2𝑛𝑚2.(3.16), that is,𝐾Any element in the matrices 𝐺𝑇 and 𝐾𝑘,𝑛𝑤=𝐴𝑗𝐺𝑘𝑛,(3.17)T𝑘,𝑛=Ω𝑁T𝑘𝐿1𝑁𝑛dΩ,(3.18) could be expressed, respectively, as𝑁𝑘Γ where 𝑢=𝑢.(3.19) is the Lagrange basic function.

As the definition domain of wavelet transformation is a double infinite interval, the wavelet coefficients are large near the endpoints of the bounded signal, which increases the computational error. To obtain a high computational precision, the additional condition of scalar fonctionelle in boundary 𝐿1𝑁𝑛=𝑖𝑁𝑖𝑢𝑖,(3.20) could be changed to𝑢𝑖 Thus,𝐾0𝐺𝐶+𝑅=𝐾𝐺𝑇0{𝑎𝑏𝑃0}{}=0.(3.21) where 𝐺T𝑘,𝑛=𝑁T𝑘𝑁𝑛.(3.22) is a known value. Recomputing (3.13) with this additional condition, it could be rewritten as follows:𝐍 Correspondingly, (3.18) changed to the form𝑁 The interpolation properties of 𝑎 and 𝑏 imply that 𝑢 and 𝑥𝑖 are the vectors consisting of the approximate values of [𝑚,𝑛] in discrete points 𝑅=2𝑗+1 and extended discrete points outside the domain of [𝑚,𝑛], respectively. Typically, periodic, symmetry, and zero extension methods could be employed. Assuming that there are 𝑥0,𝑥1,,𝑥2𝑗 discrete points in 𝐿, that is, [𝑚,𝑛], and 𝑥𝐿1,𝑥𝐿,,𝑥1 extended points in both outsides of the domain 𝑥𝑅,𝑥𝑅+1,,𝑥𝑅+𝐿, respectively, that is, 𝑢 and 𝑢𝑥1𝑥=𝑢0𝑢𝑥2𝑥=𝑢1𝑢𝑥1𝑖𝑥=𝑢𝑖𝑢𝑥𝐿1𝑥=𝑢𝐿,𝑢𝑥𝑅𝑥=𝑢𝑅1𝑢𝑥𝑅+1𝑥=𝑢𝑅2𝑢𝑥𝑅+𝑖𝑥=𝑢𝑅𝑖1𝑢𝑥𝑅+𝐿𝑥=𝑢𝑅𝐿1(3.23), the function 𝐴{𝑢𝐺1𝐺2𝐿1𝑢𝐿𝑢1𝑢0𝑢1𝑢2𝑗𝑢𝑅𝑢𝑅+1𝑢𝑅+𝐿𝑃}={0},(3.24) can be extended outside of the definition domain using the symmetry method, period method, or the zero extension method. Here, we take the symmetry method as example to illustrate the interval wavelet construction method as follows: 𝐴𝑤𝐴=𝑗𝐴𝑤(00)𝑗𝑤(01)𝐴𝑗02𝑗𝐴𝑤𝑗𝐴𝑤(10)𝑗𝑤(11)𝐴𝑗12𝑗𝐴𝑤𝑗2𝑗𝐴𝑤0𝑗2𝑗𝑤1𝐴𝑗2𝑗2𝑗,𝑎𝐺1=𝐿1,𝐿1𝑤𝑗(0)𝑎𝐿1,𝐿𝑤𝑗(1)𝑎𝐿1,1𝑤𝑗𝑎(𝐿)𝐿,𝐿1𝑤𝑗(1)𝑎𝐿,𝐿𝑤𝑗(0)𝑎𝐿,1𝑤𝑗𝑎(1𝐿)1,𝐿1𝑤𝑗(𝐿)𝑎1,𝐿𝑤𝑗(𝐿1)𝑎1,1𝑤𝑗,𝑏(0)𝐺2=𝑅,𝑅𝑤𝑗(0)𝑏𝑅,𝑅+1𝑤𝑗(1)𝑏𝑅,𝑅+𝐿𝑤𝑗𝑏(𝐿)𝑅+1,𝑅𝑤𝑗(1)𝑏𝑅+1,𝑅+1𝑤𝑗(0)𝑏𝑅+1,𝑅𝑤𝑗𝑏(1𝐿)𝑅+𝐿,𝑅𝑤𝑗(𝐿)𝑏𝑅+𝐿,𝑅+1𝑤𝑗(𝐿1)𝑏𝑅+𝐿,𝑅+𝐿𝑤𝑗.(0)(3.25) Substituting (3.23) into (3.21) results in𝐺1 where𝐺2 In the matrices 𝑎𝑛𝑘 and 𝑏𝑛𝑘, 𝑎𝑛𝑘=𝑙1𝑗𝑘𝑥𝑗𝑛,𝑏𝑛𝑘=𝑙2𝑗𝑘𝑥𝑗𝑛𝑙,(3.26)1𝑗,𝑘=1𝑖=𝐿1𝑖𝑘𝑥𝑥𝑗,𝑖𝑥𝑗,𝑘𝑥𝑗,𝑖,𝑙2𝑗,𝑘=2𝑗+1+𝐿𝑖=2𝑗+1𝑖𝑘𝑥𝑥𝑗,𝑖𝑥𝑗,𝑘𝑥𝑗,𝑖.(3.27) and 𝐺𝑇𝑏=0 could be calculated, respectively, as {𝑢𝐺1+𝐴11𝐴12𝐴13𝐴21𝐴22𝐴23𝐴31𝐴32𝐺2+𝐴330𝑢1𝑢2𝑗𝑃0}={}.(3.28) According to 𝑤𝑗𝑘2=𝑤𝑗+𝑥𝑘1𝑛=𝑁+1𝑎𝑛𝑘𝑤2𝑗𝑤𝑥𝑛,𝑘=0,,𝐿,𝑗𝑘2=𝑤𝑗𝑥𝑘,𝑘=𝐿+1,,2𝑗𝑤𝐿1,𝑗𝑘2=𝑤𝑗+𝑥𝑘2𝑗+𝑁1𝑛=2𝑗+1𝑏𝑛𝑘𝑤2𝑗𝑥𝑛,𝑘=2𝑗𝐿,,2𝑗,(3.29) and the relational expression (3.23), (3.24) could be condensed in the following way:𝐿 Subsequently, as interval interpolation basic functions, the quasi-Shannon interval wavelet can be obtained from (3.28) as follows:𝑁 where sup𝜙=[𝑁,𝑁] is the number of the external points, 𝑤𝑗,𝑘(𝑥) is the support domain of the wavelet function, that is, 𝐴(𝑥)𝑢(𝑥)+𝐵(𝑥)𝑢(𝑥)+𝐶(𝑥)𝑢(𝑥)=𝑓(𝑥),𝑥𝑎,𝑏.(4.1).

It is easy to know that the quasi-Shannon interval wavelet is a linear combination of the quasi-Shannon scaling function 𝑢𝑗(𝑥)=2𝑗𝑛=0𝑢𝑗𝑥𝑛𝑤2𝑗𝑥𝑛,𝑛=0,1,2,,2𝑗.(4.2). Therefore, the quasi-Shannon interval wavelet function possesses all the properties of the quasi-Shannon scaling function.

It should be noted that the interval wavelet function (3.29) is similar to the result obtained in [5] using different methods. But we can construct different interval wavelets possessed similar format based on different extension technologies by this method. Otherwise, this method reveals the close relationship between the restricted variational principle and the interval interpolation wavelet.

4. Interval Wavelet Numerical Method for Ordinary Deferential Equations (ODEs)

4.1. Interval Wavelet Discrete Formulation of ODEs [6, 7]

Considering the following ODE:2𝑗𝑛=0𝑢𝑗(𝑥𝑛)𝐴𝑥𝑘22𝑗𝑤𝑥(𝑘𝑛)+𝐵𝑘2𝑗𝑤𝑥(𝑘𝑛)+𝐶𝑘𝑥𝑤(𝑘𝑛)=𝑓𝑘,𝑘=0,1,2,,2𝑗.(4.3) According to the collocation method, its approximate solution can be expressed as𝑥𝑘 Substituting (4.2) into (4.1), the system of algebraic equations can be obtained as follows:𝑀𝑈=𝐹,(4.4) The solution of (4.3) is an approximate one of (4.1) at the collocation point 𝑈. Equation (4.3) could be rewritten in the new form below as follows:𝐹 where 𝑀 is the vector of solution, 𝑀 is right side vector, and 𝑚𝑘,𝑛𝑥=𝐴𝑘22𝑗𝑤𝑥(𝑘𝑛)+𝐵𝑘2𝑗𝑤𝑥(𝑘𝑛)+𝐶𝑘𝑤(𝑘𝑛),(4.5) is the stiffness matrix. Any element of the matrix 𝑤2(𝑘𝑛)=𝑗cos𝜋(𝑘𝑛)exp(𝑘𝑛)2/(2𝑟2)𝑤(𝑘𝑛)(𝑏𝑎),𝑘𝑛0,𝑘=𝑛,2(𝑘𝑛)=2𝑗+1cos𝜋(𝑘𝑛)exp(𝑘𝑛)2/2𝑟2(𝑏𝑎)21(𝑘𝑛)2+1𝑟22,𝑘𝑛2𝑗3+𝜋2𝑟23𝑟2(𝑏𝑎)2,𝑘=𝑛.(4.6) could be expressed as follows:𝑆 where 𝑢(𝑥)|𝑆=𝛽(𝑥),𝑥𝑆.(4.7) In the solution of the above differential equations, there are several types of boundary conditions commonly encountered, such as Dirichlet and Neumann boundary conditions. For simple, we only consider the Dirichlet boundary in the following discussions.

Dirichlet boundary conditions of differential equations specify the value of function on a surface 𝑢(𝑥):𝑆 That is, the solution [0,1] on the boundary 𝑢𝑗(𝑥)=𝑢𝑗2(0)𝑤𝑗𝑥+2𝑗1𝑛=1𝑢𝑗𝑥𝑛𝑤2𝑗𝑥𝑛+𝑢𝑗2(1)𝑤𝑗𝑥2𝑗.(4.8) is known. Assuming that the definition interval of the solution is 𝜀𝑢(𝑥)+𝑢(𝑥)+𝑢(𝑥)=1,0<𝑥<1,𝑢(0)=𝑢(1)=0.(4.9), the approximate solution will be𝑢(𝑥)=1𝑒𝜆2𝑒𝜆2𝑒𝜆1𝑒𝜆1𝑥1𝑒𝜆1𝑒𝜆2𝑒𝜆1𝑒𝜆2𝑥+1,(4.10) Substituting (4.8) into (4.4), the system of algebraic equations satisfying the boundary conditions could be obtained.

4.2. Numerical Example

Considering the following two-point boundary problem as a convection equation,𝜆1 Its analytical solution as shown in Figure 1 is𝜆2 where 𝜆1=1+1+4𝜀2𝜀,𝜆2=11+4𝜀2𝜀.(4.11) and 𝑥=1 are Eigenvalues:𝜀 Obviously, the boundary layer phenomenon is evident near the boundary 2𝑗1𝑛=1𝑢𝑗𝑥𝑛22𝑗𝜀𝑤(𝑘𝑛)+2𝑗𝑤(𝑘𝑛)+𝑤(𝑘𝑛)=1,𝑘=1,2,,2𝑗1.(4.12), which would increase with the decrease of the parameter 𝑤. This could be used to test the effectiveness of the interval wavelet.

629253.fig.001
Figure 1: Analytical solution of the convection (4.9).

Substituting (4.5) into (4.9) results in its wavelet discrete formulation as the following: 𝜎 The weight function 𝐿,𝜀,and𝑗 in (4.12) can be set as the Shannon, the quasi-Shannon, or the interval wavelet functions. Their numerical solutions and corresponding errors are shown in Figure 2, where the count of collocation points is 65, and the points are distributed evenly in space. Obviously, the calculation precision of the interval wavelet is higher than others.

fig2
Figure 2: The comparison of the Shannon (a), the quasi-Shannon (b), and the interval wavelet (c), and their errors of numerical results (e), (f), and (g), respectively.

Besides, the regularized width parameter 𝐿 which has been discussed in [8], the parameters 𝑢 affect the calculation precision significantly. Therefore, appropriate values should be selected when implementing the numerical method developed.

4.2.1. Total Number of External Collocation Points 𝐿

In fact, the definition of interval wavelet in (3.29) indicates that the values of 𝐿 in external collocation points are obtained by Lagrange interpolation method. So, the parameter 𝐿=1 is relative with the smoothness of the solution, which is difficult to be determined. With the increasing of 𝐿=3, the Runge phenomena could occur in the Lagrange interpolation [9]. As a result, it could decrease the computational precision, which is proved by the numerical results shown in Table 1. When 𝐿=13, the maxima of both the relative error and the absolute error are apparently less than that of the other two methods. However, opposite results were obtained, when 𝐿(𝑗=6,𝜀=1). Actually, the solution of ODE with the Dirichlet boundary condition is not smooth on the boundary in most cases, and so the bigger L cannot make any contribution to improve the computational precision except increasing the time complexity of the numerical method. Consequently, 𝐿 might be the best choice.

tab1
Table 1: Influence of the total number of external collocation points (103).
4.2.2. Parameter 𝜀=1

Table 2 shows the influence of parameter 𝜀 on the calculation precision. It can be observed that the solution is smooth as 𝜀=0.001. Although the quasi-Shannon wavelet has local property, its computational precision is lower than that of the Shannon wavelet without local property. With 𝑗=6 decreasing, the superiority of quasi-Shannon wavelet method to Shannon wavelet method exhibits gradually. When 𝜀0.1, the maximum of the relative error of the quasi-Shannon wavelet method is only 0.33, while that of Shannon wavelet method is as high as 26.6. The precision of the interval wavelet method keeps higher than the other two methods in the whole process. When 𝜀 and 𝑗, the maxima of both the relative error and the absolute error of the interval wavelet method are the lowest among the three methods.

tab2
Table 2: Influence of parameter 𝜀.

However, when the parameter 𝑗 further decreases, the relative error of the interval wavelet method is slightly larger than that of quasi-Shannon wavelet method, and the absolute error of the interval wavelet method approximately identical to that of quasi-Shannon wavelet method. This is due to the boundary layer phenomena of the two-point boundary problem appearing at the points adjacent to the boundary 𝑗. The boundary layer phenomena result in the large gradient of solution, and thus the Gibbs phenomena [10] are unavoidable. Since the construction of the interval wavelet is based on the Lagrange continuation, the Gibbs phenomena would decrease the precision of continuation. Hence, the computational precision of quasi-Shannon interval wavelet method is decreased.

4.2.3. Influence of 𝑗

The parameter 2𝑗+1 is used to control the total number of collocation points in the discrete formula of wavelet collocation method. The number of collocation points is the exponential function of the parameter 𝜀, that is, 𝑗=10. According to the explanation in Section 4.2.2, if there is no Gibbs phenomena or the Gibbs phenomena are weak enough, the precision of the interval wavelet method should be slightly higher than the quasi-Shannon wavelet method, even the value of 𝑗(𝜀=0.001) is small. The Gibbs phenomena would be continuously weakening up to none with the increment of the collocation points. Table 3 indicates that calculation errors of all three methods decrease with number of the collocation points increasing. When 𝑗, the maximum absolute error of the interval wavelet method is less than the other two methods. In addition, the maximum relative error of the interval wavelet method is less than the error of the quasi-Shannon wavelet method, but larger than the error of Shannon wavelet method.

tab3
Table 3: Influence of (103).
4.3. Comparison of Calculation Precision Obtained by Shannon Interval Wavelet Method and Quasi-Shannon Interval Wavelet Method

In this work, the Shannon interval wavelet is also constructed to further evaluate the performance of the quasi-Shannon interval wavelet for numerical analysis. Instead of (3.16), the following weight function is used to construct the Shannon interval wavelet:𝜀 The first and second derivatives of (4.13) are expressed, respectively, as follows: 𝜀0.01𝑗=10 The collocation method of Shannon interval wavelet method could be obtained by substituting (4.13)–(4.15) into (4.5).

Table 4 indicates that the maxima of both the relative error and the absolute error obtained by the Shannon interval wavelet method are less than those obtained by the quasi-Shannon interval wavelet method when the solution curve of the two-point problem is smooth. As 𝑗 decreasing, the precision of quasi-Shannon interval method become higher than that of the Shannon interval wavelet method. Especially when 𝜀, the calculation precision obtained by quasi-Shannon interval wavelet method is higher than that obtained by Shannon interval wavelet method in the whole process with the collocation number varying from a smaller number to 1025 (i.e., (105)).

tab4
Table 4: Comparison of calculation precision obtained by Shannon and quasi-Shannon interval wavelet methods.

Figure 3 also shows that when the equation is solved by the quasi-Shannon interval wavelet, the numerical error mainly lies in the boundary, while the errors of other parts are very little. On the contrast, when the Shannon interval wavelet method is utilized, the numerical error almost evenly distributes in the all domain. Consequently, the sum of the error is much larger than that of the quasi-Shannon interval wavelet.

fig3
Figure 3: Comparison of the quasi-Shannon interval wavelet and the Shannon interval wavelet (𝐿,𝜀, 𝑗, and ).

5. Conclusions

A novel general method for construction of interval wavelet based on the restricted variational principle is proposed in this paper. Both Shannon and quasi-Shannon interval wavelets could be constructed using this new method. With appropriate values of parameters , and selected, the numerical results for a differential equation show that the quasi-Shannon interval wavelet outperforms than Shannon and quasi-Shannon wavelets. Furthermore, the capability of quasi-Shannon interval wavelet for numerical analysis is better than that of Shannon interval wavelet.

It should be noted that combining the corresponding interval wavelet with the classical Galerkin method can obtain an excellent wavelet-Galerkin method for engineering problems. Since quasi-Shannon wavelet is orthogonal, this guarantees the stability of the matrix equations in Galerkin method.

Acknowledgment

This work is supported by the National Natural Science Foundation of China (no. 60772038).

References

  1. A. Cohen, I. Daubechies, and P. Vial, “Wavelets on the interval and fast wavelet transforms,” Applied and Computational Harmonic Analysis, vol. 1, no. 1, pp. 54–81, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  2. E. Quak and N. Weyrich, “Decomposition and reconstruction algorithms for spline wavelets on a bounded interval,” Applied and Computational Harmonic Analysis, vol. 1, no. 3, pp. 217–231, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  3. C. R. Ortloff, “Restricted variational principle methods for boundary value problems of kinetic theory,” Physics of Fluids, vol. 10, no. 1, pp. 230–231, 1967. View at Publisher · View at Google Scholar
  4. G. W. Wei, “Quasi wavelet and quasi interpolating wavelets,” Chemical Physics Letters, vol. 296, no. 3-4, pp. 215–222, 1998. View at Publisher · View at Google Scholar
  5. O. V. Vasilyev and S. Paolucci, “A dynamically adaptive multilevel wavelet collocation method for solving partial differential equations in a finite domain,” Journal of Computational Physics, vol. 125, no. 2, pp. 498–512, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  6. S.-L. Mei, Q.-S. Lu, S.-W. Zhang, and L. Jin, “Adaptive interval wavelet precise integration method for partial differential equations,” Applied Mathematics and Mechanics, vol. 26, no. 3, pp. 364–371, 2005. View at Publisher · View at Google Scholar · View at MathSciNet
  7. S.-L. Mei, Q.-S. Lu, and S.-W. Zhang, “An adaptive wavelet precise integration method for partial differential equations,” Chinese Journal of Computational Physics, vol. 21, no. 6, pp. 523–530, 2004. View at Google Scholar
  8. D.-C. Wan and G.-W. Wei, “The study of quasi wavelets based numerical method applied to Burgers' equations,” Applied Mathematics and Mechanics, vol. 21, no. 10, pp. 1099–1110, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  9. J. J. Leader, Numerical Analysis and Scientific Computation, Addison-Wesley, Reading, Mass, USA, 2004.
  10. I. Daubechies, Ten Lectures on Wavelets, vol. 61 of CBMS-NSF Regional Conference Series in Applied Mathematics, SIAM, Philadelphia, Pa, USA, 1992. View at Zentralblatt MATH · View at MathSciNet