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客户端

回复此楼

» 猜你喜欢

反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有291人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

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

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

ZXR319

新虫 (小有名气)

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

送红花一朵

引用回帖:

8楼: Originally posted by 独孤神宇 at 2024-09-11 21:58:39
假设 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.9592328836049 ...

老师，能不能麻烦您给个源代码？我自己弄的有问题，拟合出来反应速率和浓度都是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客户端