"""
利用正规方程求解向量θ
θ为n+1个参数组成的(列)向量
θ=[θ_0,θ_1,θ_2,...,θ_n].T T:转置
"""
import numpy as np
def normalEqn(X, y):
"""
numpy.linalg模块包含线性代数的函数。
使用这个模块,可以计算逆矩阵、求特征值、解线性方程组以及求解行列式等
:param X: 特征矩阵
:param y: 结果向量
:return: 参数向量θ
"""
theta = np.linalg.inv(X.T @ X) @ X.T @ y # inv(A):inverse of matrix A.矩阵A的逆 ;X.T @ X等价于X.T.dot(X)
return theta
if __name__ == '__main__':
XX = np.array([[1, 2], [3, 4]])
yy = np.array([1, 1])
k = normalEqn(XX, yy)
print(k)