在 Python 上处理错误时遇到问题

Question

好吧，经过一番研究，我发现我不知道我的代码有什么问题。这是我正在研究的牛顿法优化代码，我想在另一个程序的 for 循环中的几个函数上测试它，以及其他几种优化方法，但是某个函数出现了问题不要让代码继续测试一切。

问题是牛顿法要求矩阵不是奇异的才能解决系统问题，并且在这个有问题的函数上，该矩阵最终会在某些时候成为奇异矩阵，我收到一条错误消息告诉我并且停止整个过程。到目前为止，我对 try 和除了错误处理一无所知，所以经过一些研究，我弄清楚了它并在我的代码上实现它以尝试避免这个错误，但是，它仍然会因为错误而停止整个程序，会有人帮我看看吗？

这是问题部分的代码：

while np.dot(grad(x),grad(x)) > 10**(-4):
        if f(x+a*p) <= f(x)+c1*a*np.dot(grad(x),p):
            x = x+a*p
            try:
                p = -np.linalg.solve(hess(x),grad(x))
            except:
                break
            k = k+1
        else:
            a = 0.9*a
        if k >= 100:
            break
return [x, f(x),grad(x)]

hess(x) 在某个时候结束为奇异矩阵并给我一个异常错误，但代码不会尝试不执行 try 块并在这种情况下执行 except 块，因为我认为它应该，什么正在进行吗？

我尝试过一个较小的问题，例如

i = 0
while i < 10:
    if 2 < 3:
        try:
            -np.linalg.solve(A,b)
        except:
            print(9)
            break
    i = i+1

A 是一个奇异矩阵，它工作得很好，只需打印一个“9”，那么为什么在主程序上这根本不会发生呢？也许这与我在另一个程序上调用此函数或类似的事实有关？

我可以尝试通过 if 语句预先测试矩阵是否为单数并在是单数时中断循环来解决这个问题，但我觉得这对于迭代来说成本很高，所以我想避免它为了以后更好地比较这些方法，所以我想尽可能地避免这种情况。

这是我执行代码时得到的结果：

Traceback (most recent call last):

  File "<ipython-input-1-31bd7312c08f>", line 1, in <module>
    runfile('C:/Users/CLIENTE/Desktop/Coisas/Estudos 2.0/Ñ Linear/Testes.py', wdir='C:/Users/CLIENTE/Desktop/Coisas/Estudos 2.0/Ñ Linear')

  File "C:\Users\CLIENTE\Anaconda3\lib\site-packages\spyder_kernels\customize\spydercustomize.py", line 704, in runfile
    execfile(filename, namespace)

  File "C:\Users\CLIENTE\Anaconda3\lib\site-packages\spyder_kernels\customize\spydercustomize.py", line 108, in execfile
    exec(compile(f.read(), filename, 'exec'), namespace)

  File "C:/Users/CLIENTE/Desktop/Coisas/Estudos 2.0/Ñ Linear/Testes.py", line 52, in <module>
    effer = método(dado)

  File "C:\Users\CLIENTE\Desktop\Coisas\Estudos 2.0\Ñ Linear\Método_de_Newton.py", line 27, in MétodoDeNewton
    p = -np.linalg.solve(hess(x),grad(x))

  File "C:\Users\CLIENTE\Anaconda3\lib\site-packages\numpy\linalg\linalg.py", line 394, in solve
    r = gufunc(a, b, signature=signature, extobj=extobj)

  File "C:\Users\CLIENTE\Anaconda3\lib\site-packages\numpy\linalg\linalg.py", line 89, in _raise_linalgerror_singular
    raise LinAlgError("Singular matrix")

LinAlgError: Singular matrix

Answer 1

您是否尝试过捕获异常的完全限定名称？

while np.dot(grad(x),grad(x)) > 10**(-4):
        if f(x+a*p) <= f(x)+c1*a*np.dot(grad(x),p):
            x = x+a*p
            try:
                p = -np.linalg.solve(hess(x),grad(x))
            # HEY! Look here.
            except np.linalg.LinAlgError:
                break
            k = k+1
        else:
            a = 0.9*a
        if k >= 100:
            break
return [x, f(x),grad(x)]