【发布时间】:2017-09-20 08:09:04
【问题描述】:
我有一个代码如下:
def xgauss(self):
m, n = self.m, self.n
M = self.copy()
for k in range(n):
for i in range(k + 1, m):
yield M
if not likezero(M[i][k]):
lam = M[i][k] / M[k][k]
M[i] = M[i] - lam * M[k]
yield M
现在足以说如果值小于精度(用于数值比较),likezero 返回 true ;-)
我的目标是返回高斯消除算法的每次迭代以供学习(课堂)。
我有一些测试输入:
In [77]: A
Out[77]:
Matrix([
[ 1, 4, 1],
[ 1, 6, -1],
[ 2, -1, 2]
])
In [78]: Ab
Out[78]:
Matrix([
[ 1, 4, 1, 7],
[ 1, 6, -1, 13],
[ 2, -1, 2, 5]
])
奇怪与否,当我通过 Matrix Ab 调用生成器时(通过 A 我也得到相同的行为),我得到:
In [76]: list(Ab.xgauss())
Out[76]:
[Matrix([
[ 1, 4, 1, 7],
[ 0, 2, -2, 6],
[ 0, 0, -9, 18]
]), Matrix([
[ 1, 4, 1, 7],
[ 0, 2, -2, 6],
[ 0, 0, -9, 18]
]), Matrix([
[ 1, 4, 1, 7],
[ 0, 2, -2, 6],
[ 0, 0, -9, 18]
]), Matrix([
[ 1, 4, 1, 7],
[ 0, 2, -2, 6],
[ 0, 0, -9, 18]
])]
那是正确答案(应该只是最后一次迭代),但是我看不到每一步,生成器在所有迭代中返回矩阵结果。我不知道会发生什么。
【问题讨论】:
标签: python generator linear-algebra