| 查看: 1581 | 回复: 7 | ||
xiaohsu2006木虫 (著名写手)
|
[求助]
求1stopt拟合动力学参数
|
|
小弟没有高版本的，求高手帮忙，感激不尽！ 反应动力学： A+B - 2.5C+3.5D+0.5F+G A+D - C+F 3C+D - E+A r1=k1*CB r2=k2*CD r3=k3*CD dCAdt=-r1-r2+r3 dCBdt=-r1 dCCdt=2.5*r1+r2-3*r3 dCDdt=3.5*r1-r2-r3 dCEdt=r3 dCFdt=0.5r1+r2 dCG=r1 初始条件 CA=4.0799, CB=0.0707，CC=CD=CE=CF=CG=0 1 t/s CC CD CE CF 30 2.35E-04 3.36E-03 1.15E-04 4.53E-03 60 1.78E-03 8.62E-03 9.35E-04 6.27E-03 120 2.92E-03 1.20E-02 2.27E-03 6.74E-03 180 4.43E-03 1.54E-02 2.93E-03 6.99E-03 2 t/s CC CD CE CF 30 7.33E-04 5.87E-03 4.50E-04 5.01E-03 60 3.07E-03 1.36E-02 2.67E-03 6.25E-03 120 5.07E-03 1.66E-02 4.42E-03 6.62E-03 180 5.86E-03 1.71E-02 5.30E-03 7.29E-03 3 t/s CC CD CE CF 30 1.67E-03 9.37E-03 1.08E-03 4.77E-03 60 5.55E-03 1.68E-02 4.33E-03 5.91E-03 120 8.29E-03 1.77E-02 6.02E-03 7.50E-03 180 1.00E-02 1.68E-02 6.28E-03 8.12E-03 分别拟合出三个温度下的反应速率常数k1,k2,k3. |
回复此楼» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有136人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有11人回复
-
散金币祈福
已经有94人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 拜求！matlab拟合动力学参数
已经有14人回复
- 最小二乘解算拟合方程
已经有12人回复
- 求助用matlab拟合动力学方程
已经有9人回复
- matlab 拟合方法
已经有4人回复
- matlab 拟合反应动力学参数
已经有8人回复
- 1stopt拟合问题求助
已经有7人回复
- rmse，sse计算求助
已经有3人回复
- 如何用matlab对四维数据进行非线性拟合？
已经有5人回复
- 通过实验数据拟合，求解公式中的参数
已经有14人回复
- 用1stopt拟合非线性方程结果与其他软件拟合结果差异大
已经有5人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- 用1stopt或者origin非线性拟合磁导率
已经有5人回复
- matlab数据拟合问题求助
已经有7人回复
- matlab求解动力学 方程中的参数
已经有13人回复
- 动力学参数拟合
已经有26人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- MATLAB拟合工具如何引入误差棒？
已经有11人回复
- 学1stopt要多久？
已经有12人回复
- 求可以做曲面拟合的软件
已经有10人回复
- 【求助】拟合动力学方程求助
已经有13人回复
- 用1stopt进行数据拟合时，无法运行sharedmodel ,求助高手
已经有7人回复
xiaohsu2006
木虫 (著名写手)
- 应助: 23 (小学生)
- 金币: 2520.4
- 散金: 381
- 红花: 4
- 帖子: 1362
- 在线: 1137.1小时
- 虫号: 1012350
- 注册: 2010-05-06
- 性别: GG
- 专业: 多相流热物理学
2楼2014-03-24 10:23:41
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
xiaohsu2006: 金币+90, ★★★★★最佳答案 2014-03-24 11:33:21
感谢参与，应助指数 +1
xiaohsu2006: 金币+90, ★★★★★最佳答案 2014-03-24 11:33:21
|
代码及第一组数据结果： 均方差(RMSE):0.00206339294816877 残差平方和(SSE):6.81214473368416E-5 相关系数(R): 0.857496435844223 相关系数之平方(R^2): 0.735300137485547 确定系数(DC): -2.62619250526438 F统计(F-Statistic): 0.302254888086507 参数 最佳估算 -------------------- ------------- k1 0.000571090248522173 k2 0.00237024445720707 k3 0.00354220354969495 ====== 结果输出 ====== 文件: 数据文件-1 No t 目标 cc 计算 cc 目标 cd 计算 cd 目标 ce 计算 ce 目标 cf 计算 cf 1 30 0.000235 0.00250994279834437 0.00336 0.00385071637240616 0.000115 0.000211282929925092 0.00453 0.000741861291302684 2 60 0.00178 0.0041032640107171 0.00862 0.00701015349688254 0.000935 0.000793940277216935 0.00627 0.00172202530354545 3 120 0.00292 0.00514480375607497 0.012 0.0116906247177385 0.00227 0.00281539943089081 0.00674 0.00422532664962184 4 180 0.00443 0.00410080160481611 0.0154 0.0147450856032192 0.00293 0.00564846729296284 0.00699 0.00723295135707574 |
» 本帖已获得的红花（最新10朵）
xiaohsu20063楼2014-03-24 10:57:31
xiaohsu2006
木虫 (著名写手)
- 应助: 23 (小学生)
- 金币: 2520.4
- 散金: 381
- 红花: 4
- 帖子: 1362
- 在线: 1137.1小时
- 虫号: 1012350
- 注册: 2010-05-06
- 性别: GG
- 专业: 多相流热物理学
4楼2014-03-24 11:35:56
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
xiaohsu2006: 金币+200, ★★★★★最佳答案 2014-03-24 13:41:29
xiaohsu2006: 金币+200, ★★★★★最佳答案 2014-03-24 13:41:29
|
第二组： 均方差(RMSE):0.00231720632138809 残差平方和(SSE):8.59111221740947E-5 相关系数(R): 0.868731725547655 相关系数之平方(R^2): 0.754694810973007 确定系数(DC): -2.87359587747822 F统计(F-Statistic): 0.553612983479765 参数 最佳估算 -------------------- ------------- k1 0.000745071719863344 k2 0.00228132894455383 k3 0.00353913908564415 ====== 结果输出 ====== 文件: 数据文件-1 No t 目标 cc 计算 cc 目标 cd 计算 cd 目标 ce 计算 ce 目标 cf 计算 cf 1 30 0.000733 0.00325877285433546 0.00587 0.00501713253581565 0.00045 0.000275173337049803 0.00501 0.000958759952701426 2 60 0.00307 0.00529425321594842 0.0136 0.00911952654539591 0.00267 0.00103303889142659 0.00625 0.00221139144787796 3 120 0.00507 0.00650713396046731 0.0166 0.0151522454277829 0.00442 0.00365541671324483 0.00662 0.00537970201979195 4 180 0.00586 0.00495183891595699 0.0171 0.0190254293617908 0.0053 0.00731587838390978 0.00729 0.00915254679655734 第三组： 均方差(RMSE):0.00290108831059639 残差平方和(SSE):0.000134661014174064 相关系数(R): 0.87612822607638 相关系数之平方(R^2): 0.767600668527745 确定系数(DC): -1.04816539387532 F统计(F-Statistic): 2.70389954875217 参数 最佳估算 -------------------- ------------- k1 0.000850324466622732 k2 0.00300797213925827 k3 0.00339594133977107 ====== 结果输出 ====== 文件: 数据文件-1 No t 目标 cc 计算 cc 目标 cd 计算 cd 目标 ce 计算 ce 目标 cf 计算 cf 1 30 0.00167 0.00381899279623351 0.00937 0.00566811232455994 0.00108 0.000299316725528757 0.00477 0.00115548563185028 2 60 0.00555 0.00643102513322945 0.0168 0.0102027330692526 0.00433 0.00111646934774825 0.00591 0.00274722131198251 3 120 0.00829 0.00889289413916752 0.0177 0.0166430893826131 0.00602 0.00390344884962953 0.0075 0.0068866477146232 4 180 0.01 0.00874715868387847 0.0168 0.0205465263667531 0.00628 0.00772723662236323 0.00812 0.0118613247042733 |
5楼2014-03-24 13:17:56
xiaohsu2006
木虫 (著名写手)
- 应助: 23 (小学生)
- 金币: 2520.4
- 散金: 381
- 红花: 4
- 帖子: 1362
- 在线: 1137.1小时
- 虫号: 1012350
- 注册: 2010-05-06
- 性别: GG
- 专业: 多相流热物理学
|
你好，能不能帮我再跑个代码，写好了！麻烦你了，先谢过！ ConstStr r1=k1*CB,r2=k2*CD,r3=k3*CD; Variable t,CC,CD,CE,CF; InitialODEValue t=0, CA=4079.9, CB=70.7,CC=0, CD=0, CE=0, CF=0, CG=0; ODEFunction CA'=-r1-r2+r3; CB'=-r1; CC'=2.5*r1+r2-3*r3; CD'=3.5*r1-r2-r3; CE'=r3; CF'=0.5*r1+r2; CG'=r1; Data; //t/s CC CD CE CF 30 0.235008523 3.363239986 0.115121201 4.531327956 60 1.780010877 8.616751392 0.934615204 6.268066726 120 2.923654401 12.0012815 2.266319353 6.737997984 180 4.433302192 15.40323797 2.929517334 6.994715676 |
6楼2014-04-12 19:02:41
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
|
均方差(RMSE):2.06284664311107 残差平方和(SSE):68.0853803679137 相关系数(R): 0.858287409589619 相关系数之平方(R^2): 0.736657277460058 确定系数(DC): -2.62363520042484 F统计(F-Statistic): 0.302815847593867 参数 最佳估算 -------------------- ------------- k1 0.000571141146818135 k2 0.00237051259796401 k3 0.00354131349016907 ====== 结果输出 ====== 文件: 数据文件-1 No t 目标 cc 计算 cc 目标 cd 计算 cd 目标 ce 计算 ce 目标 cf 计算 cf 1 30 0.235008523 2.51033673349595 3.363239986 3.85109149173353 0.115121201 0.211249848393764 4.531327956 0.74194373907573 2 60 1.780010877 4.10425865484917 8.616751392 7.01089022534024 0.934615204 0.793820390023803 6.268066726 1.72224301149238 3 120 2.923654401 5.14738702699264 12.0012815 11.6919857722601 2.266319353 2.81501817652468 6.737997984 4.22595972144955 4 180 4.433302192 4.10534301010109 15.40323797 14.746973106494 2.929517334 5.64773674450559 6.994715676 7.23413221207872 |
7楼2014-04-12 19:52:02
xiaohsu2006
木虫 (著名写手)
- 应助: 23 (小学生)
- 金币: 2520.4
- 散金: 381
- 红花: 4
- 帖子: 1362
- 在线: 1137.1小时
- 虫号: 1012350
- 注册: 2010-05-06
- 性别: GG
- 专业: 多相流热物理学
8楼2014-05-16 14:14:51
