Ka-*_*Yip 6 matlab matrix numerical-integration ode
在不使用质量矩阵的情况下,ode 求解器ode45可以求解 y'=f(t,y)。
但是对于涉及“质量”矩阵的问题,在 ode 求解器中有一个质量矩阵选项,M(t,y)y'=f(t,y)。
“质量”矩阵究竟是什么?这个术语是否来自质量弹簧系统振荡的质量?我在文档中找不到关于此的示例代码。此外,似乎我可以在 y'=f(t,y) 等式中的 f(t,y) 中对有关 t 和 y 的信息进行编码。在什么情况/示例中会出现 M(t,y)y'=f(t,y) 需要 M(t,y) 的地方?
文档中对此进行了完整解释odeset有完整解释。是的,它可以与二阶系统的质量/惯性项相关,但它也可以表示具有相同形式的其他系统中的不同参数。此选项可以在某些情况下提高效率并处理奇异(不可逆)质量矩阵(物理系统不常见)。
Shampine 等人的《Solving ODEs with Matlab》 (PDF)一书。提供了更多详细信息和一个很好的示例(第 2.3.2 节,第 105 页) \xe2\x80\x93 请参阅batonodeMatlab。