各位高手,请教一个问题,这是我做数值模拟时的程序,当拟合出k的值后发现,其置信区间很大,这种情况怎么处理?请求大家的帮忙!
function HZSM_8
clear all;clc
expdata=[4.28 0.0724 0.0614 0.0075 0.0013 0.0021;
8.56 0.1126 0.0859 0.0164 0.0039 0.0064;
16.69 0.1529 0.1034 0.0278 0.0079 0.0138;
33.76 0.1928 0.1134 0.0412 0.0132 0.0251;
69.72 0.2305 0.1158 0.0555 0.0192 0.04;];
yexp=expdata(:,2:6);
x0=[0 0 0 0 0];
k0=[1.862 0.578 1.321 0.8 9.813 1.018 0.2409 0.626 4.5301 18.174 40.311 314.1];
options=optimset('Algorithm','Levenberg-Marquardt','LargeScale','off');
[k,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(@objfunc,k0,[],[],options,x0,yexp);
ci=nlparci(k,residual,jacobian);
fprintf('\t k1=%.4f±%.4f\n',k(1),ci(1,2)-k(1))
fprintf('\t k2=%.4f±