第三章函数逼近
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第七章 函数逼近
用简单的函数p (x )近似地代替函数f (x ),是计算数学中最基本的概念和方法之一。近
似代替又称为逼近,函数f (x )称为被逼近的函数,p (x )称为逼近函数,两者之差
)()()(x p x f x R -=
称为逼近的误差或余项
在计算数学里,所谓简单的函数主要是指可以用加、减、乘、除四则运算进行计算的函
数,如有理分式函数、多项式等。由于多项式最简单,计算其值只需用到加、减与乘三种运算,且求其微分和积分都很方便,所以常用它来作为逼近函数,而被逼近的函数f (x )一般是一个比较复杂的不易计算的函数或以表格形式给出的函数。
第六章介绍的插值法实际上也是函数逼近的一种方法。不过,它要求函数p (x )与f (x )
在节点处具有相同的函数值 (甚至要求有相同的导数值),但在非节点处,p (x ) 虽然有可能很好地逼f (x ),但也可能使逼近f (x ) 的误差很大,如果实际问题要求p (x )在区间[a , b ] 上每一点都“很好”地逼近的话,用插值多项式p (x ) 去逼近f (x )有时就要失败,所谓龙格现象,就是典型一例。
大家知道,用f (x )的泰勒(Taylor)展开式
)()()!
1()()(!)()(!2)()
)(()()(010)1(00)(200000之间与在x x x x n f x x n x f x x x f x x x f x f x f n n n n ξξ++-++-++-''+-'+=
的部分和去逼近函数f (x ),也是常用的方法。这种方法的特点是:x 越接近于x 0,误差就越小,x 越偏离x 0,误差就越大。若要使这种逼近在整个所讨论的区间上都达到精度要求,则需取很多项,这样,计算工作量就大大增加。因此,如何在给定精度下,求出计算量最小的
近似式,这就是函数逼近要解决的问题,这个问题的一般提法是:
对于函数类A 中给定的函数f (x ),要求在另一类较简单的且便于计算的函数类B (⊂ A )
中寻找一个函数p (x ),使p (x )与f (x )之差在某种度量意义下最小。 一般,最常见的函数A 是区间[a , b ]上的连续函数,记作C [a , b ]。 最常用的函数类B 有代数多项式、三角多项式以及有理分式函数等。 最常用的度量标准有两种: (一) 一致逼近
以函数f (x )和p (x )的最大误差
)()(max ]
,[x p x f b a x -∈
作为度量误差f (x ) - p (x )的“大小”的标准,在这种意义下的函数逼近称为一致逼近或均匀逼近,讲得更具体一点,也即对于任意给定的一个小正数ε >0,如果存在函数p (x ),使不等式
ε<-<<)()(max x p x f b
x a
成立,则称该函数p (x )在区间[a, b ]上一致逼近或均匀逼近于函数f (x )。 (二)平方逼近: 如果我们采用
dx x p x f b
a
⎰
-2)]()([
作为度量误差)()(x p x f -的“大小”的标准,在这种意义下的函数逼近称为平方逼近或均方逼近。这种方法要比一致逼近的相应问题简单得多。
本章主要介绍在这两种度量标准下用代数多项式p (x )去逼近区间[a, b ]上的连续函数,
也就是介绍函数的最佳一致逼近多项式和最佳平方逼近多项式。
由于正交多项式是函数逼近的重要工具,因此,下面先介绍几种常见的正交多项式。
§1 正交多项式 一、正交函数系的概念
高等数学中介绍傅立叶(Fourier)级数时,证明过函数系;
1, cos x ,sin x ,cos2x ,sin2x ,…,con nx ,sin nx ,… (7.1) 中任何两个函数的乘积在区间[-π ,π ]上的积分都等于0。我们称这个函数中任何两个函数在[-π ,π ]上是正交的,并且称这个函数为一个正交函数系。若对(7.1)中的每一个函数再分别乘以适当的数,使之成为:
nx nx x x sin 1
,
cos 1
,
,,sin 1
,
cos 1
,
21π
π
π
π
π
(7.2)
那么这个函数系在[-π ,π ]上不仅保持正交的性质,而且还地标准化的(规范的),亦即每
一个函数自乘之积,在[-π ,π ]上的积分是1。
为了使讨论更具有一般性,先要介绍一些基本概念。 1.权函数的概念 定义7.1 设ρ (x )定义在有限或无限区间[a , b ]上,如果具有下列性质: (1) ρ (x ) ≥0,对任意x ∈[a , b ], (2) 积分
dx x x n
b
a
)(ρ⎰
存在,(n = 0, 1, 2, …),
(3) 对非负的连续函数g (x ) 若
⎰
=b
a
dx x x g 0)()(ρ。
则在(a , b )上g (x ) ≡ 0,我们就称ρ (x )为[a , b ]上的权函数。 在正交多项式的讨论中,会遇到各种有意义的权函数,常用的权函数有: 1)(],1,1[],[=-=x b a ρ;
2
11)(],1,1[],[x
x b a -=
-=ρ
x e x b a -=∞=)(],,0[],[ρ
2
)(],,[],[x e x b a -=∞+-∞=ρ
等等。 2.内积的概念
定义7.2 设f (x ),g (x ) ∈ C [a , b ],ρ (x )是[a , b ]上的权函数,则称
⎰=b
a
dx x g x f x g f )()()(),(ρ
为f (x )与g (x )在[a , b ]上以ρ (x )为权函数的内积。 内积有如下性质:
(1) (f , f )≥0,且(f , f )=0 ⇔ f = 0; (2) (f , g ) = (g , f );
(3) (f 1 + f 2, g ) = (f 1, g ) + (f 2, g ); (4) 对任意实数k ,(kf , g ) = k (f , g )。 这些性质,由内积的定义不难得到证明。 3.正交性的概念
定义7.3 设f (x ),g (x ) ∈C [a , b ]若
⎰==b
a
dx x g x f x g f 0)()()(),(ρ
则称f (x )与g (x )在[a , b ]上带权ρ (x )正交。
定义7.4 设在[a , b ]上给定函数系{
} ),(,),(),(10x x x n ϕϕϕ,若满足条件
())
(),1,0,(,0,0)(),((是常数k k k
j A k j k
j A k
j x x ⎩⎨⎧==>≠= ϕϕ
则称函数系{ϕk (x )}是[a , b ]上带权ρ (x )的正交函数系,特别地,当A k ≡ 1时,则称该函数系为标准正交函数系。若定义7.4中的函数系为多项式函数系{} )(),(10x p x p ,则称{}
)(x p k