matlab最小二乘法曲线拟合 m文件错误在哪x=10000:1000:23000;y=[24.0,26.5,29.8,32.4,34.7,41.1,42.8,44.6,47.3,65.8,87.5,137.8,174.2];dy=1.5;for n=1:6[a,S]=polyfit(x,y,n);A{n}=a;da=dy*sqrt(diag(inv(S.R'*S.R)));DA{n}= da';freedom(n)=S.df;[ye,
来源：学生作业帮助网 编辑：六六作业网 时间：2024/12/26 12:46:41
matlab最小二乘法曲线拟合 m文件错误在哪x=10000:1000:23000;y=[24.0,26.5,29.8,32.4,34.7,41.1,42.8,44.6,47.3,65.8,87.5,137.8,174.2];dy=1.5;for n=1:6[a,S]=polyfit(x,y,n);A{n}=a;da=dy*sqrt(diag(inv(S.R'*S.R)));DA{n}= da';freedom(n)=S.df;[ye,
matlab最小二乘法曲线拟合 m文件错误在哪
x=10000:1000:23000;
y=[24.0,26.5,29.8,32.4,34.7,41.1,42.8,44.6,47.3,65.8,87.5,137.8,174.2];
dy=1.5;
for n=1:6
[a,S]=polyfit(x,y,n);
A{n}=a;
da=dy*sqrt(diag(inv(S.R'*S.R)));
DA{n}= da';
freedom(n)=S.df;
[ye,delta]=polyval(a,x,S);
YE{n}=ye;
D{n}=delta;
chi2(n)=sum((y-ye).^2)/dy/dy;
end
Q=1-chi2cdf(chi2,freedom);
clf,shg
subplot(1,2,1),plot(1:6,abs(chi2-freedom),'b')
xlabel('阶次'),title('chi2与自由度')
subpolt(1,2,2),plot(1:6,Q,'r',1:6,ones(1,6)*0.5)
xlabel('阶次'),title('Q与0.5线')
nod=input('根据图形选择适当的阶次（请输入数值）');
clf,shg,
plot(x,y,'kx');
xlabel('x'),ylabel('y');
axis([8000,23000,20.0,174.2]);hold on
error(x,YE{nod},D{nod,'r'});hold off
title('较适合阶次的拟合')
text(10000,150.0,['chi2='num2str(chi2(nod))'~1'int2str(freedom(nod))])
text(10000,140.0.['freedom='int2str(freedom(nod))])
text(2000,40.0,['Q='num2str(Q(nod))'~0.5'])
disp()
disp('拟合多项式系数'),diap(A{nod})
disp('拟合多项式离差'),disp(DA{nod})
第28行出错了,
matlab最小二乘法曲线拟合 m文件错误在哪x=10000:1000:23000;y=[24.0,26.5,29.8,32.4,34.7,41.1,42.8,44.6,47.3,65.8,87.5,137.8,174.2];dy=1.5;for n=1:6[a,S]=polyfit(x,y,n);A{n}=a;da=dy*sqrt(diag(inv(S.R'*S.R)));DA{n}= da';freedom(n)=S.df;[ye,
你程序里的28行是这个吗text(10000,150.0,['chi2='num2str(chi2(nod))'~1'int2str(freedom(nod))])
看不出来有没有空格啊text(10000,150.0,['chi2=' num2str(chi2(nod)) '~1' int2str(freedom(nod))])
报的什么错误,输入的nod没有越界吧