No Arabic abstract
In applications spaning from image analysis and speech recognition, to energy dissipation in turbulence and time-to failure of fatigued materials, researchers and engineers want to calculate how often a stochastic observable crosses a specific level, such as zero. At first glance this problem looks simple, but it is in fact theoretically very challenging. And therefore, few exact results exist. One exception is the celebrated Rice formula that gives the mean number of zero-crossings in a fixed time interval of a zero-mean Gaussian stationary processes. In this study we use the so-called Independent Interval Approximation to go beyond Rices result and derive analytic expressions for all higher-order zero-crossing cumulants and moments. Our results agrees well with simulations for the non-Markovian autoregressive model.
The study of correlated time-series is ubiquitous in statistical analysis, and the matrix decomposition of the cross-correlations between time series is a universal tool to extract the principal patterns of behavior in a wide range of complex systems. Despite this fact, no general result is known for the statistics of eigenvectors of the cross-correlations of correlated time-series. Here we use supersymmetric theory to provide novel analytical results that will serve as a benchmark for the study of correlated signals for a vast community of researchers.
In this study, we analytically formulated the path integral representation of the conditional probabilities for non-Markovian kinetic processes in terms of the free energy of the thermodynamic system. We carry out analytically the time-fractional kinetic equations for these processes. Thus, in a simple way, we generalize path integral solutions of the Markovian to the non-Markovian cases. We conclude that these pedagogical results can be applied to some physical problems such as the deformed ion channels, internet networks and non-equilibrium phase transition problems.
The persistence of a stochastic variable is the probability that it does not cross a given level during a fixed time interval. Although persistence is a simple concept to understand, it is in general hard to calculate. Here we consider zero mean Gaussian stationary processes in discrete time $n$. Few results are known for the persistence $P_0(n)$ in discrete time, except the large time behavior which is characterized by the nontrivial constant $theta$ through $P_0(n)sim theta^n$. Using a modified version of the Independent Interval Approximation (IIA) that we developed before, we are able to calculate $P_0(n)$ analytically in $z$-transform space in terms of the autocorrelation function $A(n)$. If $A(n)to0$ as $ntoinfty$, we extract $theta$ numerically, while if $A(n)=0$, for finite $n>N$, we find $theta$ exactly (within the IIA). We apply our results to three special cases: the nearest neighbor-correlated first order moving average process where $A(n)=0$ for $ n>1$, the double exponential-correlated second order autoregressive process where $A(n)=c_1lambda_1^n+c_2lambda_2^n$, and power law-correlated variables where $A(n)sim n^{-mu}$. Apart from the power-law case when $mu<5$, we find excellent agreement with simulations.
A nonparametric method to predict non-Markovian time series of partially observed dynamics is developed. The prediction problem we consider is a supervised learning task of finding a regression function that takes a delay embedded observable to the observable at a future time. When delay embedding theory is applicable, the proposed regression function is a consistent estimator of the flow map induced by the delay embedding. Furthermore, the corresponding Mori-Zwanzig equation governing the evolution of the observable simplifies to only a Markovian term, represented by the regression function. We realize this supervised learning task with a class of kernel-based linear estimators, the kernel analog forecast (KAF), which are consistent in the limit of large data. In a scenario with a high-dimensional covariate space, we employ a Markovian kernel smoothing method which is computationally cheaper than the Nystrom projection method for realizing KAF. In addition to the guaranteed theoretical convergence, we numerically demonstrate the effectiveness of this approach on higher-dimensional problems where the relevant kernel features are difficult to capture with the Nystrom method. Given noisy training data, we propose a nonparametric smoother as a de-noising method. Numerically, we show that the proposed smoother is more accurate than EnKF and 4Dvar in de-noising signals corrupted by independent (but not necessarily identically distributed) noise, even if the smoother is constructed using a data set corrupted by white noise. We show skillful prediction using the KAF constructed from the denoised data.
Recent experiments using fluorescence spectroscopy have been able to probe the dynamics of conformational fluctuations in proteins. The fluctuations are Gaussian but do not decay exponentially, and are therefore, non-Markovian. We present a theory where non-Markovian fluctuation dynamics emerges naturally from the superposition of the Markovian fluctuations of the normal modes of the protein. A Rouse-like dynamics of the normal modes provides very good agreement to the experimentally measured correlation functions. We provide simple scaling arguments rationalising our results.