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 -variations of the observed process. The method has been recently employed in [11] for the case of quadratic variations, i.e. . As for the fBm, it was shown that the standard quadratic variation estimator is not asymptotically normal when the Hurst index becomes bigger than and this is inconvenient for statistical applications. In order to avoid this restriction and to get an estimator which is asymptotically Gaussian for every , we will use the generalized -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
denotes the solution to the wave equation with fractional-white noise, we define the (centered) generalized
-variation statistics ( integer),(1) |
where represents the spatial increment of the solution at along a filter of power (order) and length (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 (in the spirit of [4]) whenever and in this way the restriction can be avoided by choosing a filter of order , 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 and to -variations of any order and in addition we show that in the special case and 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 -variation estimator for is consistent and asymptotically normal. Moreover, we provide a numerical analysis of the estimator when 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 and we obtain the rate of convergence when is even via the Stein-Malliavin theory. In Section 4, we show a non-central limit limit in the case and for filters of order . 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 -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
where is the Laplacian on , , and is a fractional-white Gaussian noise which is defined as a real valued centered Gaussian field , over a given complete filtered probability space , with covariance function given by
(2) |
where is the covariance of the fractional brownian motion
We denoted by the class of bounded Borel subsets of and we will assume throughout this work
The solution of the equation (2.1) is understood in the mild sense, that is, it is defined as a square-integrable centered field defined by
(3) |
where is the fundamental solution to the wave equation and the integral in (3) is a Wiener integral with respect to the Gaussian process . Recall that for (we will later restrict to this situation in our work) we have, for and ,
(4) |
We know (see e.g. [3]) that the solution (3) is well-defined if and only if
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
(5) | |||||
with . When and , this expression reduces to
(6) |
We will fix for the rest of the work and we will associate to the process its canonical Hilbert space which is defined as the closure of the linear space generated by the indicator functions with respect to the inner product
We will denote by the multiple stochastic integral of order with respect to the Gaussian process and by the Malliavin derivative with respect to this process. We refer to the Appendix for the basic elements of the Malliavin calculus.
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 and , a vector
with the convention .
For a filter of length and of order we define the space-filtered process (or the spatial increment of the process along the filter )
(7) |
We denote for
From the covariance formula (5) we can write
(8) | |||||
with
and and . We write for further use
We will need the below technical lemma to etablish the asymptotical equivalent of and similar expressions. The proof of the lemma is based on a Taylor expansion, see [6] or [10].
Lemma 1
Let and , be filters of lengths , and of orders respectively. Then
with , where .
In the sequel, we write to indicate that the sequences have the same behavior as .
3 Central limit theorem for the spatial -variations
In this section we focus on the asymptotic behavior in distribution of the -variation in space of the solution to the fractional-white wave equation, defined via a filter of power . In the first step we show the -variation satisfies a CLT when . Next, by taking 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 , we prove a multidimensional CLT, which is needed for the estimation of the Hurst parameter.
3.1 Central Limit Theorem
Fix and let be a filter of length and of power as in Definition 1. Let be given by (3). For any integer we define the centered spatial -variations of the process by
(10) |
with 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 -variation sequence . Noticing that the filtered process
as a linear combination of centered Gaussian random variables is a centered Gaussian process, we get
(11) |
where denotes the
-th absolute moment of a standard Gaussian variable given by
. We introduce the variable(12) |
It is clear that is a standard Gaussian variable and , where denotes the correlation coefficient. Using (11) and (12) we can write as follows:
Recall the expansion of the development in Hermite polynomials of the function given in Lemma 2 of [6],
where for , for and denotes the j-th Hermite polynomial defined by
Observing that for
we have from (8) that we can express as an integral with respect to the process since the increment can be expressed as (recall that represents the multiple integral of order 1 with respect to the Gaussian process ) for every :
Since we have for we get
Hence, we obtain the following chaotic expansion of the -variation statistics
(13) |
with
(14) |
Let us start by analyzing the asymptotic behavior of the mean square of each kernel .
Lemma 2
For , let be given by (14). Then
for (i.e. for and for ), where we use the notation
(15) |
Moreover, For ,
(16) |
and thus
(18) |
Using the expression
with
(19) |
we can write, with and given by (15) and (17) respectively,
(20) |
and remark that due to Lemma 1 for large enough
(21) |
where does not depend on .With this notation we can write
with
(22) |
Clearly, by the dominated convergence theorem
which by Lemma 1 is finite if , , and for all in the other cases. For , ,
converges to a positive constant and thus (16) is obtained.
In order to conclude it remains to show that the rest term (22) converges to as , for every . From (22), using the bound (21) and Lemma 1, we have the estimate
and for each ,
When the series converges we get
and when the series diverges,
if , or and . If , , the quantity
will also converge to zero for using again (21) and Lemma 1.
The fact that the series for follows from the study of the -variations of the fractional Brownian motion, see [6] or [14].
We will consider the renormalized -variation sequence
(23) |
From the above Lemma 2, it follows that
with given in the statement of Lemma 2. We will show that the sequence (23) satisfies a central limit theorem.
Theorem 1
Proof: Notice that from (13), we can write
(24) |
with given by (14). Our main tool to prove the asymptotic normality of (24) is Theorem 6.3.1 from [14]. According to it, for it suffices to show that
-
and
-
for every and ,
-
and for , ,
-
,
-
for every and ,
-
.
Point 1. in both cases follows from Lemma 2. Let us check point 2. By the definition of the contraction (48), we have for and
Comments
There are no comments yet.