zigongzhb 发表于 2008-5-26 12:33

关于ode子函数编写

我的振动程为 MX''+KX=F.用0de45求解，在编写子函数时，有没有什么简单的方法？也就是直接利用质量矩阵和刚度矩阵。一个一个方程的写太麻烦了。

gh688 发表于 2008-5-26 22:59

回复 楼主 的帖子

以前求一阶微分方程X‘=AX时候直接定义了一个子函数myfun，这个函数包含了矩阵A
function xdot=myfun(t,x,A)
xdot=A*x
二阶的微分方程你变一下，原理也是相同的
页: [1]
查看完整版本: 关于ode子函数编写