20160104--数值复习(典型例题)题稿.pptx

计算方法复习 典型概念例题 零 绪论 误差及算法 误差 算法 分类 度量 传播 舍入 截断 绝对 相对 有效数字 一元函数 n元函数 一 插值与逼近 所以，关于a０，a１为未知数的法方程组为 求g(x)=??x 在P1[0,1]中的最佳平方逼近元 解法一 这是C[0,1]上的最佳平方逼近问题. 取?０＝１, ?１＝x, P1[0,1]＝span{1,x} 记 p1(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 为P1[0,1]中对g(x)= ??x的最佳平方逼近元. 即p1(x)=4/5x+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为观测点数 定义残差向量: 所以: 令: 所求运动方程为: 二 数值积分 数值积分 基本概念 Gauss求积公式 代数精度 插值型求积公式 收敛及稳定性 数值求积思想 N-C公式 Romberg求积公式及外推加速 梯形公式 辛普森公式 例2 试确定常数A,B,C及α,使求积公式: 解 代数精确度尽可能高，并确定上述公式的代数精确度。是否为高斯型求积公式. 令: 整理得: 所以代数精确度为5次. 因为代数精确度为2×3=5次,是高斯型求积公式. 三 线性方程组 直接法 选主元Gauss消去法 矩阵三角分解法 向量和矩阵范数 追赶法 矩阵条件数 三 线性方程组 迭代法 基本概念 雅可比迭代 迭代收敛速度 高斯-塞德尔迭代 迭代格式 收敛条件 SOR迭代 例3 解 设线性方程组 的系数矩阵为: (1)写出Jacobi 迭代法的迭代格式 (2)确定a的取值范围，使方程组对应的Gauss-Seidel迭代收敛。 (1) 线性方程组 Jacobi 迭代 (2) 线性方程组 Gauss-Seidel迭代矩阵: 令 得 四 非线性方程求根 求根法 二分法 不动点迭代法及收敛性理论 牛顿迭代法 插值型迭代 弦截法 抛物线法 考虑方程 x = g(x), g(x)?C[a, b], 若 ( I ) 当 x?[a, b] 时， g(x)?[a, b]； ( II )在[a,b]上成立不等式: |g(x1)-g(x2)|≤?|x1-x2| 则（1）g在[a,b]上存在惟一不动点x* （2）任取 x0?[a, b], 由 xk+1 = g(xk) 得到的序列 {xk}（?[a,b]） 收敛于x* （3）k次迭代所得到的近似不动点xk与精确不动 点x*有有误差估计式： 定理1 局部收敛性定理 设x*为g的不动点, g(x)与g’(x)在包含x*的某 邻域U(x*) (即开区间)内连续, 且|g’(x*)|1, 则存在 ?0,当x0∈[x*-? , x*+ ? ]时, 迭代法产生的序列 {xk}?[x*-? , x*+ ? ]且收敛于x*. 定理2 用一般迭代法求方程x-lnx＝2在区间(2, ?)内的根， 要求|xk-xk-1|/|xk|=10-8 令f(x)=x-lnx-2 f(2)0, f(4)0,故方程在(2,4)内至少有一个根 因此f(x)=0在(2, ?)内仅有一个根x* 将方程化为等价方程：x＝2＋lnx 例5 解 因此，? x0?(2, ?), xk+1＝2＋lnxk产生的序列? xk ? 收敛于x* 取初值x0＝3.0,计算结果如下： k xi 0 3.000000000 1 3.098612289 2 3.130954362 3 3.141337866 4 3.144648781 5 3.145702209 6 3.146037143 7 3.146143611 8 3.146177452 9 3.146188209 10 3.146191628 11 3.146192714 12 3.146193060 13 3.146193169 14 3.146193204 另一种迭代格式 0 3.000000000 1 3.147918433 2 3.146193441 3 3.146193221 五 常微分方程数值解 数值解法 单步法(Euler法, 改进Euler法) 线性多步法 方程组与高阶方程 重要概念 重要构造方法 局部截断误差 方法精度 差分构造 泰勒展式构造 积分构造 例5 解 给定求解常微