这是我的示例代码。这是一个迭代过程的高斯赛德尔(矩阵解算器)。本质上,当误差足够小时,它就会脱离while循环。i=1
while (i>0):
x_past = x_present
j=0
while(j<3):
value=0
k=0
while(k<3):
if(k!=j):
if(i==1):
if(k>j):
value=value+0
else:
value=value+x_present[k]*eqn[j][k]
else:
value=value+x_present[k]*eqn[j][k]
else:
value=value+eqn[j][k]
k=k+1
x_present[j:j+1]=[value]
j=j+1
print "X_PAST"
print x_past
print "X_PRESENT"
print x_present
if(error(x_past, x_present)<10**-2):
break;
i=i+1
我简化了代码,使其更易于管理。如果你不明白它在做什么对解决这个问题来说并不重要。
这就是问题所在。每次x_present[j:j+1]=[value]
是运行的,x_过去等于x_现在。我不知道为什么会这样,因为我把x戋over设置为x戋present的唯一地方是在循环的顶部。如果我拿走x_past=x_present
一句话,过去永远不等于现在。这使我认为是这两种说法的某种结合导致了这个问题。
这是一个大问题,因为如果每次x_past=x_出现错误=0,并且循环在第一次迭代后中断。代码确实有效,例如,如果我告诉代码运行8次迭代,然后中断它,它就会给出它应该给出的答案。
在过去的4个小时里,我一直在想办法解决这个问题,但我完全被难住了。我使用python的时间不长,所以我在语法方面的故障排除技能并不是很好。任何帮助都将不胜感激!!