本文介绍了使用Python的Numpy库进行矩阵的基本运算
##创建全0矩阵1
2
3# 创建3x5的全0矩阵
myZero = np.zeros([3, 5])
print myZero
输出结果：1
2
3[[ 0. 0. 0. 0. 0.]
[ 0. 0. 0. 0. 0.]
[ 0. 0. 0. 0. 0.]]
##创建全1矩阵1
2
3# 创建3x5的全1矩阵
myOnes = np.ones([3, 5])
print myOnes
输出结果：1
2
3[[ 1. 1. 1. 1. 1.]
[ 1. 1. 1. 1. 1.]
[ 1. 1. 1. 1. 1.]]
##创建0~1之间的随机矩阵1
2
3# 3x4的0~1之间的随机数矩阵
myRand = np.random.rand(3, 4)
print myRand
输出结果为：1
2
3[[ 0.26845651 0.26713961 0.12632736 0.69840295]
[ 0.92745819 0.44091417 0.21733213 0.76135785]
[ 0.97161283 0.13570203 0.07819361 0.72129986]]
##创建单位矩阵1
2
3# 3x3的单位矩阵
myEye = np.eye(3)
print myEye
输出结果为：1
2
3[[ 1. 0. 0.]
[ 0. 1. 0.]
[ 0. 0. 1.]]
##矩阵求和1
print myZero + myOnes
输出结果为：1
2
3[[ 1. 1. 1. 1. 1.]
[ 1. 1. 1. 1. 1.]
[ 1. 1. 1. 1. 1.]]
##矩阵求差1
print myZero - myOnes
输出结果为：1
2
3[[-1. -1. -1. -1. -1.]
[-1. -1. -1. -1. -1.]
[-1. -1. -1. -1. -1.]]
##创建矩阵1
2myMatrix = np.mat([[1, 2, 3], [4, 5, 6], [7, 8, 9]])
print myMatrix
输出结果为：1
2
3[[1 2 3]
[4 5 6]
[7 8 9]]
##矩阵乘以常数1
print 10 * myMatrix
输出结果为：1
2
3[[10 20 30]
[40 50 60]
[70 80 90]]
##矩阵所有元素求和1
print np.sum(myMatrix)
输出结果为45
##矩阵乘法
当维度相同的时候，为各个位置对应元素的乘积
当矩阵的维度不同时，会根据一定的广播规则将维数扩充到一致的形式1
2
3
4
5
6myMatrix2 = 2 * myMatrix
print myMatrix2
print np.multiply(myMatrix, myMatrix2)
myMatrix3 = np.mat([[1], [2], [3]])
print(myMatrix3)
print myMatrix * myMatrix3
输出结果为：1
2
3
4
5
6
7
8
9
10
11
12[[ 2 4 6]
[ 8 10 12]
[14 16 18]]
[[ 2 8 18]
[ 32 50 72]
[ 98 128 162]]
[[1]
[2]
[3]]
[[14]
[32]
[50]]
##矩阵的幂1
print np.power(myMatrix, 2)
输出结果为：1
2
3[[ 1 4 9]
[16 25 36]
[49 64 81]]
##矩阵的转置1
2
3print myMatrix
print myMatrix.T
print myMatrix.transpose()
输出结果为：1
2
3
4
5
6
7
8
9[[1 2 3]
[4 5 6]
[7 8 9]]
[[1 4 7]
[2 5 8]
[3 6 9]]
[[1 4 7]
[2 5 8]
[3 6 9]]
##矩阵的其他操作：行列数、切片、复制、比较1
2
3
4
5
6
7
8
9
10
11[m, n] = myMatrix.shape
print "矩阵的行列数为:", m, n
# 按照行切片（输出矩阵的行）
print myMatrix[0]
#按照列切片 (输出矩阵的列)
print myMatrix.T[0]
#矩阵的复制
myMatrixCopy = myMatrix.copy()
print myMatrixCopy
#矩阵的比较
print myMatrix < myMatrix.T
输出结果为：1
2
3
4
5
6
7
8
9矩阵的行列数为: 3 3
[[1 2 3]]
[[1 4 7]]
[[1 2 3]
[4 5 6]
[7 8 9]]
[[False True True]
[False False True]
[False False False]]
本文章全部代码为：1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52# -*- coding: utf-8 -*-
import numpy as np
# 创建3x5的全0矩阵
myZero = np.zeros([3, 5])
print myZero
# 创建3x5的全1矩阵
myOnes = np.ones([3, 5])
print myOnes
# 3x4的0~1之间的随机数矩阵
myRand = np.random.rand(3, 4)
print myRand
# 3x3的单位矩阵
myEye = np.eye(3)
print myEye
# 矩阵求和
print myZero + myOnes
# 矩阵求差
print myZero - myOnes
# 创建矩阵
myMatrix = np.mat([[1, 2, 3], [4, 5, 6], [7, 8, 9]])
print myMatrix
# 矩阵乘以常数
print 10 * myMatrix
# 矩阵所有元素求和
print np.sum(myMatrix)
# 矩阵乘法
# 当维度相同的时候，为各个位置对应元素的乘积
# 当矩阵的维度不同时，会根据一定的广播规则将维数扩充到一致的形式
myMatrix2 = 2 * myMatrix
print myMatrix2
print np.multiply(myMatrix, myMatrix2)
# 矩阵的幂
print np.power(myMatrix, 2)
myMatrix3 = np.mat([[1], [2], [3]])
print(myMatrix3)
print myMatrix * myMatrix3
#矩阵的转置
print myMatrix
print myMatrix.T
print myMatrix.transpose()
##矩阵的其他操作：行列数、切片、复制、比较
[m, n] = myMatrix.shape
print "矩阵的行列数为:", m, n
# 按照行切片（输出矩阵的行）
print myMatrix[0]
#按照列切片 (输出矩阵的列)
print myMatrix.T[0]
#矩阵的复制
myMatrixCopy = myMatrix.copy()
print myMatrixCopy
#矩阵的比较
print myMatrix < myMatrix.T