24小时热门版块排行榜    

查看: 1836  |  回复: 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
回复此楼
已阅   关注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 CrisMessi 2026-08-18 4/200 2026-08-19 00:53 by gatelove
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 27/1350 2026-08-18 21:44 by Noways
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 9/450 2026-08-18 20:09 by 奶牛小黑
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳变了,能看出什么问题? +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 重要消息,中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[基金申请] 时间戳又变了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 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 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
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见