最小二乘法拟合指数函数（Matlab编程）,着急,拟合函数形式为：y=a×exp（-bt）；数据比较少：y=[50 40 30 20 10 ]; t=[55 63 73 100 121],用matlab编程程序该怎么写啊,需要求出未知系数a和b的值
来源：学生作业帮助网 编辑：六六作业网 时间：2024/12/25 01:31:49
最小二乘法拟合指数函数（Matlab编程）,着急,拟合函数形式为：y=a×exp（-bt）；数据比较少：y=[50 40 30 20 10 ]; t=[55 63 73 100 121],用matlab编程程序该怎么写啊,需要求出未知系数a和b的值
最小二乘法拟合指数函数（Matlab编程）,着急,
拟合函数形式为：y=a×exp（-bt）；数据比较少：y=[50 40 30 20 10 ]; t=[55 63 73 100 121],用matlab编程程序该怎么写啊,
需要求出未知系数a和b的值
最小二乘法拟合指数函数（Matlab编程）,着急,拟合函数形式为：y=a×exp（-bt）；数据比较少：y=[50 40 30 20 10 ]; t=[55 63 73 100 121],用matlab编程程序该怎么写啊,需要求出未知系数a和b的值
%方法一
y=[50 40 30 20 10 ]; t=[55 63 73 100 121];
yp=log(y);
p = polyfit(t,yp,1);
b=-p(1)
a=exp(p(2))
yf=a*exp(-b*t);
yf-y
plot(t,y,'r+',t,yf,'b-')
legend('原始点','拟合线')
%方法二
%% Fit:'exp1'.
[xData,yData] = prepareCurveData( t,y );
% Set up fittype and options.
ft = fittype( 'exp1' );
opts = fitoptions( ft );
opts.StartPoint = [145.2 -0.3];
% Fit model to data.
[fitresult,gof] = fit( xData,yData,ft,opts );
% Plot fit with data.
figure( 'Name','untitled fit 1' );
h = plot( fitresult,xData,yData );
legend( h,'y vs.t','exp1','Location','NorthEast' );
% Label axes
xlabel( 't' );
ylabel( 'y' );
fitresult
gof
%方法三
y=[50 40 30 20 10 ]';
yp=log(y);
t=[55 63 73 100 121]';
tl=ones(size(t));
t1=[tl t];
p=t1\yp;
b=-p(2)
a=exp(p(1))
yf=a*exp(-b*t);
yf-y
plot(t,y,'r+',t,yf,'b-')
legend('原始点','拟合线')
%方法四
regress
clc;
clear all;
y=[50 40 30 20 10 ]; x=[55 63 73 100 121];
a=[145.2345667367 -0.02061401470466 ];%设置初值，很关键。
f=@(a,x)a(1)*exp(-a(2)*x);
nlinfit(x,y,f,a)
%可以利用toolbox里面的拟合函数来确定初值...
全部展开
clc;
clear all;
y=[50 40 30 20 10 ]; x=[55 63 73 100 121];
a=[145.2345667367 -0.02061401470466 ];%设置初值，很关键。
f=@(a,x)a(1)*exp(-a(2)*x);
nlinfit(x,y,f,a)
%可以利用toolbox里面的拟合函数来确定初值：
（先在窗口输入x,y值）
command window输入cftool，然后再data里面选x，y数值。create data set。
Fitting里面选择new fit, type of fit 里面选择exponential,单击apply。
在result里面就可以看到a，b的值。（还有不同的组合）
然后你点击File，Generate M-file.
打开你生成的M-File，里面有st_就是matlab设置的初值。
希望可以帮到你。
收起