计算方法（潘林强）NA02d.pptVIP

第二章8节 曲线拟合与函数逼近 /* Approximation Theory */ 仍然是已知 x1 … xm ; y1 … ym, 求一个简单易算的近似函数 P(x) ? f(x)。 但是 ① m 很大； ② yi 本身是测量值，不准确，即 yi ? f (xi) 这时没必要取 P(xi) = yi , 而要使 P(xi) ? yi 总体上尽可能小。 常见做法： ? 使 最小 /* minimax problem */ 太复杂? ? 使 最小 不可导，求解困难? ? 使 最小 /* Least-Squares method */ §1 最小二乘拟合多项式 /* L-S approximating polynomials */ 确定多项式 ，对于一组数据(xi, yi) (i = 1, 2, …, m) 使得 达到极小，这里 n m。 n a a a 1 0 ? 实际上是 a0, a1, …, an 的多元函数，即 [ ] ? = - + + + = m i i n i n i n y x a x a a a a a 1 2 1 0 1 0 ... ) , ... , , ( j 在 ? 的极值点应有 k i m i n j i j i j x y x a ? ? = = - = 1 0 ] [ 2 - = ? ? ? = = = + n j m i k i i m i k j i j x y x a 0 1 1 2 记 ? ? = = = = m i k i i k m i k i k x y c x b 1 1 , 法方程组(或正规方程组) /* normal equations */ 回归系数 /* regression coefficients */ §1 L-S Approximating Polynomials 定理 L-S 拟合多项式存在唯一 (n m)。 证明：记法方程组为 Ba = c . 对任意 ，必有 。 若不然，则 存在一个 使得 … 即 是 n 阶多项式 的根 则 ? B为正定阵，则非奇异，所以法方程组存在唯一解。 Wait a second! You only gave me a critical point, but it’s not necessarily a minimum point ! 则有 , 其中 §1 L-S Approximating Polynomials 定理 Ba = c 的解确是 ? 的极小点。即：设 a 为解，则任意 b = (b0 b1 … bn )T 对应的多项式 必有 ? = = n j j j x b x F 0 ) ( ? ? = = = - ? - = m i m i i i i i b y x F y x P a 1 1 2 2 ) ( ] ) ( [ ] ) ( [ ) ( j j 证明： ? ? = = - - - = - m i i i m i i i y x P y x F a b 1 2 1 2 ] ) ( [ ] ) ( [ ) ( ) ( j j ? ? = = - - - + - = m i i i m i i i i i y x P y x P x P x F 1 2 1 2 ] ) ( [ ] ) ( ) ( ) ( [ ? ? = = - - + - = m i i i i i m i i i y x P x P x F x P x F 1 1 2 ] ) ( )][ ( ) ( [ 2 )] ( ) ( [ 0 注：? L-S method 首先要求设定 P(x) 的形式。若设n=m?1，则可取 P(x) 为过 m 个点的m?1阶插值多项式，这时? = 0。 ? P(x) 不一定是多项式，通常根据经验确定。 §1 L-S Approximating Polynomials 例： x y (xi , yi) , i = 1, 2, …, m 方案一：设 b ax x x P y + = ? ) ( 求 a 和 b 使得