scipy.integrate.quad() 不使用 lambda

Question

我目前正在制作一个涉及整改的程序。唯一可以进行复杂集成的模块是带有scipy.integrate.quad() 命令的 Scipy。但是，在我的代码中，我需要有一个表示函数本身的变量，因为最终会在校正过程中使用更多不同的函数，这些函数是通过使用前面的等式得出的。我看到的关于该命令的每个来源要么涉及 lambda，要么创建一个输出方程的定义。您手动输入方程。所以，我的问题是，有没有办法在不这样做的情况下整合它？这是我的代码：

import scipy                                                   
import sympy
from scipy.integrate import quad

def rectify(smallestterm, nextsmallestterm, poly, coef, exp):
   test = rectification(coef,exp)
   a = quad(test, smallestterm, nextsmallestterm)
   print(a)
   a = a[0]
   dist = a/2
   return dist     
def test(x):
   return rectification(coef, exp)
def rectification(coef, exp):
   u = Poly(lint(coef,exp)) #Generated Equation
   b = u.all_coeffs()
   poly = b
   c = len(poly) - 1
   exponents = []
   while c + 1 > 0:
       exponents.append(c)
       c = c - 1
   poly = function(poly,exponents)   #Generated Equation in a form that can actually be used and identified as a function.
   return sqrt(1+(pow(diff(poly),2)))

其中coef 是多项式前导系数的列表，exp 是多项式前导指数的列表。本质上，它们都将间接组合在另一个定义中，function(coef, exp)（代码未显示）输出多项式（变量“poly”）。

即

function([2,4,5],[1,6,0]) 输出

4*x**6 + 2*x + 5

此代码（函数代码上方）不起作用，因为它不允许我使用变量“a”来表示整个函数，因为它只将“a”识别为函数本身。 因此，lambda 在我的情况下不起作用。我不能简单地做这样的事情：

import scipy
import sympy
from scipy.integrate import quad
poly = 2*x**5 + 5*x**4 - 4*x**2 + 10   #Some complicated polynomial
d = diff(poly) #differential of polynomial
a = sqrt(1+(pow(d,2)))
e = quad(a, 1, 5)
e = e[1]/2
return e

如果您需要查看我的完整代码以了解此代码中的任何其他功能，请提出要求，我将很乐意提供！

Answer 1

据我了解，您的代码使用 SymPy 生成符号表达式。这些不是 Python 意义上的 函数，因为它们不能被调用。所以它们不能直接用作quad的第一个参数。 SymPy 提供了lambdify 方法，将表达式包装成一个函数：

quad(lambdify(x, expr), 0, 1)

这里expr 可以是任何符号表达式，其中包含变量x，例如expr = 4*x**6 + 2*x + 5。