研究生数值分析(20)函数的最佳平方逼近.pptVIP

* * §6 函数的最佳平方逼近 一、最佳平方逼近的概念与解法 用简单函数 p (x)去近似一个给定区间[a，b]上的连续函数 f (x)，是函数逼近要研究的问题。度量逼近误差标准有许多种，这里介绍一种称为平方逼近的函数逼近。 设函数组 都是[a，b] 上连续的函数，并且在 [a，b] 上线性无关。 则 Hn中的任一个元素为 对空间C[a , b] 中的任意两个函数 f 和 g，定义内积 以此为基底，生成空间 C[a，b] 的一个子空间 其中ρ(x)是（a，b）上的一个权函数。 定义 对于给定的函数 f (x)∈ C [a，b] ， 若 满足 为子空间Hn中对于 f (x) 的最佳 则称 平方逼近元素。 特别地，如果 则满足条件(I)的 可称为函数 f (x) 在区间[a ,b] 上带权ρ(x)的n次最佳平方逼近 多项式。具体问题中ρ(x)是给定的, 如果不指明,则ρ(x)≡1，这时内积定义为 定理1 设 f (x)∈ C [a，b] ， 是子空间 H n中对于f (x) 的最佳平方逼近元素的充分必要条件是 或对任一个 ，总有 定理2 设 f (x)∈ C [a，b] ，则在子空间 Hn 中对于 f(x) 的最佳平方逼近元素是唯一的。 因 就是求它所含的系数 求最佳平方逼近元素 故由条件（II），可表达为 这是一个以 为未知数的 n+1 元线性方程组，称该式为法方程组和正规方程。 它的系数矩阵为 由于 在[a，b] 上连续且线性无关。因而矩阵G非奇异，故法方程的解 存在且唯一。 记 称δ为最佳平方 逼近误差， 称为均方误差。 由于 ，所以有 例1 定义内积 试在 中寻求对于 的最佳平方逼近元素 p (x)。 解: 法方程为 解得 所求的最佳平方逼近元素为 二、正交函数系在最佳平方逼近中的应用 对于一般的基底 以及求解法方程的计算量都是很大的。 当 n 较大时，计算法方程中的 作基底，当ρ(x)≡1时， 的法方程系数矩阵G在 n稍大时是病态矩阵，在计算机上求解法方程，其结果不太可靠。为避免这些弊端，应采用正交基底。 虽然 容易计算，但由此形成 设 是区间[a，b]上带权 ρ(x)的正交函数系，即 那么，法方程的解为 Legendre多项式的应用 对给定的函数 f (x)∈ C [a，b] ，在区间[a，b]上求f (x)的n次最佳平方逼近多项式 这个问题，相当于在内积为 的情形下，在子空间 中寻求对f(x)的最佳平方逼近元素 现对该 Hn另取一组基底，即 其中 Lj (x) 是 j 次Legendre多项式。此时，法方程的解可直接得到，就是 所求的n次最佳平方逼近多项式为 如果所给的区间不是[-1，1]，而是一般的有限区间[a，b]，那么，可以通过变量置换将它转化为区间-1≤t≤1来处理，即令 例2 求 在[-1，1]上的三次最佳 平方逼近多项式。 解 采用Legendre多项式 作次数不高于3的多项式空间的基底。由 由公式（III）得 所求最佳平方逼近多项式为 例3 设 在区间[0，1]上的一次最佳 平方逼近多项式。 解 令 得 先求 在区间[-1，1]上的一次最佳平方 逼近多项式 ，由 可知 把 t =2x-1代人 得 在区间[0，1]上的一次最佳平方 逼近多项式 例4 测得铜导线在温度 t j 时的电阻 r j 数据如表所示，求出电阻 r 与温度 t 的近似表达式。 三、最小二乘原理 先考查一个例子。 85.1 83.9 82.35 80.8 79.25 77.8 76.3 50.0 45.1 40.0 36.0 30.1 25.0 19.1 7 6 5 4 3 2 1 把这7个点画在图上 可以看出它近似地在一条直线上 设此直线方程为 （a，b待定） 从图上可知， 不是严格在一条直线上， 上，也就是误差 一般不都全为零。我们希望选择a、b ，使 R j的平方和尽可能地小。 因此，不论怎样选择a，b总是不能使所有的 点都落在直线 即求 ，使 取最小值。 二乘原理。求得的 按最小二乘原理选择近似函数的方法称为最小二乘法。 的原理，称为最小 用这种方法求得 称为拟合函数 或者称为经验公式。 对于给定数据 ，要求在 ，使 其中 为函数类Φ中任意函数。 1、最小二乘法的提法 用最小二乘法求近似函数的问题可以归结为： 因此，用最小二乘法解决实际问题包含如下2个基本环节： 某个函数类Φ中寻求一个函数 （1）确定函数类Φ ，即确定 这不是一个单纯的数学问题，还与其它领域的一些专门知识有关。在数学上，通常的作法是将数据点