Python fsolve 对象回火

Question

我用 Python 编写了一个脚本，它使用 fsolve 找到一个相当复杂的函数的零。它的工作方式如下。有一个类只是简单地存储函数的参数。该类有一个 evaluate 方法，它根据存储的参数返回一个值，另一个方法(反转)找到函数获取提供的输出的参数。

反转方法会在每次迭代时更新函数的参数，并一直这样做，直到评估方法返回的值与提供的值之间的不匹配为零。<​​/p>

我遇到的问题是，虽然反转方法返回的值是正确的，但作为对象一部分的参数在反转方法终止后始终为 0。奇怪的是，如果我使用 root 而不是 fsolve，这个问题就会消失。据我所知，fsolve 只是 root 的包装器，其中包含求解器算法的一些设置和其他一些强制执行的设置。

这是 fsolve 的已知问题还是我在这里做了一些愚蠢的事情？下面的脚本演示了我在正弦函数上遇到的问题。

from scipy.optimize import fsolve, root
from math import sin, pi

class invertSin(object):
    def __init__(self,x):
        self.x = x

    def evaluate(self):
        return sin(self.x)

    def arcsin_fsolve(self,y):

        def errorfunc(xnew):
            self.x = xnew
            return self.evaluate() - y

        soln = fsolve(errorfunc, 0.1)
        return soln

    def arcsin_root(self,y):

        def errorfunc(xnew):
            self.x = xnew
            return self.evaluate() - y

        soln = root(errorfunc, 0.1, method = 'anderson')
        return soln

myobject = invertSin(pi/2)
x0 = myobject.arcsin_fsolve(0.5) #find x s.t. sin(x) = 0.5 using fsolve
print(x0) #this prints pi/6
x0obj = myobject.x
print(x0obj) #this always prints 0 no matter which function I invert

myobject2 = invertSin(pi/2)
x1 = myobject2.arcsin_root(0.5) #find x s.t. sin(x) = 0.5 using root
print(x1) #this prints pi/6
x1obj = myobject2.x
print(x1obj) #this prints pi/6

Answer 1

如果您在 errorfunc 中为 xnew 添加打印语句，那么您将看到 fsolve 使用列表(一个元素)。这意味着以这种方式重新解释函数，而不是原始函数。不知何故，类型信息在退出求解器后丢失，因此对该列表的地址/引用被解释为 float 据，从而给出错误的值。

在那里设置 self.x = xnew[0] 可以恢复所需的行为。