代码如下:
#encoding:utf-8 import numpy as np from numpy.linalg import inv from numpy import dot from numpy import mat X=mat([1,2,3]).reshape(3,1) Y=2*X#此实验的函数模型(最简单的模型) print(X) #theta=(X.T*X)^-1*X.T*Y theta=dot(dot(inv(dot(X.T,X)),X.T),Y) print(theta)
代码如下:
#encoding:utf-8 import numpy as np from numpy.linalg import inv from numpy import dot from numpy import mat X=mat([1,2,3]).reshape(3,1) Y=2*X#此实验的函数模型(最简单的模型) print(X) #theta=(X.T*X)^-1*X.T*Y theta=dot(dot(inv(dot(X.T,X)),X.T),Y) print(theta)