首先根據計算t是弧度制的,這樣的話有以下擬合結果:
t=[0,0.2,0.4,1.3,5.6]';
y=[0.8,0.7,0.6,0.5,0.2]';
f = @(c,t)c(1) + c(2).*exp(-t) + c(3).*cos(t);
c = [1;1;1]; %c是系數向量,對其賦任意初值
c = lsqcurvefit(f,c,t,y); ? %最小二乘擬合函數
y1 = f(c,t);
plot(t,y,t,y1,'o')
c =
0.3836 %c1
0.6606 %c2
-0.2396 %c3