Python库Numpy的代数计算

Numpy库，可以用二维数据代替矩阵，举例如下：

x = np.array([[1., 2., 3.], [4., 5., 6.]])
y = np.array([[6., 23.], [-1, 7], [8, 9]])
x.dot(y) # 等价于np.dot(x, y)，x的轴1和y的轴0的shape相等可相乘，所以x.dot(ones(3))是可行的，结果是array([6., 15.])

from numpy.linalg import inv, qr
x = np.random.randn(3, 3); mat = x.T @ x
inv(mat) # 求方阵mat的逆（逆矩阵）
q, r = qr(mat) # QR分解
np.mat(x) # 生成矩阵，建议在numpy中使用数组运算，如果转为矩阵经常需要转来转去，因为处理结果往往都是数组。