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.

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

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

where is some sequence in. Taking the Fourier transform on both sides of (2), we obtain

Substituting the value of from (1) into (3), we have

.

This gives

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

with support that generates and consequently all the wavelet spaces. To verify that is in, we need the spline identity

So, substituting (6) into (5), we have the two scale relation

where,

Let

with the corresponding two scale sequence. If is a wavelet, then there exists another called the dual wavelet of such that

For the scaling function, we define its dual by

such that

Now, we have

Taking the Fourier transform of (13), we have

where,

A necessary and sufficient condition for the duality relationship (12) is that and are dual two scale symbols in the sense that

A proof of this statement is given in ( [15] , Theorem 5.22). Also from (7) and (9), we have

where,

We observe that

See ( [15] , Section 5.3).

If is an orthogonal scaling function, then

We say that is orthogonal (o.n) B-wavelet function associated with orthogonal scaling function if

and is an orthonormal basis of, so we have

Lemma 1 Let. Then is an orthonormal family if and only if

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

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.

Following Coifman and Meyer [6] [7] , we introduce two sequences of functions and defined by

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

where,

So, (24) can be written as

Similarly, taking the Fourier transformation on both sides of (26), we have

where,

Using these conditions we can write

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

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

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

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

holds for all.

Proof In order to prove the two scale relation, we need the following identity, see ([15] , Lemma 7.9)

Taking the right-hand side of (38), and applying the identity (39), we have

This completes the proof of the theorem.

Next, we discuss the duality properties between the wavelet packets and.

Lemma 2 For all and,

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

Thus, we have

This shows that (41) also holds for.

Lemma 3 For all and, and, with,

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