problem with leasqr fitting

mrodrig via Help-octave <[email protected]>
Newsgroups gmane.comp.gnu.octave.general
Message-ID <[email protected]>
Hello,

I am trying to perform a fitting and I am obtaining the following results
pfit1 =

   375.11773
    83.26345
     0.29860

kvg = 1
iter =  3
corp1 =

    NA    NA    NA
    NA    NA    NA
    NA    NA    NA

covp1 =

    NA    NA    NA
    NA    NA    NA
    NA    NA    NA

I attach the data for consideration

thanks

datam= xlsread ('data1.xlsx');

temp=datam(1:173,1)+273;
[intensity, lambda] = regdatasmooth (temp, datam(1:173,2));

k=8.6170e-005;

[pks, loc]=findpeaks (intensity,"DoubleSided");

imaxi=intensity(loc(index=find(pks==max(pks))));
tmaxi=temp(loc(index=find(pks==max(pks))));


eacti1= 0.2;

parTLi1 = [tmaxi,imaxi,eacti1];

F1=@(temp,parTLi1)(imaxi*exp(1+parTLi1(3)./(k*temp).*(temp-parTLi1(1))/parTLi1(1)-temp.^2/parTLi1(1).^2*(1-(2*k*parTLi1(1))/parTLi1(3)).*exp(parTLi1(3)./(k*temp).*(temp-parTLi1(1))/parTLi1(1))-(2*k*parTLi1(1)/parTLi1(3))));

wt_tl=ones(size(temp));
iter_tl=800;
tol_tl=0.1;
dFdp_tl="dfdp";
dp_tl=0.0001*ones(size(parTLi1));

%tolerance
tmax_max=tmaxi+2;
tmax_min=tmaxi-2;
imax_max=imaxi+100;
imax_min=imaxi-100;
eact_max=1.3;
eact_min=0.5;

bounds= [tmax_min,tmax_max; imax_min, imax_max; eact_min,eact_max];


f2=0;

[f2, pfit1, kvg, iter, corp1, covp1, covr1,r2] = leasqr (temp, intensity,
parTLi1, F1,tol_tl, iter_tl,wt_tl, dp_tl, dFdp_tl, bounds)

data1.xlsx <https://octave.1599824.n4.nabble.com/file/t372706/data1.xlsx>  




--
Sent from: https://octave.1599824.n4.nabble.com/Octave-General-f1599825.html
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.