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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 理论与计算方法 » 【求助】参数拟合；用lsqnonlin拟合气固催化动力学模型如何得到最优初值、结果？
151/1返回列表
查看: 2192  |  回复: 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的回帖

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的回帖

liufeisheng

新虫 (初入文坛)
大神们，我理解的是，我这个问题就是解非线性超定方程组，要个最优解，不知对不对
赞一下 回复此楼
9楼2021-11-17 21:30:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
★
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenhuanlong

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

liufeisheng

新虫 (初入文坛)
引用回帖:
12楼: Originally posted by chenhuanlong at 2021-11-20 17:07:42
gytik

这是什么？

发自小木虫Android客户端
回复此楼
13楼2021-11-20 23:25:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
xiaoqiuqiu202楼
2021-11-17 18:33   回复  
liufeisheng(金币+1): 谢谢参与
yuekong4楼
2021-11-17 19:08   回复  
liufeisheng(金币+1): 谢谢参与
发自小木虫Android客户端
zzyxiang20126楼
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): 谢谢参与
相关版块跳转 我要订阅楼主 liufeisheng 的主题更新
151/1返回列表
提示: 如果您在30分钟内回复过其他散金贴，则可能无法领取此贴金币
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭