2

我正在尝试在 Python 中使用 for 循环生成一个简单的 Toeplitz 矩阵

a = np.array([[3,4,6]])
b = np.zeros((np.size(a),np.size(a)))

b[0,:] = a
b[:,0] = a

for i in range(1,np.size(a)):
    for j in range(1,np.size(a)):
        a[i,j] = a[i-1,j-1]

应该工作，除非我遗漏了什么，但给出了这个错误：

    a[i,j] = a[i-1,j-1]

IndexError: index 1 is out of bounds for axis 0 with size 1

如何在没有 scipy 的内置 toeplitz() 函数的情况下使它成为 Toeplitz 你能帮忙吗？

4

1 回答 1

0

我相信你的意思是

    b[i,j] = b[i-1,j-1]

不是 a，它显然只有一个维度。

于 2021-03-17T03:30:13.357 回答