问Python如何实现一个系数矩阵的拟合？

一个矩阵存有1200个4*4小图块的灰度直方图，每个小图块对应8个灰度级，储存在a.shape=(1200,8)里。

其中有一部分像素点是前景1，一部分是背景0，储存在c.shape=(1200,1)里。

现在需要求一个b.shape=(8,1)的矩阵，使np.dot(a,b)=c，但是我把a和c带入curve_fit(fun,a,c)拟合函数后报错，不知道是不是函数输入只能是一维矩阵。

求大佬指点，这个该怎么写呀？