实验名称 拟合法3.doc

  1. 1、本文档共14页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验名称 拟合法3

探索实验6 曲线拟合法 实验目的 了解最小二乘法与曲线拟合问题及用法,理解并掌握多项式拟合及线性模型曲线拟合算法的构造和计算,学习用计算机求近似函数的一些科学计算方法和简单的编程技术。 二、概念与结论 1. 最小二乘法则: 设f为函数y=f(x)的一组实验点(xk f(xk)), k=0,1,2,…,n 组成的,在一个函数集合P中选择一个函数 (*(x) 作为函数f(x)的近似表达式,使 这里(k 0为已知的一组数值, yk= f(xk)。函数集合P称为拟合函数类。 2.线性模型拟合: 如果用最小二乘法确定函数 (*(x) 作为函数f(x)的近似函数时,使用的拟合函数类P是由已知函数(0(x)、(1(x)、…、(m(x)的线性组合构成的线性空间 P={((x)| ((x)= a0(0(x)+ a1(1(x)+ a2(2(x)+ …+ am(m(x)}=Span[(0 ,(1 , …,(m ] 中所求出的近似函数。这里已知函数(0(x)、(1(x)、…、(m(x)称为基函数,它们彼此是线性无关的。 3.多项式拟合 在函数f(x) 的线性拟合中,如果取基函数 (j(x)=x j ,即,拟合函数类是 M =Span[1,x , x2,…,xm] 构成的m次多项式空间时得到的近似函数。这个拟合函数称为最小二乘多项式。 4.Haar条件 设(0(x)、(1(x)、…、(m(x)是闭区间[a,b]上的连续函数,如果它们的任意线性组合在点集{xj}={xj, j=0,1,2,…,n},n(m 至多有m个不同的零点 ,则称(0(x)、(1(x)、…、(m(x)是在点集{xj}上满足Haar条件。 三、程序中Mathematica语句解释: 1.Plot[f[t],{t,a,b}, PlotLabel-{m****}] 画出带有变量m的值和数字串“****”表识的图形 2. v=Plot[f[t],{t,a,b},DisplayFunction-Identity] 将函数f(t)在区间[a,b]上的图形文件存在变量v中,但不画出该图形。 3. Show[v, DisplayFunction-$DisplayFunction] 将存在变量v中的图形文件画出。 4.p={f1[x],f2[x],…,fn[x]}/.x-a 将一组函数值{f1(a),f2(a),…,fn(a)}赋值给变量p 四、方法与程序 曲线拟合也是函数逼近的重要方法之一, 它是求近似函数的一种方法,有着广泛的应用。插值法关心插值函数与被插函数在插值节点上函数值相同,曲线拟合关心的是在数据点上的误差平方和最小,它可以克服测量数据误差对确定近似函数的影响。曲线拟合法也有很多种,其中常用的有多项式拟合和线性模型拟合。 多项式拟合 如果拟合函数函数类是一个多项式集合,则称为多项式拟合,求出的拟合多项式称为最小二乘多项式,多项式拟合是最简单的拟合形式。 1.1 多项式拟合的构造过程: 令所求的拟合函数((x)= a0+ a1x+ a2x2+ …+ amxm 为确定系数a0,a1, …,am ,考虑平方和 于是求解归结为 m+1元函数S(a0,a1,…,a m)的极值问题,由多元函数极值的必要条件,有 整理后,得 这是关于a0,a1,…,a m的线性方程组,引进函数内积符号: 于是得到如下线性方程组 这个线性方程组称为法方程组或正规方程组,它的系数矩阵是对称矩阵,这个方程组是正定方程组,它有唯一解。设此解为a0*,a1*,…,a m*,则函数 (*(x)= a0*+ a1*x+ a2*x2+ …+ am*xm 既为所求的多项式拟合。 注意到正规方程组的系数具有规律,引入符号 于是正规方程可以写为更简单的形式 1.2 多项式拟合算法: 输入n+1个拟合点: (xi, yi),i=0,1,…,n 根据散点图确定拟合多项式的次数m 计算相应正规线性方程组的系数和右端项 解正规正规线性方程组,得解:a0*,a1*,…,a m* 写出拟合多项式(*(x)= a0*+ a1*x+ a2*x2+ …+ am*xm 1.3 求m次多项式拟合程序: Clear[xi,xx,yi]; xi=Input[xi=] yi=Input[yi=] n=Length[xi]; h=ListPlot[Table[{xi[[i]],yi[[i]]},{i,1,n}],PlotStyle-PointSize[0.04]] m=Input[多项式次数m=] s=Table[Sum[xi[[k]]^i,{k,1,n}],{i,0,2m}]; a=Table[s[[i+j-1]],{i,1

文档评论(0)

wyjy + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档