请教用rk方法求解变系数微分方程组的matlab程序

shjjuanzi

我建立的方程是二自由度机构的动力学方程 他的系数矩阵跟运动状态有关 每个时刻的系数矩阵都得通过运动分析获得  能用matlab编程实现吗？运动分析已经作出来 就是不知道怎么用rk法实现一步一步的积分 请高手指点

不知道说的够清楚吗？
我的思路是进行一次运动学分析获得这个时刻的系数矩阵 这时候就把系数矩阵看成常数矩阵 用龙哥库塔方法积分求解出下个时刻的数值解 然后再进行运动学分析 以此循环下去
现在我就是不知道 用龙哥库塔方法 要怎么求出某个时刻的一个数值解

[ 本帖最后由 ChaChing 于 2009-8-24 22:41 编辑 ]

无水1324

我的问题是，你这个运动分析是怎么计算出来的？是一个代数方程吗？，你可以不求出来直接方程微分方程里面，就把它当着一个随时间变化的参数就可以了阿？
另外求出来我还觉得这个是不是随时间变化的呢？