主要求解過程:
執行過程:
tdata = linspace(100,1000,10);
cdata = 1e-05.*[454 499 535 565 590 610 626 639 650 659];
x0 = [0.2,0.05,0.05];
x = lsqcurvefit('curvefun',x0,tdata,cdata)
f = curvefun(x,tdata)
plot(tdata,cdata,'o',tdata,f,'r-')
函數文件curvefun.m
function f = curvefun(x,tdata)
f = x(1)+x(2)*exp(-0.02*x(3)*tdata);
結果: