24小时热门版块排行榜    

查看: 879  |  回复: 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 的主题更新
信息提示
请填处理意见