From N parameter fractional Brownian motions to N parameter multifractional Brownian motions, , November 2002

one is isotropic while the other is not. For each of these processes, a moving average representation, a harmonizable representation, and the covariance structure are given.

The Hölder regularity is then studied. In particular, the case of an irregular exponent function

H

is investigated. In this situation, the almost sure pointwise and local Hölder exponents of the multi-parameter mBm are proved to be equal to the correspondent exponents of

H

. Eventually, a local asymptotic self-similarity property is proved. The limit process can be another process than fBm.

AMS classification: 62 G 05, 60 G 15, 60 G 17, 60 G 18.

Keywords: fractional Brownian motion, Gaussian processes, Hölder regularity, local asymptotic self-similarity, multi-parameter processes.

1 Introduction

In many applications, fractional Brownian motion (fBm) seems to fit very well to random phenomena. Recall that it can be defined by one of the four following properties. Let

H \in (0, 1)

(

H

is sometimes called the Hurst parameter).

∙ $B^{H}$ is a centered Gaussian process such that $\forall s, t \in R_{+}; E [B_{s}^{H} B_{t}^{H}] = \frac{1}{2} [s^{2 H} + t^{2 H} - | t - s |^{2 H}]$
∙ the process $B^{H}$ such that $\forall t \in R_{+}; B_{t}^{H} = \int_{- \infty}^{0} [(t - u)^{H - \frac{1}{2}} - (- u)^{H - \frac{1}{2}}] . W (d u) + \int_{0}^{t} (t - u)^{H - \frac{1}{2}} . W (d u)$ is a fBm,
∙ the process $B^{H}$ such that $\forall t \in R_{+}; B_{t}^{H} = \int_{R} \frac{e^{i t ξ} - 1}{{| ξ |}^{H + \frac{1}{2}}} . \hat{W} (d ξ)$ is a fBm,
∙ $B^{H}$ is the unique self-similar Gaussian process with stationary increments.

Its efficiency has already been shown in simulation of traffic on Internet or in finance. This induced some recent progress such as stochastic integration against fBm.

However, the main limitation of fBm is that the Hölder regularity is constant along the paths.

Multifractional Brownian motion (mBm) has been independently introduced in [4] and [13] . This process is a generalization of fractional Brownian motion where the Hurst parameter

H

is substituted by a function

t \mapsto H (t)

. As a consequence the Hölder exponent is allowed to vary along trajectories.

The different definitions by the two groups of authors provided two different representations of mBm.

Peltier and Levy-Vehel ([13] ) defined the mBm from the moving average definition of the fractional Brownian motion

X_{t} = \int_{- \infty}^{0} [(t - u)^{H (t) - \frac{1}{2}} - (- u)^{H (t) - \frac{1}{2}}] . W (d u) + \int_{0}^{t} (t - u)^{H (t) - \frac{1}{2}} . W (d u)

where

t \mapsto H (t)

is a Hölder function.

Benassi, Jaffard and Roux ([4] ) defined the mBm from the harmonizable representation of the fBm

X_{t} = \int_{R} \frac{e^{i t ξ} - 1}{{| ξ |}^{H (t) + \frac{1}{2}}} . \hat{W} (d ξ)

These two definitions were proved to be equivalent up to a multiplicative deterministic function ([6] ).

Moreover, in [3] the covariance function of this Gaussian process has been proved to be

E [X_{s} X_{t}] = D (H (s), H (t)) [| s |^{H (s) + H (t)} + | t |^{H (s) + H (t)} - | t - s |^{H (s) + H (t)}]

where

D

is a known deterministic function.

The goal of this paper is to study some multi-parameter extension of the multifractional Brownian motion, ie a stochastic process indexed by

R_{+}^{N}

, which is an mBm when

N = 1

. One extension has already been considered in [4] .

2D extension of fractional Brownian motion has been already used in various applications such as underwater terrain modeling ([14] ). It may be more realistic to allow local regularity to vary at each point : our extension of mBm in

R^{2}

may be used for this kind of application.

2 Multi-parameter extension of the fractional Brownian motion

Since multifractional Brownian motion is an extension of fractional Brownian motion, we start with a review of the existing extensions of fBm. Most of the results in this section are well-known, but we give new proofs based only on the covariance functions.

In the same way as Brownian motion has two main multi-parameter extensions:

Levy Brownian motion and Brownian sheet, two different multi-parameter extensions of fractional Brownian motion have been defined.

2.1 Levy fractional Brownian motion

This process can be seen as an isotropic extension of the fractional Brownian motion. Indeed, for the fBm, we have for all

s, t \in R_{+}

E {[X_{t} - X_{s}]}^{2} = | t - s |^{2 H}

A natural idea to extend this process for a set of index

T \subset R_{+}^{N}

is to substitute the absolute value by a norm. We get the Levy fractional Brownian motion, which is defined to be a centered Gaussian process of covariance function

\begin{matrix} E [X_{s} X_{t}] = \frac{1}{2} [∥ s ∥^{2 H} + ∥ t ∥^{2 H} - ∥ t - s ∥^{2 H}] \end{matrix}

(1)

There are several definitions of this process by its trajectories. Among these, it can be defined as integral against white noise. Lindstrom stated the following (see [9] ).

Proposition 1 The process defined by

\begin{matrix} X_{t} = \int_{R^{N}} [∥ t - u ∥^{H - \frac{N}{2}} - ∥ u ∥^{H - \frac{N}{2}}] W (d u) \end{matrix}

(2)

is a Levy fractional Brownian motion up to a multiplicative constant.

Proof

This process is obviously Gaussian and centered. Thus we only have to show that the covariance function is of the form ( 1 ). We have

\begin{matrix} E {[X_{s} - X_{t}]}^{2} & = & \int_{R^{N}} {[∥ t - u ∥^{H - \frac{N}{2}} - ∥ s - u ∥^{H - \frac{N}{2}}]}^{2} . d u \end{matrix}

\begin{matrix} = & \int_{R^{N}} {[∥ t - s - u ∥^{H - \frac{N}{2}} - ∥ u ∥^{H - \frac{N}{2}}]}^{2} . d u \end{matrix}

We consider the change of variables from

R^{N}

into itself,

v = φ (u)

, where

φ

is the linear application which maps the canonic basis of

R^{N}

to the orthonormal basis

(e_{1} = \frac{t - s}{∥ t - s ∥}, e_{2}, . . ., e_{N})

. The differential of

φ

in any

u \in R^{N}

is itself and the Jacobian

\begin{matrix} J φ_{(v_{1}, . . ., v_{N})}^{- 1} = | det (φ^{- 1}) | = 1 \end{matrix}

because the matrix of

φ

is orthogonal.

We have

\begin{matrix} ∥ t - s - u ∥^{2} & = & (t - s - u)^{2} = ∥ t - s ∥^{2} - 2 < t - s, u > + ∥ u ∥^{2} \end{matrix}

\begin{matrix} = & ∥ t - s ∥^{2} - 2 ∥ t - s ∥ . v_{1} + ∥ v ∥^{2} \end{matrix}

\begin{matrix} = & {(∥ t - s ∥ . ε_{1} - v)}^{2} \end{matrix}

We obtain

\begin{matrix} E {[X_{s} - X_{t}]}^{2} & = & \int_{R^{N}} {[∥ ∥ t - s ∥ . ε_{1} - v ∥^{H - \frac{N}{2}} - ∥ v ∥^{H - \frac{N}{2}}]}^{2} . d v \end{matrix}

\begin{matrix} = & ∥ t - s ∥^{2 H - N} \int_{R^{N}} {[∥ ε_{1} - \frac{v}{∥ t - s ∥} ∥^{H - \frac{N}{2}} - ∥ \frac{v}{∥ t - s ∥} ∥^{H - \frac{N}{2}}]}^{2} . d v \end{matrix}

and after the second change of variables,

v = ∥ t - s ∥ . w = ∥ t - s ∥ I d . w

we get

\begin{matrix} E {[X_{s} - X_{t}]}^{2} = ∥ t - s ∥^{2 H} {\underset{︸}{\int_{R^{N}} {[∥ ε_{1} - w ∥^{H - \frac{N}{2}} - ∥ w ∥^{H - \frac{N}{2}}]}^{2} . d w}}_{K_{N, H}} \end{matrix}

therefore

\begin{matrix} E [X_{s} X_{t}] = K_{N, H} [∥ s ∥^{2 H} + ∥ t ∥^{2 H} - ∥ t - s ∥^{2 H}] \end{matrix}

□

The harmonizable representation of fractional Brownian motion can also be generalized. Before that, let's recall briefly definitions of white noise and its Fourier transform.

In the following, we will denote

L_{C}^{2} (R^{N})

the set of functions

f : R^{N} \to C

such that

\int_{R^{N}} | f (u) |^{2} d u < \infty

Definition 1 The complex isonormal process is defined to be a centered Gaussian process

W = {W (f); f \in L_{C}^{2} (R^{N})}

such that

\forall f, g \in L_{C}^{2} (R^{N}); E [W (f) \bar{W (g)}] = \int_{R^{N}} f (u) \bar{g (u)} . d u

Then, white noise

W

can be defined by

W (E) = W (1_{E})

Definition 2 A Gaussian process

{\hat{W} (f); f \in L_{C}^{2} (R^{N})}

is said to be the Fourier transform of a complex isonormal process

{W (f); f \in L_{C}^{2} (R^{N})}

if for all

f \in L_{C}^{2} (R^{N})

\hat{W} (f) = W (\hat{f})

where

\hat{f}

if the Fourier transform of the function

f

The Fourier transform of white noise is defined in the same way.

This complex measure is usually used to define the harmonizable representation of fractional Brownian motion

B_{t}^{H} = \int_{R} \frac{e^{i t ξ} - 1}{{| ξ |}^{H + \frac{1}{2}}} . \hat{W} (d ξ)

that can be generalized in the following.

Proposition 2 The process defined by

\begin{matrix} X_{t} = \int_{R^{N}} \frac{e^{i 〈 t, ξ 〉} - 1}{∥ ξ ∥^{H + \frac{N}{2}}} . \hat{W} (d ξ) \end{matrix}

(3)

where

\hat{W}

is the Fourier transform of white noise in

R^{N}

, is a Levy fractional Brownian motion up to a multiplicative constant.

Proof

As will be done for multifractional Brownian field, the Fourier transform of the kernel of representation ( 2 ) could be directly computed. But as this representation defines a real centered Gaussian process, it is enough to show that the covariance function has the form ( 1 ).

For all

t \in R^{N}

, let's denote by

f_{t}

the function

ξ \mapsto \frac{e^{i < t, ξ >} - 1}{∥ ξ ∥^{H + \frac{N}{2}}}

and consider the centered Gaussian process

X = {X_{t} = \hat{W} (f_{t}); t \in R_{+}^{N}}

First of all, let's show that, almost surely,

\hat{W} (f_{t}) \in R

In fact, using

\hat{W} (f_{t}) = W ({\hat{f}}_{t})

, showing that

{\hat{f}}_{t} \in R

is sufficient.

Indeed, by

E {[I m (W (f))]}^{2} = \int_{R^{N}} {(I m (f))}^{2} = 0

f \in R

imply

W (f) \in R

almost surely, and

\begin{matrix} I m ({\hat{f}}_{t} (x)) & = & \int_{R^{N}} I m (e^{- i < x, u >} f_{t} (u)) . d u \end{matrix}

\begin{matrix} = & \int_{R^{N}} \frac{sin < t - x, u >}{∥ u ∥^{H + \frac{N}{2}}} . d u \end{matrix}

\begin{matrix} = & 0 \end{matrix}

by parity.

The process

X

is therefore real and its covariance function is

\begin{matrix} E [X_{s} X_{t}] & = & E [\hat{W} (f_{s}) \bar{\hat{W} (f_{t})}] \end{matrix}

\begin{matrix} = & \int_{R^{N}} \frac{(e^{i < s, ξ >} - 1) (e^{- i < t, ξ >} - 1)}{∥ ξ ∥^{2 H + N}} . d ξ \end{matrix}

\begin{matrix} = & \int_{R^{N}} \frac{e^{i < s - t, ξ >} - e^{i < s, ξ >} - e^{- i < t, ξ >} + 1}{∥ ξ ∥^{2 H + N}} . d ξ \end{matrix}

Then we have to consider

3

integrals of the form

\int_{R^{N}} \frac{1 - e^{i < t, ξ >}}{∥ ξ ∥^{2 H + N}} . d ξ

As in proposition 1 , for

t \in R^{N}

fixed, consider the change of variables from

R^{N}

into itself,

u = φ (ξ)

where

φ

is the linear application which maps the canonic basis of

R^{N}

to the orthonormal basis

(e_{1} = \frac{t}{∥ t ∥}, e_{2}, . . ., e_{N})

Then, we get

\int_{R^{N}} \frac{1 - e^{i < t, ξ >}}{∥ ξ ∥^{2 H + N}} . d ξ = \int_{R^{N}} \frac{1 - e^{i ∥ t ∥ . u_{1}}}{∥ u ∥^{2 H + N}} . d u

After the second change of variables

\begin{matrix} v = ∥ t ∥ . u = ∥ t ∥ I d . u \end{matrix}

\begin{matrix} d v = ∥ t ∥^{N} . d u \end{matrix}

we get

\begin{matrix} \int_{R^{N}} \frac{1 - e^{i < t, ξ >}}{∥ ξ ∥^{2 H + N}} . d ξ = \frac{∥ t ∥^{2 H + N}}{∥ t ∥^{N}} {\underset{︸}{\int_{R^{N}} \frac{1 - e^{i v_{1}}}{∥ v ∥^{2 H + N}} . d v}}_{C_{N, H} > 0} \end{matrix}

Proceeding the same way for the

2

other integrals, we can conclude

\begin{matrix} E [X_{s} X_{t}] = C_{N, H} [∥ s ∥^{2 H} + ∥ t ∥^{2 H} - ∥ t - s ∥^{2 H}] \end{matrix}

which shows that the process

{\frac{1}{\sqrt{C_{N, H}}} \hat{W} (f_{t}), t \in R_{+}^{N}}

is a Levy fractional Brownian motion.

□

2.2 Fractional Brownian sheet

On the contrary to the Levy fractional Brownian motion, this process is not isotropic. In particular, we can have different Hurst parameters in each of the

N

directions.

For the fBm, we have for all

s, t \in R_{+}

E [X_{s} X_{t}] = \frac{1}{2} [s^{2 H} + t^{2 H} - | t - s |^{2 H}]

As in the definition of Brownian sheet, another way to generalize fBm is to set the covariance equal to the tensor product of one dimensional covariances. Then, fractional Brownian sheet (fBs) is defined to be a centered Gaussian process of covariance function

\begin{matrix} E [X_{s} X_{t}] =^{N} \prod_{i = 1} \frac{1}{2} (s_{i}^{2 H_{i}} + t_{i}^{2 H_{i}} - | t_{i} - s_{i} |^{2 H_{i}}) \end{matrix}

(4)

As in the isotropic case, this process has two different representations by its trajectories.

Proposition 3 The process defined by

X_{t} = {\int_{R^{N}}}^{N} \prod_{i = 1} [| t_{i} - u_{i} |^{H_{i} - \frac{1}{2}} - | u_{i} |^{H_{i} - \frac{1}{2}}] W (d u)

is a fractional Brownian sheet, up to a multiplicative constant.

Remark 1 In [8] , Pontier/Leger introduced another moving average representation of fractional Brownian sheet.

X_{t} = {\int_{R^{N}}}^{N} \prod_{i = 1} [(t_{i} - u_{i})_{+}^{H_{i} - \frac{1}{2}} - (- u_{i})_{+}^{H_{i} - \frac{1}{2}}] W (d u)

Proof

This process is obviously Gaussian and centered. Thus, we only need to show that its covariance function has the expected form. We compute

\begin{matrix} E [X_{s} X_{t}] & = & {\int_{R^{N}}}^{N} \prod_{i = 1} [| s_{i} - u_{i} |^{H_{i} - \frac{1}{2}} - | u_{i} |^{H_{i} - \frac{1}{2}}] [| t_{i} - u_{i} |^{H_{i} - \frac{1}{2}} - | u_{i} |^{H_{i} - \frac{1}{2}}] . d u \end{matrix}

\begin{matrix} = & ^{N} \prod_{i = 1} \int_{R} [| s_{i} - u_{i} |^{H_{i} - \frac{1}{2}} - | u_{i} |^{H_{i} - \frac{1}{2}}] [| t_{i} - u_{i} |^{H_{i} - \frac{1}{2}} - | u_{i} |^{H_{i} - \frac{1}{2}}] . d u_{i} \end{matrix}

We can see that the factor corresponding to each

i

, is the covariance of a fBm with Hurst parameter

H_{i}

(or a Levy fractional Brownian motion with

N = 1

Then we have

\begin{matrix} E [X_{s} X_{t}] =^{N} \prod_{i = 1} K_{1, H_{i}} [| s_{i} |^{2 H_{i}} + | t_{i} |^{2 H_{i}} - | t_{i} - s_{i} |^{2 H_{i}}] \end{matrix}

□

This process also has an harmonizable representation, using the Fourier transform of the white noise in

R^{N}

as in the previous paragraph.

Proposition 4 For all

t = (t_{i})

, consider the function

φ_{t}

such that for all

ξ = (ξ_{i})

φ_{t} (u) =^{N} \prod_{m = 1} \frac{e^{i t_{m} ξ_{m}} - 1}{| ξ_{m} |^{H_{m} + \frac{1}{2}}}

The process defined by

X_{t} = \hat{W} (φ_{t}) = {\int_{R^{N}}}^{N} \prod_{m = 1} \frac{e^{i t_{m} ξ_{m}} - 1}{| ξ_{m} |^{H_{m} + \frac{1}{2}}} \hat{W} (d ξ)

is a fractional Brownian sheet, up to a multiplicative constant.

Proof

As in the previous proposition, let's compute the covariance function of this process.

\begin{matrix} E [X_{s} X_{t}] & = & {\int_{R^{N}}}^{N} \prod_{m = 1} \frac{(e^{i s_{m} ξ_{m}} - 1) (e^{- i t_{m} ξ_{m}} - 1)}{| ξ_{m} |^{2 H_{m} + 1}} . d ξ \end{matrix}

\begin{matrix} = & ^{N} \prod_{m = 1} \int_{R} \frac{(e^{i s_{m} ξ_{m}} - 1) (e^{- i t_{m} ξ_{m}} - 1)}{| ξ_{m} |^{2 H_{m} + 1}} . d ξ_{m} \end{matrix}

\begin{matrix} = & ^{N} \prod_{m = 1} C_{1, H_{m}} [| s_{m} |^{2 H_{m}} + | t_{m} |^{2 H_{m}} - | t_{m} - s_{m} |^{2 H_{m}}] \end{matrix}

using the same argument of the previous proposition.

□

Remark 2 The processes defined in propositions 3 and 4 are proved to have the same law. In fact, as a particular case of proposition 10 , they are indistinguishable.

2.3 Stationarity of increments and self similarity

Let us start by recalling the notion of increments in

R_{+}^{N}

For a function

f : [0, 1]^{N} \to R

and

h \in R

, one usually define the progressive difference in direction

ε_{i}

\begin{matrix} Δ_{h, i} f (x) = {\begin{matrix} f (x + h ε_{i}) - f (x) & if x, x + h ε_{i} \in [0, 1]^{N} \\ 0 & either \end{matrix} \end{matrix}

and for

h \in R^{N}

and

A = (i_{1}, . . ., i_{k})

Δ_{h, A} f = Δ_{h_{i_{1}}, i_{1}} f \circ \cdot \cdot \cdot \circ Δ_{h_{i_{k}}, i_{k}} f

Despite the temptation to define the increments by

X_{t} - X_{s}

as in one dimension, it is better to set

\begin{matrix} Δ X_{s, t} & = & Δ_{t - s, (1, . . ., N)} X_{s} \end{matrix}

\begin{matrix} = & \sum_{r \in {0, 1}^{N}} (- 1)^{N - \sum_{l} r_{l}} X_{{[s_{i} + r_{i} (t_{i} - s_{i})]}_{i}} \end{matrix}

(5)

If there exists

i \in {1, . . ., N}

such that

s_{i} = t_{i}

, we have

Δ X_{s, t} = 0

. Then, we consider

I = {i = 1, . . ., N; s_{i} \neq t_{i}}

and

Δ_{t - s, I} X_{s} = \sum_{r \in {0, 1}^{# I}} (- 1)^{# I - \sum_{l} r_{l}} X_{{[s_{i} + r_{i} (t_{i} - s_{i})]}_{i \in I}}

2.3.1 Isotropic case

In the isotropic case, the following extension of fBm's properties are well known (see [9] ).

Proposition 5 Let

X = {X_{t}; t \in R_{+}^{N}}

be a Levy fractional Brownian motion.

We have the two following properties for all

h \in R_{+}^{N}

and

a > 0

\begin{matrix} X_{t + h} - X_{h} & \overset{(d)}{=} & X_{t} - X_{0} \end{matrix}

\begin{matrix} X_{a t} & \overset{(d)}{=} & a^{H} X_{t} \end{matrix}

where

\overset{(d)}{=}

means equality of finite dimensional distributions.

Proof

For all

s

and

t

R_{+}^{N}

, we have

\begin{matrix} E [(X_{s + h} - X_{h}) (X_{t + h} - X_{h})] & = & \frac{1}{2} (E {[X_{s + h} - X_{h}]}^{2} + E {[X_{t + h} - X_{h}]}^{2} - E {[X_{t + h} - X_{s + h}]}^{2}) \end{matrix}

\begin{matrix} = & \frac{1}{2} (∥ s ∥^{2 H} + ∥ t ∥^{2 H} - ∥ t - s ∥^{2 H}) \end{matrix}

\begin{matrix} = & E [X_{s} X_{t}] \end{matrix}

For self-similarity, we compute

\begin{matrix} E [X_{a s} X_{a t}] & = & \frac{1}{2} (E {[X_{a s}]}^{2} + E {[X_{a t}]}^{2} - E {[X_{a t} - X_{a s}]}^{2}) \end{matrix}

\begin{matrix} = & \frac{1}{2} (∥ a s ∥^{2 H} + ∥ a t ∥^{2 H} - ∥ a t - a s ∥^{2 H}) \end{matrix}

\begin{matrix} = & a^{2 H} E [X_{s} X_{t}] \end{matrix}

\begin{matrix} = & E [a^{H} X_{s} a^{H} X_{t}] \end{matrix}

□

Proposition 5 implies the stationarity of increments ( 5 ).

Proposition 6 The increments of Levy fractional Brownian are stationary, ie for all

h \in R_{+}^{N}

Δ X_{h, t + h} \overset{(d)}{=} Δ X_{0, t}

Proof

We fix

h \in R_{+}^{N}

and write

\begin{matrix} Δ X_{h, t + h} = \sum_{r \in {0, 1}^{N} - {0}} (- 1)^{N - \sum_{l} r_{l}} (X_{{[h_{i} + r_{i} t_{i}]}_{i}} - X_{h}) \end{matrix}

then in the development of

E [Δ X_{h, s + h} Δ X_{h, t + h}]

, we only have terms of the form

\begin{matrix} E [(X_{{[h_{i} + r_{i} s_{i}]}_{i}} - X_{h}) (X_{{[h_{i} + ρ_{i} t_{i}]}_{i}} - X_{h})] = E [X_{{[r_{i} s_{i}]}_{i}} X_{{[ρ_{i} t_{i}]}_{i}}] \end{matrix}

using the previous proposition. Therefore we have

E [Δ X_{h, s + h} Δ X_{h, t + h}] = E [Δ X_{0, s} Δ X_{0, t}]

□

2.3.2 Non-isotropic case

In the non-isotropic case, the properties of self-similarity and stationarity of increments have been stated by Léger/Pontier (cf [8] ). Here, we give another proof based on the covariance function rather than the moving average representation.

Proposition 7 Let

X = {X_{t}; t \in R_{+}^{N}}

be a fractional Brownian sheet. We have the two following properties for all

h \in R_{+}^{N}

and

a > 0

Δ X_{h, t + h} \overset{(d)}{=} Δ X_{0, t}

X_{a t} \overset{(d)}{=} a^{\sum_{i} H_{i}} X_{t}

Proof

We consider

N

independent fBm

X^{(1)}, . . . X^{(N)}

of Hurst parameter

H_{i}

, and the process

Y = {Y_{t}; t \in R_{+}^{N}}

such that

Y_{t} = \prod_{i = 1}^{N} X_{t_{i}}^{(i)}

. We can see easily that

X

and

Y

have the same covariance function. The same result follows for the increments

{Δ X_{h, t + h}; t \in R_{+}^{N}}

and

{Δ Y_{h, t + h}; t \in R_{+}^{N}}

. As a consequence, from

\begin{matrix} Δ Y_{h, t + h} & = & \sum_{r \in {0, 1}^{N}} (- 1 {)^{N - \sum_{l} r_{l}}}^{N} \prod_{i = 1} X_{h_{i} + r_{i} t_{i}}^{(i)} \end{matrix}

\begin{matrix} = & ^{N} \prod_{i = 1} [X_{t_{i} + h_{i}}^{(i)} - X_{h_{i}}^{(i)}] \end{matrix}

then we have

\begin{matrix} E [Δ X_{h, s + h} Δ X_{h, t + h}] & = & E [Δ Y_{h, s + h} Δ Y_{h, t + h}] \end{matrix}

\begin{matrix} = & ^{N} \prod_{i = 1} {\underset{︸}{E [(X_{s_{i} + h_{i}}^{(i)} - X_{h_{i}}^{(i)}) (X_{t_{i} + h_{i}}^{(i)} - X_{h_{i}}^{(i)})]}}_{E [X_{s_{i}}^{(i)} X_{t_{i}}^{(i)}]} \end{matrix}

\begin{matrix} = & E [Δ X_{0, s} Δ X_{0, t}] \end{matrix}

For self-similarity, we compute for all

a > 0

\begin{matrix} E [X_{a s} X_{a t}] & = & ^{N} \prod_{i = 1} \frac{1}{2} ((a s_{i})^{2 H_{i}} + (a t_{i})^{2 H_{i}} - | a t_{i} - a s_{i} |^{2 H_{i}}) \end{matrix}

\begin{matrix} = & {a^{2 \sum_{i} H_{i}}}^{N} \prod_{i = 1} \frac{1}{2} (s_{i}^{2 H_{i}} + t_{i}^{2 H_{i}} - | t_{i} - s_{i} |^{2 H_{i}}) \end{matrix}

\begin{matrix} = & a^{2 \sum_{i} H_{i}} E [X_{s} X_{t}] \end{matrix}

\begin{matrix} = & E [a^{\sum_{i} H_{i}} X_{s} a^{\sum_{i} H_{i}} X_{t}] \end{matrix}

□

Therefore, we can conclude that both extensions of fBm satisfy the properties of self-similarity and stationarity of increments.

3 The multifractional Brownian motion's case

Once again, we can consider two different kinds of multi-parameter extension of mBm : isotropic and anisotropic extension. Note, first of all, that mBm already has a multi-parameter extension. Indeed, the formulation of Benassi/Jaffard/Roux in [4] was done for

t \in R^{N}

. We will see that it can be considered as an isotropic extension.

3.1 Isotropic extension

To define an isotropic extension of the mBm, the natural way is to substitute the constant

H

of the moving average representation of the Levy fractional Brownian motion, with a function.

Definition 3 Let

H : R^{N} \to (0, 1)

be a measurable function. The process

{X_{t}; t \in R_{+}^{N}}

such that

\begin{matrix} X_{t} = \int_{R^{N}} [∥ t - u ∥^{H (t) - \frac{N}{2}} - ∥ u ∥^{H (t) - \frac{N}{2}}] W (d u) \end{matrix}

(6)

is called multifractional Brownian field.

We will show that this process is the same as the process defined by Benassi/Jaffard/Roux. This result generalizes on the equivalence stated in the case

N = 1

in [6] .

Proposition 8 Let

H : R^{N} \to (0, 1)

be a measurable function. The process defined by

\begin{matrix} X_{t} = \int_{R^{N}} \frac{e^{i 〈 t, ξ 〉} - 1}{∥ ξ ∥^{H (t) + \frac{N}{2}}} . \hat{W} (d ξ) \end{matrix}

(7)

is indistinguishable, up to a multiplicative deterministic function, from the process defined by ( 6 ). This formulation is the harmonizable representation of the multifractional Brownian field.

Proof

First of all, let us compute the Fourier transform of the function

∥ . ∥^{α}

\begin{matrix} 〈 T ∥ . ∥^{α}, φ 〉 & = & 〈 ∥ . ∥^{α}, \hat{φ} 〉 \end{matrix}

\begin{matrix} = & \int_{R^{N}} ∥ t ∥^{α} (\int_{R^{N}} e^{- i < w, t >} φ (w) . d w) . d t \end{matrix}

we consider the change of variables

\begin{matrix} R^{N} \times R^{N} & \to & R^{N} \times R^{N} \end{matrix}

\begin{matrix} (w, t) & \mapsto & (w, λ = φ (t)) \end{matrix}

where

φ

is the linear application which maps the canonic basis of

R^{N}

to the orthonormal basis

(e_{1} = \frac{w}{∥ w ∥}, e_{2}, . . ., e_{N})

. We get

\begin{matrix} 〈 T ∥ . ∥^{α}, φ 〉 & = & \int_{R^{N}} \int_{R^{N}} ∥ λ ∥^{α} e^{i λ_{1} ∥ w ∥} φ (w) . d w . d λ \end{matrix}

\begin{matrix} = & \int_{R^{N}} \int_{R^{N}} \frac{∥ u ∥^{α}}{∥ w ∥^{α}} e^{- i u_{1}} φ (w) \frac{d w . d u}{∥ w ∥^{N}} \end{matrix}

using the change of variables

(w, λ) \mapsto (w, u = ∥ w ∥ λ)

. Then we have

\begin{matrix} 〈 T ∥ . ∥^{α}, φ 〉 = {\underset{︸}{(\int_{R^{N}} ∥ u ∥^{α} e^{- i u_{1}} . d u)}}_{λ_{α}} \int_{R^{N}} \frac{1}{∥ w ∥^{α + N}} φ (w) . d w \end{matrix}

Thus,

T ∥ . ∥^{α} (w) = \frac{λ_{α}}{∥ w ∥^{α + N}}

We use this result to calculate the Fourier transform of

∥ t - . ∥^{α} - ∥ . ∥^{α}

. We will use the following property : if

g (u) = f (u - α)

then

\hat{g} = e^{- i < α, v >} \hat{f} (v)

T [∥ t - . ∥^{α} - ∥ . ∥^{α}] (v) = [e^{- i < t, v >} - 1] \frac{λ_{α}}{∥ v ∥^{α + N}}

We deduce from this

T [∥ t - . ∥^{H (t) - \frac{N}{2}} - ∥ . ∥^{H (t) - \frac{N}{2}}] (v) = λ_{H (t)} \bar{(\frac{e^{i < t, v >} - 1}{∥ v ∥^{H (t) + \frac{N}{2}}})}

and

\forall t \in R^{N}

, we have almost surely

\int_{R^{N}} [∥ t - u ∥^{H (t) - \frac{N}{2}} - ∥ u ∥^{H (t) - \frac{N}{2}}] W (d u) = λ_{H (t)} \int_{R^{N}} \frac{e^{i 〈 t, ξ 〉} - 1}{∥ ξ ∥^{H (t) + \frac{N}{2}}} . \hat{W} (d ξ)

using the fact we saw previously that the second integral is almost surely real.

Therefore, by an argument of continuity, the result follows.

□

This process is obviously a centered Gaussian process. It is thus of interest to study its covariance function. The following proposition is an extension of the case

N = 1

stated in [3] .

Proposition 9 Let

{X_{t}; t \in R_{+}^{N}}

be a multifractional Brownian field. There exists a deterministic function

D_{N}^{f} : R \to R

such that the covariance function of

X

can be written

\begin{matrix} E [X_{s} X_{t}] = D_{N}^{f} (H (s) + H (t)) [∥ s ∥^{H (s) + H (t)} + ∥ t ∥^{H (s) + H (t)} - ∥ t - s ∥^{H (s) + H (t)}] \end{matrix}

(8)

Proof

The easiest way to show this result is to use the harmonizable representation.

By definition of

\hat{W}

, we have

E [X_{s} X_{t}] = \int_{R^{N}} \frac{(e^{i < s, ξ >} - 1) (e^{- i < t, ξ >} - 1)}{∥ ξ ∥^{H (s) + H (t) + N}} . d ξ

This integral has already been calculated for a Levy fractional Brownian motion with a parameter

H = \frac{H (s) + H (t)}{2}

. Then we have

E [X_{s} X_{t}] = {\underset{︸}{(\int_{R^{N}} \frac{1 - e^{i u_{1}}}{∥ u ∥^{H (s) + H (t) + N}} . d u)}}_{D_{N}^{f} (H (s) + H (t))} [∥ s ∥^{H (s) + H (t)} + ∥ t ∥^{H (s) + H (t)} - ∥ t - s ∥^{H (s) + H (t)}]

with

D_{N}^{f} (x) = \int_{R^{N}} \frac{1 - e^{i u_{1}}}{∥ u ∥^{x + N}} . d u

□

3.2 Non isotropic extension

Another way to extend the multifractional Brownian motion for a set of index included in

R_{+}^{N}

, is to copy the definition of the Brownian sheet.

Definition 4 Let

H : R_{+}^{N} \to (0, 1)^{N}

be a measurable function. The process

{X_{t}; t \in R_{+}^{N}}

such that

X_{t} = {\int_{R^{N}}}^{N} \prod_{i = 1} [| t_{i} - u_{i} |^{H_{i} (t) - \frac{1}{2}} - | u_{i} |^{H_{i} (t) - \frac{1}{2}}] W (d u)

where

W

is the white noise, is called multifractional Brownian sheet (mBs).

As in the case of the isotropic extension, there also exists a harmonizable representation of the mBs.

Proposition 10 Let

H : R_{+}^{N} \to (0, 1)^{N}

be a measurable function. For all

t = {(t_{i})}_{i \in {1; . . .; N}}

, we consider the function

φ_{t}

such that for all

ξ = (ξ_{i})

φ_{t} (u) =^{N} \prod_{m = 1} \frac{e^{i t_{m} ξ_{m}} - 1}{| ξ_{m} |^{H_{m} (t) + \frac{1}{2}}}

The process defined by

X_{t} = \hat{W} (φ_{t}) = {\int_{R^{N}}}^{N} \prod_{m = 1} \frac{e^{i t_{m} ξ_{m}} - 1}{| ξ_{m} |^{H_{m} (t) + \frac{1}{2}}} \hat{W} (d ξ)

is indistinguishable, up to a multiplicative deterministic function, from the process defined previously. This formulation is the harmonizable representation of the multifractional Brownian sheet.

Proof

We have already seen that for each

m \in {1, . . ., N}

T [| t_{m} - . |^{H_{m} (t) - \frac{1}{2}} - | . |^{H_{m} (t) - \frac{1}{2}}] (ξ_{m}) = λ_{H_{m} (t)} \bar{(\frac{e^{i t_{m} ξ_{m}} - 1}{| ξ_{m} |^{H_{m} (t) + \frac{1}{2}}})}

Moreover, we compute

\begin{matrix} T (^{N} \prod_{m = 1} [| t_{m} - . |^{H_{m} (t) - \frac{1}{2}} - | . |^{H_{m} (t) - \frac{1}{2}}]) & (ξ) \end{matrix}

\begin{matrix} = & \int_{R^{N}} {e^{- i < ξ, x >}}^{N} \prod_{m = 1} [| t_{m} - x_{m} |^{H_{m} (t) - \frac{1}{2}} - | x_{m} |^{H_{m} (t) - \frac{1}{2}}] . d x \end{matrix}

\begin{matrix} = & {\int_{R^{N}}}^{N} \prod_{m = 1} e^{- i ξ_{m} x_{m}} [| t_{m} - x_{m} |^{H_{m} (t) - \frac{1}{2}} - | x_{m} |^{H_{m} (t) - \frac{1}{2}}] . d x \end{matrix}

\begin{matrix} = & ^{N} \prod_{m = 1} T [| t_{m} - . |^{H_{m} (t) - \frac{1}{2}} - | . |^{H_{m} (t) - \frac{1}{2}}] (ξ_{m}) \end{matrix}

Therefore

{\underset{︸}{(^{N} \prod_{i = 1} λ_{m} (t))}}_{λ (t)} \hat{W} \bar{(^{N} \prod_{m = 1} \frac{e^{i t_{m} .} - 1}{| . |^{H_{m} (t) + \frac{1}{2}}})} = W (^{N} \prod_{m = 1} [| t_{m} - . |^{H_{m} (t) - \frac{1}{2}} - | . |^{H_{m} (t) - \frac{1}{2}}])

We use the same arguments as in proposition 8 to conclude.

□

The following proposition shows that the covariance structure of multifractional Brownian sheet, is a generalization of the fBs's one.

Proposition 11 Let

{X_{t}; t \in R_{+}^{N}}

be a multifractional Brownian sheet. There exists a deterministic function

D^{s} : R^{N} \to R

such that

\begin{matrix} E [X_{s} X_{t}] = D^{s} {(H (s) + H (t))}^{N} \prod_{m = 1} [| s_{m} |^{H_{m} (s) + H_{m} (t)} + | t_{m} |^{H_{m} (s) + H_{m} (t)} - | t_{m} - s_{m} |^{H_{m} (s) + H_{m} (t)}] \end{matrix}

(9)

Proof

As usually, we use the harmonizable representation of the process

\begin{matrix} E [X_{s} X_{t}] & = & {\int_{R^{N}}}^{N} \prod_{m = 1} \frac{(e^{i s_{m} ξ_{m}} - 1) (e^{- i t_{m} ξ_{m}} - 1)}{| ξ_{m} |^{H_{m} (s) + H_{m} (t) + 1}} . d ξ \end{matrix}

\begin{matrix} = & ^{N} \prod_{m = 1} \int_{R} \frac{(e^{i s_{m} ξ_{m}} - 1) (e^{- i t_{m} ξ_{m}} - 1)}{| ξ_{m} |^{H_{m} (s) + H_{m} (t) + 1}} . d ξ_{m} \end{matrix}

We remark that the factor corresponding to each

m

, is the covariance of a multifractional Brownian motion, with has already been calculated. Therefore we have

\begin{matrix} E [X_{s} X_{t}] =^{N} \prod_{m = 1} D_{1}^{f} (H_{m} (s) + H_{m} (t)) [| s_{m} |^{H_{m} (s) + H_{m} (t)} + | t_{m} |^{H_{m} (s) + H_{m} (t)} - | t_{m} - s_{m} |^{H_{m} (s) + H_{m} (t)}] \end{matrix}

□

Remark 3 The form of the previous covariance function gives the idea to consider the process

Y = {Y_{t}; t \in R_{+}^{N}}

defined from

N

independent multifractional Brownian motions

X^{(i)}

with parameter

H_{i}

Y_{t} = X_{t^{(1)}}^{(1)} . . . X_{t^{(N)}}^{(N)}

Although

Y

is not a Gaussian process, it is easily seen that it has the same covariance function as a multifractional Brownian sheet. This remark will be often used in the following.

4 Regularity

A lot of properties are known about the regularity of the trajectories of Brownian motion and fractional Brownian motion. As we will see, in the case of the multi-parameter extension of the mBm, we have to make some assumptions about the regularity of

H

before studying the continuity of trajectories. In the definitions of mBm (cf [1] and [4] ), the function

H

is supposed to be Hölder continuous.

4.1 Continuity of the two extensions

We first recall the Kolmogorov's criterion.

Theorem 1 (Kolmogorov) Consider a process

X = {X_{t}; t \in R_{+}^{N}}

such that there exists

C > 0

p > 0

and

γ > N

\begin{matrix} \forall s, t \in R_{+}^{N}, E [| X_{t} - X_{s} |^{p}] \leq C ∥ t - s ∥^{γ} \end{matrix}

(10)

Then, there exists a modification

Y = {Y_{t}; t \in R_{+}^{N}}

X

that is Hölder continuous of any order

q \in (0, \frac{γ - N}{p})

As usually, the quantity

E [| X_{t} - X_{s} |^{2}]

is studied for

s, t \in [a, b]

where

a ≼ b

and then, a patching argument is used to extend to

s, t \in R_{+}^{N}

4.1.1 Isotropic case

Lemma 1 For all

η

and

μ

such that

0 < η < μ < 1

, the multiplicative factor

D_{N}^{f}

of covariance function in ( 9 ), is positive and belongs to

C^{\infty} ([η, μ])

Moreover, its order

n

derivative is given by

\begin{matrix} {D_{N}^{f}}^{(n)} (x) = \int_{R^{N}} \frac{1 - e^{i u_{1}}}{∥ u ∥^{x + N}} {ln}^{n} \frac{1}{∥ u ∥} . d u \end{matrix}

(11)

Proof

As the integral of a positive function,

D_{N}^{f}

is positive. By an argument of uniform convergence of integrals ( 11 ) on

[η, μ]

D_{N}^{f}

C^{\infty} ([η, μ])

and the derivatives are obtained by derivations of the integrand.

□

Proposition 12 For all

s, t \in [a, b]

, we have

\begin{matrix} \frac{1}{2} E {[X_{t} - X_{s}]}^{2} & = & D [H (s) + H (t)] \times ∥ t - s ∥^{H (s) + H (t)} \end{matrix}

\begin{matrix} + \frac{1}{2} [\frac{\partial^{2} φ}{\partial x^{2}} (H (s) + H (t); ∥ s ∥) + \frac{\partial^{2} φ}{\partial x^{2}} (H (s) + H (t); ∥ t ∥)] \times {(H (t) - H (s))}^{2} \end{matrix}

\begin{matrix} + O_{a, b} [(H (t) - H (s)) (∥ t ∥ - ∥ s ∥)] + o_{a, b} {(H (t) - H (s))}^{2} \end{matrix}

(12)

where

φ (x, y) = D (x) y^{x}

Proof

Using the covariance function of the multifractional Brownian field, we have

\begin{matrix} \frac{1}{2} E [| X_{s} - X_{t} |^{2}] & = & D [2 H (s)] ∥ s ∥^{2 H (s)} - D [H (s) + H (t)] ∥ s ∥^{H (s) + H (t)} \end{matrix}

\begin{matrix} + D [2 H (t)] ∥ t ∥^{2 H (t)} - D [H (s) + H (t)] ∥ t ∥^{H (s) + H (t)} \end{matrix}

\begin{matrix} + D [H (s) + H (t)] ∥ t - s ∥^{H (s) + H (t)} \end{matrix}

(13)

We have to get a second order expansion of this expression.

We introduce the function

φ

defined by

φ (x, y) = D (x) y^{x}

We can write

\begin{matrix} \frac{1}{2} E [| X_{s} - X_{t} |^{2}] & = & φ (2 H (s), ∥ s ∥) - φ (H (s) + H (t), ∥ s ∥) \end{matrix}

\begin{matrix} + φ (2 H (t), ∥ t ∥) - φ (H (s) + H (t), ∥ t ∥) \end{matrix}

\begin{matrix} + D [H (s) + H (t)] ∥ t - s ∥^{H (s) + H (t)} \end{matrix}

(14)

We use the second order expansion

\begin{matrix} φ (2 H (s), ∥ s ∥) - φ (H (s) + H (t), ∥ s ∥) & = & (H (s) - H (t)) \times \frac{\partial φ}{\partial x} (H (s) + H (t), ∥ s ∥) \end{matrix}

\begin{matrix} + \frac{{(H (s) - H (t))}^{2}}{2} \times \frac{\partial^{2} φ}{\partial x^{2}} (H (s) + H (t), ∥ s ∥) \end{matrix}

\begin{matrix} + o_{a, b} {(H (s) - H (t))}^{2} \end{matrix}

An inversion of roles between

s

and

t

provides the expansion of

φ (2 H (t), ∥ t ∥) - φ (H (s) + H (t), ∥ t ∥)

Then ( 14 ) becomes

\begin{matrix} \frac{1}{2} E [| X_{s} - X_{t} |^{2}] & = & (H (t) - H (s)) \times [\frac{\partial φ}{\partial x} (H (s) + H (t), ∥ t ∥) - \frac{\partial φ}{\partial x} (H (s) + H (t), ∥ s ∥)] \end{matrix}

\begin{matrix} + \frac{{(H (t) - H (s))}^{2}}{2} \times [\frac{\partial^{2} φ}{\partial x^{2}} (H (s) + H (t), ∥ s ∥) + \frac{\partial^{2} φ}{\partial x^{2}} (H (s) + H (t), ∥ t ∥)] \end{matrix}

\begin{matrix} + D [H (s) + H (t)] ∥ t - s ∥^{H (s) + H (t)} + o_{a, b} {(H (t) - H (s))}^{2} \end{matrix}

Since

(H (t) - H (s)) \times [\frac{\partial φ}{\partial x} (H (s) + H (t), ∥ t ∥) - \frac{\partial φ}{\partial x} (H (s) + H (t), ∥ s ∥)]

O_{a, b} [(H (t) - H (s)) (∥ t ∥ - ∥ s ∥)]

, the result follows.

□

Corollary 1 For all

s, t \in [a, b]

, we have

\begin{matrix} \frac{1}{2} E {[X_{t} - X_{s}]}^{2} & = & D [2 H (t)] \times ∥ t - s ∥^{2 H (t)} \end{matrix}

\begin{matrix} + \frac{\partial^{2} φ}{\partial x^{2}} (2 H (t); ∥ t ∥) \times {(H (t) - H (s))}^{2} \end{matrix}

\begin{matrix} + o_{a, b} {(H (t) - H (s))}^{2} + o_{a, b} (∥ t - s ∥^{2 H (t)}) \end{matrix}

(15)

where

φ (x, y) = D (x) y^{x}

Proof

Using the expansion of

D [H (s) + H (t)]

and

\begin{matrix} ∥ t - s ∥^{H (s) + H (t)} = ∥ t - s ∥^{2 H (t)} - (H (t) - H (s)) ∥ t - s ∥^{2 H (t)} ln ∥ t - s ∥ + o_{a, b} {(H (t) - H (s))}^{2} \end{matrix}

we get

\begin{matrix} D [H (s) + H (t)] \times ∥ t - s ∥^{H (s) + H (t)} & = & D [2 H (t)] \times ∥ t - s ∥^{2 H (t)} \end{matrix}

(16)

\begin{matrix} + o_{a, b} (∥ t - s ∥^{2 H (t)}) + o_{a, b} {(H (t) - H (s))}^{2} \end{matrix}

Moreover as

H (t) < 1

for all

t \in [a, b]

, we have

ε = 1 - H (t) > 0

and

\begin{matrix} 2 (H (t) - H (s)) (∥ t ∥ - ∥ s ∥) & = & 2 (H (t) - H (s)) {(∥ t ∥ - ∥ s ∥)}^{\frac{ε}{2}} \times {(∥ t ∥ - ∥ s ∥)}^{1 - \frac{ε}{2}} \end{matrix}

\begin{matrix} \leq & {(H (t) - H (s))}^{2} {(∥ t ∥ - ∥ s ∥)}^{ε} + {(∥ t ∥ - ∥ s ∥)}^{2 - ε} \end{matrix}

that implies

\begin{matrix} (H (t) - H (s)) (∥ t ∥ - ∥ s ∥) = o_{a, b} {(H (t) - H (s))}^{2} + o_{a, b} (∥ t - s ∥^{2 H (t)}) \end{matrix}

(17)

We conclude by ( 12 ), ( 16 ) and ( 17 ) using first order expansion of

\frac{\partial^{2} φ}{\partial x^{2}}

x

and

y

□

Using the continuity of

D

D^{'}

and

D^{''}

, we can state from the previous proposition

Corollary 2 There exist positive constants

K

and

L

such that

\begin{matrix} \forall s, t \in [a, b]; E {[X_{t} - X_{s}]}^{2} \leq K ∥ t - s ∥^{2 H (t)} + L {| H (t) - H (s) |}^{2} \end{matrix}

(18)

Corollary 3 Suppose

H

β

-Hölder continuous. There exists a constant

M

such that

\begin{matrix} \forall s, t \in [a, b]; E {[X_{t} - X_{s}]}^{2} \leq M ∥ t - s ∥^{2 (β \land H (t))} \end{matrix}

(19)

4.1.2 Non-isotropic case

Lemma 2 There exists positive constants

K

and

L

such that

\begin{matrix} \forall s, t \in [a, b]; E [| X_{t} - X_{s} |^{2}] \leq K ∥ t - s ∥^{2 {min}_{i} H_{i} (t)} + L ∥ H (t) - H (s) ∥^{2} \end{matrix}

(20)

Proof

By remark 3 , we have

\begin{matrix} E {[X_{s} - X_{t}]}^{2} & = & E {[^{N} \prod_{i = 1} X_{s^{(i)}}^{(i)} -^{N} \prod_{i = 1} X_{t^{(i)}}^{(i)}]}^{2} \end{matrix}

\begin{matrix} = & E [(^{N} \prod_{i = 1} X_{s^{(i)}}^{(i)} - X_{t^{(1)}}^{(1)} \prod_{i > 1} X_{s^{(i)}}^{(i)}) + (X_{t^{(1)}}^{(1)} \prod_{i > 1} X_{s^{(i)}}^{(i)} - X_{t^{(1)}}^{(1)} X_{t^{(2)}}^{(2)} \prod_{i > 2} X_{s^{(i)}}^{(i)}) \end{matrix}

\begin{matrix} + \cdot \cdot \cdot + {((^{N - 1} \prod_{i = 1} X_{t^{(i)}}^{(i)}) X_{s^{(N)}}^{(N)} -^{N} \prod_{i = 1} X_{t^{(i)}}^{(i)})]}^{2} \end{matrix}

By the inequality of convexity

{(\sum a_{i})}^{2} \leq n \sum a_{i}^{2}

, we get

\begin{matrix} E {[X_{s} - X_{t}]}^{2} \leq N {E {[\prod_{i > 1} X_{s^{(i)}}^{(i)}]}^{2} E {[X_{s^{(1)}}^{(1)} - X_{t^{(1)}}^{(1)}]}^{2} + \dots + E {[X_{s^{(N)}}^{(N)} - X_{t^{(N)}}^{(N)}]}^{2} E {[^{N - 1} \prod_{i = 1} X_{s^{(i)}}^{(i)}]}^{2}} \end{matrix}

Since there exists a constant

M = M_{a, b}

such that

\forall t \in [a, b], \forall i; E {[X_{t^{(i)}}^{(i)}]}^{2} \leq M

we get

\begin{matrix} E {[X_{s} - X_{t}]}^{2} \leq N M^{n - 1} \sum_{i = 1}^{N} E {[X_{s^{(i)}}^{(i)} - X_{t^{(i)}}^{(i)}]}^{2} \end{matrix}

(21)

Using

E {[X_{s^{(i)}}^{(i)} - X_{t^{(i)}}^{(i)}]}^{2} \leq K_{i} | s^{(i)} - t^{(i)} |^{2 H_{i} (t)} + L_{i} {(H_{i} (s) - H_{i} (t))}^{2}; \forall i = 1, . . ., N

( 21 ) implies

\begin{matrix} E {[X_{t} - X_{s}]}^{2} \leq N M^{n - 1} [(\sum_{i = 1}^{N} K_{i}) ∥ t - s ∥^{2 {min}_{i} H_{i} (t)} + (\sum_{i = 1}^{N} L_{i}) ∥ H (t) - H (s) ∥^{2}] \end{matrix}

□

Corollary 4 Suppose

H

β

-Hölder continuous. There exists a positive constant

M

such that

\begin{matrix} \forall s, t \in [a, b]; E {[X_{t} - X_{s}]}^{2} \leq M ∥ t - s ∥^{2 (β \land {min}_{i} H_{i} (t))} \end{matrix}

(22)

4.1.3 Existence of a continuous modification

In both isotropic and anisotropic cases, under Hölder regularity assumptions for

H

, we have an inequality

E {[X_{t} - X_{s}]}^{2} \leq K ∥ t - s ∥^{α}

But to use the Kolmogorov criterion, we need to have

α > N

As the random variable

X_{t} - X_{s}

is Gaussian, we can write, for each integer

n

E {[X_{t} - X_{s}]}^{2 n} \leq λ_{n} K ∥ t - s ∥^{n . α}

and choose

n

such that

n . α > N

We conclude by a classical patching argument. For

a

and

b

, Kolmogorov's theorem gives a continuous process

Y^{a, b} = {Y_{t}^{a, b}; t \in [a, b]}

. Consider

a^{'}

and

b^{'}

such that

[a, b] \subset [a^{'}, b^{'}]

. The processes

Y^{a, b}

and

Y^{a^{'}, b^{'}}

coincide on

[a, b]

Thus

\forall t \in [a, b]; P {Y_{t}^{a, b} = Y_{t}^{a^{'}, b^{'}}} = 1

and, by continuity

P {Y_{t}^{a, b} = Y_{t}^{a^{'}, b^{'}}; \forall t \in [a, b]} = 1

Then we can define a process

Y

R_{+}^{N}

who coincides with

Y^{a, b}

[a, b]

and we can seen easily that this process is continuous.

4.2 Hölder exponents

The notion of Hölder function is well known. It is interesting to consider a localized version of this notion.

For the paths of a process

X

, one usually define two kinds of exponent (see [1] , [2] ):

∙ the pointwise Hölder exponent

\begin{matrix} α (t_{0}) & = & sup {α; {lim}_{h \to 0} \frac{| X_{t_{0} + h} - X_{t_{0}} |}{∥ h ∥^{α}} = 0} \end{matrix}

\begin{matrix} = & sup {α; {limsup}_{ρ \to 0} \frac{{sup}_{s, t \in B (t_{0}, ρ)} | X_{t} - X_{s} |}{ρ^{α}} < \infty} \end{matrix}

∙ the local Hölder exponent

\begin{matrix} \tilde{α} (t_{0}) = sup {α; {limsup}_{ρ \to 0} {sup}_{s, t \in B (t_{0}, ρ)} \frac{| X_{t} - X_{s} |}{∥ t - s ∥^{α}} < \infty} \end{matrix}

We can see easily that for all

t_{0}

, we have

\begin{matrix} \tilde{α} (t_{0}) \leq α (t_{0}) \end{matrix}

(23)

A study of these exponents, in the case of

1

D mBm, is made in [2] .

Remark 4 If

H

β

-Hölder continuous, then the local Hölder exponent

\tilde{β} (t)

H

at every point is not smaller than

β

Conversely, suppose that the local Hölder exponent of

H

at every point of a compact

[a, b]

is positive. Then

H

β

-Hölder continuous on

[a, b]

with

β = {inf}_{t \in [a, b]} \tilde{β} (t)

In the following, we suppose that

H

admits positive local Hölder exponent

\tilde{β} (t_{0})

at every point

t

Proposition 13 Let

X = {X_{t}; t \in R_{+}^{N}}

be a multifractional Brownian field.

For all

t_{0} \in R_{+}^{N}

, the local Hölder exponent of

X

t_{0}

is almost surely given by

\begin{matrix} \tilde{α} (t_{0}) = \tilde{β} (t_{0}) \land H (t_{0}) \end{matrix}

(24)

and the pointwise Hölder exponent of

X

t_{0}

satisfies almost surely

\begin{matrix} α (t_{0}) = β (t_{0}) \land H (t_{0}) \end{matrix}

(25)

where

β (t_{0})

and

\tilde{β} (t_{0})

denote the pointwise and local Hölder exponents of

H

t_{0}

As a consequence of this result, if

H

satisfies

\forall t \in R_{+}^{N}; β (t) < H (t)

the Hölder regularity of multifractional Brownian field of parameter function

H

is given by the regularity of

H

(and not by the value of

H

). This point is developed in [7] .

Proposition 14 Let

X = {X_{t}; t \in R_{+}^{N}}

be a multifractional Brownian sheet. For all

t_{0} \in R_{+}^{N}

, the local Hölder exponent of

X

t_{0}

is almost surely given by

\begin{matrix} \tilde{α} (t_{0}) = \tilde{β} (t_{0}) \land {min}_{i} H_{i} (t_{0}) \end{matrix}

(26)

and the pointwise Hölder exponent of

X

t_{0}

satisfies almost surely

\begin{matrix} α (t_{0}) = β (t_{0}) \land {min}_{i} H_{i} (t_{0}) \end{matrix}

(27)

where

β (t_{0})

and

\tilde{β} (t_{0})

denote the pointwise and local Hölder exponents of

H

t_{0}

The proofs of propositions 13 and 14 are detailed in the three following paragraphs.

4.2.1 Lower bound for the local Hölder exponent

A lower bound for the local Hölder exponent is directly given by Kolmogorov's theorem. Indeed, for

X

a multifractional Brownian field or a multifractional Brownian sheet indexed by

[a, b]

, for all

n \in N

, there exists

λ_{n} > 0

such that

E {[X_{t} - X_{s}]}^{2 n} \leq λ_{n} ∥ t - s ∥^{n . α}

with

α = 2 {inf}_{[a, b]} (\tilde{β} \land H)

α = 2 {inf}_{[a, b]} (\tilde{β} \land {min}_{i} H_{i})

Kolmogorov's theorem states that there exists a modification of

X

, which is

q

-Hölder continuous for all

q \in (0, \frac{α}{2} - \frac{N}{2 n})

. Then, for all

t_{0} \in R_{+}^{N}

and all

a, b \in R_{+}^{N}

such that

a ≺ b

and

t_{0} \in (a, b)

, we have

\forall n \in N; \tilde{α} (t_{0}) \geq \frac{α}{2} - \frac{N}{2 n}

and therefore, taking the limit

n \to \infty

\tilde{α} (t_{0}) \geq \frac{α}{2}

H

is continuous, we can take the limit

(a, b) \to (t_{0}, t_{0})

and we get

∙ in the isotropic case,

\begin{matrix} \tilde{α} (t_{0}) \geq \tilde{β} (t_{0}) \land H (t_{0}) \end{matrix}

(28)

∙ in the non-isotropic case,

\begin{matrix} \tilde{α} (t_{0}) \geq \tilde{β} (t_{0}) \land {min}_{i} H_{i} (t_{0}) \end{matrix}

(29)

4.2.2 Lower bound for the pointwise Hölder exponent

By ( 23 ), paragraph 4.2.1 provides a lower bound for the pointwise Hölder exponent. However, it can be improved in the case

\tilde{β} (t_{0}) < β (t_{0})

Let

X = {X_{t}; t \in R_{+}^{N}}

be a multifractional Brownian field. By corollary 2 , there exist positive constants

K

and

L

such that for all

s, t \in R_{+}^{N}

E {[X_{t} - X_{s}]}^{2} \leq K ∥ t - s ∥^{2 H (t)} + L | H (t) - H (s) |^{2}

and by corollary 3 , there exists positive constants

α

and

M

such that

\forall s, t \in [a, b]; E {[X_{t} - X_{s}]}^{2} \leq M ∥ t - s ∥^{α}

Therefore, using Kolmogorov's criterion, there exists a modification of

X

, which is

ν

-Hölder continuous for all

ν \in] 0, \frac{α}{2} [

. In the following, we consider such a

ν

with

\frac{1}{ν} \in N

For all

ε > 0

, there exist

ρ_{0} > 0

and

M > 0

such that for all

ρ < ρ_{0}

and all

t \in B (t_{0}, ρ)

E {[\frac{X_{t} - X_{t_{0}}}{ρ^{β (t_{0}) \land H (t_{0}) - ε}}]}^{2} \leq M ρ^{ε}

Then, setting

γ = β (t_{0}) \land H (t_{0}) - ε

, for all

p \in N^{*}

\begin{matrix} P {| X_{t} - X_{t_{0}} | > ρ^{γ}} \leq E {[\frac{X_{t} - X_{t_{0}}}{ρ^{γ}}]}^{2 p} \leq M_{p} ρ^{p ε} \end{matrix}

Let

ρ = 2^{- n}

and for all

m \in N

D_{m} = {t_{0} + k . 2^{- (n + m)}; k \in {0, \pm 1, . . ., \pm 2^{m}}^{N}}

let us compute

\begin{matrix} P {{max}_{k \in {\pm 1, . . ., \pm 2^{m}}^{N}} \frac{| X_{t_{0} + k . 2^{- (m + n)}} - X_{t_{0}} |}{2^{- γ n}} > 1} \end{matrix}

\begin{matrix} \leq \sum_{k \in {\pm 1, . . ., \pm 2^{m}}^{N}} P {| X_{t_{0} + k . 2^{- (n + m)}} - X_{t_{0}} | > 2^{- γ n}} \end{matrix}

\begin{matrix} \leq M_{p} 2^{(m + 1) N} 2^{- p ε n} \end{matrix}

Let us take

m = \frac{1 + ⌊ γ ⌋}{ν} n = κ n

and

p \in N

such that

N \frac{1 + ⌊ γ ⌋}{ν} - p ε < 0

. By the Borel-Cantelli lemma, there exists a finite random variable

n^{*}

such that almost surely,

\begin{matrix} \forall n \geq n^{*}; {max}_{k \in {0, . . ., \pm 2^{κ n}}^{N}} | X_{t_{0} + k . 2^{- (1 + κ) n}} - X_{t_{0}} | \leq 2^{- γ n} \end{matrix}

(30)

From ( 30 ), we show that, almost surely, for all

m \in N

, we have

\begin{matrix} \forall t \in D_{m}; | X_{t} - X_{t_{0}} | \leq C 2^{- γ n} \end{matrix}

(31)

∙ if $0 \leq m \leq κ n$ , ( 31 ) follows directly from ( 30 )

∙ if

m > κ n

, for

t \in D_{m}

, let

C_{t_{0}, t}^{κ n} = {x \in D_{κ n}; \forall i, (t_{0})_{i} \leq x_{i} \leq t_{i}}

Then consider

\hat{t} \in B (t, 2^{- (1 + κ) n}) \cap C_{t_{0}, t}^{κ n}

As the paths of

X

are

ν

-Hölder continuous, we have

\begin{matrix} | X_{\hat{t}} - X_{t} | \leq \tilde{C} 2^{- ν (1 + κ) n} \leq \tilde{C} 2^{- γ n} \end{matrix}

and by ( 30 ),

| X_{\hat{t}} - X_{t_{0}} | \leq 2^{- γ n}

Using the triangular inequality, the result follows.

Therefore, ( 31 ) leads to

\begin{matrix} \forall m \in N; \forall s, t \in D_{m}; | X_{t} - X_{s} | \leq 2 C 2^{- γ n} \end{matrix}

Using the continuity of

X

and

m \to + \infty

, we get

{sup}_{s, t \in B (t_{0}, 2^{- n})} | X_{t} - X_{s} | \leq 2 C 2^{- γ n}

and therefore, almost surely,

\begin{matrix} {limsup}_{ρ \to 0} {sup}_{s, t \in B (t_{0}, ρ)} \frac{| X_{t} - X_{s} |}{ρ^{γ}} < + \infty \end{matrix}

(32)

By ( 32 ), for all

ε > 0

, almost surely

α (t_{0}) \geq β (t_{0}) \land H (t_{0}) - ε

Taking

ε \in Q_{+}

, we have almost surely

\begin{matrix} α (t_{0}) \geq β (t_{0}) \land H (t_{0}) \end{matrix}

(33)

For a multifractional Brownian sheet

X

, by lemma 2 , we get in the same way that, almost surely

\begin{matrix} α (t_{0}) \geq β (t_{0}) \land H_{i} (t_{0}) \end{matrix}

(34)

for all

i = 1, . . ., N

4.2.3 Upper bound for the pointwise Hölder exponent

The main result getting the upper bound for the Hölder exponents, is the following lemma, a direct consequence of proposition 12 using continuity of

D

D^{'}

and

D^{''}

Lemma 3 Let

X = {X_{t}; t \in R_{+}^{N}}

be a multifractional Brownian field. For all

[a, b] \subset R_{+}^{N}

, there exist positive constants

k_{1}

k_{2}

l_{1}

l_{2}

such that

\begin{matrix} \forall s, t \in [a, b]; & E {[X_{t} - X_{s}]}^{2} \geq k_{1} ∥ t - s ∥^{2 H (t)} - l_{1} {(H (t) - H (s))}^{2} \end{matrix}

(35)

\begin{matrix} E {[X_{t} - X_{s}]}^{2} \geq k_{2} {(H (t) - H (s))}^{2} - l_{2} ∥ t - s ∥^{2 H (t)} \end{matrix}

(36)

Proof

We only have to study the multiplicative factors of

∥ t - s ∥^{2 H (t)}

and

{(H (t) - H (s))}^{2}

in ( 12 )

∙ Let $k_{1} = {inf}_{t \in [a, b]} D [2 H (t)]$ and $l_{2} = {sup}_{t \in [a, b]} D [2 H (t)]$ .
By continuity of $t \mapsto D [2 H (t)]$ on the compact $[a, b]$ and as the function $D$ is positive (lemma 1 ), for all $t \in [a, b]$ $0 < k_{1} \leq D [2 H (t)] \leq l_{2} < + \infty$

∙ and let

Φ (t) = ∥ t ∥^{2 H (t)} \times {D [2 H (t)] {ln}^{2} ∥ t ∥ - 2 D^{'} [2 H (t)] ln ∥ t ∥ + D^{''} [2 H (t)]}

By lemma 1 ,

\begin{matrix} Φ (t) = \int_{R^{N}} {\underset{︸}{\frac{1 - e^{i u_{1}}}{∥ u ∥^{2 H (t) + N}} {(ln ∥ t ∥ - ln ∥ u ∥)}^{2}}}_{\geq 0} . d u \end{matrix}

Let

k_{2} = {inf}_{t \in [a, b]} Φ (t)

and

l_{1} = {sup}_{t \in [a, b]} Φ (t)

As previously, for all

t \in [a, b]

, we have

0 < k_{2} \leq Φ (t) \leq l_{1} < + \infty

□

Lemma 4 Let

X = {X_{t}; t \in R_{+}^{N}}

be a multifractional Brownian sheet. For all

[a, b] \subset R_{+}^{N}

, there exist positive constants

k_{1}

k_{2}

l_{1}

l_{2}

such that

\begin{matrix} \forall s, t \in [a, b]; & t - s \in R_{+} . ε_{i} \end{matrix}

\begin{matrix} E {[X_{t} - X_{s}]}^{2} \geq k_{1} ∥ t - s ∥^{2 H_{i} (t)} - l_{1} {(H_{i} (t) - H_{i} (s))}^{2} \end{matrix}

(37)

\begin{matrix} E {[X_{t} - X_{s}]}^{2} \geq k_{2} {(H_{i} (t) - H_{i} (s))}^{2} - l_{2} ∥ t - s ∥^{2 H_{i} (t)} \end{matrix}

(38)

Proof

For all

s

t

such that

t - s \in R_{+} . ε_{i}

, using lemma 3 , we have

\begin{matrix} E {[X_{t} - X_{s}]}^{2} & = & E {[X_{t^{(i)}}^{(i)} - X_{s^{(i)}}^{(i)}]}^{2} \prod_{j \neq i} E {[X_{t^{(j)}}^{(j)}]}^{2} \end{matrix}

\begin{matrix} \geq & k_{1} | t_{i} - s_{i} |^{2 H_{i} (t)} - l_{1} {(H_{i} (t) - H_{i} (s))}^{2} \end{matrix}

and

\begin{matrix} E {[X_{t} - X_{s}]}^{2} \geq k_{2} {(H_{i} (t) - H_{i} (s))}^{2} - l_{2} | t_{i} - s_{i} |^{2 H_{i} (t)} \end{matrix}

□

From this result, the upper bound for the pointwise exponent is a consequence of the following lemma whose proof is the same as the case

N = 1

(see [1] )

Lemma 5 Let

X = {X_{t}; t \in R_{+}^{N}}

be a Gaussian process. Assume there exists

μ \in (0, 1)

such that for all

ε > 0

, there exist a sequence

{(h_{n})}_{n \in N}

{(R_{+}^{N})}^{*}

converging to 0, and a constant

c > 0

such that

\forall n \in N; E {[X_{t + h_{n}} - X_{t}]}^{2} \geq c ∥ h_{n} ∥^{2 μ + ε}

Then we have almost surely

α (t) \leq μ

Let

X = {X_{t}; t \in R_{+}^{N}}

be a multifractional Brownian field (resp. multifractional Brownian sheet). Let

β (t_{0})

be the pointwise Hölder exponent of

H

t_{0}

We consider the two cases :

∙ if

H (t_{0}) < β (t_{0})

(resp.

H_{i} (t_{0}) < β (t_{0})

), by definition of

β (t_{0})

, we have

\begin{matrix} {lim}_{h \to 0} \frac{∥ H (t_{0} + h) - H (t_{0}) ∥}{∥ h ∥^{H (t_{0})}} = 0 \end{matrix}

Hence, by ( 35 ) (resp. ( 37 )), there exists a positive constant

C

such that

\begin{matrix} E {[X_{t_{0} + h} - X_{t_{0}}]}^{2} \geq C ∥ h ∥^{2 H (t_{0})} \end{matrix}

Then, by lemma 5

\begin{matrix} α (t_{0}) \leq H (t_{0}) (resp. H_{i} (t_{0})) \end{matrix}

(39)

∙ if

H (t_{0}) > β (t_{0})

(resp.

H_{i} (t_{0}) > β (t_{0})

), we consider

α \in (β (t_{0}); H (t_{0}))

(resp.

α \in (β (t_{0}); H_{i} (t_{0}))

). There exists a positive constant

C

and a sequence

{(h_{n})}_{n \in N}

converging to

0

such that

\begin{matrix} \forall n \in N; ∥ H (t_{0} + h_{n}) - H (t_{0}) ∥ > C ∥ h_{n} ∥^{α} \end{matrix}

Then, by ( 36 ) (resp. ( 38 ))

\begin{matrix} \forall n \in N; E {[X_{t_{0} + n} - X_{t_{0}}]}^{2} & > & k_{2} C ∥ h_{n} ∥^{2 α} - l_{2} ∥ h_{n} ∥^{2 H (t_{0})} \end{matrix}

\begin{matrix} \geq & C^{'} ∥ h_{n} ∥^{2 α} \end{matrix}

hence, by lemma 5

α \geq α (t_{0})

and therefore

\begin{matrix} α (t_{0}) \leq β (t_{0}) \end{matrix}

(40)

We can restate the upper bounds ( 39 ) and ( 40 ) of the pointwise Hölder exponent of

X

t_{0}

\begin{matrix} α (t_{0}) \leq β (t_{0}) \land H (t_{0}) (resp. β (t_{0}) \land H_{i} (t_{0})) \end{matrix}

(41)

4.2.4 Upper bound for the local Hölder exponent

By ( 23 ), any upper bound for the pointwise Hölder exponent is an upper bound for the local Hölder exponent. But we can improve on this result in the case

\tilde{β} (t_{0}) < H (t_{0})

. We first give an analogous of lemma 5 for the local exponent

Lemma 6 Let

X = {X_{t}; t \in R_{+}^{N}}

be a Gaussian process. Assume there exists

μ \in (0, 1)

such that for all

ε > 0

, there exist two sequences

{(h_{n})}_{n \in N}

and

{(l_{n})}_{n \in N}

{(R_{+}^{N})}^{*}

converging to 0, and a constant

c > 0

such that

\forall n \in N; E {[X_{t_{0} + h_{n}} - X_{t_{0} + l_{n}}]}^{2} \geq c ∥ h_{n} - l_{n} ∥^{2 μ + ε}

Then we have almost surely

\tilde{α} (t_{0}) \leq μ

Proof

Let

ε > 0

and consider two sequences

(h_{n})_{n \in N}

and

(l_{n})_{n \in N}

as in the statement.

For all

n \in N

, the law of the random variable

\frac{X_{t_{0} + h_{n}} - X_{t_{0} + l_{n}}}{∥ h_{n} - l_{n} ∥^{μ + ε}}

N (0, σ_{n}^{2})

From the assumption, we have

σ_{n} \to + \infty

n \to + \infty

Then, for all

λ > 0

\begin{matrix} P {\frac{∥ h_{n} - l_{n} ∥^{μ + ε}}{| X_{t_{0} + h_{n}} - X_{t_{0} + l_{n}} |} < λ} & = & P {\frac{| X_{t_{0} + h_{n}} - X_{t_{0} + l_{n}} |}{∥ h_{n} - l_{n} ∥^{μ + ε}} > \frac{1}{λ}} \end{matrix}

\begin{matrix} = & \int_{| x | > \frac{1}{λ}} \frac{1}{\sqrt{2 π} σ_{n}} exp (- \frac{x^{2}}{2 σ_{n}^{2}}) . d x \end{matrix}

\begin{matrix} = & \frac{1}{2 π} \int_{| x | > \frac{1}{λ σ_{n}}} exp (- \frac{x^{2}}{2}) . d x \overset{n \to + \infty}{⟶} 1 \end{matrix}

Therefore the sequence

{(\frac{∥ h_{n} - l_{n} ∥^{μ + ε}}{| X_{t_{0} + h_{n}} - X_{t_{0} + l_{n}} |})}_{n \in N}

converges to

0

in probability.

then there exists a subsequence which converges to

0

almost surely. Then we have almost surely

\tilde{α} (t_{0}) \leq μ + ε

. Taking

ε \in Q_{+}

, the result follows.

□

Let

α \in (\tilde{β} (t_{0}); H (t_{0}))

(resp.

α \in (\tilde{β} (t_{0}); H_{i} (t_{0}))

). As

\begin{matrix} {limsup}_{ρ \to 0} {sup}_{s, t \in B (t_{0}, ρ)} \frac{| H (t) - H (s) |}{∥ t - s ∥^{α}} = + \infty \end{matrix}

for all

M > 0

, there exists

ρ_{0} > 0

such that

\begin{matrix} \forall ρ < ρ_{0}; \exists s, t \in B (t_{0}, ρ); | H (t) - H (s) | > M ∥ t - s ∥^{α} \end{matrix}

Therefore we can construct two sequences

(h_{n})

and

(l_{n})

converging to

0

such that

\begin{matrix} \forall n \in N; | H (t_{0} + h_{n}) - H (t_{0} + l_{n}) | > M ∥ h_{n} - l_{n} ∥^{α} \end{matrix}

By lemma 6 , we can deduce

\begin{matrix} \tilde{α} (t_{0}) \leq \tilde{β} (t_{0}) \end{matrix}

(42)

4.3 Directional Hölder exponents

One may also define directional pointwise and local Hölder exponents in the direction

u \in U = {u \in R^{N}; ∥ u ∥ = 1}

\begin{matrix} α_{u} (t_{0}) & = & sup {α; {lim}_{ρ \to 0} \frac{| X_{t_{0} + ρ . u} - X_{t_{0}} |}{ρ^{α}} = 0} \end{matrix}

and

\begin{matrix} {\tilde{α}}_{u} (t_{0}) = sup {α; {limsup}_{ρ \to 0} {sup}_{\binom{s, t \in B (t_{0}, ρ)}{s, t \in t_{0} + R . u}} \frac{| X_{t} - X_{s} |}{∥ t - s ∥^{α}} < \infty} \end{matrix}

As previously, for all

u \in U

, we have

\begin{matrix} {\tilde{α}}_{u} (t_{0}) \leq α_{u} (t_{0}) \end{matrix}

(43)

Moreover, we can see easily that for all

u \in U

, we have

\begin{matrix} α (t_{0}) \leq α_{u} (t_{0}) and \tilde{α} (t_{0}) \leq {\tilde{α}}_{u} (t_{0}) \end{matrix}

(44)

Proposition 15 Let

X = {X_{t}; t \in R_{+}^{N}}

be a multifractional Brownian field.

For all

t_{0} \in R_{+}^{N}

and all

u \in U

, the local Hölder exponent of

X

t_{0}

in the direction

u

is almost surely given by

\begin{matrix} {\tilde{α}}_{u} (t_{0}) = {\tilde{β}}_{u} (t_{0}) \land H (t_{0}) \end{matrix}

(45)

and the pointwise Hölder exponent of

X

t_{0}

in the direction

u

satisfies almost surely

\begin{matrix} α_{u} (t_{0}) = β_{u} (t_{0}) \land H (t_{0}) \end{matrix}

(46)

where

β_{u} (t_{0})

and

{\tilde{β}}_{u} (t_{0})

denote the pointwise and local Hölder exponents of

H

t_{0}

in the direction

u

Proof

Let

t_{0} \in R_{+}^{N}

u \in U

and consider the stochastic process

\tilde{X} = {{\tilde{X}}_{ρ} = X_{t_{0} + ρ . u}; ρ > 0}

By definition,

α_{u} (t_{0})

and

{\tilde{α}}_{u} (t_{0})

are respectively the pointwise and local Hölder exponents of

\tilde{X}

0

Let

\tilde{H} (ρ) = H (t_{0} + ρ . u)

. We have

\begin{matrix} \frac{1}{2} E {[{\tilde{X}}_{ρ} - {\tilde{X}}_{η}]}^{2} & = & \frac{1}{2} E {[X_{t_{0} + ρ . u} - X_{t_{0} + η . u}]}^{2} \end{matrix}

\begin{matrix} = & D [2 \tilde{H} (ρ)] | ρ - η |^{2 \tilde{H} (ρ)} + \frac{\partial^{2} φ}{\partial x^{2}} (2 \tilde{H} (ρ); ∥ t_{0} + ρ . u ∥) \times {(\tilde{H} (ρ) - \tilde{H} (η))}^{2} \end{matrix}

\begin{matrix} + o {(\tilde{H} (ρ) - \tilde{H} (η))}^{2} + o (| ρ - η |^{2 \tilde{H} (ρ)}) \end{matrix}

by corollary 1 .

Then, using the same method as in proposition 13 , the result follows.

□

Proposition 16 Let

X = {X_{t}; t \in R_{+}^{N}}

be a multifractional Brownian sheet. For all

t_{0} \in R_{+}^{N}

, the local Hölder exponent of

X

t_{0}

in the direction

ε_{i}

is almost surely given by

\begin{matrix} {\tilde{α}}_{ε_{i}} (t_{0}) = {\tilde{β}}_{ε_{i}} (t_{0}) \land H_{i} (t_{0}) \end{matrix}

(47)

and the pointwise Hölder exponent of

X

t_{0}

in the direction

u

satisfies almost surely

\begin{matrix} α_{ε_{i}} (t_{0}) = β_{ε_{i}} (t_{0}) \land H_{i} (t_{0}) \end{matrix}

(48)

where

β_{ε_{i}} (t_{0})

and

{\tilde{β}}_{ε_{i}} (t_{0})

denote the pointwise and local Hölder exponents of

H

t_{0}

in the direction

ε_{i}

Proof

As in the proof of lemma 2 , there exists a constant

M > 0

and a one-parameter mBm

X^{(i)}

such that

\begin{matrix} E {[X_{t_{0} + ρ . ε_{i}} - X_{t_{0} + η . ε_{i}}]}^{2} = M . E {[X_{t_{0}^{(i)} + ρ . ε_{i}}^{(i)} - X_{t_{0}^{(i)} + η . ε_{i}}^{(i)}]}^{2} \end{matrix}

Then, using the same method as in proposition 14 , the result follows.

□

4.4 Application of Dudley's theory

Another way to study the regularity of our processes is to examine the behavior around zero of the modulus of continuity

ω_{X, T} (δ) = {sup}_{s, t \in T; d (s, t) \leq δ} | X_{s} - X_{t} |

When the process studied is Gaussian, it is convenient to consider the pseudo-metric

d (s, t) = E {[X_{s} - X_{t}]}^{2}

As usually, we define the ball of radius

r > 0

about

t \in T

ℬ_{d} (t, r) = {s \in T; d (s, t) < r}

and we say that

(T, d)

is totally bounded if for all

ε > 0

, there exists

t_{1}, . . ., t_{m} \in T

such that

T \subset^{m} ⋃_{i = 1} ℬ_{d} (t_{i}, ε)

When

(T, d)

is totally bounded, we can define the metric entropy

ε \mapsto D (ε, T, d)

where

D (ε, T, d)

is minimum number of balls of radius

ε

required to cover

T

The following theorem allows to improve on the results of the previous paragraph.

Theorem 2 (Dudley's Theorem) Consider a centered Gaussian process

X = {X_{t}; t \in R_{+}^{N}}

indexed by the pseudo-metric space

(T, d)

. If

(T, d)

is totally bounded and if

\int_{0}^{1} \sqrt{ln D (r, T, d)} d r < \infty

, then

X

has a continuous modification

Y = {Y_{t}; t \in R_{+}^{N}}

. Moreover, there exists a universal constant

C > 0

such that

{limsup}_{δ \to 0^{+}} \frac{ω_{Y, T} (δ)}{\int_{0}^{δ} \sqrt{ln D (\frac{r}{2}, T, d)} d r + C δ \sqrt{ln ln \frac{1}{δ}}} \leq 24

To apply this result, we first need to verify the assumptions about the metric entropy.

Lemma 7 Let

T \subset R_{+}^{N}

measurable and

d

a pseudo-metric on

T

. If there exists

C > 0

and

α > 0

such that

\forall s, t \in T; d (s, t) \leq C ∥ s - t ∥^{α}

then there exists

r_{0} > 0

such that

\forall r \in [0, r_{0}]

D (r, T, d) \leq C^{\frac{N}{α}} . L e b (T) . r^{- \frac{N}{α}}

We saw previously that the

2

multi-parameter extensions of the mBm, satisfy the assumption of this lemma with

T = [a, b]

. Then there exists

C_{a, b} > 0

α = α_{a, b} > 0

and

r_{0} > 0

such that for all

r \in [0, t_{0}]

D (r, [a, b], d) \leq C^{\frac{N}{α}} . L e b ([a, b]) . r^{- \frac{N}{α}}

As a consequence,

([a, b], d)

is totally bounded and in the neighborhood of

0

, we have

\sqrt{ln D (r, T, d)} \leq \sqrt{K - \frac{N}{α} ln r}

therefore the integral

\int_{0}^{1} \sqrt{ln D (r, T, d)} d r

is finite and we can apply Dudley's theorem. We get

Proposition 17 Let

X = {X_{t}; t \in R_{+}^{N}}

be one of the multi-parameter extension of the mBm. For all

a ≺ b

, there exists

C_{a, b} > 0

and

α = α_{a, b} > 0

such that

{limsup}_{ε \to 0} \frac{{sup}_{∥ s - t ∥ \leq ε} | X_{s} - X_{t} |}{ε^{α} \sqrt{ln \frac{1}{ε}}} \leq C_{a, b} \sqrt{N}

Proof

First of all, we study the quotient

\frac{\int_{0}^{δ} \sqrt{ln D (\frac{r}{2}, T, d)} d r + C δ \sqrt{ln ln \frac{1}{δ}}}{δ \sqrt{ln \frac{1}{δ}}} = \frac{N (δ)}{D (δ)}

The derivative of the numerator is

\begin{matrix} N^{'} (δ) & = & \sqrt{ln D (\frac{δ}{2}, T, d)} + C [\sqrt{ln ln \frac{1}{δ}} - \frac{1}{2 ln \frac{1}{δ} \sqrt{ln ln \frac{1}{δ}}}] \end{matrix}

\begin{matrix} \leq & \sqrt{K - \frac{N}{α} ln δ} + C \sqrt{ln ln \frac{1}{δ}} [1 - \frac{1}{2 ln \frac{1}{δ} ln ln \frac{1}{δ}}] \end{matrix}

and the derivative of the denominator is

\begin{matrix} D^{'} (δ) & = & \sqrt{ln \frac{1}{δ}} [1 - \frac{1}{2 ln \frac{1}{δ}}] \end{matrix}

\begin{matrix} \sim & \sqrt{ln \frac{1}{δ}} \end{matrix}

Then we have

{limsup}_{δ \to 0^{+}} \frac{N^{'} (δ)}{D^{'} (δ)} \leq \sqrt{\frac{N}{α}}

and by a L'Hopital's rule type argument,

{limsup}_{δ \to 0^{+}} \frac{N (δ)}{D (δ)} \leq \sqrt{\frac{N}{α}}

Then we have

{limsup}_{δ \to 0} \frac{{sup}_{d (s, t) \leq δ} | X_{s} - X_{t} |}{δ \sqrt{ln \frac{1}{δ}}} \leq 24 \sqrt{\frac{N}{α}}

The problem is now to transform

{sup}_{d (s, t) \leq δ}

into

{sup}_{∥ s - t ∥ \leq δ}

To do this, we write

\forall s, t \in [a, b]; d (s, t) \leq C_{a, b} ∥ s - t ∥^{α}

then

∥ s - t ∥ \leq ε = \frac{δ^{\frac{1}{α}}}{C_{a, b}^{\frac{1}{α}}}

implies

d (s, t) \leq δ

and we get

{limsup}_{δ \to 0} \frac{{sup}_{∥ s - t ∥ \leq ε} | X_{s} - X_{t} |}{C_{a, b} ε^{α} \sqrt{α ln \frac{1}{ε}}} \leq 24 \sqrt{\frac{N}{α}}

which gives the expected result.

□

This result is more powerful than knowledge of Hölder exponents. It gives the behavior of

| X_{s} - X_{t} |

in a ball around

t_{0}

5 Locally asymptotic self-similarity

Extending fBm into multifractional Brownian motion implies the loss of the two properties of self-similarity and stationarity of increments. However, a weak form of self-similarity remains, called locally asymptotic self-similarity (see [1] , [4] ). As we will see, this property still holds for the two kinds of extension of mBm in

R^{N}

Theorem 3 Let

X = {X_{t}; t \in R_{+}^{N}}

be a multifractional Brownian field.

For all

t_{0} \in R_{+}^{N}

, the law of the process

Y^{α} (ρ) = {Y_{u}^{α} (ρ) = \frac{X_{t_{0} + ρ u} - X_{t_{0}}}{ρ^{α}}; u \in R_{+}^{N}}

converges weakly if one of the following two conditions holds

1. $α = H (t_{0})$ and $H (t_{0}) < {inf}_{u, v} β_{u v} (t_{0})$ where $β_{u v} (t_{0}) = sup {α; {lim}_{ρ \to 0} \frac{| H (t_{0} + ρ u) - H (t_{0} + ρ v) |}{ρ^{α}} = 0}$ .
Then, the limit measure is the law of a fractional Brownian field with parameter $H (t_{0})$ .
2. $α = {inf}_{u, v} β_{u v} (t_{0})$ , $H (t_{0}) > {inf}_{u, v} β_{u v} (t_{0})$ and for all $u, v \in R_{+}^{N}$ , the following limit exists ${lim}_{ρ \to 0} \frac{| H (t_{0} + ρ u) - H (t_{0} + ρ v) |}{ρ^{{inf}_{u, v} β_{u v} (t_{0})}} = Γ (u, v)$ with $(u, v) \mapsto \frac{Γ (u, v)}{∥ u - v ∥^{2 β}}$ bounded on $[a, b]^{2}$ for some $β > 0$ .
The limit measure is the law of a Gaussian process $Y^{{inf}_{u, v} β_{u v} (t_{0})}$ such that $E {[Y_{u}^{{inf}_{u, v} β_{u v} (t_{0})} - Y_{v}^{{inf}_{u, v} β_{u v} (t_{0})}]}^{2} = K_{t_{0}} {[Γ (u, v)]}^{2}$

Remark 5 As in the Levy fBm's case in proposition 6 , the same result as theorem 3 can be stated for the increments

Δ X

defined in section 2.3 . The law of the process

Y^{α} (ρ) = {Y_{u}^{α} (ρ) = \frac{Δ X_{t_{0}, t_{0} + ρ u}}{ρ^{α}}; u \in R_{+}^{N}}

converges weakly under the same assumptions.

In the case

N = 1

, for all

u, v \in R_{+}

, we have

β_{u v} (t_{0}) = β (t_{0})

. Therefore, theorem 3 has a simpler statement. The two cases to be considered, depend of the comparison between

H (t_{0})

and the pointwise exponent

β (t_{0})

H

The following example shows that the limit considered in the second case, can be non trivial.

Example 1 In the case

N = 1

, let

H (t) = \frac{3}{4} + t^{\frac{1}{2}}

for

t \in [0, \frac{1}{4}]

For

t_{0} = 0

, we compute, for all

u, v

and

ρ > 0

\frac{| H (ρ . u) - H (ρ . v) |}{ρ^{\frac{1}{2}}} = | u^{\frac{1}{2}} - v^{\frac{1}{2}} | < | u - v |^{\frac{1}{2}}

The limit measure is the law of a centered Gaussian process

Y

such that

E {[Y_{u} - Y_{v}]}^{2} = K_{0} {(u^{\frac{1}{2}} - v^{\frac{1}{2}})}^{2}

E [Y_{u} Y_{v}] = K_{0} u^{\frac{1}{2}} v^{\frac{1}{2}}

Theorem 4 Let

X = {X_{t}; t \in R_{+}^{N}}

be a multifractional Brownian sheet.

The law of the process

Y^{α} (ρ) = {Y_{u}^{α} (ρ) = \frac{Δ X_{t_{0}, t_{0} + ρ u}}{ρ^{\sum_{i} α_{i}}}; u \in R_{+}^{N}}

converges weakly if for all

i \in {1, . . ., N}

, one of the following two conditions holds

1. $α_{i} = H_{i} (t_{0})$ and $H_{i} (t_{0}) < {inf}_{u, v} β_{u v}^{i} (t_{0})$ where $β_{u v}^{i} (t_{0}) = sup {α; {lim}_{ρ \to 0} \frac{| H_{i} (t_{0} + ρ u) - H_{i} (t_{0} + ρ v) |}{ρ^{α}} = 0}$ .
2. $α_{i} = {inf}_{u, v} β_{u v}^{i} (t_{0})$ , $H_{i} (t_{0}) > {inf}_{u, v} β_{u v}^{i} (t_{0})$ and ${lim}_{ρ \to 0} \frac{| H_{i} (t_{0} + ρ u) - H_{i} (t_{0} + ρ v) |}{ρ^{{inf}_{u, v} β_{u v}^{i} (t_{0})}} = Γ_{i} (u, v)$ with $(u, v) \mapsto \frac{Γ_{i} (u, v)}{∥ u - v ∥^{2 β_{i}}}$ bounded on $[a, b]^{2}$ for some $β_{i} > 0$ .

As usually, the proof of weak convergence proceeds in two steps. First, we need to show finite dimensional convergence, and then, use a tightness argument.

Lemma 14.2 and theorem 14.3 in [10] , for instance, allow then to conclude.

5.1 Finite dimensional convergence

As the considered processes are Gaussian, we only have to show the convergence of covariance functions.

5.1.1 Multifractional Brownian field

By ( 12 ), we compute

\begin{matrix} ρ^{2 α} E {[Y_{u}^{α} (ρ) - Y_{v}^{α} (ρ)]}^{2} & = & E {[X_{t_{0} + ρ u} - X_{t_{0} + ρ v}]}^{2} \end{matrix}

\begin{matrix} = & D [H (t_{0} + ρ u) + H (t_{0} + ρ v)] \times ∥ ρ . (u - v) ∥^{H (t_{0} + ρ u) + H (t_{0} + ρ v)} \end{matrix}

\begin{matrix} + \frac{\partial^{2} φ}{\partial x^{2}} (2 H (t_{0} + ρ u); ∥ t_{0} + ρ u ∥) \times {(H (t_{0} + ρ u) - H (t_{0} + ρ v))}^{2} \end{matrix}

\begin{matrix} + o (∥ ρ . (u - v) ∥^{2}) + o {(H (t_{0} + ρ u) - H (t_{0} + ρ v))}^{2} \end{matrix}

(49)

To show that

ρ^{H (t_{0} + ρ u) + H (t_{0} + ρ v)} \sim ρ^{2 H (t_{0})}

in the neighborhood of

ρ = 0

, we study

\begin{matrix} [H (t_{0} + ρ u) + H (t_{0} + ρ v) - 2 H (t_{0})] ln ρ & = & \frac{H (t_{0} + ρ u) - H (t_{0})}{∥ ρ . u ∥^{α}} \times ∥ ρ . u ∥^{α} ln ρ \end{matrix}

\begin{matrix} + \frac{H (t_{0} + ρ v) - H (t_{0})}{∥ ρ . v ∥^{α}} \times ∥ ρ . v ∥^{α} ln ρ \end{matrix}

for

α < β (t_{0})

(u; ρ) \mapsto ∥ ρ . u ∥^{α} ln ρ

is bounded on

[a, b] \times [0, 1]

and

\forall u \in [a, b]; \frac{H (t_{0} + ρ u) - H (t_{0})}{∥ ρ . u ∥^{α}} \overset{ρ \to 0}{⟶} 0

we have

[H (t_{0} + ρ u) + H (t_{0} + ρ v) - 2 H (t_{0})] ln ρ \overset{ρ \to 0}{⟶} 0

Therefore, in the neighborhood of

ρ = 0

, the first term of ( 49 ) is equivalent to

D [2 H (t_{0})] ∥ u - v ∥^{2 H (t_{0})} \times ρ^{2 H (t_{0})}

and the second to

\frac{\partial^{2} φ}{\partial x^{2}} (2 H (t_{0}); ∥ t_{0} ∥) \times {(H (t_{0} + ρ u) - H (t_{0} + ρ v))}^{2}

Let

β_{u v} (t_{0}) = sup {α; {lim}_{ρ \to 0} \frac{| H (t_{0} + ρ u) - H (t_{0} + ρ v) |}{ρ^{α}} = 0}

. We have to distinguish the two following cases

∙ if

H (t_{0}) < {inf}_{u, v} β_{u v} (t_{0})

, by definition of

β_{u v} (t_{0})

\begin{matrix} \forall u, v \in R_{+}^{N}; {lim}_{ρ \to 0} \frac{| H (t_{0} + ρ u) - H (t_{0} + ρ v) |}{ρ^{H (t_{0})}} = 0 \end{matrix}

Therefore

\begin{matrix} \forall u, v \in R_{+}^{N}; E {[Y_{u}^{H (t_{0})} (ρ) - Y_{v}^{H (t_{0})} (ρ)]}^{2} \overset{ρ \to 0}{⟶} {\underset{︸}{D [2 H (t_{0})] ∥ u - v ∥^{2 H (t_{0})}}}_{E {[B_{u}^{H (t_{0})} - B_{v}^{H (t_{0})}]}^{2}} \end{matrix}

where

B^{H (t_{0})}

denotes fractional Brownian field of parameter

H (t_{0})

∙ if

H (t_{0}) > {inf}_{u, v} β_{u v} (t_{0})

, for all

α < {inf}_{u, v} β_{u v} (t_{0})

, as

\begin{matrix} \forall u, v \in R_{+}^{N}; {lim}_{ρ \to 0} \frac{| H (t_{0} + ρ u) - H (t_{0} + ρ v) |}{ρ^{α}} = 0 \end{matrix}

we have

\begin{matrix} \forall u, v \in R_{+}^{N}; \frac{1}{ρ^{2 α}} E {[X_{t_{0} + ρ u} - X_{t_{0} + ρ v}]}^{2} \overset{ρ \to 0}{⟶} 0 \end{matrix}

Moreover, since there exists

u, v \in R_{+}^{N}

such that

H (t_{0}) > β_{u v} (t_{0})

, we can consider

α \in (β_{u v} (t_{0}); H (t_{0}))

. The limit

{limsup}_{ρ \to 0} \frac{| H (t_{0} + ρ u) - H (t_{0} + ρ v) |}{ρ^{α}} = + \infty

implies

\begin{matrix} \exists u, v \in R_{+}^{N}; {limsup}_{ρ \to 0} \frac{1}{ρ^{2 α}} E {[X_{t_{0} + ρ u} - X_{t_{0} + ρ v}]}^{2} = + \infty \end{matrix}

Therefore

E {[Y_{u}^{α} (ρ) - Y_{v}^{α} (ρ)]}^{2}

admits a limit for all

u, v \in R_{+}^{N}

when

ρ \to 0

if and only if

\begin{matrix} {\begin{matrix} α = {inf}_{u, v} β_{u v} (t_{0}) \\ and \\ {lim}_{ρ \to 0} \frac{| H (t_{0} + ρ u) - H (t_{0} + ρ v) |}{ρ^{{inf}_{u, v} β_{u v} (t_{0})}} = Γ (u, v) \in R_{+}^{*} \end{matrix} \end{matrix}

In that case, we have for all

u

v

R_{+}^{N}

\begin{matrix} E {[Y_{u}^{α} (ρ) - Y_{v}^{α} (ρ)]}^{2} \overset{ρ \to 0}{⟶} \frac{\partial^{2} φ}{\partial x^{2}} (2 H (t_{0}); ∥ t_{0} ∥) {[Γ (u, v)]}^{2} \end{matrix}

Remark 6 We can see easily that

\begin{matrix} β_{\frac{u}{∥ u ∥}} (t_{0}) \land β_{\frac{v}{∥ v ∥}} (t_{0}) \leq β_{u v} (t_{0}) \end{matrix}

(50)

hence

\begin{matrix} {inf}_{u \in U} β_{u} (t_{0}) \leq {inf}_{u, v} β_{u v} (t_{0}) \end{matrix}

(51)

Conversely, assume there exist

u, v \in U

such that

β_{u} (t_{0}) < β_{v} (t_{0})

, and let

α \in (β_{u} (t_{0}); β_{v} (t_{0}))

. The inequality

\frac{| H (t_{0} + ρ u) - H (t_{0}) |}{ρ^{α}} \leq \frac{| H (t_{0} + ρ u) - H (t_{0} + ρ v) |}{ρ^{α}} + \frac{| H (t_{0} + ρ v) - H (t_{0}) |}{ρ^{α}}

implies

{limsup}_{ρ \to 0} \frac{| H (t_{0} + ρ u) - H (t_{0} + ρ v) |}{ρ^{α}} = + \infty

and therefore

α > β_{u v} (t_{0})

. Then

{inf}_{u, v} β_{u v} (t_{0}) \leq {inf}_{u \in U} β_{u} (t_{0})

, which gives

\begin{matrix} {inf}_{u, v} β_{u v} (t_{0}) = {inf}_{u \in U} β_{u} (t_{0}) \end{matrix}

(52)

5.1.2 Multifractional Brownian sheet

In the non-isotropic case, using remark 3 , consider

N

independent mBm

X^{(i)}

with parameter function

H_{i}

\begin{matrix} E [Δ X_{t_{0}, t_{0} + ρ u} Δ X_{t_{0}, t_{0} + ρ v}] =^{N} \prod_{i = 1} E [(X_{t_{0}^{(i)} + ρ u^{(i)}}^{(i)} - X_{t_{0}^{(i)}}^{(i)}) (X_{t_{0}^{(i)} + ρ v^{(i)}}^{(i)} - X_{t_{0}^{(i)}}^{(i)})] \end{matrix}

As in the isotropic case, for all

i \in {1, . . ., N}

, consider

β_{u v}^{i} (t_{0}) = sup {α; {lim}_{ρ \to 0} \frac{| H_{i} (t_{0} + ρ u) - H_{i} (t_{0} + ρ v) |}{ρ^{α}} = 0}

Each process

X^{(i)}

is locally asymptoticly self-similar, therefore

E [\frac{X_{t_{0}^{(i)} + ρ u^{(i)}}^{(i)} - X_{t_{0}^{(i)}}^{(i)}}{ρ^{α_{i}}} \times \frac{X_{t_{0}^{(i)} + ρ v^{(i)}}^{(i)} - X_{t_{0}^{(i)}}^{(i)}}{ρ^{α_{i}}}] \overset{ρ \to 0}{⟶} E [Y_{u}^{α_{i}} Y_{v}^{α_{i}}]

where

Y^{α_{i}}

denotes

∙ fractional Brownian motion of parameter $α_{i} = H_{i} (t_{0})$ , in the case $H_{i} (t_{0}) < {inf}_{u, v} β_{u v}^{i} (t_{0})$ ,
∙ the centered Gaussian process such that $E {[Y_{u}^{α_{i}} - Y_{v}^{α_{i}}]}^{2} = K_{t_{0}} {[Γ_{i} (u, v)]}^{2}$ where $α_{i} = {inf}_{u, v} β_{u v}^{i} (t_{0})$ , in the case $H_{i} (t_{0}) > {inf}_{u, v} β_{u v}^{i} (t_{0})$ and ${lim}_{ρ \to 0} \frac{| H_{i} (t_{0} + ρ u) - H_{i} (t_{0} + ρ v) |}{ρ^{{inf}_{u, v} β_{u v}^{i} (t_{0})}} = Γ_{i} (u, v)$ with $Γ_{i}$ bounded on $[a, b]^{2}$ .

Then we conclude

E [\frac{Δ X_{t_{0}^{(i)}, t_{0}^{(i)} + ρ u^{(i)}}^{(i)}}{ρ^{\sum_{i} α_{i}}} \times \frac{Δ X_{t_{0}^{(i)}, t_{0}^{(i)} + ρ v^{(i)}}^{(i)}}{ρ^{\sum_{i} α_{i}}}] \overset{ρ \to 0}{⟶} E [Y_{u}^{α} Y_{v}^{α}]

where

Y^{α} = \prod_{i = 1}^{N} Y^{α_{i}}

5.2 Tightness of laws

The study of weak convergence is well-known for stochastic processes indexed by

R_{+}

. A comprehensive review was made by Billingsley (cf [5] ) for a compact set of index (

[0, 1]

). In ([11] ), Karatzas and Shreeve stated the same kind of results for the whole

R_{+}

. The case of

R_{+}^{N}

can be found in ([10] ) whose corollary 14.9 provides

Proposition 18 Consider a sequence of continuous processes

{(X^{(n)})}_{n \in N}

with

X^{(n)} = {X_{t}^{(n)}; t \in R_{+}^{N}}

(Ω, ℱ, P)

such that

1. there exists a positive constant $ν$ such that ${sup}_{n \geq 1} E {| X_{0}^{(n)} |}^{ν} < \infty$
2. for all $T > 0$ and all $s$ , $t$ in $[0, T]^{N}$ , there exist positive constants $α$ , $β$ and $C_{T}$ such that ${sup}_{n \geq 1} E {| X_{t}^{(n)} - X_{s}^{(n)} |}^{α} \leq C_{T} ∥ t - s ∥^{N + β}$

Then the probability measures

P_{n} \overset{Δ}{=} P . {(X^{(n)})}^{- 1}

(C (R_{+}^{N}), ℬ (C (R_{+}^{N})))

form a tight sequence.

We verify the conditions of proposition 18 , in the case of mBm, in the following sections.

5.2.1 Multifractional Brownian field

By ( 18 ), there exist positive constants

K_{T}

and

L_{T}

such that for all

u

v

[0, T]^{N}

\begin{matrix} ρ^{2 α} E {[Y_{u}^{α} (ρ) - Y_{v}^{α} (ρ)]}^{2} & = & E {[X_{t_{0} + ρ u} - X_{t_{0} + ρ v}]}^{2} \end{matrix}

\begin{matrix} \leq & K_{T} ∥ ρ . (u - v) ∥^{2 H (t_{0} + ρ u)} \end{matrix}

\begin{matrix} + L_{T} {| H (t_{0} + ρ u) - H (t_{0} + ρ v) |}^{2} \end{matrix}

Therefore,

\begin{matrix} E {[Y_{u}^{α} (ρ) - Y_{v}^{α} (ρ)]}^{2} \leq K_{T}^{'} ρ^{2 (H (t_{0}) - α)} . ∥ (u - v) ∥^{2 H (t_{0})} + L_{T} \frac{{| H (t_{0} + ρ u) - H (t_{0} + ρ v) |}^{2}}{ρ^{2 α}} \end{matrix}

∙ In the case

H (t_{0}) < {inf}_{u, v} β_{u v} (t_{0})

, there exists

M_{T} > 0

such that

\begin{matrix} E {[Y_{u}^{H (t_{0})} (ρ) - Y_{v}^{H (t_{0})} (ρ)]}^{2} \leq M_{T} ∥ u - v ∥^{2 H (t_{0})} \end{matrix}

∙ In the case

H (t_{0}) > {inf}_{u, v} β_{u v} (t_{0})

, under the assumption

{lim}_{ρ \to 0} \frac{| H (t_{0} + ρ u) - H (t_{0} + ρ v) |}{ρ^{{inf}_{u, v} β_{u v} (t_{0})}} = Γ (u, v)

with

(u, v) \mapsto \frac{Γ (u, v)}{∥ u - v ∥^{2 β}}

bounded on

[a, b]^{2}

, there exists

M_{T} > 0

such that

\begin{matrix} E {[Y_{u}^{{inf}_{u, v} β_{u v} (t_{0})} (ρ) - Y_{v}^{{inf}_{u, v} β_{u v} (t_{0})} (ρ)]}^{2} \leq M_{T} ∥ u - v ∥^{2 (β \land H (t_{0}))} \end{matrix}

Since the process

Y^{α}

is Gaussian, we get an exponent greater than

N

in the usual way. Then we can conclude by proposition 18 that the laws of

Y^{α}

are tight.

5.2.2 Multifractional Brownian sheet

In the same way as in paragraph 4.1.2

\begin{matrix} E {[Y_{u}^{(ρ)} - Y_{v}^{(ρ)}]}^{2} & = & \frac{1}{ρ^{2 \sum_{i} α_{i}}} E [^{N} \prod_{i = 1} (X_{t_{0}^{(i)} + ρ u^{(i)}}^{(i)} - X_{t_{0}^{(i)}}^{(i)}) -^{N} \prod_{i = 1} (X_{t_{0}^{(i)} + ρ v^{(i)}}^{(i)} - X_{t_{0}^{(i)}}^{(i)})] \end{matrix}

\begin{matrix} \leq & K \sum_{i} E {[\frac{X_{t_{0}^{(i)} + ρ u^{(i)}}^{(i)} - X_{t_{0}^{(i)} + ρ v^{(i)}}^{(i)}}{ρ^{α_{i}}}]}^{2} \end{matrix}

then, under the assumptions of theorem 4 , there exists a positive constant

M_{T}

such that

\begin{matrix} E {[Y_{u}^{(ρ)} - Y_{v}^{(ρ)}]}^{2} \leq M_{T} ∥ u - v ∥^{2 {min}_{i} α_{i}} \end{matrix}

We conclude as in the isotropic case.

Acknowledgement The author thanks Jacques Lévy-Véhel for all their fruitful discussions, especially about the Hölder regularity.

References

A. Ayache and J. Lévy Véhel, Generalized multifractional Brownian motion: definition and preliminary results, in: Fractals: theory and application in engineering, M. Dekking, J. Levy-Véhel, E. Lutton and C. Tricot, Springer, 1999.
A. Ayache and J. Lévy Véhel, Generalized multifractional Brownian motion, SISP, 3, 1/2, 7-18, 2000.
A. Ayache, S. Cohen and J. Lévy Véhel, The covariance structure of multifractional Brownian motion, with application to long range dependence. ICASSP, 2000.
A. Benassi, S. Jaffard and D. Roux, Elliptic Gaussian random processes, Rev. Mat. Ibe., vol. 13, 19-89, 1998.
P. Billingsley, Convergence of probability measures, Wiley series in probability and statistics, 2nd edition, 1999.
S. Cohen, From self-similarity to local self-similarity : the estimation problem, in: Fractals : theory and application in engineering, M. Dekking, J. Levy-Véhel, E. Lutton and C. Tricot, Springer, 1999.
E. Herbin and J. Lévy Véhel, Fine analysis of the regularity of Gaussian processes: Stochastic $2$ -microlocal analysis, preprint, 2004.
S. Leger and M. Pontier, Drap brownien fractionnaire, Note aux CRAS, Paris, t.329, série I, mathématiques, 893-898, 1999.
T. Lindstrom, Fractional Brownian fields as integrals of white noise, Bull. London Math. Soc., 25, 83-88, 1993.
O. Kallenberg, Foundations of modern probability, Springer, 1997.
I. Karatzas, S. Shreve, Brownian motion and stochastic calculus, Springer, 1991.
D. Khoshnevisan, Multiparameter processes, an introduction to random fields, Springer, 2002.
R. Peltier and J. Lévy-Véhel, Multifractional Brownian motion: definition and preliminary results, Rapport de recherche INRIA 2645, 1995.
B. Pesquet-Popescu, Modélisation bidimensionnelle de processus non stationnaires et application à l'étude du fond sous-marin, thèse de l'ENS Cachan, 1998.

From $N$ parameter fractional Brownian motions to $N$ parameter multifractional Brownian motions

Erick Herbin INRIA, Domaine de Voluceau, Rocquencourt, BP 105, 78153 Le Chesnay Cedex, France erick.herbin@inria.fr and Dassault Aviation, 78 quai Marcel Dassault, 92552 Saint-Cloud Cedex, France erick.herbin@dassault-aviation.fr

November 2002