Thesis Banouh 2
Thesis Banouh 2
net/publication/352005870
CITATIONS READS
0 466
1 author:
Hicham Banouh
Université de Bouira
10 PUBLICATIONS 26 CITATIONS
SEE PROFILE
All content following this page was uploaded by Hicham Banouh on 31 May 2021.
Faculté de Mathématiques
Thèse
Présentée pour l’obtention du diplôme de DOCTORAT 3ème cycle
En : MATHEMATIQUES
Spécialité : Contrôle Optimal des EDP
Par : BANOUH Hicham
Sujet
Thesis
rd
Presented to obtain a 3 cycle DOCTORATE degree
In : MATHEMATICS
SPECIALTY : Optimal Control of PDE
By : BANOUH Hicham
Subject
Dissertation
by
Hicham BANOUH
Thesis
Submitted in partial fulfilment of the
requirements for the degree of
Doctor of Mathematics
2020
© Copyright by
Hicham BANOUH
2020
Approved by
First Reader: BEN MABROUK Anouar, Advisor
Associate-Professor, Higher institute of Applied Mathematics and Computer
Sciences, University of Kairaouan, Tunisia,
-------------------------------------
-------------------------------------
-------------------------------------
-------------------------------------
-------------------------------------
-------------------------------------
-------------------------------------
Acknowledgments
Firstly, I would like to thank my advisor, Dr. Anouar Ben Mabrouk, for his
unlimited patience, knowledge, guidance, support, and encouragement. You
have inspired me and helped me to grow as a mathematician. I cannot thank you
enough.
I thank my co-advisor Professor M’hamed Kesri for his great help during those
three last years and for the very usefull conversations we had.
I am very thankful to the chairman and the the members of the Thesis committee
for thier reports and thier presence.
To Mom and my late Dad, thank you for instilling in me a love of learning and
discovery. It has taken me far and served me well. Thank you for your
unconditional love, encouragement, and for putting up with my complaints.
Finally to Mina, Mazigh and Mayas : thank you for your constant love, support
and patience. Your unwavering belief in me kept me going when things were
toughest. I could not have done this without you.
Hicham BANOUH
KESRI M’hamed,
ABSTRACT
The project falls within the general framework of harmonic analysis on Clifford
algebras. We propose more precisely to extend the famous Heisenberg uncertainty principle to
the context of Clifford algebras by applying the so-called Clifford wavelets.
In the present work, our aim is to study and establish a new Heisenberg uncertainty
principle based on Clifford wavelet transform. We recall that the majority of uncertainty
principles in their different forms are based essentially on Fourier transform and the
wavelet uncertainty principles already established in the literature did not been extended
to general Clifford algebra framework.
ⴰⴳⵣⵓⵍ
ⴰⵙⴻⵏⴼⴰⵔ ⴰⴳⵉ ⵢⴻⵜⵜⵉⵍⵉ ⴷⴻⴳ ⵙⵙⴻⴳⵔⴻⵔ ⴰⵎⴰⵜⵓ ⵏ ⵜⴰⵙⵍⴻⴹⵜ ⵜⴰⵎⵜⴰⵡⴰⵏⵜ ⴷⴻⴳ ⵓⵊⴱⴰⵔⴻⵏ ⵏ
Clifford. ⵏⵙⵙⵓⵎⴻⵔ ⴷⴻⴳ ⵓⵇⴷⵉⵛ ⴰⴳⵉ ⴰⵏⵙⵎⴻⵖⵔⴻⴷ ⴰⵣⵡⵉⵔ ⵏ ⵜⴻⴳⵍⴰⴳⵉⵍⵜ ⵏ Heisenberg ⴰⵎⵓⵍⴰⵏ
ⴷⴻⴳ ⵓⵙⵙⴻⴳⵔⴻⵔ ⵏ ⵓⵊⴱⴰⵔⴻⵏ ⵏ Clifford ⵙ ⵓⵙⵏⴻⵙⵏ ⵏ ⵜⵉⵎⴷⴻⵙⵡⴰⵍ ⵜⵉⵎⴻⵥⵢⴰⵏⵉⵏ ⵏ Clifford.
ⴷⴻⴳ ⵜⵙⴻⵔⵙⵉⵜ ⴰⴳⵉ ⵏⴻⵜⵜⵏⴰⴷⵉ ⴰ ⵏⴻⵍⵎⴻⴷ ⴰⴽⴻⴷ ⴰⴷⵏⴰⴼ ⴰⵣⵡⵉⵔ ⵏ ⵜⴻⴳⵍⴰⴳⵉⵍⵜ ⵏ Heisenberg
ⴰⵎⴰⵢⵏⵓⵜ ⵢⴻⵜⵜⵡⴰⴱⵏⴰⵏ ⵖⴻⴼ ⵜⴰⴱⴻⴷⴷⵉⵍⵜ ⵏ ⵜⵉⵎⴻⵥⵢⴰⵏⵉⵏ ⵏ Clifford. ⴰⴷ ⵏⴻⵙⵎⴻⴽⵜⵉ ⴱⴻⵍⵍⵉ ⵜⵉⴳⴷⵉ
ⵏ ⵉⵣⵡⵉⵔⴻⵏ ⵏ ⵜⴻⴳⵍⴰⴳⵉⵍⵜ ⴷⵉ ⵜⴰⵍⵖⵉⵡⵉⵏ ⵏⵙⴻⵏ ⵉⵎⴳⴰⵔⴰⴷⴻⵏ ⵜⵜⵡⴰⴱⵏⴰⵏⵜ ⵖⴻⴼ ⵜⴱⴻⴷⵉⵍⵜ ⵏ Fourier
ⴻⵔⵏⵓ ⵉⵣⵡⴰⵔⴻⵏ ⵏ ⵜⴻⴳⵍⴰⴳⵉⵍⵜ ⵢⴻⵜⵜⵡⴰⴼⴻⵏ ⴷⵉ ⵜⵙⴻⴽⵍⴰ ⵓⵔ ⵜⵜⵡⵙⵎⴻⵖⵔⴻⵏ ⴰⵔⴰ ⵖⴻⵔ
ⵙⴰⵙⵙⴻⴳⵔⴻⵔ ⵏ ⵓⵊⴱⴰⵔⴻⵏ ⵏ Clifford ⴰⵎⴰⵜⵓ.
ملخص
نقترح على نحو أكثر دقة توسيع.يندرج المشروع في اإلطار العام للتحليل التوافقي في جبر كليفورد
بمويجات كليفورد. مبدأ عدم اليقين الشهير لهايزنبيرغ في سياق جبر كليفورد من خالل تطبيق ما يسمى
تم تطوير تحليل كليفورد كامتداد للتحليل التوافقي الكالسيكي حيث تم تمديد،في الفيزياء الرياضية
مفاهيم مثل تحويالت فورييه و المويجات في حالة جبر كليفورد
هدفنا هو دراسة وإنشاء مبدأ جديد لهايزنبيرغ لعدم اليقين يستند إلى تحويل مويجات,في العمل الحالي
وأن مبادئ، ونذكر أن معظم مبادئ عدم اليقين في أشكالها المختلفة تقوم أساسا على تحويل فورييه.كليفورد
عدم اليقين التي وضعت بالفعل في المؤلفات لم تمتد إلى إطار الجبر العام لكليفورد.
RESUME
Le projet s'inscrit dans le cadre général de l'analyse harmonique sur les algèbres de
Clifford. Nous proposons plus précisément d'étendre le fameux principe d'incertitude
d’Heisenberg au contexte des algèbres de Clifford en appliquant les ondelettes de
Clifford.
Dans ce travail, notre objectif est d'étudier et d'établir un nouveau principe d'incertitude
d’Heisenberg basé sur la transformation en ondelettes de Clifford. Nous rappelons que la
majorité des principes d'incertitude sous leurs différentes formes sont basés
essentiellement sur la transformation de Fourier et que les principes d’incertitudes relatifs
aux ondelettes n’ont pas été etendus aux algèbres de Clifford générales.
Contents
1.1 Introduction
One can never know with perfect accuracy both of those two important
factors which determine the movement of one of the smallest particles
— its position and its velocity. It is impossible to determine accurately
both the position and the direction and speed of a particle at the same
instant.
next Plancherel like rules. These findings have been applied next to prove different
uncertainty principles including Heisenberg-Weyl’s, Hardy’s, Beurling’s and loga-
rithmic ones in the case of the new quaternionic offset linear canonical transform.
Recently El-Haoui and Fahlaoui established in [78] several uncertainty inequalities
in the real Clifford algebra Rp,q such as Hausdorf-Young inequality and qualitative
uncertainty principles of Donoho-Stark.
In [81] expansions of signals with respect to Gabor wavelets and short time
Fourier transform have been investigated. Using Heisenberg group techniques sta-
ble iterative algorithms for signal analysis and synthesis have been developed.
These algorithms have been shown to be convergent for a variety of norms and
compatibility with the time-frequency localization of signals has been proved.
While the Clifford geometric algebra Fourier Transform is global, the author
in [107] introduced the local Clifford geometric algebra wavelet concept using the
similitude group SIM (n). As an explicit example, the author introduced Clifford
Gabor wavelets. In [108] the same author derived a new directional uncertainty
principle for quaternion valued functions by means of quaternion Fourier transfor-
mation generalized to the case of Clifford geometric algebras.
In [106] basic concept multivector functions and their vector derivative in geo-
metric algebra have been introduced. Concepts of Fourier transform and Clifford
and some useful properties have been also investigated in the same framework of
geometric algebras. An uncertainty principle has been next developed for many
cases of Clifford wavelets and shown to be useful for signal processing.
In [111] a generalization of the Fourier transform in some Clifford geometric
algebras has been extended and adopted for real Clifford geometric algebra Fourier
transform. This has been applied next to to define and prove the uncertainty
Outline
In [177], the authors gave a numerical method for solving non-linear optimal
control problems with inequality constraints by Legendre wavelet approximations.
In his PhD thesis, the author in [92] studied the interaction between problems
in control theory for partial differential equations and inequalities of the uncer-
tainty principle type. He also investigated the connection between compactness
of localization operators and uncertainty principles from an abstract harmonic
analysis perspective and gave results applied to the wavelet transform.
In [129], the authors investigated performance limitations and trade-offs in the
control design for linear time-invariant systems and showed that control specifi-
cations in time domain and in frequency domain are always mutually exclusive
determined by uncertainty relations.
For the moment, a theory of controllability and observability of Clifford algebra
valued differential systems is mainly studied and applied for the special case of the
algebra of quaternions H. For example in [216, 217] and [215] the authors used
quaternion valued equations to determine the best attitude for atmospheric entry
(see also [1, 121]).
In this work we obtain anew results stating that the product of variances be-
tween Clifford-Fourier transform denoted by fb and the Clifford wavelet transform
Tψ [f ] (w.r.t an admissible Clifford algebra-valued mother wavelet ψ) of a square
integrable function f is lower bounded.
Z Z
da 1
q n o
( kbk Tψ [f ] (a, ·, s)k22 ds) 2 ξk fb ≥ 2n+1 π n Aψ kf k22 + 2 |hf1 , f2 i|
Spin(n) R+ an+1 2
where
1 da
R R R
ψ a,b,s (x)∂bk Tψ [f ] (a, b, s) an+1
f1 (x) =
dV (b)ds
Aψ Spin(n) Rn R+
1 da
R R R
ψ a,b,s (x)bk Tψ [f ] (a, b, s) an+1
f2 (x) =
dV (b)ds
Aψ Spin(n) Rn R+
These results will be detailed in Chapter 4. Moreover, we recall that they have
been partially published in the following papers [19, 18]:
We intend that this result will be extended to stronger inequalities. Some future
directions are exposed in the last chapter.
2.1 Introduction
the signal and the time duration of that segment. Low-frequency pieces tend to
last a long time, whereas high frequencies occur in general for a short moment only.
For example, vowels have low frequency and last for a long period of time, whereas
consonants short bursts of high frequency. Fourier transform (Fourier 1878) gives
us only informations about the frequency domain (symbolized by the variable ξ)
but no information on time localization (the variable x). For this reason, we need a
time-frequency representations : transform the signal from a one variable function
to a function of two variables : time and frequency as the Short Fourier Transform
or Gabor Transform (Denis Gabor in [86]) and the wavelet transform.
We obtain wavelets by starting with a function ψ of the real variable x. This
function is called a mother wavelet if it is well localized and oscillating. (It resem-
bles to a wave because it oscillates, and it is a wavelet because it is localized). The
localization condition is expressed in the usual way by saying that the function
decreases rapidly to zero as |x| −→ ∞. The second condition suggests that ψ
vibrates like a wave. Mathematically, we require that the integral of ψ be zero and
that preferably the other first N -moments of ψ also vanish.
2.2 Wavelets on R
To deal with and/or to conduct wavelet analysis of functions we usually need and
start with one source function which will be called next the mother wavelet and
which plays the role of the analysing source. Such function should satisfy several
assumptions to be able to analyse functions next.
2.2.1 Admissibility
The analysing function must be square integrable But to have the CWT well
defined we add a condition on that analysis wavelet. This condition (admissibility
property) assures that the CWT can be inverted and so we can reconstruct the
signal again.
This will filter the polynomial components (at least of order ≤ N ) which is usually
the smoothest part of the signal and will show only the singularities and the sharp
parts. For example, if the analysing wavelet has vanishing moment of first order,
then the linear aspects of the signal will be ignored.
A next example is the Mexican Hat wavelet or Marr wavelet [144] obtained as the
2nd derivative of a Gaussian function
d2 − x2 x2
ψM H (x) = − 2 (e 2 ) = (1 − x2 )e− 2 .
dx
An important example is also due to Morlet [208, 93] and is based on the mother
wavelet
x2
ψM (x) = eic0 x− 2 .
This wavelet is closely related to human perception, both in the processes of au-
dition and vision.
a > 1, the support of ψ a,b will be wider than that of ψ and will correspond to the
√
lower frequencial part of the signal. The factor 1/ a ensures the conservation of
the L2 -norm of both ψ and ψ a,b and plays the role of a normalization constant.
All these copies ψ a,b have the same L2 -norm as the mother wavelet ψ. Indeed,
ψ a,b L2 (R,C,dx)
= kψkL2 (R,C,dx) .
which yields that all these daughter wavelets are admissible also.
where ψ
fa (x) = √1 ψ( −x ).
a a
It is interpreted as a filter with a function of zero
momentum.
We may remark that if the wavelet ψ has N -vanishing moments, the CWT
Tψ [f ] (a, b) of f will have the same order of magnitude as aN +1+1/2 (See [68, Page
102]). Indeed,
Z
Tψ [f ](a, b) = f (x)[ψ a,b (x)]c dx
R
x−b c x−b c
Z Z
1
= √ { f (b)[ψ( )] dx + f 0 (b)(x − b)[ψ( )] dx
a R a R a
Z 00 Z (N )
f (b) x−b c f (b) x−b c
+ [ψ( )] dx + · · · + (x − b)N [ψ( )] dx + · · · }
R 2 a R N! a
We see that the N first terms are zero. Consequently, putting x = at + b in the
(N + 1)-integral we get an approximation
|Tψ [f ]| ∼ aN +1+1/2 .
From the definition of the CWT we may easily notice the following properties. Let
ψ and φ be two admissible wavelets, f, g ∈ L2 (R, C, dx) and α, β ∈ C.
• Linearity rule
• Anti-linearity
c
• Tψ [f ] (a, b) = Tψ [f ] ( a1 , −b .
a
Besides as in Fourier analysis, we have also the possibility to reconstruct the anal-
ysed functions even non periodic one with analogues of Dirichlet, Parseval and
Plancherel rules by using the CWT.
Tψ [f ] (a, b) = f, ψ a,b
D E
= f, ψ
b da,b , by Parseval formula for Fourier transform.
√
Z
= fb(ξ)[ ae−iξb ψ(aξ)]
b c
dξ
ZR
√ iξb c
= ae fb(ξ)[ψ(aξ)]
b dξ
R
and so
√ −iξb
Z
c
[Tψ [g](a, b)] = ae [b g (ξ)]c ψ(aξ)dξ.
b
R
This yields that
Z Z +∞ Z Z +∞ Z
c da
√ iξb c
Tψ [f ](a, b)[Tψ [g](a, b)] 2 db = { ae fb(ξ)[ψ(aξ)]
b dξ}
R 0 a R 0 R
√ −iξb
Z
da
×{ ae [b g (ξ)]c ψ(aξ)dξ}
b db (2.5)
R a2
Z Z +∞ Z Z
da
= { fb(ξ)[bg (ξ)]c [ψ(aξ)]
b cb
ψ(aξ)dξdξ} db.
R 0 R R a
Now, we characterize the image of L2 (R, C, dx) by a given CWT and introduce
thus some reproducing kernels associated to the CWT. Let Hψ = L2 (R+ dadb
∗ ×R, Aψ a2 ).
2.3 Conclusion
In this chapter, we presented the wavelets tool used for analysing functions. We
saw some proprieties of those wavelets and proved that under certain conditions,
the continuous wavelet transform is a an invertible operator. We will use similar
characteristics in the Clifford algebra framework which will be the subject of the
next chapter.
3.1 Introduction
ied for the analysis of higher dimensional signals such as colour images based on
Clifford Hermite wavelets. In [88] a sophisticated model based on Marr wavelet
kernel has been developed and applied on samples of intensity values for each pixel
in an image to estimate the probability density function of the pixel intensity. Marr
wavelet has been also applied in [211] to detect and characterize two-dimensional
vortex for a synthetic flow and propeller wake. See also [87], [202].
In [39], the authors studied the historical development of quaternion and Clif-
ford Fourier transforms and wavelets. Basic concepts have been revisited and
mathematical formulations has been enlightened. Hypercomplex Fourier trans-
forms and wavelets has been revisited with overviews on quaternion Fourier trans-
forms, Clifford Fourier transforms, quaternion and Clifford wavelets.
In [48], [49], the authors introduced new definition for general geometric Fourier
transform covering some Clifford cases. They showed necessary constraints to
obtain linearity, scaling and shift theorem. Applications in image/signal processing
and mathematical imaging vision in general have been discussed.
In this chapter we propose to state the most relevant results on Clifford alge-
bras and the theory of functions taking value in them. In the first section, we give
some preliminary concepts on Clifford algebra as their definitions, the matrix rep-
resentation of Clifford algebras and the Spin group (which will be used to describe
rotations in Rn ). In the forth section, we give a review on monogenic functions
which are counterparts of the holomorphic functions on the complex plan. In the
two final sections, we present the generalizations of the classical Fourier transform
and the wavelet transform in the Clifford algebras framework.
In the literature, there are different ways to introduce Clifford algebras. We will
present in this part some of them. The readers may be referred to the list of
references provided in this document for more details and other constructions.
It starts from the Grassmann exterior algebra [91] Rn of the linear space Rn ,
V
ei ei = e2i = 1.
In [179] another construction of Clifford algebras has been provided. Let (V, Q) be
a n−dimensional quadratic space and A an associative algebra with the following
rules on addition and multiplication :
x2 = Q(x)
xy + yx = 2B(x, y)
ei ej + ej ei = 2B(ei , ej ).
ei ei = e2i = 1
∀1 ≤ i ≤ n.
This last construction resembles to the constructions formulated in [71] and [140].
The first and simplest example is the Clifford algebra R0,1 where the elements are
written on the form x = x0 1 + x1 e1 where x0 , x1 ∈ R and e21 = −1. We have in
fact By making the identity e1 = i, the imaginary unit we have the isomorphism
R0,1 ' C.
The second example is the well known quaternions algebra [100] denoted R0,2 with
its elements q = q0 1 + q1 e1 + q2 e2 + q12 e12 where qi ∈ R for i = 0, 1, 2, q12 ∈ R and
where e2i = −1, and ei ej = −ej ei for i 6= j. Choosing e1 = i, e2 = j, e12 = k the
elements of the canonical basis of the Euclidean space R3 we obtain the Hamilton
algebra constructed in [99] known also as the quaternion algebra which satisfies
H ' R2,0 .
In the later we mainly consider the Clifford algebra Rn of the euclidean space Rn .
Any element of the Clifford algebra Rn can be written as
X
a= aA eA
A
where A is an arbitrary ordered multi-indices
= a0 + (a1 e1 + a2 e2 + a3 e3 )
+ a123 e123
and the a’s are real numbers. Let k be the cardinality of A (|A| = k), we have
X Xn X
a= aA eA = aA e A .
A k=0 k=|A|
The subspace Rkn = SpanR {eA | |A| = k} will be called subspace of grade k.
For example, the subspace of grade 0 is the field R whose elements are called
scalars, the one of grade 1 is the vector space Rn composed of vectors, the elements
of one of grade 2 are called bivectors and finally, the one dimensional subspace of
So that
a = a∅ 1 + (a1 e1 + .... + an en )
scalars vectors
+ · · · + (a123...n e123....n ).
pseudo−scalar
n
X
We may write a = [a]k where [a]k is the projector of a on Rkn . The operator
k=0
[.]k : Rn −→ Rkn for k = 0, 1, · · · , n satisfy for a, b ∈ Rn and λ ∈ R
[[a]k ]k = [a]k
A aA eA ∈ Cn or λ = a + ib ∈ Cn with a, b ∈ Rn .
P
The Clifford algebra Rn is Z2 −grader, which means that it is the direct sum of
an even and odd subspaces
M M
−
Rn = R+
n ⊕ Rn = Rkn ⊕ Rkn
k even k odd
and any Clifford number a can be written as
a = [a]+ + [a]−
where [a]± ∈ R±
n . These two components
satisfy to the
following inclusions [163]
−
R+ +
R+
n Rn n Rn
⊂ Rn and
+
⊂ R−
n.
R− R
n n
−
R−n Rn
+
Also, as Rn is the direct sum of the even and odd subspaces and being of dimension
2n we have
−
dim(R+
n ) = dim(Rn ) = 2
n−1
.
Definition 3.1. The centre of the Clifford algebra Rn is the set of elements which
commute with all the other elements of the algebra. It will be noted by Z
We have
for n even
R
Z(Rn ) =
for n odd
R ⊕ Re123···n
One remarkable concept in Clifford algebra is its attempt to mix between the
Grassmann exterior algebra and the Hamilton quaternion one. Let x and y be two
vectors. We define the Clifford product by
xy = x · y + x ∧ y
Xn
x·y = − xj yj = − < x, y >Rn
j=1
X
x∧y = ei ej (xi yj − xj yi )
i<j
Observing that y ∧ x = −x ∧ y we get
yx = y · x + y ∧ x = x · y − x ∧ y.
So that
1 1
x · y = (xy + yx) and x ∧ y = (xy − yx).
2 2
and
1
x ∧ Ak = [xAk ]k+1 = (xAk + (−1)k Ak x). (3.2)
2
In this case the Clifford product will be
xAk = x · Ak + x ∧ Ak .
We can even expand to the Clifford product of two k−grade,Ak , and l−grade, Bl ,
elements (see [105, page 6])
and
Ak ∧ Bl = [Ak Bl ]|k+l|
On the real Clifford algebra Rn we may also define different types of involutions.
A first type is the main involution. It extends the vectorial reflection through
the origin to the whole algebra. Sometimes called grade involution, it is denoted
e and has following proprieties : for λ ∈ R, x ∈ Rn we have
a=a
e
e and a = [a]+ − [a]−
e
(ab)
g =e aeb and ^
(a + b) = e
a + eb.
Next, the Clifford conjugation is defined as the composition of the main invo-
lution and the reversion as follows,
a∗ .
ae∗ = e
It corresponds to the complex and quaternion conjugation in the case of R0,1 ' C
and R0,2 ' H respectively. The superposition will be denoted by so we have
a∗ .
a = ae∗ = e
So that
|A|(|A|+1)
λ = λ, x = −x, ej = −ej , eA = (−1) 2 eA
(ab) = ba, (a + b) = a + b, a = a
Let us now consider the two following groups, the formed by the products of
invertible vectors (known as the versor group)
k
Y
Γ2p,q := { vi : v1 , v2 , · · · , vk ∈ (Rp,q )× }
i=1
and
Γ1p,q := {b ∈ R× e −1 for all vector x ∈ Rp,q }.
p,q ; bxb
N : Rp,q −→ Rp,q
a 7−→ aa
Definition 3.2. We define the P in(p, q) group as the subgroup of Γ1p,q of elements
b for which N (b) = ±1. It can be shown that it is also the kernel of the mapping
N : Γ1p,q −→ R∗ .
The same way we define Spin(p, q) as the subgroup of P in(p, q) of product of
an even number of elements which the spinor norm equals ±1 or
Spin(p, q) = P in(p, q) ∩ R+
p,q
2l
Y
= {s ∈ Rp,q ; s = ω j with ω 2j = ±1, 1 ≤ j ≤ 2l}. (3.5)
j=1
ea eb = ea+b
In the latter, we consider only functions with values on a Clifford algebra but
where the variable is a vector in Rn taken as a part of the Clifford algebra.
Xn
We put x = (x0 , x) ∈ Rn+1 and where x = ei xi ∈ Rn so we can write
i=1
n
X Xn
x = e 0 x0 + e i xi = ei xi
i=1 i=0
The vector space Rn can be seen as the hyper-plan {x = (x0 , x) ∈ Rn+1 : x0 = 0}.
Consider functions defined on the (n + 1)-dimensional vector space Rn+1 :
f : R ⊕ Rn ' Rn+1 −→ Rn ( or Cn ).
It may be expressed as
X
f (x) = eA fA (x), (3.7)
A
where fA are real-valued functions and A ⊂ {1, 2, · · · , n}. Its conjugate f is given
by
X
f (x) = eA fA (x)
A
We mainly use the complex Clifford conjugation (see 3.3. The inner product (3.8)
satisfies the Cauchy-Schwartz inequality
< f, g >L2 (Rn ,Rn ,dV (x)) ≤ kf kL2 (Rn ,Rn ,dV (x)) kgkL2 (Rn ,Rn ,dV (x)) . (3.10)
(i) (a + b) · ∇f = a · ∇f + b · ∇f for a, b ∈ Rn .
∂x : C (r) −→ C (r−1)
n
∂f X
f 7−→ ∂x f = ek (3.12)
k=1
∂ xk
and the generalized Cauchy-Riemann operator or Weyl operator as
∂
Dx = + ∂x (3.13)
∂x0
also known as the Fueter-Delanghe operator (see [84],[69]). For the sake of sim-
plicity we denote from now on
∂
= ∂i
∂xi
so the operators (3.12) and (3.13) become
Xn
∂x = ei ∂i and Dx = ∂0 + ∂x
i=1
and their conjugates will be
n
X
∂x = ei ∂i and Dx = ∂0 − ∂x .
i=1
They have an action from the left
X
∂x f (x) = ei eA ∂i fA (x)
i,A
and right
X
f ∂x (x) = eA ei ∂i fA (x).
i,A
(see the formulas (3.2)). By the same way
Dx f = ∂0 f + ∂x f and
f Dx = f ∂0 + f ∂x .
parts
1
div f = ∂x · f = (∂x f + f ∂x )
2
and
1
curl f = ∂x ∧ f = (∂x f − f ∂x ).
2
We consider the action of the group Spin(n) on a Clifford-valued function f (see
[89, Thm 3.6]) given by
∂x Ls = Ls ∂x .
Proof. We have
n
X
∂x Ls f (x) = ei ∂i {sf (sxs)s}
i=1
Xn
= ei s∂i {f (sxs)} s
i=1
n
X
= s( ei ∂i {f (sxs)})s
i=1
= s∂x f (sxs)s
= Ls ∂x f (x).
Remark 3.8. We remark (see [168, pp. 139]) that the Dirac operator ∂x maps even
parts to odd parts and odd parts to even parts (see the splitting in (3.3))
and
∂x [f (x)]− = [∂x f (x)]+
As such, f is monogenic if and only if ∂x [f (x)]+ = 0 and ∂x [f (x)]− = 0. More
generally, the equation
∂x f = g
where f and g are two Clifford algebra-valued functions defined on Rn , can be
splitted into a system of two equations
∂x [f ]+ = [g]−
∂x [f ]− = [g]+
resp.
f Dx = 0.
In the special case where the function f takes values only on the vector space
Rn taken as a component of the real Clifford algebra Rn , so it can be written as
X n
f (x) = ei fi (x)
i=0
where, by analogy with (3.7), the (n + 1) functions fi are real valued. In this case,
the monogenicity condition (3.14) becomes
X
ei ej ∂i fj (x) = 0.
i,j
Or as f is vector-valued ∂x f = ∂x · f + ∂x ∧ f = 0 so we have
n
X
∂x · f = − ∂xj fj = −div(f )
j=1
X
∂x ∧ f = eij ∂xi fj − ∂xj fi = rot(f )
i<j
known as the Riesz system which describes an irrotational flow [203]. One can
show that
where ∆n and ∆n+1 are the Laplacian in Rn and Rn+1 respectively. This means
that a monogenic function is also a harmonic one and hence infinitely differentiable;
even more, it is an multivector-valued analytic function in Ω ⊂ Rn and so each
of its components fA is real-analytic in its domain. Now, we will try to write the
previous operators in terms of spherical coordinates. For that, we recall that we
can write x = rη with η = ni=0 ei ηi and ηi = xri for i = 0, 1, · · · , n. More explicitly
P
where ∂η spherical Dirac operator acts on the sphere S n . We have [37, Page 49]
n ∂η
X ∂θi
· ∂θi
∂η = 2 and ∂η = x ∧ ∂x .
∂ηi
i=1
∂θi
Putting
Consequently,
1 1 f?
∂x = η ∂r + Γη = (∂r + Γ )η
r r η
and
1f
∂x = ∂r + Γη η
r
Observing that
∆n = ∂x ∂x = ∂x ∂x (3.16)
and that ηη = ηη = 1 we obtain
1 fη )∂r + 1 (Γ
∆n = ∂r2 + (Γη + Γ fη Γη − Γη ).
r r
By the same way we get
1 f?η )∂r + 1 (Γ
∆n = ∂r2 + (Γ?η + Γ f? Γ? − Γ?η ).
r r η η
This means that
1 fη )∂r + 1 (Γ
fη Γη − Γη ) = ∂ 2 + 1 (Γ? + Γ f? )∂r + 1 (Γ
∂r2 + (Γη + Γ r η η
f? Γ? − Γ? )
r r r r η η η
tions we get
∆η = (n − 1) I − Γ?η Γ?η .
n
X
Putting E = xi ∂i called the Euler operator, we obtain E = r∂r . As a result,
i=0
we get the following identities
• xDx = E + Γη ,
• x∂x + ∂x x = −2E − n,
• Γη = −
P
i<j eij (xi ∂j − xj ∂i ).
Proposition 3.11. [37, Page 48] If Ω is open and star-shaped with respect to the
origin and u : Ω → R is harmonic in Ω then the function
Z 1 Z 1
n−1 n−1
g(x) = f (x)e0 + t ∂x f (tx)xdt − t ∂x f (tx)dt
0 0 0
is left monogenic in Ω and its scalar part is precisely the function f .
As for the previous result, we seek a monogenic function starting with an ana-
lytic one. For this, we use the Cauchy-Kowalevski extension (CK-extension). The
idea behind the CK-extension is to characterize solutions of a system of partial dif-
ferential equations by the restriction of some of their derivatives to a sub-manifold
of co-dimension one.
Let f be an analytic function on an open set U of Rn , then we can always
monogenically extend it to an open set Ω of Rn+1 where U = Rn ∩ Ω (known as
the Cauchy-Kowalevski extension (see for instance the original works [51, 133] and
more recently [38, 40, 199, 51][37, Page 111] and [71, Page 151] ), then the solution
F of the system
∂x F (x0 , x) = −∂x F (x0 , x) in Rm+1
0
F (0, x) = f (x)
=0
Using the following result we can extend a holomorphic complex function into
a monogenic one.
Theorem 3.4.1. [189, 84] Let n be odd and let z = x + iy 7−→ f (x + iy) be
holomorphic function. Then the function
n−1
(∆n+1 ) 2 f (x0 + x)
is monogenic.
The authors in [171] proved a similar result for n even. Now, we present another
method of extending a function on Rn into a monogenic one. This is due to
[132, 172]
if and only if
Supp(fb) ⊂ B(0, R)
where
B(0, R) = {x ∈ Rn : |x| < R}
In this case we have
Z
1
F (x0 , x) = e(x0 , x, ξ)fb(ξ)dξ,
(2π)n Rn
where
and
ξ
1
χ± (ξ) = 1+i
2 |ξ|
Now, we give some results about integration of Clifford-valued function which are
similar to those that are well known in the theory of complex-valued functions as
a Stokes formula and a Cauchy representation formula.
dσx = η(x)dS(x)
As a consequence we have.
Proof. We only prove the result for the left-monogenic function f . We consider the
sphere S n−1 (y, r) for r > 0 chosen small enough such that the disc whose boundary
f (x) − f (y)
Z Z
= n−1 dσx + f (y) dσx
S n−1 (y,r) x − y S n−1 (y,r)
Theorem 3.4.5 – Mean Value Theorem. Let y ∈ U and R > 0 and consider
D(y, R), the disk of centre y and radius R such that D(y, R) ⊂ U then
Z
1
f (y) = n f (x)dV (x).
R an D(y,R)
G(ξ)
b = G(ξ).
Proposition 3.14. The two operators H and exp(−i π2 H) are Fourier invariant
in the sense that
[) = H(fb)
H(f
and that
\ π π
exp(−i H)(f ) = exp(−i H)(fb).
2 2
Proof. We have
−1
Z
1
[)(ξ) =
H(f n (∆n f (x) + x2 f (x) + nf (x))e−ix·ξ dV (x)
(2π) Rn 2
2
−1 1
Z Z
−ix·ξ 1
= [ n
2
∂ f (x)e dV (x) + n x2 f (x)e−ix·ξ dV (x)
2 (2π) 2 Rn x (2π) 2 Rn
Z
1
+n n f (x)e−ix·ξ dV (x)]
(2π) 2 Rn
−1 h 2 b i
= ξ f (ξ) + ∂ξ2 fb(ξ) + nfb(ξ)
2
= H(fb)(ξ).
Lets calculate H
\ k (f ). We have from the first assertion :
\
H \
k (f ) = H {H k−1 (f )}
= HH\
k−1 (f )
\
= HH {H k−2 (f )}
= H2 H\
k−2 (f )
..
.
= Hk H
\ 0 (f )
= Hk fb.
This gives us
∞ k
\ π X −i π
exp(−i H)(f )(ξ) = 2
Hk fb(ξ)
2 k=0
k!
π
= exp(−i Hfb).
2
The idea behind the extension of the classical Fourier transform to Clifford
algebra-valued one resides in the generalization of the Hermite operator into a
multivector-valued one. We use the factorization of the Laplace operator by the
Dirac operator to obtain two new Clifford-valued operators that factorize the Her-
mite operator. This method has been developed in [35, 66, 30, Ch. 12]. For that
purpose we introduce the following operators (see [32])
1
O1 = (∂x − x) (∂x + x)
2
1
O2 = (∂x + x) (∂x − x)
2
This operators has the following properties
Proposition 3.15.
O1 = H + Γx
O2 = H − Γx + n
= 2 H + n2
O1 + O2
= 2 Γx − n2
O1 − O2
And
1
O2 = (∂x + x) (∂x − x)
2
1 2
∂x − x2 − n + x∂x − ∂x x + n
=
2
1
= H + (x∂x − ∂x x + 2n − n)
2
= H − Γx + n.
We have also
O1 + O2 = H + Γx + H − Γx + n
n
=2 H+
2
O1 − O2 = H + Γx − H + Γx − n
n
= 2 Γx − .
2
For the two operators O1 and O2 to be used in the definition of the new Fourier
As
O
1 = H + Γx
O2
= H − Γx + n
then
O
1f
d = Hf
c +Γ
dx f = O1 f
b
O
c2 = H{c −Γ dx f + nf = O2 f .
b b
Now, we are able to define the Clifford-Fourier transform
H+ + H− = O1 + O2 − n.
So we set
n
H
+ = O1 − 2
n
H−
= O2 − 2
hence
= H + Γx − n2
H
+
(3.21)
= H − Γx − n2 .
H−
One alternative for the operators H+ and H− is given in [66, Sec. 12.3.2] and in
[35, Def. 4.2]
H+
= O1
H−
= O2 − n
and in this case we have
H
+ = H + Γx
H−
= H − Γx .
Keeping our operators as in (3.21) we obtain in terms of the exponential operator
π
F+ = exp(−i H+ )
2
π h n i
= exp −i H + Γx −
2 2 (3.22)
π n π
= exp −i Γx − exp −i H
2 2 2
π n
= exp −i Γx − F
2 2
and in the same way
π n
F− = exp i Γx − F. (3.23)
2 2
and
Z
1 π n −ix·ξ
F− [f ](ξ) = n exp i Γξ − e f (x)dV (x). (3.24)
(2π) 2 Rn 2 2
We have defined the Clifford-Fourier transform as a pair of operators satisfying
F 2 = F+ F−
so we can write
1 1
F = F+2 F−2
where the square root is the Fractional Fourier Transform (see [170] and [66,
Ch.11]) given by
1 π
F+ = exp −i H+
2
4
1 π
F−2 = exp −i H− .
4
We can factorize the classical Fourier transform as the products
π π
F = exp −i H+ exp −i H−
4 4
π π
= exp −i H− exp −i H+
4 4
This transform can be inverted. Using (3.22) and (3.23) we have
π
−1
F+ = exp i H+
2
and
π
F−−1 = exp i H−
2
which gives
π n −1
F+−1 = exp i Γx − F
2 2
and
π n −1
F−−1 = exp −i Γx − F
2 2
The Clifford-Fourier transform has properties similar to those of the classical
Fourier transform.
F+ [f a + gb] = F+ [f ] a + F+ [g] b
and
F− [f a + gb] = F− [f ] a + F− [g] b.
n n
X ∂ X ∂
∂y = ei = ei λ = λ∂x .
i=1
∂yi i=1
∂xi
So finally,
x x
Γ λx = − ∧ ∂ λx = − ∧ λ∂x = Γx .
λ λ
and
Proof. We know from the proprieties of the classical Fourier transform that
Then
π n
F+ [•f (•)] (ξ) = exp −i Γξ − F [•f (•)] (ξ)
2 2
π n
= i exp −i Γξ − ∂ξ F [f ] (ξ)
(∞ 2 2 )
X (−i π )k n k
=i 2
Γξ − ∂ξ F [f ] (ξ). (3.25)
k=0
k! 2
We know from [162, Thm 2.4] that
∂ξ Γξ + Γξ ∂ξ = (n − 1)∂ξ .
Then
h i
Γξ ∂ξ = ∂ξ n − 1 − Γξ ,
and
n n
∂ξ = ∂ξ .
2 2
So
n n−2
Γξ − ∂ξ = ∂ξ − Γξ .
2 2
Applying Γξ − n
2
gives us
n
n n n−2
Γξ − Γξ − ∂ξ = Γξ − ∂ξ − Γξ
2 2 2 2
n−2 n−2
= ∂ξ − Γξ − Γξ .
2 2
and
and
(ξ ∧ x)
r
π −1/2
K− (x, ξ) = |x ∧ ξ| (1 − x · ξ)J1/2 (|x ∧ ξ|) + J3/2 (|x ∧ ξ|)(x · ξ)
2 |x ∧ ξ|
and for n = 6 we have
pπ −1/2 (x·ξ)2
K+ (x, ξ) = 2 {|x ∧ ξ| J1/2 (|x ∧ ξ|) + |x∧ξ| J3/2 (|x ∧ ξ|)
2
+|x ∧ ξ|−3/2 2(x · ξ)J3/2 (|x ∧ ξ|) + (ξ ∧ x)J3/2 (|x ∧ ξ|) + (x·ξ)
|x∧ξ|
(ξ ∧ x)J 5/2 (|x ∧ ξ|) }
and
pπ −1/2 (x·ξ)2
K− (x, ξ) = 2
{|x ∧ ξ| J1/2 (|x ∧ ξ|) + J (|x
|x∧ξ| 3/2
∧ ξ|)
(x·ξ)2
−|x ∧ ξ|−3/2 2(x · ξ)J3/2 (|x ∧ ξ|) + (ξ ∧ x)J3/2 (|x ∧ ξ|) + |x∧ξ|
(ξ ∧ x)J5/2 (|x ∧ ξ|) }.
In [65], the authors found a general expression for the kernels for all even dimen-
sions.
We introduce the concept of the Clifford-wavelet transform and some of its impor-
tant properties to be used later. Besides of the translation and dilation, we will use
rotation using the action of the spin group which is a double-cover of the special
orthogonal group in Rn . In this context, a function ψ ∈ L1 ∩L2 (Rn , Rn , dV (x)) will
< ∞.
Example 3.22. From [9] we know that for 0 < t < 1−n−2α
2
we have
Z
α
xk Gn,t,α+t (x) 1 + |x|2 dV (x) = 0
Rn
where Gn,t,α+t are the Clifford-Gegenbauer polynomials and specially
Z Z
2 α
α+t
(−1)t ∂xt 1 + |x|2
Gn,t,α+t (x) 1 + |x| dV (x) = dV (x)
Rn Rn
=0
and this is just the admissibility condition 3.6. So the functions x 7−→ ψn,t,α (x) =
α+t
(−1)t ∂xt (1 + |x|2 ) can be taken as mother wavelets. We call them Clifford-Gegenbauer Wavelets.
Example 3.23. In [28] the authors defined the so called Clifford-Laguerre Wavelets
as
(−1)k ∂xk exp(−|x|)|x|α+2k P +
< ψ a,b,s , f >L2 (Rn ,Rn ,dV (x)) = 0, ∀a > 0, b ∈ Rn and s ∈ Spin(n).
we obtain
Since
Z h i†
n
[
<ψ a,b,s , fb > 2 n
L (R ,Rn ,dV (x)) = a 2 eib·ξ s ψ(asξs)
b sfb(ξ)dV (ξ) = 0,
Rn
then necessarily
h i†
s ψ(asξs) sfb(ξ) = 0, ∀ξ ∈ Rn .
b
Recall now that for a fixed ξ 6= 0 in Rn (see [66], pages 48 and 49)
It results that
fb = 0 and so f = 0.
As for the real case, we define the Clifford Wavelet Transform as the projection
of the signal f on the set of admissible Clifford wavelets
• Covariance by translation
• Covariance by Dilation
1 • a b
Tψ n f( ) (a, b, s) = Tψ [f ] ( , , s).
λ2 λ λ λ
Proof. We have
†
s(x − b)s
Z
1
Tψ [f (• − c] (a, b, s) = n s ψ sf (x − c)dV (x).
a 2 Rn a
Put y = x − c =⇒ x − b = y − (b − c) and dV (x) = dV (y) so
s(y − (b − c))s †
Z
1
Tψ [f (• − c] (a, b, s) = n s ψ sf (y)dV (y)
a 2 Rn a
= Tψ [f ] (a, b − c, s).
s(λy − b)s †
•
Z
1 λ n
Tψ n f( ) (a, b, s) = ( ) 2 s ψ sf (y)dV (y)
λ2 λ a Rn a
" !#†
s(y − λb )s
Z
1
= a n s ψ a sf (y)dV (y)
( λ ) 2 Rn λ
a b
= Tψ [f ] ( , , s).
λ λ
Finally, for the case of the Spin rotation first recall the action of the Spin
group (3.4.1.2)
Ls : f (x) → sf (s̄xs)s
Then we have
†
s(x − b)s
Z
1
Tψ [Lt f ] (a, b, s) = n s ψ stf (txt)tdV (x)
a 2 Rn a
put y = txt =⇒ x = tyt and dV (x) = tdV (y)t. So
" !#†
s(tyt − b)s
Z
1
Tψ [Lt f ] (a, b, s) = n s ψ stf (tyt)tdV (y),
a 2 Rn a
having tt = tt = 1, then
" !#†
s(tyt − b)s
Z
1
Tψ [Lt f ] (a, b, s) = n s ψ stf (tyt)tdV (y)
a2Rn a
" !#†
st(y − tbt)ts
Z
1
=t n ts ψ {st} f (tyt)tdV (y)t
a 2 Rn a
Knowing all that, we can now introduce a result that permits us to invert the
Clifford wavelet transform
Put
h i†
Φψ (a, s, ξ) [f ] (−b) = ψ(asξs) sf (ξ) (−b)
b b
and similarly
h i†
Φψ (a, s, ξ) [g] (−b) = ψ(asξs) sb
b g (ξ) (−b).
We obtain
n n
Tψ [f ] (a, b, s) = a 2 s(2π) 2 Φψ\
(a, ξ, s) [f ] (−b)
and
n n
Tψ [g] (a, b, s) = a 2 s(2π) 2 Φψ\
(a, ξ, s) [g] (−b).
Applying Parseval formula we get
D E
\
Φψ (a, \
•, s) [f ], Φψ (a, •, s) [g] = hΦψ (a, •, s) [f ] , Φψ (a, •, s) [g]i .
And so
Z Z Z
1 † da
[Tψ [f ] , Tψ [g]] = (Φψ (a, ξ, s) [f ] (ξ)) Φψ (a, ξ, s) [g] (ξ)dV (b) ds
(2π)n Aψ Rn a
Spin(n) R+
Z Z (Z h †
† h
†
)
1 i i da
= n
( ψ(asξs)
b sfb(ξ) ψ(asξs)
b sb
g (ξ)dV (ξ) ds
(2π) Aψ Rn a
Spin(n) R+
Z Z Z h i† i†
1 h da
= n
f (ξ) sψ(asξs) ψ(asξs) sb
b b b g (ξ)dV (ξ) ds
(2π) Aψ Spin(n) Rn a
R+
Z h †
Z Z † da
1 i h i
= f
b(ξ) s ψ(asξs)
b ψ(asξs)
b s ds gb(ξ)dV (ξ).
(2π)n Aψ Rn Spin(n) a
R+
Observing now that
Aψ
Z Z h i† da
sψ(asξs)
b ψ(asξs)
b s ds = , (3.30)
Spin(n) a (2π)n
R+
we get immediately
Z h i†
[Tψ [f ] , Tψ [g]] = fb(ξ) gb(ξ)dV (ξ)
Rn
We can also say that it is an isometry between the two spaces of integrable
functions L2 (Rn , Rn , dV (x)) and L2 (R+ × Rn × Spin(n), A−1
ψ a
−(n+1)
dadV (b)ds).
This is an analogue to Parseval’s formula and as a result we have also a Plancherel’s
formula
Z Z Z
da
|Tψ [f ] (a, b, s)|2 n+1
dV (b)ds = Aψ kf k22 . (3.32)
a
Spin(n) Rn R+
As a result of the last Proposition and as in the real case, we have here a Clifford-
wavelet reconstruction formula.
In other words, the Clifford wavelet transform decomposes the signal f in terms
of the analysing wavelets ψ a,b,s with coefficients Tψ [f ].
Proof. Let f and g two square integrable Clifford-valued functions with Clifford
wavelet transforms (with respect to a mother wavelet ψ) Tψ [f ] and Tψ [g] respec-
kernel.
3.7 Conclusion
4.1 Introduction
The uncertainty principle is central for information processing and quantum physics.
In the Clifford algebras framework, the uncertainty principle provides us data
about the way a multivector valued function and its Clifford-Fourier transform are
related. We will see that the uncertainty principle for the Clifford wavelet trans-
form establishes a lower bound of the product of the variances of Clifford wavelet
transform of a square integrable multivector-valued function and its Clifford-Fourier
transform.
First, we recall some results concerning the classical Heisenberg uncertainty prin-
ciple. For more backgrounds on the uncertainty principle, its variants, Fourier and
wavelet transforms on the Euclidean space Rn the readers may be referred also to
[122], [166], [191] and [201]. Mathematically, the uncertainty principle states that a
non-zero function and its Fourier transform cannot both be sharply localized. The
next theorem formally summarizes the Heisenberg’s Uncertainty Principle.
Which reads as
1
kAf k2 kBf k2 ≥ |h[A, B] f, f i| . (4.2)
2
which completes the proof.
Lets apply theorem (4.2.1) to the case of the Fourier transform. We mainly review
the results of [83] and [118]. Let f ∈ L1 ∩ L2 (Rn , Rn , dV (x)) We define for k ∈
{1, 2, · · · , n} the two families of operators
Ak f (x) = xk f (x)
whence
∂d
xk f = ξk fb .
2 2
kAk f k2 = kxk f k2
and
kBk f k2 = ξk fˆ
2
Applying theorem (4.2.1) to both operators Ak and Bk yields :
1
kAk f k2 kBk f k2 ≥ |< [Ak , Bk ] f, f >| .
2
As the commutator
[Ak , Bk ] f = Ak Bk f − Bk Ak f
= −f (x)
we obtain finally
1
kxk f k2 ξk fb ≥ kf k22 . (4.3)
2 2
Theorem 4.2.2. For A and B two symmetric operators on a Hilbert space H and
for f ∈ L2 (H, dx) we have
1
q
2
kAf k2 kBf k2 ≥ |h[A, B] f, f i|2 + [A, B]+ f, f (4.4)
2
where [A, B]+ = AB + BA is the anti-commutator.
As for Theorem 4.3.1, the first step is to apply this result to the Fourier transform.
So we have
Theorem 4.2.3. [80] Let f ∈ L1 ∩ L2 (Rn , Rn , dV (x)) and define Ak f (x) = xk f (x)
and Bk f (x) = 1i ∂xk f (x) then
√
kxk f k2 ξk fb ≥ 2(kf k22 + |2 hxk ∂xk f, f i|).
2
So
k∂xk f k2 = ∂d
xk f = ξk fb .
2 2
Then (4.5) becomes
1
q
kxk f k2 ξk fb ≥ kf k42 + |h2xk ∂xk f + f, f i|2
2 2q
1 2
= kf k42 + kf k22 + 2 hxk ∂xk f, f i
2q
1
≥ kf k42 + 4 |hxk ∂xk f, f i|2 .
2
q √ √
Since for a, b > 0, a +b then 2(a + b). So
2 2 a+b 1
≥ a 2 + b2 ≥
2 2 2
1
q
kxk f k2 ξk fb ≥ kf k42 + 4 |hxk ∂xk f, f i|2
2 2
√
≥ 2(kf k22 + |2 hxk ∂xk f, f i|).
In this section, we establish the main result of the thesis : a new Heisenberg
uncertainty principle for the Clifford wavelet transform. This result have been
published in [19].
Proof. As the daughter wavelet ψ a,b,s has the following Fourier expression
n
[
ψ a,b,s (ξ) = a 2 e−i<b,ξ> sψ(asξs)s,
b
we get
h i†
n
−1
Tψ [f ] (a, b, s) = a F
2 s ψ(as
b • s) sfb(•) (b)
and in the frequency domain
n
h i†
\
Tψ [f ](a, ξ, s) = a s ψ(asξs) sfb(ξ).
2 b (4.6)
Therefore by (3.8)
Z
\ 2
|ξk Tψ [f ](a, ξ, s)| dV (ξ)
Rn
Z h i†
n
= |ξk a 2 s ψ(asξs)
b sfb(ξ)|2 dV (ξ)
Rn
Z h i† † h i†
n n
= ξk a 2 s ψ(asξs) sf (ξ) ξk a 2 s ψ(asξs)
b b b sfb(ξ)dV (ξ)
Rn
Z h i† h i† †
= ξk2 an s ψ(asξs)
b sfb(ξ) s ψ(asξs)
b sfb(ξ) dV (ξ)
Rn
Z h i† h i†
= ξk2 an s ψ(asξs) sf (ξ) f (ξ) sψ(asξs)sdV
b b b b (ξ)
Rn
Z h i† h i†
= ξk2 an ψ(asξs) f (ξ) f (ξ) ψ(asξs)dV
b b b b (ξ)
Rn
Z h i† h i†
= ξk2 an ψ(asξs)
b ψ(asξs)
b fb(ξ) fb(ξ) dV (ξ). (4.7)
Rn
Using (4.7) we obtain
Z Z Z
\ 2 da
|ξk T ψ [f ](a, ξ, s)| dV (ξ) n+1 ds
a
Spin(n) R+ Rn
Z Z Z h i† i†
2 n
h da
= ξk a ψ(asξs)
b ψ(asξs)
b fb(ξ) fb(ξ) dV (ξ) n+1 ds
a
Spin(n) R+ Rn
Z Z Z h i† i†
2
h da
= ξk ψ(asξs) ψ(asξs)
b b f (ξ) f (ξ)
b b dV (ξ) ds
a
Spin(n) R+ Rn
h i†
Z Z Z ψ(asξs)
b ψ(asξs)
b
h i†
2
= dads ξk fb(ξ) fb(ξ) dV (ξ).
a
Rn +
Spin(n) R
According to (3.30), we get finally
Z Z Z
\ 2 da Aψ 2
|ξk Tψ [f ](a, ξ, s)| dV (ξ) n+1 ds = ξk f
b .
a (2π)n 2
Spin(n) R+ Rn
Z Z
da 2 da
kbk Tψ [f ] (a, •, s)k22 ds × \
ξk Tψ [f ](a, •, s) ds
Spin(n)×R+ an+1 Spin(n)×R+ 2 an+1
Z !2
1 da
≥ |Tψ [f ] (a, b, s)|2 dV (b) n+1 ds . (4.8)
2 Spin(n)×R+ ×Rn a
Now, using Lemma 4.1 and the fact that the Clifford wavelet transform is an isometry,
we get by (3.32)
Z Z Z
da
|Tψ [f ] (a, b, s)|2 dV (b)ds = Aψ kf k22 .
an+1
Spin(n) Rn R+
1
2
Z Z 1
da Aψ 2 2 1
kbk Tψ [f ] (a, •, s)k22 ds ξk fb ≥ Aψ kf k22 . (4.9)
a n+1 (2π)n 2 2
Spin(n) R+
Hence, we obtain
1
2
Z Z n
da (2π) 2 p
kbk Tψ [f ] (a, •, s)k22 n+1 ds ξk fb ≥ Aψ kf k22 .
a 2
2
Spin(n) R+
ple
In the present section we state and prove the second main result which concerns
a sharper formulation of the Clifford-wavelet uncertainty principle.
Z Z
da 1
q n o
( kbk Tψ [f ] (a, •, s)k22 n+1
ds) 2 ξk fb ≥ 2n+1 π n Aψ kf k22 + 2 |hf1 , f2 i|
Spin(n) R+ a 2
where
1 da
R R R
ψ a,b,s (x)∂bk Tψ [f ] (a, b, s) an+1
f1 (x) =
dV (b)ds
Aψ Spin(n) Rn R+
1 da
R R R
ψ a,b,s (x)bk Tψ [f ] (a, b, s) an+1
f2 (x) =
dV (b)ds
Aψ Spin(n) Rn R+
√ Z Z
da
≥ 2 (kTψ [f ] (a, •, s)k22 +|2 hbk ∂bk Tψ [f ] (a, •, s), Tψ [f ] (a, •, s)i|) ds.
Spin(n) R+ an+1
By the inequality of Cauchy-Schwartz
Z Z
da
kbk Tψ [f ] (a, •, s)k2 ξk Tψ [f\ ] (a, •, s) ds
Spin(n) R+ 2 an+1
Z Z
da 1
≤( kbk Tψ [f ] (a, •, s)k22 n+1 ds) 2
Z Spin(n)Z R+ a
2 da 1
×( ξk Tψ [f\
] (a, •, s) n+1
ds) 2 .
Spin(n) R+ 2 a
Then
Z Z Z Z
da 1 2 da 1
( kbk Tψ [f ] (a, •, s)k22 ds) × (
2 ξk Tψ [f\
] (a, •, s) ds) 2
Spin(n) R+ an+1 Spin(n) R+ 2 a n+1
√ Z Z
≥ 2 (kTψ [f ] (a, •, s)k22
Spin(n) R+
da
+ |2 hbk ∂bk Tψ [f ] (a, •, s), Tψ [f ] (a, •, s)i|) ds
an+1
√ Z Z
da
= 2 (kTψ [f ] (a, •, s)k22 ds
Spin(n) R+ an+1
√ Z Z
da
+2 2 |hbk ∂bk Tψ [f ] (a, •, s), Tψ [f ] (a, •, s)i|) ds.
Spin(n) R+ an+1
Knowing that
2 Aψ 2
da
R R
\
ξk Tψ [f ](a, •, s) n+1 ds = ξk fb
Spin(n) R+ (2π)n
2 a 2
Then we have
Z Z
da 1 Aψ 2 √
( kbk Tψ [f ] (a, •, s)k22 ds) 2 × ( ξ k
b ) 12 ≥ 2Aψ kf k2 .
f 2
Spin(n) R+ an+1 (2π)n 2
√ Z Z
da
+2 2 |hbk ∂bk Tψ [f ] (a, •, s), Tψ [f ] (a, •, s)i|) ds.
Spin(n) R+ an+1
So
Z Z
da 1
q
( kbk Tψ [f ] (a, •, s)k22 n+1
ds) 2 ξk fb ≥ 2n+1 π n Aψ kf k22
Spin(n) R+ a 2
s
2n+3 π n
Z Z
da
+ |hbk ∂bk Tψ [f ] (a, •, s), Tψ [f ] (a, •, s)i|) ds
Aψ Spin(n) R+ an+1
q
= 2n+1 π n Aψ kf k22
s
2n+3 π n
Z Z
da
+ |hbk ∂bk Tψ [f ] (a, •, s), Tψ [f ] (a, •, s)i|) ds.
Aψ Spin(n) R+ an+1
We know that
Z Z
da
|hbk ∂bk Tψ [f ] (a, •, s), Tψ [f ] (a, •, s)i|) ds
Spin(n) R+ an+1
Z Z Z
da
= [∂bk Tψ [f ] (a, b, s)]† bk Tψ [f ] (a, b, s)dV (b) ds
Spin(n) R+ Rn an+1
Z Z Z
da
≥ [∂bk Tψ [f ] (a, b, s)]† bk Tψ [f ] (a, b, s)dV (b) ds
Spin(n) R+ Rn an+1
Then we have
Z Z
da 1
q
( kbk Tψ [f ] (a, •, s)k22 ds) 2 ξk fb ≥ 2n+1 π n Aψ kf k22
Spin(n) R+ an+1 2
q
+ 2n+3 π n Aψ |[∂bk Tψ [f ] , bk Tψ [f ]]| .
we may write
Z Z
da 1
q n o
( kbk Tψ [f ] (a, •, s)k22 ds) 2 ξk fb ≥ 2n+1 π n Aψ kf k22 + 2 |hf1 , f2 i|
Spin(n) R+ an+1 2
where
1 da
R R R
ψ a,b,s (x)∂bk Tψ [f ] (a, b, s) an+1
f1 (x) =
dV (b)ds
Aψ Spin(n) Rn R+
1 da
R R R
ψ a,b,s (x)bk Tψ [f ] (a, b, s) an+1
f2 (x) =
dV (b)ds
Aψ Spin(n) Rn R+
This result is evidently more powerful than the one given in Theorem (4.3.1) (see
[20, 19]). For that, we notice that there is a slit increase in the lower bound :
q n o √ np
n o
2n+1 π n Aψ kf k22 + 2 |hf1 , f2 i| = 2(2π) 2 Aψ kf k22 + 2 |hf1 , f2 i|
np
(2π) 2 Aψ
≥ kf k22
2
[149, 151] .
4.5 Conclusion
In this chapter, we could formulate and prove a new result on Clifford wavelet
uncertainty principle stating that we can’t know simultaneously the values of the
Clifford-Fourier transform and the Clifford wavelet transform of a square integrable
multi-vector valued function. The results are based on the generalizations of the
uncertainty principle to Clifford-Fourier transform.
5.1 Conclusion
In this thesis, we answered the main problematic which is that the Clifford wavelet
transform and the Clifford-Fourier transform of multivector valued function can’t
be both sharp an uncertainty principle associated with the continuous wavelet
transform in the Clifford algebra’s settings has been formulated and proved. Start-
ing from the definition of real Clifford algebra and the real continuous wavelet
transform, we have presented a continuous Clifford wavelet transform, displayed its
properties and formulated an associated uncertainty principle. This research aimed
to state a new uncertainty principle for the Clifford wavelet transform. Based
on proprieties of Clifford algebra-valued monogenic admissible mother wavelets
and harmonic analysis of the Clifford-Fourier transform, we concluded on the im-
possibility for a Clifford wavelet transform of a function and its Clifford-Fourier
transform to be simultaneously sharply concentrated. which expresses the limita-
tions on the simultaneous concentration of Tψ [f ], and fb. This results have been
published in
5.2 Perspectives
We attend to apply the results of this thesis to some concrete Clifford wavelets such
as Clifford-Hermite, Clifford-Bessel, Clifford-Laguerre and Clifford-Gegenbauer
wavelets [28, 45, 42, 12, 11, 9, 14, 29] and other.
wavelet transform
Following [78], [2] and [113] we may try to extend the Dohono-Stark uncertainty
principle for ε-concentrated Clifford wavelet transforms. We recall that for Ω ⊂ Rn
the function f : Ω −→ Rn is ε-concentrated in the Lp norm on Ω if there exists
Ω > 0 such
Z p1
p
|f (x)| dV (x) ≤ εΩ k f k p .
Rn \Ω
As for the continuous wavelet transform, the shearlet transform [135, 60, 134, 61]
can be derived from a square-integrable group representation of a specific group
namely the shear group R∗ × Rn−1 × Rn with the operators defined the following
way : we set for a ∈ R∗ and s = (s1 , s2 , · · · , sn−1 ) ∈ Rn−1 the dilation matrix
a ... 0
. .
.. . . ...
Aa =
0 ... a
and the shear matrix
1 s1 s2 . . . sn−1
..
.
0 1 s1 s2
.. .. .. .. ..
Ss =
. . . . . .
.. ..
. . 1 s1
0 ... ... 0 1
We may try to extend this transform into the Clifford algebra Rn and deduce a new
uncertainty principle (see for example [192, 59, 206, 46, 167, 155]) and formulate
new uncertainty principles .
transform
We may find a generalised result for the Lp -variance instead of the square integrable
one and a logarithmic uncertainty principle based on [155] and [154].
2013.
[7] B. Amri and L. T. Rachdi. Uncertainty Principle in Terms of Entropy for
the Riemann–Liouville Operator. Bulletin of the Malaysian Mathematical
Sciences Society, 39(1):457–481, 2015.
[8] J.-P. Antoine, R. Murenzi, P. Vandergheynst, and S. T. Ali. Two-
Dimensional Wavelets and their Relatives. Cambridge University Press,
2004.
[9] S. Arfaoui. Some Monogenic Gegenbauer-Jacobi Polynomials and Associ-
ated Spheroidal Wavelets. PhD thesis, Faculté Des Sciences De Monastir,
Universite De Monastir, 2017.
[10] S. Arfaoui and A. Ben Mabrouk. Some Generalized Clifford-Jacobi Polyno-
mials and Associated Spheroidal Wavelets. arXiv preprint arXiv:1704.03513,
April 2017.
[11] S. Arfaoui and A. Ben Mabrouk. Some Old Orthogonal Polynomials Revis-
ited and Associated Wavelets : Two-Parameters Clifford-Jacobi Polynomi-
als and Associated Spheroidal Wavelets. Acta Applicandae Mathematicae,
155(1):177–195, December 2017.
[12] S. Arfaoui and A. Ben Mabrouk. Some Ultraspheroidal Monogenic Clifford-
Gegenbauer-Jacobi Polynomials and Associated Wavelets. Advances in Ap-
plied Clifford Algebras, 27:2287–2306, May 2017.
[13] S. Arfaoui, A. Ben Mabrouk, and C. Cattani. New Type of Gegenbauer-
Hermite Monogenic Polynomials and Associated Clifford Wavelets. Journal
of Mathematical Imaging and Vision, October 2019.
[14] S. Arfaoui, A. Ben Mabrouk, and C. Cattani. New Type of Gegenbauer-
Jacobi-Hermite Monogenic Polynomials and Associated Continuous Clifford
1882-1974), 2018.
[52] C. K. Chui. An Introduction to Wavelets. Academic Press Professional, Inc.,
San Diego, CA, USA, 1992.
[53] W. K. Clifford. Applications of Grassmann’s Extensive Algebra. American
Journal of Mathematics, 1(4):350, 1878.
[54] W. K. Clifford. On the Classification of Geometric Algebras. Mathematical
Papers, pages 397–401, 1882.
[55] J. Cnops. The Wavelet Transform in Clifford Analysis. Computational Meth-
ods and Function Theory, 1(2):353–374, November 2001.
[56] J. Cnops. An Introduction to Dirac Operators on Manifolds. Birkhäuser
Boston, 2002.
[57] J. M. Combes, A. Grossmann, Ph. Tchamitchian, and Allan D. Pierce.
Wavelets Time-Frequency Methods and Phase Space. Springer Berlin Hei-
delberg, 1989.
[58] A. K. Common and F. Sommen. Axial monogenic functions from holo-
morphic functions. Journal of Mathematical Analysis and Applications,
179(2):610–629, 1993.
[59] S. Dahlke, G. Kutyniok, P. Maass, C. Sagiv, H.-G. Stark, and G. Teschke.
The Uncertainty Principle Associated with the Continuous Shearlet Trans-
form. International Journal of Wavelets, Multiresolution and Information
Processing, 6(02):157–181, 2008.
[60] S. Dahlke, G. Steidl, and G. Teschke. The continuous shearlet transform in
arbitrary space dimensions. Journal of Fourier Analysis and Applications,
16(3):340–364, October 2009.
[61] S. Dahlke and G. Teschke. The continuous shearlet transform in higher
[120] B. Jefferies and A. Mcintosh. The Weyl Calculus and Clifford Analysis.
Bulletin of the Australian Mathematical Society, 57(2):329–341, April 1991.
[121] B. Jiang, Y. Liu, K.-I. Kou, and Z. Wang. Controllability and Observability
of Linear Quaternion-valued Systems. Acta Mathematica Sinica, 06 2020.
[122] P. Jorgensen and F. Tian. Non-commutative analysis. World Scientific, 2017.
[123] G. Kaiser. A Friendly Guide to Wavelets. Birkhauser Boston Inc., Cam-
bridge, MA, USA, 1994.
[124] M. Kesri. Solutions globales du système intégro-différentiel non linéaire de
Liouville–Boltzmann–Poisson et l'équation spectrale. Colloquium Mathe-
maticum, 106(2):265–281, 2006.
[125] M. Kesri. An Explicit Method for Optimal Control Problems. Applied Math-
ematics and Computation, 219(15):8213–8221, apr 2013.
[126] M. Kesri. On the Class of Infinite Horizon Optimal Control Problems Affine,
Which Are Structurally Stable. Mathematical Methods in the Applied Sci-
ences, 39(11):2799–2807, oct 2015.
[127] M. Kesri. Continuous Dependence of the Value Function on a Class of One-
Dimensional Autonomous Optimal Control Problems and Existence of Solu-
tions. Indian Journal of Pure and Applied Mathematics, 50(1):269–281, mar
2019.
[128] A. Kheirabadi, A. M. Vaziri, and S. Effati. Solving Optimal Control Prob-
lem Using Hermite Wavelet. Numerical Algebra, Control & Optimization,
9(1):101–112, 2019.
[129] J. King. Uncertainty Principle in Control Theory, Part I: Analysis of Per-
formance Limitations. ArXiv pre-print : 1402.5468v1.
[130] V. Kisil. Wavelets in Banach Spaces. Acta Applicandae Mathematicae,
59(1):79–109, 1999.
[131] K. I. Kou, J.-Y. Ou, and J. Morais. On uncertainty principle for quaternionic
linear canonical transform. Abstract and Applied Analysis, 2013:1–15, 2013.
[132] K.-I. Kou and T. Qian. The Paley-Wiener Theorem in Rn with the Clifford
Analysis Setting. Journal of Functional Analysis, 189(1):227–241, feb 2002.
[133] S. V. Kovalevskaja. Zur Theorie Der Partiellen Differentialgleichungen. J.
Fur Die Reine Und Angew. Mathem., 80:1–32, 1875.
[134] D. Kumar. Multi dimensional continuous shearlet transform and multi scale
wavelet shrinkage with integro differential equations. American Review of
Mathematics and Statistics, 4(2):39–47, December 2016.
[135] G. Kutyniok and D. Labate, editors. Shearlets. Birkhäuser Boston, 2012.
[136] P. Van Lancker. Clifford analysis on the sphere. In Clifford Algebras and
Their Application in Mathematical Physics, pages 201–215. Springer Nether-
lands, 1998.
[137] G. Laville. Some Topics in Clifford Analysis. Advances in Applied Clifford
Algebras, 19(3-4):721–775, nov 2009.
[138] J.-W. Lee, L.-W. Liu, H.-K. Hong, and J.-T. Chen. Applications of the Clif-
ford Algebra Valued Boundary Element Method to Electromagnetic Scatter-
ing Problems. Engineering Analysis with Boundary Elements, 71:140–150,
2016.
[139] Z. Li, H. Yang, and Y. Qiao. A New Cauchy Integral Formula in the Complex
Clifford Analysis. Advances in Applied Clifford Algebras, 28(4), jul 2018.
[140] P. Lounesto. Clifford Algebras and Spinors. Cambridge University Press,
2001.
[141] W. Luo and J. Du. Generalized Cauchy Theorem in Clifford Analysis and
2013.
[189] M. Sce. Osservazioni Sulle Serie Di Potenze Nei Moduli Quadratici. Atti
Accad. Naz.Lincei. Rend. Cl. Sci. Fis. Mat. Nat., 8(23):220–225, 1957.
[190] K. K. Selig. Uncertainty Principles Revisited. Electron. Trans. Numer. Anal.,
14:165–177, 2002.
[191] D. Sen. The Uncertainty Relations in Quantum Mechanics. Current Science,
107(2):203–2018, July 2014.
[192] J. Sharma and S. K. Singh. Clifford Valued Shearlet Transform. Advances
in Applied Clifford Algebras, 30(3), jun 2020.
[193] S. Shirrell and R. Walter. Hermitian Clifford Analysis and Its Connec-
tions with Representation Theory. Complex Variables and Elliptic Equations,
62(9):1329–1342, March 2017.
[194] H. Skibbe and M. Reisert. Spherical tensor algebra: A toolkit for 3d image
processing. J. Math. Imaging Vis., 58(3):349–381, 2017.
[195] I. N. Sneddon. Fourier Transforms. Dover Publications Inc., 2003.
[196] F. Sommen. Special Functions in Clifford Analysis and Axial Symmetry.
Journal of Mathematical Analysis and Applications, 130(1):110–133, Febru-
ary 1988.
[197] F. Sommen. Clifford analysis and integral geometry. In Clifford Algebras
and their Applications in Mathematical Physics, pages 293–311. Springer
Netherlands, 1992.
[198] F. Sommen. Clifford Analysis in Two and Several Vector Variables. Appli-
cable Analysis, 73(1-2):225–253, oct 1999.
[199] F. Sommen and H. De Schepper. Introductory Clifford Analysis. Operator
Theory, pages 1–27, 2015.