今天在阅读数学建模的时候看到了差分那章
其中有一个用matlab求线性的代码,这里我贴出来
这里我送上 Python代码
In [39]: import numpy as np
...: from scipy.optimize import nnls
...: x = np.array([[1,2,3,4,5],[1,1,1,1,1]])
...: x = x.T
...: y = np.array([11,12,13,15,16])
...: nnls(x,y)
...:
Out[39]: (array([ 1.3, 9.5]), 0.5477225575051656)
In [40]: np.linalg.lstsq(x,y)
Out[40]: (array([ 1.3, 9.5]), array([ 0.3]), 2, array([ 7.69121313, 0.91936964]))
In [41]: