象征性地计算插值函数的输入？

Question

我有一个相当复杂的函数H(x)，我正在尝试求解x 的值，使得H(x) = constant。我想使用从离散间隔生成的插值对象和 H(interval) 的相应输出来执行此操作，其中其他输入保持不变。我表示插值对象f。

我的问题是插值对象的调用函数接受一个array_like，所以传递一个符号给f(x)来使用sage的求解器方法是不可能的。关于如何解决这个问题的任何想法？

我有插值函数f。我想解方程f(x) == sageconstant forx.

    from scipy.interpolate import InterpolatedUnivariateSpline as IUspline
    import numpy as np

    #Generating my interpolation object
    xint = srange(30,200,step=.1)
    val = [H(i,1,.1,0,.2,.005,40) for i in srange(30,299,step=.1)]
    f = IUspline(xint,val,k=4)

    #This will yield a sage constant
    eq_G(x) = freeB - x 

    #relation that I would like to solve
    eq_m(x) = eq_G(39.9) == f(x)
    m = solve(eq_m(x),x)

上面的代码（更具体的f(x)）生成

"TypeError: 无法将数组数据从 dtype('0') 转换为 dtype('float64') 根据“安全”规则。

编辑：任何函数H(x) 都会导致相同的错误，因此H(x) 是什么并不重要。为简单起见（当我说H 很复杂时，我不是在开玩笑），试试H(x) = x。然后该块将读取：

    from scipy.interpolate import InterpolatedUnivariateSpline as IUspline
    import numpy as np

    #Generating my interpolation object
    xint = srange(30,200,step=.1)
    H(x) = x
    val = [H(i) for i in srange(30,299,step=.1)]
    f = IUspline(xint,val,k=4)

    #This will yield a sage constant
    eq_G(x) = freeB - x 

    #relation that I would like to solve
    eq_m(x) = eq_G(39.9) == f(x)
    m = solve(eq_m(x),x)

Answer 1

使用 numpy 和 scipy 时，更喜欢 Python 类型而不是 Sage 类型。

使用 Python 整数和浮点数代替 Sage 整数和实数。

也许你可以像这样修复你的代码。

from scipy.interpolate import InterpolatedUnivariateSpline as IUspline
import numpy as np

# Generate interpolation object
xint = srange(30,200,step=.1)
xint = [float(x) for x in xint]
val = [float(H(i,1,.1,0,.2,.005,40)) for i in srange(30,299,step=.1)]
f = IUspline(xint,val,k=4)

# This will yield a Sage constant
eq_G(x) = freeB - x 

# relation that I would like to solve
eq_m(x) = eq_G(39.9) == f(x)
m = solve(eq_m(x),x)