多项式曲面拟合numpy

Question

如何z=f(x,y)在具有完整交叉项的numpy中使用多项式拟合2D曲面？

Answer 1

这本质上是数字病态的,但你可以做这样的事情:

import numpy as np

x = np.random.randn(500)
y = np.random.randn(500)
z = np.random.randn(500) # Dependent variable

v = np.array([np.ones(500), x, y, x**2, x * y, y**2])

coefficients, residues, rank, singval = np.linalg.lstsq(v.T, z)

你添加的术语越多,数字越糟糕.你确定要多项式插值吗？

对于多项式还有其他基础,其值的矩阵没有那么严格的条件,但我不记得它们被称为什么; 但是,任何大学水平的数值分析教科书都会有这种材料.