24小时热门版块排行榜    

查看: 2462  |  回复: 14

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
独孤神宇: 金币+1 2018-07-05 14:30:10
l1003785517: 金币+20, ★★★很有帮助 2018-07-05 15:28:26
CODE:
Title "kinetic values";
Constant R=8.314;
VarConstant T=[275,282,287,292];
ConstStr a1=k1*exp(-Ea1/R/T),a2=k2*exp(-Ea2/R/T);
Variables x(1:3),y(1:2);
Parameters a(3:5);
ParameterDomain = [0,];
ShareModel;
Function y1=a1*x1/(1+a3*x1+a4*x2+a5*x3);
         y2=(a2*x2-a1*x1)/(1+a3*x1+a4*x2+a5*x2);
Data; //T=275
6.31317        0.029302 0.007644 0.0574 0.0345
6.29412        0.044590 0.011466 0.0407 0.0298
6.27507        0.059878 0.015288 0.0241 0.0251
6.26237        0.069433 0.018473 0.0278 0.0203
6.24713        0.080899 0.022295 0.0315 0.0208
6.23316        0.091091 0.026117 0.0268 0.0194
6.22109 0.099372 0.029939 0.0220 0.0160
6.21093        0.107016 0.032487 0.0172 0.0127
Data; //T=282
6.14807        0.133133        0.069433        0.2894        0.1089
6.062345        0.186004        0.102557        0.1519        0.0714
6.021705        0.202566        0.126763        0.0144        0.0338
5.991225        0.212758        0.147147        0.0425        -0.0038
5.968365        0.220402        0.162435        0.0705        0.0217
5.937885        0.234416        0.178997        0.0575        0.0297
5.91566        0.244608        0.1911        0.0445        0.0201
5.892165        0.254163        0.205114        0.0316        0.0105
Data; //T=287
6.0452        0.187278        0.118482        0.4615        0.1918
5.89026        0.277095        0.184093        0.273        0.1397
5.80136        0.326781        0.223587        0.0845        0.0876
5.74675        0.356083        0.249067        0.0829        0.0356
5.69976        0.377741        0.274547        0.0813        0.0243
5.67182        0.386022        0.294294        0.0677        0.0176
5.648325        0.390481        0.313404        0.0541        0.0152
5.614035        0.400673        0.33761        0.0404        0.0129
Data; //T=292
5.90677        0.220402        0.224224        0.6236        0.1761
5.741035        0.303212        0.307671        0.4216        0.1584
5.559425        0.383474        0.409591        0.2197        0.1408
5.434965        0.437619        0.480298        0.2059        0.1231
5.311775        0.497497        0.543998        0.1921        0.0677
5.26669        0.514696        0.572026        0.1219        0.0332
5.229225        0.527436        0.596869        0.0817        0.0195
5.201285        0.53508        0.617253        0.0585        0.0058

结果不好,自己好好检查下公式和数据:

Root of Mean Square Error (RMSE): 0.0829320200375007
Sum of Squared Residual: 0.440174076640027
Correlation Coef. (R): 0.0189691245411682
R-Square: 0.000359827685858348
Adjusted R-Square: 0.51886698386634
Determination Coef. (DC): -2.88389955816035
F-Statistic: 0.226197729277559

Parameter                  Best Estimate
--------------------        -------------
k1        4.60347494177944E35
ea1        166179.343552407
k2        9.73280643547253E23
ea2        93748.2902053327
a3        65737.9703181681
a4        4.29079677286443E-8
a5        49386056.7270752
11楼2018-07-04 23:05:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

引用回帖:
6楼: Originally posted by l1003785517 at 2018-07-04 12:55:03
是的,作为一个求助者真的不应该出现这种错误,很抱歉。
我在代码中给出的公式,Function y1=a1*x1/(1+a3*x1+a4*x2+a5*x3);
                                                    y2=(a2*x2-a1*x1)/(1+a3*x1+a4 ...

你是四组数据对应四组a,k和ea吧。
一团和气
12楼2018-07-05 11:21:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

l1003785517

新虫 (小有名气)

引用回帖:
11楼: Originally posted by dingd at 2018-07-04 23:05:52
Title "kinetic values";
Constant R=8.314;
VarConstant T=;
ConstStr a1=k1*exp(-Ea1/R/T),a2=k2*exp(-Ea2/R/T);
Variables x(1:3),y(1:2);
Parameters a(3:5);
ParameterDomain = ;
...

非常感谢您的帮助,我想问一下:我把您给的代码在1.5版上运行的结果和您给的结果为什么不一样呢?而且每一次运行结果都不同,是因为软件版本太低的原因吗?
您说的结果不好应该是公式的问题,我自己推了几个可能的公式。您能分别帮我算一下这四个公式吗?谢谢!
(1):Function y1=a1*x1/(1+a3*x1+a4*x2+a5*x3);
                           y2=a2*x1/(1+a3*x1+a4*x2+a5*x3);
(2):Function y1=a1*a3*x1/(1+a3*x1+a4*x2+a5*x3)^4;
                           y2=a2*a3*x1/(1+a3*x1+a4*x2+a5*x3)^6;
(3):Function y1=a1*a3*x1/(1+a3*x1+a4*x2+a5*x3)^7;
                           y2=a2*a3*x1/(1+a3*x1+a4*x2+a5*x3)^11;
(4):Function y1=a1*a4*x2/(1+a3*x1+a4*x2+a5*x3);
                           y2=a2*a4*x3/(1+a3*x1+a4*x2+a5*x3);
13楼2018-07-05 15:41:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

l1003785517

新虫 (小有名气)

引用回帖:
12楼: Originally posted by lipenggg at 2018-07-05 11:21:36
你是四组数据对应四组a,k和ea吧。...

我是希望通过这四组数据得到一组参数的,用matlab拟合动力学方程需要给出参数的初值,只能由1stopt计算了,不知道还有什么办法可以实现!非常感谢您的回复
14楼2018-07-05 15:42:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zy_1019

新虫 (初入文坛)

引用回帖:
11楼: Originally posted by dingd at 2018-07-04 23:05:52
Title "kinetic values";
Constant R=8.314;
VarConstant T=;
ConstStr a1=k1*exp(-Ea1/R/T),a2=k2*exp(-Ea2/R/T);
Variables x(1:3),y(1:2);
Parameters a(3:5);
ParameterDomain = ;
...

你好,我想问一下,算出来的结果怎么看好不好?
15楼2018-07-25 10:19:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 l1003785517 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 看来今天不会放榜了? +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 紫杉醇
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +15 医学老男孩 2026-08-20 18/900 2026-08-21 13:46 by 医学老男孩
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 7/350 2026-08-21 12:15 by weiyin
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 今天基金会出结果吗?20260819 +15 kkkl_v 2026-08-19 16/800 2026-08-21 11:22 by 365372687
[基金申请] 科研孤儿太难了 +15 我4大白菜 2026-08-20 16/800 2026-08-21 11:07 by yming1319
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳今天,20号变了 +4 archvillain 2026-08-20 4/200 2026-08-21 08:10 by gatelove
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 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 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 今天维护系统维护 祝所有人 高中 +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 蔡棒棒菂
信息提示
请填处理意见