24小时热门版块排行榜    

查看: 1098  |  回复: 4

l1003785517

新虫 (小有名气)

[求助] 动力学方程参数计算 已有1人参与

Variables x(1:3),y;
Parameters a(1:6);
Constants R=8.314, T=563.15;
Function y=-a1*exp(-a2/R/T)*(x1^a3)/(1+a4*(x2^a5)*(x3^a6));
Data;
6.34514148        0.0255312        0.01212732        0.007
6.10642476        0.2489292        0.02744604        0.0068
5.97876876        0.3702024        0.03382884        0.0055
5.89260096        0.446796                0.04340304        0.0048
5.7796254        0.5489208        0.0542538        0.0041
5.67941544        0.63828                0.06510456        0.0027
5.59899216        0.7148736        0.06893424        0.002
5.51601576        0.7914672        0.07531704        0.002
5.3711262        0.9318888        0.079785                0.0019
5.25751236        1.0403964        0.08489124        0.0018
5.14708992        1.148904                0.08680608        0.0017
动力学方程中有6个待求参数,a1为指前因子,a2为活化能。
为什么我用这个代码每次运行得到的结果都不一样啊,从书上知道是出现过拟合现象,但是不知道怎么加限制条件,还望各位虫友帮忙!谢谢@月只蓝@beefly
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
仅从数学角度看“a1*exp(-a2/R/T)”项就可用“a1”代替。

各参数的范围得楼主根据自己实际问题的物理意义来给定,他人无法而知。不加范围约束肯定会是多解。

» 本帖已获得的红花(最新10朵)

2楼2017-04-04 11:33:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

l1003785517

新虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by dingd at 2017-04-04 11:33:32
仅从数学角度看“a1*exp(-a2/R/T)”项就可用“a1”代替。

各参数的范围得楼主根据自己实际问题的物理意义来给定,他人无法而知。不加范围约束肯定会是多解。

谢谢回复,这是一个动力学方程,a1是指前因子,a2是活化能,我也知道你说的这一项是可以合并的,但是合并之后我就无法得到这两个参数的具体数值了,所以还想问问您有没有什么办法可以解决这个问题?

发自小木虫IOS客户端
3楼2017-04-08 14:34:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

自己根据各参数的实际物理意义加范围约束。

[ 发自手机版 http://muchong.com/3g ]
4楼2017-04-08 17:12:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

l1003785517

新虫 (小有名气)

引用回帖:
4楼: Originally posted by dingd at 2017-04-08 17:12:57
自己根据各参数的实际物理意义加范围约束。

1.5版本的是不是做不了约束问题?

发自小木虫IOS客户端
5楼2017-04-09 07:19:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 l1003785517 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 05:38 by k06BKNblrGRH
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 04:52 by k06BKNblrGRH
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-17 04:26 by k06BKNblrGRH
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-17 03:52 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-17 03:16 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-17 02:38 by Tql5LQhh5rLK
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 7/350 2026-09-17 02:02 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 6/300 2026-09-17 01:50 by Tql5LQhh5rLK
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-17 01:50 by Tql5LQhh5rLK
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +8 QUjhNVAcOSff 2026-09-13 11/550 2026-09-16 18:37 by rRvxEj3zorCy
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 16:51 by i3KBRD19ERkG
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 15:48 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-16 07:12 by LH5NkK5Ud4bh
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:39 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
信息提示
请填处理意见