作业帮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/05/09 11:40:18
![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,](/uploads/image/z/5924745-9-5.jpg?t=matlab%E6%9C%80%E5%B0%8F%E4%BA%8C%E4%B9%98%E6%B3%95%E6%9B%B2%E7%BA%BF%E6%8B%9F%E5%90%88+m%E6%96%87%E4%BB%B6%E9%94%99%E8%AF%AF%E5%9C%A8%E5%93%AAx%3D10000%3A1000%3A23000%3By%3D%5B24.0%2C26.5%2C29.8%2C32.4%2C34.7%2C41.1%2C42.8%2C44.6%2C47.3%2C65.8%2C87.5%2C137.8%2C174.2%5D%3Bdy%3D1.5%3Bfor+n%3D1%3A6%5Ba%2CS%5D%3Dpolyfit%28x%2Cy%2Cn%29%3BA%7Bn%7D%3Da%3Bda%3Ddy%2Asqrt%28diag%28inv%28S.R%27%2AS.R%29%29%29%3BDA%7Bn%7D%3D+da%27%3Bfreedom%28n%29%3DS.df%3B%5Bye%2C)
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没有越界吧