特定函数的数据拟合

%指定函数拟合
clear
clc
syms t %再声明一个自变量
x = [1;2;3;4;5;6;7];
y = [88.3;88.0;70.0;65.4;65.0;60.0;50.5];
plot(x,y,’ro’); % 画出散点图
%假设函数是f=acos(kt)e^wt
f = fittype(‘a*cos(k*t)*exp(w*t)’,’independent’,’t’,’coefficients’,{‘a’,’k’,’w’});
cfun = fit(x,y,f) %显示拟合函数
xi = 0:.1:20;
yi = cfun(xi);
plot(x,y,’r*’,xi,yi,’b-‘);