资源描述：
《数值分析作业-曲线拟合的最小二乘法》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、学号：专业：学院：一、实验课题：曲线拟合的最小二乘法从随机的数据中找岀其规律性，给出其近似表达式的问题，在生产实践和科学实验中大量存在，通常利用数据的最小二乘法求得拟合曲线。二、理论意义和实用价，如果己知函数f（x）在若干点xi（i=l,2,…，n）处的值yb便可根据插值原理来建立插值多项式作为f00的近似。但在科学实验和生产实践中，往往会遇到这样一种情况，即节点上的函数值并不是很精确的，这些函数值是由实验或观测得到的数据，不可避免地带有测量误差，如果耍求所得的近似函数曲线精确无误地通过所有的点（Xi，yi），就会使曲线保留着一些测试误差。当个别数据的误差较人时，插值效果显然是不理想的。此外
2、，由实验或观测提供的数据个数往往很多,如果用插值法,势必得到次数较高的插值多项式，这样计算起来很烦琐。所以我们设想:在大量的随机数据X（X：、X2、X3......Xn＞与Y（yi、y2、......yn）.从看似无规律的这两组离散数据中，找到一条一条曲线Y=F（x），使数据点均在离此曲线的上方或下方不远处,它既能反映数据的总体分布，又不至于出现局部较大的波动，更能反映被逼近函数的特性，使求得的逼近函数与已知函数从总体上來说其偏差按某种方法度量达到最小,这就是曲线拟合最小二乘法。在对给出的实验（或观测）数据作曲线拟合时，一般希望各实验（或观测）数据与拟合曲线的偏差的平方和最小,这就是最小二乘
3、原理。曲线拟合不要求曲线通过所有已知点，而是要求得到的近似函数能反映数据的基本关系。此外，由于实验或观测提供的数据个数往往很多，如果用插值法，势必得到次数较高的插值多项式，这样计算起来很麻烦，缺乏实用价值，所以从某些意义上来说，在解决实际问题的过程中，曲线拟合更具有实用价值。三、计算过程将给定数据作散点图，始图所示，选择形如5々）11^吻沾他山3）作为拟合曲线，这里0(x)=span(l,x,x2,x3)(4.8)设/(x)eC[a,b],O=span｛识0，…，(pn｝。C[a，b]，妒()，…，识,，是正交函数族，则ao-(/，％)-•衾al♦♦=(/，仍)•••K♦X.參=(/，^*)
4、《P(，识々)，‘<(又)=S“x)=E(/，DttoII(pkII2(Pk(4.9)均方误差II弋⑴II2=II/(幻—<(x)
5、
6、X)-S:(X)，/(x)))-(E0II(PkII2=II/II2-Z根据以上公式求出:ao=O.000a,=-0.0052,a2=0.2634a3=0.0178rToII(PkII2y所求的拟合曲线为SjXh-O.0052X+0.2634x2+0.0178x3t(分)0510152025303540455055Y*(104)01.272.162.863.443.874.154.374.514.584.024.64Si(X)0.01781.20872.164
7、62.91133.47453.88004.15364.32114.40824.44074.44444.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+ai(x)+a2(x2)求出aQ=-0.0024a,=0.2037,a2=0.2305所求的S2(X)=-O.0024+0.2037x+0.2305x
8、2tm0510152025303540455055Y*(104)01.272.162.863.443.874.154.374.514.584.024.64秘)0.23051.18942.02932.75023.35213.83494.19874.44354.56934.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.
9、19690.1661Ul曲线拟合合破量与时间的关系五、结构程序设计x=0:5:55y=[01.272.162.863.443.874.154.374.514.584.024.64];xnum=length(x);plot(x,y/o’)，holdon;pi=polyfit(x,y,3)pj=polyfit(x,y,2)newx=0:0.1:55;newyi=polyval(pi,newx);newyj=pol
