1.快速構造一個矩陣

• 構造一個3*5的全為0的矩陣

a = np.zeros(shape=(3, 5), dtype=int)
print(a)
[[0 0 0 0 0]
 [0 0 0 0 0]
 [0 0 0 0 0]]

• 構造一個2*4的全為1的矩陣

a = np.ones(shape=(2, 4), dtype=int)
print(a)
[[1 1 1 1]
 [1 1 1 1]]

• 構造一個全為666的矩陣

a = np.full((3, 5), 666)
print(a)
[[666 666 666 666 666]
 [666 666 666 666 666]
 [666 666 666 666 666]]

• 構造一個3*5，元素是[0透绩，10）之間隨機數(shù)的矩陣

a = np.random.randint(0, 10, (3, 5))
print(a)
[[7 5 0 9 3]
 [7 3 9 8 7]
 [9 1 6 9 7]]

• 構造一個3*5，所有元素符合均值為0壁熄，方差為1的正態(tài)分布的隨機數(shù)的矩陣

a = np.random.normal(0, 1, (3, 5))
print(a)
[[ 0.43121389 -1.45724221  0.84369408 -1.62208387  0.42111614]
 [-2.70974994  0.46920864  1.48373216 -1.72006643 -0.29006381]
 [-0.70056842  0.22743593  0.6276454  -0.78630736 -1.17294585]]

• 生成一個對角陣

b = np.mat(np.eye(3, 3, dtype=int))
print(b)
# 生成一個3*3的對角矩陣
[[1 0 0]
 [0 1 0]
 [0 0 1]]

• 生成一個對角線為指定數(shù)的對角矩陣

b = np.mat(np.diag([1, 2, 3]))
print(b)
[[1 0 0]
 [0 2 0]
 [0 0 3]]

2.矩陣的常用操作

a = np.mat([[1, 1, 1], [2, 3, 6], [4, 5, 7]])
[[1 1 1]
 [2 3 6]
 [4 5 7]]

• 常用的查詢方法

print(a.ndim)  # 查詢幾維數(shù)組 2
print(a.shape)  # 查詢幾行幾列 (3, 3)
print(a.size)  # 查詢元素個數(shù) 9
print(a[1, 1])  # 輸出1行1列（0開始算起）元素 3
print(a.max())  # 計算矩陣中最大元素 7
print(max(a[:, 1]))  # 計算第一列中的最大值,得到的是矩陣 [[5]]
print(a[1, :].max())  # 計算第二行的最大值帚豪，得到一個數(shù) 6
print(np.max(a, 0))  # 計算所有列的最大值
[[4 5 7]]
print(np.max(a, 1))  # 計算所有行的最大值
[[1]
 [6]
 [7]]
print(np.argmax(a, 0))  # 求所有列的最大值的索引
[[2 2 2]]
print(np.argmax(a[1, :]))  # 計算第二行的最大值在該行的索引 2

• 矩陣的拆分

print(a[:2, :2])  # 從原矩陣中分割出（0-2）*（0-2）的子矩陣
[[1 1]
 [2 3]]

• 矩陣的合并

d1 = np.mat(np.ones((2, 2)))
d2 = np.mat(np.eye(2))
print(d1, d2)
[[1. 1.]
 [1. 1.]]
 [[1. 0.]
 [0. 1.]]
d3 = np.vstack((d1, d2))  # 兩個矩陣按列合并
print(d3)
[[1. 1.]
 [1. 1.]
 [1. 0.]
 [0. 1.]]
d4 = np.hstack((d1, d2))
print(d4)  # 按行合并
[[1. 1. 1. 0.]
 [1. 1. 0. 1.]]

3.矩陣的運算

• 一個1n矩陣與一個n1矩陣相乘

b1 = np.mat([1, 2, 3])
b2 = np.mat([[3], [2], [1]])
print(b1, b2)
b = b1 * b2
print(b)
[[1 2 3]]
 [[3]
 [2]
 [1]]
[[10]]

• 矩陣的點乘

b1 = np.mat([1, 2])
b2 = np.mat([3, 4])
print(b1, b2)
b = np.multiply(b1, b2)
print(b)
[[1 2]] [[3 4]]
[[3 8]]

• 矩陣與數(shù)相乘

b1 = np.mat([1, 1])
b = b1 * 2
print(b)
# 矩陣與數(shù)的點乘
[[2 2]]

• 使用mat方法構建的矩陣求逆矩陣

c1 = np.mat(np.eye(2, 2)*0.5)  # 構建一個2行2列的對角矩陣，元素為0.5
print(c1)
c = c1.I
print(c)
[[0.5 0. ]
[0.  0.5]]
[[2. 0.]
[0. 2.]]

• 使用array方法構建的矩陣求逆矩陣

c1 = np.array([[1, 2], [3, 4]])
c = np.linalg.inv(c1)
print(c)
[[-2.   1. ]
 [ 1.5 -0.5]]

• 使用mat方法構建的矩陣求轉置

c1 = np.mat([[1, 2], [3, 4]])
c = c1.T
print(c)
[[1 3]
 [2 4]]

• 使用array方法構建的矩陣求轉置

c1 = np.array([[1, 2], [3, 4]])
c = c1.transpose()
print(c)
[[1 3]
 [2 4]]

• 求行列式的值

c = np.array([[1, 2], [3, 4]])
print(np.linalg.det(c))
-2.0000000000000004

• 求矩陣的特征值與特征向量

c = np.array([[1, 2], [3, 4]])
print(np.linalg.eig(c))
# 所得的元組中草丧，第一個為特征值元組狸臣，第二個為相對應的特征向量
(array([-0.37228132,  5.37228132]), array([[-0.82456484, -0.41597356],
       [ 0.56576746, -0.90937671]]))

• 求解線性方程組

c = np.array([[1, 2], [3, 4]])
d = np.array([[5], [10]])
print(np.linalg.solve(c, d))
"""
求解線性方程組
1X + 2Y = 5
3X + 4Y = 10
"""
[[0. ]
 [2.5]]

• 求列之和與行之和

c = np.mat([[1, 1], [2, 3], [4, 5]])
c1 = c.sum(axis=0)
print(c1)
# 求列和，得到1*2的矩陣
c2 = c.sum(axis=1)
print(c2)
# 求行和昌执，得到3*1的矩陣
[[7 9]]
[[2]
 [5]
 [9]]