插值函数与自定义函数乘积的数值积分报错及解决求助

你遇到的TypeError: unsupported operand type(s) for *: 'function' and 'FunctionClass'是因为你直接尝试把两个函数对象相乘了——Python里函数本身不能用*运算符做乘法操作，你需要定义一个新的函数，把两个函数的计算结果相乘后再传入积分函数。

具体解决方案步骤： #

1. 定义组合积分函数
   首先要确保你的f1是可以接收积分变量ktsq的函数（如果f1依赖其他非积分变量，可以通过闭包或者额外参数传递）。比如：

   • 如果f1没有额外参数：

     def integrand(ktsq):
         # 先分别计算两个函数的结果，再相乘
         return alphasinterp(ktsq) * f1(ktsq)
     

   • 如果f1有非积分变量（比如常数a、b）：

     def integrand(ktsq, a, b):
         return alphasinterp(ktsq) * f1(ktsq, a, b)
     

2. 调用积分函数并传入组合函数
   使用sciquad（即scipy.integrate.quadrature）时，它支持通过args参数传递额外的非积分变量：

   • 无额外参数的情况：

     result = sciquad(integrand, 0, 21400)
     

   • 有额外参数的情况（比如a=2、b=3）：

     result = sciquad(integrand, 0, 21400, args=(2, 3))
     

完整示例代码 #

结合你原有的插值代码，这里给出一个可运行的示例（假设f1是一个简单的平方根函数）：

import numpy
from numpy import loadtxt
from scipy.interpolate import *
from scipy.integrate import quadrature as sciquad

# 原有插值逻辑
x_given = numpy.arange(1,21400,21400/25000)
y_given = loadtxt("Yvalues.txt", comments="#", delimiter=",", unpack=False)
func_spline = interp1d(x_given,y_given, kind='cubic')
def alphasinterp(ktsq): return func_spline(ktsq)

# 自定义任意函数f1（替换成你实际的函数逻辑）
def f1(ktsq):
    return numpy.sqrt(ktsq)

# 定义组合积分函数
def integrand(ktsq):
    return alphasinterp(ktsq) * f1(ktsq)

# 执行积分
result = sciquad(integrand, 0, 21400)
print(result)

为什么原来的写法不行？ #

alphasinterp*f1这个表达式是在尝试把两个函数对象本身相乘，而不是让它们先计算出对应ktsq的结果再相乘。Python的函数对象不支持这种乘法操作，所以必须通过包装函数的方式，在函数内部完成结果的相乘逻辑。

内容的提问来源于stack exchange，提问作者CAF