《数学函数逼近》课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《《数学函数逼近》课件.ppt》由用户(晟晟文业)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学函数逼近 数学 函数 逼近 课件
- 资源描述:
-
1、EAST CHINA INSTITUTE OF TECHNOLOGY理学院 2.1 数据拟合数据拟合(最小二乘法最小二乘法)实例:考察某种纤维的强度与其拉伸倍数的关系,下表是实际测定的24个纤维样品的强度与相应的拉伸倍数 通过观察或测量到一组离散数据序列,当所得数据比较准确的时候,可构造插值函数,构造的原则是要求插值的函数通过这些插值节点。但通常所测得的数据有误差,如果数据系列无法同时满足某特定函数,即插值无法满足要求,在这种情况下只能通过逼近函数最优的靠近样点,即数据数据拟合。拟合。EAST CHINA INSTITUTE OF TECHNOLOGY理学院编 号 拉伸倍数 强 度编 号 拉伸
2、倍数 强 度11.91.41355.5221.3145.2532.11.81565.542.52.5166.36.452.72.8176.5662.72.5187.15.373.531986.583.52.72087944218.98.51043.52298114.54.2239.58.1124.63.524108.1iiyxiiyxEAST CHINA INSTITUTE OF TECHNOLOGY理学院1234567891012345678912345678910123456789纤维强度随拉伸倍数增加而增加系要关系应是线性关的主与拉伸倍数因此可以认为强度xy并且24个点大致分布在一条直线
3、附近xxy10)(为待定参数其中10,-(1)EAST CHINA INSTITUTE OF TECHNOLOGY理学院01 ()()(,)iiy xxx y我们希望与所有的数据点 样本点越接近越好。必须找到一种度量标准来衡量什么曲线最接近所有数据点。一、最小二乘法的基本概念iiiyxy)(令一般使用mii0222在回归分析中称为残差miiiyxy02)(准偏离程度大小的度量标与数据点作为衡量),()(iiyxxy称为平方误差EAST CHINA INSTITUTE OF TECHNOLOGY理学院在回归分析中称为残差平方和。mii0222miiiyxy02)(注意(1)式是一条直线。,x y
4、但的关系并不一定是线性关系。因此将问题一般化:EAST CHINA INSTITUTE OF TECHNOLOGY理学院)(,xSyyx的关系为设来自函数类其中)(xS来自线性函数类中如)()1(xy为给定的一组数据设),1,0)(,(miyxii),1,0)(nixi的基函数为设函数类nm一一般般要要求求即生成的函数集是由也称,),1,0)(nixi)(,),(),(10 xxxspannmii0222miiiyxS02)(仍然定义平方误差njjjxaxS0)()(EAST CHINA INSTITUTE OF TECHNOLOGY理学院我们选取的度量标准是:)(*xS中选取一个函数在函数类
5、njjjxaxS0*)()(*)(*)(*)(*1100 xaxaxann22*使得miiiyxS02)(*(miiixSyxS02)()(min22)(minxS0()()njjjS xax其中为中的任意函数。-(2)-(3)EAST CHINA INSTITUTE OF TECHNOLOGY理学院),1,0(,)(njaxSj如何求拟合系数后在确定了拟合函数呢?满足拟合条件使得)3()()(*0*njjjxaxS*0(3)*()()njjjSxax称满足条件的求函数的方法为数据拟合的最小二乘法.*0*()()njjjSxax为最小二乘解.0()(),(0,1,)njjjjS xaxajn为
6、拟合函数为拟合系数.22*称为最小二乘解的平方误差.EAST CHINA INSTITUTE OF TECHNOLOGY理学院 miinjijjyxa020)(miiiyxS02)(二、法方程组22njjjxaxS0)()(由(0,1,)jajn为拟合系数的函数.可知因此可假设),(10naaa miinjijjyxa020)(因此求最小二乘解转化为EAST CHINA INSTITUTE OF TECHNOLOGY理学院的问题点极小值的最小值求*,*,*,)(),(1010nnaaaaaa由多元函数取极值的必要条件0),(10knaaaank,1,0)()(200ikmiinjijjxyxa
7、ka0得即miikimiiknjijjxyxxa000)()()(0)()()(00 ikmiinjikijjxyxxa200()mnjjiiijaxy EAST CHINA INSTITUTE OF TECHNOLOGY理学院miikimiiknjijjxyxxa000)()()(miikinjjikmiijxyaxx000)()()(nk,1,0-(4)miikiikmiinnikmiiikmiixyxxaxxaxxa00011000)()()()()()()(nk,1,0即EAST CHINA INSTITUTE OF TECHNOLOGY理学院元线性方程组的是一个关于显然1,)4(10
8、naaan引入记号)(,),(),(10mrrrxxxr),(10myyyf)()(),(0ijmiikjkxx则由内积的概念可知imiikkyxf0)(),(-(5)-(6),(jk),(kj显然内积满足交换律EAST CHINA INSTITUTE OF TECHNOLOGY理学院方程组(4)便可化为),(),(),(),(1100faaaknknkknk,1,0-(7)的线性方程组常数项为这是一个系数为),(),(fkjk将其表示成矩阵形式naaa10),(),(),(10fffn),(),(),(01000n),(),(),(11101n),(),(),(10nnnn-(8)EAST
9、CHINA INSTITUTE OF TECHNOLOGY理学院上的法方程组在点式为函数序列称mnxxxxxx,)(,),(),()8(1010的基为函数类由于)(,),(),(10 xxxn必然线性无关因此)(,),(),(10 xxxn并且其系数矩阵为对称阵所以法方程组的系数矩阵非奇异,即0),det(nnji根据Cramer法则,法方程组有唯一解*,*,*,1100nnaaaaaaEAST CHINA INSTITUTE OF TECHNOLOGY理学院*),*,*,(10naaa miinjijjyxa020)(),(10naaa即是的最小值22*miiiyxS02)(*(miiixS
展开阅读全文