Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#8 |
高级会员
注册日期: 2007-06-24
年龄: 70
帖子: 188
声望力: 21 ![]() |
![]()
你提供的模型可能不恰当,请用该模型进行拟合:
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 |
![]() |
![]() |