如果您可以将其编写为质量矩阵形式，则 Julia 中的 DifferentialEquations.jl 可以做到这一点。您不会在教程中找到它，但您可以提供质量矩阵作为SDEProblem. 一些僵硬的求解器可以解决这个问题（我看到它还没有很好地记录哪些，但它是辛和隐含的欧拉形式）。

我会提醒这是相当新的，所以我会在非 DAE 上测试它，即使用质量矩阵功能，然后添加$M^{-1}$到您的 SDE 方程并检查以确保它们计算相同的东西。如果您发现任何差异，请提出问题。此处隐式处理的架构基于经过良好测试的 ODE 求解器，但由于我看到我们在持续集成套件中没有明确地进行质量矩阵 SDE 测试，但我只会在信任它之前进行此测试.

目前还没有适用于 SDAE 的自适应时间步长算法（仅适用于 SDE）……嗯……已发布。有一个用于质量矩阵形式的加性噪声​​ SDAE 的自适应时间步长刚性求解器，它正在生产中，但要等到出版后才会发布。希望几个月后我可以编辑这篇文章，使其更加可靠，但是现在我们所拥有的至少应该可以解决问题。