资源描述:
《数值计算方法-第6章-函数逼近ppt课件.ppt》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、第6章函数逼近实际问题中,通过测量或数值计算得到一批离散的数据,希望通过某种函数(曲线)来描述它,且使得它在某种意义下最“贴近”这批数据,这就是数据拟合,也称为函数逼近.一组实验数据:函数逼近的概念函数逼近的例子i1234Xi2468yi1.12.84.97.2从图形上可看出,数据分布接近于直线:如何选取a,b,使得直线“最好”地贴近于数据点?记,残差评判残差大小的标准?衡量残差大小的标准使残差的绝对值之和最小,即:第6章函数逼近这一标准虽然简单,但使用上不太方便.使残差绝对值最大的分量达到最小,即:这一方法称为最佳一致逼近.使
2、残差的平方和达到最小,即:这一方法称为最佳平方逼近,通常也称为曲线(数据)拟合的最小二乘法.该方法较简单,是应用中常用的一种方法.一.数据拟合的最小二乘法多项式拟合(xi,yi)第6章函数逼近最小二乘法的基本思想一组给定的数据点(i=0,1,…n)选取近似函数类H,寻求函数使得最小.即,要求F(a0,…,am)极小H通常采用比较简单的函数类,如低阶多项式,指数函数等.方法数据点:m次多项式:残差平方和:确定拟合多项式的系数第6章函数逼近j=0j=1j=m称为正则方程组,或法方程组.记为,可以证明,该方程组有唯一解(a0,a1,…
3、,am),从而得Pm(x).但要注意,系数矩阵A通常是病态的,其条件数cond(A)非常大.多项式拟合的例子第6章函数逼近用二次多项式拟合这批数据.i123456Xi012345yi0.3560.8051.0050.9420.6680.325解.二次拟合多项式:正则方程(n=6):拟合曲线图指数拟合方法取极小关于x的线性函数称为非线性最小二乘问题.即:第6章函数逼近当数据点接近于指数曲线分布时,可用指数函数进行拟合:改进的方法数据点拟合多项式:正则方程:线性最小二乘问题分式线性拟合当数据点分布接近于函数时第6章函数逼近作变换:对
4、数据点进行线性最小二乘拟合.参数a,b拟合曲线当数据点分布接近于函数时作变换:对数据点进行线性最小二乘拟合.参数a,b拟合曲线上面两种拟合中,参数a,b满足以下正则方程:最小二乘拟合的一般步骤第6章函数逼近通过观测数据点的形态,确定拟合函数的形式;对于一些简单的非线性问题,通过合适的变换化成线性问题;由最小二乘的正则方程确定拟合函数中的参数;对于非线性问题,将线性拟合函数再反变换成非线性函数.观测数据点取极小几点说明对于指数函数拟合问题:直接拟合(非线性)非线性拟合结果数据变换变换后的数据点线性拟合参数非线性拟合函数线性拟合反变
5、换实际问题中,各个数据点的重要性可能不相等,定义误差函数:加权系数线性最小二乘法的一般形式方法要求F极小,即第6章函数逼近数据点记:取拟合函数:线性无关函数组:误差函数:正则方程:由于线性无关,系数矩阵非奇异,方程有唯一解.下面证明就是所求的最小二乘函数.为权系数.定理6.1(最小二乘解)证明.对证毕#第6章函数逼近设(a0,a1,…,am)是满足上述正则方程的解,则是数据点(xi,yi)(i=1,…,n)的最小二乘函数.证明II即,正则方程的第k个方程正交多项式方法的提出第6章函数逼近正则方程一般是病态的,尤其是m比较大的情况
6、,这样得到的系数(a0,a1,......,am)误差较大.正则方程的系数矩阵为对角阵.若选取的函数组满足:最小二乘函数:若函数组满足上述条件,定义(正交函数族)则称为以为权,关于点集(x1,x2,…,xn)的正交函数族.正交多项式的构造第6章函数逼近其中,证明过程类似于后面的定理6.4.利用正交多项式求拟合多项式的过程第6章函数逼近得拟合多项式:构造正交函数系:解正则方程,得系数:a0,a1,…,am.i123456Xi012345yi0.3560.8051.0050.9420.6680.325例.利用正交函数求下列数据的最小
7、二乘二次(m=2)拟合多项式解.n=6,取构造正交多项式系:例(续)第6章函数逼近解正则方程,得拟合函数的各个系数.拟合多项式:与多项式直接拟合结果完全一致.二.正交多项式基本概念第6章函数逼近定义6.2(正交函数系)特例:正交函数系的例子若函数系满足:则称函数系在[a,b]上关于权函数的正交函数系.是区间上关于权函数的正交函数系.三角函数系:称为标准正交函数系.正交函数系的性质——由下面两个定理给出:定理6.2(线性无关性)证明[采用反证法]假设线性相关,即存在不全为0的系数ci,使第6章函数逼近不妨设由正交性矛盾.故线性无关
8、.正交函数系必定是线性无关的.定理6.3(正交多项式的充要条件)设是最高次系数非零的k次多项式,则是[a,b]上关于权函数的正交多项式的充要条件是,对任意次数不超过k-1次的多项式,都有:定理6.3(续)证毕#第6章函数逼近证明[充分性]若对成立.取即另外,故为