求matlab编程求动力学方程的方法 - 计算模拟 - 其他 - 第2页小木虫论坛-学术科研互动平台

化学工程

11楼2016-02-19 15:37:55

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

引用回帖:

11楼: Originally posted by cheng.jie at 2016-02-19 15:37:55
请问去哪可以下到1stopt软件？这个软件我在网上搜了好久没有安装包...

你是做微分方程参数拟合，得用这个版本：http://muchong.com/bbs/viewthread.php?tid=8050245

1stOpt软件常见问题与解答：http://muchong.com/bbs/viewthread.php?tid=9176191

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

12楼2016-02-19 15:47:40

cheng.jie

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 136.8
帖子: 33
在线: 18.5小时
虫号: 2347007
注册: 2013-03-14
性别: GG
专业: 资源化工

引用回帖:

12楼: Originally posted by 月只蓝 at 2016-02-19 15:47:40
你是做微分方程参数拟合，得用这个版本：http://muchong.com/bbs/viewthread.php?tid=8050245

1stOpt软件常见问题与解答：http://muchong.com/bbs/viewthread.php?tid=9176191...

好的，谢谢！请问我这个问题用matlab好解决吗

化学工程

13楼2016-02-19 16:11:06

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

引用回帖:

13楼: Originally posted by cheng.jie at 2016-02-19 16:11:06
好的，谢谢！请问我这个问题用matlab好解决吗...

不算难。
不过这种问题，MATLAB效果不一定有1stOpt好。

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

14楼2016-02-19 16:16:11

cheng.jie

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 136.8
帖子: 33
在线: 18.5小时
虫号: 2347007
注册: 2013-03-14
性别: GG
专业: 资源化工

引用回帖:

14楼: Originally posted by 月只蓝 at 2016-02-19 16:16:11
不算难。
不过这种问题，MATLAB效果不一定有1stOpt好。...

能否帮忙看一下我的matlab程序，我做了俩，一个是积分法求方程（文件已上传），一个是微分法，微分法比较简单，但是拟合效果太差，积分法我不会，刚用matlab，有些东西不是很理解，万分感谢

化学工程

15楼2016-02-19 16:26:06

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

引用回帖:

15楼: Originally posted by cheng.jie at 2016-02-19 16:26:06
能否帮忙看一下我的matlab程序，我做了俩，一个是积分法求方程（文件已上传），一个是微分法，微分法比较简单，但是拟合效果太差，积分法我不会，刚用matlab，有些东西不是很理解，万分感谢...

你应该参照的是黄华江那本书上的代码。
可成功运行的MATLAB代码，可参考：http://muchong.com/bbs/viewthread.php?tid=6425538

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

16楼2016-02-19 16:37:21

cheng.jie

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 136.8
帖子: 33
在线: 18.5小时
虫号: 2347007
注册: 2013-03-14
性别: GG
专业: 资源化工

引用回帖:

16楼: Originally posted by 月只蓝 at 2016-02-19 16:37:21
你应该参照的是黄华江那本书上的代码。
可成功运行的MATLAB代码，可参考：http://muchong.com/bbs/viewthread.php?tid=6425538...

是啊，我用您推荐的1stopt，怎么算不出来，您看问题出在哪？麻烦老师了

ds.png

化学工程

17楼2016-02-19 16:56:27

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

引用回帖:

17楼: Originally posted by cheng.jie at 2016-02-19 16:56:27
是啊，我用您推荐的1stopt，怎么算不出来，您看问题出在哪？麻烦老师了http://

ds.png
...

文本形式的代码

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

18楼2016-02-19 16:57:37

cheng.jie

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 136.8
帖子: 33
在线: 18.5小时
虫号: 2347007
注册: 2013-03-14
性别: GG
专业: 资源化工

引用回帖:

18楼: Originally posted by 月只蓝 at 2016-02-19 16:57:37
文本形式的代码...

哦哦，是我复制的问题吗，我直接粘过去的，不是输的
Title "Type your title here";
Parameters b1,b2,b3;
Variable t,CA,H2;
ODEFunction CA'=-b1*CA^b2*H2^b3;
Data;
0.2       0.291       0.023
0.3       0.186       0.015
0.4       0.134       0.015
0.5       0.112       0.008
0.6       0.101       0.007

化学工程

19楼2016-02-19 17:04:32

cheng.jie

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 136.8
帖子: 33
在线: 18.5小时
虫号: 2347007
注册: 2013-03-14
性别: GG
专业: 资源化工

引用回帖:

4楼: Originally posted by dingd at 2016-02-19 10:09:08
1stopt求解：

Variable t,CA,H2;
ODEFunction CA'=-b1*CA^b2*H2^b3;
Data;
0.2 0.291 0.023
0.3 0.186 0.015
0.4 0.134 0.015
0.5 0.112 0.008
0.6 0.101 0.007

结果：
均方差(RMSE):0.00100891050 ...

请问您用的1stopt是哪个版本，在网上我下载了一个，按照下面的输入算不出来，是什么原因呢，谢谢！
Title "Type your title here";
Parameters b1,b2,b3;
Variable t,CA,H2;
ODEFunction CA'=-b1*CA^b2*H2^b3;
Data;
0.2       0.291       0.023
0.3       0.186       0.015
0.4       0.134       0.015
0.5       0.112       0.008
0.6       0.101       0.007

化学工程

20楼2016-02-19 22:17:02