应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 理论与计算方法 » 【求助】参数拟合；用lsqnonlin拟合气固催化动力学模型如何得到最优初值、结果？
51/1返回列表
查看: 2054  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 liufeisheng 的 18 个金币 ，回帖就立即获得 1 个金币，每人有 1 次机会
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

liufeisheng

新虫 (初入文坛)

[交流] 【求助】参数拟合；用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
感谢回答:handS前几天请教过版主，但我调了参数后还是很拉跨，麻烦了各位
回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2021-11-17 17:59:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenhuanlong

禁虫 (文学泰斗)
★
liufeisheng(金币+1): 谢谢参与
本帖内容被屏蔽
12楼2021-11-20 17:07:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

liufeisheng

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by xiaoqiuqiu20 at 2021-11-17 18:33:40

感谢大佬回复，还望具体指导啊
赞一下 回复此楼
3楼2021-11-17 18:48:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liufeisheng

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by yuekong at 2021-11-17 19:08:00

大佬回复太言简意赅啊，敢问如何搞定？
赞一下 回复此楼
5楼2021-11-17 19:18:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liufeisheng

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by zzyxiang2012 at 2021-11-17 19:24:46
h

这惜字如金啊

发自小木虫Android客户端
回复此楼
7楼2021-11-17 19:53:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
yuekong4楼
2021-11-17 19:08   回复  
liufeisheng(金币+1): 谢谢参与
发自小木虫Android客户端
nono200910楼
2021-11-17 23:53   回复  
liufeisheng(金币+1): 谢谢参与
。 发自小木虫Android客户端
zzyxiang20126楼
2021-11-17 19:24   回复  
liufeisheng(金币+1): 谢谢参与
h 发自小木虫IOS客户端
查看全部 15 个回答
提示: 如果您在30分钟内回复过其他散金贴，则可能无法领取此贴金币
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕323求调剂 +6 李白26 2026-03-07 6/300 2026-03-14 06:05 by JourneyLucky
[考研] 346分材料求调剂 +3 snow_反季节版 2026-03-07 3/150 2026-03-14 05:55 by JourneyLucky
[考研] 求调剂，一志愿山东大学，数二英一，302 +3 nini412 2026-03-08 3/150 2026-03-14 04:23 by JourneyLucky
[考研] 一志愿哈工大材料 初试成绩323 +3 手机用户 2026-03-08 3/150 2026-03-14 03:27 by JourneyLucky
[考研] 313分生物学求调剂 +6 Yyt杨1 2026-03-09 8/400 2026-03-14 03:00 by JourneyLucky
[考研] 293求调剂 +5 上班不着吉 2026-03-09 5/250 2026-03-14 02:37 by JourneyLucky
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[考研] 一志愿北京化工大学材料与化工296分求调剂 +16 稻妻小编 2026-03-09 18/900 2026-03-14 02:00 by JourneyLucky
[考研] 328，0703考生求调剂，一志愿为东北师范大学 +4 观素律 2026-03-09 5/250 2026-03-14 01:24 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 一志愿西南交大，材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
[考研] 研究生招生 +3 徐海涛11 2026-03-10 7/350 2026-03-12 14:26 by 徐海涛11
[考研] 一志愿江南大学085701环境工程专硕总分287求调剂 +5 18266118446 2026-03-09 5/250 2026-03-11 16:51 by 2020015
[考博] 找博导 +4 小呆呆熊 2026-03-07 4/200 2026-03-10 11:08 by 天道酬勤2026686
[考研] 0307化学求调剂 +6 0ok0k 2026-03-07 6/300 2026-03-07 20:10 by pies112
信息提示
关闭
请填处理意见
关闭