python 稀疏矩阵的计算

1 稀疏矩阵的学习

1. SciPy教程 - 稀疏矩阵库scipy.sparse
2. 稀疏矩阵运算库
3. 官方文档

2 稀疏矩阵的实现

code

# coding = utf-8
from scipy.sparse import rand
'''生成随机稀疏矩阵'''
A = rand(4, 4, density=0.25, format="csc", random_state=42)
B = rand(4, 4, density=0.25, format="csc", random_state=40)
print(A.todense())
print(B.todense())
print((A.todense())*(B.todense()))

'''生成csc型的稀疏矩阵'''
import numpy as np
from scipy.sparse import csc_matrix,find
row = np.array([0,2,2,0,1,2])
col = np.array([0,0,1,2,2,2])
data1 = np.array([1,2,3,4,5,6])
data2 = np.array([10,20,30,40,50,60])
X = csc_matrix((data1, (row, col)), shape=(3, 3))
Y = csc_matrix((data2, (row, col)), shape=(3, 3))
print(X.toarray())
QQ = find(X) #返回非零元素的位置和值
print(QQ)
print((X-(X.T)).toarray())
print(Y.toarray())
print('稀疏矩阵中矩阵的乘法')
print((X*Y).toarray())  # 稀疏矩阵中矩阵的乘法
print('稀疏矩阵中对应元素的乘法')
print((csc_matrix.multiply(X,Y).toarray())) # 稀疏矩阵中对应元素的乘法
print('稀疏矩阵中矩阵的乘法')
print((csc_matrix.dot(X,Y)).toarray())  # 稀疏矩阵中矩阵的乘法
print('稀疏矩阵中对应元素的加法')
print((X+Y).toarray())  # 稀疏矩阵中对应元素的加法
print('稀疏矩阵中对应元素的2次方')
print((csc_matrix.power((X-Y),2)).toarray())    # 稀疏矩阵中对应元素的2次方
print('计算出来的结果是exp(x)-1,稀疏矩阵中对应元素的计算')
print(( csc_matrix.expm1(X)).toarray()) # 计算出来的结果是exp(x)-1,稀疏矩阵中对应元素的计算

result