| 查看: 1982 | 回复: 13 | |||
Jason_Cindy铁虫 (小有名气)
|
[求助]
四阶龙格库塔求解微分方程组求救 已有1人参与
|
||
 如图所示系数矩阵均是关于Vm Vgm及P和α表达式，如何利用四阶龙格库塔法对上述微分方程组进行求解，望大家不吝赐教，不胜感激！ 先说下我目前的想法：通过线性方程组的矩阵运算可以求出dα/dz的表达式，进而利用龙格库塔方法进行求解（不进行迭代，只算一次），可得α1,然后利用（α1-α）/dz即可得到这步下的dα/dz，让后回到线性方程组的求解中，得到dp/dz的表达式，以此类推进行求解，不知道这样是否合理。正确的方法又该如何求解？  微分方程组.png |
回复此楼
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3881
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
3楼2015-10-22 10:36:46
Jason_Cindy
铁虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 195.1
- 散金: 151
- 红花: 1
- 帖子: 102
- 在线: 23.8小时
- 虫号: 2977176
- 注册: 2014-02-19
- 专业: 多相流热物理学
2楼2015-10-22 08:39:21
Jason_Cindy
铁虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 195.1
- 散金: 151
- 红花: 1
- 帖子: 102
- 在线: 23.8小时
- 虫号: 2977176
- 注册: 2014-02-19
- 专业: 多相流热物理学
4楼2015-10-22 10:38:44
Jason_Cindy
铁虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 195.1
- 散金: 151
- 红花: 1
- 帖子: 102
- 在线: 23.8小时
- 虫号: 2977176
- 注册: 2014-02-19
- 专业: 多相流热物理学
5楼2015-10-23 14:35:27
破晓林海雪原
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 2045.1
- 散金: 122
- 红花: 4
- 帖子: 236
- 在线: 159.8小时
- 虫号: 4048270
- 注册: 2015-09-02
- 专业: 计算机软件
|
我觉得你这种是正确的，求解的时候可能要都求出来，因为下一次迭代会用到相关数据。。做过一个二阶的，迭代两千次左右，结果不好，猜测是误差累积以及初始值选择不当的问题，因为前面的结果可以，后来的才逐渐变坏。。 发自小木虫Android客户端 |
6楼2015-10-23 16:16:08
破晓林海雪原
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 2045.1
- 散金: 122
- 红花: 4
- 帖子: 236
- 在线: 159.8小时
- 虫号: 4048270
- 注册: 2015-09-02
- 专业: 计算机软件
7楼2015-10-23 16:17:22
Jason_Cindy
铁虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 195.1
- 散金: 151
- 红花: 1
- 帖子: 102
- 在线: 23.8小时
- 虫号: 2977176
- 注册: 2014-02-19
- 专业: 多相流热物理学
8楼2015-10-23 17:48:49
破晓林海雪原
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 2045.1
- 散金: 122
- 红花: 4
- 帖子: 236
- 在线: 159.8小时
- 虫号: 4048270
- 注册: 2015-09-02
- 专业: 计算机软件
9楼2015-10-23 19:57:42
破晓林海雪原
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 2045.1
- 散金: 122
- 红花: 4
- 帖子: 236
- 在线: 159.8小时
- 虫号: 4048270
- 注册: 2015-09-02
- 专业: 计算机软件
10楼2015-10-23 20:14:33
