| 查看: 2112 | 回复: 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前几天请教过版主，但我调了参数后还是很拉跨，麻烦了各位 |
回复此楼» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
UJN物理学专业调剂
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有263人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名
已经有0人回复
-
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
一志愿物理所凝聚态物理352求调剂
已经有0人回复
» 抢金币啦！回帖就可以得到:
-
齐齐哈尔大学材料科学与工程学院 诚招调剂生
+1/978
-
重庆工商大学环境与资源学院/废油资源化技术与装备教育部工程研究中心接收研究生调剂
+5/330
-
安徽师范大学有机化学专业调剂
+1/282
-
上海海洋大学与中国水产科学研究院联培 学术型研究生 水产 (A+学科）有调剂名额
+2/144
-
双一流南京医科大学招计算机、AI、统计、生物信息等方向26年9月入学博士
+1/84
-
浙江理工大学蒋仲庆 物理/新能源/材料 招9名
+2/82
-
东莞理工土木工程和土木水利专业硕士接受调剂
+1/80
-
2026-大连大学～六盘水师范学院联合培养硕士研究生，招收化学、环境等相关专业
+1/41
-
安徽工业大学有调剂名额
+1/39
-
天津理工大学功能晶体研究院（晶体材料全国重点实验室）杰青团队招收2026年硕士研究生
+1/37
-
大连工业大学 高分子材料 接收考研调剂
+1/36
-
江西师大化学工程学院接收2026调剂生
+1/15
-
烟台大学化学化工学院班庆福课题组招收调剂生
+1/11
-
哈工大化工学院李蒙刚课题组招收2026年秋季入学博士生
+1/8
-
三峡大学 招收调剂（数学、物理学、力学、工学）
+1/7
-
【博士招生】广东工业大学国家优青课题组招收2026年博士生（环境/化学/材料/微生物）
+1/5
-
IL-4蛋白在肿瘤免疫抑制性骨髓生成中的作用机制研究
+1/4
-
招聘智慧建造或城乡规划背景硕士生
+1/3
-
调剂招生化学化工（0817、0856）
+1/3
-
武汉纺织大学材料学院包海峰教授课题组招收2026年材料类调剂考生
+1/3
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): 谢谢参与
