python数据分析——numpy数组学习（3）

矩阵的计算

1、矩阵乘法

扩展库numpy中，支持矩阵的乘法运算，可直接计算

例：

import numpy as np
x=np.matrix([1,2,3],[4,5,6])
y=np.matrix([[1,2],[3,4],[5,6]])
print(x*y)

结果：

[[22 28]
[49 64]]

2、计算相关系数矩阵

扩展库numpy提供了corrcoef（）计算相关系数矩阵

例：

import numpy as np
print(np.corrcoef([1,2,3,4],[4,3,2,1]))
print(np.corrcoef([1,2,3,4],[8,3,2,1]))
print(np.corrcoef([1,2,3,4],[1,2,3,4]))
print(np.corrcoef([1,2,3,4],[1,2,3,5]))

结果：

[[ 1. -1.]
[-1. 1.]]
[[ 1. -0.91350028]
[-0.91350028 1. ]]
[[1. 1.]
[1. 1.]]
[[1. 0.98270763]
[0.98270763 1. ]]

3、计算逆矩阵

扩展库numpy的线性代数子模块linalg中提供了用来计算逆矩阵的函数inv（）

例：

import numpy as np
x=np.matrix([[1,2,3],[4,5,6,],[7,8,9]])
y=np.linalg.inv(x)
print(y)
print(x*y)
print(y*x)

结果：

[[ 3.15251974e+15 -6.30503948e+15 3.15251974e+15]
[-6.30503948e+15 1.26100790e+16 -6.30503948e+15]
[ 3.15251974e+15 -6.30503948e+15 3.15251974e+15]]
[[ 0. 1. -0.5]
[ 0. 2. -1. ]
[ 0. 3. 2.5]]
[[ 12. 16. 12. ]
[-10. -16. -6. ]
[ 1.5 4. 2.5]]