nlifit2.m

这是天津大学刘则毅老师编写的那本书中的源代码

lab=[.228 .821 .495 .678 .727 .458 .002 .231 .631 .811 ...
     .050 .483 .116 .456 .358 .258 .767 .487 .748 .695 ...
     .981 .429 .664 .059 .758 .823 .682 .440 .342 .162];
cap=[.802 .771 .758 .452 .845 .084 .295 .546 .017 .223 ...
     .161 .836 .930 .185 .485 .249 .511 .425 .817 .958 ...
     .021 .277 .129 .906 .145 .006 .521 .495 0.092 .934]; 
logq=[-1.359 .193 -.165 -.473 -.563 -2.218 -5.586 -1.315 -3.879 -1.377 ...
   -2.539 -.324 -1.530 -1.151 -.951 -1.695 -.649 -.270 .031 -.125 ...
   -3.633 -.773 -1.678 -2.301 -2.270 -5.150 -.253 -.614 -2.089 -1.275];
x=[lab' cap'];
b0=[1 -1 .5 -1]
[b,r,j]= nlinfit(x,logq,'ces',b0);
b
[yp,ci]=nlpredci('ces',x,b,r,j);
yy=[x logq' yp r yp-ci yp+ci]
nlinto(x, logq, 'ces',b0);