Chap4Sec3插值与数据拟合3.pptVIP

机电工程学院 样条插值与最小二乘拟合 三次样条插值背景 三次样条插值的分类 三次样条插值函数的构造 最小二乘拟合意义 矛盾方程组的最小二乘解 多项式最小二乘拟合方法 三次样条插值背景 分段低阶插值，收敛性好，但光滑性不够理想。在工业设计中，对曲线光滑性要求高，如：流线型机翼、船体设计 设想这样一插值曲线：次数不高于3次，整个曲线2阶连续可导，称为三次样条函数插值。 §4.6 三次样条插值 定义4.6.1 　　　 设 。三次样条函数 , 且在每个 上为三次多项式 。若它同时还满足 ，则称它为 f 的三次样条插值函数。 样条插值与最小二乘拟合 三次样条插值背景 三次样条插值的分类 三次样条插值函数的构造 最小二乘拟合意义 矛盾方程组的最小二乘解 多项式最小二乘拟合方法 三次样条插值的分类 n+1个节点n个子区间 n+1个条件 3n-3个条件 三次样条插值的分类 三次样条插值的分类 由数据特征决定 样条插值与最小二乘拟合 三次样条插值背景 三次样条插值的分类 三次样条插值函数的构造 最小二乘拟合意义 矛盾方程组的最小二乘解 多项式最小二乘拟合方法 三次样条插值函数的构造 通常记为S[i](x) 拉格朗日插值 三次样条插值函数的构造 在区间[xi-1, xi]积分2次， 可得 S(x) 和 S(x) ： 三次样条插值公式 同时得到， 其中x ∈[xi-1, xi] 同时得到， 其中x ∈[xi-1, xi] 其中x ∈[xi-1, xi]且xi-0= xi 其中x ∈[xi, xi+1]且xi+0= xi 插值区间内点 共n-1个方程, n+1个变量 还有什么条件没有用到呢？ 三次样条插值函数的构造 同样补充了两个方程 三次样条逼近f(x)的步骤 ① 计算 ?j , ?j , gj ; ② 计算 Mj (追赶法等) ; ③ 找到 x 所在区间 ( 即找到相应的 i ) ; ④ 由该区间上的 S[i](x) 算出 f(x) 的近似值。 样条插值与最小二乘拟合 三次样条插值背景 三次样条插值的分类 三次样条插值函数的构造 最小二乘拟合意义 矛盾方程组的最小二乘解 多项式最小二乘拟合方法 §4.7 曲线拟合的最小二乘法 仍然是已知 x1 … xm ; y1 … ym, 求一个简单易算的近似函数 P(x) ? f(x)。 但是 ① m 很大； ② yi 本身是测量值，不准确，即 yi ? f (xi)。 这时没必要取 P(xi) = yi , 而要使 P(xi) ? yi 总体上尽可能小。 常见做法： ? 使 最小 /* 最大最小问题 */ 太复杂? ? 使 最小 不可导，求解困难? ? 使 最小 /* 最小二乘法 */ 　　考虑一般的线性无关函数族?={?0(x), ?1(x), … , ?n(x), … }，其有限项的线性组合 称为广义多项式。 常见多项式： ? { ?j(x) = x j } 对应代数多项式。 ? { ?j(x) = cos jx }、{ ?j(x) = sin jx } ? { ?j(x), ?j(x) }对应三角多项式。 ? { ?j(x) = e kj x , ki ? kj } 对应指数多项式。 对于一组数据 (xi, yi) (i = 1, 2, …, n) 使得 达到极小，这里 m n。 设： 样条插值与最小二乘拟合 三次样条插值背景 三次样条插值的分类 三次样条插值函数的构造 最小二乘拟合意义 矛盾方程组的最小二乘解 多项式最小二乘拟合方法 使Q取得极小值的解x称为矛盾方程组Ax=b的最优近似解。 该方程称为矛盾方程组Ax=b的法方程或正规方程。法方程的解就是矛盾方程组的最优近似解。 样条插值与最小二乘拟合 三次样条插值背景 三次样条插值的分类 三次样条插值函数的构造 最小二乘拟合意义 矛盾方程组的最小二乘解 多项式最小二乘拟合方法 本次课结束！ 谢谢大家