Generalized k-variations and Hurst parameter estimation for the fractional wave equation via Malliavin calculus

1 Introduction

Since several decades, the statistical inference in stochastic (partial) differential equations (S(P)DE in the sequel) constitutes an intensive research direction in probability theory and mathematical statistics. Nowadays, a particular case of wide interest is represented by the S(P)DE driven by fractional Brownian motion (fBm) and related processes, due to the vast area of application of such stochastic models. Many recent works concern the estimation of the drift parameter for stochastic equation driven by fractional Brownian motion (we refer, among many others, to

[1], [9], [12], [24]), while fewer works deal with the estimation of the Hurst parameter in such stochastic equations.

In this paper, we consider the stochastic wave equation driven by an additive Gaussian noise which behaves as a fractional Brownian motion in time and as a Wiener process in space (we call it fractional-white noise). Our purpose is to construct and analyze an estimator for the Hurst parameter of the solution to this SPDE based on the observation of the solution at a fixed time and at a discrete number of points in space. The wave equation with fractional noise in time and/or in space has been studied in several works, such as [3], [7], [8], [11], [20] etc. We will use a standard method to construct estimators for the Hurst parameter, which is based on the $k$ -variations of the observed process. The method has been recently employed in [11] for the case of quadratic variations, i.e. $k = 2$ . As for the fBm, it was shown that the standard quadratic variation estimator is not asymptotically normal when the Hurst index becomes bigger than $\frac{3}{4}$ and this is inconvenient for statistical applications. In order to avoid this restriction and to get an estimator which is asymptotically Gaussian for every $H$ , we will use the generalized $k$ -variations, which basically means that the usual increment of the process is replaced by a higher order increment. The idea comes from the reference [10] and since it has been used by many authors (see e.g. [6] or [5]). More precisely, if $(u (t, x), t \geq 0, x \in R)$

denotes the solution to the wave equation with fractional-white noise, we define the (centered) generalized

k

-variation statistics (

k \geq 1

integer),

V_{N} (k, α) = \frac{1}{N - l} N \sum i = l ⎡ ⎢ ⎢ ⎢ ⎣ \frac{{∣ ∣ U^{α} (\frac{i}{N}) ∣ ∣}^{k}}{E {∣ ∣ U^{α} (\frac{i}{N}) ∣ ∣}^{k}} - 1 ⎤ ⎥ ⎥ ⎥ ⎦,

(1)

where $U^{α} (\frac{i}{N})$ represents the spatial increment of the solution $u$ at $\frac{i}{N}$ along a filter $α$ of power (order) $p \geq 1$ and length $l + 1 \geq 1$ (see the next section for the precise definition).

By using chaos expansion and the recent Stein-Malliavin calculus we show that the sequence (1) satisfies a Central Limit Theorem (CLT) as $N \to \infty$ (in the spirit of [4]) whenever $p > H + \frac{1}{4}$ and in this way the restriction $H < \frac{3}{4}$ can be avoided by choosing a filter of order $p \geq 2$ , i. e. by replacing, for example, the usual increment by a higher order increment. We will obtain the rate of convergence under the Wasserstein distance for this convergence in law and we also prove a multidimensional CLT. So we generalize the findings in [11] to filters of any power $p \geq 1$ and to $k$ -variations of any order $k \geq 1$ and in addition we show that in the special case $p = 1$ and $H > \frac{3}{4}$ a non-Gaussian limit theorem occurs with limit distribution related to the Rosenblatt distribution (but more complex than it).

These theoretical results are then applied to the estimation of Hurst index of the solution the the fractional-white wave equation. Based on the behavior of the sequence (1), we prove that the associated $k$ -variation estimator for $H$ is consistent and asymptotically normal. Moreover, we provide a numerical analysis of the estimator when $k = 2$ by analysing its performance on various filters and for several values of the Hurst parameter and confirming via simulation the theoretical results.

We organized the paper as follows. Section 2 contains some preliminaries. We present in this part the basic facts concerning the solution to the fractional-white wave equation, we introduce the filters and the increment of the solution along filters. In Section 3, we prove a CLT for the sequence (1) for any integer $k \geq 1$ and we obtain the rate of convergence when $k$ is even via the Stein-Malliavin theory. In Section 4, we show a non-central limit limit in the case $k = 2, H > \frac{3}{4}$ and for filters of order $p = 1$ . Section 5 concerns the estimation of the Hurst parameter of the solution to the fractional-white wave equation (2.1). We included here theoretical results related to the behavior of the $k$ -variations estimators for the Hurst index as well as simulations and numerical analysis for the performance of the estimators. Section 6 (the Appendix) contains the basic tools from Malliavin calculus needed in the paper.

2 Preliminaries

We introduce here the fractional-white heat equation and its solution and we present the basic definitions and the notation concerning the filters used in our work.

2.1 The solution to the wave equation with fractional-colored noise

The object of our study will be the solution to the following stochastic wave equation

⎧⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎨⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎩∂2u∂t2(t,x)=Δu(t,x)+˙WH(t,x),t≥0,x∈Rd,d≥1,\omit\rm span\omit\rm span\@@LTX@noalign\vskip5.690551pt\omitu(0,x)=0,x∈Rd,\omit\rm span\omit\rm span\@@LTX@noalign\vskip5.690551pt\omit∂u∂t(0,x)=0,x∈Rd,

where $Δ$ is the Laplacian on $R^{d}$ , $d ⩾ 1$ , and $W^{H}$ is a fractional-white Gaussian noise which is defined as a real valued centered Gaussian field $W_{H} = {W_{t}^{H} (A); t \in [0, T], A \in B_{b} (R^{d})}$ , over a given complete filtered probability space $(Ω, F, (F_{t})_{t ⩾ 0}, P)$ , with covariance function given by

E (W_{t}^{H} (A) W_{s}^{H} (B)) = R_{H} (t, s) λ (A \cap B), \forall A, B \in B_{b} (R^{d}),

(2)

where $R_{H}$ is the covariance of the fractional brownian motion

R_{H} (t, s) = \frac{1}{2} (t^{2 H} + s^{2 H} - | t - s |^{2 H}), s, t \geq 0.

We denoted by $B_{b} (R^{d})$ the class of bounded Borel subsets of $R^{d}$ and we will assume throughout this work $H \in (\frac{1}{2}, 1) .$

The solution of the equation (2.1) is understood in the mild sense, that is, it is defined as a square-integrable centered field $u = (u (t, x); t \in [0, T], x \in R^{d})$ defined by

u (t, x) = \int_{0}^{t} \int_{R^{d}} G_{1} (t - s, x - y) W^{H} (d s, d y), t \geq 0, x \in R^{d},

(3)

where $G_{1}$ is the fundamental solution to the wave equation and the integral in (3) is a Wiener integral with respect to the Gaussian process $W^{H}$ . Recall that for $d = 1$ (we will later restrict to this situation in our work) we have, for $t \geq 0$ and $x \in R$ ,

G_{1} (t, x) = \frac{1}{2} 1_{{| x | < t}} .

(4)

We know (see e.g. [3]) that the solution (3) is well-defined if and only if

d < 2 H + 1

and it is self-similar in time and stationary in space. Other properties of the solution can be found in [3], [7] or [23]. In particular, the spatial covariance of the solution can be expressed as follows

	$E (u (t, x) u (t, y))$	$=$	$\frac{1}{2} (c_{H} \| y - x \|^{2 H + 1} - t \frac{\| y - x \|^{2 H}}{2} + \frac{t^{2 H + 1}}{2 H + 1}) 1_{{\| y - x \| < t}}$		(5)
			$+ \frac{(2 t - \| y - x \|)^{2 H + 1}}{8 (2 H + 1)} 1_{{t \leq \| y - x \| < 2 t}}$		(5)

with $c_{H} = \frac{4 H - 1}{4 (2 H + 1)}$ . When $t > 1$ and $x, y \in [0, 1]$ , this expression reduces to

E (u (t, x) u (t, y))

=

\frac{1}{2} (c_{H} | y - x |^{2 H + 1} - t \frac{| y - x |^{2 H}}{2} + \frac{t^{2 H + 1}}{2 H + 1}) for x, y \in [0, 1] .

(6)

We will fix for the rest of the work $t > 1$ and we will associate to the process $(u (t, x), x \in [0, 1])$ its canonical Hilbert space $H$ which is defined as the closure of the linear space generated by the indicator functions ${1_{[0, x]}, x \in [0, 1]}$ with respect to the inner product

⟨ 1_{[0, x]}, 1_{[0, y]} ⟩_{H}

=

E (u (t, x) u (t, y)) .

We will denote by $I_{q}$ the multiple stochastic integral of order $q \geq 1$ with respect to the Gaussian process $(u (t, x), x \in [0, 1])$ and by $D$ the Malliavin derivative with respect to this process. We refer to the Appendix for the basic elements of the Malliavin calculus.

We will also use in Section 4 multiple stochastic integrals with respect to the fractional-white noise $W^{H}$ with covariance (2). We use the notation $I_{q}^{W}$ to indicate the multiple integral of order $q \geq 1$ with respect to $W^{H}$ .

2.2 Filters

In this paragraph we will define the filters and the increments of the solution to (2.1) along filters. We start with several definitions and notations needed along this paper.

Definition 1

Given $l \in N$ and $p \in N^{*}$

, a vector

α = (α_{0}, . . ., α_{l})

is called a filter of length

l + 1 \geq 1

and order (or power)

p \geq 1

such that

⎧ ⎨ ⎩ \begin{matrix} \sum_{q = 0}^{l} α_{q} q^{r} & = 0, & 0 \leq r \leq p - 1, \sum_{q = 0}^{l} α_{q} q^{p} & \neq 0 \end{matrix}

with the convention $0^{0} = 1$ .

For a filter $α = (a_{0}, a_{1}, . ., a_{l})$ of length $l + 1 \geq 1$ and of order $p \geq 1$ we define the space-filtered process (or the spatial increment of the process $u$ along the filter $α$ )

U^{α} (\frac{i}{N}) = l \sum r = 0 a_{r} u (t, \frac{i - r}{N}) for i = l, . ., N .

(7)

We denote for $j \geq 1$

π_{H}^{α, N} (j) := E [U^{α} (\frac{i}{N}) U^{α} (\frac{i + j}{N})] .

From the covariance formula (5) we can write

	$π_{H}^{α, N} (j)$	$=$	$l \sum r_{1}, r_{2} = 0 a_{r_{1}} a_{r_{2}} E [u (t, \frac{i - r_{1}}{N}) u (t, \frac{i + j - r_{2}}{N})]$		(8)
		$=$	$k_{1} \frac{1}{N^{2 H}} Φ_{H, α} (j) + k_{2} \frac{1}{N^{2 H + 1}} Φ_{H + \frac{1}{2}, α} (j)$		(8)

with

Φ_{H, α} (j) = l \sum r_{1}, r_{2} = 0 a_{r_{1}} a_{r_{2}} | j + r_{1} - r_{2} |^{2 H}, j \geq 0

and $k_{1} = - \frac{t}{4}$ and $k_{2} = \frac{c_{H}}{2} = \frac{4 H - 1}{8 (2 H + 1)}$ . We write for further use

c_{1} (H) = \frac{- t}{4} l \sum q, r = 0 α_{q} α_{r} | q - r |^{2 H}, and c_{2} (H) = \frac{c_{H}}{2} l \sum q, r = 0 α_{q} α_{r} | q - r |^{2 H + 1} .

(9)

In particular, from (8)

	$π_{H}^{α, N} (0)$	$=$	$E {[U^{α} (\frac{i}{N})]}^{2} = k_{1} \frac{1}{N^{2 H}} Φ_{H, α} (0) + k_{2} \frac{1}{N^{2 H + 1}} Φ_{H + \frac{1}{2}, α} (0)$
		$=$	$c_{1} (H) \frac{1}{N^{2 H}} + c_{2} (H) \frac{1}{N^{2 H + 1}}$

We will need the below technical lemma to etablish the asymptotical equivalent of $Φ_{H, α}$ and similar expressions. The proof of the lemma is based on a Taylor expansion, see [6] or [10].

Lemma 1

Let $H \in R^{+} ∖ N$ and $α^{(1)}$ , $α^{(2)}$ be filters of lengths $l_{1} + 1$ , $l_{2} + 2$ and of orders $p_{1}, p_{2} \geq 1$ respectively. Then

l_{1} \sum q = 0 l_{2} \sum r = 0 α_{q}^{(1)} α_{r}^{(2)} | q - r + k |^{2 H} \sim k \to \infty

κ_{H} k^{2 H - 2 p}

with $κ_{H} = \sum_{q = 0}^{l_{1}} \sum_{r = 0}^{l_{2}} α_{q}^{(1)} α_{r}^{(2)} \frac{2 H (2 H - 1) \dots (2 H - 2 p + 1)}{2 p!} (q - r)^{2 p}$ , where $p = min (p_{1}, p_{2})$ .

In the sequel, we write $a_{k} \sim_{k \to \infty} b_{k}$ to indicate that the sequences $a_{k}, b_{k}$ have the same behavior as $k \to \infty$ .

3 Central limit theorem for the spatial $k$ -variations

In this section we focus on the asymptotic behavior in distribution of the $k$ -variation in space of the solution to the fractional-white wave equation, defined via a filter of power $p \geq 1$ . In the first step we show the $k$ -variation satisfies a CLT when $p > H + \frac{1}{4}$ . Next, by taking $k$ to be an even integer, we derive a Berry-Esséen type bound for this convergence in distribution via the Stein-Malliavin calculus. Restricting ourselves in addition to $k = 2$ , we prove a multidimensional CLT, which is needed for the estimation of the Hurst parameter.

3.1 Central Limit Theorem

Fix $t > 1$ and let $α$ be a filter of length $l + 1 \geq 1$ and of power $p \geq 1$ as in Definition 1. Let $u$ be given by (3). For any integer $k \geq 1$ we define the centered spatial $k$ -variations of the process $(u (t, x), x \in R)$ by

V_{N} (k, α) = \frac{1}{N - l} N \sum i = l ⎡ ⎢ ⎢ ⎢ ⎣ \frac{{∣ ∣ U^{α} (\frac{i}{N}) ∣ ∣}^{k}}{E {∣ ∣ U^{α} (\frac{i}{N}) ∣ ∣}^{k}} - 1 ⎤ ⎥ ⎥ ⎥ ⎦

(10)

with $U^{α} (\frac{i}{N})$ given by (7). We will show that the sequence (10) satisfies a CLT. In order to do this we will use a criterion based on Malliavin calculus. The first step is to expand in chaos the $k$ -variation sequence $V_{N} (k, α)$ . Noticing that the filtered process $U^{α}$

as a linear combination of centered Gaussian random variables is a centered Gaussian process, we get

E (U^{α} {(\frac{i}{N})}^{k}) = E_{k} E {(U^{α} {(\frac{i}{N})}^{2})}^{\frac{k}{2}},

(11)

where $E_{k}$ denotes the $k$

-th absolute moment of a standard Gaussian variable given by

E_{k} = \frac{2^{\frac{k}{2}} Γ (\frac{k + 1}{2})}{Γ (\frac{1}{2})}

. We introduce the variable

Z^{α} (\frac{i}{N}) = \frac{U^{α} (\frac{i}{N})}{(π_{H}^{α, N} (0))^{1 / 2}} .

(12)

It is clear that $Z^{α} (\frac{i}{N})$ is a standard Gaussian variable and $Corr (Z^{α} (\frac{i}{N}), Z^{α} (\frac{j}{N})) = Corr (U^{α} (\frac{i}{N}), U^{α} (\frac{j}{N}))$ , where $Corr$ denotes the correlation coefficient. Using (11) and (12) we can write $V_{N}$ as follows:

V_{N} (k, α)

=

\frac{1}{N - l} N \sum i = l ⎡ ⎣ \frac{| U^{α} (\frac{i}{N}) |^{k}}{E | U^{α} (\frac{i}{N}) |^{k}} - 1 ⎤ ⎦ = \frac{1}{N - l} N \sum i = l ⎡ ⎣ \frac{| Z^{α} (\frac{i}{N}) |^{k}}{E_{k}} - 1 ⎤ ⎦ .

Recall the expansion of the development in Hermite polynomials of the function $H^{k} (t) = \frac{| t |^{k}}{E_{k}} - 1$ given in Lemma 2 of [6],

H^{k} (t) = \infty \sum j = 1 c_{j}^{k} H_{j} (t),

where $c_{2 j + 1}^{k} = 0$ for $j ⩾ 0$ , $c_{2 j}^{k} = \frac{1}{(2 j)!} \prod_{i = 0}^{j - 1} (k - 2 i)$ for $j ⩾ 1$ and $H_{j} (t)$ denotes the j-th Hermite polynomial defined by

H_{j} (t) = [\frac{j}{2}] \sum a = 0 (- 1)^{a} \frac{a!}{(j - 2 a)! a!} 2^{- a} t^{j - 2 a} .

Observing that for

C_{i, α} := l \sum q = 0 α_{q} 1_{[0, \frac{i - q}{N}]}

we have from (8) that ${∥ ∥ ∥ \frac{C_{i, α}}{(π_{H}^{α, N} (0))^{1 / 2}} ∥ ∥ ∥}_{H} = 1$ we can express $Z^{α} (\frac{i}{N})$ as an integral with respect to the process $(u (t, x), x \in R)$ since the increment $u (t, y) - u (t, x)$ can be expressed as $I_{1} (1_{[x, y]})$ (recall that $I_{1}$ represents the multiple integral of order 1 with respect to the Gaussian process $(u (t, x), x \in [0, 1])$ ) for every $x < y$ :

Z^{α} (\frac{i}{N})

=

I_{1} (\frac{C_{i, α}}{(π_{H}^{α, N} (0))^{1 / 2}}) .

Since we have $H_{q} (I_{1} (h)) = \frac{1}{q!} I_{q} (h^{\otimes q})$ for $∥ h ∥_{H} = 1$ we get

$V_{N} (k, α)$	$=$	$\frac{1}{N - l} N \sum i = l H^{k} (Z^{α} (\frac{i}{N})) = \frac{1}{N - l} \sum q ⩾ 1 c_{2 q}^{k} N \sum i = l H_{2 q} (Z^{α} (\frac{i}{N}))$
	$=$	$\frac{1}{N - l} \sum q ⩾ 1 c_{2 q}^{k} N \sum i = l H_{2 q} (I_{1} (\frac{C_{i, α}}{(π_{H}^{α, N} (0))^{1 / 2}}))$
	$=$	$\frac{1}{N - l} \sum q ⩾ 1 \frac{c_{2 q}^{k}}{(2 q)!} N \sum i = l I_{2 q} ⎛ ⎝ {(\frac{C_{i, α}}{(π_{H}^{α, N} (0))^{1 / 2}})}^{\otimes 2 q} ⎞ ⎠ .$

Hence, we obtain the following chaotic expansion of the $k$ -variation statistics

V_{N} (k, α) = \frac{1}{N - l} N \sum i = l \infty \sum q = 1 \frac{c_{2 q}^{k}}{(2 q)!} I_{2 q} ⎛ ⎝ \frac{C_{i, α}^{\otimes 2 q}}{(π_{H}^{α, N} (0))^{q}} ⎞ ⎠ = \sum q \geq 1 I_{2 q} (f_{N, 2 q})

(13)

with

f_{N, 2 q} = \frac{c_{2 q}^{k}}{(2 q)!} \frac{1}{N - l} N \sum i = l \frac{C_{i, α}^{\otimes 2 q}}{(π_{H}^{α, N} (0))^{q}} .

(14)

Let us start by analyzing the asymptotic behavior of the mean square of each kernel $f_{N, 2 q}$ .

Lemma 2

For $N, q \geq 1$ , let $f_{N, 2 q}$ be given by (14). Then

(N - l) (2 q)! ∥ f_{N, 2 q} ∥_{H^{\otimes 2 q}}^{2} \to_{N \to \infty} \frac{(c_{2 q}^{k})^{2}}{(2 q)!} \sum v \in Z (φ_{H, α} (v))^{2 q} := σ_{2 q}^{2}

for $H < p - \frac{1}{4 q}$ (i.e. $H < 1 - \frac{1}{4 q}$ for $p = 1$ and $H \in (\frac{1}{2}, 1)$ for $p \geq 2$ ), where we use the notation

φ_{H, α} (v) = \frac{Φ_{H, α} (v)}{Φ_{H, α} (0)} .

(15)

Moreover, $σ^{2} := \sum_{q \geq 1} σ_{2 q}^{2} < \infty .$ For $p = q = 1$ , $H = 3 / 4$

\frac{N - l}{log (N - l)} 2! ∥ f_{N, 2} ∥_{H^{\otimes 2}}^{2} \to_{N \to \infty} c^{2} := \frac{(c_{2}^{k})^{2}}{2} lim N log N \sum | v | \leq N (ρ_{H, α} (v))^{2} < \infty .

(16)

Proof: From (14), we get

	$(2 q)! ∥ f_{N, 2 q} ∥_{H^{\otimes 2 q}}^{2}$	$=$	$\frac{(c_{2 q}^{k})^{2}}{(2 q)!} \frac{1}{(N - l)^{2}} N \sum i, j = l \frac{⟨ C_{i, α}, C_{j, α} ⟩_{H}^{2 q}}{(π_{H}^{α, N} (0))^{2 q}}$
		$=$	$\frac{1}{(N - l)^{2}} \frac{(c_{2 q}^{k})^{2}}{(2 q)!} N \sum i, j = l {(ρ_{H}^{α, N} (j - i))}^{2 q},$

where we used the notation

ρ_{H}^{α, N} (v) = \frac{π_{H}^{α, N} (v)}{π_{H}^{α, N} (0)} for v \in Z .

(17)

Next, we write

\frac{1}{N - l} N \sum i, j = l {(ρ_{H}^{α, N} (j - i))}^{2 q} = \sum v \in Z {(ρ_{H}^{α, N} (v))}^{2 q} 1_{{| v | \leq N - l}} \frac{N - | v | - l}{N - l}

and thus

(N - l) (2 q)! ∥ f_{N, 2 q} ∥_{H^{\otimes 2 q}}^{2}

=

\frac{(c_{2 q}^{k})^{2}}{(2 q)!} \sum v \in Z {(ρ_{H}^{α, N} (v))}^{2 q} 1_{{| v | \leq N - l}} \frac{N - | v | - l}{N - l} .

(18)

Using the expression

ρ_{H}^{α, N} (v) = \frac{k_{1} Φ_{H, α} (v) N^{- 2 H} + k_{2} Φ_{H + \frac{1}{2}, α} (v) N^{- 2 H - 1}}{k_{1} Φ_{H, α} (0) N^{- 2 H} + k_{2} Φ_{H + \frac{1}{2}, α} (0) N^{- 2 H - 1}} = \frac{Φ_{H, α} (v) + a_{N} (v)}{Φ_{H, α} (0) + a_{N} (0)}

with

a_{N} (v) = \frac{k_{2}}{k_{1} N} Φ_{H + \frac{1}{2}, α} (v)

(19)

we can write, with $φ_{H, α}$ and $ρ_{H}^{α, N}$ given by (15) and (17) respectively,

b_{N, H} (v) := ρ_{H}^{α, N} (v) - φ_{H, α} (v)

(20)

and remark that due to Lemma 1 for $v$ large enough

| b_{N, H} (v) | \sim | v | \to \infty ∣ ∣ ∣ a_{N} (v) \frac{1}{Φ_{H, α} (0) + a_{N} (0)} ∣ ∣ ∣ \leq C \frac{1}{N} v^{2 H + 1 - 2 p},

(21)

where $C > 0$ does not depend on $N, v$ .With this notation we can write

$(N - l) (2 q)! ∥ f_{N, 2 q} ∥_{H^{\otimes 2 q}}^{2}$	$=$	$\frac{(c_{2 q}^{k})^{2}}{(2 q)!} \sum v \in Z {(φ_{H, α} (v) + b_{N, H} (v))}_{H, α}^{2 q} 1_{{\| v \| \leq N - l}} \frac{N - \| v \| - l}{N - l}$
	$=$	$\frac{(c_{2 q}^{k})^{2}}{(2 q)!} \sum v \in Z 2 q \sum m = 0 (\frac{2 q}{m}) φ_{H, α} (v)^{m} (b_{N, H} (v))^{2 q - m} 1_{{\| v \| \leq N - l}} \frac{N - \| v \| - l}{N - l}$
	$=$	$\frac{(c_{2 q}^{k})^{2}}{(2 q)!} \sum v \in Z φ_{H, α} (v)^{2 q} 1_{{\| v \| \leq N - l}} \frac{N - \| v \| - l}{N - l} + r_{N, q, 1}$

with

r_{N, q, 1} = \frac{(c_{2 q}^{k})^{2}}{(2 q)!} \sum v \in Z 2 q - 1 \sum m = 0 (\frac{2 q}{m}) φ_{H, α} (v)^{m} (b_{N, H} (v))^{2 q - m} 1_{{| v | \leq N - l}} \frac{N - | v | - l}{N - l} .

(22)

Clearly, by the dominated convergence theorem

\frac{(c_{2 q}^{k})^{2}}{(2 q)!} \sum v \in Z φ_{H, α} (v)^{2 q} 1_{{| v | \leq N - l}} \frac{N - | v | - l}{N - l} \to_{N \to \infty} σ_{2 q}^{2},

which by Lemma 1 is finite if $p = 1$ , $H < 1 - \frac{1}{4 q}$ , and for all $H \in (1 / 2, 1)$ in the other cases. For $q = p = 1$ , $H = 3 / 4$ ,

\frac{1}{log (N - l)} \sum v \in Z φ_{H, α} (v)^{2} 1_{{| v | \leq N - l}} \frac{N - | v | - l}{N - l}

converges to a positive constant and thus (16) is obtained.

In order to conclude it remains to show that the rest term $r_{N, q, 1}$ (22) converges to $0$ as $N \to \infty$ , for every $q \geq 1$ . From (22), using the bound (21) and Lemma 1, we have the estimate

| r_{N, q, 1} |

\leq

C 2 q - 1 \sum m = 0 (\frac{2 q}{m}) \frac{1}{N^{2 q - m}} \sum 1 \leq v \leq N - l | v |^{(2 H - 2) m} | v |^{(2 H + 1 - 2 p) (2 q - m)} := 2 q - 1 \sum m = 0 r_{N, q, 1, m}

and for each $m = 0, . ., 2 q - 1$ ,

r_{N, q, 1, m} \leq \frac{C}{N^{2 q - m}} \sum 1 \leq v \leq N - l | v |^{(2 H - 2 p) 2 q + 2 q - m} .

When the series $\sum_{v \in Z} | v |^{(2 H - 2 p) 2 q + 2 q - m}$ converges we get

r_{N, q, 1, m} \leq C \frac{1}{N^{2 q - m}} \leq C \frac{1}{N} \to_{N \to \infty} 0

and when the series diverges,

r_{N, q, 1, m} \leq C \frac{1}{N^{2 q - m}} N^{(2 H - 2 p) 2 q + 2 q - m + 1} \leq C N^{(2 H - 2 p) 2 q + 1} \to_{N \to \infty} 0

if $p = 1, H < 1 - \frac{1}{4 q}$ , or $p \geq 2$ and $H \in (\frac{1}{2}, 1)$ . If $p = q = 1$ , $H = 3 / 4$ , the quantity

\frac{1}{log (N - l)} r_{N, q, 1, m}

will also converge to zero for $m = 0, 1$ using again (21) and Lemma 1.

The fact that the series $σ^{2} = \sum_{q \geq 1} σ_{2 q}^{2} < \infty$ for $H < p - \frac{1}{4 q}$ follows from the study of the $k$ -variations of the fractional Brownian motion, see [6] or [14].

We will consider the renormalized $k$ -variation sequence

G_{N} (k, α) = \sqrt{N - l} V_{N} (k, α) .

(23)

From the above Lemma 2, it follows that

E {[G_{N} (k, α)]}_{N}^{2} \to_{N \to \infty} σ^{2}

with $σ^{2}$ given in the statement of Lemma 2. We will show that the sequence (23) satisfies a central limit theorem.

Theorem 1

For a filter $α$ of order $p \geq 1$ and of length $l + 1 \geq 1$ , with $p > H + \frac{1}{4}$ , let $G_{N} (k, α)$ be given by (23). Then the sequence ${(G_{N} (k, α))}_{N \geq 1}$ converges in distribution, as $N \to \infty$ , to the Gaussian law $N (0, σ^{2})$ . Moreover, for $p = 1$ , $H = 3 / 4$ , the sequence ${(\frac{1}{\sqrt{log (N - l)}} G_{N} (k, α))}_{N \geq 1}$ converges in distribution to $N (0, c^{2})$ . The constants $σ^{2}, c^{2}$ are those appearing in Lemma 2.

Proof: Notice that from (13), we can write

G_{N} (k, α) = \sum q \geq 1 I_{2 q} (g_{N, 2 q}) with g_{N, 2 q} = \sqrt{N - l} f_{N, 2 q}

(24)

with $f_{N, 2 q}$ given by (14). Our main tool to prove the asymptotic normality of (24) is Theorem 6.3.1 from [14]. According to it, for $p > H + 1 / 4$ it suffices to show that

$(2 q)! ∥ g_{N, 2 q} ∥_{H^{\otimes 2 q}}^{2} \to_{N \to \infty} σ_{2 q}^{2}$ and $σ^{2} := \sum_{q \geq 1} σ_{2 q}^{2} < \infty,$
for every $q \geq 1$ and $r = 1, . ., 2 q - 1$ , $∥ g_{N, 2 q} \otimes_{r} g_{N, 2 q} ∥_{H^{\otimes 4 q - 2 r}} \to_{N \to \infty} 0,$

and for $p = 1$ , $H = 3 / 4$ ,

$\frac{1}{log (N - l)} (2 q)! ∥ g_{N, 2 q} ∥_{H^{\otimes 2 q}}^{2} \to_{N \to \infty} 1_{{q = 1}} c^{2}$ ,
for every $q \geq 1$ and $r = 1, . ., 2 q - 1$ , $\frac{1}{log (N - l)} ∥ g_{N, 2 q} \otimes_{r} g_{N, 2 q} ∥_{H^{\otimes 4 q - 2 r}} \to_{N \to \infty} 0,$
${lim}_{M \to \infty} {sup}_{N \geq 1} \sum_{q \geq M + 1} \frac{1}{log (N - l)} (2 q)! ∥ g_{N, 2 q} ∥_{H^{\otimes 2 q}}^{2} = 0$ .

Point 1. in both cases follows from Lemma 2. Let us check point 2. By the definition of the contraction (48), we have for $q \geq 1$ and $r = 1, . ., 2 q - 1$

g_{N, 2 q} \otimes_{r} g_{N, 2 q}

=

\frac{1}{N - l} \frac{(c_{2 q}^{k})^{2}}{(2 q)!} N \sum i, j = l ⟨ C_{i <}

Generalized k-variations and Hurst parameter estimation for the fractional wave equation via Malliavin calculus

Authors

On quadratic variations of the fractional-white wave equation

Statistical analysis of discretely sampled semilinear SPDEs: a power variation approach

Hurst index estimation in stochastic differential equations driven by fractional Brownian motion

Estimation Of all parameters in the Fractional Ornstein-Uhlenbeck model under discrete observations

Parameter estimation for discretely sampled stochastic heat equation driven by space-only noise revised

Parameter estimation for SPDEs based on discrete observations in time and space

Multi-dimensional parameter estimation of heavy-tailed moving averages

1 Introduction

2 Preliminaries

2.1 The solution to the wave equation with fractional-colored noise

2.2 Filters

Definition 1

Lemma 1

3 Central limit theorem for the spatial $k$ -variations

3.1 Central Limit Theorem

Lemma 2

Theorem 1

Generalized k-variations and Hurst parameter estimation for the fractional wave equation via Malliavin calculus

Authors

Related Research

On quadratic variations of the fractional-white wave equation

Statistical analysis of discretely sampled semilinear SPDEs: a power variation approach

Hurst index estimation in stochastic differential equations driven by fractional Brownian motion

Estimation Of all parameters in the Fractional Ornstein-Uhlenbeck model under discrete observations

Parameter estimation for discretely sampled stochastic heat equation driven by space-only noise revised

Parameter estimation for SPDEs based on discrete observations in time and space

Multi-dimensional parameter estimation of heavy-tailed moving averages

This week in AI

1 Introduction

2 Preliminaries

2.1 The solution to the wave equation with fractional-colored noise

2.2 Filters

Definition 1

Lemma 1

3 Central limit theorem for the spatial k-variations

3.1 Central Limit Theorem

Lemma 2

Theorem 1

3 Central limit theorem for the spatial $k$ -variations