(压缩)8计算方法第八章函数逼近与曲线拟合(2011-2012学年).pdf

第八章 函数逼近与曲线拟合 8.1 引言 8.2 正交多项式 8.3 最佳平方逼近 8.4 曲线拟合的最小二乘法 8.5 最佳平方三角逼近与离散傅立叶变换 计算方法 8.1 引言 在科学计算与工程应用中，经常要计算函数值和描绘函数曲 线．当函数仅在某区间中的有限点集上给定函数值时，要在该区 间上给出近似于该函数的简单表达式，或者当函数表达式相对复 杂时，要得到与它相近的简单函数，这些构造区间[a, b] 上的简单 函数去近似已知函数的问题，就是函数逼近问题．上章讨论的插 值法是函数逼近问题的一种，插值函数在结点处的值等于已知函 数值，而逼近函数则可以不然． 本章讨论的函数逼近问题是指：对函数类 A 中给定的函数 f (x) ，即f (x) ∈A ，要得到在另一类简单的便于计算的函数类 B 中的函数p (x) ，即p (x) ∈B ，使p (x) 与f (x) 的误差在某种 度量意义下是最小的． 计算方法 函数类A 通常是无限维的，例如区间[a, b] 上的全体连续 函数．函数类B 通常是有限维的，例如次数不超过 n （n 为 正整数）的全体多项式．函数类B 还可以是某个有理函数类、 分段低次多项式类等． 一般地，设有一组在 C[a, b] 中的线性无关的函数集合 n { } ϕ (x) ，则元素 i i 0 ϕ(x) ∈Φ span{ϕ (x), ϕ (x), ,ϕ (x)} ⊂C[a, b] 0 1 n 可表示为 ϕ(x) a ϕ (x) +a ϕ (x) ++a ϕ (x) 0 0 1 1 n n 常见的函数逼近问题就是：对任给f (x) ∈C[a, b] ，在子 * * 空间Φ中找一个元素ϕ (x) ∈Φ，使f (x) −ϕ (x) 在某种度 量标准下最小．最常用的度量标准有两种： 计算方法 1）以函数f (x) 和ϕ(x) 的最大误差 max f (x) −ϕ(x) a≤x ≤b 作为 f (x) 和ϕ(x) 的误差大小的度量标准，在这种意义下的 函数逼近称为一致逼近． 2 ）用积分 b 2 ∫ [f (x) −ϕ(x) ]dx a 作为 f (x) 和ϕ(x) 的误差大小的度量标准，在这种意义下的 函数逼近称为平方逼近． 关于一致逼近的理论，本课程不做详细讨论，只叙述著 名的魏尔斯特拉斯（Weierstrass ）定理如下： 计算方法 ε 定理 8.1 若f (x) ∈C[a, b] ，则对任何 0 ，总 存在多项式p (x) ，使得对于一切a ≤x ≤b 都有 |f (x) −p ( x) |ε ． 本章主要研究多项式 p (x) 与区间[a, b] 上的连续函数 f (x) 的最佳平方逼近问题，还要介绍曲线拟合，讨论最佳平 方三角逼近和快速傅立叶变换，最后简述有理逼近．下面介 绍函数逼近中常用的几种正交多项式． 计算方法 8.2 正交多项式 8.2.1 正交函数系 正交多项式是函数逼近的重要工具，在下一章将讨论的 数值积分中也有重要应用．下面先给出权函数、带权内积和