函数逼近.doc

上传人:本田雅阁 文档编号:2056093 上传时间:2019-02-09 格式:DOC 页数:28 大小:630.51KB
返回 下载 相关 举报
函数逼近.doc_第1页
第1页 / 共28页
函数逼近.doc_第2页
第2页 / 共28页
函数逼近.doc_第3页
第3页 / 共28页
亲,该文档总共28页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述

《函数逼近.doc》由会员分享,可在线阅读,更多相关《函数逼近.doc(28页珍藏版)》请在三一文库上搜索。

1、第七章 函数逼近用简单的函数p(x)近似地代替函数f (x),是计算数学中最基本的概念和方法之一。近似代替又称为逼近,函数f (x)称为被逼近的函数,p (x)称为逼近函数,两者之差 称为逼近的误差或余项在计算数学里,所谓简单的函数主要是指可以用加、减、乘、除四则运算进行计算的函数,如有理分式函数、多项式等。由于多项式最简单,计算其值只需用到加、减与乘三种运算,且求其微分和积分都很方便,所以常用它来作为逼近函数,而被逼近的函数f (x)一般是一个比较复杂的不易计算的函数或以表格形式给出的函数。第六章介绍的插值法实际上也是函数逼近的一种方法。不过,它要求函数p (x)与f (x)在节点处具有相同

2、的函数值 (甚至要求有相同的导数值),但在非节点处,p (x) 虽然有可能很好地逼f (x),但也可能使逼近f (x) 的误差很大,如果实际问题要求p (x)在区间a, b 上每一点都“很好”地逼近的话,用插值多项式p (x) 去逼近f (x)有时就要失败,所谓龙格现象,就是典型一例。大家知道,用f (x)的泰勒(Taylor)展开式的部分和去逼近函数f (x),也是常用的方法。这种方法的特点是:x越接近于x0,误差就越小,x越偏离x0,误差就越大。若要使这种逼近在整个所讨论的区间上都达到精度要求,则需取很多项,这样,计算工作量就大大增加。因此,如何在给定精度下,求出计算量最小的近似式,这就是

3、函数逼近要解决的问题,这个问题的一般提法是:对于函数类A中给定的函数f (x),要求在另一类较简单的且便于计算的函数类B( A)中寻找一个函数p (x),使p (x)与f (x)之差在某种度量意义下最小。一般,最常见的函数A是区间a, b上的连续函数,记作Ca, b。最常用的函数类B有代数多项式、三角多项式以及有理分式函数等。最常用的度量标准有两种:(一) 一致逼近以函数f (x)和p (x)的最大误差作为度量误差f (x) - p (x)的“大小”的标准,在这种意义下的函数逼近称为一致逼近或均匀逼近,讲得更具体一点,也即对于任意给定的一个小正数e 0,如果存在函数p (x),使不等式成立,则

4、称该函数p (x)在区间a, b上一致逼近或均匀逼近于函数f (x)。(二)平方逼近:如果我们采用作为度量误差的“大小”的标准,在这种意义下的函数逼近称为平方逼近或均方逼近。这种方法要比一致逼近的相应问题简单得多。本章主要介绍在这两种度量标准下用代数多项式p (x)去逼近区间a, b上的连续函数,也就是介绍函数的最佳一致逼近多项式和最佳平方逼近多项式。由于正交多项式是函数逼近的重要工具,因此,下面先介绍几种常见的正交多项式。1 正交多项式一、正交函数系的概念高等数学中介绍傅立叶(Fourier)级数时,证明过函数系;1, cosx,sinx,cos2x,sin2x,connx,sinnx, (

5、7.1)中任何两个函数的乘积在区间-p ,p 上的积分都等于0。我们称这个函数中任何两个函数在-p ,p 上是正交的,并且称这个函数为一个正交函数系。若对(7.1)中的每一个函数再分别乘以适当的数,使之成为:(7.2)那么这个函数系在-p ,p 上不仅保持正交的性质,而且还地标准化的(规范的),亦即每一个函数自乘之积,在-p ,p 上的积分是1。为了使讨论更具有一般性,先要介绍一些基本概念。1权函数的概念定义7.1 设r (x)定义在有限或无限区间a, b上,如果具有下列性质:(1) r (x) 0,对任意x a, b,(2) 积分存在,(n = 0, 1, 2, ),(3) 对非负的连续函数

6、g (x) 若。则在(a, b)上g (x) 0,我们就称r (x)为a, b上的权函数。在正交多项式的讨论中,会遇到各种有意义的权函数,常用的权函数有:;等等。2内积的概念定义7.2 设f (x),g (x) C a, b,r (x)是a, b上的权函数,则称为f (x)与g (x)在a, b上以r (x)为权函数的内积。内积有如下性质:(1) (f, f )0,且(f, f )=0 f = 0;(2) (f, g) = (g, f );(3) (f1 + f2, g ) = (f1, g) + (f2, g);(4) 对任意实数k,(kf, g) = k (f, g )。这些性质,由内积的

7、定义不难得到证明。3正交性的概念定义7.3 设f (x),g(x) C a, b若则称f (x)与g (x)在a, b上带权r (x)正交。定义7.4 设在a, b上给定函数系,若满足条件则称函数系jk (x)是a, b上带权r (x)的正交函数系,特别地,当Ak 1时,则称该函数系为标准正交函数系。若定义7.4中的函数系为多项式函数系,则称为以r (x)为权的在a, b上的正交多项式系。并称pn(x)是a, b上带权r (x)的n次正交多项式。例1 验证多项式: 在上带权r (x) = 1两两正交。解 容易验证而 由定义7.4,结论成立。有了以上的基本概念,下面我们介绍几个常用的正交多项式。

8、二、常用的正交多项式1切比雪夫()多项式切比雪夫多项式具有很多重要性质,是函数逼近的重要工具,并且有广泛的应用。定义7.5 称多项式(7.3)为n次的切比雪夫多项式(第一类)。切比雪夫多项式Tn (x)具有以下性质:(1) 正交性:由 Tn (x)所组成的序列 Tn (x)是在区间-1, 1上带权的正交多项式序列。且(7.4)证 因为,令 , 则,于是 (2) 递推关系相邻的三个切比雪夫多项式具有三项递推关系式:(7.5)证 显然,n = 0时,时,当n1时,令x = cosq ,则由三角恒等式即得移项就得上述递推关系(7.5)。由三项递推关系式可依次写出如下常用的前面几个切比雪夫多项式的表达

9、式:可见Tn (x)也是普通的n次多项式。(3) 奇偶性:切比雪夫多项式Tn (x),当n为奇数时为奇函数;n为偶数时为偶函数。这是因为(4) Tn (x)在区间-1, 1上有n个不同的零点。证 由于 令 有 所以在区间0q p 上有n个值使 即在0,p中有n个不同的零点,且由于Tn (x)是n次多项式,所以至多有n个零点,现已找到n个不同的零点,则每一个xk都是Tn (x)的单重零点。显然,Tn (x)在-1, 1中的零点都是实的、互异的,且全部在-1, 1内。(5) Tn (x) 在-1, 1上有n + 1个不同的极值点使Tn (x)轮流取得最大值1和最小值-1。证 因为在区间0q p 上

10、有n + 1个点使cosnq 顺次取+1及-1,因此,n + 1个点 使Tn (x)顺次为+1及-1,即由于cosnq 之最大值是+1,最小值是-1,因此,我们把这n + 1个点叫做Tn (x)在-1, +1上的极值点,也称它为Tn (x)的交错点组。这是切比雪夫多项式的一个重要性质。如果将Tn (x)的零点xk和极值点按大小排列,则有(6)切比雪夫多项式的极值性质显然,由三项递推关系式容易验证Tn (x)的最高次项系数为2n-1 (n = 1, 2, )。譬如 T1(x) = x, 最高次项系数为21-1, T2(x) = 2x2 1,最高次项系数为22-1。由三项递推关系可看到,每升高一次

11、幂要乘2,则是一个首项系数为1的n次多项式,又由性质(5)知在点顺次达到它在-1, 1上的极值即若用Hn (x)表示首项系数为1的n次多项式的全体,则有定理7.1 在-1x 1上,在首项系数为1的一切n次多项式Hn (x)中与零的偏差最小,且其偏差为即,对于任何,有证明 (用反证法)假若存在首项系数为1的另一n次多项式它对零的偏差比与零的偏差还小,即令 由于,p (x)均属于Hn (x),则q (x)是一个次数不超过n 1的多项式。在Tn (x)的交错点组处,由于即 在处轮流取到因而则有即,q (x)在n + 1个交错点处轮流取正负值,由连续函数的性质知,q (x)应当具有n个零点,但已知q

12、(x)至多是n 1次多项式,因而q (x) 0,所以,。这个结果与假设矛盾,从而定理得证。从这个定理知,所有首项系数为1的n次多项式在区间-1, 1上的最大值满足定理7.1 称为切比雪夫多项式的极性,这种极性也是切比雪夫多项式的一种重要性质。作为应用,介绍一下多项式插值余项的极小化。设在-1, 1上给定n + 1个互异的节点x0, x1, , xn,函数f (x)在-1, 1上具有n + 1阶连续导数f (n+1) (x),对f (x)作多项式插值时,由拉格朗日插值的余项表示式其中,若得到,则有显然,余项的大小,取决于因子的大小。现在我们提出一个问题:怎样选取节点xj ( j = 0, 1,

13、2, , n)能使 尽可能小?由于是一个最高次项系数为1的n + 1次多项式,由Tn(x)的极性讨论知,当xj满足时,取得极小,亦即只要插值节点xk取成n + 1次比雪夫多项式的零点则插值公式的余项在全区间-1, 1上的最大绝对值为极小,此时,有余项公式:注意:如果插值区间是a, b,而不是-1, 1,总可以作变换把函数变换成其中 -1t 1,即可将定义在区间a,b上的函数f (x)化为新变量t的定义在区间-1, 1上的函数g (t)。以上讨论在数值积分中将要用到。2勒让德(Legendre)多项式定义7.6 多项式(7.6)称为n次勒让德多项式。显然 p0 (x) = 1pn (x)的最高次

14、幂xn项的系数是勒让德多项式具有以下一些重要性质:(1) 正交性勒让德多项式序列pn(x)是在-1, 1上带权r (x) = 1的正交多项式序列。即(7.7)(2) 递推关系相邻的三个勒让德多项式具有三项递推关系式:(7.8)由此可得,(3) 奇偶性:当n为偶数时,pn (x)为偶函数;当n为奇数时,pn (x)为奇函数。即这是由于pn (x)的表达式中的(x2 - 1)n是偶次多项式,经过偶次求导仍为偶次多项式,经过奇次求导则为奇次多项式,故n为偶数时,pn (x)为偶函数;n为奇数时,pn (x)为奇函数。(4) pn (x)的n个零点都是实的、相异的,且全部在区间-1, 1内部。3其它常

15、用的正交多项式一般而言,如果给定的区间a, b不同且权函数不同,正交多项式也就不同,我们再介绍几个常用的正交多项式。(1) 第二类切比雪夫多项式定义7.7 称(7.9)为第二类切比雪夫多项式。可以证明un(x)是在区间-1, 1上带权函数的正交多项式序列。即(7.10) 相邻的三项具有递推关系式:(7.11)(2)拉盖尔(Laguerre)多项式定义7.8 称多项式(7.12)为拉盖尔多项式。可证 Ln(x)是在区间0, +上带权r (x) = e-x的正交多项式序列。即(7.13) 相邻的三项具有递推关系式:(7.14)(3)埃尔米特(Hermite)多项式定义7.9 称多项式(7.15)为

16、埃尔米特多项式。可以证明 Hn(x)是在区间(-, +)上带权函数的正交多项式序列。即(7.16) 相邻的三项具有递推关系式:(7.17)2 最佳一致逼近一、最佳一致逼近的概念定义7.10 设函数f (x)是区间a, b上的连续函数,对于任意给定的e 0,如果存在多项式p (x),使不等式成立,则称多项式p (x)在区间a, b上一致逼近(或均匀逼近)于函数f (x)。那么,对于在区间a, b上的连续函数f (x),是否存在多项式p(x)一致逼近于f (x)呢?这个问题有许多人研究过。德国数学家维尔斯特拉斯(Weierstrass)在1885年曾给出下述著名定理。维尔斯特拉斯定理 若f (x)

17、是区间a, b上的连续函数,则对于任意e 0,总存在多项式p (x),使对一切a x b有证明从略。维尔斯特拉斯定理表明,连续函数f (x)可以用多项式p (x)逼近到任意精确程度,但维尔斯特拉斯定理只在理论上肯定了闭区间上的连续函数可以用多项式以任意精确度来逼近,并没有给出确定逼近得最快的多项式的方法。事实上,如果精确度要求较高,则用来逼近的多项式的次数一般也很高,这就增加了计算工作量。因而,在实际计算时,我们总量希望在一定的精确度要求下,逼近多项式的次数越低越好。切比雪夫从这样的观点去研究一致逼近问题,他不让逼近多项式的次数n趋于无穷大,而是先把n加以固定。对于给定的a, b上的连续函数f

18、 (x),他提出在次数不超过n的多项式的集合pn中去寻找一个多项式,使它在a, b上“最佳地逼近”f (x)。这里最佳逼近的意思是指对f (x)的偏差。和其它任一p(x) pn 对f (x)的偏差比较时是最小的,也就是说(7.18)这就是通常所谓的最佳一致逼近问题,也称为切比雪夫逼近问题。若这样的存在,则称是函数f (x)在区间a, b上的n次最佳一致逼近多项式,简称为最佳逼近多项式。现在要问:最佳逼近多项式是否存在?是否唯一?如何构造?我们不妨设n次多项式显然应与p (x)的系数a0, a1, ,an有关。若记 则j 应是n + 1个系数a0, a1, ,an的正值连续函数。我们称多元函数j

19、 (a0, a1, ,an)的最小值 (7.19)为f (x)与p (x)在a, b上最小偏差。对照公式(7.18)可知,寻求f (x)在a, b上的n次最佳一致逼近多项式的问题就归结为求多元函数j (a0, a1, ,an)的最小值问题。可以证明,存在唯一的能使j 成立,也即存在唯一的满足关系式这就说明,对于a, b上的任意连续函数f (x),其n次最佳逼近多项式是存在且唯一的,(因证明较繁,这里我们不予讨论)。下面我们主要介绍的构造。二、最佳一致逼近多项式的求法定理7.2 是a, b上连续函数f (x)的n次最佳一致逼近多项式的充分必要条件是在区间a, b上至少有n + 2个点使得其中m

20、是f (x)与在a, b上的偏差,即(7.20)s 为“1”或“-1”。因证明较复杂,这里从略。点集x1, x2, ,xn+2称为切比雪夫交错点组。其中每一个xk (k = 1, 2, ,n+2)称为交错点。定理7.2表明,若用最佳一致逼近多项式 来近似代替f(x),则误差在a, b上的分布是十分均匀的。定理7.2 常称为切比雪夫定理,由此定理还可推得下面两个有用的推论。推论1 设则f (x)在pn中的最佳一致逼近多项式,就是f (x)在a, b上的某个n次拉格朗日插值多项式。证明 设是f (x)在a, b上的n次最佳一致逼近多项式,根据切比雪夫定理可知,连续函数在区间a, b上至少有n +

21、2个点使其交替变号,这就是说,方程=0在a, b上至少有n + 1个根,即存在n + 1个点 (k = 0, 1, 2, , n ),使,(k = 0, 1, 2, , n )所以,实际上就是以为插值节点的f (x)的n次拉格朗日插多项式。推论1表明:如果能适当选择插值节点,求出使偏差为最小的n次拉格朗日插值多项式,那么就能求得f (x)的n次最佳一致逼近多项式。推论2 如果函数f (x)在区间a, b上有n + 1阶导数,且f (n+1)(x)在a, b上恒为正(或负),那么区间a, b的端点a和b都属于的交错点组。证明(反证法)设a (或b)不属于的交错点组,那么函数R (x) = 在开区

22、间(a, b)内至少有n +1个点,使其取得最大值和最小值,则由取得极值的必要条件,必有,(k = 1, 2, , n+1 )反复用罗尔定理可知,在(a, b)由至少存在一点x ,使但是 ,从而有这与在a, b上恒为正(或负)的已知条件矛盾。切比雪夫定理不仅给出了最佳一致逼近多项式的特性,并从理论上给出了寻找最佳一致逼近多项式的方法。但通常,在讨论连续函数的多项式逼近中,要为某一连续函数求取最佳一致逼近多项式往往是十分困难的。下面我们仅介绍比较简单的线性最佳一致逼近多项式的求法以及切比雪夫多项在求最佳一致逼近时的应用。1线性最佳一致逼近多项式的求法设函数f (x)在a, b上有二阶导数,且在a

23、, b上不变号(即恒为正或负),则可按下面方法求f (x)在a, b上的线性最佳一致逼近多项式。设,则根据切比雪夫定理,在a, b上至少存在三个点:,使(7.21)其中s = 1,k = 1, 2, 3。由于在a, b上不变号,根据推论2可知,区间a, b的两个端点a, b都属于的交错点组,即有则另一交错点x2必位于a, b的内部,且它是的极值点,所以,而从而(即等于在点x2处的斜率)。又因为在a, b上不变号,所以一阶导数在a, b上严格单调,因此它只能取一次,这就证明了在a, b内除了点x2外不能再有其它极值点,也即在(a, b)内有,且只有一个极值点,它就是交错点x2。于是由(7.21)

24、式得:解此方程组得把代入的表示式,即求得f (x)在a, b上线性最佳一致逼近多项式为(7.24)线性最佳一致逼近的几何意义是:直线与弦AB平行且过线段AC的中点D,其方程为见图7-1,其中图7-1例7.2 求函数在上的线性最佳一致逼近多项式。解 显然在上不变号,故,x3 = 1,由(7.22)式得由 及(7.22)式得,所以 再由(7.23)式得故函数这就是函数在区间1/4, 1上的线性最佳一致逼近多项式。2切比雪夫多项式在函数逼近中的应用在1中我们已经讨论了切比雪夫多项式的极值性质,也即在区间-1, 1上,所有最高次项系数为1的一切n次多项式中,n次切比雪夫多项项Tn (x)与1/2n-1

25、的乘积与零的偏差最小,且其偏差为1/2(n-1)。利用这个极值性质,切比雪夫多面式就成为-1, 1上逼近其它函数f (x)的重要工具。下面分两种情况讨论(1) 若f (x)是n次多项式,则它的n 1次最佳一致逼近多项式能精确求出。例7.3 已知,求其在-1, 1上的二次最佳一致逼近多项式。解 令它是首项系数为1的3次多项式,要它在-1, 1上的最大值最小,由切比雪夫多项式的极值性质知,在-1, 1上,p3中与零偏差最小的首项系数为1的三次多项式是,所以令其中则就是f (x)在-1, 1上的二次最佳一致逼近多项式。(2) 若f (x)不是多项式,也可利用切比雪夫多项式的极值性质,求出f (x)的

26、近似最佳一致逼近多项式。常用的有“切比雪夫插值法”和“缩减幂级数法”两种近似方法。我们将简单地介绍用切比雪夫多项式来降低逼近多项式的次数。我们常要在一定区间上求一个函数在一定误差范围内的逼近多项式,自然,我们希望此逼近多项式的次数越低越好。由于xk也可用切比雪夫多项式Tn(x)表示,即若把普通n次多项式中的所有用上述切比雪夫多项式去代替,则pn(x)可改写为在满足误差要求的情况下,可利用切比雪夫多项式的极值性质,把pn(x)的高次幂的项缩减下来,使它成为m (mn - 1)次多项式,这个m次多项式可以作为在已给精度下的f (x)在-1, 1上的近似的最佳一致逼近多项式。下面通过简单的例子来说明

27、这种方法。例7.4 设在区间-1, 1上,要计算f (x) = ex,欲找一个近似多项式p (x)近似代替ex,使误差e 0.01。解 先对f (x) = ex在x = 0处作泰勒展开,有若取前六项之和作为ex的近似,这时截断误差R5(x)满足即误差确实满足所提要求,但注意p5(x)是一个5次多项式。若利用xk和Tk (x)的关系式,p5(x)也可表示为这里,我们看到Tk的足标k越小,Tk的系数就越小,由于在上,(k = 0, 1, 2, ),故可以略去次数Tk (x)项,而这就意味着逼近多项式的次数降低,这正是我们希望达到的。此题中,若略去其最后的两项,则所增添的误差为从而在上,若用近似代替

28、ex的总误差不超过0.0038 + 0.0058 = 0.0096 0.01亦即用y(x)代替ex 也满足所提要求。再利用Tk (x)和xk之间的关系,把y (x)重写为通常的x的乘幂形式,则这个y (x)是一个3次多项式,比原来的p5 (x)降低2次,从而用y (x)求值的计算工作量比用 p5(x)少。这个方法为缩减幂级数法。上述方法既使逼近多项式的次数降低了,而且还可以使误差的分布更为均匀,因此可以作为f (x)的近似最佳一致逼近多项式。对于任一有限区间a, b上的逼近问题,可以通过变量替换把a, b区间转化为-1, 1上的逼近问题类似讨论。3 最佳平方逼近一般而言,在a, b上对给定的函

29、数求它的一致逼近函数比较困难,下面我们介绍在a, b上较易计算的另一种逼近方法最佳平方逼近。一、预备知识1函数系的线性关系定义7.11若函数,在区间a, b上连续,如果关系式当且仅当时才成立,则称函数在a, b上是线性无关的,否则称线性相关。如果函数系jk (x)(k = 0, 1, 2, )中的任何有限个函数线性无关,则称函数系jk (x)为线性无关函数系,例如1, x, , xn, 就是在区间a, b上的线性无关函数系。设是a, b上线性无关的连续函数a0, a1, , an是任意实数,则的全体是Ca, b的一个子集,记为并称是生成集合的一个基底。例如 Pn = Span 1, x, x2

30、, , xn表示由基底1, x, , xn, 生成的多项式集合。下面给出判断函数系jk (x)(k = 0, 1, 2, n)线性无关的一个充要条件。定理7.3 连续函数在a, b上线性无关的充分必要条件是它们的克莱姆(Gram)行列式Gn 0,其中(7.25)证明 设a0, a1, , an是一组实数,使现分别用乘上式,然后在a, b上积分,于是得方程组显然,上述齐次线性代数方程组只有零解的充要条件是它的系数行列式。证毕。根据本定理即可断定:在区间a, b上带权r (x)的正交函数系是线性无关的函数系。2广义多项式我们已知,n次代数多项式是函数1, x, x2, , xn的线性组合。现在我们

31、推广多项式的概念。设函数系,线性无关,则其有限项的线性组合(7.26)称为广义多项式。例如三角多项式就是一个广义多项式。二、函数的最佳平方逼近我们已知,用平方误差取得最小,作为度量标准研究函数的逼近函数p (x),就是最佳平方逼近问题。若设表示任意一个不高于n次的多项式函数,Pn = Span 1, x, , xn,则最佳平方逼近多项式的定义可叙述为:定义7.12 对于给定的函数,若n次多项式满足关系式(7.27)则称S*(x)为f (x)在区间a, b上的n次最佳平方逼近多项式。实用上,为用使问题的讨论更有一般性,我们可以把定义进行更一般的推广,若把xj推广为一般的连续函数jj(x),且要求

32、j0, j1, , jn线性无关,并且积分可带权函数,则定义7.12可改为定义7.13 对于给定的函数如果存在使(7.28)则称S*(x)是f (x)在集合F中的最佳平方逼近函数。显然定义7.12是定义(7.13)中取r (x) = 1, ji (x) = xj 的特殊情况,作这样的推广,使我们可以考虑更一般的,象之类在实际问题中有用的带权函数的平方逼近问题。在时,满足条件(7.28)的S* (x),就是函数f (x)的n次最佳平方逼近多项式。显然,求最佳平方逼近函数的问题可归结为求它的系数,使多元函数取得极小值,也即点()是I (a0, ,an)的极点。由于I (a0, a1, ,an)是关

33、于a0, a1, ,an的二次函数,利用多元函数取得极值的必要条件, (k = 0, 1, 2, , n)即(7.29)得方程组 (7.30)如采用函数内积记号那么,方程组可以简写为 这是一个包含n + 1个未知元a0, a1, , an的n + 1阶线性代数方程组,写成矩阵形式为(7.31)此方程组叫做求aj (j = 0, 1, 2, , n)的法方程组。显然,其系数行列式就是克莱姆行列式Gn = Gn (j0, j1, , jn)。由于j0, j1, , jn线性无关,故Gn 0,于是上述方程组存在唯一解。从而肯定了函数f (x)在F中如果存在最佳平方逼近函数,则必是例7.5 。求(不超

34、过二次的多项式全体),使最小。解:设,即取方程组(7.31)成为其中,于是得解得所以偏差由于在实际问题中方程组(7.31)的系数矩阵常是病态的或高度病态的,为了得到有效的结果,我们将选为带权的正交系,则由(7.30)得所以(7.32)注意到由(7.12)算出的ak就是f (x)的广义福里哀级数的系数,而得到的逼近函数就是广义福里哀级数的部分和。这样,不仅计算ak容易,而且不受n的影响。因为要确定方程组必须先确定n。n变了,方程组也随之改变。现在,只需把a0, a1, , an算好,那么,对任意的n, 就是我们要找的多项式(级数的部分和)。我们仍以f (x) = x4为例。利用福里哀一切比雪夫级数的部分和来求x4在H2中的最佳平方逼近多项式。已知于是所以于是所求多项式为即

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 其他


经营许可证编号:宁ICP备18001539号-1