1

我正在尝试使用 numpy 矩阵列表,但遇到了一个烦人的问题。

假设我从一个包含 10 个 2x2 零矩阵的列表开始

para=[numpy.matrix(numpy.zeros((2,2)))]*(10)

我像这样访问单个矩阵

para[0]
para[1]

等等。到现在为止还挺好。

现在,我只想修改第二个矩阵的第一行,其余的保持不变。所以我这样做

para[1][0]=numpy.matrix([[1,1]])

第一个索引指向列表中的第二个矩阵,第二个索引指向该矩阵中的第一行,将其替换为 [1,1]。

但奇怪的是,这个命令将列表中所有十个矩阵的第一行更改为 [1,1] 而不是我想要的第二个。是什么赋予了?

4

1 回答 1

5

当您将初始列表乘以 10 时,您最终会得到一个包含 10 个 numpy 数组的列表,这些数组实际上是对相同底层结构的引用。修改一个将修改所有这些,因为实际上只有一个 numpy 数组,而不是 10 个。

如果您需要证明,请查看 REPL 中的此示例:

>>> a = numpy.zeros(10)
>>> a = [numpy.zeros(10)]*10
>>> a[0] is a[1]
True
>>> 

is运算符检查两个对象是否实际上相同(而不是它们的值是否相等)。

您应该做的是使用列表推导来生成初始数组而不是乘法,如下所示:

para=[numpy.matrix(numpy.zeros((2,2))) for i in range(10)]

这将调用numpy.matrix()十次而不是一次并生成 10 个不同的矩阵。

于 2013-02-23T22:37:09.063 回答