Scipy odeint 没有参数或括号的方法

Question

请帮助 - 我无法理解我自己的代码！哈哈 我是 python 的新手，经过多次试验和错误，我的代码可以正常工作，但其中有一个特定部分我不明白。

在下面的代码中，我通过 scipy 的 odeint 函数解决了一个相当基本的 ODE。我的目标是在此蓝图的基础上构建更复杂的系统。

我的问题：如何在没有任何参数且没有右括号的情况下调用方法.reaction_rate_simple？这在 python 中是什么意思？我应该在某处使用静态方法吗？

如果有人对此有任何反馈 - 也许这是一段糟糕的代码，但有更好的解决方法！

非常感谢任何回复和帮助！

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint

class batch_reator:

    def __init__(self, C_init, volume_reactor, time_init, time_end):
        self.C_init = C_init
        self.volume = volume_reactor
        self.time_init = time_init
        self.time_end = time_end
        self.operation_time = (time_end - time_init)
    
    def reaction_rate_simple(self, concentration_t, t, stoch_factor, order, rate_constant):
        reaction_rate = stoch_factor * rate_constant * (concentration_t ** order)
        return reaction_rate

    def equations_system(self, kinetics):
        dCdt = kinetics
        return dCdt


C_init = 200
time_init, time_end = 0, 1000 
rate_constant, volume_reactor, order, stoch_factor = 0.0001, 10, 1, -1
time_span = np.linspace(time_init, time_end, 100)
 
Batch_basic = batch_reator(C_init, volume_reactor, time_init, time_end)       
kinetics = Batch_basic.reaction_rate_simple

sol = odeint(Batch_basic.equations_system(kinetics), Batch_basic.C_init, time_span, args=(stoch_factor, order, rate_constant))

plt.plot(time_span, sol)
plt.show() 

Answer 1

我假设你指的是这条线

kinetics = Batch_basic.reaction_rate_simple

您不是在调用它，而是将方法保存为变量，然后将该方法传递给equations_system(...)，它只是返回它。我不熟悉 odeint，但根据文档，它接受一个可调用对象，这就是你给它的。

在python函数中，lambdas、类都是可调用的，可以赋值给变量，传递给函数，根据需要调用。

在这种特殊情况下，来自odeint docs 的回调定义说

func callable(y, t, ...) 或 callable(t, y, ...) 计算 y 在 t 处的导数。如果签名是可调用的(t, y, ...)，则参数 tfirst 必须设置为 True。

所以前两个参数是由 odeint 传入的，另外三个来自你指定的参数。