Automatic Determination of Set of Multivariate Basis Polynomials Based on Recursion and Application of LSPCR to High-dimensional Uncertainty Quantification of Multi-conductor Transmission Lines

Weiwei Chen,Yuhang Ji,Ziyang Qi,Liping Yan,Xiang Zhao
DOI: https://doi.org/10.1109/access.2023.3348132
IF: 3.9
2024-01-01
IEEE Access
Abstract:Uncertainty quantification (UQ) in polynomial chaos expansion (PCE) suffers from the curse of dimensionality, which is fundamentally reflected in that the number of the PCE coefficients to be found grows rapidly as the number of random inputs increases, and thus the number of samples required increases dramatically. Regardless of the approach taken to alleviate this disaster, we are usually faced with the technical challenge of automatically and programmatically determining the basis functions during the PCE implementation, which are composed of a large number of multivariate polynomials. To address this problem, this paper proposes an algorithm based on the recursive idea for automatically determining the basis functions, and combines this algorithm with the latest weight-based regression called least squares polynomial chaos regression (LSPCR), proposes the implementation of the LSPCR method in the medium- and high-dimensional cases, and evaluates its performance in the electromagnetic coupling problem of the transmission lines (TLs) in this case. Subsequently, the performance of six algorithms for LSPCR are compared in the medium- and high-dimensional, low-order cases, each algorithm consisting of a sampling strategy selected from asymptotic sampling (AS), standard sampling (SS), and coherence-optimal sampling (COS), paired with a norm problem chosen from either least squares optimization (LSO) or ℓ 1 -minimization (ℓ 1 -M) problems. Numerical experiments demonstrate the excellence of the proposed algorithm. Furthermore, the results also show that only the algorithms that employ the SS (SSs) and the algorithms based on the COS (COSs) require no more than 1.4% of the total Monte Carlo (MC) computation time when producing similarly accurate results to the MC, regardless of the chosen norm problem.
What problem does this paper attempt to address?