24小时热门版块排行榜    

查看: 1072  |  回复: 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的回帖

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的回帖
查看全部 5 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[教师之家] 跳槽后在研项目怎么办? +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
请填处理意见