24小时热门版块排行榜    

查看: 1068  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天基金会出结果吗?20260819 +13 kkkl_v 2026-08-19 13/650 2026-08-19 22:37 by 宝贝虫子
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 aOn84t7OD1JN 2026-08-19 4/200 2026-08-19 22:31 by BqkHsOAGISWU
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 PnssoZglMj3v 2026-08-19 5/250 2026-08-19 22:19 by BqkHsOAGISWU
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 PnssoZglMj3v 2026-08-19 4/200 2026-08-19 22:16 by BqkHsOAGISWU
[基金申请] 时间戳变了,能看出什么问题? +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +54 医学老男孩 2026-08-13 121/6050 2026-08-19 20:28 by 医学老男孩
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 PnssoZglMj3v 2026-08-19 5/250 2026-08-19 17:42 by oCMhZAsgNFoO
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
信息提示
请填处理意见