Available online at www.sciencedirect.com
Physics Procedia 33 (2012) 1346 – 1353
2012 International Conference on Medical Physics and Biomedical Engineering
Construction of a type of Biorthogonal Binary Wave Functions by Multiresolution Analysis Approach Delin Hua School of Education Nanyang Institute of Technology Nanyang 473000, China
[email protected]
Abstract Wavelet analysis is nowadays a widely used tool in applied mathematics. The concept of vector-valued binary wavelets with two-scale dilation factor associated with an orthogonal vector-valued scaling function is introduced. The existence of orthogonal vector-valued wavelets with two-scale is discussed. A necessary and sufficient condition is provided by means of vector-valued multiresolution analysis and paraunitary vector filter bank theory. An algorithm for constructing a sort of orthogonal vector-valued wavelets with compact support is proposed, and their orthogonal properties are investigated.
©2012 2011Published Published by Elsevier Ltd. Selection and/or peer-review under responsibility of [name Committee. organizer] © by Elsevier B.V. Selection and/or peer review under responsibility of ICMPBE International Index Terms—Dilation equation; vector-valued multiresol-ution analysis; binary scaling functions; binary wavelets; matrix theory; orthogonal; frame wavelets
Introduction Construction of wavelet bases is an important aspect of wavelet analysis, and there are many kinds of scalar scaling functions and scalar wavelet functions. The main advantage of wavelets is their timefrequency localization property. Already they have led to exciting applications in signal analysis [1], fractals [2], image processing [3] and so on. Sampling theorems play a basic role in digital signal processing. They ensure that continuous signals can be processed by their discrete samples. Vector-valued wavelets are a class of generalized multiwavelets [4]. Chen and Cheng [5] introduced the notion of vectorvalued wavelets and showed that multiwavelets can begenerated from the component functions in vectorvalued wavelets. Vector-valued wavelets and multi-wavelets are different in the following sense. In real life, Video images are vector-valued signals. Vector-valued wavelet transforms have been recently studied
1875-3892 © 2012 Published by Elsevier B.V. Selection and/or peer review under responsibility of ICMPBE International Committee. doi:10.1016/j.phpro.2012.05.221
1347
Delin Hua / Physics Procedia 33 (2012) 1346 – 1353
for image coding by W. Li. Chen and Cheng studied orthogonal compactly supported vector-valued wavelets with 2-scale. Inspired by [5-7], we are about to investigate the construction of a class of orthogonal compactly supported vector-valued wavelets with three-scale. Similar to uni-wavelets, it is more complicated and meaningful to investigate vector-valued wav-elets with 4-scale. Based on an observation in [5,6], another purpose of this article is to introduce the notion of orthogonal vector-valued wavelet packets with two-scale and investigate their properties. Notations and multiresolution analysis Set s be a constant and 2 s Z , where Z and Z den--ote all non-negative integers and all integers, respectively. By L2 ( R 2 , C v ) , we denote the aggregate of arbitrary vector-valued functions F (t ) ,i.e, L2 ( R 2 , C s ) : { F (t ) ( f1 (t ), f 2 (t ), , f s (t ) )T : f (t ) L2 ( R 2 ), 1, 2, , s }, where T means the transpose of a vector. For example, video images an -d digital films are examples of vectorth column at time t. For F (t ) L2 ( R 2 , C v ), valued functions where g (t ) denotes the pixel on the u 2 1/ 2 F denotes the norm of vector-valued function F (t ), i.e., F : ( | f (t ) | dt ) , and its 1 R integration is defined as R F (t ) dt : ( R2 f1 (t ) dt , f 2 (t ) dt , , R T f ( t ) dt ) .The Fourier transform of F (t ) is defined by s 2 2
2
2
R
F( ) :
R
2
F (t ) e
it
dt .
For two vector-valued functions F , G F ( ), G ( ) : F (t ) G (t )* dt , 2
L2 ( R 2 , C v ) , their symbol inner product is defined by
R
(1)
where * means the transpose and the complex conjugate, and I u denotes the v v identity matrix. A sequence {Fl (t )}l Z U L2 ( R 2 , C u ) is called an orthonormal set of the subspace U , if the following condition is satisfied
Fj ( ), Fk ( )
I , j, k
j ,k v
Z,
(2)
where j , k is the Kronecker symbol, i.e., j , k 1 as j k and j , k 0 otherwise. Definition 1. We say that F (t ) Y L2 ( R 2 , C s ) is an orthogonal vector-valued function of Y if its translates { F (t v)}v Z 2 is an orthonormal set of the subspace Y , i.e.,
F(
n), F (
v)
I , n, v Z 2 .
n ,v u
(3)
Definition 2[5]. A sequence {Fv (t )}v Z 2 U L2 ( R 2 , C s ) is called an orthonormal basis of Y , if it satisfies (2), and for any G (t ) U , there exists a unique sequence of s s constant matrices { M k }k Z 2 such that
G (t ) v Z2
M v Fv (t ).
(4)
1348
Delin Hua / Physics Procedia 33 (2012) 1346 – 1353 [5]
Definition 3 . A vector-valued multiresolution analysis of L2 ( R 2 , C s ) is a nested sequence of closed 2 2 s subspaces {Vl }l z such that (i) Vl Vl 1 , l Z ; (ii) l z Vl {0}; l z Vl is dense in L ( R , C ) , 2 2 s l where 0 is the zero vector of L ( R , C ) ; (iii) F (t ) V0 if and only if F (2 t ) Vl , ; (iv) there is (t ) V0 such that the sequence { (t v ), v Z 2 } is an orthonormal basis of subspace V0 . On the basis of Definition 2 and Definition 3, we obtain (t ) satisfies the following equation
(t ) u Z2
Pv (2t u ),
(5)
constant matrices, i.e., {Pk }k z has only Where {Pk }k Z 2 , is a finite supported sequence of finite non-zero terms, and the others are zero matrices. By taking the Fouries transform for the both sides of (5), and assuming ( ) is continuous at zero, we have
(2 )
R2 ,
( ) ( ),
4 ( ) u Z2
(6)
Pu exp{ iu }.
(7)
Let W j ( j Z ) denote the orthocomplement subspace of V j in V j 1 and there exist three vectorvalued functions Gs (t ) L2 ( R 2 , C s ), r 1, 2,3 such that their translations and dilations form a Riesz basis of W j , i.e., j 2 W j clos (span{Gs (2 t u ) : s 1, 2, 3; u Z }), j Z . ( 8) L ( R ,C ) Since r (t ) W0 V1 , r 1, 2,3 there exist three finitely supported sequences {Bk( s ) }k z of v v constant matrices such that 2
r
u
(t ) u Z2
4
(s)
Bu( r ) (2t u ), r 1, 2,3.
( ) u Z
2
(9)
Bu( s ) exp{ iu }.
(10)
Then, the refinement equation (10) becomes thefollowing
(2 )
r
(r )
( ) ( ), s 1, 2,3,
R2 .
(11)
If (t ) L2 ( R 2 , C s ) is an orthogonal vector-valued scaling function, then it follows from (3) that
( ),
(
v)
v Z 2.
I ,
0,u v
(12)
We say that r (t ) L2 ( R 2 , C s ), r onal vector-valued wavelet functions associated with the vector-valued scaling function (t ) , if they satisfy
(
n),
and the family { r
( ),
(
r r
(
n)
Lemma 1 . Let F (t )
F( k Z
O, r 1, 2,3, n, v Z ,
(13)
(t v), r 1, 2,3, v Z } is an orthonomal basis of W0 . Thus we have
[6]
2
v)
1, 2,3 are orthog
2k ) F (
r,
I , r,
0, n u 2
2
1, 2,3; n Z .
(14)
L ( R , C ) . Then F (t ) is an orthogonal vector-valued function if and only if
2k )*
s
Is ,
R2 .
(15)
1349
Delin Hua / Physics Procedia 33 (2012) 1346 – 1353
Lemma 2. If (t ) L2 ( R 2 , C s ) , defined by (5), is an orthogonal vector-valued scaling function, then for v Z , we have the following equalities,
P (P
4v
)*
4
I.
(16)
0, v s
Z 3
(
) (
)
*
(17)
2
Is ,
R .
0
Proof. By substituting equation (5) into the relation (12), for
(
I
0, k s
k ), ( ) 2
l Z u Z
1 4
2
l Z u Z
2
(
Pl
2k l ), (
R
2
k
Z , we obtain that
Pl (2t 2k l ) (2t u )* ( Pu )* dt
u ) ( Pu )*
1 4u
Z
2
Pu ( Pu
4k
)* .
Thus, both Theorem 1 and formulas (16), (23) and (24) provide an approach to design a class of compactly supported orthogonal vector-valued wavelets. Design of vector-valued wavelets
In the following, we begin with considering the existence of a class of orthogonal vector- valued wavelets.
compactly
supported
(t ) L2 ( R 2 , C s ) defined by (5), be an orthogonal vector-valued scaling function. Assume r (t ) L2 ( R 2 , C s ), r 1, 2,3, and ( ) and (r ) ( ) are defined by (7) and (10), respectively. Then r (t ) are orthogonal vector-valued wavelet functions associated with (t ) if and only if
Theorem 1. Let
3
(
(r )
)
)*
(
O, r 1, 2,3
(18)
0 3
(r )
(
( )
)
)*
(
r,
Is ,
(19)
0
where r , s 1, 2, l Z
l Z
2
2
Pl ( Bl( r 2) u )* Bl( r ) ( Bl(
R . or equivalently,
O, r 1, 2,3, u Z 2 ;
) * 2u
)
4
I r , s 1, 2, 3, u Z 2 .
r,
0,u v ,
(20) (21)
Proof. Firstly, we prove the necessity. By Lemma 1 and (6), (11) and (13), we have
(2
O u Z
u )
s
(2
u )*
2
(
u ) (
u )
u Z2
u )*
( 3
( 0
(r )
/ 2)
u )*
( (r )
(
/ 2)* .
1350
Delin Hua / Physics Procedia 33 (2012) 1346 – 1353
It follows from formula (14) and Lemma 1 that
Iv
r,
r u Z
(2
2u )
2u )*
(2
2
(r )
(
u ) (
u )
2
u Z
u )*
( 3
(r )
( )
(
u )*
( ( )
)
)* .
(
0
Next, the sufficiency of the theorem will be proven. From the above calculation, we have
(2 u Z
u )
r
(2
u )*
2
3
(
/ 2)
(r )
/ 2)*
(
I .
0, r u
0
Furthermore
( ), r
r
r
(
1
k)
2
2u )* e 4ik d 1 ( k) 2
(2
( ),
[0, / 2]2
(2
2u )
u Z2
Z2
O, r 1, 2,3, k [0, / 2]2
r
(2
2u )
u Z2
2u )* e 4ik d Z2 0, k r , I s , k 1, 2,3} are a family of Thus, (t ) and r (t ), r 1, 2,3 are mutually orthogonal, and { (t ), orthogonal vector-valued functions. This shows the orthogonality of { s ( v ), r 1, 2,3}v z . Similar to [7,Proposition 1], we can prove its completeness in W0 . Theorem 2. Let (t ) L2 ( R 2 , C s ) be a 5-coefficient compactly supported orthogonal vector-valued (2
scaling functions satisfying the following refinement equation:
(t )
P0 (2t ) P1 (2t 1)
P4 (2t 4).
Assume there exists an integer 4 , such that (4 I v P ( P )* ) 1 P ( P )* is a positive definite matrix. Define Qr ( r 1, 2,3) to be two essentially distinct Hermitian matrice, which are all invertible and satisfy
(Qr ) 2
P ( P )* ] 1 P ( P )* .
[4 I s
(22)
Define B (j r ) B
Then r
Qs Pj ,
(r) j
(t )
j
1
(Qs ) Pj , r
j
r 1, 2,3, j {0,1, 2,3, 4}.
(t ) ( r 1, 2,3) , defined by (24), are orthogonal vector-valued wavelets associated with (r)
B0
(2t )
(r)
B1
(2t 1)
(r)
B4
(2t
4)
(23)
(t ) : (24)
1351
Delin Hua / Physics Procedia 33 (2012) 1346 – 1353
1 . By Lemma 2, (20) and (21), it suffices to show that {B0( r ) , B1( r ) ,
Proof. For convenience, let
B2( r ) , B3( r ) , B4( r ) , r 1, 2,3} satisfy the following equations: P0 ( B4( r ) )*
O,
(r ) * 4 0 (r ) * 0 0 (r ) ( ) * 0 4 (r ) (r ) * 0 0
P (B ) P (B )
B (B )
B (B )
r 1, 2,3,
O,
(25)
r 1, 2,3, (r ) * 1
P1 ( B ) O,
P4 ( B )
r,
(r ) 1
(26) (r ) * 4
O,
(27)
{1, 2,3} ,
(s) * 1
(28) (s) 4
B (B )
(s) * 4
B (B )
4I s .
(29)
If {B0( r ) , B1( r ) , B2( r ) , B3( r ) , r 1, 2} are given by (23), then equations (26), (27) and (29) follow from (16). For the proof of (28) and (30), it follows from (16) and (27) that
P0 ( B0( r ) )*
P1 ( B1( r ) )*
P2 ( B2( r ) )*
[ P0 ( P0 )*
P2 ( P2 )*
( P1 ( P1 )*
P1 ( P1 )* )( Bs )
B0r ( B0r )*
P3 ( P3 )*
B1r ( B1r )*
1
P3 ( B3( r ) )*
P4 ( P4 )* ] Qs
P4 ( B4( r ) )*
P1 ( P1 )* (Qs )
1
O.
B2r ( B2r )*
B3r ( B3r )*
B4r ( B4r )*
Qr {P1 (P1 )* [P1 (P1 )* ] 1[4Is P1 (P1 )* ]P1 (P1 )*}(Qr ) 1 Qr {P1 ( P1 )*
4I s
P1 ( P1 )* }(Qr )
1
4I s .
So, (28), (30) follow. This completes the proof of Thm 2. Example1. Let
(t )
L2 ( R 2 , C 3 ) be 5-coefficient orthogonalvector valued scaling function satisfy the
following equation:
(t )
P0 (2t ) P1 (2t 1)
Where P3
P4
O, P0 ( P4 )*
P4 (2t 4). O,
P0 ( P0 )*
P0
P1 ( P1 )* 2 2 1 2
2 2 2 3
0
0
P2 ( P2 )*
P3 ( P3 )*
0 1 2 3 3
, P1
P4 ( P4 )*
4I3 .
1
0
0
0
2 6
0 ,
0
0
3 3
1352
Delin Hua / Physics Procedia 33 (2012) 1346 – 1353
P2
2 Q1
2 2 2 3
0
0
0
2 3 3 2
.
0
53
0
4
0
53
0
, Q2
2
0
0
2
0
53
0
1 ,
0
53
0
0
1 . By using (22), we can choose
Suppose 2
2 2 1 2
2
0
4
By applying formula (24), we get that 1
1
2
2 53
(1)
B0
106 0
0
106
53
159
53
2 (1)
, B1
106
0
0
0
6
B2(1)
1 2 53 106 0
B0(2)
0
0 2
6 3
1 2 106 159 0
1 2 53 106
0
6
6
0
0
1 2 106 159 0
0 53 , 53 6 6 0 53 . 53 6 6
Delin Hua / Physics Procedia 33 (2012) 1346 – 1353
2 B1(2)
B2(2)
0
0
0
106 6
0
0
0
1 2 53 106
1 2 106 159
0
0
,
2 6 3 0 53 . 53 6 6 ( )
(2t ) Applying Theorem 2, we obtain that G (t ) B0 B1( ) (2t 1) B4( ) (2t 4), 1, 2,3 are orthogonal vector-valued wavelet functions associated with the orthogonal vector-valued scaling function Conclusion
A necessary and sufficient condition on the existence of a class of orthogonal vector-valued wavelets is presented. An algorithm for constructing a class of compactly supported orthogonal vectorvalued wavelets is proposed. References [1]
N. Zhang X.Wu,``Lossless of color masaic images". IEEE Trans Image Delivery}, 2006, 15(6), pp. 1379-1388
[2]
S. Efromovich, J. akey , M. Pereyia, N. Tymes, ``Data-Diven and Optimal Denoising of a Signal and Recovery of its
Derivation Using Multiwavelets", IEEE Trans Signal Processing, 2,004, 52(3), pp.628-635 [3]
Z. Shen , Nontensor product wavelet packets in\ L2 ( R s ) [J]. SIAM Math Anal.1995, 26(4): 1061-1074.
[4]
X. G. Xia, B. W. Suter, `Vector-valued wavelets and vector filter banks", IEEE Trans Signal Processing}, 1996, 44(3),
pp. 508-518. [5] Solitons
Q. Chen, Z..Cheng, ``A study on compactly supported orthogonal vector-valued wavelets andwavelet packets”, Chaos,
& Fractals}, 2007, 31(4), pp.1024-1034.
[6]
C.K.Chui, J.Lian, ``A study of orthonormal multiwavelets". Appli Numer Math, 1996, 20(1), pp. 273-298
[7]
S. Yang. Z.Cheng, H. Wang, ``Construction of biorthogonal multiwavelets", Math Anal Appl}, 2002, 276(1), pp.1-12.
[8]
S. Li and H. Ogawa., “Pseudoframes for Subspaces with Applications,” J. Fourier Anal. Appl.., Vol. 10,No 4, pp.409-
431, 2004. [9]
S. Yang , Biorthogonal two-direction refinable function and two-direction wavelet. Appl Math & Comput . 2006 , 1822
(3) : 1717-1724 [10]
Q. Chen , A . Huo , “The research of a class of biorthogonal compactly supported vector-valued wavelets , Chaos ,
Solitons & Fractals , 41 (2 ), pp . 951-961, August 2009”
1353