24小时热门版块排行榜    

查看: 2221  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学319求调剂 +4 锦鲤0909 2026-03-17 4/200 2026-03-17 18:21 by 重科小霸王
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考博] 26申博 +4 八旬速览 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[基金申请] 国自科面上基金字体 +6 iwuli 2026-03-12 7/350 2026-03-16 21:18 by sculhf
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
信息提示
请填处理意见