1. Introduction
Spline wavelet is one of the most important wavelets in the wavelet family. In both applications and wavelet theory, the spline wavelets are especially interesting because of their simple structure. All spline wavelets are linear combination of B-splines. Thus, they inherit most of the properties of these basis functions. The simplest example of an orthonormal spline wavelet basis is the Haar basis. The orthonormal cardinal spline wavelets in were first constructed by Battle [2] and Lemarié [3] . Chui and Wang [4] found the compactly supported spline wavelet bases of and developed the duality principle for the construction of dual wavelet bases [1] [5] .
Wavelets are a fairly simple mathematical tool with a variety of possible applications. If is an orthonormal basis of, then is called a wavelet. Usually a wavelet is derived from a given multiresolution analysis of. The construction of wavelets has been discussed in a great number of papers. Now, considerable attention has been given to wavelet packet analysis as an important generalization of wavelet analysis. Wavelet packet functions consist of a rich family of building block functions and are localized in time, but offer more flexibility than wavelets in representing different kinds of signals. The main feature of the wavelet transform is to decompose general functions into a set of approximation functions with different scales. Wavelet packet transform is an extension of the wavelet transform. In wavelet transformation signal decomposes into approximation coefficients and detailed coefficients, in which further decomposition takes place only at approximation coefficients whereas in wavelet packet transformation, detailed coefficients are decomposed as well which gives more wavelet coefficients for further analysis.
For a given multiresolution analysis and the corresponding orthonormal wavelet basis of, wavelet packets were constructed by Coifman, Meyer and Wickerhauser [6] [7] . This construction is an important generalization of wavelets in the sense that wavelet packets are used to further decompose the wavelet components. There are many orthonormal bases in the wavelet packets. Efficient algorithms for finding the best possible basis do exist. Chui and Li [8] generalized the concept of orthogonal wavelet packets to the case of nonorthogonal wavelet packets. Yang [9] constructed a scale orthogonal multiwavelet packets which were more flexible in applications. Xia and Suter [10] introduced the notion of vector valued wavelets and showed that multiwavelets can be generated from the component functions in vector valued wavelets. In [11] , Chen and Cheng studied compactly supported orthogonal vector valued wavelets and wavelet packets. Other notable generalizations are biorthogonal wavelet packets [12] , non-orthogonal wavelet packets with r-scaling functions [13] .
The outline of the paper is as follows. In, we introduce some notations and recall the concept of B-splines and wavelets. In, we discuss the B-spline wavelet packets and the corresponding dual wavelet packets.
2. Preliminaries
In this Section, we introduce B-spline wavelets (or simply B-wavelets) and some notions used in this paper.
Every mth order cardinal spline wavelet is a linear combination of the functions. Here the function is the mth order cardinal B-spline. Each wavelet is constructed by spline multiresolution analysis. Let m be any positive integer and let denotes the mth order B-spline with knots at the set of integers such that
The cardinal B-splines are defined recursively by the equations
We use the following convention for the Fourier transform,
The Fourier transform of the scaling function is given by
(1)
For each, we set, and for each, let denotes the -closure of the algebraic span of. Then is said to generate spline multiresolution analysis if the following conditions are satisfied.
1)
2);
3),
4) for each, is a Riesz basis of.
Following Mallat [14] , we consider the orthogonal complementary subspaces that is;
5).
6).
7).
These subspaces, are called the wavelet subspaces of relative to the B-spline. Since and, we have
, (2)
where is some sequence in. Taking the Fourier transform on both sides of (2), we obtain
. (3)
Substituting the value of from (1) into (3), we have
.
This gives
(4)
So, (2) can be written as
,
which is called the two scale relation for cardinal B-splines of order.
Chui and Wang [1] , introduced the following mth order compactly supported spline wavelet or B-wavelet
, (5)
with support that generates and consequently all the wavelet spaces. To verify that is in, we need the spline identity
. (6)
So, substituting (6) into (5), we have the two scale relation
, (7)
where,
. (8)
Let
, (9)
with the corresponding two scale sequence. If is a wavelet, then there exists another called the dual wavelet of such that
. (10)
For the scaling function, we define its dual by
, (11)
such that
. (12)
Now, we have
(13)
Taking the Fourier transform of (13), we have
(14)
where,
(15)
A necessary and sufficient condition for the duality relationship (12) is that and are dual two scale symbols in the sense that
. (16)
A proof of this statement is given in ( [15] , Theorem 5.22). Also from (7) and (9), we have
(17)
where,
(18)
We observe that
(19)
See ( [15] , Section 5.3).
If is an orthogonal scaling function, then
. (20)
We say that is orthogonal (o.n) B-wavelet function associated with orthogonal scaling function if
, (21)
and is an orthonormal basis of, so we have
. (22)
Lemma 1 Let. Then is an orthonormal family if and only if
. (23)
Proof See ([15] , page no. 75].
Theorem 1 Let defined by (13) is an orthonormal scaling function. Assume that whereas and are defined by (15) and (18) respectively. Then is an orthonormal wavelet function associated with if and only if
(24)
Proof Let us suppose that is an orthonormal wavelet function associated with. By Lemma 1 and (21), we have
Again by Lemma 1 and (22), we have
On the other hand, let (24) holds.
Now,
Also,
Thus, and are orthogonal and is an orthonormal wavelet function associated with.
3. B-Spline Wavelet Packets and Their Duals
Following Coifman and Meyer [6] [7] , we introduce two sequences of functions and defined by
(25)
(26)
where
When and, we have
and for and, we have
We call the sequence of B-spline wavelet packets induced by the wavelet and its corresponding scaling function whereas denotes the corresponding sequence of dual wavelet packets. By applying the Fourier transformation on both sides of (25), we have
(27)
where,
(28)
(29)
So, (24) can be written as
(30)
Similarly, taking the Fourier transformation on both sides of (26), we have
(31)
where,
(32)
(33)
Using these conditions we can write
(34)
We are now in a position to investigate the properties of B-spline wavelet packets.
Theorem 2 Let be any orthonormal scaling function and its corresponding family of B-spline wavelet packets. Then for each, we have
(35)
Proof Since satisfies (35) for. We may proceed to prove (35) by induction. Suppose that (35) holds for all, where, a positive integer and. We have
, where denote the largest integer not exceeding. By induction hypothesis and Lemma 1, we have
(36)
By using (27), (30) and (36), we obtain
Hence, by Lemma 1, (35) follows.
Theorem 3 Let be a B-spline wavelet packet with respect to the orthonormal scaling function. Then for every, we have
(37)
Proof By (27), (30) and (36), for we have
For the family of B-spline wavelet packets corresponding to some orthonormal scaling function, consider the family of subspaces
generated by. We observe that
where is the MRA of generated by, and is the sequence of orthogonal complementary (wavelet) subspaces generated by the wavelet. Then the orthogonal decomposition
may be written as
.
A generalization of the above result for other values of can be written as
.
Theorem 4 For the B-spline wavelet packets, the following two scale relation
(38)
holds for all.
Proof In order to prove the two scale relation, we need the following identity, see ([15] , Lemma 7.9)
(39)
Taking the right-hand side of (38), and applying the identity (39), we have
(40)
This completes the proof of the theorem.
Next, we discuss the duality properties between the wavelet packets and.
Lemma 2 For all and,
. (41)
Proof We will prove (41) by induction on. The case is the same as our assumption (12) on the dual scaling functions and. Suppose that (41) holds for all where,
where is a positive integer. Then for, we can write for some,
according to the proof of Theorem 7.24 in [15] . From the Fourier transform formulations of equations (25) and (26) and using (34) we have
Since, it follows from the induction hypothesis that for all, and this is equivalent to
(42)
Thus, we have
This shows that (41) also holds for.
Lemma 3 For all and, and, with,
(43)
Proof By applying the Fourier transform formulations of Equations (25) and (26) and using (42) and (34), we have as in the proof of Lemma 2 that