计算方法 第五章第三节最优一致逼近
数值分析之最小二乘法与最佳一致逼近

这里 ( x) 0是 [a, b]上的权函数,它表示不同点 ( xi , f ( xi ))
处的数据比重不同.
5
S ( x) a00 ( x) a11 ( x) ann ( x) (n m) 用最小二乘法求拟合曲线的问题,就是在
S ( x )中求一函数 y S * ( x), 使误差取得最小.
23
结果如下:
24
2
用正交多项式做最小二乘拟合
用最小二乘法得到的法方程组,其系数矩阵
n . G是病态的
(
(k 0,1,, n). (5.6) j 0 如果 0 ( x), 1 ( x), , n ( x) 是关于点集
k
, j )a j d k
{xi } (i 0,1,, m) 带权 ( xi ) (i 0,1,, m)
4 2.00 8.46 2.135
21
( 0 , y ) yi 9.404,
i 0 4
4
(1 , y ) xi yi 14.422.
i 0
故有法方程
5 A 7.50b 9.404, 7.50 A 11.875b 14.422.
解得
A 1.122, b 0.505, a e A 3.071.
使误差平方和
* 2 [ S ( x ) y ] i i i 0 2 i i 0 m m m
min
S ( x )
2 [ S ( x ) y ] , i i i 0
这里
S ( x) a00 ( x) a11 ( x) ann ( x)
(n m).
y S * ( x) 与所给数据{( xi , yi ), i 0,1,, m} 拟合.
计算方法最佳一致逼近多项式切比雪夫多项式ppt课件

0,1,2,… , n)
轮流取得最大值1和最小值 1,{xk }称为交错点组。
- 1 x4
x 3
x2 0
x 1
x0 1
证: 将xk
cos
kπ n
,
(k
1,2,… , n)
代入Tn(x)的表达式,得到
Tn(x)
cos[narccos(cos
kπ )] n
cos[kπ]
(1)k
1
T2(x) T1(x)
多项式,且 max | f(x)
-1 x 1
Ln(x)
|
1 2n(n 1)!
||
f (n 1) (x)
||
证明:
max
-1 x 1
|
f(x)
Ln(x)
|
(n
1 1)!
||
f(n1)(x)
||||
(x
x0 )(x
x1) … (x
xn)
||
(n
1 1)!
||
f (n 1) (x)
||||
1 2n
xn )
|
要使 max 1 x 1
|
(x
x0 )(x
x1) … (x
xn )
|
取极小值, 只需令:
(x x0 )(x x1) … (x xn)
1 2n
Tn1(x),
最佳一致 逼近0的 多项式
而上式成立的充分必要条件是x0, x1,…xn是切比雪夫 多项式的0点。
将Lagrange插值多项式Ln(x)的节点取为Tn1(x) 的0点 :
cos[(2k
1)π] 2
0 (k
1,2, … , n)
图为T11(x)的零点,一共有11个
计算方法最佳一致逼近多项式切比雪夫多项式

路漫漫其修远兮, 吾将上下而求索
2020年4月11日星期六
内容
1. 函数逼近的基本概念 2. 切比雪夫多项式 3. 最佳一致逼近多项式 4. 切比雪夫多项式在函数逼近中的应用 5. 利用切比雪夫多项式的0点构造最佳逼近多
项式的例子
路漫漫其修远兮, 吾将上下而求索
y
y=L (x)
路漫漫其修远兮, 吾将上下而求索
一致逼近的几何意义
x Home
切比雪夫多项式
路漫漫其修远兮, 吾将上下而求索
切比雪夫(Chebyshev)多项式
• 切比雪夫多项式在逼近理论中有重要的应用 • 。切比雪夫多项式的0点可以用于构造具有最佳
一致逼近性质的插值多项式。
切比雪夫多项式的(简单)定义:
三、切比雪夫多项式在函数逼近中的应用
希望构造最高次幂xn 系数为1 的多项式:
路漫漫其修远兮, 吾将上下而求索
…
三、切比雪夫多项式在函数逼近中的应用
证明比较复杂,省略。
路漫漫其修远兮, 吾将上下而求索
这个定理的 结论非常重要
怎样才能使得拉格朗日插值多项式成为最佳逼近 ?
…
偏差估计
路漫漫其修远兮, 吾将上下而求索
吾将上下而求索
(5)切比雪夫多项式的极值点 …
路漫漫其修远兮, 吾将上下而求索
…
1
T2(x )
T1(x
)
-1
1
T3(x ) 路漫漫其修远兮,
吾将上下而求索
T4(x )
-1
T3(x)有3个0值点,4个极值点
总结: Tn(x)具有很好的性质。
y
x
Tn(x)是n阶多项式,具有n个0点,n+1个极值点;有 界[-1, 1]; T1(x), T3(x),…只含x的奇次项,是奇函数
数值分析(22)连续函数的最佳一致逼近.ppt

故
e1 e0
b
e 1 1.7183
10
由 f '( x1 ) e x1 e 1
求出 x1 ln(e 1) 0.5413
e0 e0.5413
0.5413
a
1.7183
0.[0,1]上的最佳一致逼近多项式为
P( x) 0.8940 1.7183 x
22
定理2(Chebyshev定理) 设f ( x) C[a, b], Pn ( x) H,则Pn( x)是f ( x)的
最佳一致逼近的充分必要条件是f ( x) Pn ( x)在[a, b] 上至少有n 2交错点组成的交错点组。
对n 1, f ( x) P1( x)有n 2 3个交错点。
a0
f (a) f ( x1 ) 2
f
(b)
f (a)
a
x1
ba
2
这样就得到f ( x)的线性最佳一致逼近多项式为
P1( x) a0 a1 x
数值分析
数值分析
例:选取常数a, b,使max | ex (a bx) | 达到最小。 0 x1
解:设P( x) a bx为f ( x) e x在[0,1]上的最佳一致 逼近多项式。
数值分析
数值分析
定义 设函数f ( x) C[a, b], 称点集
{ xk }kn0 { x0, x1, xn } 是f ( x)在[a, b]的交错点组,当且仅当满足
f ( xk ) (1)k
f (x)
(k 0,1, 2
, n)
其中 取1或 1。
例 f ( x) sin x 在[0, 2]的交错点组{1 , 3}。
x[ 1,1]
T n( x)
计算方法最佳平方逼近-最小二乘法

i0
i0
为最小。
这种要求误差(偏差)平方和最小的拟合称为 曲线拟合的最小二乘法。
(1)直线拟合
设已知数据点 (xi , yi ), i 1, 2, … , m 分布大致为 一条直线。作拟合直线 y(x) a0 a1x, 该直线
不是通过所有的数据点 (xi , yi ) ,而是使偏差平方和
F(a0 , a1 )
计算方法 (Numerical Analysis)
第5次 最佳平方逼近与曲线拟合的最小二乘法
主要内容
• 最佳平方逼近 • 曲线拟合的最小二乘法
最佳平方逼近
函数逼近的类型
• 最佳一致逼近:使用多项式对连续函数进行一致 逼近。逼近误差使用范数
||
f(x)
-
s(x)
||
max
a x b
|
f(x)
-
1(1 x2 )dx 1(0.934 0.426x) 1 x2dx
0
0
0.0026
最大值误差 :
同学们自己求一下
|| δ(x) || max | 1 x2 (0.934 0.426x) | 0.066
例题 求f(x) x在[1/4,1]上的在Φ span{1, x} 中的关于ρ(x) 1的最佳平方逼近多项式。
10 27
88 x 135
平方误差 :|| δ(x) ||22
1xdx
1
( 10 27
7 12
31 80
) 88
135
4
1.02
p1* (x)
10 27
88 x. 135
1
f(x) x
平方误差 : || δ(x) ||22 0.0001082.
最佳一致逼近多项式

§3最佳一致逼近多项式2-1 最佳一致逼近多项式的存在性切比雪夫从另一观点研究一致逼近问题,他不让多项式次数n 趋于无穷,而是固定n ,记次数小于等于n 的多项式集合为n H ,显然],[b a C H n ⊂。
记{1,,,}n n H span x x =L , n x x ,,,1L 是],[b a 上一组线性无关的函数组,是n H 中的一组基。
n H 中的元素)(x P n 可表示为01()n n n P x a a x a x =+++L ,其中n a a a ,,,10L 为任意实数。
要在n H 中求)(*x P n 逼近],[)(b a C x f ∈,使其误差)()(max min )()(max *x P x f x P x f n bx a H P n b x a n n −=−≤≤∈≤≤ 这就是通常所谓最佳一致逼近或切比雪夫逼近问题。
为了说明这一概念,先给出以下定义。
定义1 ],[)(,)(b a C x f H x P n n ∈∈,称)()(max ),(x P x f P f P f n bx a nn −=−=∆≤≤∞ 为)(x f 与)(x P n 在],[b a 上的偏差。
显然),(,0),(n n P f P f ∆≥∆的全体组成一个集合,记为)},({n P f ∆,它有下界0。
若记集合的下确界为,)()(max inf )},({inf x P x f P f E n b x a H P n H P n n n n n −=∆=≤≤∈∈ 则称之为)(x f 在],[b a 上最小偏差。
定义2 假定],[)(b a C x f ∈,若存在n n H x P ∈)(*,n n E P f =∆),(*, 则称)(*x P n 是)(x f 在],[b a 上的最佳一致逼近多项式或最小偏差逼近多项式,简称最佳逼近多项式。
注意,定义并未说明最佳逼近多项式是否存在,但可证明下面的存在定理。
最佳一致和平方逼近ppt课件

7
三、 Ca,b 上的最佳一致逼近的特征
引理4.1
设 f x 是区间a,b 上的连续函数,Pn* x 是 f x 的n次最佳一致逼近多项式,则 f x Pn* x 必同时
min f
x Pn* x
Pn xHn
f x Pn x
其中,H n代表由全体代数多项式构成的集合。
4
§2 最佳一致逼近多项式
一、最佳一致逼近多项式的存在性
定理4.9
对任意的 f xCa,b, 在 H n 中都存在对 f x 的最佳一致逼近多项式,记为 pn* x ,使得
f (x)
存在正负偏差点。
8
y
Oa
y f x En
y f x
y f x En
bx
9
定理 4.10( Chebyshev定理)
设 f x 是区间 a,b 上的连续函数,则 Pn* x 是 f x 的n次最佳一致逼近多项式的充要条件是: f x Pn* x 在区间a,b 上存在一个至少有 n 2 个交错偏差点组成,
注: 显然, f , Pn 0 , f , Pn 的全体组成一个
集合,记作 f , Pn ,它有下界0。
6
2、偏差点
定义
设 f xCa,b, PxHn, 若在 x x0 上有
P x0 f x0 max P x f x , a xb
则称 x0 是 P x f (x) 的偏差点。
由推论1,f x P1 x 在 a,b 上恰好有3个点构成的交错
组,且区间端点 a, b 属于这个交错点组,设另一个交错点为 x2 ,
6.2最佳一致逼近

注: 该性质又被称为Chebyshev多项式的最小模性质.
利用这个极值性质,切比雪夫多面式就成为[-1, 1]上逼 近其它函数f (x)的重要工具。 2.切比雪夫多项式在函数逼近中的应用 下面分两种情况讨论
多项式
2、 C a, b上最佳一致逼近多项式的存在性
定理6.2(Borel定理)
对任意的 f x C a, b , 在
Hn
中都存在对
* f x 的最佳一致逼近多项式,记为 pn x ,使得
f x pn* x
pn x H n
inf
且
引理6.2
limBn ( f , x ) f ( x ),
n
x [0,1] 一致成立
引理6.3 若函数r阶可导
limBn
n (m)
( f , x ) f ( m ) ( x ), x [0,1]
m<=r
引理6.4 如果函数 f C[0,1] 并且是凸函数,则
Bn ( f , x)
6.2 最佳一致逼近
最佳一致逼近
1、算例
e
x
e
1 2 1 x x 2
0.8
1 x
最佳一致逼近
1 2 1 3 1 4 1 x x x x 2 6 24
1 2 1 3 1 x x x 2 6
最佳一致逼近
e
1 2 1 x x 2
x
1 x
最佳一致逼近
1 2 1 3 1 4 1 x x x x 2 6 24 1 2 1 3 1 x x x 2 6
也是凸函数
6.2.2 最佳一致逼近多项式 相关的概念 定义6.8 若 P x Hn , f x C a, b, n
数值分析(22)连续函数的最佳一致逼近

插值逼近的性质
插值逼近的误差
插值逼近的误差取决于插值多项式的阶数和插值点的选择,一般 来说,阶数越高,误差越小。
插值逼近的稳定性
插值逼近的稳定性取决于插值多项式的选择和计算方法,选择合适 的插值多项式和计算方法可以提高稳定性。
插值逼近的应用
插值逼近在数值分析、数学建模、信号处理等领域近
多项式逼近是一种常用的逼近方法,通 过将函数表示为一系列多项式的和,来 逼近原函数。多项式逼近具有精度高、 适用范围广等优点,但计算量大、稳定 性差。
VS
插值法
插值法是一种常用的多项式逼近方法,通 过构造一个多项式来逼近原函数。插值法 具有数学基础扎实、计算稳定等优点,但 需要解决插值节点过多导致计算量大、数 值不稳定性等问题。
最佳一致逼近的误差通常用范数表示,常用的范数有L∞范数、 L2范数和L1范数等。
逼近的数学模型
01
逼近问题通常可以转化为求解一个 泛函极值问题,即寻找一个多项式 p(x),使得它在给定区间[a, b]上与 目标函数f(x)的误差最小。
02
逼近问题的数学模型可以表示为 求解一个极值条件下的优化问题 ,常用的方法有梯度法、牛顿法 、拟牛顿法等。
深入研究逼近定理
进一步探索逼近定理的内在机制,为逼近理论的 发展提供理论支持。
逼近误差分析
对逼近误差进行深入分析,建立更加精确的逼近 误差估计,提高逼近精度。
推广逼近理论
将逼近理论应用于更广泛的领域,如微分方程、 积分方程等,推动相关领域的发展。
逼近在实际问题中的应用拓展
数值计算
利用最佳一致逼近方法进行数值计算,提高计算精度和效率。
CHAPTER
最佳一致逼近的方法
线性逼近的方法
最佳一致逼近

主讲 孟纯军
函数逼近和函数空间
回忆一下向量空间的定义. 多项式空间 C[a,b], 连续函数空间
定义:设S是线性空间,x1,...., xn S 若存在不全为零的数a1,...., an ,使得 a1x1 an xn 0 称x1,...., xn线性相关,否则,线性无关。
||
f
(x)
pˆ n (x) ||
min
p( x)n
||
f
(x)
pn (x) ||
其中 n 表示次数不超过n的多项式全体。
称pˆn (x)为f (x)在[a,b]的最佳逼近n次多项式。
最佳逼近多项式一定存在。
定义:给定f (x) C[a,b], p(x) n
若在x0 [a,b]处有:
函数的内积
函数空间C[a,b], (x)为给定的权函数,
对任何f (x), g(x) C[a,b],
b
( f (x), g(x)) a (x) f (x)g(x)dx
为函数f (x), g(x)的内积。
由函数的内积导出范数:
1
|| f (x) ||2 ( f (x), f (x))2
为u1,, un线性无关。
证明:设k1,, kn为n个数,则 u1 ,, un线性无关等价于 k1u1 knun 0 (1) 只有零解,即k1 kn 0
将方程(1)两边用ui做内积,得到
(u1, u1) (u1, u2 ) (u1, un ) k1 0
0.0571 0.0604
0
其中第一列为自变量x 的值。
4.1091 2.1951
0 0
1.3811 0
最佳一致逼近多项式

( f , p n ) E n,
*
( 3 .3 ) 或
则称 p n ( x ) 是 f ( x ) 在 [ a , b ]上的 n 次 最佳一致逼近多项式 最小偏差逼近多项式 ,简称 最佳逼近多项式
*
*
.
定理 2 若 f ( x ) C [ a , b ],
*
则总存在 p n ( x ) H n , 使得 。
证明:令 ( x ) | P ( x ) f ( x ) |, 则 ( x ) 连续,因而可以达到最 即存在 x 0 , 使得 ( x 0 ) max ( x ) || P ( x ) f ( x ) || 。
a xb
大值,
这说明 x 0 是 P ( x ) 的一个偏差点,不妨设 由于 P ( x ) 是最佳逼近多项式,则
三、最佳一致逼近多项式
1.零次最佳一致逼近多项式 对于n=0的P0(x)有: P0(x) =(M+m)/2 其中M、m分别为f (x) 的最大值和最小值。 ∵f(x)C[a,b],由闭区间上连续函数性质;在[a,b]上存在两点x1,x2 使f (x1)=M, f (x2)=m, 即:x1,x2为偏差点(负,正)使:
axb
f (x)
n
(x)
即在H中 (x)与f(x)之差的绝对值的最大值是最小的,H中 任一ψ (x)与f(x)之差的绝对值都比它大,这样的 (x)为 f(x)在H中的最佳一致逼近函数。
定义1
设 f ( x ) C [ a , b ],
pn ( x ) H n , 称
a xb
逼近多项式
推论2 设f(x)C[a,b],则f(x)在Hn中的最佳一致逼近多项 式Pn(x),就是f (x)在[a,b]上的某个n次Lagrange插 值多项式。 证明∵Pn(x)有n+2个偏差点,亦即使f (x) -Pn (x)在[a,b]上至少 有n+2个点交替换正负号,亦就是说f(x) Pn(x)=0在[a,b]上有n+1 个根存在n+1个点:a x0<…< xn b使f (xi) Pn (xi)=0 即:f (xi)=Pn(xi) (i =0,1,2,…,n) , 所以,以此作为插值条件可得 到Pn(x),因此,Pn(x)就是以x0,x1,…,xn为插值节点的n次值多项 式。 切比雪夫定理不仅给出了最佳一致逼近多项式的特征, 并从理论上给出了寻找最佳一致逼近多项式的方法:
数值分析5.3 最佳平方及一致逼近

为 f(x) 与Pn(x)在[a, b]上的偏差. 显然 ( f , Pn ) 0 , ( f , Pn ) 的全体组成一个集合,记 为{ ( f , Pn ) },它有下界0. 若记集合的下确界为
*5.4 最佳平方逼近
*5.4.1 最佳平方逼近及其计算
最佳平方逼近问题:
对f(x)∈C[a, b]及C[a, b]中的一个子集
C [a , b ]
span{ 0 ( x ), 1 ( x ), n ( x )}
若存在函数 S ( x ) ,使下式成立
f ( x ) S ( x ) min f ( x ) S ( x ) 2 min
注2:由于勒让德多项式{Pk(x)}是在[-1, 1]上由多项
式基{1,x,…,xn,…}正交化得到的,因此利用函数的勒让
德展开得到最佳平方逼近多项式与由
S ( x) a a x a x .
n
0
1
n
直接通过解法方程得到Hn中的最佳平方逼近多项式
是一致的,但是当n较大时法方程出现病态,计算误
1 2
得法方程组
解出
1 1 2
1 2 1 3
a0 1.147 , a 1 0.609
a0 0.934, a1 0.426,
故得所求的一次最佳平方逼近多项式为 平方误差为 2 ( x ) 2 ( f , f ) ( S1 , f )
1 k j
( f ( x ), k ( x )) f ( x )x k dx d k
0
1
若用H表示 Gnn G( j ,k )nn 对应的矩阵,即
3.3最佳一致逼近多项式

( 1, k 1,2,3).
12
由于 f ( x) 在 [a, b]上不变号, 故 f ( x) 单调, f ( x) a1 在 (a, b)内只有一个零点,记为 x2, 于是
且点 xk cos
k (k 0,1, , n) 是 Tn ( x)的切比雪夫交错点组, n
8
由定理5可知,区间 [1, 1] 上 x n 在 H n 1 中最佳逼近多项式
为 Pn*1 ( x), 即 ( x) 是与零的偏差最小的多项式. n
定理得证.
9
例3 求 f ( x) 2 x3 x 2 2 x 1 在 [1, 1]上的最佳2次逼 近多项式. 解 由题意,所求最佳逼近多项式 P2* ( x) 应满足
5
使 P( x), 用反证法,若存在 Q( x) H n , Q( x)
f ( x) Q( x)
f ( x ) P( x ) .
由于
P( x) Q( x) [ P( x) f ( x)] [Q( x) f ( x)]
在点 x1 , x2 ,, xn 2 上的符号与 P( xk ) f ( xk )(k 1,, n 2) 一致, 故 P( x) Q( x) 也在 n 2 个点上轮流取“+”、“-”号. 由连续函数性质,它在 [a, b] 内有 n 1 个零点,但因
于是得 1 x 2 的最佳一次逼近多项式为
P 1 ( x) 0.955 0.414 x,
即
1 x 2 0.955 0.414 x, 0 x 1;
最佳一致和平方逼近

~ 由于首项系数为1的 n + 1次Chebyshev多项式 Tn+1 ( x) 由于首项系数为1 Chebyshev多项式
无穷范数最小, 故有 无穷范数最小,
f ( x ) − Pn ( x ) ~ = Tn+1 ( x) a0
于是 Pn ( x ) = f ( x ) − a0 Tn+1 ( x)
a ≤ x ≤b
∆ ( f , Pn ) = f − Pn
= max f ( x ) − Pn ( x )
上的偏差 偏差。 为 f ( x ) 与 P ( x) 在 [ a, b ] 上的偏差。 n
注: 显然, ( f , Pn ) ≥ 0 ,{∆ ( f , Pn )} 的全体组成一个 显然, ∆
一、最佳一致逼近的概念
定义 对于任意 设函数 f ( x) 是区间 [ a, b ] 上的连续函数, 上的连续函数, 给定的 ∀ε > 0 ,如果存在多项式 P ( x ) ,使不等式
max f ( x ) − P ( x ) < ε
a ≤ x ≤b
成立, 则称多项式 P ( x ) 在区间 [ a , b ] 上一致逼近 成立, (或均匀逼近)于函数 f ( x ) 。 均匀逼近)
(1)定义 (1)定义 Chebyshev多项式 称 Tn = cos ( n arccos x ) , x ≤ 1 为n次Chebyshev多项式.
[注] 令 θ = arccos x , 则 cos θ = x
n 2 n−2 2 4 n−4 4 而 cos nθ = cos θ − Cn cos θ sin θ + Cn cos θ sin θ − L
y = f ( x)
数值分析(西北工业大学出版社) 最佳一致逼近

n 1,2,
性质2 : Tn ( x )是n次代数多项式
性质3 : Tn ( x )的最高次幂 n的系数为 n1 x 2
性质4 : | Tn ( x ) | 1, | x | 1
性质5 : n为奇数, Tn ( x )为奇函数 ; n为偶数, Tn ( x )为偶函数.
性质6: Tn ( x )是[1,上带权 ( x ) 1] 1 1 x2 的正交多项式
1 1 x
2
的正交多项式
证:由
1
Tn ( x ) cos n
n m
得到
1
( x )T ( x )T
0
( x )dx
1
1
Tm ( x )Tn ( x ) 1 x
2
dx
x cos
cos m cos n sin ( sin )d , m n 0 cos m cos nd , m n 0 2 0 0, m n
Tn ( x ) cos(n arccos x ), 1 x 1 Tn ( x ) cos n , 0
有了明确的认识。 在这里给出契贝晓夫多项式是为了解决前面提出的 最佳一致逼近问题:
max | f ( x ) ci* x i | max | f ( x ) ci x i |
2、正交多项式
关于多项式族:
gk ( x ) Ak x k Ak 1 x k 1 A1 x A0
如果存在权函数 ( x ), 使得
b
Ak 0, k 0,1,2,
mn 0, b 2 ( x ) gm ( x ) gn ( x )dx ( x ) gn ( x )dx 0, a a
03(2)-最佳一致逼近

§2 最佳一致逼近一、最佳一致逼近的概念定义3.10 设函数f (x )是区间[a , b ]上的连续函数,对于任意给定的ε >0,如果存在多项式p (x ),使不等式ε<-<<)()(max x p x f bx a 成立,则称多项式p (x )在区间[a , b ]上一致逼近(或均匀逼近)于函数f (x )。
那么,对于在区间[a , b ]上的连续函数f (x ),是否存在多项式p (x )一致逼近于f (x )呢?这个问题有许多人研究过。
德国数学家维尔斯特拉斯(Weierstrass)在1885年曾给出下述著名定理。
维尔斯特拉斯定理 若f (x )是区间[a , b ]上的连续函数,则对于任意ε >0,总存在多项式p (x ),使对一切a ≤x ≤b 有ε<-)()(x p x f证明从略。
维尔斯特拉斯定理表明,连续函数f (x )可以用多项式p (x )逼近到任意精确程度,但维尔斯特拉斯定理只在理论上肯定了闭区间上的连续函数可以用多项式以任意精确度来逼近,并没有给出确定逼近得最快的多项式的方法。
事实上,如果精确度要求较高,则用来逼近的多项式的次数一般也很高,这就增加了计算工作量。
因而,在实际计算时,我们总量希望在一定的精确度要求下,逼近多项式的次数越低越好。
切比雪夫从这样的观点去研究一致逼近问题,他不让逼近多项式的次数n 趋于无穷大,而是先把n 加以固定。
对于给定的[a , b ]上的连续函数f (x ),他提出在次数不超过n 的多项式的集合p n 中去寻找一个多项式)(*x p n ,使它在[a , b ]上“最佳地逼近”f (x )。
这里最佳逼近的意思是指)(*x p n 对f (x )的偏差。
)()(max *x p x f n bx a -<< 和其它任一p (x ) ∈ p n 对f (x )的偏差)()(max x p x f bx a -<<比较时是最小的,也就是说{})()(max min )()(max )(*x P x f x p x f bx a p x p n b x a n-=-<<∈<<(3.18)这就是通常所谓的最佳一致逼近问题,也称为切比雪夫逼近问题。
计算方法 第五章第三节最优一致逼近

次最优一致逼近 多项式
y P (x) n
可见Pn(x) 是 f (x)的 某一个插 值多项式
0
y f (x)
y f (x)
y f (x)
x
由切比雪夫定理可推出: Pn(x) f (x) 在定义域上
至少变号 n+1 次,故至少有n+1 个根。
二、切比雪夫多项式的性质
下面我们求最好的直线所满 足的直线方程。设该方程为
p(x) a0 a1x.
由图示知它与 f (x) tan1 x的误差 R(x) f (x) p(x)在 x 0,
x 和x 1处依顺序变号,且绝对值达到最大,即有
R(0) f (0) p(0) , R( ) f ( ) p( ) , R(1) f (1) p(1) .
上式等价于求实数 a0,a1,...,an,使得多项式
pn (x) a0 a1x ... an xn
满足
S( a0,a1,...,an ) :
f
pn
inf pPn
f p
定理5.3.2 (存在性定理)
f (x) C [a,b],总存在最优一致逼近多项式
满足
pn (x) c0 c1x ... cn xn
k 0,1,..., n 1,
其中 1或 1,称这样的点组为 p(x) f (x)的(切比雪夫)交错点组。
推论
设 f (x) Cn1 [a,b]且 f n1(x)在[a,b]上保号,p(x)
Pn span{1, x,.., xn}为 f (x)在区间[a,b]上的n次最优一致逼近多 项式,则 p(x) f (x)在区间[a,b]上恰好存在n+2个交错点,且
计算方法最佳一致逼近多项式-切比雪夫多项式专题培训课件

cos(nθθs )c in o(snθθ )si
co s1 () n θ cos(θ n s θi) nc(o θ n
Tn(x)在 1,[1]上有的 n个 零不 点同 xk co(s22k 1 n)π , (k1,2 …,,n)
证:将xk
cos(2k 1)π, 2n
(k
1,2,…, n)
代入Tn(x)的表达式,得到
Tn(x)
cos[narccso(cos(2k 1)π)] 2n
cos[(2k 1)π] 2
计算方法最佳一致
逼近多项式-切比 雪夫多项式
内容
1. 函数逼近的基本概念 2. 切比雪夫多项式 3. 最佳一致逼近多项式 4. 切比雪夫多项式在函数逼近中的应用 5. 利用切比雪夫多项式的0点构造最佳逼近多
项式的例子
函数逼近的基本概念
第3章 函数逼近与曲线拟合
§1 函数逼近的基本概念
一、函数逼近与函数空间
实际应用需要使用简单函数逼近已知复杂函数。
函数逼近问题: 对于函数类A中给定函的数
f(x),要求在另一类较简单便的于计算的函
数类
BA
B
A
中找一个函数p(,x使) p(x)与f(x的) 误差在某
种度量意义下达到最. 小
定1 理(Weaisesrf)s(若 tx rC ) [b a],则 ,ε0, 多项式 使p得 (x),
得知:情况a)如 为果 奇n数,则n2(xT)只含n的偶, 次方 Tn1(x)只含x的偶方 数, 次从而左n端 1(xT)只含x的偶; 次 情况b)如果n为,偶则数2xn(Tx)只含x的奇, 次Tn方 1(x) 只含x的奇次方,左从端而 nT1(x)只含x的奇次方
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
称 x0为 p(x) f (x)的一个“负”的偏差点。
直接构造最优一致逼近多项式的确比较困难,不妨换 个角度,先考察它应该具备的性质。有如下结论:
定理5.3.4 f (x) C[a,b],p(x) Pn span{1, x,.., xn}为 f (x)在区间[a,b] 上的n次最优一致逼近多项式,则 p(x) f (x)同时存在正负偏差点。
性质3.
xk
cos
2k 1 ,
2n
k 1, 2,..., n.
Tn (x) 1且Tn (x)在区间[1,1]上恰好有n 1个交错点
k
xk cos n , k 0,1, 2,..., n.
性质4.
设 pn (x)是首项系数为1的n次多项式,则
pn
21n
Tn
21n.
性质4称为切比雪夫多项式的极性,这种极性是我们构造近似 最优一致逼近的依据。
三、近似最优一致逼近多项式
(一)切比雪夫插值多项式
设 pn (x)是函数 f (x) Cn1[1,1]以区间[1,1]上n 1点x0,x1,..., xn 为节点的n次Lagrange多项式,则公式误差为
f (x) pn(x)
次最优一致逼近 多项式
y P (x) n
可见Pn(x) 是 f (x)的 某一个插 值多项式
0
y f (x)
y f (x)
y f (x)
x
由切比雪夫定理可推出: Pn(x) f (x) 在定义域上
至少变号 n+1 次,故至少有n+1 个根。
二、切比雪夫多项式的性质
以 x 0,1为节点作一次Lagrange插值: tan1 x x,
0.0493是否为在区间[0,1]上计算 4 其函误数差tRa2n(x1)x的ta最n好1 x近 似4 x一在次两式端呢点?附近很小,
但在x /4-1处达到最大值 0.0711.
用最佳平方一次多项式拟合: tan1 x 0.0429 0.7918x, 其误差的绝对值在x 1处达到最大值 0.0493.
f
( (n1) x
(n 1)!
)
(
x
x0
)(
xx1Fra bibliotek)...(
x
x
n
),
x [1,1].
记 M n1
max |
1 x1
f
(n1) (x) | .
则
|
f
(x)
pn (x) |
M n1 | (x (n 1)!
x0)(x
x1)...(x
xn) |.
由切比雪夫多项式的极性,当 x0,x1,...,xn取为n 1次切比雪夫 多项式的零点
| p(x0 ) f (x0) |
f p max | f (x) p(x) | axb
则称 x0为 p(x) f (x)的一个偏差点。 而且,若
p(x0 ) f (x0 ) ,
称 x0为 p(x) f (x)的一个“正”的偏差点。 若
p(x0 ) f (x0 ) ,
由于x p(x为 ) 误0.差03函56数R0(.7x)8的54极x,值点 ,0故.0还35应6.有
1 0.2146, 2 0.0711, 3 0.0493,
R '() f '() p '() 0.
0.0356.
整理得
a0 ,
a1 / 4 0.7854,
下面我们求最好的直线所满 足的直线方程。设该方程为
p(x) a0 a1x.
由图示知它与 f (x) tan1 x的误差 R(x) f (x) p(x)在 x 0,
x 和x 1处依顺序变号,且绝对值达到最大,即有
R(0) f (0) p(0) , R( ) f ( ) p( ) , R(1) f (1) p(1) .
S( c0,c1,...,cn ) :
f
pn
inf pPn
f p
定理5.3.3 (惟一性定理) f (x) C [a,b],其最优一致逼近多项式
pn (x) c0 c1x ... cn xn 是惟一的。
定义5.3.1 设 f (x) C[a,b],p(x) Pn span{1, x,.., xn},
k 0,1,..., n 1,
其中 1或 1,称这样的点组为 p(x) f (x)的(切比雪夫)交错点组。
推论
设 f (x) Cn1 [a,b]且 f n1(x)在[a,b]上保号,p(x)
Pn span{1, x,.., xn}为 f (x)在区间[a,b]上的n次最优一致逼近多 项式,则 p(x) f (x)在区间[a,b]上恰好存在n+2个交错点,且
n
f (x) akTk (x) Sn (x) akTk (x)
k=0
k=0
f (x) Sn (x) akTk (x) an T 1 n1(x) k=n+1
由于Tn1(x)有n 2个交错点,故可以近似地说,f (x) Sn (x) 也有n 2个交错点。于是根据切比雪夫定理,Sn (x)是函数 f (x)在区间[1,1]上的一种近似n次最优一致逼近多项式。
证 令 arc cos x, 则 x cos,
cos(n 1) cos(n 1) 2cos cos n,
Tn1(x) Tn1(x) 2xTn (x)
性质2.
Tn(x)为 n 次多项式,首项系数为 2n1,T2n(x)只含 x 的偶次幂, T2n+1(x)只含x 的奇次幂。且Tn(x)在区间 [0,1]上有 n个零点:
Runge 现象说明了这一点。如果用最佳平方逼近,多项式的摆 动现象也说明高次最佳平方逼近多项式拟合不一定就会达到好 的效果。那么,对于在区间 [a, b] 上连续的函数 f (x) ,是否必 存在多项式序列 {Pn(x)},使得在区间 [a, b] 上一致地逼近函数 f (x)呢?
定理5.3.1 维尔斯特拉斯(Weierstrass)定理
xˆk
1 2
(b a) cos
2k 1
2(n 1)
a
b ,
k 0,1,..., n.
作变换:x a b a b t. 22
例5.3.2 P192
(二)截断切比雪夫级数
设函数 f (x) C [a,b]. 当正交函数取为切比雪夫多项式族时,得
xk
cos 2k 1 ,
2(n 1)
k 0,1,..., n,
时,
max
1 x1
|
(x
x0
)( x
x1 )...( x
xn
)
|
达到最小值。 此时
| f (x) pn (x) | 2nM(nn11)!.
我们称这样的插值多项式为切比雪夫插值多项式。
当 f (x) C [a,b]时,切比雪夫插值多项式的插值节点变为
f p max | f (x) p(x) | axb
称为 f (x)与 p(x)的偏差。
因此, 最优一致逼近问题等价于求多项式 pn (x) c0 c1x ... cn xn ,
使之与函数 f (x)的偏差达到最小。
定义5.3.2 设 f (x) C[a,b],p(x) Pn span{1, x,.., xn},若在 x0 处满足
为 f (x)在区间[a,b]上的n次最优一致逼近多项式的充要条件是 p(x) f (x)
在区间[a, b] 上至少存在n+2个轮流取“正”、“负”的偏差点。即有n+2
个点 a x0 x1 ... xn1 b,使之成立
p(xk ) f (xk ) (1)k
f p ,
f : max | f (x) |, f (x) C[a,b] axb
称为无穷范数或者一致范数
设函数 f (x)在区间[a,b]上连续,pn (x) Pn span{1, x,.., xn}使得
f
pn
inf pPn
f p
成立,则称函数 pn (x)为函数 f (x)在区间[a,b]上的n次最优一致 逼近多项式。
2
(
f
,Tk )
2
1 f (x)Tk (x)dx
1
1 x2
2
f (cos ) cos k d ,
0
k 1, 2,...
因此,函数 f (x)的切比雪夫级数即为函数 f (cos )的Fourier级数。
当 f (x)在区间[-1,1]上分段光滑时,切比雪夫级数在此区间一致收敛于 f (x).
第五章
最佳逼近
§3. 最优一致逼近
一、 最优一致逼近的概念与求法
例5.3.1 试求一次多项式在区间 [0,1] 上逼近函数 f (x) tan1 x.
解 Taylor公式: 在零点展开,tan1 x x, 其误差R1(x) tan1 x x在x 0附近很小, 在x 1处达到最大值 0.2146.
上式等价于求实数 a0,a1,...,an,使得多项式
pn (x) a0 a1x ... an xn
满足
S( a0,a1,...,an ) :
f
pn
inf pPn
f p