作者:淋殁 | 来源:互联网 | 2023-07-30 08:07
%自定义拟合函数f(t)a*cos(k*t)*exp(w*t)clc,clearsymstx[0;0.4;1.2;2;2.8;3.6;4.4;5.2;6;7.2;8;9.2;
%自定义拟合函数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-');
结果:
cfun =
General model:
cfun(t) = a*cos(k*t)*exp(w*t)
Coefficients (with 95% confidence bounds):
a = 0.9987 (0.9836, 1.014)
k = 1.001 (0.9958, 1.006)
w = -0.2066 (-0.2131, -0.2002)