Abstracts
One reminds for all function f : <img border=0 width=18 height=19 src="../../../../../img/revistas/tema/v15n3/r_bastao.jpg">n → <img border=0 width=17 height=21 src="../../../../../img/revistas/tema/v15n3/r_bastao_barra.jpg"> the so-called (min, +)-wavelets which are lower and upper hulls build from (min, +) analysis [12, 13]. One shows that this analysis can be applied numerically to the Weierstrass and Weierstrass-Mandelbrot functions, and that one recovers their theoretical Hölder exponents and fractal dimensions.
(min, +)-wavelets; Hölder exponents; Weierstrass functions
Lembrando que para todas as funções f : <img border=0 width=21 height=21 src="../../../../../img/revistas/tema/v15n3/r_bastao.jpg">n → <img border=0 width=19 height=23 src="../../../../../img/revistas/tema/v15n3/r_bastao_barra.jpg">, as chamadas (min, +)-wavelets são construções do fecho inferior e superior, vindos da análise (min, +) [12, 13]. Mostra-se que esta análise pode ser aplicada numericamente às funções de Weierstrass e Weierstrass-Mandelbrot, e que recupera-se os seus expoentes de Hölder teóricos e dimensões fractais.
(min, +)-wavelets; expoentes de Hölder; funções de Weierstrass
Numerical calculations of Hölder exponents for the Weierstrass functions with (min, +)-wavelets
M. GondranI; A. KenoufiII, * * Corresponding author: Abdel Kenoufi
IUniversity Paris Dauphine, Lamsade, Paris, France. E-mail: michel.gondran@polytechnique.org
IIScientific COnsulting for Research & Engineering (SCORE), Strasbourg, France. E-mail: kenoufi@s-core.fr
ABSTRACT
One reminds for all function f : n → the so-called (min, +)-wavelets which are lower and upper hulls build from (min, +) analysis [12, 13]. One shows that this analysis can be applied numerically to the Weierstrass and Weierstrass-Mandelbrot functions, and that one recovers their theoretical Hölder exponents and fractal dimensions.
Keywords: (min, +)-wavelets, Hölder exponents,Weierstrass functions.
RESUMO
Lembrando que para todas as funções f : n → , as chamadas (min, +)-wavelets são construções do fecho inferior e superior, vindos da análise (min, +) [12, 13]. Mostra-se que esta análise pode ser aplicada numericamente às funções de Weierstrass e Weierstrass-Mandelbrot, e que recupera-se os seus expoentes de Hölder teóricos e dimensões fractais.
Palavras-chave: (min, +)-wavelets , expoentes de Hölder, funções de Weierstrass.
1 INTRODUCTION
Genesis of wavelets theory started in 1946 with D. Gabor [9], who introduced the Windowed Fourier Transform (WFT)
for the local spectral analysis of radar signals. The localization is reached due to fast decaying window function g(x) 0. Even if WFT exhibits many powerful and practical features, there are some defects compared to Fourier Transform. The transform (1.1) can not resolve efficiently wavelengths longer than the window g(x) width. Conversely, for signal with high frequencies, short decomposition needs a broad window with a large number of periods. Thus, signal reconstruction in this case adds a large number of terms with comparable amplitudes and hence becomes numerically unstable. Finally, one needs a scheme with a wide window for low frequency signals and a narrow window for high frequency ones. Such a scheme, was independently suggested as a tool for geophysical studies by several authors at the beginning of 1980s [19, 23]. Wavelets Theory (WT) word was introduced in analysis by J. Morlet [17, 11]. It is considered nowadays as a preferable alternative to the Fourier analysis, used where and when the signals are random and comprised of fluctuations of different scales, such as in turbulence phenomena [3]. WT has been immediately followed by several applications in science and engineering, such as signal processing and detection, fractals, self-similar objects, self-similar random processes, like turbulence and Brownian motion [2]. WT was then mathematically formalized by Grossman & Morlet [11], Goupillaud et al. [10], Daubechies [6] and some other authors. Practically, WT is a separate convolution of the signal in question with a family of functions obtained from some basic one, the basic wavelet called mother wavelet or analysing function, by shifts τ and dilatations a:
An interesting point of view is to consider WT as a realization of a function decomposition with respect to the representation of the affine group [2]: x
ax + b.We refer the reader to useful and interesting articles [1, 8, 20, 4, 7, 5] and both theoretical and practical book [2].
Nevertheless, wavelets decompositions are limited by their linear features. This present article aims to apply for pathological functions such as Weierstrass functions, a non-linear transform, called (min, +) transform, which has been already defined within (min, +) analysis [12, 13, 18]. This one consists to replace in the scalar product definition of two real-valued functions f and g defined on a domain X, the real number field (, +, ×) with the (min, +) dioid (∪{+ ∞}, min, +). The classical scalar product 〈f, g〉 = ∫x∈Xf(x)g(x)dx becomes then the (min, +) scalar product [12]:
The demonstration that it is a scalar product within the (min, +) dioid is straightforward and easy excepted for its linearity.
One has to show that 〈f, g〉min + is distributive according to min, which means
〈f, min{g1, g2)〉(min,+) = min{〈f, g1〉(min,+), 〈f, g2〉(min,+)},
and linear according to the addition of a scalar λ : 〈f (x), λ + g(x)〉(min,+) = λ + 〈f, g〉(min,+). The linearity is obvious since infx∈X{f (x) + λ + g(x)} = λ + infx∈X{f(x) + g(x)}. Distributivity is obtained in two steps. One has first to prove this equality with mean of two inequalities. We start first with the simple relations:
〈f, g1〉(min,+)< f(x) + g1(x), and 〈f, g2〉(min,+)< f(x) + g2(x), ∀x.
This gives min{〈f, g1〉(min,+), 〈f, g2〉(min,+)} < min{f(x) + g1(x), f(x) + g2(x)}∀x. And since
min{f(x) + g1(x), f(x) + g2(x)} = f (x) + min{g1(x), g2(x)},
one has min{〈f, g1〉(min,+), 〈f, g2〉(min,+)} < f(x) + min{g1(x), g2(x)}∀x, which yields to the inequality
In a second step, one can write
〈f, min{g1, g2}〉(min,+)< f(x) + min{g1(x), g2(x)} < f(x) + g1(x) ∀x,
which becomes
and in the same manner
〈f, min{g1, g2}〉(min,+)< f(x) + min{g1(x), g2(x)} < f(x) + g2(x) ∀x,
giving now
and then from (1.4) and (1.5)
From relations (1.3) and (1.6), one deduces finally the equality and thus the distributivity.
With this (min, +) scalar product, one obtains a distribution-like theory : the operator is linear and continuous according the dioid structure (∪{+ ∞}, min, +), non-linear and continuous according to the classical structure (, +, ×). The non-linear distribution δmin(x) defined as
δ(min,+)(x) = {0 if x = 0, +∞ else}
is similar in (min, +) analysis to the classical Dirac distribution. Then, one has
In (min, +) analysis, the Legendre-Fenchel transform which permits to get Hamiltonian from Lagrangian and which has an important role in physics is similar in (min, +) analysis to the Fourier transform in the classical one [18].
In this article, we explore how (min, +)-wavelets decomposition and reconstruction could be an interesting signal processing tool, since (min, +) transforms can be applied to a larger class of functions than the functions treated with classical wavelet transforms, especially to lower semi-continuous functions [12], such as x g(x) · Floor(x) for instance, where g is a continuous function.
In this paper, one focus on Weierstrass and Weierstrass-Mandelbrot functions which are classical examples of functions continuous everywhere but differentiable nowhere [22].
One introduces in the following Section 2, the (min, +)-wavelets decomposition and reconstruction of a signal with mean of (min, +) scalar product. in Section 3, we show how the (min, +)wavelets allow a characterisation of Hölder functions. in Section 4, we apply these results to numerical calculations of Hölder exponents of Weierstrass-like functions and compare them to the theoretical values [21]. This permits to deduce immediately their fractal dimensions.
2 (min, +)-WAVELETS
The usual wavelet transform of a function f from n to is a linear transform defined for all scales a ∈+ and points b ∈n, which can be computed according to the equation (1.2):
where Ψ is a function called mother wavelet or analysing function. It has to be zero average and exhibiting oscillations until a certain order p. This can be written as
In (min, +) analysis, a set of non-linear transforms has been introduced for lower semi-continuous functions [12, 18], the so-called (min, +)-wavelets transforms which are defined for a function f : n → and for all a ∈+ and b ∈n such as:
where h is a basis analysing function (upper semi-continuous and inf-compact verifying h(0) = 0), like the following functions:
Since Tf (a, x) < f (x) for all a > 0, Tf (a, x) is a lower hull of f(x). For any lower bounded and lower semi-continuous function, one has a reconstruction formula like in the linear wavelets theory [2]:
which can be simplified within the (min, +) theory in
The (min, +)-wavelets analysis will be based on simultaneous analysis of lower hulls Tf (a, b), and upper hulls of f represented by Tf +(a, b) defined by:
For the upper hulls Tf +(a, b), we have a reconstruction formula which is symmetric to lower hulls Tf (a, b) (2.9, 2.8):
which simplifies as well as:
For each analysing function h, one has [13]:
because Tf(a, x) (respectively Tf +(a, x)) are functions decreasing with scales (respectively increasing) and converging to f*(x) (respectively f *(x)), the lower semi-continuous closure of f (respectively upper semi-continuous closure) when the scale tends to 0.
Remark 1. We use the word "wavelet" by analogy with linear wavelets since the decomposition and reconstruction formula are very similar and since one just replaces the usual real number field (, +, ×) with the (min, +) dioid (∪{+ ∞}, min, +). Another name can be (min, +) pen or (min, +) hulls.
Remark 2. The shift and scale parameters have the same meaning as in Linear Wavelet Theory: for high frequencies, one needs small scales, and the inverse as well. But the relation between them is not simply proportionally inverse as in linear theory, because it depends on the choice of analyzing function hα, and this introduces non-linear dependency between scale and frequency. This leads to a relation such as ν = γ(a, α), where ν, and a are respectively the frequency and the scale, and γ a non-linear function decreasing with a.
Definition 1. (min, +)-wavelet is defined as the couple {Tf (a, x), Tf +(a, x)}. For all
+,the a-oscillation of f is defined:
In the case of analysing function h∞, one has
and ΔTf (a, x) = sup|xy|<a f (y) inf|xz|<af (z) corresponds to the a-oscillation defined in one dimension by Tricot [21]: oscaf (x) = supy,z∈[xa,x+a][ f (y) f (z)].
3 CHARACTERISATION OF HÖLDERIAN FUNCTIONS WITH (min, +)-WAVELETS ANALYSIS
The calculations of oscillations according to the analysing function and the scale will permit to study the global and local regularity of a function.
First, let's start with the case of global regularity of a Hölderian function for which it exists H (0 < H < 1) and a constant K such as
It is a sufficient but not necessary condition for a function to be continuous. In the case of fractal function, K is related to its fractal dimension.
Theorem 1. The function f is Hölderian with exponent H, 0 < H < 1, if and only if it exists a constant C such as for all a, one of the following condition is verified:
Demonstration:
-
Demonstration for the case of analysing function h∞ is classic [21]: if f verifies (3.17), let's consider some x and y in n. One can assume that f (x) > f (y). Then, one has
this yields to
|f (x) f (y)| < ΔTf(a, x) < KaH < K|x y|H .
Conversely, let's assume that |f (x) f (y)| < K |x y|H for all y. Let y1 such as f (y1) =sup|xz|<af (z) and y2 such as f (y2) =inf|xz|<af (z). One has then
ΔTf(a, x) = f (y1) f (y2) = f (y1) f (y) + f (y) f (y2),
which yields to
ΔTf(a, x) < |f (y1) f (y)| + |f (y) f (y2)| < 2KaH .
In the case of analysing functions hα, α > 1, let's suppose first that f verifies (3.18). We consider x and y in n with f (x) > f (y). Reconstruction equation (2.12) of f (x) can be written as
and the equation of reconstuction for f (y) (2.8)
One deduces
thus
and the optimisation on the scale a implies that it exists K such as |f (x) f (y)| < K |x y|H .
Conversely, let's assume that |f (x) f (y)| < K |x y|H for all x and y. Using (2.8) and (2.11), one has
We deduce that
whose optimisation gives (3.18). □
Let's consider now the case of local irregularity at x0 where the function is Hölderian: it exists H (0 < H < 1) and a constant K such as
Theorem 2. The function f is Hölderian at point x0, with exponent H, 0 < H < 1, if and only if it exists a constant C such as for all a, one has one of the following conditions:
Demonstration:
-
In the case of the analysing function h∞, if f verifies (3.20) for all x, let a = |x x0|. One has then inequations sup|xx0|<af (z) > f (x) > f (x0) > inf|xz| < af (z) or sup|xx0 |<af (z) > f (x0) > f (x) >inf|xz|<af (z). In both cases one gets
|f (x) f (x0)| < ΔTf(a, x) < 2C|x x0|H .
Conversely, let suppose |f (x) f (x0)| < K |x x0|H for all x, and y1 such as f (y1) = sup|xz|<af (z) and y2 such as f (y2) =inf|xz|<af (z); One has then
that means
ΔTf(a, x) < K(|y1 x|H + |x x0|H + |y2 x|H + |x x0|H ),
this yields to
ΔTf(a, x) < 2K(aH + |x x0|H).
-
For analysing functions hα, α > 1, we assume first that f verifies (3.21). Let's consider x in n and the two cases, f (x) > f (x0),and f (x) < f (x0). In the first case, one uses the reconstruction equations
and
For the second case, one uses a symmetric reconstruction method. This yields to
which gives
This implies that it exists a constant K such as
|f (x) f (x0)| < K|x x0|H .
Conversely, let suppose that |f (x) f (x0)| < K|x x0|H for all x. With mean of (2.8) and (2.11), one has
Since
f (x) f (y) = f (x) f (x0) + f (x0) f (y),
one deduces
which yields to
whose optimisation gives (3.21). □
One gets here a reciprocal relation which is not fully obtained with linear wavelets [15].
4 HÖLDER EXPONENTS CALCULATION FOR WEIERSTRASS FUNCTIONS
We exhibit an application of the (min, +)-wavelets analysis to the Weierstrass function in order to compute its Hölder exponent H and its fractal dimension D. This one is a typical example of function continuous everywhere but nowhere differentiable [22]. One consider the general form of Weierstrass functions on [0, 2π]
with ωH > 1and {φm}m>0, constant or randomly distributed variable.
Those functions are Hölderian (and anti-Hölderian) with coefficient H and fractal dimension [16, 14]:
One calculates for all scales s = k · scalemin with k an integer from 1 to 10 and scalemin =102, the following function of scale for h2 and h∞
For the Weierstrass function, the upper bound of the sum is replaced with a finite constant M =15 which is sufficient for our tests. Thus, the truncated Weierstrass function can be written as
and is represented with its (min, +)-wavelets decomposition on Figure 1 for φm= 0.
We made numerical calculations to determine Hölder exponents. The fractal dimension is then directly given by equation (4.25). Computations were performed for H ∈ , ω =2, for analysing functions h∞ and h2, for both cases of zero and random φmwith a uniform probability measure in [0, 2π].
The slope of the linear part of curves for small scales gives the value of Hölder exponent.
Hölder exponents calculations for random phase Weierstrass functions are summarized on Tables 1 and 2. According to equations (3.17, 3.18), the slopes and Hölder exponents are very close to the theoretical value H = for h∞ and for h2 [21, 14]. The fractal dimension is then given by D = 2 H = . Same result for H = with a slope of H = for h∞ and for h2. They confirm that the Hölder exponents and fractal dimensions of Weierstrass function remain the same in the case of a uniform random phase [16, 14].
The Weierstrass function W(t) = (ωH)m cos(ωmt), verifies
which is not a scaling invariance property [21]. In order to circumvent that, one can build the Weierstrass-Mandelbrot function
Since
the change of variable m' = m + 1 leads to WM(ωt) = ωHWM(t), which has scaling invariance property. Hölder exponents calculations for a truncated version of this function are exhibited on Figures 2 and 3, confirming thus the validity of (min, +)-wavelets decomposition for its Hölder exponents computation.
5 CONCLUSION AND PERSPECTIVES
We have presented in this article a promising tool to determine numerically Hölder exponents of Weierstrass-like functions which are exhibiting fractal properties. It is based on (min, +) analysis and proposes a signal decomposition using the (min, +) scalar product. By analogy with Linear Wavelet Theory, this permits to define (min, +)-wavelets , which are lower and upper hulls of a signal at a certain scale.
ACKNOWLEDGEMENTS
We thank Mikhail Altaisky, Thierry Lehner and René Voltz for useful and helpful discussions about turbulence and wavelets.
Received on November 9, 2013
Accepted on August 29, 2014
- [1] A. Arneodo, F. Argoul, E. Bacry, J. Elezgaray & J.-F. Muzy. Ondelettes, multifractales et turbulence Diderot, Paris, (1995).
- [2] M.V. Altaisky. Wavelets, Theory, Applications and Implementation Universities Press, (2005).
- [3] M.V. Altaisky. Multiscale theory of turbulence in wavelet representation. Number 410 (3). Doklady Akademii Nauk, (2006).
- [4] N.M. Astaf'eva. Wavelet analysis: basic theory and some applications. Uspekhi fizicheskih nauk, 166(11) (1994), 11461170, (in Russian).
- [5] C.K. Chui. An Introduction to Wavelets Academic Press Inc., (1992).
- [6] I. Daubechies. Orthonormal bases of compactly supported wavelets. Comm. Pure. Apl. Math., 41 (1988), 909996.
- [7] I. Daubechies. Ten lectures on wavelets S.I.A.M., Philadelphie, (1992).
- [8] M. Farge. Wavelets and their application to turbulence. Annual review of fluid mechanics, 24 (1992), 395457.
- [9] D. Gabor. Theory of communication. Proc. IEE, 93 (1946), 429457.
- [10] P. Goupillaud, A. Grossmann & J. Morlet. Cycle-octave and related transforms in seismic signal analysis. Geoexploration, 23 (1984/85), 85102.
- [11] A. Grossmann & J. Morlet. Decomposition of Hardy functions into square-integrable wavelets of constant shape. SIAM J. Math. Anal., 15(4) (1984), 723736.
- [12] M. Gondran. Analyse minplus. C.R. Acad. Sci. Paris, 323 (1996), 371375.
- [13] M. Gondran. Convergences de fonctions à valeurs dans Rk et analyse minplus complexe. C.R. Acad. Sci. Paris, 329 (1999), 783788.
- [14] B.R. Hunt. The hausdorff dimension of graphs of weierstrass functions. Proceedings of the American Mathematical Society, 126 (1998), 791800.
- [15] S. Jaffard. Décompositions en ondelettes, "developments of mathematics 1950-2000" (J-P. Pier ed.). pages 609634, (2000).
- [16] J. Mallet-Paret J.L. Kaplan & J.A. Yorke. The lyapounov dimension of nowhere differentianle attracting torus. Ergod. Th. Dynam. Sys., 4 (1984), 261281.
- [17] J. Morlet, G. Arens & I. Fourgeau. Wave propagation and sampling theory. Geophysics, 47 (1982), 203226.
- [18] M. Minoux & M. Gondran. Graphs, Dioids and Semirings Springer, (2008).
- [19] J. Morlet. Sampling theory and wave propagation. In: Proc. 51st Annu. Meet. Soc. Explor. Geophys., Los-Angeles, (1981).
- [20] K. Shneider & M. Farge. Wavelet approach for modelling and computing turbulence, volume 1998-05 of Advances in turbulence modelling. Von Karman Institute for Fluid Dynamics, Bruxelles, (1998).
- [21] C. Tricot. Courbes et dimension fractale Springer-Verlag, (1993).
- [22] K. Weierstrass. Über continuirliche funktionen eines reellen arguments, die f¨ur keinen werth des letzteren einen bestimmten differentialquotienten besitzen. Karl Weiertrass Mathematische Werke, Abhandlungen II, Gelesen in der Königl. Akademie der Wissenchaften am 18 Juli 1872, 1967.
- [23] V.D. Zimin. Hierarchic turbulence model. Izv. Atmos. Ocean. Phys., 17 (1981), 941946.
Publication Dates
-
Publication in this collection
09 Feb 2015 -
Date of issue
Dec 2014
History
-
Received
09 Nov 2013 -
Accepted
29 Aug 2014