第四章 最佳逼近
最佳逼近
若p S是最佳逼近元,则f p S,即
( f - g , 1 ( x)) 0, ( f - g , 2 ( x)) 0,..., ( f - g , m ( x))
称为法方程或者正规方程。
( f ( x ) - p( x )) ( x ) [ y (c ( x ) ... c
函数的最佳逼近
主讲 孟纯军
插值法是用多项式近似的表示函数,并要 求在他们的某些点处的值相拟合. 最佳逼近(或者曲线拟和)也是用简单 函数逼近复杂函数(或未知函数),但 是,逼近的原则和插值的原则不一样。
离散情形
最小二乘拟合直线 最小二乘拟合多项式 非线性拟合
Hilbert空间中的投影定理
1的基为1, x, 则g 1是f ( x)在最佳逼近元的充要条件为
(f g ,1)= ( f ( xi ) g ( xi )) 1 ( f ( xi ) (a bxi )) 1 0
i 1 i 1
n
n
(f g ,x)= ( f ( xi ) g ( xi )) xi ( f ( xi ) (a bxi )) xi 0
2 i 1
n
即如下最佳逼近问题:
1.子空间为 m,即次数不超过n的多项式, 取它的基函数为 1, x,...x ;
m
2. 在 m中找一个元素p( x),使它与给定函数 f ( x)最靠近,即 ( p( xi) f ( xi )) 2 min 。
i 1 n
p m是f ( x)在最佳逼近元的充要条件为
b=0.9068
最小二乘拟合直线为y= 0.0147 +0.9068x
最佳平方逼近
n
因为 f p*, p * p cj c j f p*, j 0 及
j0
( p * p, p * p) 0, 故 ( f p, f p) ( f p*, f p*).
2 则 f (x) 1 1 t g(t), 1 x 1
2 先求g(t)在区间 [-1,1] 旳一次最佳平方逼近多项式.
由
c0 *
1 2
(g,
L0 )
1 2
1 1
1 2
1 tdt 2 , 3
c1
3 2
(g,
L1 )
3 2
1 1
t 2
1 tdt 2 . 5
可知
2
2
22
q1(t) 3 L0 (x) 5 L1(x) 3 5 t,
例6 定义内积 ( f , g)
1
f (x)g(x)dx
0
试在H1=Span{1,x}中谋求对于f(x)= x 旳最佳平方逼近
元素p(x).
解 法方程为
1 12
1 2
13
c0
c1
2 2
3 5
解得
c0
4, 15
c1
4 5
所求的最佳平方逼近元素为 p(x) 4 4 x. 0 x 1 15 5
(n
,n
)
cn*
( f ,n )
因为 0x,1x, , nx 线性无关, 能够推得上系数阵是
非奇异旳. 故 (5. 82) 有唯一解 { c*j }.
四、最佳平方逼近旳误差
记 ( f p*, f p*), 称其为最佳平方逼近误差, 利用
第四章 3最佳平方逼近(1)
§3 最佳平方逼近摘要:介绍内积赋范线性空间中的最佳平方逼近的特征,最佳逼近元存在唯一性及求解;2,L a b ρ⎡⎤⎣⎦中的最佳平方逼近问题的讨论。
一、赋范线性空间定义1 线性空间X 称为赋范线性空间,如果其上赋予一个满足如下3条性质的函数:(1)0,00,;x x x x X ≥=⇔=∈ (2),,;α=∈∈ax a x x X(3),,.x y x y x y X +≤+∈ 则称是线性空间X 的范数。
例1 欧氏空间(欧几里德空间(Euclid ))n:12221,,=⎛⎫=∈⎪⎝⎭∑nnk k x x x范数称为欧氏范数或2-范数。
例2 1-范数赋范线性空间n:11,,==∈∑nnk k x x x称为1-范数。
定义2 赋范线性空间中的最佳逼近:若Y 是赋范线性空间X 的一个线性子空间,x X ∈,则称量(),inf y Yx Y x y∈∆=-为子空间Y 对元素x 的最佳逼近,而使上式成立的元素*y 称为最佳逼近元,且Y 称为逼近子空间。
二、内积空间定义3 假设X 是一线性空间,如果其上赋予一个满足如下4条性质的二元函数(),:()()(1)(,)(,),,;(2)(,)(,),,,;(3)(,)(,)(,),,,;(4),0,;,00,ααα=∀∈=∀∈∈+=+∀∈≥∀∈=⇔=x y y x x y X x y x y x y X x y z x z y z x y z X x x x X x x x则称X 为内积空间。
例3 欧几里得空间n: (),,,=∈Tnx y x y x y内积→范数:2x ,2x 满足范数的3条性质。
内积空间→赋范线性空间定义4 内积空间中的最佳逼近:假设(1,2,,)ϕ=i i n 是内积空间X 中的n 个线性无关的元素,f X ∈,则子集{}12,,,ϕϕϕΦ=n n span对f 的最佳平方逼近定义为()2,min ϕϕ∈Φ∆Φ=-nn f f . (1) 使(1)成立的那个元素称为最佳逼近元素。
高代课程论文——最佳逼近
高等代数课程论文——最佳逼近:学号:班级:一.摘要欧几里德空间(Euclidean Space) 简称为欧氏空间,在数学中是对欧几里德所研究的2维和3维空间的一般化。
这个一般化把欧几里德对于距离、以及相关的概念长度和角度,转换成任意数维的坐标系。
这是有限维、实和积空间的“标准”例子。
欧氏空间是一个特别的度量空间,它使得我们能够对其的拓扑性质,例如紧性加以调查。
积空间是对欧氏空间的一般化。
积空间和度量空间都在泛函分析中得到了探讨。
欧几里德空间在对包含了欧氏几何和非欧几何的流形的定义上发挥了作用。
一个定义距离函数的数学动机是为了定义空间中围绕点的开球。
这一基本的概念正当化了在欧氏空间和其他流形之间的微分。
微分几何把微分,会同导入机动性手法,局部欧氏空间,探讨了非欧氏流形的许多性质。
当一个线性空间定义了积运算之后它就成为了欧几里德空间。
二.关键字欧式空间最佳逼近函数构造三、问题的阐述1、欧式空间的定义设V是实数域R上的线性空间或称为向量空间,若V上定义着正定对称双线性型g(g称为积),则V称为(对于g的)积空间或欧几里德空间(有时仅当V是有限维时,才称为欧几里德空间)。
具体来说,g是V上的二元实值函数,满足如下关系:(1)g(x,y)=g(y,x);(2)g(x+y,z)=g(x,z)+g(y,z);(3)g(kx,y)=kg(x,y);(4)g(x,x)>=0,而且g(x,x)=0当且仅当x=0时成立。
这里x,y,z是V中任意向量,k是任意实数(1)<ξ,η>=<η,ξ>(2)<ξ+η,ζ>= <ξ,ζ>+<η, ζ>(3)<aξ,η>= a<ξ,η>(4)当ξ≠0 时<ξ,ξ>0这里ξ,ζ,η是V的任意向量a是任意实数,那么V叫作对这个积来说的一个欧几里的空间,简称欧式空间。
2、举例说明例1:在Rn里对于任意两个向量ξ=(x1,x2,…,Xn)η=(y1,y2,…,Yn)规定<ξ,η>=x1y1+x2y2+…+xnyn容易验证关于积的公理被满足,因而Rn对于这样定义的积来说作成一个欧式空间。
4章§3 最佳平方逼近
定理6
ϕ0 (x),ϕ1(x),Lϕn−1(x), 在[a,b]上线性无关的充分必要条
件是它的克来姆(Gramer)行列式 G ≠ 0, ,其中 n−1
(ϕ0 ,ϕ0 ) Gn−1 = G(ϕ0,ϕ1, Lϕn−1 = (ϕ1,ϕ0 ) L
(ϕ0 ,ϕ1) L (ϕ0 ,ϕn−1) (ϕ1,ϕ1) L (ϕ1,ϕn−1) L L L
用{1,x,…,xn)做基,求最佳平方逼近多项式,当n较大时,系数 矩阵(3.16)式是高度病态的(病态矩阵概念见第七章),求法方程 (3.13)的解,舍人误差很大,这时要用正交多项式做基,才能求得 最小平方逼近多项式(见§5).
(ϕn−1,ϕ0 ) (ϕn−1,ϕ1) L (ϕn−1,ϕn−1)
(3.10) 定理的证明由读者完成.
一、函数的最佳平方逼近
现在我们研究在区间[a,b]上一般的最佳平方逼近问题. 对 f (x) ∈C[a,b]及C[a,b]中的一个子集ϕ = span{ 0 ,ϕ1 ,L,ϕn ) , ϕ 若存在S*(x) ∈ϕ ,使
(3.2)
则在(a,b)上g(x)≡0,就称 ρ(x) 为区间(a,b)上的权函数. 定义5 设f(x),g(x)∈[a,b],是[a,b]上的权函数,积分
( f , g) = ∫ ρ(x) f (x)g( x)dx
a
b
称为函数f(x)与g(x)在[a,b]上的内积.
容易验证这样定义的内积满足下列四条公理; 1) ( f , g)=(g, f ) ; 2) (cf , g)=c(f , g) ,c为常数; 3)
于是有
∑(ϕ ,ϕ )a
j=0 k, j
n
j
= ( f ,ϕk )
(k = 0,1,L, n)
最佳平方逼近
所求的
应该使下式达极小:
由
整理得到
计算积分后,得法方程组
解之得 从而得到最佳平方逼近一次多项式
三、正交基函数的选择 如果我们选择子空间
正交,即 则法方程
简化为
即 容易求得 并得到最佳平方逼近
例3.2. 已知
在区间 [-1,1]上两两正交,试求
在这个
区间上的最佳平方逼近二次多项式,并给出误差估计。
应该使
整体达最小。 通过这种度量标准求得拟合曲线y=f(x)的方法,
就称作曲线拟合的最小二乘法。 按照以上思想来求出f(x)的拟合曲线,首先需
要确定出f(x)所属的函数类,然后进一步求出具体 函数,具体按照以下步骤进行。
二、最小二乘法拟合曲线的步骤
第一步:根据如下已知点的坐标,在坐标系里描点
第二步:根据图示,确定曲线所属的函数类型,例 如多项式函数类、三角函数类、指数函数 类、对数函数类等。假设所确定的函数类 的基函数为
而n+1元函数
在区间
上具有一阶连续导函数,因此根据
极值原理,在最小值点
处:
而 于是 即
利用内积 可以得到 这是一个含有n+1个变量的方程组,具体形式为:
再写成 矩阵形式为
这是关于n+1个变量
的线性方程组,并称
其为法方程组,或者正规方程组。
解此方程组,就可以得到 了f(x) 的最佳平方逼近:
,也就得到
同时,还需要给出连续函数
空间上的一个度量标准,下面通过内积给出平方范数。
二、连续函数的平方范数
已知所有连续函数构成的集合C[a,b]是一个线性
空间,对于C[a,b]中的任意函数
、 ,定义
实数
最佳逼近PPT课件.ppt
称线性空间Y为内积空间,(f,g)为内积。
有可能给的条件个数n大于多项式
P(x)的待定系数个数,如,10个插值条件求
5次多项式,该问题是无解的。
有时我们所需的近似函数不一定是多项式。
在实际问题中,往往并不要求近似函
数φ(x)所表示的曲线通过这些观测点,而只
要求由已知数据(xi,yi)(i=0,1,…,n)找出x,y之 间的依赖关系,使得近似函数φ(x)能充分地
a0(n 1) a1 n xi
i0
a0
n i0
n
xi a1
i0
xi2
n
n
am xim
yi
i0
i0
n
n
am
x m 1 i
xi yi
i0
i0
a0
n i0
n
xim a2
i0
x m 1 i
n
n
m
xi2m
xim yi
i0
i0
(4―73)
写成矩阵形式
(4-74)
例8 设有一组数据表
故在基点x0,x1,x2,…,xn上φ(x)与f(x)有
误差
ri=φ(xi)-yi, i=0,1,2,…,n
(4―69)
称ri为用φ(x)拟合f(x)的偏差。
设 函 数 关 系 y=f(x) 的 一 组 观 测 数 据 为
(xi,yi)(i=0,1,2,…,n),欲求一个m(m<n)次多项 式
Pm(x)=α0+α1x+…+αmxm
最佳逼近与最优求积问题
最佳靠近与最优求积问题是数值分析中的经典问题,这些问题在工程学、科学和数学中都有广泛的应用。
本文将介绍最佳靠近和最优求积两个问题的基本观点与相关算法,包括泰勒展开、拉格朗日插值法、牛顿插值法、切比雪夫靠近、勒让德积分等。
本文还将探讨最佳靠近和最优求积的应用,比如在数据拟合、信号处理、图像压缩、插值、微积分求解等方面的应用。
最后,我们还将介绍最佳靠近和最优求积的优缺点及其将来的进步方向。
关键词:最佳靠近、最优求积、泰勒展开、拉格朗日插值法、牛顿插值法、切比雪夫靠近、勒让德积分、数据拟合、信号处理、图像压缩、插值、微积分求解一、引言是数值分析中的经典问题之一。
最佳靠近是指在数学上寻找一个函数或曲线,使其在某种范数下与实际数据或函数的距离最小。
最优求积则是指在一个区间内,寻找一个带权重的多项式函数,使其与实际函数的误差最小。
这两个问题屡屡在工程学、科学和数学中出现,并且有着广泛的应用。
本文将介绍最佳靠近和最优求积问题的基本观点与相关算法,以及它们在实际应用中的作用和不足之处。
二、最佳靠近最佳靠近的主要目标是在给定点集中找到一个函数,使得该函数与实际数据的差距最小。
在数学上,我们可以利用不同的靠近算法,比如泰勒展开、拉格朗日插值法、牛顿插值法、切比雪夫靠近等等。
其中,泰勒展开用于对光滑函数进行靠近,拉格朗日插值法和牛顿插值法则常用于数据拟合,而切比雪夫靠近则可用于对非光滑函数进行靠近。
三、最优求积最优求积的主要目标是在给定区间内,寻找一个带权重的多项式函数,使得其与实际函数的误差最小。
最优求积的解决方法浩繁,比较常见的方法包括勒让德积分、拉盖尔积分、傅里叶变换等等。
这些方法可用于求解微积分问题、信号处理、图像压缩等问题。
四、最佳靠近和最优求积的应用最佳靠近和最优求积在实际应用中被广泛使用。
它们可用于数据拟合、信号处理、图像压缩、插值、微积分求解等方面。
比如在数据拟合中,我们可以利用最小二乘法对数据进行拟合,得到一条最佳靠近的直线或曲线。
最佳一致逼近
构造C[0,1]上W=&(1,x,…,x9)到f(x)=e x上的最佳逼近学院:数学与计算机科学学院班级:2011级数学与应用数学姓名:学号:指导教师:日期:2012.06.20构造C[0,1]上W=&(1,x ,…,x 9)到f(x)=e x上的最佳逼近(西北民族大学数学与应用数学专业,兰州 730124)指导教师摘要: 本文通过对最佳逼近的研究,着重分析其构造方法,从而使得对知识的掌握更加连贯及牢固。
通过对它的研究,我们对其有了更深的了解。
关键词:最佳逼近,正射影,傅里叶系数最佳平方逼近一般而言,在[a , b ]上对给定的函数求它的一致逼近函数比较困难,下面我们介绍在[a , b ]上较易计算的另一种逼近方法――最佳平方逼近。
一、预备知识1.函数系的线性关系定义1 若函数)(,),(),(10x x x n ϕϕϕ ,在区间[a , b ]上连续,如果关系式0)()()()(221100=++++x a x a x a x a n n ϕϕϕϕ 当且仅当0210=====n a a a a 时才成立,则称函数)(,),(),(10x x x n ϕϕϕ 在[a , b ]上是线性无关的,否则称线性相关。
如果函数系{ϕk (x )}(k = 0, 1, 2, …)中的任何有限个函数线性无关,则称函数系{ϕk (x )}为线性无关函数系,例如{1, x , …, x n , …}就是在区间[a , b ]上的线性无关函数系。
设)(,),(),(10x x x n ϕϕϕ 是[a , b ]上线性无关的连续函数a 0, a 1, …, a n 是任意实数,则)()()()(1100x a x a x a x S n n ϕϕϕ+++=的全体是C [a , b ]的一个子集,记为},,,{Span 10n ϕϕϕ =Φ并称)(,),(),(10x x x n ϕϕϕ 是生成集合的一个基底。
最佳平方逼近
4 15 4 5 x. 0 x 1
所求的最佳平方逼近元
素为 p ( x )
5.6.2 正交系在最佳平方逼近中的应用
当 0x,1x, , nx, 是正交系时,求解最佳平方逼 近式(5.82)中的系数非常容易. 目标: 求下面的最佳平方逼近式中的系数
( 5.84)
函数 f 的 L-最佳平方逼近函数为
pL ( x )
n k 0
c k L k ( x ),
(L)
1 x 1
(5.85)
遇到区间[a,b], 通过下面的变换把问题转化到[-1,1]上处理.
x a b 2 b a 2 t.
函数 f (x) 的 Legendre 无穷级数
三、最佳平方逼近函数的求解
利用 ( f p*, j )=0, 可求出最佳平方逼近函数 p*. 设
p * (x)
n
故
0
ck k ( x)
*
n * k
k 0
f
p *,
j
f , c
j k 0 j
k
,
j
(5.82)
n
ck k ,
1
( f , g)
1
f ( x ) g ( x )dx
L-正交多项式为 L0x, L1x, , Lnx, 用(5.83), 有
ck
(L)
( f , Lk ) ( Lk , Lk )
2k 1 2
1
1
L
k
( x ) f ( x )dx ,
k 0,1, 2, ..., n
最佳一致逼近
主讲 孟纯军
函数逼近和函数空间
回忆一下向量空间的定义. 多项式空间 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
计算方法 最佳一致逼近多项式-切比雪夫多项式精编版
22
22
22
22
xk
cos (2k 1)π , (k 22
1,2,… ,11)
接近-1和1的地方越密。过这些0点作平行于y轴的直
线,这些直线与上半单位元的交点形成了一个关于圆
弧的等距的点的集合。
(5)切比雪夫多项式的极值点
Tn(x)在[1,1]上有n 1个不同的极值点
x k
cos kπ , (k n
π,
m n 0.
(2.12)
证:令x cosθ ,则
1
1
1
1
x2
Tm(x)Tn(x)dx
0 cos(mθ )cos(nθ )dcosθ
π
1 cos2θ
π cos(mθ )cos(nθ )dθ 0
根据积化
1 2
3)则对n 1的情况,由递推公式 Tn1(x) 2xTn(x) Tn1(x)
得知:情况a)如果n为奇数,则2xTn(x)只含n的偶次方, Tn1(x)只含x的偶数次方,从而左端Tn1(x)只含x的偶次方; 情况b)如果n为偶数,则2xTn(x)只含x的奇次方,Tn1(x) 只含x的奇次方,从而左端Tn1(x)只含x的奇次方
切比雪夫多项式的前几项:
T0(x) cos(0) 1 T1(x) cos(arccosx) x T2(x) cos(2arccosx) 2x2 1 T3(x) cos(3arccosx) 4x3 3x
课堂练习:推出T4(x)
切比雪夫多项式的性质
(1)基本递推关系
计算方法 (Numerical Analysis)
第4次 最佳一致逼近多项式
最佳逼近
n
n
这两个方程称为法方程,或者正规方程。
将法方程组整理,得到
n n na ( xi )b yi i 1 i 1 n n n ( x )a ( x 2 )b x y i ii i i 1 i 1 i 1
解法方程组,就可得到a, b, 从而得到最小二乘拟合直线y a bx.
线性最小二乘拟合
给定函数y f ( x)在点x1,...xn处的函数值y1,... yn ,
给定子空间S,设S的基函数为1 ( x),....m ( x)(m n),
求p( x) c11 ( x) .... cmm ( x) S ,
使得 ( p( xi ) yi ) min
解 子空间S的基函数为 2 , 拟合函数为a bx 2, 1,x
则法方程组为:
5 5 2 5a ( xk )b yk k 1 k 1 5 5 5 ( x 2 )a ( x 4 )b= x 2 y k k k k k 1 k 1 k 1
解法方程组,得到c1 ,..., cm ,
从而得到最小二乘拟合: p( x) c11 ( x) .... cmm ( x)
例 用y a bx 最小二乘拟合如下数据:
2
25 31 38 44 x 19 y 19.0 32.3 49.0 73.3 97.8
i 1 i i 1 i i 1 i 1 1 i m
n
n
m
( xi ))]1 ( xi ) 0
(f(x )-p(x )) ( x ) [y (c ( x ) ... c
i 1 i i 2 i i 1 i 1 1 i m
数值逼近课程设计最佳逼近
数值逼近课程设计最佳逼近一、课程目标知识目标:1. 让学生掌握数值逼近的基本概念,理解最佳逼近的原理及其在数值计算中的应用。
2. 使学生能够运用不同的数值方法进行数据逼近,并分析其优缺点。
3. 帮助学生建立误差分析的基本框架,培养他们评价逼近效果的能力。
技能目标:1. 培养学生运用数学软件或编程工具实现数值逼近算法的能力。
2. 提高学生解决实际问题时选择合适数值逼近方法的能力,并能进行相应的参数调优。
3. 培养学生通过团队协作,共同解决复杂数值计算问题的能力。
情感态度价值观目标:1. 培养学生对数学科学的兴趣,特别是对数值逼近这一领域的探索热情。
2. 增强学生的实证思维,培养他们严谨的科学态度和精益求精的学术追求。
3. 通过数学建模和问题解决,激发学生的创新意识,增强他们面对挑战时的自信心和坚持到底的决心。
本课程设计针对高中年级学生,考虑到他们已具备一定的数学基础和逻辑思维能力,课程性质为理论与实践相结合。
在教学过程中,要求教师注重启发式教学,鼓励学生主动探究和动手实践,通过案例分析、小组讨论等形式,提高学生的问题解决能力和团队合作精神。
课程目标的设定旨在使学生不仅掌握数值逼近的相关知识,而且能够将这些知识应用于实际问题中,培养他们的综合素养。
二、教学内容本章节教学内容围绕以下三个方面展开:1. 数值逼近基本概念:- 介绍逼近的概念、逼近的误差及度量方法。
- 解释最佳逼近的定义及其判定标准。
2. 数值逼近方法:- 分析常用的数值逼近方法,如插值法、最小二乘法、样条插值等。
- 详述各种方法的原理、步骤和适用范围。
教学大纲:a. 插值法:Lagrange插值、Newton插值、Hermite插值等。
b. 最小二乘法:线性最小二乘法及其应用。
c. 样条插值:线性样条、二次样条和三次样条插值。
3. 数值逼近应用及误差分析:- 结合实际案例,展示数值逼近方法在实际问题中的应用。
- 分析逼近过程中的误差来源,探讨如何降低误差。
hilbert空间的最佳逼近定理的几何意义
hilbert空间的最佳逼近定理的几何意义引言Hilbert空间是数学中一个重要的概念,它是一个完备的内积空间,常常用于描述物理现象和工程问题。
在Hilbert空间中,最佳逼近定理是一个重要的结果,它揭示了在Hilbert空间中,我们可以通过选择合适的元素来最佳逼近一个给定的目标元素。
本文将深入探讨最佳逼近定理的几何意义,以及其在几何学领域中的应用。
Hilbert空间概述在介绍最佳逼近定理之前,先来了解一下Hilbert空间的基本概念。
Hilbert空间是一个实或复的向量空间,配以一个内积,它是一个完备的度量空间。
在Hilbert空间中,我们可以定义向量的长度、角度和距离,这使得Hilbert空间成为了研究几何性质和进行几何分析的理想工具。
最佳逼近定理的表述最佳逼近定理是Hilbert空间理论中的一个重要结果,它描述了如何通过选择合适的元素来最佳逼近一个给定的元素。
具体而言,最佳逼近定理表明,在Hilbert空间中,对于任意一个给定的元素,总存在一个最佳逼近序列,使得在所有逼近中,这个序列收敛到目标元素,并且存在一个收敛的逼近序列能够达到最佳逼近。
最佳逼近定理可以用数学公式表示如下:定理:设H为Hilbert空间,f是H中的一个元素,E是H中的子空间。
则存在一个最佳逼近序列{en},使得:1.对于任意n,en属于E;2.对于任意e属于E,||f-en|| <= ||f-e||,其中||.||表示H中的范数;3.对于序列{en}的每一个子序列{en_k},都存在项ek使得||f-ek||是最小的。
最佳逼近定理的几何意义最佳逼近定理的几何意义十分重要,它从几何的角度解释了Hilbert空间中的最佳逼近现象。
在Hilbert空间中,我们可以将元素看作空间中的点,而子空间可以看作空间中的平面或曲面。
最佳逼近定理告诉我们,在给定一个点的情况下,我们可以选择一个平面或曲面,使得这个点到平面或曲面的距离最小。
计算方法最佳一致逼近多项式切比雪夫多项式
(x) pn(x)|
(3.2)
称为f(x)在 b][上 a,的最小 。偏差
定义8设f(x) C[a,b], 若存在n*p(x) Hn, 使得
Δ(f,pn*) En (最小偏差),(3.3) 则称pn*(x)是f(x)在a,[b]上的n次最佳一 逼致 近多项式,简称最近 佳多 逼项式。
(4)切比雪夫多项式的零点
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(2k2n1)π)]
cos[(2k 1)π] 2
而上式成立的充分必要条件是x0, x1,…xn是切比雪夫 多项式的0点。
将Lagran值g多 e插 项n(式 x)L的节点n取 1(x为 )
的0点 :
x k
cos(2k-1)π, k 2(n1)
1, … ,n,n1
此时 Ln(,x)具有近致 似逼 最近 佳的 一性
定理7设f(x) Cn1[1,1]L,n(x)为插值多项式, 其插值节点 0, x1,...,xn取为切比雪夫多项 n1式T 的0点,则 n(xL)是f(x)-在 1,1[ ]上的最佳一致逼近 多项m - 1 式 x 1 , |f a 且( x L x n () x |2 ) n (1 n 1|)f |(! n 1 ()x | |)
T0(x ) cos(01) T1(x )cos(ax rc)cxos T 2(x)cos(2 sa xr )2 c2 x c 1 o
T 3 (x )cos(3 sa xr )4 c3 x c 3 ox
4最佳逼近
第四章 最佳逼近1解:作变换)()(a b t a t x -+==ϕ,则当],[b a x ∈时,]1,0[∈t ,记:]1,0[)),(())(()(∈-+==t a b t a f t f t g ϕ,则其Bernstain 多项式为:,2,1,)1()()1()(0=-⎪⎭⎫ ⎝⎛-+=-⎪⎭⎫ ⎝⎛=-=-=∑∑n t t C a b n i a f t t C n i g g B in i n i ni i n i n i ni n再将ab a x t --=代入上式即得)(x f 在],[b a 上的Bernstain 多项式:,2,1,)1()()(0=---⎪⎭⎫⎝⎛--⎪⎭⎫ ⎝⎛-+=-=∑n a b a x a b a x C a b n i a f f B i n in i ni n 4.假设],[b a C f ∈,证明f 关于0P 的最佳一致逼近多项式为:2m M +,其中:)(max ),(min ],[],[x f M x f m b a x b a x ∈∈==。
证明:因)(x f 在],[b a 上连续,故存在],[,21b a x x ∈使:)(min )(],[1x f m x f b a x ∈==,)(max )(],[2x f M x f b a x ∈==。
(1) 若M m =,则)(x f 为常数,显然2m M c +=就是)(x f 在],[b a 上的0次最佳一致逼近多项式。
(2) 若M m ≠则21x x ≠,且记2m M c +=,由于当M x f m b x a ≤≤≤≤)(,时,于是:22)(m M mM c x f --≤-≤-,即2)(mM c x f -∞=-,又∞--=-cx f c x f )()(1,∞-=-c x f c x f )()(2,即21,x x 为误差曲线c x f -)(的两个正负相间的偏差点,由契比雪夫定理知,2mM c +=就是)(x f 在],[b a 上的0次最佳一致逼近多项式。
第4章 函数最优逼近法
n
n
n
2 2
由多元函数极值的必要条件,S取得最小必有 S 0, (k 0,1, , n) ck
S ( c j j , crr ) 2( c j j , f ) f
j 0 r 0 j 0
n n S 又 cr (k ,r ) c j (k ,i ) 2(k , f ) ck r 0 j 0
令 = max xi , 于是 0有 S
b a
w( x )[ p ( x ) f ( x )] 2 dx
(4 3)
选取c0 , c1 ,
, cn , 使得S最小的p( x)称为f ( x)在[a, b]上相对于
权函数w( x)的最优平方逼近函数.
4.1.1 最优平方逼近函数之内积
(0 , n ) c0 (0 , y ) (1 , n ) c1 (1 , y ) (n ,n ) cn (n , y )
n S 由 0得, (k ,i )c j (k , f ), k 0,1, ck j 0
n
n
n
2 2
,n
(0 ,0 ) (0 ,1 ) (1 ,0 ) (1 ,1 ) (n ,0 ) (n ,1 )
简写为 Ac b
(0 ,n ) c0 (0 , f ) (1 ,n ) c1 (1 , f ) (n ,n ) cn (n , f )
0
1
所以
0.06277
对离散数据的曲线拟合最小二乘法
问题(回顾):在科学实验中,得到函数y=f(x)的一 组实验数据: ( xi , yi ) (i 1,2, , m) ,求曲线
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第四章最佳逼
近
学习目标:掌握最佳一致逼近和最佳平方逼近的基本理论和
方法、以及最小二乘法常用
的正交多项式以及正交多项
式的性质。
重点为最佳一致
逼近和最佳平方逼近的特征
性质(如契比雪夫定理等)
以及最佳一致逼近和最佳平
方逼近多项式的计算方法。
§1 C[a ,b ]上的最佳一致逼近
不难验证,[a ,b ]上所有连续函数的全体构成一无限维线性空间,
简记为C[a,b]。
为描述方便,引进符号函数 ,称为C[a,b]
上的一致范数或契比雪夫(Chebyshev )范数,其定义为
∞⋅],[]
,[,)(max b a b a x C f x f f ∈∀=∈∞考虑所有n 次代数多项式的全体形成的集合
. 不难验证,P n 是C [a ,b ]上的n+1维线性子空间。
{
}n n x x span P ,,,1 =
对给定的函数f (x )∈C [a ,b ]称量:
)
()(min ),(x p x f P f n P p n -=∆∈为f (x )关于P n 的最佳一致逼近,简称最佳逼近,也称为契比雪夫逼近。
满足上式的多项式p *(x )称为f (x )在[a ,b ]上的最佳逼近多项式,而线性空间 P n 也称为逼近子空间。
围绕这一问题,人们马上会问:最佳逼近多项式是否存在?是否唯一?如果存在,如何寻找或构造它?对这些问题的回答构成了最佳一致逼近研究的中心内容。
定理(契比雪夫定理) 对任意
是f 的最佳一致逼近多项式的充要条件是f - p 在[a ,b ]上存在的至少有n +2个点组成的交错点组。
n b a p p C f ∈∈,],[推论1 如果 ,那么在
中存在唯一的元素为f 的最佳一致逼近多项式
],[b a C f ∈n p 推论 2 如果f 在[a ,b ]上有n +1阶导数,且 在
(a ,b )上保号(恒正或恒负),那么契比雪夫交
错组唯一,且区间[a ,b ]的端点属于契比雪夫交错组。
)1(+n f
§4.3 最佳平方逼近
假设X是一线性空间,如果其上赋予一个满足如下4条性质的二元函数(·,·):
(1)(x, y)=(y, x),x, y∈X;
(2)(αx, y)=α(x, y),x, y∈X;α∈R;
(3)(x+y, z)=(x, z)+(y, z),x, y, z∈X;
(4)(x, y)≥0,x, y∈X, (x, y)=0↔x=0
则称X为内积空间。
最简单的内积空间是欧几里德空间,它是线性空间R n,按内积
(x, y)=x T y, x, y∈R n (20)
所构成,其中T表示向量的转置。
容易验收证,由(20)式定义的内积满足内积的上述4条性质。
利用内积可以引进范数,若X 是一内积空间,则容易
验证函数 满足范数的3条性质。
于是内
积空间按上述范数构成赋范线性空间,由此可以引出内积空间中的最佳逼近问题:假设φi (i =1,2,…,n )是内积空间X 中n 个线性无关的元素,f ∈X ,则子集 Φn=span{φ1,φ2,…,φn}
对f 的最佳平方逼近定义为:
2min ),(ϕϕ-=Φ∆Φ∈f f n n (21)
特别使(21)式成立的那个元素称为最佳逼近元。
),(2x x x =
类似地,可以提出如下问题:最佳平方逼近元是否存在?如果存在,是否唯一?其特征又如何?对此我们首先建立如下的结论:
定理12 假设X 是内积空间,f ∈X ,则φ*∈Φn 为f 的最佳逼近元的充分必要条件是
(22) n i f i ,,2,1,0),(* ==-ϕϕ。