使用 ode45 在 MATLAB 中求解具有离散值的微分方程

Question

我有一个微分方程-

L'(x) = F1(x,L(x))

使用 ode45，我得到了 L(x) 的解。我有一个 L(x) 值数组，用 L_val 表示。使用这个解决方案，我打算解决另一个微分方程。

w'(x)=L(x)/x

如何求解 w(x)？特别是因为 L(x) 不是 x 的函数，而是离散值的数组。

Answer 1

使用累积梯形积分函数： https://fr.mathworks.com/help/matlab/ref/cumtrapz.html

或者，您可以使用其他更准确的规则（Simpson 规则，在Matlab central 中搜索函数文件）。另一种实用的方法是计算函数的高阶插值并直接积分。