使用numpy/scipy的Toeplitz矩阵

Question

在Octave或Matlab中，有一种简洁、简洁的方法来创建大的Toeplitz矩阵，例如：

T = toeplitz([1,-0.25,zeros(1,20)])

这节省了大量的时间，否则，通过使用额外的代码行来填充矩阵中的几十个或几百个零。在

然而，我似乎不能对scipy或numpy执行相同的操作，尽管这两个库都有toeplitz（）和zeros（）函数。 有没有类似的方法可以做到这一点呢？还是我必须自己制定一个程序来完成（不是一个大问题，但仍然是一个麻烦）？在

谢谢

F

Answer 1

目前我认为你能做的最好的是：

from numpy import concatenate, zeros
from scipy.linalg import toeplitz

toeplitz(concatenate([[1., -.25], zeros(20)]))

但是，对于python 3.5，我们将拥有：

^{pr2}$

所以这是值得期待的。在