| 查看: 1511 | 回复: 10 | |||
| 【悬赏金币】回答本帖问题，作者ZXR319将赠送您 40 个金币 | |||
[求助]
请会1stopt 的大佬帮帮我已有2人参与
|
|||
|
我在做 lh 机理动力学，已经有一部分数据了，目前已经用 matlab 处理了，拟合初值是随便给的，我现在想用1stopt 软件拟合一下，拟合的结果作为 matlab 的参数初值，以便更准确。求大佬帮我拟合一下，最好有源代码，金币我可以尽量给，有多少都可以给，好商量。 数据如下：t ca cb 0 0.1679 0 15 0.1511 0.0164 30 0.1322 0.0376 45 0.115 0.0512 60 0.1037 0.06 90 0.1003 0.069 120 0.09486 0.07368 还需要提供什么可随时联系我，拜谢各位！ 图2是我用 matlab 拟合的思路，目标函数就是计算值和实验测量值的差（比如残差平方和）最小。   @月只蓝 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
求2026年在台湾举行的物理和材料领域国际学术会议信息
已经有0人回复
-
求国际会议网站
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有94人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有19人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31014.8
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 855.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
2楼2024-09-09 13:31:42
|
老师您的意思是数据里每一个时间下应该写出来对应的反应速率吗？我通过实验只能得到不同时间的反应物浓度和生成物浓度，那在1stopt 中可以通过代码得到不同时刻的反应速率吗？我在 Matlab 中是通过 ODE45函数求反应速率的 发自小木虫Android客户端 |
3楼2024-09-09 14:57:27
|
Function r = (k * KA * cA * sqrt(KH * PH2 / Kbc)) / ((1 + KA * cA + KB * cB + sqrt(KH * PH2 / Kbc))^2); 这是模型方程的代码形式 发自小木虫Android客户端 |
4楼2024-09-09 15:28:45
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 18111.5
- 红花: 53
- 帖子: 2879
- 在线: 603.7小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
5楼2024-09-09 22:28:55
送红花一朵|
谢谢老师我做的是加氢反应，1molA 和3mol 氢气生成1molB，A 和 B 摩尔比是1:1，所以是不是-dcA/dt=r；dcB/dt=r。不知道这样可以计算了吗？ 发自小木虫Android客户端 |
6楼2024-09-10 08:18:34
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31014.8
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 855.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
假设 dcB/dt=r 模型公式: cb' = (k*ka*ca*sqrt(kh*3/kbc))/((1+ka*ca+kb*cb+sqrt(kh*3/kbc))^2) 残差平方和(SSE): 0.000100907681851596 均方差(RMSE): 0.00410096902881087 相关系数(R): 0.978424525403226 相关系数之平方(R^2): 0.957314551910527 修正R平方(Adj. R^2): 0.928857586517545 确定系数(DC): 0.956042109592344 F统计(F-Statistic): 5.07803520426272 参数 最佳估算 ---- -------- k 3947916360160.45 ka 4.19389648022836E15 kh 0.406855254346946 kbc 8.49921244458393 kb 1.28426131584115E16 ====== 结果输出 ====== 文件: 数据文件-1 No t 目标 cb 计算 cb 1 15 0.0164 0.023012428915109 2 30 0.0376 0.0370171975730106 3 45 0.0512 0.0473206201899215 4 60 0.06 0.05546914388147 5 90 0.069 0.068244316554898 6 120 0.07368 0.0782291213083908 |
7楼2024-09-11 21:54:51
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31014.8
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 855.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
|
假设 dcA/dt=-r 模型公式: ca' = -(k*ka*ca*sqrt(kh*3/kbc))/((1+ka*ca+kb*cb+sqrt(kh*3/kbc))^2) 残差平方和(SSE): 0.000192586899315775 均方差(RMSE): 0.005665493495948 相关系数(R): 0.959232883604991 相关系数之平方(R^2): 0.920127724989146 修正R平方(Adj. R^2): 0.866879541648577 确定系数(DC): 0.91769457654023 F统计(F-Statistic): 2.54023664039648 参数 最佳估算 ---- -------- k 2837450.21546848 ka 2.67743713648235E-7 kh 67.0303631498266 kbc 0.100651239217611 kb 687.010099643109 ====== 结果输出 ====== 文件: 数据文件-1 No t 目标 ca 计算 ca 1 15 0.1511 0.143509362937322 2 30 0.1322 0.129868469396974 3 45 0.115 0.120093597136612 4 60 0.1037 0.112222162940757 5 90 0.1003 0.0997387665702967 6 120 0.09486 0.0893241064493534 |
8楼2024-09-11 21:58:39
9楼2024-09-13 09:34:45
送红花一朵|
老师，能不能麻烦您给个源代码？我自己弄的有问题，拟合出来反应速率和浓度都是0。Title "动力学模型2"; Parameters k,KA,KB,KH,Kbc; Variable t,CA,CB; Constant PH2=3; Function r = (k * K_A * c_A * sqrt(K_H * P_H2 / K_bc)) / ((1 + K_A * c_A + K_B * c_B + sqrt(K_H * P_H2 / K_bc))^2); Data; t CA CB 0 0.1679 0 30 0.1314 0.0367 60 0.0964 0.0719 90 0.0759 0.096 120 0.0575 0.1113 STARTPROGRAM var i: integer; begin for i = 1 to length(t)-1 do r = -(c_A - c_A[i-1]) / (t - t[i-1]); end ENDPROGRAM FIT 这是Ai 给出的代码，我不知道怎么改 发自小木虫Android客户端 |
10楼2024-09-15 09:34:32