用ODE23 Matlab解决牛顿冷却方法

Question

我是Matlab的新手，尝试使用DEQ来解决牛顿冷却ODE23，实际上是pathname expansion版本。我不知道如何将函数输入ODE函数。

我T_s = 19, T_c = 84 r = 0.025已经过[0, 300] seconds。{/ p>

等式为(dT_c)/dt = -r(T_c - T_s)

任何人都可以帮我使用链接功能来解决此DEQ吗？

Answer 1

Matlab支持higher order functions以及anonymous functions。您可以将f(T_c, t)传递给ode23tx 如下：

ode23tx(@(T_c, T_s)(-r*(T_c - T_s), tspan, y0, arg4, varargin)

我假设r是不变的，但你可以轻松地改变这一点。 阅读链接函数的documentation以获取另一个的详细信息 参数。