摘要 Abstract
有理函数逼近是一种简单但灵活的替代多项式逼近的方法,能够在不产生振荡伪影的情况下捕捉复杂的非线性特性。然而,由于可能产生虚假奇点,针对噪声数据使用有理函数的尝试很少。为了避免奇点的出现,我们利用Bernstein多项式及其系数的适当条件,强制分母严格为正。虽然这减少了可表达的有理多项式的范围,但它保留了有理函数的所有优点,同时在噪声数据场景中保持了多项式逼近的稳健性。我们的数值实验表明,现有的有理逼近方法在逼近域内持续产生虚假极点,而我们的方法不会在逼近域内产生极点,并且对具有多个变量的函数提供的拟合效果优于多项式逼近甚至惩罚样条。此外,在使用谱方法数值求解微分方程时,保证区间内无极点对于估计非恒定系数至关重要,这可以紧凑地表示原始微分方程,使数值求解器能够快速达到高精度,正如我们的实验所示。
Rational function approximations provide a simple but flexible alternative to polynomial approximation, allowing one to capture complex non-linearities without oscillatory artifacts. However, there have been few attempts to use rational functions on noisy data due to the likelihood of creating spurious singularities. To avoid the creation of singularities, we use Bernstein polynomials and appropriate conditions on their coefficients to force the denominator to be strictly positive. While this reduces the range of rational polynomials that can be expressed, it keeps all the benefits of rational functions while maintaining the robustness of polynomial approximation in noisy data scenarios. Our numerical experiments on noisy data show that existing rational approximation methods continually produce spurious poles inside the approximation domain. This contrasts our method, which cannot create poles in the approximation domain and provides better fits than a polynomial approximation and even penalized splines on functions with multiple variables. Moreover, guaranteeing pole-free in an interval is critical for estimating non-constant coefficients when numerically solving differential equations using spectral methods. This provides a compact representation of the original differential equation, allowing numeric solvers to achieve high accuracy quickly, as seen in our experiments.