24小时热门版块排行榜    

查看: 1741  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Jason_Cindy

铁虫 (小有名气)

[求助] 四阶龙格库塔求解微分方程组求救 已有1人参与


如图所示系数矩阵均是关于Vm Vgm及P和α表达式,如何利用四阶龙格库塔法对上述微分方程组进行求解,望大家不吝赐教,不胜感激!
先说下我目前的想法:通过线性方程组的矩阵运算可以求出dα/dz的表达式,进而利用龙格库塔方法进行求解(不进行迭代,只算一次),可得α1,然后利用(α1-α)/dz即可得到这步下的dα/dz,让后回到线性方程组的求解中,得到dp/dz的表达式,以此类推进行求解,不知道这样是否合理。正确的方法又该如何求解?

四阶龙格库塔求解微分方程组求救
微分方程组.png
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

引用回帖:
2楼: Originally posted by Jason_Cindy at 2015-10-22 08:39:21
求解答…

做拉普拉斯变换,一次求出;再做逆变换可以不?
3楼2015-10-22 10:36:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

Jason_Cindy

铁虫 (小有名气)

2楼2015-10-22 08:39:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jason_Cindy

铁虫 (小有名气)

引用回帖:
3楼: Originally posted by feixiaolin at 2015-10-22 10:36:46
做拉普拉斯变换,一次求出;再做逆变换可以不?...

用R-K法求解还需要拉氏变换吗?我不太懂怎么处理呀
5楼2015-10-23 14:35:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

破晓林海雪原

木虫 (小有名气)

我觉得你这种是正确的,求解的时候可能要都求出来,因为下一次迭代会用到相关数据。。做过一个二阶的,迭代两千次左右,结果不好,猜测是误差累积以及初始值选择不当的问题,因为前面的结果可以,后来的才逐渐变坏。。

发自小木虫Android客户端
6楼2015-10-23 16:16:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见