Python目标函数含变上限积分时的最小化优化报错解决

• 混用了sympy符号计算对象和scipy数值计算逻辑：代码里intM是sympy输出的符号表达式，直接传入仅接受数值输入的目标函数时，scipy优化器传入的浮点型t无法和符号表达式做运算，触发类型错误。
• sympy积分环节用法错误：原代码计算的是不定积分，本身带未知常数项，不符合「0到t积分」的定积分要求。
• 之前尝试quad时逻辑写错：积分上限是优化过程中每一步动态传入的t值，完全可以在目标函数内部实时计算，不需要提前固定积分限。
• 通配符导入from scipy import *会覆盖命名空间内的同名函数，额外触发不可预期的计算错误。

不需要依赖sympy，在目标函数内部每一步迭代时动态计算当前t对应的定积分即可，注意给t设置正的取值边界，避免t=0的除零错误：

import math
from scipy.optimize import minimize_scalar
from scipy.integrate import quad

# 固定参数
Cf = 100
Cp = 50
Eta = 72.66585511711865
Beta = 1.18609324

def R(t):
    return math.exp(-t/Eta)**Beta

def f(t):
    return (Beta/(Eta**Beta)) * (t**(Beta-1)) * (math.exp(-t/Eta)**Beta)

def M(t):
    return t * f(t)

def C(t):
    # 排除非正t的无效解
    if t <= 0:
        return 1e10
    # 动态计算0到当前t的定积分
    integral_M, _ = quad(M, 0, t)
    numerator = Cp * R(t) + Cf * (1 - R(t))
    denominator = t * R(t) + integral_M
    return numerator / denominator

# 采用有界优化，设置t的合理取值范围
res = minimize_scalar(C, bounds=(1e-6, 1000), method='bounded')
print(res)

如果要提前用sympy推导积分解析式，必须用sp.lambdify把符号表达式转成可接受数值输入的函数，且要计算0到t的定积分而非不定积分：

import math
import sympy as sp
from scipy.optimize import minimize_scalar

Cf = 100
Cp = 50
Eta = 72.66585511711865
Beta = 1.18609324

# 定义正实数符号变量，简化积分结果
t_sym = sp.Symbol('t', positive=True)
# 构造符号被积函数，计算0到t的定积分
M_sym = t_sym * (Beta/(Eta**Beta)) * t_sym**(Beta-1) * (sp.exp(-t_sym/Eta)**Beta)
intM_sym = sp.integrate(M_sym, (t_sym, 0, t_sym))
R_sym = sp.exp(-t_sym/Eta)**Beta

# 把符号表达式转成可输入数值的计算函数
intM_func = sp.lambdify(t_sym, intM_sym, 'math')
R_func = sp.lambdify(t_sym, R_sym, 'math')

def C(t):
    if t <= 0:
        return 1e10
    integral_M = intM_func(t)
    r_t = R_func(t)
    numerator = Cp * r_t + Cf * (1 - r_t)
    denominator = t * r_t + integral_M
    return numerator / denominator

res = minimize_scalar(C, bounds=(1e-6, 1000), method='bounded')
print(res)

两种方案运行后都能得到正确的最优解，最优t约为37.2，对应最小成本约为1.63。

• 禁止使用from xxx import *的通配符导入方式，避免命名空间污染引发的奇怪报错。
• 做带积分项的优化时，优先选数值积分动态计算的方案，避免符号积分无法得到闭式解、或转数值时出现兼容问题。
• 对物理意义为时间的变量t做优化时，必须设置正的取值边界，避免无效值进入计算触发错误。

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