飞天小男警 发表于 2013-3-13 10:24

如何用龙格库塔法解耦合的二阶微分方程组

    大家好,问题是一个方程组,mx''+cx'+kx=0,因为m为非对角矩阵,就是质量耦合,用龙格库塔法要怎么算?具体情况如下:function dy=fun(t,y)
dy(1)=y(2);
dy(3)=y(4);
dy(5)=y(6);
dy(2)=。。。。。。-e*dy(4)-i*dy(6);
dy(4)=。。。。。。。。
dy(6)=。。。。。。。。
t(end,1)
dy=
省略掉了没用的东西,在算dy(2)的时候出现了dy(4)和dy(6)由于没有初值[初值赋的y(1)到y(6)]算不下去,真心求教该怎么算,十分感谢!!

飞天小男警 发表于 2013-3-17 00:03

解决了,谢谢大家

godlinux 发表于 2013-4-23 22:20

飞天小男警 发表于 2013-3-17 00:03 static/image/common/back.gif
解决了,谢谢大家

楼主你是怎么解决的?可以说一下吗?万分感激

mxlzhenzhu 发表于 2013-7-12 23:38

楼主解决了么,对于质量耦合的振动方程组用什么算法最高效呢?

z13464238735 发表于 2013-8-6 10:07

飞天小男警 发表于 2013-3-17 00:03 static/image/common/back.gif
解决了,谢谢大家

朋友,我也遇到你这样的问题,请问怎么解决的?

TWO-SOCKS 发表于 2014-6-16 19:16

楼主你好!我最近也在解这个方程,困扰中 请问你是怎么解的?谢谢!

llffox 发表于 2014-6-30 15:59

我也想知道怎么解决的,楼主能指教一下么
页: [1]
查看完整版本: 如何用龙格库塔法解耦合的二阶微分方程组