f=@(x,xdata) 28.83*xdata(1,./(0.13+xdata(1,+xdata(1,.^2./x(1)+x(2)*xdata(2,);
xdata=[0.03276245 0.037514287 0.02777479 0.035865345 0.033207492 0.039172673 0.058890928 0.06598988 0.063845034 0.070846079 0.060763037 0.083675767 0.148866584 0.153289562 0.150439553 0.134706648 0.14891361 0.164153345 0.353145247 0.313251484 0.373761142 0.360901466 0.353697346 0.361937093;0.040560018 0.038002633 0.027442326 0.031787511 0.041754224 0.052949472 0.069234367 0.076544576 0.068018698 0.072672711 0.076786118 0.083558949 0.187125943 0.178016219 0.170146036 0.17161249 0.186164469 0.173520858 0.355085849 0.334727787 0.374815352 0.357844418 0.41236388 0.426641669];
ydata=[11.32441 10.7704 10.22368 11.11952 13.16454 14.32423 22.9482 24.10042 22.79831 20.17602 23.25512 22.65642 28.94045 29.99842 28.94679 29.09359 28.96675 26.9409 14.83138 19.4787 10.70323 9.11525 11.29219 9.5724];
a=[10,200];
[x,resnorm]=LSQCURVEFIT(f,a,xdata,ydata)
y=28.83*xdata(1,./(0.13+xdata(1,+xdata(1,.^2./x(1)+x(2)*xdata(2,);
plot3(xdata(1,,xdata(2,,ydata,'*',xdata(1,,xdata(2,,y,'r.');
sum(sum((y-ydata).^2))/size(y,2)

运行结果显示
??? Attempt to execute SCRIPT lsqcurvefit as a function.

Error in ==> lsqcurvefit at 5
[x,resnorm]=LSQCURVEFIT(f,a,xdata,ydata) @beefly 返回小木虫查看更多