Track papers across labs, topics, and trends.
Search papers, labs, and topics across Lattice.
This paper introduces coded Hankel polynomial chaos (CH-PC), a novel spectral approach for identifying dominant polynomial-chaos modes that enhances traditional sparse-regression methods. By utilizing a finite generating transform and low-rank Hankel matrices, the authors achieve exact recovery and noise stabilization in polynomial chaos expansions, even under conditions of sampling and observation errors. Numerical experiments demonstrate the method's efficacy in recovering dominant modes from both sparse Legendre benchmarks and a stochastic Darcy problem, showcasing its robustness in practical applications.
Coded Hankel polynomial chaos can recover dominant modes with unprecedented accuracy, even in the presence of noise and sampling errors.
Identification of dominant polynomial-chaos modes is usually formulated as a sparse-regression problem on a sampled multivariate polynomial dictionary. We develop coded Hankel polynomial chaos (CH-PC), a complementary spectral formulation for dominant-mode identification. A finite generating transform converts PCE coefficients into a coefficient-generating polynomial, and evaluation along a geometric phase orbit produces a finite exponential sum. Its model order and spectral nodes are encoded by low-rank Hankel matrices, while coordinate phase shifts attach root-of-unity labels from which the full polynomial multi-indices are recovered. Coordinate-shifted probes are combined as common-node snapshots, and independent phase encodings provide redundant representations when a single spectral encoding is poorly conditioned. For finite observations, population, finite-data, and observed probes are kept distinct: sampling or quadrature error and observation error enter as separate Hankel perturbations, which are then connected to spectral stability, discrete decoding, and phase voting. For tensor-product candidate sets, the generating kernel factorizes into one-dimensional sums and can be evaluated without assembling the full multivariate PCE design matrix. Numerical experiments on sparse Legendre benchmarks and a stochastic Darcy problem illustrate exact recovery, noise stabilization, unknown-order identification by phase persistence, and dominant-mode recovery for a PDE-generated quantity of interest.