samson888 发表于 2010-3-30 19:55

求助---求解状态方程中遇到的问题

小弟最近在求解状态方程时遇到问题,还请大虾们帮忙指点;
问题是:状态方程 dx/dt=Ax+Bu,其中dx/dt表示x对时间t的一阶导数,A是6X6阶矩阵,B是6X3阶矩阵,U是输入变量(常数)但是矩阵A中包含有状态变量x,也就是说方程中的系数A是状态变量的函数,我考虑用4阶龙格库塔法(自编程序)解状态方程,采用定步长,因此,在计算中,每循环一次,方程系数必须重新更新计算一次,把新的状态变量赋到A中,所以必须把与状态变量有关的系数放在循环内,但是如何编写这个循环子程序,使得每次更新的系数释放出来重新回到循环中去呢???请高手指导,万分感谢
页: [1]
查看完整版本: 求助---求解状态方程中遇到的问题