线性以及非线性方程组解法.pptVIP

3.4 非线性方程组的解法 Gauss-Seidel迭代； 3.4 非线性方程组的解法 加速系数法； 迭代 加速 合适的w取值，加速收敛，避免数值振荡； 3.4 非线性方程组的解法 牛顿迭代法； 非线性方程 f(x)=0； 构造序列 uk 逼近根 x*； 迭代式； x* xk xk-1 xk+1 3.4 非线性方程组的解法 牛顿法的收敛性； ek= xk - x*； 如果 ，那么迭代法是 p 阶收敛的； 牛顿法是两阶的； 3.4 非线性方程组的解法 非线性方程组的牛顿迭代法； 3.4 非线性方程组的解法 电力系统计算机仿真 福州大学 电气工程与自动化学院电力系 邵振国 第三章 线性和非线性方程组的解法 3.1 问题描述 3.2 线性方程组的迭代法 3.3 线性方程组的直接法 3.4 非线性方程组的解法 3.1 问题描述 线性方程组 AX=B A是非奇异系数矩阵，B是右端项； 找到X*，使得AX*=B; 3.1 问题描述 非线性方程组 含有 x 的高次项； 潮流方程； F(X) = 0； 找到X*，使得F(X*) = 0； 3.1 问题描述 线性方程组的解法； 直接法； 迭代法； 非线性方程组的解法； 迭代法； 通常无法直接求解； 3.2 线性方程组的迭代法 3.2.1 简单迭代 3.2.2 Gauss-Seidel迭代 3.2.3 松弛法 3.2.4 迭代的收敛性 等效变换； AX=B，aii ≠ 0 简单迭代法的思路； 3.2.1 简单迭代法 3.2.1 简单迭代法 迭代过程； 构造向量序列 X[k]，使其收敛到X*； 3.2.1 简单迭代法 简单迭代法的特征； X[k] = D-1B - (I - D-1A)X[k-1] 一阶； X[k+1]只和X[k]有关，和X[k-1]等无关； 线性； X[k+1]是X[k]的线性函数； 定常； 迭代式的系数是常数； 3.2.1 简单迭代法 收敛性定义； 如果对于任意选取的初始值X[0]，迭代式X[k] = D-1B -(I -D-1A)X[k-1]生成的序列X[k] 具有相同的极限 X* ，并且 X* 是 AX=B 的解，那么称迭代是收敛的。 构造迭代式要注意的问题； 不是等效变换得来的迭代式都是收敛的； 如何构造收敛的迭代式； 3.2.2 Gauss-Seidel迭代法 Gauss-Seidel迭代法； 计算用 3.2.2 Gauss-Seidel迭代法 Gauss-Seidel迭代式的构造； AX = B； A = D-L-U； 恒等变形：(D-L)X = UX + B； 构造迭代式：(D-L)Xk+1=UXk + B； Xk+1 = (D-L)-1UXk + (D-L)-1B； 优点： 节省内存； 收敛性好； D -L -U 分析用 3.2.3 松弛法 松弛法； w = 1 3.2.3 松弛法 松弛法； Xk+1 = w·(D-1LXk+1 + D-1UXk+D-1B)+(1-w)·Xk； Xk+1 = w·(D-1LXk+1 + D-1UXk+D-1B)+(1-w)·Xk； 计算用 3.2.4 迭代的收敛性 一般迭代法的收敛性分析； 简单迭代法：X[k] = D-1B -(I -D-1A)X[k-1] 记一般迭代式为：X[k] = GX[k-1]＋E 设X[k]收敛于AX=B的解X*，那么X* = GX*＋E； 记误差 ek =X[k]-X*，那么ek = Gek-1 = Gke0； 因为收敛，limek=limGke0=0，所以limGk=0； 定义：G的谱半径S(G)=max|mi|，mi是特征根； 定理：limGk=0 = S(G)1； ∴迭代过程收敛= S(G)1； 收敛速度和S(G)有关； 3.2.4 迭代的收敛性 松弛法的收敛性分析； A = Q -R = (D -wL)/w –[(1-w)D+wU]/w； QX = RX + B； Xk+1 = Q-1RXk + Q-1B； Xk+1 = (D -wL)-1[(1-w)D+wU]Xk + Q-1B； 收敛的必要条件； 0w2； 分析用 3.2.4 迭代的收敛性 如何判断迭代是否收敛； 基本定义； 如果矩阵A不能通过行的次序调换和相应列的次序调换成为 ，A11和A22是方阵，则称A为不可约。 若方阵满足 且至少对应一个i，上式中严格的不等号成立，则称矩阵A具有对角优势。 如果系数矩阵A不可约且具有对角优势，则简单迭代收敛； 如果还有0w1，则松弛法收敛； 3.