数值分析作业-曲线拟合的最小二乘法

1、数值分析实验报告 课题八： 曲线拟合的最小二乘法 姓 名: 学 号： 专 业： 学 院： 一、实验课题:曲线拟合的最小二乘法 从随机的数据中找出其规律性，给出其近似表达式的问题，在生产实践和科学实验中大量存在，通常利用数据的最小二乘法求得拟合曲线。t(分）0510152025303540455055Y*(10-4）01.272.162.863.443.874.154.374.514.584.024.在某冶炼过程中，根据统计数据的含碳量与时间关系，试求含碳量与时间t的拟合曲线。二、理论意义和实用价值。 如果已知函数f(x)在若干点xi(i=1,2,n)处 的值yi,便可根据插值原理来建立插值多项

2、式作为f(x)的近似。但在科学实验和生产实践中，往往会遇到这样一种情况，即节点上的函数值并不是很精确的，这些函数值是由实验或观测得到的数据，不可避免地带有测量误差，如果要求所得的近似函数曲线精确无误地通过所有的点(xi,yi),就会使曲线保留着一些测试误差。当个别数据的误差较大时,插值效果显然是不理想的。此外,由实验或观测提供的数据个数往往很多,如果用插值法,势必得到次数较高的插值多项式，这样计算起来很烦琐。 所以我们设想:在大量的随机数据X（X1、X2、X3Xn)与Y(y1、y2、yn),从看似无规律的这两组离散数据中，找到一条一条曲线Y=F(x),使数据点均在离此曲线的上方或下方不远处,它

3、既能反映数据的总体分布,又不至于出现局部较大的波动,更能反映被逼近函数的特性,使求得的逼近函数与已知函数从总体上来说其偏差按某种方法度量达到最小,这就是曲线拟合最小二乘法。 在对给出的实验(或观测)数据作曲线拟合时,一般希望各实验(或观测)数据与拟合曲线的偏差的平方和最小,这就是最小二乘原理。 曲线拟合不要求曲线通过所有已知点,而是要求得到的近似函数能反映数据的基本关系。此外，由于实验或观测提供的数据个数往往很多，如果用插值法，势必得到次数较高的插值多项式，这样计算起来很麻烦，缺乏实用价值，所以从某些意义上来说,在解决实际问题的过程中，曲线拟合更具有实用价值。三、计算过程 将给定数据作散点图，

4、始图所示，选择形如S1(X)=a1(x)+a2(x2)+a3(x3)作为拟合曲线，这里x=span(1,x,x2,x3)根据以上公式求出：a0=0.000 a1=-0.0052 ，a2=0.2634 a3=0.0178 所求的拟合曲线为S1(X)=-0.0052X+0.2634x2+0.0178x3t(分）0510152025303540455055Y*(10-4）01.272.162.863.443.874.154.374.514.584.024.64S1(X)0.01781.20872.16462.91133.47453.88004.15364.32114.40824.44074.4444

5、4.4450逼近误差0.0178-0.06130.00460.05130.03450.01000.0036-0.0489-0.1018-0.13930.4244-0.1950平方误差0.00030.00380.00000.00260.00120.00010.00000.000240.01040.01940.18010.0380 若取S2(X)=a0+a1(x)+a2(x2) 求出a0=-0.0024 a1=0.2037 ，a2=0.2305 所求的S2(X)=-0.0024+0.2037x+0.2305x2t(分）0510152025303540455055Y*(10-4）01.272.162

6、.863.443.874.154.374.514.584.024.64S1(X)0.23051.18942.02932.75023.35213.8349 4.19874.44354.5693 4.57604.46374.2324逼近误差0.2305-0.0806-0.1307-0.1098-0.0879-0.03510.04870.07350.0593-0.0040.4437-0.4076平方误差0.05310.00650.01710.01210.00770.00120.00240.00540.00350.00000.19690.1661四、曲线拟合图五、结构程序设计 x=0:5:55y=0

7、1.27 2.16 2.86 3.44 3.87 4.15 4.37 4.51 4.58 4.02 4.64;xnum=length(x);plot(x,y,o),hold on;pi=polyfit(x,y,3)pj=polyfit(x,y,2)newx=0:0.1:55;newyi=polyval(pi,newx);newyj=polyval(pj,newx);yi=polyval(pi,x)yj=polyval(pj,x)erri=yi-yerrj=yj-yplot(x,erri,.r),hold on;plot(x,errj,*b),hold on;S_erri=sum(erri.2)

8、S_errj=sum(errj.2)plot(newx,newyi,r,newx,newyj,b),grid off,hold off; x = 0 5 10 15 20 25 30 35 40 45 50 55pi = 0.0000 -0.0052 0.2634 0.0178pj = -0.0024 0.2037 0.2305yi = Columns 1 through 8 0.0178 1.2087 2.1646 2.9113 3.4745 3.8800 4.1536 4.3211 Columns 9 through 12 4.4082 4.4407 4.4444 4.4450yj = C

9、olumns 1 through 8 0.2305 1.1894 2.0293 2.7502 3.3521 3.8349 4.1987 4.4435 Columns 9 through 12 4.5693 4.5760 4.4637 4.2324erri = Columns 1 through 8 0.0178 -0.0613 0.0046 0.0513 0.0345 0.0100 0.0036 -0.0489 Columns 9 through 12 -0.1018 -0.1393 0.4244 -0.1950errj = Columns 1 through 8 0.2305 -0.0806 -0.1307 -0.1098 -0.0879 -0.0351 0.0487 0.0735 Columns 9 through 12 0.0593 -0.0040 0.4437 -0.4076S_erri = 0.2583S_errj = 0.4720六、结果讨论和分析 最小二乘意义下的曲线拟合,由于方程比较简单,实际中应用