Question

今天我的任务是制作一个直方图来表示A ^ n的运算，其中A是一个矩阵，但仅适用于矩阵中的特定条目。

例如，假设我有一个矩阵，其中行总和为1。第一个条目是一些特定的十进制数。但是，如果我将该矩阵提升到第二个幂，那么第一个条目就会变成别的东西，如果我将该矩阵提升到第三个幂，它会再次变化，等等 - 这就是我需要绘制的内容。

现在我的尝试是创建一个空列表，然后使用for循环将矩阵乘法产生的条目添加到列表中。但是，它所做的只是将最终矩阵乘法的结果打印到列表中，而不是在每次迭代时打印它的值。

以下是我所说的具体代码：

print("The intial probability matrix.")

print(tabulate(matrix))

baseprob = []

for i in range(1000):

     matrix_n = numpy.linalg.matrix_power(matrix, s)

     baseprob.append(matrix_n.item(0))

print(baseprob)

print("The final probability matrix.")

print(tabulate(matrix_n))

这是完整的代码，以及我得到的输出。

http://pastebin.com/EkfQX2Hu

Answer 1

当然它只打印最终值，你正在进行相同的操作，matrix^s，1000次。您需要s每次更改1000次。

如果您想计算matrix(0)位置matrix^i的所有值，其中i是从1到s的每个值（您的最终权力）：

baseprob = []

for i in range(1,s): #changed to do a range 1-s instead of 1000

     #must use the loop variable here, not s (s is always the same)
     matrix_n = numpy.linalg.matrix_power(matrix, i)

     baseprob.append(matrix_n.item(0))

然后baseprob将matrix(0)，matrix^1等等matrix^2保留matrix^s。

用矩阵条目制作直方图？

1 个答案: