24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1513  |  回复: 10
【悬赏金币】回答本帖问题,作者ZXR319将赠送您 40 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ZXR319

新虫 (小有名气)

[求助] 请会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 拟合的思路,目标函数就是计算值和实验测量值的差(比如残差平方和)最小。

请会1stopt 的大佬帮帮我


请会1stopt 的大佬帮帮我-1


@月只蓝 发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZXR319

新虫 (小有名气)

送红花一朵
引用回帖:
5楼: Originally posted by hzlhm at 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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

独孤神宇

版主 (知名作家)

数值计算
2楼2024-09-09 13:31:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZXR319

新虫 (小有名气)

引用回帖:
2楼: Originally posted by 独孤神宇 at 2024-09-09 13:31:42
模型缺少 r的数据。

老师您的意思是数据里每一个时间下应该写出来对应的反应速率吗?我通过实验只能得到不同时间的反应物浓度和生成物浓度,那在1stopt 中可以通过代码得到不同时刻的反应速率吗?我在 Matlab 中是通过 ODE45函数求反应速率的

发自小木虫Android客户端
3楼2024-09-09 14:57:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZXR319

新虫 (小有名气)

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见