Distinguishing between fractional Brownian motion with random and constant Hurst exponent using sample autocovariance-based statistics.

Fractional Brownian motion (FBM) is a canonical model for describing dynamics in various complex systems. It is characterized by the Hurst exponent, which is responsible for the correlation between FBM increments, its self-similarity property, and anomalous diffusion behavior. However, recent research indicates that the classical model may be insufficient in describing experimental observations when the anomalous diffusion exponent varies from trajectory to trajectory. As a result, modifications of the classical FBM have been considered in the literature, with a natural extension being the FBM with a random Hurst exponent. In this paper, we discuss the problem of distinguishing between two models: (i) FBM with the constant Hurst exponent and (ii) FBM with random Hurst exponent, by analyzing the probabilistic properties of statistics represented by the quadratic forms. These statistics have recently found application in Gaussian processes and have proven to serve as efficient tools for hypothesis testing. Here, we examine two statistics-the sample autocovariance function and the empirical anomaly measure-utilizing the correlation properties of the considered models. Based on these statistics, we introduce a testing procedure to differentiate between the two models. We present analytical and simulation results considering the two-point and beta distributions as exemplary distributions of the random Hurst exponent. Finally, to demonstrate the utility of the presented methodology, we analyze real-world datasets from the financial market and single particle tracking experiment in biological gels.

Testing of two-dimensional Gaussian processes by sample cross-covariance function.

In this paper, we address the issue of testing two-dimensional Gaussian processes with a defined cross-dependency structure. Multivariate Gaussian processes are widely used in various applications; therefore, it is essential to identify the theoretical model that accurately describes the data. While it is relatively straightforward to do so in a one-dimensional case, analyzing multi-dimensional vectors requires considering the dependency between the components, which can significantly affect the efficiency of statistical methods. The testing methodology presented in this paper is based on the sample cross-covariance function and can be considered a natural generalization of the approach recently proposed for testing one-dimensional Gaussian processes based on the sample autocovariance function. We verify the efficiency of this procedure on three classes of two-dimensional Gaussian processes: Brownian motion, fractional Brownian motion, and two-dimensional autoregressive discrete-time process. The simulation results clearly demonstrate the effectiveness of the testing methodology, even for small sample sizes. The theoretical and simulation results are supported by analyzing two-dimensional real-time series that describe the main risk factors of a mining company, namely, copper price and exchange rates (USDPLN). We believe that the introduced methodology is intuitive and relatively simple to implement, and thus, it can be applied in many real-world scenarios where multi-dimensional data are examined.

Divergence-Based Segmentation Algorithm for Heavy-Tailed Acoustic Signals with Time-Varying Characteristics.

Many real-world systems change their parameters during the operation. Thus, before the analysis of the data, there is a need to divide the raw signal into parts that can be considered as homogeneous segments. In this paper, we propose a segmentation procedure that can be applied for the signal with time-varying characteristics. Moreover, we assume that the examined signal exhibits impulsive behavior, thus it corresponds to the so-called heavy-tailed class of distributions. Due to the specific behavior of the data, classical algorithms known from the literature cannot be used directly in the segmentation procedure. In the considered case, the transition between parts corresponding to homogeneous segments is smooth and non-linear. This causes that the segmentation algorithm is more complex than in the classical case. We propose to apply the divergence measures that are based on the distance between the probability density functions for the two examined distributions. The novel segmentation algorithm is applied to real acoustic signals acquired during coffee grinding. Justification of the methodology has been performed experimentally and using Monte-Carlo simulations for data from the model with heavy-tailed distribution (here the stable distribution) with time-varying parameters. Although the methodology is demonstrated for a specific case, it can be extended to any process with time-changing characteristics.

Identification, Decomposition and Segmentation of Impulsive Vibration Signals with Deterministic Components-A Sieving Screen Case Study.

Condition monitoring is a well-established field of research; however, for industrial applications, one may find some challenges. They are mostly related to complex design, a specific process performed by the machine, time-varying load/speed conditions, and the presence of non-Gaussian noise. A procedure for vibration analysis from the sieving screen used in the raw material industry is proposed in the paper. It is more for pre-processing than the damage detection procedure. The idea presented here is related to identification and extraction of two main types of components: (i) deterministic (D)-related to the unbalanced shaft(s) and (ii) high amplitude, impulsive component randomly (R) appeared in the vibration due to pieces of ore falling down of moving along the deck. If we could identify these components, then we will be able to perform classical diagnostic procedures for local damage detection in rolling element bearing. As deterministic component may be AM/FM modulated and each impulse may appear with different amplitude and damping, there is a need for an automatic procedure. We propose a method for signal processing that covers two main steps: (a) related to R/D decomposition and including signal segmentation to neglect AM/FM modulations, iterative sine wave fitting using the least square method (for each segment), signal filtering technique by subtraction fitted sine from the raw signal, the definition of the criterion to stop iteration by residuals analysis, (b) impulse segmentation and description (beginning, end, max amplitude) that contains: detection of the number of impulses in a decomposed random part of the raw signal, detection of the max value of each impulse, statistical analysis (probability density function) of max value to find regime-switching), modeling of the envelope of each impulse for samples that protrude from the signal, extrapolation (forecasting) envelope shape for samples hidden in the signal. The procedure is explained using simulated and real data. Each step is very easy to implement and interpret thus the method may be used in practice in a commercial system.