用矩阵条目制作直方图？

Question

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

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

现在我的尝试是创建一个空列表，然后使用 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。