数值分析matab版第三章.docVIP

实验报告三 题目：曲线拟合的最小二乘法 摘要：在实验中往往要从一组实验数据（xi,yi)(i=0,1,2,3，...，m）中寻找自变量x与因变量y之间的函数关系y=F（x）。由于观测数据往往不准确，因为要求y=F（x）经过所有点（xi，yi），而只要求在给定点xi上误差按某种标准最小。 数学原理： 所谓的曲线拟合的最小二乘法就是使拟合多项式在各节点处的偏差P(xi）-yi的平方和达到最小。 习题22. function c =myployfit(x,y,m) n =length(x); b =zeros(1,m+1); f =zeros(n,m+1); for k=1:m+1 f(:,k)=x.^(k-1); end a=f*f; b=f*y; c=a\b; c=flipud(c); x=19:6:44; y=[19.0,32.3,49.0,73.3,97.8]; m=1; C=mypolyfit(x,y,m) C = 3.309999999999997 -48.329999999999913 plot(x,y,ok,LineWidth,2); grid on; title(最小二乘一次拟合效果图); x0=0:2:44; y0=C(1).*x0+C(2); hold on,plot