24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者ZXR319将赠送您 40 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ZXR319

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 74.6
帖子: 72
在线: 7.3小时
虫号: 33303510
注册: 2023-03-13
专业: 化学工程及工业化学

[求助] 请会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论文润色/翻译怎么收费? 已经有140人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有19人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复

1楼2024-09-09 11:24:26

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

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31014.8
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 855.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

感谢参与，应助指数 +1

引用回帖:

6楼: Originally posted by ZXR319 at 2024-09-10 08:18:34
谢谢老师我做的是加氢反应，1molA 和3mol 氢气生成1molB，A 和 B 摩尔比是1:1，所以是不是-dcA/dt=r；dcB/dt=r。不知道这样可以计算了吗？
...

假设 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