%自定义拟合函数f(t)=a*cos(k*t)*exp(w*t)
clc,clear
syms t
x=[0;0.4;1.2;2;2.8;3.6;4.4;5.2;6;7.2;8;9.2;10.4;11.6;12.4;13.6;14.4;15];%列向量
y=[1;0.85;0.29;-0.27;-0.53;-0.4;-0.12;0.17;0.28;0.15;-0.03;-0.15;-0.071;0.059;0.08;0.032;-0.015;-0.02];
f=fittype('a*cos(k*t)*exp(w*t)','independent','t','coefficients',{'a','k','w'}); %fittype是自定义拟合函数
cfun=fit(x,y,f) %根据自定义拟合函数f来拟合数据x,y
xi=0:0.1:20;
yi=cfun(xi);
plot(x,y,'r*',xi,yi,'b-');
通过改变fittype里的函数表达式来完成自定义拟合,可以通过cfun.a来查看和输出拟合系数a