飞天小男警 发表于 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]
查看完整版本: 如何用龙格库塔法解耦合的二阶微分方程组