计算方法7.pptVIP

2、第k次消元 (2≤k≤n-1) ： 设第k-1次消元已完成，有 用-mik乘第k个方程，加到第i个方程，消去第i个 方程到第n个方程的未知数xk，得到 其中 3、第n-1次消元后，有 消元过程自上而下进行，每次都是顺序选取主 对角线上元素akk(k)作为主元素，消去以下各行 相应的变元，叫顺序消去法。 二、回代过程 对上式自下而上逐步回代，得 其中： 主元素所在的行称为主行。 三、高斯消去法计算步骤总结 用k表示消元过程的次序时， 1、消元过程 2、回代过程 四、计算机计算高斯消去法框图及特点 开始 读入aij, bi,i,j=1,2,…,n 1 k aik/akk = aik i=k+1,k+2,…,n k=n-1? 输出b1,b2,…,bn 结束 k+1=k y N 选主元 aij-aikakj = aij i,j=k+1,k+2,…,n bi-aikbk = bi i=k+1,k+2,…,n 计算机计算高斯消去法特点 1、按消元规则进行运算后，对角线以下元素 为0，故运算中对角线以下元素不作计算， 以减小计算量。 2、对角线以下元素对回代求解无影响，故将 乘数放在该处，即 以节省存储单元。 3、对角线以上元素和常数项采用“原地”工作 方式，即经变换后的元素仍放在原来的位 置上以节省存储单元。 aij-aikakj = aij i,j=k+1,k+2,…,n bi-aikbk = bi i=k+1,k+2,…,n 4、回代后的值仍放在常数项存储单元以节省存储单元。 这时b1,b2,…,bn中存放的就是输出值x1,x2,… ,xn 。 例1：用计算机求解方程组 解： 用增广矩阵各元素的位置代表计算机的存储单 元，增广矩阵元素数值的变化代表存储单元所 存放数据的变化。 放入初始数据 第一次消元 第二次消元 第一次回代 第二次回代 第三次回代 * 2.2.4 迭代过程的加速 一、加权法 二、埃特金加速法和斯蒂芬森迭代法 1、埃特金加速法 2、斯蒂芬森算法 ② 说明： 埃特金算法适用于任意线性序列，斯蒂芬森算法只针对不动点迭代。 只要迭代函数φ(x)满足φ(x*)≠1。不管原迭 代法是否收敛，由它构造的斯蒂芬森迭代法至 少是平方收敛的。它是对不动点迭代法的收敛 性和收敛速度的一种改善方法。 当原来的迭代法已有二阶或更高收敛阶时， 没必要使用斯蒂芬森算法，这时加速效果已 不显著。 2.3.1 牛顿迭代公式 一、迭代公式 定理1: 2.3.3 改进牛顿迭代公式 一、简化牛顿法 其中c为常数 二、牛顿下山法 下山成功条件： 三、重根时的修正算法 则迭代过程 是线性收敛，不是平方收敛。 修正算法 是平方收敛的。 这种修正方法需要预知重根数m的值。 例1: 解: 代入牛顿迭代法，有 所以牛顿迭代法为线性收敛。 m=2。若选取迭代格式 所以牛顿迭代法为平方收敛。 例2： 解： 牛顿迭代法 重根数m=2时，修正公式 取x0=1.5 1.414216865 1.436607143 2 1.414213562 1.425497619 3 1.458333333 1.5 牛顿值xk 1.416666667 1 1.5 0 xk (m已知时) k 2、当重根数未知时，令 牛顿迭代法修正为 则该方法是二阶收敛的。 迭代格式 取x0=1.5 1.414213562 1.414216865 1.416666667 1.5 xk (m已知时) 1.414211438 1.436607143 2 1.414213562 1.425497619 3 1.458333333 1.5 牛顿值xk 1.411764706 1 1.5 0 xk (m未知时) k 例： 解： 牛顿迭代法 作业：《数值计算方法》P55：2.1-2.14题 《科学与工程计算》P67：1-17题 第三章 解线性代数方程 组的直接法 3.1 高斯消去法 3.1.1 基本思想 3.1.2 计算步骤 3.1.3 使用条件 3.1.4 计算量 3.2 列主元高斯消去法 3.2.1 数值不稳定 3.2.2 选主技术 3.2.3 计算矩阵行列式 3.2.4 高斯—约当消去法 3.3 矩阵三角分解法 3.3.1 矩阵三角分解原理 3.3.2 解线性方程组的三角分解法 3.3.3 平方根法 3.3.4 追赶法 3.4 向量和矩阵的范数 3.4.1 向量的范数 3.4.2