您好我正在尝试使用中心有限差分法解决一组转换为 ODE 和 DAE 的 pde。我已经使用 MATLAB 的“求解”命令来确定边界条件下虚拟节点的系数。我正在使用 ode15s 求解器并得到“此 DAE 的索引大于 1”daeic12 (line76) 错误代码。我已将代码作为链接附在下面。要运行的脚本文件是 main_code_live，它是一个实时脚本，解释了 pde 是如何转换为 ode 和 dae 以及边界条件的，其余都是运行主脚本所必需的函数文件。我尝试使用 Mathworks 网站上的很多建议，但仍然没有成功。任何帮助，将不胜感激！https://drive.google.com/drive/folders/1emgOcdsnP81MP5ZBERefEJ0d4IP0qbSL?usp=sharing