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

回复此楼