拟合完成。

clear ,clc
data=load('case.txt');%将数据放于工作空间,改名为case.txt
x=data(:,1);
y=data(:,3);
fx=@(a,x)a(4)*exp(-a(1)*x).*sin(a(2)*x+a(3))+a(5);%拟合目标函数
a0=[-1 1 3 0.1 4];%初始点
[ax,resnorm] = lsqcurvefit(fx,a0,x,y);%最小二乘法拟合
figure
plot(x,y,x,fx(ax,x),'r')
legend('原数据','拟合')
结果:
ax =
-0.1320 3.4604 -43.4173 -0.0390 4.0066
即:
f(x)=4.00658 - 0.0390409*exp(0.132015*x)*sin(3.4604*x - 43.4173)