多项式表面拟合numpy

Question

如何使用满交叉项的numpy多项式拟合2D表面z=f(x,y)？

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) 

您添加的更多方面，更糟糕的事情搞定，数值。你确定你想要一个多项式插值？

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