[MATLAB基础] 曲线拟合求助

[lincsu]

slgu兄可否帮忙看看，我看到你回答的上面一个醉酒酒精问题，我这个类似，但是我怎么知道初始值

[slgu]

数据文件打不开，我就很难帮你了。
另外，第一个系数以coeff(1)^2出现似乎不合常规。

[lincsu]

您好，数据文件下载后，解压缩，
然后会发现一个txt文件，吧txt文件改成.mat文件就能读出来，是一个matlab数据文件，可否给我您的邮箱什么？我传给您

非常感谢

[lincsu]

请查看这个附件，上传不了excel文件，附件是txt，吧后缀名字改成xlxs的excel文件就可以看到数据
非常感谢

[slgu]

[email protected]
最好用excel格式。

[slgu]

你提供的模型可能不恰当，请用该模型进行拟合：
fx=@(b,x)b(1)*exp(b(2)*x+b(3)*x.^2+b(4)*x.^3+b(5)*x.^4);
请运行下列小程序：
clear,clc
x=[...]';%清输入x数据
y=[...]';%请输入y数据
b=[2.5540354806*10^(-41), 2653.32513908, -30459.7660992, 153943.375849, -292083.699481]
fx=@(b,x)b(1)*exp(b(2)*x+b(3)*x.^2+b(4)*x.^3+b(5)*x.^4);
b=lsqcurvefit(fx,b,x,y)
b=nlinfit(x,y,fx,b)
n=length(x)
SSy=var(y)*(n-1)
y1=fx(b,x);
q=(y-y1)'*(y-y1)
rsquare=(SSy-q)/SSy
figure(1),clf
plot(x,y,'ko')
hold on
plot(x,y1,'linewidth',2)
legend('data','fit','location','best')
axis tight