Section: Research Program

Stochastic models

The second axis in the theoretical developments of the Regularity team aims at defining and studying stochastic processes for which various aspects of the local regularity may be prescribed.

Multifractional Brownian motion

One of the simplest stochastic process for which some kind of control over the Hölder exponents is possible is probably fractional Brownian motion (fBm). This process was defined by Kolmogorov and further studied by Mandelbrot and Van Ness, followed by many authors. The so-called “moving average” definition of fBm reads as follows:

Y t = - 0 ( t - u ) H - 1 2 - ( - u ) H - 1 2 . 𝕎 ( d u ) + 0 t ( t - u ) H - 1 2 . 𝕎 ( d u ) ,

where 𝕎 denotes the real white noise. The parameter H ranges in (0,1), and it governs the pointwise regularity: indeed, almost surely, at each point, both the local and pointwise Hölder exponents are equal to H.

Although varying H yields processes with different regularity, the fact that the exponents are constant along any single path is often a major drawback for the modeling of real world phenomena. For instance, fBm has often been used for the synthesis natural terrains. This is not satisfactory since it yields images lacking crucial features of real mountains, where some parts are smoother than others, due, for instance, to erosion.

It is possible to generalize fBm to obtain a Gaussian process for which the pointwise Hölder exponent may be tuned at each point: the multifractional Brownian motion (mBm) is such an extension, obtained by substituting the constant parameter H(0,1) with a regularity function H:+(0,1).

mBm was introduced independently by two groups of authors: on the one hand, Peltier and Levy-Vehel [29] defined the mBm {Xt;t+} from the moving average definition of the fractional Brownian motion, and set:

X t = - 0 ( t - u ) H ( t ) - 1 2 - ( - u ) H ( t ) - 1 2 . 𝕎 ( d u ) + 0 t ( t - u ) H ( t ) - 1 2 . 𝕎 ( d u ) ,

On the other hand, Benassi, Jaffard and Roux [45] defined the mBm from the harmonizable representation of the fBm, i.e.:

X t = e i t ξ - 1 ξ H ( t ) + 1 2 . 𝕎 ^ ( d ξ ) ,

where 𝕎^ denotes the complex white noise.

The Hölder exponents of the mBm are prescribed almost surely: the pointwise Hölder exponent is αX(t)=H(t)αH(t) a.s., and the local Hölder exponent is α˜X(t)=H(t)α˜H(t) a.s. Consequently, the regularity of the sample paths of the mBm are determined by the function H or by its regularity. The multifractional Brownian motion is our prime example of a stochastic process with prescribed local regularity.

The fact that the local regularity of mBm may be tuned via a functional parameter has made it a useful model in various areas such as finance, biomedicine, geophysics, image analysis, .... A large number of studies have been devoted worldwide to its mathematical properties, including in particular its local time. In addition, there is now a rather strong body of work dealing the estimation of its functional parameter, i.e. its local regularity. See http://regularity.saclay.inria.fr/theory/stochasticmodels/bibliombm for a partial list of works, applied or theoretical, that deal with mBm.

Self-regulating processes

We have recently introduced another class of stochastic models, inspired by mBm, but where the local regularity, instead of being tuned “exogenously”, is a function of the amplitude. In other words, at each point t, the Hölder exponent of the process X verifies almost surely αX(t)=g(X(t)), where g is a fixed deterministic function verifying certain conditions. A process satisfying such an equation is generically termed a self-regulating process (SRP). The particular process obtained by adapting adequately mBm is called the self-regulating multifractional process [3] . Another instance is given by modifying the Lévy construction of Brownian motion [4] . The motivation for introducing self-regulating processes is based on the following general fact: in nature, the local regularity of a phenomenon is often related to its amplitude. An intuitive example is provided by natural terrains: in young mountains, regions at higher altitudes are typically more irregular than regions at lower altitudes. We have verified this fact experimentally on several digital elevation models [8] . Other natural phenomena displaying a relation between amplitude and exponent include temperatures records and RR intervals extracted from ECG [9] .

To build the SRMP, one starts from a field of fractional Brownian motions B(t,H), where (t,H) span [0,1]×[a,b] and 0<a<b<1. For each fixed H, B(t,H) is a fractional Brownian motion with exponent H. Denote:

X ¯ ̲ α ' β ' = α ' + β ' - α ' X - min K X max K X - min K X

the affine rescaling between α' and β' of an arbitrary continuous random field over a compact set K. One considers the following (stochastic) operator, defined almost surely:

Λ α ' , β ' : 𝒞 0 , 1 , α , β 𝒞 0 , 1 , α , β Z ( . ) B ( . , g Z ( . ) ̲ ¯ α ' β '

where αα'<β'β, α and β are two real numbers, and α',β' are random variables adequately chosen. One may show that this operator is contractive with respect to the sup-norm. Its unique fixed point is the SRMP. Additional arguments allow to prove that, indeed, the Hölder exponent at each point is almost surely g(t).

An example of a two dimensional SRMP with function g(x)=1-x2 is displayed on figure 1 .

Figure 1. Self-regulating miltifractional process with g(x)=1-x2

We believe that SRP open a whole new and very promising area of research.

Multistable processes

Non-continuous phenomena are commonly encountered in real-world applications, e.g. financial records or EEG traces. For such processes, the information brought by the Hölder exponent must be supplemented by some measure of the density and size of jumps. Stochastic processes with jumps, and in particular Lévy processes, are currently an active area of research.

The simplest class of non-continuous Lévy processes is maybe the one of stable processes [56] . These are mainly characterized by a parameter α(0,2], the stability index (α=2 corresponds to the Gaussian case, that we do not consider here). This index measures in some precise sense the intensity of jumps. Paths of stable processes with α close to 2 tend to display “small jumps”, while, when α is near 0, their aspect is governed by large ones.

In line with our quest for the characterization and modeling of various notions of local regularity, we have defined multistable processes. These are processes which are “locally” stable, but where the stability index α is now a function of time. This allows to model phenomena which, at times, are “almost continuous”, and at others display large discontinuities. Such a behaviour is for instance obvious on almost any sufficiently long financial record.

More formally, a multistable process is a process which is, at each time u, tangent to a stable process [49] . Recall that a process Y is said to be tangent at u to the process Yu' if:

lim r 0 Y ( u + r t ) - Y ( u ) r h = Y u ' ( t ) , (4)

where the limit is understood either in finite dimensional distributions or in the stronger sense of distributions. Note Yu' may and in general will vary with u.

One approach to defining multistable processes is similar to the one developed for constructing mBm [29] : we consider fields of stochastic processes X(t,u), where t is time and u is an independent parameter that controls the variation of α. We then consider a “diagonal” process Y(t)=X(t,t), which will be, under certain conditions, “tangent” at each point t to a process tX(t,u).

A particular class of multistable processes, termed “linear multistable multifractional motions” (lmmm) takes the following form [11] , [10] . Let (E,,m) be a σ-finite measure space, and Π be a Poisson process on E× with mean measure m× ( denotes the Lebesgue measure). An lmmm is defined as:

Y ( t ) = a ( t ) ( 𝖷 , 𝖸 ) Π 𝖸 < - 1 / α ( t ) > | t - 𝖷 | h ( t ) - 1 / α ( t ) - | 𝖷 | h ( t ) - 1 / α ( t ) ( t ) . (5)

where x<y>:=sign(x)|x|y, a:+ is a C1 function and α:(0,2) and h:(0,1) are C2 functions.

In fact, lmmm are somewhat more general than said above: indeed, the couple (h,α) allows to prescribe at each point, under certain conditions, both the pointwise Hölder exponent and the local intensity of jumps. In this sense, they generalize both the mBm and the linear multifractional stable motion [57] . From a broader perspective, such multistable multifractional processes are expected to provide relevant models for TCP traces, financial logs, EEG and other phenomena displaying time-varying regularity both in terms of Hölder exponents and discontinuity structure.

Figure 2 displays a graph of an lmmm with linearly increasing α and linearly decreasing H. One sees that the path has large jumps at the beginning, and almost no jumps at the end. Conversely, it is smooth (between jumps) at the beginning, but becomes jaggier and jaggier as time evolves.

Figure 2. Linear multistable multifractional motion with linearly increasing α and linearly decreasing H