Next Article in Journal
Extension of Fuzzy ELECTRE I for Evaluating Demand Forecasting Methods in Sustainable Manufacturing
Next Article in Special Issue
Higher-Order Benjamin–Ono Model for Ocean Internal Solitary Waves and Its Related Properties
Previous Article in Journal
The Probabilistic Dual Hesitant Fuzzy Multi-Attribute Decision-Making Method Based on Cumulative Prospect Theory and Its Application
Previous Article in Special Issue
Existence and General Energy Decay of Solutions to a Coupled System of Quasi-Linear Viscoelastic Variable Coefficient Wave Equations with Nonlinear Source Terms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quadratic Phase Multiresolution Analysis and the Construction of Orthonormal Wavelets in L2(ℝ)

1
Department of Mathematics, IIT Patna, Patna 801106, Bihar, India
2
Department of Mathematics, Texas A&M University-Kingsville, Kingsville, TX 78363-8202, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Axioms 2023, 12(10), 927; https://doi.org/10.3390/axioms12100927
Submission received: 2 September 2023 / Revised: 23 September 2023 / Accepted: 25 September 2023 / Published: 28 September 2023

Abstract

:
The multi-resolution analysis (MRA) associated with quadratic phase Fourier transform (QPFT) serves as a tool to construct orthogonal bases of the L 2 ( R ) . Consequently, it assumes a pivotal role in facilitating potential applications of QPFT. Inspired by the sampling theorem applicable to band-limited signals in the QPFT domain, this paper formulates the development of the MRA linked with QPFT. Subsequently, we develop a method for constructing orthogonal bases for L 2 ( R ) , followed by some examples.

1. Introduction

The Fourier transform (FT) is a remarkable discovery in the field of mathematical sciences, which has had a profound impact on many branches of science and engineering [1,2]. Over time, the domain of Fourier analysis has witnessed numerous mathematical breakthroughs, leading to significant advancements and profound implications of the classical Fourier transform. Notable developments that stem from the conventional Fourier transform include the fractional Fourier transform [3,4], linear canonical transform [5], special affine Fourier transform [6], and the relatively recent quadratic-phase Fourier transform [7]. The quadratic-phase Fourier transform (QPFT) extends the classical Fourier transform, incorporating quadratic phase factors into its kernel. In the QPFT, the signal’s time-domain representation is multiplied by a quadratic phase term before computing its frequency-domain representation. This additional quadratic phase term allows for a more flexible analysis of signals with time-varying frequency content. The QPFT provides a unified framework for handling transient and non-transient signals, making it particularly useful for analyzing signals with time-varying properties. It has found applications in various fields, including signal processing, time-frequency analysis, and communication systems. The mathematical expression for the QPFT involves five real parameters ( A , B , C , D , E ) that control characteristics of the quadratic phase term. Adjusting these parameters can tailor the QPFT to suit specific signal processing requirements. Overall, the quadratic-phase Fourier transform enhances the traditional Fourier transform’s capabilities, enabling a more versatile and powerful analysis of signals in the time-frequency domain. Shah et al. studied short-time quadratic-phase Fourier transform as well as quadratic-phase wavelet transform (QPWT) with many applications in [8,9]. Also, the quadratic phase Fourier wavelet transform was explored by Prasad and Sharma in [10]. Over the past decades, various integral transforms, including the Fourier, fractional Fourier, and linear canonical transforms, have been extensively explored in time-frequency analysis.
On the other hand, multiresolution analysis (MRA) is a powerful mathematical framework used in signal and image processing introduced by Mallat [11], particularly in the field of wavelet analysis. It provides a systematic and hierarchical approach to analyzing signals or images at different levels of detail or resolution. The concept of MRA is rooted in the idea of representing a signal or an image in terms of a series of subspaces, each capturing different levels of frequency or scale information. Madych [12] established elementary properties of MRA in L 2 ( R n ) with scaling functions represented as characteristic functions. Subsequently, Zhang [13] explored scaling functions and wavelets in standard MRA, providing a characterization of the support of the Fourier transform of these scaling functions. Malhotra and Vashisht [14] contributed to understanding scaling functions on Euclidean spaces. The MRA associated with FrWT was also introduced in [15], where FrWT analyzes signals in the time-frequency-FrFD domain. Ahmad [16] studied fractional MRA and associated scaling functions in L 2 ( R ) . Dai et al. [17] proposed a novel fractional wavelet transform (FRWT) and studied MRA associated with the developed FRWT, together with the construction of the orthogonal fractional wavelets. Shah and Lone [18,19] studied special affine MRA and the construction of orthonormal wavelets in L 2 ( R ) and studied Shannon’s sampling theorem for the quadratic-phase Fourier transform, which serves as a comprehensive sampling theorem applicable to a broad range of integral transforms. Shah and Tantary [20] formulated the sampling theorem for the QPFT and developed a novel convolution structure for efficient filtering in the quadratic-phase Fourier domain and also gave the advantages of the proposed convolution structure and its integration with the Wigner distribution to filter out undesired signal components. As a generalization of FT, QPFT can represent adaptively signals in both time and FT domains. Therefore, QPFT not only breaks through the limitation of FT in time-Fourier domain analysis but also overcomes the limitation of FT in indicating the signal’s characteristics. QPFT successfully inherits the advantages of MRA for FT. The MRA and the construction of orthogonal wavelets associated with QPFT are crucial in its perspective applications. Thus, detecting the MRA and the construction of the orthogonal wavelets related to QPFT is necessary. Therefore, our primary concern is introducing the notion of quadratic phase MRA, which allows a smoother construction of orthonormal bases simply and insightfully.
The main contributions of this article are as follows:
  • To give an alternative proof of Shannon’s sampling theorem associated with the quadratic phase Fourier transform.
  • Inspired by the sampling theorem of band-limited signals in the QPFT domain, the MRA associated with quadratic phase wavelet transform is developed.
  • Discuss the construction of the orthonormal basis of L 2 ( R ) starting from a given scaling function.
  • To give examples of quadratic phase wavelets from given scaling functions.
The rest of the article follows this structure: in Section 2, we offer a comprehensive introduction to the basics, covering the QPFT and also obtain some of its fundamental properties that are new in the literature. Moving on to Section 3, we give an alternative proof of the sampling theorem for the band-limited theorem in the QPFT domain. Based on this sampling theorem, we define a novel MRA and discuss constructing an orthonormal basis for L 2 ( R ) , followed by some examples in Section 4. Finally, in Section 5, we conclude our paper.

2. Preliminaries

In time-frequency analysis, a recent signal processing tool that has garnered attention is the quadratic-phase Fourier transform (QPFT), introduced by Castro et al. [21]. This transformative tool offers a unified approach to handling transient and non-transient signals.
Definition 1.
Given a parameter set Λ = ( A , B , C , D , E ) , the QPFT of f L 2 ( R ) is denoted as Q Λ f and is defined by
( Q Λ f ) ( ξ ) = 1 2 π R f ( t ) K Λ ( t , ξ ) d t ,
where K Λ ( t , ξ ) is a quadratic-phase kernel and is given by
K Λ ( t , ξ ) = e i ( A t 2 B t ξ + C ξ 2 + D t + E ξ )
and the corresponding inversion formula and Parseval’s formula for the QPFT reads
f ( t ) = R ( Q Λ f ) ( ξ ) K Λ ( t , ξ ) ¯ d ξ ,
and
f , g = | B | Q Λ f , Q Λ g f , g L 2 ( R ) ,
where A , B , C , D , E R , B 0 .

Some Properties Associated with the Quadratic Phase Fourier Transform

The lemma below gives some formula for the QPFT, that will be used later in this paper.
Lemma 1.
Let f L 2 ( R ) , a 0 , then the following holds:
(1) 
( Q Λ { f ( a t ) } ) ( ξ ) = 1 | a | e i C ξ 2 ξ a 2 + E ξ ξ a Q Λ e i A t a 2 A t 2 + D t a D t f ( t ) ξ a .
(2) 
( Q Λ { e i A a t 2 A t 2 + D a t D t f ( a t ) } ) ( ξ ) = 1 | a | e i C ξ 2 ξ a 2 + E ξ ξ a Q Λ f ( t ) ξ a .
Proof. 
Using the definition of the quadratic phase Fourier transform, we have
( Q Λ { f ( a t ) } ) ( ξ ) = 1 2 π R f ( a t ) e i ( A t 2 B t ξ + C ξ 2 + D t + E ξ ) d t = 1 2 π | a | R f ( t ) e i { A t a 2 B a t ξ + C ξ 2 + D a t + E ξ } d t = 1 2 π | a | R e i A t a 2 A t 2 + D t a D t f ( t ) e i A t 2 B t ξ a + C ξ a 2 + D t + E ξ a × e i C ξ 2 ξ a 2 + E ξ ξ a d t = 1 | a | e i C ξ 2 ξ a 2 + E ξ ξ a R e i A t a 2 A t 2 + D t a D t f ( t ) × 1 2 π e i A t 2 B t ξ a + C ξ a 2 + D t + E ξ a d t = 1 | a | e i C ξ 2 ξ a 2 + E ξ ξ a Q Λ e i A t a 2 A t 2 + D t a D t f ( t ) ξ a ,
i.e.,
( Q Λ { f ( a t ) } ) ( ξ ) = 1 | a | e i C ξ 2 ξ a 2 + E ξ ξ a Q Λ e i A t a 2 A t 2 + D T a D t f ( t ) ξ a .
Now,
( Q Λ { ϕ ( t ) f ( a t ) } ) ( ξ ) = 1 2 π R ϕ ( t ) f ( a t ) e i ( A t 2 B t ξ + C ξ 2 + D t + E ξ ) d t = 1 2 π | a | R ϕ t a f ( t ) e i A t a 2 B t a ξ + C ξ 2 + D t a + E ξ d t = 1 | a | e i C ξ 2 ξ a 2 + E ξ ξ a Q Λ e i A t a 2 A t 2 + D t a D t ϕ t a f ( t ) ξ a .
If we take ϕ t = e i A t 2 A a t 2 + D t D a t , then
( Q Λ { e i A t 2 A a t 2 + D t D a t } ) ( ξ ) = 1 | a | e i C ξ 2 ξ a 2 + E ξ ξ a Q Λ f ( t ) ξ a ,
i.e.,
( Q Λ { e i A a t 2 A t 2 + D a t D t f ( a t ) } ) ( ξ ) = 1 | a | e i C ξ 2 ξ a 2 + E ξ ξ a Q Λ f ( t ) ξ a .
This finishes the proof. □
The following lemma is an important tool in proving the Shannon’s sampling theorem in the QPFT domain. It says that if a function f is band-limited in the quadratic phase Fourier transform domain, then there is a function g depending on f such that it is band-limited in the Fourier domain. Note that, from here on, we take the value B in Λ as positive.
Lemma 2.
Assume a signal f ( t ) is band limited to Ω Λ in quadratic phase Fourier domain with parameter Λ = ( A , B , C , D , E ) and B > 0 . Let
g ( t ) = ( Q Λ f ) ( ξ ) e i [ B ξ t + C ξ 2 + E ξ ] d ξ ,
then g ( t ) is a signal band-limited in ( B Ω Λ , B Ω Λ ) in the Fourier domain.
Proof. 
Given,
g ( t ) = ( Q Λ f ) ( ξ ) e i [ B ξ t + C ξ 2 + E ξ ] d ξ .
Taking the Fourier transform, we have
F { g ( t ) } ( ω ) = ( Q Λ f ) ( ξ ) e i [ B ξ t + C ξ 2 + E ξ ] d ξ e i t ω 2 π d t = ( Q Λ f ) ( ξ ) e i [ C ξ 2 + E ξ ] e i [ B ξ t ] e i t ω 2 π d t d ξ = ( Q Λ f ) ( ξ ) e i [ C ξ 2 + E ξ ] 1 2 π e i ( B ξ ω ) t d t d ξ = ( Q Λ f ) ( ξ ) e i C ξ 2 + E ξ 2 π δ ( B ξ ω ) d ξ = ( Q Λ f ) ( ξ ) e i ( C ξ 2 + E ξ ) 2 π δ ( B ξ ω ) d ξ = 2 π B ( Q Λ f ) ξ B e i C ξ 2 B 2 + E ξ B 2 π δ ( ξ ω ) d ξ = 2 π B ( Q Λ f ) ω B e i C ω 2 B 2 + E ω B , i . e . , F { g ( t ) } ( ω ) = 2 π B ( Q Λ f ) ω B e i ( C ω 2 B 2 + E ω B ) .
Since ( Q Λ f ) ω is band-limited to Ω Λ , so s u p p ( Q Λ f ) [ Ω Λ , Ω Λ ] , i.e.,
( Q Λ f ) ω 0 a . e . , | ω | > Ω Λ ( Q Λ f ) ω B 0 a . e . , ω B > Ω Λ ( Q Λ f ) ω B 0 a . e . , | ω | > B Ω Λ ,
i.e., supp ( F { g ( t ) } ) [ B Ω Λ , B Ω Λ ] Thus, ( F { g ( t ) } ) ( ω ) 0 a . e . , | ω | > B Ω Λ . Therefore, g ( t ) is band-limited to B Ω Λ in the Fourier domain. □

3. Sampling Theorem for Band Limited Signal in QPFT Domain

The sampling theorem, also known as the Nyquist–Shannon sampling theorem, is a fundamental principle in signal processing and digital signal theory. It provides guidance on how to accurately reconstruct a continuous-time analog signal from its discrete samples. The theorem states that, to avoid aliasing and to perfectly reconstruct the original signal, the sampling rate (i.e., the number of samples taken per second) must be at least twice the highest frequency in the analog signal. Mathematically, if a band-limited signal contains a range of k frequencies, it can be accurately reconstructed by taking 2 k evenly spaced samples. Taking additional samples would prove redundant, whereas fewer samples would lead to a loss in signal quality. The sampling theorem can be expressed as follows: If a continuous-time signal is band-limited, meaning it contains no frequencies higher than a certain maximum frequency (known as the Nyquist frequency), then the signal can be completely reconstructed from its samples if the sampling rate is greater than or equal to twice the Nyquist frequency.
Inspired by the sampling theorem of band-limited signal in QPFT domain, in this section, the MRA associated with QPFT is established in the next section. The sampling theorem of a band-limited signal associated with QPFT is given by the following theorem.
Theorem 1.
Let signal f ( t ) be band-limited to Ω Λ in QPFT-domain having parameter Λ = ( A , B , C , D , E ) , B > 0 . Then, the following sampling theorem expansion for f ( t ) holds:
f ( t ) = e i ( A t 2 + D t ) n Z f ( n T ) e i A ( n T ) 2 + D ( n T ) s i n c B Ω Λ ( t n T ) π ,
where T is the sampling period and satisfies T = π B Ω Λ and is called as the Nyquist rate of sampling theorem associated with the quadratic phase Fourier transform.
Proof. 
We have
f ( t ) = 1 2 π R ( Q Λ f ) ξ e i ( A t 2 B t ξ + C ξ 2 + D t + E ξ ) d ξ = 1 2 π R ( Q Λ f ) ξ e i ( B t ξ + C ξ 2 + E ξ ) e i ( A t 2 + D t ) d ξ = R ( Q Λ f ) ξ e i ( B t ξ + C ξ 2 + E ξ ) d ξ 1 2 π e i ( A t 2 + D t ) = 1 2 π e i ( A t 2 + D t ) g ( t ) .
Since g ( t ) is band-limited to ( B Ω Λ , B Ω Λ ) in the Fourier domain, by applying the classical Shannon’s sampling theorem, we get
g ( t ) = n Z g ( n T ) s i n c B Ω Λ ( t n T ) π ,
where T = π B Ω Λ is the sample period. Therefore,
f ( t ) = 1 2 π e i ( A t 2 + D t ) n Z g ( n T ) s i n c B Ω Λ ( t n T ) π = e i ( A t 2 + D t ) n Z e i { A ( n t ) 2 + D ( n t ) } f ( n T ) s i n c B Ω Λ ( t n T ) π ,
i.e.,
f ( t ) = e i ( A t 2 + D t ) n Z f ( n T ) e i { A ( n t ) 2 + D ( n t ) } s i n c B Ω Λ ( t n T ) π .
This finishes the proof. □

4. Multiresolution Analysis Associated with QPFT

This section is devoted to the MRA associated with the QPFT. To introduce the definition, we first start with the following discussion, which has mainly to do with the Shannon’s sampling theorem discussed before. It motivated us to define an MRA associated with the QPFT. In what follows, the results also show the existence of the so-developed MRA.
When Ω Λ = π B , the set of band-limited signals in QPFT domain is denoted by V 0 Λ , i.e.,
V 0 Λ = f L 2 ( R ) : ( Q Λ f ) u = 0 , | u | Ω Λ = π B ,
where sampling period T = 1 . Therefore, from the sampling theorem, for all f ( t ) V 0 Λ , we get
f ( t ) = n Z e i ( A t 2 + D t ) f ( n ) e i { A n 2 + D n } s i n c B Ω Λ ( t n ) π .
Since Ω Λ = π B , we have
f ( t ) = n Z e i { A ( t 2 n 2 ) + D ( t n ) } f ( n ) s i n c ( t n ) ,
i.e.,
f ( t ) = n Z f ( n ) ϕ Λ , 0 , n ( t ) ,
where ϕ Λ , 0 , n ( t ) = e i { A ( t 2 n 2 ) + D ( t n ) } s i n c ( t n ) .
Combining with the orthogonality of { ϕ Λ , 0 , n } n Z , we can further obtain that { ϕ Λ , 0 , n = e i { A ( t 2 n 2 ) + D ( t n ) } s i n c ( t n ) } n Z forms an orthonormal basis of V 0 Λ .
When Ω Λ = 2 π B , T = 1 2 , the set of band-limited signal in the QPFT domain is denoted by V 1 Λ , i.e.,
V 1 Λ = f L 2 ( R ) : ( Q Λ f ) u = 0 , | u | Ω Λ = 2 π B .
Therefore, according to sampling theorem Equation (3),
g V 0 Λ ( Q Λ f ) ( u ) = 0 , | u | π B ( Q Λ f ) ( u ) = 0 , | u | 2 π B g V 1 Λ ,
i.e., V 0 Λ V 1 Λ . Also for T = 1 2 , Ω Λ = 2 π B , f V 1 Λ can be written as
f ( t ) = e i ( A t 2 + D t ) n Z f n 2 e i A n 2 2 D n 2 s i n c 2 t n 2 = n Z f n 2 e i A t 2 n 2 2 + D t n 2 s i n c 2 t n = n Z 1 2 f n 2 ϕ Λ , 1 , n ( t ) ,
where ϕ Λ , 1 , n ( t ) = 2 1 2 e i A t 2 n 2 2 + D t n 2 s i n c ( 2 t n ) .
It can be further obtained that { ϕ Λ , 1 , n = 2 1 2 e i A t 2 n 2 2 + D t n 2 s i n c ( 2 t n ) } n Z forms an orthonormal basis of V 1 Λ . For a = 2 , Equation (1) can be written as
( Q Λ { e i A 2 t 2 A t 2 + D 2 t D t f ( 2 t ) } ) ( ξ ) = 1 2 e i C ξ 2 ξ 2 2 + E ξ ξ 2 Q Λ f ξ 2 .
This implies that f V 0 Λ if and only if e i A 2 t 2 A t 2 + D 2 t D t f ( 2 t ) V 1 Λ . This is because
f V 0 Λ ( Q Λ f ) ( ξ ) = 0 , | ξ | Ω Λ = π B ( Q Λ f ) ξ 2 = 0 , ξ 2 π B ( Q Λ f ) ξ 2 = 0 , ξ 2 π B 1 2 e i C ξ 2 ξ 2 2 + E ξ ξ 2 Q Λ f ( t ) ξ 2 = 0 , ξ 2 π B ( Q Λ { e i A 2 t 2 A t 2 + D 2 t D t f ( 2 t ) } ) ( ξ ) = 0 , ξ 2 π B e i A 2 t 2 A t 2 + D 2 t D t f ( 2 t ) V 1 Λ ,
i.e., f V 0 Λ e i A 2 t 2 A t 2 + D 2 t D t f ( 2 t ) V 1 Λ .
Generally, let Ω Λ = 2 k π B , T = 1 2 k ,
V k Λ = { f L 2 ( R ) : ( Q Λ f ) ( ξ ) = 0 , | u | 2 k π B } .
Now, f V k Λ , we have
f ( t ) = n Z f n 2 k e i A t 2 n 2 k 2 + D t n 2 k s i n c 2 k t n = n Z 1 2 k 2 f n 2 k ϕ Λ , k , n ( t ) ,
where ϕ Λ , k , n ( t ) = 2 k 2 e i A t 2 n 2 k 2 + D t n 2 k s i n c 2 k t n . Thus, { ϕ Λ , k , n ( t ) } n Z forms an orthonormal basis of V k Λ .
Thus, k Z , we have
(a)
V k Λ V k + 1 Λ ;
(b)
f ( t ) V k Λ e i A 2 t 2 A t 2 + D 2 t D t f ( 2 t ) V k + 1 Λ ;
(c)
k Z V k Λ = { 0 } and k Z V k Λ ¯ = L 2 ( R ) .
To put it briefly, the sampling theorem for band-limited signals in the QPFT domain serves as the inspiration to establish an orthonormal MRA associated with QPFT.
In this section, our focus lies on introducing the concept of a quadratic phase MRA within the space L 2 ( R ) . This MRA will hold significant importance in developing the quadratic phase orthonormal wavelet basis for L 2 ( R ) . Initially, we present the formal definition of a special affine MRA in L 2 ( R ) .
Definition 2.
An orthogonal MRA associated with QPFT is defined as a sequence of closed subspace V k Λ , k Z such that
(A) 
V k Λ V k + 1 Λ , k Z .
(B) 
f ( t ) V k Λ e i A 2 t 2 A t 2 + D 2 t D t f ( 2 t ) V k + 1 Λ k Z .
(C) 
k Z V k Λ = { 0 } and k Z V k Λ ¯ = L 2 ( R ) .
(D) 
There exist a function ϕ L 2 ( R ) such that { ϕ Λ , 0 , n ( t ) = e i { A ( t 2 n 2 ) + D ( t n ) } s i n c ( t n ) } n Z is an orthonormal basis of the subspace V 0 Λ , where ϕ is called the scaling function of the given MRA.
Lemma 3.
The family { ϕ Λ , 0 , n : n Z } , given by the above, constitute an orthonormal system in L 2 ( R ) iff
k Z ( F ϕ ) ( ξ + 2 k π ) 2 = 1 2 π .
Proof. 
We have
( Q Λ { ϕ Λ , 0 , n } ) ( ξ ) = R ϕ ( t n ) e i { A ( t 2 n 2 ) + D ( t n ) } 1 2 π e i ( A t 2 B t ξ + C ξ 2 + D t + E ξ ) d t = R ϕ ( t n ) e i { A ( n 2 ) + D ( n ) } 1 2 π e i ( B t ξ + C ξ 2 + E ξ ) d t = 1 2 π e i { A n 2 + D n B n ξ + C ξ 2 + E ξ } R ϕ ( t ) e i t ( B ξ ) d t ,
i.e., ( Q Λ { ϕ Λ , 0 , n } ) ( ξ ) = e i { A n 2 + D n B n ξ + C ξ 2 + E ξ } ( F ϕ ) ( B ξ ) . Now since,
ϕ Λ , 0 , n , ϕ Λ , 0 , l = B Q Λ { ϕ Λ , 0 , n } , Q Λ { ϕ Λ , 0 , l } = B R ( Q Λ { ϕ Λ , 0 , n } ) ( ξ ) ( Q Λ { ϕ Λ , 0 , l } ) ( ξ ) ¯ d ξ = B R e i ( A n 2 + D n B n ξ + C ξ 2 + E ξ ) e i ( A l 2 + D l B l ξ + C ξ 2 + E ξ ) ( F ϕ ) ( B ξ ) 2 d ξ = B R e i ( A ( n 2 l 2 ) + D ( n l ) B ( n l ) ξ ) ( F ϕ ) ( B ξ ) 2 d ξ = B e i ( A ( n 2 l 2 ) + D ( n l ) R e i B ( n l ) ξ ) ( F ϕ ) ( B ξ ) 2 d ξ = e i ( A ( n 2 l 2 ) + D ( n l ) R e i ( n l ) ξ ) ( F ϕ ) ( ξ ) 2 d ξ .
Since { ϕ Λ , 0 , n } n Z forms an orthonormal basis of V 0 Λ , we have
ϕ Λ , 0 , n , ϕ Λ , 0 , l = δ n , l k , l R .
This implies e i ( A ( n 2 l 2 ) + D ( n l ) R e i ( n l ) ξ ) ( F ϕ ) ( ξ ) 2 d ξ = δ n , l . Set n l = q then
e i [ A ( q + l ) 2 l 2 ] + D q R e i q ξ ( F ϕ ) ( ξ ) 2 d ξ = δ q + l , l R e i q ξ ( F ϕ ) ( ξ ) 2 d ξ = δ q , 0 , q Z k Z 2 k π 2 ( k + 1 ) π e i q ξ ( F ϕ ) ( ξ ) 2 d ξ = δ q , 0 k Z 0 2 π e i q ( η + 2 k π ) ( F ϕ ) ( η + 2 k π ) 2 d η = δ q , 0 k Z 0 2 π e i q η ( F ϕ ) ( η + 2 k π ) 2 d η = δ q , 0 , i . e . , 0 2 π e i q η k Z ( F ϕ ) ( η + 2 k π ) 2 d η = δ q , 0 .
Let F ( u ) = k Z ( F ϕ ) ( η + 2 k π ) 2 . Then
F ( u + 2 π ) = u Z ( F ϕ ) ( u + 2 π + 2 k π ) 2 = u Z ( F ϕ ) ( u + 2 ( k + 1 ) π ) 2 = r Z ( F ϕ ) ( u + 2 r π ) 2 = F ( u ) ,
i.e., F ( u ) is 2 π periodic. Therefore,
0 2 π e i q η F ( η ) d η = δ q , 0 1 2 π 0 2 π e i q η F ( η ) d η = 1 2 π δ q , 0 F ( η ) = 1 2 π ,
i.e., k Z ( F ϕ ) ( ξ + 2 k π ) 2 = 1 2 π .
Conversely, let k Z ( F ϕ ) ( ξ + 2 k π ) 2 = 1 2 π , then
ϕ Λ , 0 , n , ϕ Λ , 0 , l = e i ( A ( n 2 l 2 ) + D ( n l ) R e i ( n l ) ξ ) ( F ϕ ) ( ξ ) 2 d ξ = e i ( A ( n 2 l 2 ) + D ( n l ) k Z 2 k π 2 ( k + 1 ) π e i ( n l ) ξ ( F ϕ ) ( ξ ) 2 d ξ = e i ( A ( n 2 l 2 ) + D ( n l ) k Z 0 2 π e i ( n l ) ξ ( F ϕ ) ( ξ + 2 k π ) 2 d ξ = e i ( A ( n 2 l 2 ) + D ( n l ) 0 2 π e i ( n l ) ξ k Z ( F ϕ ) ( ξ + 2 k π ) 2 d ξ = e i ( A ( n 2 l 2 ) + D ( n l ) 1 2 π 0 2 π e i ( n l ) ξ d ξ = 1 , n = l 0 , n l = δ n , l ,
i.e., the system is orthonormal. Hence the conclusion follows. □
Let { V k Λ } k Z be an orthonormal MRA of L 2 ( R ) . Since ϕ Λ , 0 , 0 ( t ) V 0 Λ V 1 Λ , and { ϕ Λ , 1 , n } n Z forms an orthonormal basis of V 1 Λ , so there exists { h n } n Z such that
ϕ Λ , 0 , 0 ( t ) = n Z h n ϕ Λ , 1 , n ( t ) = n Z h n 2 ϕ ( 2 t n ) e i A t 2 n 2 2 + D ( t n 2 ) .
Equation (4) is called the quadratic phase refinement equation. Here,
h n = 2 R ϕ Λ , 0 , 0 ( t ) ϕ ( 2 t n ) ¯ e i A t 2 n 2 2 + D ( t n 2 ) d t = 2 R ϕ ( t ) ϕ ( 2 t n ) ¯ e i A n 2 2 + D n 2 d t .
Now, from Equation (4), we have
ϕ ( t ) e i ( A t 2 + D t ) = 2 n Z h n ϕ ( 2 t n ) e i A t 2 n 2 2 + D ( t n 2 ) .
Taking QPFT on both sides we get
R ϕ ( t ) e i ( A t 2 + D t ) e i ( A t 2 B t ξ + C ξ 2 + D t + E ξ ) d t = 2 R n Z h n ϕ ( 2 t n ) e i A t 2 n 2 2 + D ( t n 2 ) e i ( A t 2 B t ξ + C ξ 2 + D t + E ξ ) d t .
This gives
R ϕ ( t ) e i ( B t ξ + C ξ 2 + E ξ ) d t = 2 n Z h n R ϕ ( 2 t n ) e i A n 2 2 D ( n 2 ) e i ( B t ξ + C ξ 2 + E ξ ) d t i . e . , e i ( C ξ 2 + E ξ ) R ϕ ( t ) e i B t ξ d t = e i ( C ξ 2 + E ξ ) 2 n Z h n R ϕ ( 2 t n ) e i A n 2 2 + D ( n 2 ) e i B t ξ d t .
Thus,
R ϕ ( t ) e i t ( B ξ ) d t = 2 n Z h n e i A n 2 2 + D ( n 2 ) R ϕ ( 2 t n ) e i t ( B ξ ) d t = n Z h n e i A n 2 2 + D ( n 2 ) R 2 ϕ ( 2 t n ) e i t ( B ξ ) d t = n Z h n e i A n 2 2 + D ( n 2 ) R 1 2 ϕ ( t n ) e i t B ξ 2 d t = n Z h n e i A n 2 2 + D ( n 2 ) R 1 2 ϕ ( t ) e i ( t + n ) B ξ 2 d t = n Z 1 2 h n e i A n 2 2 + D ( n 2 ) n B ξ 2 R ϕ ( t ) e i t B ξ 2 d t ,
i.e., ( F ϕ ) ( B ξ ) = Λ 0 B ξ 2 ( F ϕ ) B ξ 2 , where Λ 0 B ξ 2 = 1 2 n Z h n e i A n 2 2 + D ( n 2 ) n B ξ 2 .
Equivalently,
( F ϕ ) ( ξ ) = Λ 0 ξ 2 ( F ϕ ) ξ 2 ,
where Λ 0 ξ = 1 2 n Z h n e i A n 2 2 + D ( n 2 ) n ξ .
Since { ϕ Λ , 0 , n } n Z is an orthonormal basis or orthonormal system of V 0 Λ , so by Lemma 3, we have
k Z ( F ϕ ) ( ξ + 2 k π ) 2 = 1 2 π .
Hence, we can write
1 2 π = k Z ( F ϕ ) ( ξ + 2 k π ) 2 = k Z Λ 0 ξ + 2 k π 2 ( F ϕ ) ξ + 2 k π 2 2 = k Z Λ 0 ξ + 2 k π 2 2 ( F ϕ ) ξ + 2 k π 2 2 = k Z Λ 0 ξ 2 + 2 k π 2 ( F ϕ ) ξ 2 + 2 k π 2 + k Z Λ 0 ξ 2 + ( 2 k + 1 ) π 2 ( F ϕ ) ξ 2 + ( 2 k + 1 ) π 2 .
Observe that
Λ 0 ξ = 1 2 n Z h n e i A n 2 2 + D ( n 2 ) n ξ Λ 0 ξ + 2 π = 1 2 n Z h n e i A n 2 2 + D ( n 2 ) n ( ξ + 2 π ) = 1 2 n Z h n e i A n 2 2 + D ( n 2 ) n ξ e 2 i n π = 1 2 n Z h n e i A n 2 2 + D ( n 2 ) n ξ i . e . , Λ 0 ( ξ + 2 π ) = Λ 0 ( ξ ) .
Also,
Λ 0 ( ξ + 2 k π ) = Λ 0 ( ξ ) .
Hence, from (7)
1 2 π = k Z Λ 0 ξ 2 + 2 k π 2 ( F ϕ ) ξ 2 + 2 k π 2 + k Z Λ 0 ξ 2 + ( 2 k + 1 ) π 2 ( F ϕ ) ξ 2 + ( 2 k + 1 ) π 2 = Λ 0 ξ 2 2 k Z ( F ϕ ) ξ 2 + 2 k π 2 + Λ 0 ξ 2 + π 2 k Z ( F ϕ ) ξ 2 + ( 2 k + 1 ) π 2 .
Now using Equation (6), we have
Λ 0 ξ 2 2 · 1 2 π + Λ 0 ξ 2 + π 2 · 1 2 π = 1 2 π Λ 0 ξ 2 2 + Λ 0 ξ 2 + π 2 = 1 .
Given an orthogonal MRA { V n Λ } n Z , we define another sequence { W n Λ } n Z , of closed subspaces of L 2 ( R ) by V j + 1 Λ = V j Λ W j Λ , j Z . Followed by a definition, these subspace inherit the scaling property of { V j Λ } j Z , namely
f ( t ) W j Λ f ( 2 t ) e i [ A { ( 2 t ) 2 t 2 } + D ( 2 t t ) ] W j + 1 Λ , j Z .
Moreover, the subspaces W j Λ are mutually orthogonal with the following decomposition formula
L 2 ( R ) = j Z W j Λ .
Note that condition (10) means that any orthonormal basis for L 2 ( R ) can be constructed by finding out an orthonormal basis for the subspace W j Λ . On the other hand, condition (9) implies that the quadratic phase basis can be constructed as long as the orthonormal basis for W 0 Λ is found. Therefore, our main concern is to construct a mother function ψ Λ , 0 , 0 in W 0 Λ such that { ψ Λ , 0 , k : k Z } forms an orthonormal basis of W 0 Λ .
Suppose ψ Λ , 0 , 0 ( t ) W 0 Λ V 1 Λ , there exists { d k } k Z such that
ψ Λ , 0 , 0 ( t ) = 2 k Z d k ϕ ( 2 t k ) e i A t 2 k 2 2 + D t k 2 .
Equation (11) is called quadratic phase wavelet equation. Taking QPFT on both sides, we get
( F ψ ) ( B ξ ) = Λ 1 B ξ 2 ( F ϕ ) B ξ 2 , i . e . , ( F ψ ) ( ξ ) = Λ 1 ξ 2 ( F ϕ ) ξ 2 ,
where
Λ 1 ( ξ ) = 1 2 n Z d n e i A n 2 2 + D n 2 n ξ .
Since W 0 Λ and V 0 Λ , are orthogonal in V 1 Λ , we have
0 = ϕ Λ , 0 , k , ϕ Λ , 0 , l = B Q Λ { ϕ Λ , 0 , k } , Q Λ { ϕ Λ , 0 , l } .
Observe that
( Q Λ { ϕ Λ , 0 , k } ) ( ξ ) = 1 2 π R ϕ Λ , 0 , k ( t ) e i ( A t 2 B t ξ + C ξ 2 + D t + E ξ ) d t = 1 2 π R ϕ ( t k ) e i { A k 2 + D k } e i ( B t ξ + C ξ 2 + E ξ ) d t = 1 2 π e i { A k 2 + D k } e i ( C ξ 2 + E ξ ) R ϕ ( t k ) e i t ( B ξ ) d t = 1 2 π e i { A k 2 + D k } e i ( C ξ 2 + E ξ ) R ϕ ( t ) e i ( t + k ) B ξ d t = e i { A k 2 + D k } e i ( C ξ 2 + E ξ ) e i k ( B ξ ) 1 2 π R ϕ ( t ) e i t ( B ξ ) d t i . e . , Q Λ ( ϕ Λ , 0 , k ) ( ξ ) = e i { A k 2 + D k + C ξ 2 + E ξ k B ξ } ( F ϕ ) ( B ξ ) .
Similarly,
Q Λ ( ψ Λ , 0 , l ) ( ξ ) = e i { A l 2 + D l + C ξ 2 + E ξ l B ξ } ( F ψ ) ( B ξ ) .
Therefore, using Equations (15) and (16) in (14), we get
0 = R e i { A k 2 + D k + C ξ 2 + E ξ k B ξ } ( F ϕ ) ( B ξ ) e i { A l + D l + C ξ 2 + E ξ l B ξ } ( F ψ ) ( B ξ ) ¯ d ξ = e i { A ( k 2 l 2 ) + D ( k l ) } R e i ( k l ) B ξ ( F ϕ ) ( B ξ ) ( F ψ ) ( B ξ ) ¯ d ξ = 1 B e i { A ( k 2 l 2 ) + D ( k l ) } R e i ( l k ) ξ ( F ϕ ) ( ξ ) ( F ψ ) ( ξ ) ¯ d ξ = 1 B e i { A ( k 2 l 2 ) + D ( k l ) } R e i ( l k ) ξ Λ 0 ξ 2 Λ 1 ξ 2 ¯ ( F ϕ ) ξ 2 2 d ξ = 1 B e i { A ( k 2 l 2 ) + D ( k l ) } m Z 4 m π 4 ( m + 1 ) π e i ( l k ) ξ Λ 0 ξ 2 Λ 1 ξ 2 ¯ ( F ϕ ) ξ 2 2 d ξ = 1 B e i { A ( k 2 l 2 ) + D ( k l ) } m Z 0 4 π e i ( l k ) ( ξ + 4 m π ) Λ 0 ξ 2 + 2 m π Λ 1 ξ 2 + 2 m π ¯ × ( F ϕ ) ξ 2 + 2 m π 2 d ξ = 1 B e i { A ( k 2 l 2 ) + D ( k l ) } m Z 0 4 π e i ( l k ) ξ Λ 0 ξ 2 Λ 1 ξ 2 ¯ ( F ϕ ) ξ 2 + 2 m π 2 d ξ = 1 B e i { A ( k 2 l 2 ) + D ( k l ) } m Z [ 0 2 π e i ( l k ) ξ Λ 0 ξ 2 Λ 1 ξ 2 ¯ ( F ϕ ) ξ 2 + 2 m π 2 d ξ + 2 π 4 π e i ( l k ) ξ Λ 0 ξ 2 Λ 1 ξ 2 ¯ ( F ϕ ) ξ 2 + 2 m π 2 d ξ ] = 1 B e i { A ( k 2 l 2 ) + D ( k l ) } m Z [ 0 2 π e i ( l k ) ξ Λ 0 ξ 2 Λ 1 ξ 2 ¯ ( F ϕ ) ξ 2 + 2 m π 2 d ξ + 0 2 π e i ( l k ) ( ξ + 2 π ) Λ 0 ξ + 2 π 2 Λ 1 ξ + 2 π 2 ¯ ( F ϕ ) ξ + 2 π 2 + 2 m π 2 d ξ ] = 1 B e i { A ( k 2 l 2 ) + D ( k l ) } [ 0 2 π e i ( l k ) ξ Λ 0 ξ 2 Λ 1 ξ 2 ¯ m Z ( F ϕ ) ξ 2 + 2 m π 2 d ξ + 0 2 π e i ( l k ) ( ξ + 2 π ) Λ 0 ξ 2 + π Λ 1 ξ 2 + π ¯ m Z ( F ϕ ) ξ 2 + ( 2 m + 1 ) π 2 d ξ ] = 1 B e i { A ( k 2 l 2 ) + D ( k l ) } [ 0 2 π e i ( l k ) ξ Λ 0 ξ 2 Λ 1 ξ 2 ¯ d ξ + 0 2 π e i ( l k ) ( ξ + 2 π ) Λ 0 ξ 2 + π Λ 1 ξ 2 + π ¯ d ξ ] = 1 B e i { A ( k 2 l 2 ) + D ( k l ) } 0 2 π e i ( l k ) ξ Λ 0 ξ 2 Λ 1 ξ 2 ¯ + Λ 0 ξ 2 + π Λ 1 ξ 2 + π ¯ d ξ .
Therefore,
0 2 π e i ( l k ) ξ Λ 0 ξ 2 Λ 1 ξ 2 ¯ + Λ 0 ξ 2 + π Λ 1 ξ 2 + π ¯ d ξ = 0 .
From Equation (17), we conclude that
Λ 0 ξ 2 Λ 1 ξ 2 ¯ + Λ 0 ξ 2 + π Λ 1 ξ 2 + π ¯ = 0 .
This implies
Λ 0 ξ Λ 1 ξ ¯ + Λ 0 ξ + π Λ 1 ξ + π ¯ = 0 .
Equation (18) can be written in the matrix form as
M M * = I 2 × 2 ,
where M * denotes the conjugate transpose of M, I 2 × 2 is the identity matrix, and
M = Λ 0 ( ξ ) Λ 0 ( ξ + π ) Λ 1 ( ξ ) Λ 1 ( ξ + π ) .
Since Λ 0 ( ξ ) and Λ 0 ( ξ + π ) cannot vanish together on a set of non-zero measures due to the orthogonal property, there exists a 2 π - periodic function λ ( ξ ) such that
Λ 1 ( ξ ) , Λ 1 ( ξ + π ) = λ ( ξ ) Λ 0 ( ξ + π ) ¯ , λ ( ξ ) Λ 0 ( ξ ) ¯ .
Since,
Λ 1 ( ξ ) = λ ( ξ ) Λ 0 ( ξ + π ) ¯ Λ 1 ( ξ + π ) = λ ( ξ + π ) Λ 0 ( ξ + 2 π ) ¯ .
Using Equation (8),
Λ 1 ( ξ + π ) = λ ( ξ + π ) Λ 0 ( ξ ) ¯ .
Therefore, we have
λ ( ξ ) Λ 0 ( ξ ) ¯ = λ ( ξ + π ) Λ 0 ( ξ ) ¯ { λ ( ξ + π ) + λ ( ξ ) } Λ 0 ( ξ ) ¯ = 0 λ ( ξ + π ) + λ ( ξ ) = 0 .
Therefore, λ ( ξ ) is 2 π - periodic, it can be expressed as
λ ( ξ ) = k Z c k e i k ξ ,
where
c k = 1 2 π 0 2 π λ ( ξ ) e i k ξ d ξ = 1 2 π 0 π λ ( ξ ) e i k ξ d ξ + π 2 π λ ( ξ ) e i k ξ d ξ = 1 2 π 0 π λ ( ξ ) e i k ξ d ξ + 0 π λ ( ξ + π ) e i k ( ξ + π ) d ξ = 1 2 π 0 π λ ( ξ ) e i k ξ d ξ + 0 π ( 1 ) k λ ( ξ + π ) e i k ξ d ξ = 1 2 π 0 π λ ( ξ ) e i k ξ d ξ + ( 1 ) k 0 π λ ( ξ ) e i k ξ d ξ = 1 2 π { 1 ( 1 ) k } 0 π λ ( ξ ) e i k ξ d ξ .
Therefore, λ ( ξ ) = k Z c k e i k ξ , where
c k = 1 2 π { 1 ( 1 ) k } 0 π λ ( ξ ) e i k ξ d ξ c k = 0 k = 2 m , m Z .
Therefore,
λ ( ξ ) = l Z c 2 l + 1 e i ( 2 l + 1 ) ξ = l Z c 2 l + 1 e 2 i l ξ e i ξ = e i ξ l Z c 2 l + 1 e 2 i l ξ = e i ξ γ ( 2 ξ ) ,
where γ ( ξ ) = l Z c 2 l + 1 e 2 i l ξ . Now,
Λ 0 ( ξ ) = 1 2 n Z h n e i A n 2 2 + D n 2 n ξ Λ 0 ( ξ + π ) = 1 2 n Z h n e i A n 2 2 + D n 2 n ( ξ + π ) = 1 2 n Z h n e i A n 2 2 + D n 2 n ( ξ + π ) = 1 2 n Z h n e i A n 2 2 + D n 2 n ξ e i n ξ = 1 2 n Z ( 1 ) n h n e i A n 2 2 + D n 2 n ξ .
Therefore,
Λ 0 ( ξ + π ) ¯ = 1 2 n Z ( 1 ) n h n ¯ e i A n 2 2 + D n 2 n ξ .
Thus, from (20)
Λ 1 ( ξ ) = λ ( ξ ) Λ 0 ( ξ + π ) ¯ = e i ξ γ ( 2 ξ ) Λ 0 ( ξ + π ) ¯ .
In particular, for γ ( 2 ξ ) = 1 , using (13), we have
Λ 1 ( ξ ) = e i ξ γ ( 2 ξ ) Λ 0 ( ξ + π ) ¯ .
This implies
1 2 n Z d n e i A n 2 2 + D n 2 n ξ = e i ξ 1 2 n Z ( 1 ) n h n ¯ e i A n 2 2 + D n 2 n ξ i . e . , n Z d n e i A n 2 2 + D n 2 e i n ξ = n Z ( 1 ) n h n ¯ e i A n 2 2 + D n 2 e i ( n 1 ) ξ i . e . , n Z d n e i A n 2 2 + D n 2 e i n ξ e i k ξ = n Z ( 1 ) n h n ¯ e i A n 2 2 + D n 2 e i ( n 1 ) ξ e i k ξ .
Integrating both sides, we get
n Z d n e i A n 2 2 + D n 2 R e i ( n k ) ξ d ξ = n Z ( 1 ) n h n ¯ e i A n 2 2 + D n 2 R e i ( 1 n k ) ξ d ξ n Z d n e i A n 2 2 + D n 2 δ n k = n Z ( 1 ) n h n ¯ e i A n 2 2 + D n 2 δ 1 n k k Z d k e i A k 2 2 + D k 2 = ( 1 ) 1 k h 1 k ¯ e i A 1 k 2 2 + D 1 k 2 d k = ( 1 ) 1 k h 1 k ¯ e i A 1 k 2 2 + D 1 k 2 + A k 2 2 + D k 2 d k = ( 1 ) 1 k h 1 k ¯ e i A ( 1 k ) 2 + k 2 4 + D 2 .
Therefore, equivalently, we can write the wavelet coefficients d k of Equation (11) as
d k = ( 1 ) 1 k h 1 k ¯ e i A ( 1 k ) 2 + k 2 4 + D 2 , k Z .
The above discussion can be summarized in the following theorem.
Theorem 2.
If { V n Λ } n Z is the quadratic phase MRA associated with the scaling function ϕ , then there exists a function ψ such that
ψ Λ , 0 , 0 ( t ) = 2 k Z d k ϕ ( 2 t k ) e i A t 2 k 2 2 + D t k 2 ,
where d k is given by (21) with h k = 2 e i A k 2 2 + D k 2 R ϕ ( t ) ϕ ( 2 t k ) ¯ d t , i.e., the system ψ Λ , k , n , k , n Z is an orthonormal basis of L 2 ( R ) .
Example 1.
It is observed in the earlier discussion that the function s i n c ( t ) = sin ( π t ) π t is a scaling function for the quadratic phase MRA { V n Λ } n Z , where { ϕ Λ , 0 , n ( t ) = e i { A ( t 2 n 2 ) + D ( t n ) ) } s i n c ( t n ) } n Z is an orthonormal basis of the subspace V 0 Λ . Hence,
h n = 1 2 , n = 0 2 sin π n 2 e i A n 2 4 + D n 2 π n , n 0 ,
which results in
d n = 1 2 e i A 4 + D 2 , n = 1 2 ( 1 ) n cos π n 2 e 1 4 i n ( A n + 2 D ) π ( n 1 ) , n 1 .
Thus, the quadratic phase wavelet corresponding to the scaling function s i n c ( t ) is given by
ψ Λ , 0 , 0 ( t ) = 2 d 1 e i A t 2 1 2 2 + D t 1 2 ϕ ( 2 t 1 ) + 2 n Z , n 1 d n ϕ ( 2 t n ) e i A t 2 n 2 2 + D t n 2 = e i t ( A t + D ) ϕ ( 2 t 1 ) + 2 n Z , n 1 2 ( 1 ) n cos π n 2 e i t ( A t + D ) π ( n 1 ) ϕ ( 2 t n ) = e i t ( A t + D ) sin ( ( 2 t 1 ) π ) ( 2 t 1 ) π + 2 n Z , n 1 2 ( 1 ) n cos π n 2 e i t ( A t + D ) π ( n 1 ) sin ( ( 2 t n ) π ) ( 2 t n ) π .
The plots of the real and imaginary part of ψ Λ , 0 , 0 are given below for the particular choice of the parameter Λ = ( 0 , 1 , 0 , 0 , 0 ) and Λ = ( 1 3 , 1 , 0 , 1 5 , 0 )
Example 2.
Let ϕ ( t ) = χ [ 0 , 1 ) ( t ) , where χ [ 0 , 1 ) ( t ) is a characteristic function on [ 0 , 1 ) . It is a matter of simple verification that the set ϕ Λ , 0 , n ( t ) : n Z forms an orthonormal system. Hence it forms an orthonormal basis of the set V 0 Λ , thus is a scaling function associated with the MRA { V n Λ } n Z . Thus,
h n = 1 2 , n = 0 1 2 e i A 1 2 2 + D 1 2 , n = 1 0 , o t h e r w i s e
resulting in
d n = 1 2 , n = 0 1 2 e i A 1 2 2 + D 1 2 , n = 1 0 , o t h e r w i s e .
Thus, the quadratic phase wavelet corresponding to the scaling function χ [ 0 , 1 ) ( t ) is given by
ψ Λ , 0 , 0 ( t ) = e i t ( A t + D ) , 0 t < 1 2 e i t ( A t + D ) , 1 2 t < 1 .
The plots of the real and imaginary parts of ψ Λ , 0 , 0 are given below for the particular choice of the parameter Λ = ( 0 , 1 , 0 , 0 , 0 ) and Λ = ( 1 3 , 1 , 0 , 1 5 , 0 )
Remark 1.
By virtue of Lemma 3 and the Definition 2 of MRA we can say that any function that serves as a scaling function in the classical MRA will also serve as a scaling function for the MRA given by Definition 2. But, depending on the choice of parameters Λ, we can have different quadratic phase wavelets and thus different families of orthonormal bases of L 2 ( R ) . In particular, for the choice of the parameter Λ = ( 0 , 1 , 0 , 0 , 0 ) , we get the classical wavelets and the quadratic phase wavelets for Λ = ( 1 3 , 1 , 0 , 1 5 , 0 ) (see Figure 1, Figure 2, Figure 3 and Figure 4). The flexibility in the choice of the parameters results in the development of some novel families of orthonormal bases of L 2 ( R ) corresponding to the same scaling function.

5. Conclusions

The MRA and the construction of orthogonal wavelets for QPFT play a vital role in facilitating prospective applications of QPFT. In this paper, we gave an alternative proof of the Shannon’s sampling theorem applicable to the band-limited signal in the QPFT. Inspired by the theorem, we developed an MRA associated with QPFT. Subsequently, we discussed the construction of the quadratic phase wavelets for a given scaling function, followed by some of its examples.

Author Contributions

Conceptualization, B.G. and A.K.V.; methodology, B.G. and N.K.; software, B.G.; validation, B.G., A.K.V. and R.P.A.; formal analysis, B.G. and N.K.; investigation, B.G. and A.K.V.; resources, R.P.A.; writing—original draft preparation, B.G. and N.K.; writing—review and editing, B.G. and N.K.; visualization, A.K.V. and R.P.A.; supervision, A.K.V. and R.P.A.; project administration, A.K.V.; funding acquisition, B.G., A.K.V. and N.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work is partly supported by UGC File No. 16-9 (June 2017)/2018(NET/CSIR), DST SERB FILE NO. MTR/2021/000907 and CSIR File No. 09/1023(0035)/2020-EMR-I), New Delhi, India.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
MDPIMultidisciplinary Digital Publishing Institute
DOAJDirectory of open access journals
TLAThree letter acronym
LDLinear dichroism

References

  1. Debnath, L.; Shah, F.A. Wavelet Transforms and Their Applications; Birkhäuser: New York, NY, USA, 2015. [Google Scholar]
  2. Debnath, L.; Shah, F.A. Lecture Notes on Wavelet Transforms; Birkhäuser: Boston, MA, USA, 2017. [Google Scholar]
  3. Namias, V. The fractional order Fourier transform and its application to quantum mechanics. IMA J. Appl. Math. 1980, 25, 241–265. [Google Scholar] [CrossRef]
  4. Almeida, L.B. The fractional Fourier transform and time-frequency representations. IEEE Trans. Signal Process. 1994, 42, 3084–3091. [Google Scholar] [CrossRef]
  5. Xu, T.Z.; Li, B.Z. Linear Canonical Transform and Its Applications; Science Press: Beijing, China, 2013. [Google Scholar]
  6. Abe, S.; Sheridan, J.T. Generalization of the fractional Fourier transformation to an arbitrary linear lossless transformation an operator approach. J. Phys. Math. Gen. 1994, 27, 4179. [Google Scholar] [CrossRef]
  7. Castro, L.P.; Haque, M.R.; Murshed, M.M.; Saitoh, S.; Tuan, N.M. Quadratic Fourier transforms. Ann. Funct. Anal. 2014, 5, 10–23. [Google Scholar] [CrossRef]
  8. Shah, F.A.; Lone, W.Z.; Tantary, A.Y. Short-time quadratic-phase Fourier transform. Optik 2021, 245, 167689. [Google Scholar] [CrossRef]
  9. Shah, F.A.; Lone, W.Z. Quadratic-phase wavelet transform with applications to generalized differential equations. Math. Methods Appl. Sci. 2022, 45, 1153–1175. [Google Scholar] [CrossRef]
  10. Prasad, A.; Sharma, P.B. The quadratic-phase Fourier wavelet transform. Math. Methods Appl. Sci. 2020, 43, 1953–1969. [Google Scholar] [CrossRef]
  11. Mallat, S.G. Multiresolution approximations and wavelet orthonormal bases of L2(ℝ). Trans. Am. Math. Soc. 1989, 315, 69–87. [Google Scholar]
  12. Madych, W.R. Some elementary properties of multiresolution analyses of L2(ℝn). In Wavelets: A Tutorial in Theory and Applications; Chui, C.K., Ed.; Academic Press: Boston, MA, USA, 1992; pp. 259–277. [Google Scholar]
  13. Zhang, Z. Supports of Fourier transforms of scaling functions. Appl. Comput. Harmon. Anal. 2007, 22, 141–156. [Google Scholar] [CrossRef]
  14. Malhotra, H.K.; Vashisht, L.K. On scaling functions of non-uniform multiresolution analysis in L2(ℝ). Int. J. Wavelets Multiresolut. Inf. Process. 2020, 18, 1950055. [Google Scholar] [CrossRef]
  15. Shi, J.; Liu, X.; Zhang, N. Multiresolution analysis and orthogonal wavelets associated with fractional wavelet transform. Signal Image Video Process. 2015, 9, 211–220. [Google Scholar] [CrossRef]
  16. Ahmad, O.; Sheikh, N.A.; Shah, F.A. Fractional multiresolution analysis and associated scaling functions in L2(ℝ). Anal. Math. Phys. 2021, 11, 47. [Google Scholar] [CrossRef]
  17. Dai, H.; Zheng, Z.; Wang, W. A new fractional wavelet transform. Commun. Nonlinear Sci. Numer. Simul. 2017, 44, 19–36. [Google Scholar] [CrossRef]
  18. Shah, F.A.; Lone, W.Z. Special affine multiresolution analysis and the construction of orthonormal wavelets in L2(ℝ). Appl. Anal. 2023, 102, 2540–2566. [Google Scholar] [CrossRef]
  19. Lone, W.Z.; Shah, F.A. Shift-invariant spaces and dynamical sampling in quadratic-phase Fourier domains. Optik 2022, 260, 169063. [Google Scholar] [CrossRef]
  20. Shah, F.A.; Tantary, A.Y. Sampling and multiplicative filtering associated with the quadratic-phase Fourier transform. Signal Image Video Process. 2023, 17, 1745–1752. [Google Scholar] [CrossRef]
  21. Castro, L.P.; Minh, L.T.; Tuan, N.M. New convolutions for quadratic-phase Fourier integral operators and their applications. Mediterr. J. Math. 2018, 15, 13. [Google Scholar] [CrossRef]
Figure 1. Plot of the real part and imaginary parts of ψ Λ , 0 , 0 corresponding to Λ = ( 0 , 1 , 0 , 0 , 0 ) .
Figure 1. Plot of the real part and imaginary parts of ψ Λ , 0 , 0 corresponding to Λ = ( 0 , 1 , 0 , 0 , 0 ) .
Axioms 12 00927 g001
Figure 2. Plot of the real part and imaginary parts of ψ Λ , 0 , 0 corresponding to Λ = ( 1 3 , 1 , 0 , 1 5 , 0 ) .
Figure 2. Plot of the real part and imaginary parts of ψ Λ , 0 , 0 corresponding to Λ = ( 1 3 , 1 , 0 , 1 5 , 0 ) .
Axioms 12 00927 g002
Figure 3. Plot of the real part and imaginary parts of ψ , 0 , 0 corresponding to Λ = ( 0 , 1 , 0 , 0 , 0 ) .
Figure 3. Plot of the real part and imaginary parts of ψ , 0 , 0 corresponding to Λ = ( 0 , 1 , 0 , 0 , 0 ) .
Axioms 12 00927 g003
Figure 4. Plot of the real part and imaginary parts of ψ , 0 , 0 corresponding to Λ = ( 1 3 , 1 , 0 , 1 5 , 0 ) .
Figure 4. Plot of the real part and imaginary parts of ψ , 0 , 0 corresponding to Λ = ( 1 3 , 1 , 0 , 1 5 , 0 ) .
Axioms 12 00927 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gupta, B.; Kaur, N.; Verma, A.K.; Agarwal, R.P. Quadratic Phase Multiresolution Analysis and the Construction of Orthonormal Wavelets in L2(ℝ). Axioms 2023, 12, 927. https://doi.org/10.3390/axioms12100927

AMA Style

Gupta B, Kaur N, Verma AK, Agarwal RP. Quadratic Phase Multiresolution Analysis and the Construction of Orthonormal Wavelets in L2(ℝ). Axioms. 2023; 12(10):927. https://doi.org/10.3390/axioms12100927

Chicago/Turabian Style

Gupta, Bivek, Navneet Kaur, Amit K. Verma, and Ravi P. Agarwal. 2023. "Quadratic Phase Multiresolution Analysis and the Construction of Orthonormal Wavelets in L2(ℝ)" Axioms 12, no. 10: 927. https://doi.org/10.3390/axioms12100927

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