数值计算六讲.pptx

数值分析与算法 (6) Numerical Analysis Algorithms 清华大学 计算机系 喻文健 Wenjian Yu 2 第六章 函数逼近与函数插值 用较简单的函数近似表示未知函数、或已知的复杂函数 逼近：整体上近似 (整体误差最小) 插值：离散自变量点上函数值相等 (误差为0) 本章内容 函数逼近的基本概念 连续函数的最佳平方逼近 曲线拟合的最小二乘法 多项式插值 (拉格朗日, 牛顿) 分段多项式插值 样条函数插值 Wenjian Yu 3 函数逼近问题的例子 傅里叶变换, 信号的频谱分析 表格函数, 曲线拟合 Wenjian Yu 4 函数逼近的基本概念 Wenjian Yu 5 函数逼近的基本概念 Wenjian Yu 6 函数逼近的基本概念 A 复内积则不同 Wenjian Yu 7 函数逼近的基本概念 充分性的证明 … 矛盾! Wenjian Yu 8 函数逼近的基本概念 (非负) (多项式可积) (无局部恒为零) (广义的2-范数) 不能做权函数 Wenjian Yu 9 函数逼近的基本概念 Wenjian Yu 10 连续函数的最佳平方逼近 Wenjian Yu 11 最佳平方逼近 n元二次函数 法方程方法 Wenjian Yu 12 最佳平方逼近 充分性的证明看书 Wenjian Yu 13 最佳平方逼近 Hilbert矩阵 6.2.2小节 Wenjian Yu 14 正交函数族与正交多项式 正交 适合于任意线性空间基的正交化 Wenjian Yu 15 几种正交多项式 法方程的求解变得很简单, 直接计算 其他定义域的正交多项式 有关习题 Wenjian Yu 16 曲线拟合的最小二乘法 Wenjian Yu 17 曲线拟合问题 最佳平方逼近问题! 实例演示Matlab, Excel Wenjian Yu 18 线性最小二乘 Wenjian Yu 19 线性最小二乘 – 法方程法 算法6.2 Wenjian Yu 20 线性最小二乘 – 法方程法 (用多项式拟合比较放心) 用加权内积技术求解① ① ② Wenjian Yu 21 线性最小二乘 – 法方程法 例6.5, 6.6 ① ② Wenjian Yu 22 线性最小二乘 – 正交变换法 例6.7 正交阵 1965, G. Golub提出 Wenjian Yu 23 线性最小二乘 – 正交变换法 mn mm mn n列 n行 取等号条件: 数值稳定, 实用 Wenjian Yu 24 线性最小二乘 Wenjian Yu 25 多项式插值 Wenjian Yu 26 插值与多项式插值 插值的基本概念 为离散点配曲线，并要求曲线通过各个离散点 一种特殊的“逼近” 目的与用途 MS Word的曲线绘制功能 图形学/CAD: 画一条通过离散点的光滑曲线 对表达式未知的表格函数, 估算中间点函数值 快速、方便地计算复杂数学函数的函数值 其他数值方法的基础: 简单函数近似复杂的或未知函数(非线性方程、数值积分与微分、微分方程数值解法) ？ Wenjian Yu 27 插值与多项式插值 定义6.6 DFT Wenjian Yu 28 插值与多项式插值 解存在、唯一吗? (Vandermonde阵) Th6.6 Wenjian Yu 29 Lagrange插值法 “两点式”直线公式: 不便于计算、以及理论分析 Wenjian Yu 30 Lagrange插值法 (一个常用的技巧) Lagrange插值函数 Wenjian Yu 31 多项式插值误差估计 Rolle定理 至少有n+1个互不相同的 Wenjian Yu 32 Lagrange插值小结 (例6.9) Wenjian Yu 33 Newton插值 插值多项式 “点斜式”直线公式: n个插值节点: 牛顿插值公式: 怎么算? Wenjian Yu 34 Newton插值 (为了方便地计算Newton插值系数) (递归定义) Th6.8 Th6.9 Wenjian Yu 35 Newton插值 有k-1个节点相同 (例6.10) Wenjian Yu 36 Newton插值 (多项式插值余项的另一种形式) 插值余项 Wenjian Yu 37 Newton插值 计算复杂, n较大时矩阵有病态性 便于动态增、减插值节点 思考题: 比较两种插值法计算未知点处函数值的计算量 Wenjian Yu 38 分段多项式插值 Wenjian Yu 39 插值与多项式插值 插值的基本概念 为离散点配曲线，并要求曲线通过各个离散点 一种特殊的“逼近” 目的与用途 MS Word的曲线绘制功能 图形学/CAD: 画一条通过离散点的光滑曲线 对表达式未知的表格函数, 估算中间点函数值 快速、方便地