24小时热门版块排行榜    

查看: 1275  |  回复: 4

瑞aaa

新虫 (正式写手)

[求助] 拟合动力学参数已有1人参与

在下目前在搞反应动力学,由于实验室无人会,只能自己摸索,目前知道了动力学方程和一些数据,想拟合出参数,并求一个什么MRE和SRS值,求各位大佬帮忙matlab或1stop拟合一下,缺什么东西可以告诉我,另外有搞动力学的同学可以一起讨论一下,真是太难了,麻烦各位了!!!
参数:k1  K1 KA KB KC KD KE
r=k1(CA*CB-CC*CD/K1)/(1+KA*CA+KB*CB+KC*CC+KD*CD+KE*CE)^2
CA                   CB                    CC                    CD                    CE
3.914545455        2.146060606        4.249745353        2.737478991        0.588533232       
3.828484848        2.232121212        4.308152457        2.706725388        0.56087973       
3.767272727        2.293333333        4.340183215        2.703875994        0.531698367       
3.711515152        2.349090909        4.37840252        2.683194959        0.514160096       
4.119544848        1.941061212        4.053019871        2.925930562        0.596807143       
3.995503636        2.065102424        4.15042417        2.855163176        0.57017023       
3.930325455        2.130280606        4.185406506        2.850376686        0.539974384       
3.867855152        2.192750909        4.230895022        2.821869957        0.522992597       
4.242424242        1.818181818        3.929393939        3.05030303        0.596060606       
4.122233333        1.938372727        4.028575792        2.972130233        0.57505155       
4.036520606        2.024085455        4.090222032        2.934550481        0.550985063       
3.939393939        2.121212121        4.164242424        2.883636364        0.527878788       
4.363636364        1.696969697        3.822139403        3.143599982        0.610018191       
4.181818182        1.878787879        3.973000571        3.023695828        0.579061177       
4.109090909        1.951515152        4.022277992        2.997868258        0.555611325       
4.060606061        2                        4.044611444        3.001686202        0.529459929       
4.44                        1.620606061        3.731246671        3.249021809        0.595489096       
4.303030303        1.757575758        3.84742957        3.153625708        0.574702298
4.181818182        1.878787879        3.930166596        3.109363778        0.536227202       
4.121212121        1.939393939        3.97701748        3.076268069        0.522472026
回复此楼

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

» 猜你喜欢

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

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
r 的数据没有,另  1stOpt是不区分大小写的。
数值计算
2楼2019-03-24 19:49:27
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

瑞aaa

新虫 (正式写手)

送红花一朵
参数:KG  K1 KA KB KC KD KE
r=kG(CA*CB-CC*CD/K1)/(1+KA*CA+KB*CB+KC*CC+KD*CD+KE*CE)^2

r                       CA                   CB                    CC                    CD                    CE
3.2707       3.914545455        2.146060606        4.249745353        2.737478991        0.588533232        
3.2707       3.828484848        2.232121212        4.308152457        2.706725388        0.56087973        
3.2707       3.767272727        2.293333333        4.340183215        2.703875994        0.531698367        
3.2707       3.711515152        2.349090909        4.37840252        2.683194959        0.514160096        
4.5701       4.119544848        1.941061212        4.053019871        2.925930562        0.596807143        
4.5701       3.995503636        2.065102424        4.15042417        2.855163176        0.57017023        
4.5701       3.930325455        2.130280606        4.185406506        2.850376686        0.539974384        
4.5701       3.867855152        2.192750909        4.230895022        2.821869957        0.522992597        
5.938        4.242424242        1.818181818        3.929393939        3.05030303        0.596060606        
5.938         4.122233333        1.938372727        4.028575792        2.972130233        0.57505155        
5.938         4.036520606        2.024085455        4.090222032        2.934550481        0.550985063        
5.938         3.939393939        2.121212121        4.164242424        2.883636364        0.527878788        
6.0373       4.363636364        1.696969697        3.822139403        3.143599982        0.610018191        
6.0373       4.181818182        1.878787879        3.973000571        3.023695828        0.579061177        
6.0373       4.109090909        1.951515152        4.022277992        2.997868258        0.555611325        
6.0373       4.060606061        2                        4.044611444        3.001686202        0.529459929        
6.5347       4.44                        1.620606061        3.731246671        3.249021809        0.595489096        
6.5347       4.303030303        1.757575758        3.84742957        3.153625708        0.574702298
6.5347       4.181818182        1.878787879        3.930166596        3.109363778        0.536227202        
6.5347       4.121212121        1.939393939        3.97701748        3.076268069        0.522472026
麻烦您了!
3楼2019-03-24 20:19:26
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
瑞aaa: 金币+10, ★★★★★最佳答案 2019-03-29 09:47:17
引用回帖:
3楼: Originally posted by 瑞aaa at 2019-03-24 20:19:26
参数:KG  K1 KA KB KC KD KE
r=kG(CA*CB-CC*CD/K1)/(1+KA*CA+KB*CB+KC*CC+KD*CD+KE*CE)^2

r                       CA                   CB                    CC                    CD                    ...

均方差(RMSE): 0.261169780173993
残差平方和(SSR): 1.36419308152264
相关系数(R): 0.975770230455897
相关系数之平方(R^2): 0.952127542643954
修正R平方(Adj. R^2): 0.930032562325779
确定系数(DC): 0.952110049532069
卡方系数(Chi-Square): 0.130267014611591
F统计(F-Statistic): 42.7168754119675

参数        最佳估算
----------        -------------
kg        -0.163893312002954
k1        1.28342164757612
ka        54.4962573347966
kb        -84.1449692576913
kc        67.1452055737064
kd        -71.3880858976849
ke        -209.837213092912

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

No.        实测值 r        计算值 r
1        3.2707        3.4113409520742
2        3.2707        3.16478637059528
3        3.2707        3.31479775083469
4        3.2707        3.1364858910573
5        4.5701        5.06508488701176
6        4.5701        4.72130539974879
7        4.5701        4.95833678567805
8        4.5701        4.86620414657382
9        5.938        5.80949663720816
10        5.938        5.54919114864549
11        5.938        5.53030156430345
12        5.938        5.39069859262174
13        6.0373        6.11525742685167
14        6.0373        5.81601394103538
15        6.0373        5.89662035754303
16        6.0373        6.22021510593907
17        6.5347        6.59899589023096
18        6.5347        6.4584076607694
19        6.5347        6.68931109384878
20        6.5347        6.71179872058994
数值计算
4楼2019-03-24 20:51:56
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

瑞aaa

新虫 (正式写手)

太感谢了,您又让我看到了我那被拒n次后发四区的希望,我还有几个模型需要请您计算,明天发给您,您有时间就动动手,我不急,真的感谢你

发自小木虫Android客户端
5楼2019-03-24 21:12:40
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 瑞aaa 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 博后基金,博管会会提前知道消息吗? +4 yuyiang 2024-06-13 4/200 2024-06-16 11:40 by yangyuzhong4
[教师之家] 每次骚扰女学生的都是院系领导,而不是普通教师,小编们要注意措辞正确 +8 zju2000 2024-06-15 10/500 2024-06-16 09:22 by Quakerbird
[考博] 34岁读博士晚吗 +26 emitdne 2024-06-13 26/1300 2024-06-16 07:16 by liyeqik
[找工作] 江西双非一本和四川双一流高校如何选择? 5+7 寒山敲钟 2024-06-12 23/1150 2024-06-16 06:52 by puterde
[基金申请] Nature 11日发文,中国著名学者们称造假迫不得已 +5 babu2015 2024-06-14 5/250 2024-06-15 20:40 by lc231001
[教师之家] 请问事业编制和年薪制冲突吗? +6 ZHONGWU_U 2024-06-14 6/300 2024-06-15 20:16 by Ermito
[论文投稿] 投稿时忘记修改一作 +7 gll123456 2024-06-13 11/550 2024-06-15 11:49 by gll123456
[基金申请] 有没有机械的前辈分享一下评上海优都是什么成果啊 +7 wulala800 2024-06-10 7/350 2024-06-15 09:33 by 晓目崇
[论文投稿] 投了一篇4区的SCI,审稿人一个拒稿,一个小修,编辑给了大修。 +9 安稳22123 2024-06-13 10/500 2024-06-14 23:45 by jurkat.1640
[食品] 食品博士导师 +6 小李醒yy 2024-06-11 9/450 2024-06-14 23:37 by 小李醒yy
[基金申请] 博士后基金需要结题吗? +8 zhouchuck 2024-06-13 8/400 2024-06-14 17:27 by liuyupu132
[基金申请] 工材E口JQ有消息了吗 +4 babyduck 2024-06-11 4/200 2024-06-14 17:23 by firepick
[考博] 申博找导师 +4 疏影横斜水清浅3 2024-06-13 6/300 2024-06-14 14:31 by zxl_1105
[有机交流] ππ堆积会发生在有机溶剂中吗 5+3 zibuyu0420 2024-06-13 4/200 2024-06-14 14:17 by 小肉干
[基金申请] 国自然基金公布的时候基金号有吗 +8 潇洒怡惜 2024-06-13 11/550 2024-06-14 11:24 by JRfei
[论文投稿] ACS Nano投稿后分配到副编辑手里12天了,能确定送审了吗? +5 潇洒怡惜 2024-06-12 10/500 2024-06-14 09:56 by 潇洒怡惜
[考博] 还能读上博士吗 +14 小伙走起 2024-06-09 24/1200 2024-06-14 06:32 by 18787065176
[硕博家园] 科研求助 +5 杲www 2024-06-12 6/300 2024-06-13 16:16 by 姓李名明
[论文投稿] 摩擦磨损论文投稿 +3 jmysan 2024-06-12 3/150 2024-06-13 08:36 by 莱茵润色
[教师之家] 公办双非,学生论文升学就业都不行。一本大学的论文升学就业没问题吧? +5 河西夜郎 2024-06-09 5/250 2024-06-10 17:59 by yyallen2003
信息提示
请填处理意见