24小时热门版块排行榜

返回列表

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 liufeisheng 的 18 个金币，回帖就立即获得 1 个金币，每人有 1 次机会

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

liufeisheng

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 134.3
帖子: 47
在线: 10.3小时
虫号: 1379323

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