今天在閱讀數學建模的時候看到了差分那章
其中有一個用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]:
