29、函数变化率与运动模拟的深入探究

函数变化率与运动模拟的深入探究

1. 随时间绘制体积变化

为了对函数进行柯里化，我们可以定义一个新函数，它接受前三个参数，并返回一个接受最后一个参数 T 的新函数：

def approximate_volume_function(q, v0, dt):
    def volume_function(T):
        return approximate_volume(q, v0, dt, T)
    return volume_function

这个函数直接从流量函数生成一个可绘制的体积函数。由于源代码中提供的体积函数 volume(0) 等于 2.3，我们将 v0 设为该值。最后，我们尝试 dt 值为 0.5，这意味着我们以半小时（30 分钟）为间隔计算体积变化。让我们看看它与原始体积函数的绘制对比：

plot_function(approximate_volume_function(flow_rate, 2.3, 0.5), 0, 10)
plot_function(volume, 0, 10)

好消息是，输出结果与原始体积函数非常接近！但 approximate_volume_function 生成的结果是锯齿状的，每 0.5 小时有一个台阶。你可能会猜测这与我们的 dt 值 0.5 有关，