数值计算-方程.pptVIP

塞德尔采用迭代公式 简单迭代公式 k 0 1 2 3 4 5 6 x1(k) 0 0.3000 0.8804 0.9840 0.9978 0.9997 1.0000 x2(k) 0 1.5600 1.9445 1.9922 1.9987 1.9999 2.0000 x3(k) 0 2.6840 2.9539 2.9938 2.9991 2.9999 3.0000 迭代法求解线性方程组一般步骤 建立迭代格式 设置初始值X(0),精度要求ε max(|x(k+1)-x(k)|) ε? 迭代计算,k=k+1 迭代结束，X*=X(k+1) 迭代格式的收敛性 x1 x2 x3 0 0 0 0 1 11 -14 -3 2 -69 81 66 3 -499 -374 -429 不可约： 对角优势： 若矩阵A不能通过行的次序调换和相应的列的次序调换成为： 则A不可约 若矩阵A满足 且自少有一个i值，上式中严格的不等号成立，则A具有对角优势 若系数矩阵A不可约且具有对角优势，简单迭代法必收敛。 定理： 雅克比 赛德尔 初值 迭代14次 雅克比 迭代11次 塞德尔 矩阵求逆 设Ti 为初等行变换，若矩阵A经一系列初等行变换后变为单位矩阵I，即： 则有： 提供了一种矩阵求逆的方法 A-1 高斯消去法 A-1 例：求A的逆阵A-1 作业 建立一个收敛的迭代格式，来求解下列方程组，并编出相应的计算程序(FORTRAN语言或C语言均可，提交纸质文档)： * Ax=b能用高斯消去法的充要条件是A的各阶顺序主子式不为零 * A的行列式不为零，称为非奇异阵 线性方程组数值解法及矩阵求逆 火箭发射速度和加速度 测定过程 Time, t Velocity, v s m/s 5 106.8 8 177.2 12 279.2 方程为： 矩阵形式： T DNS P ER CN QF SF 超重元素合成： A1=AP N阶线性方程组： 矩阵表示： 数值求解方法： 2. 迭代解法 a) 简单迭代法(雅克比) b) 赛德尔迭代法 1. 直接解法 a) 高斯消元法 b) 主元素消元法 高斯消去法 a1x1 + a2x2 + a3x3 = d b1x1 + b2x2 +b3x3 = e c1x1 + c2x2 + c3x3 = f X1 = … X2 = … X3=… 增广矩阵 矩阵行运算 矩阵初等行变换 矩阵初等行变换 互换矩阵两行位置 用非零数乘(除)矩阵某行 将矩阵某行的倍数加到矩阵的另一行上 注：只要 A 非奇异，即 A?1 存在，则可通过逐次消元及初等行变换，将方程组化为三角形方程组，求出唯一解。 第一步：首先用 除第一行所有元素，然后对下面(n–1)行中的第i行分别用 乘第一行再与第i行相加，以消去第一列中的其余元素，此时变为： 其中： 其中： 第k步，是对矩阵 进行下列运算： r1/10 r2+2r1 r3+r1 r2/9.6 r3+2.2r2 r3/4.625 如果 ，消去过程会失败 如果 ，会使计算精度降低 解决方法：主元素消去法 例： 精确解为: 8个 8个 用高斯消去法计算： 方法： 列主元消去法 行主元消去法 全主元消去法 思路：对调方程的次序或变量的排列, 使得除数最大 例： 解为: 8个 8个 例： 回代： 即： 得： 简单迭代法(雅克比) . . . . . . 算法： 即： old new 当相对误差小于某预先设定值时，停止迭代 k 0 1 2 3 4 5 6 7 8 x1(k) 0 0.3000 0.8000 0.9180 0.9716 0.9804 0.9962 0.9986 0.9995 x2(k) 0 1.5000 1.7600 1.9260 1.9700 1.9897 1.9961 1.9986 1.9995 x3(k) 0 2.0000 2.6600 2.8640 2.9540 2.9823 2.9938 2.9977 2.9992 则可得解X=(0.9995, 1.9995, 2.9992)T 塞德尔迭代法 例： 计算物理·数值计算IICOMPUTATIONAL PHYSICS·NUMERICAL COMPUTATION EQUATION SOLVIN