请教用rk方法求解变系数微分方程组的matlab程序
我建立的方程是二自由度机构的动力学方程 他的系数矩阵跟运动状态有关 每个时刻的系数矩阵都得通过运动分析获得能用matlab编程实现吗?运动分析已经作出来 就是不知道怎么用rk法实现一步一步的积分 请高手指点不知道说的够清楚吗?
我的思路是进行一次运动学分析获得这个时刻的系数矩阵 这时候就把系数矩阵看成常数矩阵 用龙哥库塔方法积分求解出下个时刻的数值解 然后再进行运动学分析 以此循环下去
现在我就是不知道 用龙哥库塔方法 要怎么求出某个时刻的一个数值解
[ 本帖最后由 ChaChing 于 2009-8-24 22:41 编辑 ]
回复 沙发 shjjuanzi 的帖子
我的问题是,你这个运动分析是怎么计算出来的?是一个代数方程吗?,你可以不求出来直接方程微分方程里面,就把它当着一个随时间变化的参数就可以了阿?另外求出来我还觉得这个是不是随时间变化的呢?
页:
[1]