2025-12-12 | | Total: 3
We develop here a stochastic framework for modeling and segmenting transient spindle- like oscillatory bursts in electroencephalogram (EEG) signals. At the modeling level, individ- ual spindles are represented as path realizations of a two-dimensional Ornstein{Uhlenbeck (OU) process with a stable focus, providing a low-dimensional stochastic dynamical sys- tem whose trajectories reproduce key morphological features of spindles, including their characteristic rise{decay amplitude envelopes. On the signal processing side, we propose a segmentation procedure based on Empirical Mode Decomposition (EMD) combined with the detection of a central extremum, which isolates single spindle events and yields a collection of oscillatory atoms. This construction enables a systematic statistical analysis of spindle features: we derive empirical laws for the distributions of amplitudes, inter-spindle intervals, and rise/decay durations, and show that these exhibit exponential tails consistent with the underlying OU dynamics. We further extend the model to a pair of weakly coupled OU processes with distinct natural frequencies, generating a stochastic mixture of slow, fast, and mixed spindles in random temporal order. The resulting framework provides a data- driven framework for the analysis of transient oscillations in EEG and, more generally, in nonstationary time series.
We approximate the first Dirichlet eigenpair of the $p$-Laplace operator for $2 \leq p < \infty$ on both Euclidean and surface domains. We emphasize large $p$ values and discuss how the $p \to \infty$ limit connects to the underlying geometry of our domain. Working with large $p$ values introduces significant numerical challenges. We present a surface finite element numerical scheme that combines a Newton inverse-power iteration with a new domain rescaling strategy, which enables stable computations for large $p$. Numerical experiments in $1$D, planar domains, and surfaces embedded in $\mathbb{R}^3$ demonstrate the accuracy and robustness of our approach and show convergence towards the $p \to \infty$ limiting behavior.
We investigate matrices satisfying the Kreiss condition $$\|(zI-T)^{-1}\|\le\cfrac{K}{|z|-1}, \hspace{0.7 cm} |z|>1, $$ with $K$ lying arbitrarily close to $1.$ We provide lower bounds for the power growth of such matrices, which complement and refine related estimates due to Nikolski and Spijker-Tracogna-Welfert. We also study operators that satisfy a variant of the above Kreiss condition where $K$ is replaced by $1+ε(|z|)$, where the positive continuous function $ε(|z|)$ tends to $0$ as $|z|\to 1^+.$ We show that, if the spectum of $T$ touches the unit circle only at a single point and the resolvent of $T$ satisfies a growth restriction along the unit circle, it is possible to choose $ε$ so that this Kreiss-type condition guarantees similarity to a contraction. At the core of our proof lies a positivity argument involving the double-layer potential operator. Counterexamples related to less restrictive choices of $ε$ are also provided.