24小时热门版块排行榜    

查看: 1073  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 什么时候开奖? +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[基金申请] 2026国自然函评费到账 +15 羊腰板 2026-08-21 16/800 2026-08-23 10:45 by process2012
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 14/700 2026-08-23 10:30 by wrm
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 07:41 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 10/500 2026-08-23 07:17 by OEbVnUOu01ol
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见