python中Chebyshev多项式的拟合

Question

我试图用多项式拟合2d数据，发现在一定程度的多项式之后，numpy给出了“等级警告”。另一方面，使用Chebyshev多项式进行拟合并不会发出这样的警告。我可以用numpy的Chebyshev来拟合同样的数据。输出为不同等级的Chebyshev多项式。

我要计算这些Chebyshev多项式的数值，以便与正规多项式拟合的结果进行比较。这是示例代码。

import numpy as np
from numpy.polynomial import Chebyshev as T
import os
from mpmath import *
mp.dps = 16 
mp.pretty = True

x = [1, 2, 3, 4, 5]
y = [0.90, 8.15, 26.84, 64.87, 124.0]

deg = 3
popt_poly = np.polyfit(x, y, deg)
popt_cheb = T.fit(x, y, deg)
ypred = np.polyval(popt_poly, x)
print(popt_poly)
print(ypred)
print(popt_cheb)

提前感谢你。

Answer 1

请使用"window=“参数。在本例中: window=1,5