[求助]最小二乘曲线拟合中一个小问题!!(有程序)

haohaoljf · 发表于 2006-6-1 16:59

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

曲线拟合的程序如下：（摘至《精通matlab6.5》） <A>% </A>被拟合的原始数据 x=0:0.1:1;y=[2.1,2.3,2.5,2.9,3.2,3.3,3.8,4.1,4.9,5.4,5.8]; dy=0.15; % 原数据 y 的标准差 for n=1:6 % 依次用 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; % 计算不同阶次的 <IMG src="http://www.ddvip.net/machine/matlab/index1/img_matlab/30/30_clip_image002_0000.gif">量。 end Q=1-chi2cdf(chi2,freedom); % 用于判断拟合良好度 % 适当度的图示 subplot(1,2,1),plot(1:6,abs(chi2-freedom),'b') xlabel(' 阶次 '),title('chi 2 与自由度 ') subplot(1,2,2),plot(1:6,Q,'r',1:6,ones(1,6)*0.5) xlabel(' 阶次 '),title('Q 与 0.5 线 ') 
<A>% </A>适当的三阶多项式拟合情况图示 clf,plot(x,y,'b+');axis([0,1,1,6]);hold on errorbar(x,YE{3},D{3},'r');hold off title(' 较适当的三阶拟合 ') text(0.1,5.5,['chi2=' num2str(chi2(3)) '~' int2str(freedom(3))]) text(0.1,5,['freedom=' int2str(freedom(3))]) text(0.6,1.7,['Q=' num2str(Q(3)) '~0.5']) 哪个高手能帮我解决一下dy=0.15是怎么的来的 原数据y的标准差根据std(y)好象不是这么多啊! 我就郁闷啊,毕业设计就要结束了,跪请哪位高手帮忙解决一下!! [em06][em06][em01][em01][em01]

ericlin · 发表于 2006-6-2 12:16

也有可能是作者笔误吧？

haohaoljf · 发表于 2006-6-3 16:04

但运行没有错误的！！

账号		自动登录	找回密码
密码			我要加入

[求助]最小二乘曲线拟合中一个小问题!!(有程序)

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（ericlin）也有可能是作者笔误吧？