2
我喜歡這個轉身稀疏矩陣
>>>import numpy as np
>>>from scipy.sparse import *
>>>A = csr_matrix((np.identity(3)))
>>>print A
(0, 0) 1.0
(1, 1) 1.0
(2, 2) 1.0
爲了更好地理解A
一些稀疏矩陣是這樣的:
>>>print A.todense()
[[ 1. 0. 0.]
[ 0. 1. 0.]
[ 0. 0. 1.]]
,我想有一個運營商(我們稱之爲op1(n)
)這樣做:
>>>A.op1(1)
[[ 0. 1. 0.]
[ 0. 0. 1.]
[ 1. 0. 0.]]
=>使最後n
列第一n
的, 所以
>>>A == A.op1(3)
true
。是否有一些內置解決方案,(EDIT :),它又返回一個稀疏矩陣? 與roll
解決辦法:
X = np.roll(X.todense(),-tau, axis = 0)
print X.__class__
回報
<class 'numpy.matrixlib.defmatrix.matrix'>
作品完美!返回一個ndarray並且沒有矩陣?有不同的嗎? – 2012-07-29 19:38:03