在 MATLAB 中求解耦合微分方程和代数方程

计算科学 matlab pde 非线性方程 线法 造型
2021-12-25 12:49:16

我想用直线法在 MATLAB 中求解一个由 7 个耦合微分方程和 1 个代数方程组成的系统。我可以为每个独立的方程做一些假设,但我不能一起解决这 8 个方程。你能帮助我吗?

1个回答

您所描述的称为微分代数方程(DAE)系统。根据系统的指标,这些问题可能很容易解决,也可能非常困难。看看: http: //www.mathworks.com/matlabcentral/fileexchange/7481-manuscript-of--solving-index-1-daes-in-matlab-and-simulink-

Matlab 的求解器ode15sode23t可以处理索引为 1 的 DAE。如果系统比较僵硬,ode15s 更合适。

如果您对它们背后的理论或什么是indexstiff感兴趣,请查看:http ://www.scholarpedia.org/article/Differential-algebraic_equations

编辑: 我没有意识到 OP 是专门询问线条的方法。这里有一个关于 DAE 的部分:http ://www.scholarpedia.org/article/Method_of_lines