最佳平方逼近.ppt

《最佳平方逼近.ppt》由会员上传分享，免费在线阅读，更多相关内容在PPT专区-天天文库。

1、1.函数逼近在数值计算中经常遇到求函数值的问题，手算时常常通过函数表求得，用计算机计算时若把函数表存入内存进行查表，则占用单元太多，不如直接用公式计算方便。因此，我们希望求出便于计算且计算量省的公式近似已知函数f(x)，例如，泰勒展开式的部分和就是f(x)的一种近似公式，用它求x0附近的函数值f(x)，误差较小，当

2、x-x0

3、较大时，误差就很大。例如f(x)=ex在[-1,1]上用：近似ex，其误差：于是误差分布如图：xy-11它在整个区间上误差较大，若在计算机上用这种方法计算ex，如精度要求较高，则需取很多项，这样即费时又多占存储单

4、元。因此，我们要求在给定精度下计算次数最少的近似公式，这就是函数逼近要解决的问题。定义近似代替又称为逼近，函数f(x)称为被逼近函数；P(x)称为逼近函数，两者之差称为逼近的误差。函数逼近问题可叙述为：对函数类A中给定的函数f(x)，需要在另一类较简单的便于计算的函数类B（B∈A）中，找一个函数P(x)，使P(x)与f(x)之差在某种度量意义下达到最小。函数类A通常是区间[a,b]上的连续函数，记作C[a,b]；函数类B通常是代数多项式，有理多项式，三角多项式，分段多项式等容易计算的函数。最常用的度量标准有两种：1、一致逼近（

5、均匀逼近）以作为度量误差f(x)-P(x)的“大小”标准。2、平方逼近（均方逼近）以作为度量误差f(x)-P(x)的“大小”标准。5.6函数的最佳平方逼近5.6.1最佳平方逼近的概念与解法一、最佳逼近的意义设{0x,1x,,nx}C[a,b],它们线性无关.又给定f(x)C[a,b],求p*(x)HnSpan{0x,1x,,nx},使得f(x)p*(x)在某种意义下最小.二、最佳平方逼近的概念定义对于给定的f(x)C[a,b],若有p*(x)Hn,使得(fp*,fp*)min

6、{(fp,fp)

7、pHn}.则称p*(x)是(在子空间Hn中)对f(x)的最佳平方逼近函数.下面用到的内积为定理5.7设f(x)C[a,b]，p*(x)Hn,在Hn中,p*(x)是对f(x)最佳平方逼近的函数(fp*,j)=0,j=0,1,…,n.其中,{0x,1x,,nx}为子空间Hn的一组基.证:()反证法,设有函数kx,使得(fp*,k)k0,令q(x)p*(x)kxk/(k,k),显然,q(x)Hn.利用内积的性质,可得k这说明,p*(x)不是对f(x)

8、最佳平方逼近的函数,矛盾.()若(fp*,j)=0,j=0,1,…,n成立,对任意的p(x)Hn,有p*在Hn中是对f的最佳平方逼近函数.证毕.即设f(x)C[a,b],p*(x)Hn,在Hn中,p*(x)是对f(x)最佳平方逼近的函数对任意的p(x)Hn,均有(fp*,p)=0.由于及定理5.8设f(x)C[a,b],在子空间Hn中,对f(x)最佳平方逼近的函数是唯一的.证明假定,在Hn中,p(x)和q(x)都是对f(x)最佳平方逼近的函数,由定理5.7的系,知这说明,p(x)q(x)于[a,b].三、最佳平方

9、逼近函数的求解利用(fp*,j)=0,可求出最佳平方逼近函数p*.设故这是一个以c*0,c*1,,c*n为未知数的线性方程组.称(5.82)为法方程或正规方程.法方程的矩阵形式是由于0x,1x,,nx线性无关,可以推得上系数阵是非奇异的.故(5.82)有唯一解{c*j}.四、最佳平方逼近的误差记(fp*,fp*),称其为最佳平方逼近误差,利用(fp*,p*)=0,可有例6定义内积试在H1=Span{1,x}中寻求对于f(x)=的最佳平方逼近元素p(x).解法方程为当0x,1x,,n

10、x,是正交系时,求解最佳平方逼近式(5.82)中的系数非常容易.目标:求下面的最佳平方逼近式中的系数变为5.6.2正交系在最佳平方逼近中的应用一、Legendre多项式的应用给定函数f(x)C[a,b],求f(x)的Legendre最佳平方逼近.Legendre多项式的权函为(x)1,故内积L-正交多项式为L0x,L1x,,Lnx,用(5.83),有函数f的L-最佳平方逼近函数为（5.85）遇到区间[a,b],通过下面的变换把问题转化到[-1,1]上处理.函数f(x)的Legendre无穷级数例7求函数

11、f(x)ex在[-1,1]上的Legendre三次最佳平方逼近多项式.解:前4个Legendre多项式为使用求系数2.3504,所求三次最佳平方逼近多项式为例8求f(x)=在区间[0,1]上的一次最佳平方逼近多项式.解