From Volterra Series to Kunchenko Stochastic Polynomials: Half a Century of Non-Gaussian Estimation Methodology

TL;DR

This paper traces the evolution of Kunchenko's semiparametric non-Gaussian estimation methods over fifty years, highlighting theoretical developments, practical algorithms, and applications in signal processing.

Contribution

It introduces Kunchenko stochastic polynomials as a unified framework for moment-based estimation, hypothesis testing, and signal decomposition, connecting historical theory with modern applications.

Findings

Kunchenko polynomials form a coherent family of moment-cumulant procedures.

The paper establishes a formal link between Volterra models and Kunchenko polynomials.

PMM efficiency depends on moment existence, matrix nondegeneracy, and variance reduction.

Abstract

This paper reconstructs the half-century evolution of the scientific school founded by Yuriy P. Kunchenko (1939--2006) as the development of a semiparametric methodology for non-Gaussian estimation. Starting with Kunchenko's 1972/1973 dissertation applying Volterra series to estimate parameters of random processes, the trajectory is followed through 2006--2026. Kunchenko stochastic polynomials are presented as a coherent family of moment-cumulant procedures: the polynomial maximization method (PMM) for parameter estimation, polynomial criteria for hypothesis testing, and decomposition in spaces with a generating element. The paper details the school's structure: a verified genealogy of 15 defended dissertations, collaborations in Poland, Slovakia, and Germany, and the R package EstemPMM. A recent 2026 paper on Volterra-based signal processing is analyzed, showing how Kunchenko's nonlinear formulation reappears in applied radio engineering. We build a formal bridge between finite Volterra models and generalized Kunchenko polynomials, while separating the MMSE/L2 criterion from PMM: the former is a covariance projection for kernel adaptation, whereas PMM is a parameter-dependent moment procedure. PMM efficiency claims are stated conditionally: gains require that moments exist, the centered correlant matrix is nondegenerate, and the variance reduction coefficient is below one. The concluding research program operationalizes the historical reconstruction into testable statistical and signal-processing tasks.