矩阵论 复习典型例题.pptVIP

计算方法复习 典型概念例题 零 绪论 误 差 及 算 法 误差 算法 分类 度量 传播 舍入 截断 绝对 相对 有效数字 一元函数 n 元函数 一 插值与逼近 插值法 工具 多项式插值 分段多项式 插值 差商 差分 插值基函数 存在唯一性 误差估计 插值公式 Hermite 插值 分段线性 分段三次 Hermite 插值 三次样条插值 函数逼近 预备知识 函数逼近方法 范数 内积 正交多项式 最佳一致逼近 最佳平方逼近 最小二乘拟合 三角函数逼近 帕德逼近 所以，关于 a ０ ， a １ 为未知数的法方程组为 求 g ( x )= ?? x 在 P 1 [0,1] 中的最佳平方逼近元 解法一 这是 C[0,1] 上的最佳平方逼近问题 . 取 ? ０ ＝１ , ? １ ＝ x , P 1 [0,1] ＝ span {1, x } 记 p 1 ( x )= a ０ ＋ a １ x ( ? ０ , ? 0 )=1 ， ( ? ０ , ? 1 )=1/2 ， ( ? 1 , ? 1 )=1/3 ( ? ０ ,g)=2/3 ， ( ? 1 ,g)=2/5. 例 1 解得 a ０ =4/15 ， a １ ＝ 4/5 为 P 1 [0,1] 中对 g ( x )= ?? x 的最佳平方逼近元 . ? ? ? ? ? ? ? ? ? 5 2 3 1 2 1 3 2 2 1 1 0 1 0 a a a a 即 p 1 ( x )=4/5 x +4/15 例 1 观测物体过原点的直线运动 , 得到所示数据 , 求运动方程 . 时间 t/s 0 0.9 1.9 3.0 3.9 5.0 距离 s/m 0 10 30 50 80 110 解 作直线模型 : at + s =0 n 为观测点数 定义残差向量 : T 1 1 2 2 ( , , , ) n n V at s at s at s ? ? ? ? 2 2 ( ) I a V ? 2 1 ( ) n i i i at s ? ? ? ? 2 1 1 2 2 n n i i i i i a t t s ? ? ? ? ? ? 1 2 ( ) n i i i i dI at s t da ? ? ? ? 所以 : 2 53.63 2 1078 dI a da ? ? ? ? 令 : 2 53.63 2 1078 0 dI a da ? ? ? ? ? 20.1007 a ? ? 所求运动方程为 : 20.1007 0 t s ? ? ? 二 数值积分 数 值 积 分 基本概念 Gauss 求积公式 代数精度 插值型求积公式 收敛及稳定性 数值求积思想 N-C 公式 Romberg 求积公式及外推加速 梯形公式 辛普森公式 例 2 试确定常数 A,B,C 及 α , 使求积公式 : 解 2 2 ( ) ( ) (0) ( ) f x dx Af Bf Cf ? ? ? ? ? ? ? ? 代数精确度尽可能高，并确定上述公式的代数精 确度。是否为高斯型求积公式 . ? ? 1 ? x f 令 : 2 2 4 dx A B C ? ? ? ? ? ? ? ? x x f ? 2 2 0 xdx A C ? ? ? ? ? ? ? ? ? ? 2 x x f ? 2 2 2 2 2 16 3 x dx A C ? ? ? ? ? ? ? ? ? 3 f x x ? 2 3 3 3 2 0 x dx A C ? ? ? ? ? ? ? ? 整理得 : A C ? 2 4 A B ? ? 2 8 3 A ? ? ? ? 4 x x f ? 2 4 4 4 2 64 5 x dx A C ? ? ? ? ? ? ? 4 32 5 A ? ? 2 8 32 3 5 ? ? 12 5 ? ? ? 10 9 A ? 16 9 B ? 10 9 C ? ? ? 5 f x x ? 2 5 5 5 2 0 x dx A C ? ? ? ? ? ? ? ? ? ? 6 f x x ? 2 2 6 6 6 2 2 128 8 12 2 7 3 5 x dx A C ? ? ? ? ? ? ? ? ? ? 所以代数精确