数值分析5_6最佳平方逼近.ppt

函 数 逼 近;f(x)=ex 在[-1,1]上用：; 它在整个区间上误差较大，若在计算机上用这种方法计算ex ，如精度要求较高，则需取很多项，这样即费时又多占存储单元。因此，我们要求在给定精度下计算次数最少的近似公式，这就是函数逼近要解决的问题。;?定义 近似代替又称为逼近，函数 f(x) 称为被逼近函数；P(x) 称为逼近函数，两者之差称为逼近的误差。;?最常用的度量标准有两种： 1、一致逼近（均匀逼近） 以 作为度量误差f(x)- P(x) 的“大小” 标准。 2、平方逼近（均方逼近） 以 作为度量误差f(x)- P(x)的“大小” 标准。;5.6 函数的最佳平方逼近;定理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 . 利用内积的性质, 可得;(?) 若 ( f ?p*, ?j )=0, j=0,1,…,n 成立, 对任意的 p(x)?Hn ,有;定理5.8 设 f (x)?C[a, b], 在子空间 Hn 中, 对 f (x) 最佳平方逼近的函数是唯一的. 证明 假定, 在Hn 中, p(x) 和 q(x) 都是对 f (x) 最佳平方逼近的函数, 由定理5.7的系,知;三、最佳平方逼近函数的求解 利用 ( f ?p*, ?j )=0, 可求出最佳平方逼近函数 p*. 设 故; 由于 ?0?x?,?1?x?, ?, ?n?x? 线性无关, 可以推得上系数阵是非奇异的. 故 (5. 82) 有唯一解 { c*j }. ;例6 定义内积 试在H1=Span{1,x}中寻求对于f(x)= 的最佳平方逼近元素p(x).; 当 ??0?x?,?1?x?, ?, ?n?x?, ? ?是正交系时,求解最佳平方逼近式(5.82)中的系数非常容易.;;一、Legendre 多项式的应用;遇到区间[a,b], 通过下面的变换把问题转化到[-1,1]上处理.;;所求三次最佳平方逼近多项式为;例8 求 f(x)= 在区间 [0,1] 上的一次最佳平方逼近多项式.;把 t=2x-1 代入 q1(x)，就得 在区间[0,1]的一次最佳平方逼近多项式: ;二、Chebyshev多项式的应用;函数 f 的 C-最佳平方逼近函数为; 例9 求 f(x)=arcsinx 按 Chebyshev 多项式展开的 n=7 的部分和.;;数学符号 x0,x1, ?, xm, ??0?x?,?1??x?, ?, ?n?x?, ? ? k=0,1, ?, n ? ??x??0 n ?0 a ? x ? b k ? 1, n?m,,, ??? ?? ? ?? ? ? ?????? ??? ??????? ΓΔΘΛΞΦΨΩ ?k ?0 f (x)?C[a, b] (xj,yj),j=0,1,…,m, a ? x0 ? x1 ? ? ? xm ?b ,, ????? ????? ????? ????? ????? ????? ????? ??????? ? ????? ??? ? ? ? ? ? ? ? ? ????? ? ? ? ? ? ?? ? ? ? ? ? ? ?? ?? ? ?( f ? p*, f ? p*), , ????????? ????? ????? ?????? ????????? ? ? ? ??????? ? ? ? ? ? ? ? ? ? ? ? ? ①②③④⑤⑥⑦⑧⑨ ?(x) ?1, , ?L0?x?, L1?x?, ?, Ln?x??,