python：numpy矩阵操作

转置矩阵

import numpy as np

a = np.array([[9, 12.802, -45],
            [12.802, 22.348, -79.057],
            [-45, -79.057, 285]])
print(a)
print("transpose matrix is\n", a.transpose())

逆矩阵

a_inverse = np.linalg.inv(a)
print("inverse matrix is\n", a_inverse)

行列式

a_determinant = np.linalg.det(a)
print("determinant is\n", a_determinant)

伴随矩阵

a_adj = a_determinant * a_inverse
print("adjoint matrix is\n", a_adj)

矩阵乘法

b = np.array([[29.289],
              [43.543],
              [-149.258]])
c = np.dot(a_inverse, b)
print(c)

用循环创建矩阵

from itertools import chain

x = []
temp = []
for i in range(1, 10):
    temp.append([1, np.log(i), -i])
x.append(temp)

print(x)
x = list(chain.from_iterable(x))
y = np.dot(x, c)
print("y", y)

在for循环中需要创建一个临时列表temp来暂时储存每一行的值，再用.append()方法添加到矩阵x中。

注意：
1.这种方法会多产生一层列表，可以用itertools.chain包去掉外层列表。
2. 产生的矩阵是python list而不是numpy array，可以用np.array()转换