24小时热门版块排行榜    

查看: 752  |  回复: 7
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

eric-sun

新虫 (初入文坛)

[交流] 【求助】请教:matlab求解微分方程组

微分方程组:
dA1/dz=f(A2,A3,z),
dA2/dz=g(A1,A3,z),
dA3/dz=h(A1,A2,z),

初始条件:A1(z=0)、A2(z=0)、A3(z=0)是三个已知三维常矩阵

如何用四阶龙格库图方法求解? 如果用ode函数求解,可否将初始条件三个矩阵作为一个整体代入
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

eric-sun

新虫 (初入文坛)

引用回帖:
Originally posted by 边际浪子 at 2009-10-24 18:30:
是因为耦合吗?如果是,可以采取迭代算法,具体的先用上一层的值求A1,然后用上一层的啊!A2、A3和这一层的A1求A2,再用上一层的A3和这一层的A1、A2求A3,反复进行多次,直到两次的误差比较小,进入下一层循环。

请问“两次误差比较小”是什么意思?
7楼2009-11-09 18:17:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

边际浪子

银虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
nono2009(金币+0,VIP+0):这好像不是matlab的风格。 10-24 20:33
是因为耦合吗?如果是,可以采取迭代算法,具体的先用上一层的值求A1,然后用上一层的啊!A2、A3和这一层的A1求A2,再用上一层的A3和这一层的A1、A2求A3,反复进行多次,直到两次的误差比较小,进入下一层循环。
3楼2009-10-24 18:30:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

边际浪子

银虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
nono2009(金币+1,VIP+0):谢谢交流!欢迎常来。 10-25 15:55
Matlab也可以写程序的:)
外面加一层循环即可!
4楼2009-10-25 15:20:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

eric-sun

新虫 (初入文坛)

引用回帖:
Originally posted by 边际浪子 at 2009-10-24 18:30:
是因为耦合吗?如果是,可以采取迭代算法,具体的先用上一层的值求A1,然后用上一层的啊!A2、A3和这一层的A1求A2,再用上一层的A3和这一层的A1、A2求A3,反复进行多次,直到两次的误差比较小,进入下一层循环。

1、学友说的是算法,四阶RK也是算法,差别在于迭代斜率的确定公式;

2、我想确定一下,有没有不用自己写算法的、现有的函数求解耦合微分方程组(ODE类函数是否可以用三维矩阵作为初始条件)?

3、如果是自己写算法,如何估计计算误差(误差算法表达式)?

谢谢回帖!
5楼2009-10-26 09:22:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 14:05 by GkCVyODlx4mR
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 13:54 by GkCVyODlx4mR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:41 by l0VvVHGBGRLv
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +8 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:40 by l0VvVHGBGRLv
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:18 by l0VvVHGBGRLv
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:07 by l0VvVHGBGRLv
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:39 by ZPa0EcMwuECS
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:28 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:17 by ZPa0EcMwuECS
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 06:55 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见