24小时热门版块排行榜    

查看: 1410  |  回复: 8

nuonuoy07

禁虫 (初入文坛)

本帖内容被屏蔽

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
简单方便就用1stOpt试试,本论文也用不少例子。如果数据和公式放上来可帮你试试。

» 本帖已获得的红花(最新10朵)

2楼2013-11-09 15:07:55
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖


感谢参与,应助指数 +1
dbb627: 金币+1, 感谢参与讨论 2013-11-11 16:39:56
一级反应的特点是lnc-t图为一直线,用excel就可以做。
1)线性拟合,再转换
拟合的结果 lnc=a*t+b
经过转换c=exp[a*t+b]=exp* exp[a*t]= k* exp[a*t]  

2)直接拟合,在excel下做关于c和t的x-y散点图,但完成后,点数局线,在弹出窗口中中点“添加趋势线”;在新窗口“类型”中点指数函数,“选项”中勾选“显示公式”。可得形如c=k*exp(a*t)的方程。
3楼2013-11-09 21:11:47
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
dbb627: 金币+1, 谢谢指导 2013-11-11 16:39:36
MATLAB拟合动力学方程,可参见:
http://muchong.com/bbs/viewthread.php?tid=6425538&authorid=1122189
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
4楼2013-11-10 09:48:39
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nuonuoy07

禁虫 (初入文坛)

本帖内容被屏蔽

5楼2013-11-11 09:20:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nuonuoy07

禁虫 (初入文坛)

送红花一朵
本帖内容被屏蔽

6楼2013-11-11 10:45:54
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


dbb627: 金币+1, 谢谢指导 2013-11-11 16:39:21
你得说清楚数据中哪个是C?,哪个是T?还有公式中“(-E)/RT”是否应该是“(-E)/(R*T)”
7楼2013-11-11 16:13:29
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nuonuoy07

禁虫 (初入文坛)

本帖内容被屏蔽

8楼2013-11-14 17:57:45
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
nuonuoy07: 金币+100, ★★★★★最佳答案, 很有帮助,谢谢 2013-11-18 15:08:38
1stOpt求解代码:
CODE:
VarConstant T1=[116,124,134,145];
Constant R=8.314;
Variable t, C;
ODEFunction C'=k0*exp((-E)/(R*T1))*C^n;
Data;
0        7.682756
8        7.674475
28        7.631414
43        7.532271
58        7.427824
73        7.331353
88        7.198061
103        7.106492
118        7.015722
133        6.919351
148        6.841528
168        6.722987
175        6.673695
188        6.565408
201        6.49893
216        6.438377
231        6.327329
246        6.24636
261        6.1754
276        6.112681
291        6.045913
306        5.9739
321        5.866216
336        5.821641
351        5.759524
363        5.70366
Data;
0        7.686871
20        7.650561
35        7.41855
50        7.232251
65        6.998809
80        6.814956
95        6.640572
110        6.508211
125        6.367512
140        6.237562
155        6.064225
170        5.951161
185        5.811696
200        5.705514
215        5.60015
230        5.486666
245        5.389096
260        5.291871
275        5.198639
290        5.098263
305        5.027874
320        4.948585
335        4.864672
350        4.72527
365        4.602336
380        4.638632
395        4.465021
410        4.295954
425        4.31433
440        4.358915
455        4.228184
470        4.208097
485        4.21552
Data;
0        7.670636
8        7.685203
23        7.628831
33        7.474568
48        7.094107
63        6.719624
78        6.372542
93        6.098449
108        5.810503
123        5.608571
138        5.400397
153        5.227872
168        5.033441
183        4.909425
198        4.754524
213        4.63621
228        4.494776
243        4.393757
258        4.288467
273        4.23253
288        4.081427
303        4.055749
318        3.928262
333        3.90762
348        3.821904
363        3.769493
378        3.666622
393        3.666622
408        3.633595
423        3.599774
Data;
0        7.682569
10        7.620369
25        7.041338
40        6.318462
55        5.799877
70        5.429969
85        5.175477
100        4.971523
120        4.772523
135        4.641462
150        4.556123
165        4.443969
180        4.354723
195        4.252046
210        4.180154

结果:
计算用时(时:分:秒:微秒): 00:00:14:387
均方差(RMSE):0.151311034827501
残差平方和(SSE):2.28950292605693
相关系数(R): 0.994268050822041
相关系数之平方(R^2): 0.98856895688546
确定系数(DC): 0.975639027938885
F统计(F-Statistic): 100.112499451528

参数                  最佳估算
--------------------        -------------
k0        -0.406739756143606
e        7607.96563296904
n        1.91466838951117

====== 结果输出 ======

文件: 数据文件-1
No        t        目标 c        计算 c
1        8        7.674475        7.62268916812342
2        28        7.631414        7.47638562144282
3        43        7.532271        7.37014086750566
4        58        7.427824        7.26674980940958
5        73        7.331353        7.16610061638949
6        88        7.198061        7.06808718430579
7        103        7.106492        6.97260877526312
8        118        7.015722        6.87956968402978
9        133        6.919351        6.78887892896556
10        148        6.841528        6.70044996538686
11        168        6.722987        6.58592121810165
12        175        6.673695        6.54671711915529
13        188        6.565408        6.47508293101024
14        201        6.49893        6.40493462856371
15        216        6.438377        6.32578182436413
16        231        6.327329        6.24848097127953
17        246        6.24636        6.17296874709095
18        261        6.1754        6.0991846633189
19        276        6.112681        6.02707090913635
20        291        6.045913        5.95657220545712
21        306        5.9739        5.8876356684358
22        321        5.866216        5.8202106816798
23        336        5.821641        5.75424877653301
24        351        5.759524        5.68970351984356
25        363        5.70366        5.63905733310901

文件: 数据文件-2
No        t        目标 c        计算 c
1        20        7.650561        7.4425897345322
2        35        7.41855        7.26892701421296
3        50        7.232251        7.10286072826668
4        65        6.998809        6.9439099034303
5        80        6.814956        6.79163305604943
6        95        6.640572        6.64562423709766
7        110        6.508211        6.50550954192637
8        125        6.367512        6.37094402237766
9        140        6.237562        6.24160894826542
10        155        6.064225        6.11720937304836
11        170        5.951161        5.99747196506465
12        185        5.811696        5.88214307119673
13        200        5.705514        5.77098698446956
14        215        5.60015        5.66378439100422
15        230        5.486666        5.56033097507194
16        245        5.389096        5.46043616382055
17        260        5.291871        5.36392199565628
18        275        5.198639        5.27062209832638
19        290        5.098263        5.18038076451655
20        305        5.027874        5.09305211429794
21        320        4.948585        5.00849933506914
22        335        4.864672        4.926593990771
23        350        4.72527        4.84721539313287
24        365        4.602336        4.77025002855984
25        380        4.638632        4.6955910350112
26        395        4.465021        4.62313772386542
27        410        4.295954        4.55279514233077
28        425        4.31433        4.48447367245401
29        440        4.358915        4.4180886632122
30        455        4.228184        4.35356009255257
31        470        4.208097        4.29081225658007
32        485        4.21552        4.22977348338702

文件: 数据文件-3
No        t        目标 c        计算 c
1        8        7.685203        7.50024248057276
2        23        7.628831        7.19956089677627
3        33        7.474568        7.01162728547921
4        48        7.094107        6.74672320184585
5        63        6.719624        6.50034137699418
6        78        6.372542        6.27062991147204
7        93        6.098449        6.05597424751773
8        108        5.810503        5.85496046922264
9        123        5.608571        5.66634519164149
10        138        5.400397        5.48903070117778
11        153        5.227872        5.32204431080594
12        168        5.033441        5.16452112077062
13        183        4.909425        5.01568954805804
14        198        4.754524        4.87485912026904
15        213        4.63621        4.7414101317294
16        228        4.494776        4.61478483917664
17        243        4.393757        4.49447993662914
18        258        4.288467        4.38004009812273
19        273        4.23253        4.27105241592326
20        288        4.081427        4.16714159286996
21        303        4.055749        4.06796577240642
22        318        3.928262        3.97321290993283
23        333        3.90762        3.88259760538259
24        348        3.821904        3.79585833017221
25        363        3.769493        3.71275499250768
26        378        3.666622        3.6330667939306
27        393        3.666622        3.55659033732925
28        408        3.633595        3.48313795272053
29        423        3.599774        3.41253621216435

文件: 数据文件-4
No        t        目标 c        计算 c
1        10        7.620369        7.33215581538417
2        25        7.041338        6.8605005331292
3        40        6.318462        6.44370722363601
4        55        5.799877        6.07284254266887
5        70        5.429969        5.74080175102825
6        85        5.175477        5.44186515531925
7        100        4.971523        5.17137741252791
8        120        4.772523        4.84846353852255
9        135        4.641462        4.63058317560207
10        150        4.556123        4.43070919615031
11        165        4.443969        4.24672613965759
12        180        4.354723        4.07683538674423
13        195        4.252046        3.9194981961732
14        210        4.180154        3.77339057576542
9楼2013-11-14 20:40:01
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nuonuoy07 的主题更新
信息提示
请填处理意见