Question

我想从一维数组构造一个(n,n) - 数组，其中每一行相对于前一行被移位（包裹）一个。以下代码执行此操作：

import numpy as np

r = np.array([1, 2, 3, 4, 5])
n = len(r)
MM = np.zeros((n, n), dtype=r.dtype)
for k in range(n):
    MM[k, :] = np.roll(r, k)

print(MM)

导致：

[[1 2 3 4 5]
 [5 1 2 3 4]
 [4 5 1 2 3]
 [3 4 5 1 2]
 [2 3 4 5 1]]

对于Numpy中的大for，有没有办法更快地执行Numpy，即避免r - 循环？

Answer 1

查看scipy.linalg.circulant

In [255]: r
Out[255]: array([1, 2, 3, 4, 5])

In [256]: circulant(r).T
Out[256]: 
array([[1, 2, 3, 4, 5],
       [5, 1, 2, 3, 4],
       [4, 5, 1, 2, 3],
       [3, 4, 5, 1, 2],
       [2, 3, 4, 5, 1]])

或scipy.linalg.toeplitz

In [257]: toeplitz(np.roll(r[::-1], 1), r)
Out[257]: 
array([[1, 2, 3, 4, 5],
       [5, 1, 2, 3, 4],
       [4, 5, 1, 2, 3],
       [3, 4, 5, 1, 2],
       [2, 3, 4, 5, 1]])

在Numpy

1 个答案: