Python稀疏矩阵计算中while循环更新值始终返回相同结果问题求助

问题排查与解决方案 #

我们先对齐你代码的逻辑预期：你使用的scipy.sparse.csr_matrix的*运算符默认对应矩阵乘法，符合你注释中r(t+1) = HS * r(t)的计算规则，手动计算前3轮的预期输出应该为：

• t=1第一轮计算后r值：[[1.], [0.5]]
• t=2第二轮计算后r值：[[1.5], [1.]]
• t=3第三轮计算后r值：[[2.5], [1.5]]

如果运行时每轮输出完全相同，可按照以下优先级排查错误： #

常见错误点修正 #

• 检查r变量的初始化位置
  如果你实际运行的代码把r = initial_PS写在了while循环内部，会导致每轮计算都用初始值重算，输出自然完全一致。你当前贴的代码里这行在循环外，确认实际代码没有放错位置。
• 检查乘法方法是否误用
  如果你实际代码用了逐元素乘法HS.multiply(r)替代*运算符，第一轮计算后结果就会固定，后续所有迭代输出都不会变化。确认你用的是矩阵乘法逻辑。
• 替换运算符规避版本兼容问题
  部分旧版本scipy的csr_matrix运算符存在逻辑bug，可以用官方推荐的dot方法替代*运算符修改代码，稳定性更高。

修正后可运行参考代码 #

import scipy.sparse as sparse

HS = [[1,1],[1,0]]
HS = sparse.csr_matrix(HS)
initial_PS = [[0.5],[0.5]]
initial_PS = sparse.csr_matrix(initial_PS)

t = 1
r = initial_PS

while t < 10:
  # 用dot方法执行矩阵乘法
  r = HS.dot(r)
  # 转为稠密矩阵打印可读性更高
  print(f"t={t}时r的值：\n{r.todense()}")
  t += 1

内容的提问来源于stack exchange，提问作者LLZ