| 查看: 2052 | 回复: 14 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 liufeisheng 的 18 个金币 ，回帖就立即获得 1 个金币，每人有 1 次机会 | |||
[交流]
【求助】参数拟合；用lsqnonlin拟合气固催化动力学模型如何得到最优初值、结果？
|
|||
 【求助】参数拟合；用lsqnonlin拟合气固催化动力学模型时不同参数初值得到不同的结果，请问如何得到最优初值→最优结果？若用其他软件，有无类似教程（案例）推荐？希望案例详细直接些，小白不胜感激 学校无1stOpt\Lingo, 有Matlab\Origin Pro\Aspen\Maple 我的Matlab程序： clc, clear; a=textread('CH4.txt'); %读取txt的数据 y0=a(:,1); %提出因变量反应速率r的数据 x0=a(:,[2:5]); %自变量数据矩阵 k0=[5 7.8 4.6 69.4]; %待拟合k(1)-k(4)参数初值 lb=[-inf -inf -inf -inf]; %k(1)-k(4)下边界 ub=[+inf +inf +inf +inf];%k(1)-k(4)上边界 fun=@(k)(k(1).*(a(:,2)).^0.5.*(a(:,3)).^0.5).*(1+k(2).*(a(:,2)).^0.5.*(a(:,3)).^0.5+k(3).*(a(:,2)).^0.5.*(a(:,3)).^(-0.5)+k(4).*a(:,5)).^(-2)-y0; options = optimoptions(@lsqnonlin,'Algorithm','trust-region-reflective'); [k,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(fun,k0,lb,ub,options); ci=nlparci(k,residual,'jacobian',jacobian); 输出结果： k = 0.006998532596973 8.341744654161287 5.000558918360854 78.381545090556088 ci= 1.0e+03 * -0.000233641132182 0.000247638197376 -0.135382312161568 0.152065801469891 -0.115657322370384 0.125658440207105 -1.746709223913331 1.903472314094444 CH4.txt 数据： 3.37298E-05 0.045483955 0.192267913 0.003064331 0.000940799 3.21878E-05 0.046131365 0.192383275 0.002923178 0.001038222 3.1274E-05 0.046183108 0.192564555 0.002841214 0.001038421 3.14155E-05 0.046821577 0.192744839 0.002852714 0.000778737 3.84628E-05 0.047128941 0.194625771 0.054185236 0.000912819 3.52457E-05 0.046876849 0.194669018 0.053872792 0.000782523 3.51931E-05 0.046017958 0.194717184 0.053899781 0.00101082 3.44983E-05 0.046459181 0.194844229 0.053840691 0.001011045 3.88481E-05 0.047733607 0.195747671 0.105901047 0.001046088 4.38977E-05 0.046653568 0.196054833 0.106435268 0.000915325 4.09879E-05 0.046313497 0.195876834 0.106032215 0.000980822 3.9251E-05 0.047132934 0.19592966 0.105862742 0.001078872 2.8978E-05 0.037387191 0.155197021 0.002659132 0.00105032 2.95116E-05 0.037516347 0.155353051 0.00270676 0.000951828 2.79537E-05 0.037481776 0.155383917 0.002564274 0.000952052 2.86171E-05 0.038352003 0.155459011 0.002624085 0.000886453 2.77704E-05 0.036271974 0.155810267 0.043244686 0.00105194 2.81628E-05 0.037183107 0.155909698 0.043267296 0.000723187 2.61699E-05 0.036920894 0.155947996 0.043074371 0.000920381 2.8119E-05 0.036925417 0.156023334 0.043234911 0.001117742 3.58298E-05 0.037319526 0.156426627 0.085190713 0.001283718 3.20259E-05 0.037363173 0.156461043 0.084808392 0.000822899 3.34526E-05 0.037150974 0.156562126 0.084963969 0.000888858 3.64802E-05 0.038366946 0.156522617 0.085222442 0.000921728 拟合模型： r_(〖CH〗_4 )=(kp_(〖CO〗_2)^0.5 p_(H_2)^0.5)/(1+〖k_1 p〗_(〖CO〗_2)^0.5 p_(H_2)^0.5+k_2 (p_(〖CO〗_2)^0.5)⁄(p_(H_2)^0.5 )+k_3 p_CO )^2 感谢回答 :hand S前几天请教过版主，但我调了参数后还是很拉跨，麻烦了各位 |
回复此楼» 猜你喜欢
帮我的英语口语老师找学生
已经有4人回复
-
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
物理学I论文润色/翻译怎么收费?
已经有96人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有3人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
» 抢金币啦！回帖就可以得到:
-
广东工业大学-木质纤维素高值化利用团队招博士研究生
+1/84
-
山东师范大学有机化学专业胡忠燕老师课题组招收2026届硕士研究生以及调剂生
+1/46
-
福建师范大学化学与材料学院杜克钊团队招生
+1/40
-
南开大学罗景山教授课题组招聘博士后（光/电催化、钙钛矿光伏方向）
+2/38
-
arXiv 预印本求推荐
+2/36
-
广东工业大学管理科学与工程专业博士招生
+1/34
-
中科院化学所 宋延林 课题组招聘合成化学方向博士后（开展打印合成化学方向研究）
+1/33
-
山东理工大学 机械工程学院 多孔金属材料课题组 招收2026年入学博士
+1/32
-
华中科技大学物理学院于涛教授课题组 2026年拟招收优秀考研生1名
+2/30
-
新疆大学招收学硕调剂
+1/29
-
国家级领军人才团队高分子材料方向2026年博士研究生招生
+1/19
-
湖南大学材料科学与工程学院招收博士研究生
+1/9
-
博士后招聘（高薪40万+）
+1/8
-
课题组招收环境及相关专业调剂硕士研究生（欢迎优秀学生加入）
+1/8
-
重庆大学药学院闫海龙课题组拟招收2026年申请考核制博士研究生
+1/6
-
请问为什么我的贴子自动到“已删除的主题”里而没有被发布？
+1/5
-
上海第二工业大学-朱大海课题组招生（过线就能调剂）
+1/5
-
2026年东北石油大学“页岩油气钻采高效井眼清洁”创新团队招硕士生
+1/3
-
哈尔滨工业大学（深圳）绿色化学团队招收2026年秋季入学博士生
+1/2
-
招收环工、环科、化工、应化、化学、生工、生物、材化等相关专业研究生，要求过国家线
+1/2
3楼2021-11-17 18:48:45
5楼2021-11-17 19:18:39
这惜字如金啊7楼2021-11-17 19:53:28
9楼2021-11-17 21:30:28
★
liufeisheng(金币+1): 谢谢参与
liufeisheng(金币+1): 谢谢参与
|
a(:,3)列数据拟合公式里没用上吧，1stOpt稳定唯一解： Root of Mean Square Error (RMSE): 3.34612849534681E-6 Sum of Squared Residual: 2.68717821776925E-10 Correlation Coef. (R): 0.694037081019363 R-Square: 0.481687469829878 Parameter Best Estimate ---------- ------------- k1 0.000137780278403927 k2 0.490633000674274 k3 -0.805989362395384 k4 -53.9028500300834 |
11楼2021-11-18 16:13:47
chenhuanlong
禁虫 (文学泰斗)
★
liufeisheng(金币+1): 谢谢参与
liufeisheng(金币+1): 谢谢参与
|
本帖内容被屏蔽 |
12楼2021-11-20 17:07:42
13楼2021-11-20 23:25:59
简单回复
2021-11-17 18:33
回复
liufeisheng(金币+1): 谢谢参与
yuekong4楼
2021-11-17 19:08
回复
2021-11-17 19:24
回复
liufeisheng(金币+1): 谢谢参与
h 发自小木虫IOS客户端
tzynew8楼
2021-11-17 20:06
回复
liufeisheng(金币+1): 谢谢参与
nono200910楼
2021-11-17 23:53
回复
liufeisheng(金币+1): 谢谢参与
。 发自小木虫Android客户端
miaojiabing14楼
2021-11-23 17:09
回复
liufeisheng(金币+1): 谢谢参与
jjijjj15楼
2021-11-24 10:04
回复
liufeisheng(金币+1): 谢谢参与
