guojing2345 发表于 2012-10-23 21:57

关于龙格库塔法的问题，高手指点

请问这个方程可以拿龙格库塔法求见？

kezairenjian 发表于 2012-10-24 08:24

这个用精细积分法去计算吧

guojing2345 发表于 2012-10-24 11:04

高手可以赐教一下吗

guojing2345 发表于 2012-10-24 11:06

由于矩阵R不可逆，所以无法求得他的逆！所以感觉无法用龙格求解？

gghhjj 发表于 2012-10-24 16:57

guojing2345 发表于 2012-10-24 11:06 static/image/common/back.gif
由于矩阵R不可逆，所以无法求得他的逆！所以感觉无法用龙格求解？

一般这类方程用差分法求解比较多

hsfy919 发表于 2012-10-26 15:34

是矩阵奇异导致不可逆吧，matlab有现成的函数求解微分代数型方程。

guojing2345 发表于 2012-10-28 16:36

可否赐教一下！具体一点

yanyongju 发表于 2012-11-7 14:59

R奇异，显然是微风代数方程，精细积分可以计算，可能花费时间比较多，
Matlab的ODE15s可以解，需要设置刚度矩阵M
建议使用DAE求解方法，比如向后差分BDF，基于时间的中心差分离散等可以稍快一点。

guojing2345 发表于 2012-11-11 17:19

变系数矩阵怎么调用ODE15s，调用的时候矩阵随着时间是变化的

guojing2345 发表于 2012-11-11 17:20

特别感谢楼上！

guojing2345 发表于 2012-11-11 17:22

现在方程有点变化！真是感激各位C:\Users\asus\Desktop

guojing2345 发表于 2012-11-11 17:23

矩阵Q随时间变化

gghhjj 发表于 2012-11-14 16:53

guojing2345 发表于 2012-11-11 17:23 static/image/common/back.gif
矩阵Q随时间变化

基本思路应该不会有大的变化
页: [1]
查看完整版本: 关于龙格库塔法的问题，高手指点