如何使用 Matlab 集成非平滑函数？可以使用现金卡普方法吗？Dormand-Prince 方法在集成非平滑解时似乎会出错。

我正在使用 MATLAB ODE45 来积分我的常微分方程。

$\{q\}$是长度为 14 的列向量。右侧的激振力以下列形式变化

在哪里$\{h_p\}$,$\{g_p\}$,$\{h_g\}$， 和$\{g_g\}$是长度为 6 的列向量

我称它为不光滑是因为$\{F_\text{excit}\}$随$F_\text{mf}$.

方程组可以在以下情况下求解$F_\text{mf}=0$. 什么时候$F_\text{mf}$是非零，我收到一个警告，实际上表明无法执行积分。

Warning/Error Message in Matlab

Warning: Failure at t=1.821477e+00.  Unable to meet
integration tolerances without reducing the step 
size below the smallest value
allowed (3.552714e-15) at time t.

我确实尝试过 ODE15S，但错误出现的时间比 ODE45 早得多。两者都给了我同样的错误

我用方括号来表示矩阵，用花括号来表示向量。