函数的最佳逼近课件.ppt
《函数的最佳逼近课件.ppt》由会员分享,可在线阅读,更多相关《函数的最佳逼近课件.ppt(75页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、关于函数的最佳逼近现在学习的是第1页,共75页01.04.20231第三章 函数的最佳逼近1 最佳逼近问题最佳逼近问题一、函数的逼近方法一、函数的逼近方法关于函数的关于函数的n次次多项式多项式逼近方法逼近方法,已知有下面的几种:,已知有下面的几种:1.Taylor展式展式如果如果误差为误差为现在学习的是第2页,共75页01.04.20232第三章 函数的最佳逼近2.2.插值多项式插值多项式 同为同为n 次多项式,哪一个逼近效果更好呢?这时可以建次多项式,哪一个逼近效果更好呢?这时可以建立一个立一个度量标准度量标准来进行度量。在所建立的度量标准之下,来进行度量。在所建立的度量标准之下,就可以给出
2、就可以给出最佳的最佳的n 次逼近多项式次逼近多项式。注:注:除了用除了用多项式多项式来逼近一个函数来逼近一个函数 f(x),也可以用其它具也可以用其它具有某种有某种共同特征的函数共同特征的函数来逼近来逼近 f(x),并求出其相应的最佳,并求出其相应的最佳逼近。逼近。例如,例如,现在学习的是第3页,共75页01.04.20233第三章 函数的最佳逼近3.3.最佳逼近问题最佳逼近问题 给定给定函数空间函数空间X X 中的一个中的一个子集合子集合 ,对于对于某一已知函数某一已知函数f(x)X,在在 中寻求一个函数中寻求一个函数p(x)作为作为函数函数f(x)关于某个关于某个度量标准度量标准下的最佳逼
3、近下的最佳逼近函数函数,称之为称之为最佳逼近问题最佳逼近问题。X 本章我们主要考虑本章我们主要考虑连续函数空间连续函数空间X=Ca,b上的最佳逼近问题,这时的子集上的最佳逼近问题,这时的子集合合可以取为由具有可以取为由具有某种共同特征某种共同特征的函数的函数组成,例如组成,例如多项式函数多项式函数、三角函数三角函数、指指数函数数函数、分式有理函数分式有理函数等。等。同时,还需要给出连续函数空间同时,还需要给出连续函数空间上的一个上的一个度量标准度量标准,下面先通过,下面先通过内积内积给出给出平方范数平方范数。p(x)从总体上更能反映从总体上更能反映f(x)的特性或总体上其偏差按某种度量达到最小
4、的特性或总体上其偏差按某种度量达到最小现在学习的是第4页,共75页01.04.20234第三章 函数的最佳逼近二、连续函数的平方范数二、连续函数的平方范数 已知所有连续函数构成的集合已知所有连续函数构成的集合Ca,b是一个是一个线性空间线性空间,对于,对于Ca,b中的任意函数中的任意函数 f(x)、g(x),定义实数定义实数可以证明此实数满足性质:可以证明此实数满足性质:这时,称这时,称(f,g)为为 f(x)与与 g(x)的的内积内积。(1).(f,g)=(g,f);(2).(f,g)=(f,g),R;(3).(f+g,h)=(f,h)+(g,h);(4).(f,f)0,当且仅当当且仅当 f
5、=0 时时(f,f)=0 现在学习的是第5页,共75页01.04.20235第三章 函数的最佳逼近为函数为函数 f(x)的的平方(欧氏)范数平方(欧氏)范数,且满足以下性质:,且满足以下性质:给出了函数的给出了函数的范数范数,便给出了函数的一个,便给出了函数的一个度量标准度量标准,在此度,在此度量标准之下,就可以找出量标准之下,就可以找出 f(x)在不同函数类中的最佳逼近。下在不同函数类中的最佳逼近。下面就来考虑面就来考虑这一最佳逼近问题的解决这一最佳逼近问题的解决。并称并称(3.1)(1)f2 0,f2=0,当且仅当当且仅当 f=0 ;(2)c f2=|c|f2;(3)f+g2 f2+g2;
6、无穷范数无穷范数现在学习的是第6页,共75页01.04.20236第三章 函数的最佳逼近柯西柯西施瓦施瓦茨不等式茨不等式现在学习的是第7页,共75页01.04.20237第三章 函数的最佳逼近基函数基函数2 函数的最佳平方逼近函数的最佳平方逼近一、公式的推导一、公式的推导 对于连续函数空间对于连续函数空间 Ca,b 中的元素中的元素 f(x)及其及其子空间子空间所谓所谓 f(x)在在 中的中的最佳平方逼近最佳平方逼近最佳平方逼近最佳平方逼近,就是存在,就是存在使得对于一切使得对于一切都有:都有:广义多广义多项式项式有限维有限维现在学习的是第8页,共75页01.04.20238第三章 函数的最佳
7、逼近不等式不等式 说明说明所求的所求的满足等式:满足等式:其中其中(3.2)由于由于pn*(x)是由其系数是由其系数c0*,c1*,cn*唯一确定的,因此,只要我唯一确定的,因此,只要我们求出了满足们求出了满足(3.2)的的 c0*,c1*,cn*,就可以求出就可以求出f(x)最佳平方逼最佳平方逼近近:投影投影现在学习的是第9页,共75页01.04.20239第三章 函数的最佳逼近(3.3)构造多元函数构造多元函数根据根据则则这时等式这时等式(3.4)意味着意味着(3.5)现在学习的是第10页,共75页01.04.202310第三章 函数的最佳逼近(3.5)(3.3)的极小值点的极小值点。(3
8、.4)也就是说,求出满足也就是说,求出满足等式等式(3.4)的的 pn*(x),等价于求出满足等价于求出满足等式等式(3.5)的的 c0*,c1*,cn*。由由(3.5)可知可知 c0*,c1*,cn*是是 n+1 元二次函数函数元二次函数函数现在学习的是第11页,共75页01.04.202311第三章 函数的最佳逼近而而n+1元函数元函数在区间在区间(-,+)上具有上具有一阶连续导函数一阶连续导函数,因此根据,因此根据极值原理极值原理,在,在最小值点最小值点 c0*,c1*,cn*处:处:而而于是于是即即现在学习的是第12页,共75页01.04.202312第三章 函数的最佳逼近利用内积利用
9、内积可以得到可以得到这是一个含有这是一个含有n+1个变量的方程组,具体形式为:个变量的方程组,具体形式为:现在学习的是第13页,共75页01.04.202313第三章 函数的最佳逼近再写成再写成矩阵形式为矩阵形式为现在学习的是第14页,共75页01.04.202314第三章 函数的最佳逼近这是关于这是关于n+1个个变量变量c0,c1,cn 的线性方程组,并称其为的线性方程组,并称其为法方程法方程组组,或者,或者正规方程组正规方程组。解此方程组,就可以得到解此方程组,就可以得到c0*,c1*,cn*,也就得到了也就得到了f(x)的的最佳平方逼近:最佳平方逼近:格拉姆格拉姆(Gram)矩阵矩阵最佳
10、平方逼近函数存在惟一最佳平方逼近函数存在惟一现在学习的是第15页,共75页01.04.202315第三章 函数的最佳逼近二、误差估计二、误差估计最佳平方逼近的平方误差为最佳平方逼近的平方误差为由方程组由方程组可得可得对于最佳逼近解对于最佳逼近解现在学习的是第16页,共75页01.04.202316第三章 函数的最佳逼近于是,最佳平方逼近于是,最佳平方逼近的平方误差为的平方误差为如果如果(3.6)则称则称(3.6)(3.6)为为 f(x)的在的在a,b上的最佳平方逼近上的最佳平方逼近n次多项式次多项式。n较大时,法方程组出现病态(第六章讲,实习题六6-3Hilbert矩阵),可取基函数为正交基函
11、数(如三角函数)现在学习的是第17页,共75页01.04.202317第三章 函数的最佳逼近*求求连续函数最佳平方逼近的步骤连续函数最佳平方逼近的步骤*1.给定给定a,b上的连续函数上的连续函数f(x),及子空间及子空间2.利用内积利用内积给出法方程组给出法方程组现在学习的是第18页,共75页01.04.202318第三章 函数的最佳逼近3.求出法方程组的解求出法方程组的解 c0*,c1*,cn*,得到最佳平方逼近得到最佳平方逼近4.求出求出平方平方误差误差称为称为均方均方误差误差现在学习的是第19页,共75页01.04.202319第三章 函数的最佳逼近 例例3.13.1求求 在在 上的上的
12、最佳平方逼近最佳平方逼近一次多项式,一次多项式,并估计误差。并估计误差。直接套用公式:直接套用公式:解:设解:设 令基函数为令基函数为 则需要求解的方程组为:则需要求解的方程组为:现在学习的是第20页,共75页01.04.202320第三章 函数的最佳逼近 这时由这时由 得到得到于是得到法方程组于是得到法方程组 现在学习的是第21页,共75页01.04.202321第三章 函数的最佳逼近解之得解之得 最佳平方逼近最佳平方逼近一次多项式为一次多项式为 关于误差,由误差估计式关于误差,由误差估计式现在学习的是第22页,共75页01.04.202322第三章 函数的最佳逼近得到得到 现在学习的是第2
13、3页,共75页01.04.202323第三章 函数的最佳逼近 例例3.2 求求 f(x)=arctanx 在在0,1 上的最佳平方逼近二次多项式,上的最佳平方逼近二次多项式,并估计误差。并估计误差。解:设解:设 P2(x)=c0+c1 x+c2x2,则则 需要写出法方程组需要写出法方程组 这时这时现在学习的是第24页,共75页01.04.202324第三章 函数的最佳逼近现在学习的是第25页,共75页01.04.202325第三章 函数的最佳逼近法方程组为法方程组为解得:解得:且且现在学习的是第26页,共75页01.04.202326第三章 函数的最佳逼近本节本节(2)小结小结1.1.何为连续
14、函数最佳平方逼近多项式?何为连续函数最佳平方逼近多项式?2.2.如何计算连续函数的最佳平方逼近如何计算连续函数的最佳平方逼近n次多项式?次多项式?3.如何估计最佳平方逼近如何估计最佳平方逼近n次多项式的误差?次多项式的误差?4.练习:试求函数练习:试求函数 f(x)=1/x 在区间在区间1,3上的最佳平方逼近一次多上的最佳平方逼近一次多项式并估计误差。项式并估计误差。现在学习的是第27页,共75页01.04.202327第三章 函数的最佳逼近3 3 离散数据拟合的最小二乘法离散数据拟合的最小二乘法 当我们得到的实验数据是准确值时,可以用代数插值的方当我们得到的实验数据是准确值时,可以用代数插值
15、的方法,求出原函数的近似表达式。法,求出原函数的近似表达式。经常由观察或测试可得到经常由观察或测试可得到 y=f(x)的一组离散数据的一组离散数据:但是,这组离散数据由观察或测试得到,往往并非完全精确,但是,这组离散数据由观察或测试得到,往往并非完全精确,如果用插值的方法来逼近,效果就不会太好如果用插值的方法来逼近,效果就不会太好。这时可以考虑用这时可以考虑用最小二乘法最小二乘法进行数据拟合,给出逼近曲线。其进行数据拟合,给出逼近曲线。其特点特点是:是:所求的逼近曲线不一定经过这些离散点,但却尽可能的靠近所求的逼近曲线不一定经过这些离散点,但却尽可能的靠近原曲线原曲线。(xi,yi),yi=f
16、(xi),i=0,1,m离散点的最佳平方逼近离散点的最佳平方逼近-几何上称为几何上称为曲线拟合曲线拟合(curve fitting)现在学习的是第28页,共75页01.04.202328第三章 函数的最佳逼近最小二乘拟合曲线最小二乘拟合曲线现在学习的是第29页,共75页01.04.202329第三章 函数的最佳逼近三次样条函数插值曲线三次样条函数插值曲线现在学习的是第30页,共75页01.04.202330第三章 函数的最佳逼近LagrangeLagrange插值曲线插值曲线现在学习的是第31页,共75页01.04.202331第三章 函数的最佳逼近一、数据拟合的最小二乘法的思想一、数据拟合的
17、最小二乘法的思想 已知离散数据:已知离散数据:(xi,yi),i=0,1,2,m,假设我们用函数假设我们用函数 逼近函数逼近函数f(x),则两个函数在每一个点,则两个函数在每一个点xi都会产生一个误差:都会产生一个误差:我们希望所求的逼近函数在每一个我们希望所求的逼近函数在每一个xi 处所产生的误差处所产生的误差i 的绝对值的绝对值|i|达最小达最小。但这样分别考虑太困难,所以我们应考虑整体误差。但这样分别考虑太困难,所以我们应考虑整体误差现在学习的是第32页,共75页01.04.202332第三章 函数的最佳逼近应该使应该使整体达最小整体达最小(误差的平方和最小)(误差的平方和最小)。通过这
18、种通过这种度量标准度量标准求得拟合曲线的方法,就称作求得拟合曲线的方法,就称作曲线拟合的最曲线拟合的最小二乘法小二乘法(最小二乘逼近最小二乘逼近)。按照以上思想求按照以上思想求 f(x)的拟合曲线(的拟合曲线(逼近函数逼近函数)时,首先需要)时,首先需要确定出确定出 f(x)所属的函数类,然后进一步求出具体函数,具体按照所属的函数类,然后进一步求出具体函数,具体按照以下步骤进行以下步骤进行。现在学习的是第33页,共75页01.04.202333第三章 函数的最佳逼近二、最小二乘法拟合曲线的步骤二、最小二乘法拟合曲线的步骤第二步:根据图示判断点第二步:根据图示判断点(xi,yi)所反映的函数类,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 函数 最佳 逼近 课件
限制150内