想通过最小二乘法拟合二元多项式z=[x,sin(x),cos(x)]B[y;sin(y);cos(y)],B是3*3的系数矩阵
def fit_func(x,y,B):
U=np.array([x,np.sin(x),np.cos(x)])
V=np.array([x,np.sin(y),np.cos(y)])
V_t=V.T
A=np.dot(U,B)
f=np.dot(A,V_t)
return f
Xi=np.array([8.19,2.72,6.39,8.71,4.7,2.66,3.78])
Yi=np.array([7.01,2.78,6.47,6.71,4.1,4.23,4.05])
Zi=np.array(np.random.rand(7))
错误信息:
A=np.dot(U,B)
ValueError: shapes (4,16) and (7,) not aligned: 16 (dim 1) != 7 (dim 0)