数值分析与方法03.pptVIP

最后得, 得到与原方程组等价的上三角形方程组。 (2)回代过程 只要ann(n)≠0即可回代求解： 3、高斯消去法的条件 高斯消去法要求aii(i)≠0(i=1,2,…,n) Th3.1 方程组Ax=b能用顺序高斯消去法求解=A的各阶顺序主子式都不为0. 若消元过程中允许对增广矩阵进行行交换,则方程组Ax=b可用消元法求解=A可逆. 4、高斯消去法的计算量 n3/3 二、选主元消去法 高斯消元法存在两个问题： (1) 若akk(k)=0, 则顺序消元过程无法实现; (2) 即使akk(k)≠0, 若它与第k行到第n行的其他元素相 比绝对值较小,这样会造成舍入误差的传播,使解失真. 1、列主元消去法 (1) 按列选主元 设方程组Ax=b的增广矩阵为 在A(1)第一列中选取绝对值最大的元素作主元. 若i1≠1，则交换增广矩阵的第一行与第i1行。 经第一次消元得 重复上述过程，设已完成第k-1次消元(1≤k≤n-1) 原方程组变为A(k)x=b(k)。其增广矩阵为： 第k次消元：先选主元及行交换：在方框内的诸元素aik(k）中选出绝对值最大者，即确定ik，使得 (若|aik,k(k)|=0, 说明Ax=b无确定解) 若ik≠k，则交换第ik行和第k行元素，然后进行消元。...这样做下去，直至k=n-1。 然后进行回代过程。 (2) 算法实现 算法组成: 1) 选主元. 2) 换行. 3) 消元. 4) 回代. 存贮方式: 二维数组A(n,n)存放矩阵A的元素 一维数组b(n)存放常向量b的元素 对k=1,2,…,n-1做到步5： 步1、按列选主元素 步2、若aik,k=0，则输出‘A奇异’；停机。 步3、若ik=k，则转步4；否则 换行：akj?aik,j bk?bik (j=k,k+1,…,n) 步4、计算乘子mik：aik←mik=aik/akk (i=k+1,…,n) 步5、消元计算： aij←aij-mikakj (i, j=k+1,…,n) bi←bi-mikbk (i=k+1,…,n) 步6、回代求解：bn←bn/ann (i=n-1,n-2,…,1) 2、全主元消去法 在(1.8)中，若每次选主元不局限于第k列的方框内，而在整个主子阵 中选取，这种消去法为全主元消去法 (1) 选主元过程 与列主元消去法相比增加的步骤为： 1）确定ik, jk使 若aik jk=0, 则‘A奇异’，停止计算。否则做2）。 2）做如下行变换、列变换 行变换： akj(k)?aik,j(k) (k≤j≤n) bk(k)?bik(k) 列变换：aik(k)?ai,jk(k) (k≤i≤n) (2) 算法实现 算法组成：1) 选主元. 2) 交换行、列, 记录未知数的次序. 3) 消元. 4) 回代. 5) 调整未知数次序. 存贮方式: 用一维数组IZ(n)记录未知数的次序. 高斯全主元消去法比列主元消去法的工作量大,但数值稳定性更好. 第二节 矩阵三角分解法 3.2.1、LU分解法 若A的顺序主子式都不为0，则 A=LU 其中， 单位下三角阵 非奇异上三角阵 则方程组Ax=b转化为 LUx=b 令 Ux=y， 则Ly=b 则求解Ax=b等价于求解 Ux=y， Ly=b 因为A=LU，由矩阵乘法得 a1j=u1j , j=1,2,…,n ai1=li1u11 i=2,…,n 可推出 u1j=a1j, j=1,2,…,n li1=ai1/u11, i=2,…,n 即求出了U的第一行和L的第一列元素 设已求出U的第一行到第r-1行元素和L的第一列到第r-1列元素, 则 矩阵A的LU分解也称为Doolittle分解。 Doolittle分解法解方程组Ax=b的算法： 计算量n3/3 3.2.2、平方根法 1、平方根法 定理3.5：若A的各阶顺序主子式非0，则A可分解为A=LDU，其中，L是单位下三角阵，U是单位上三角阵，D是对角阵，且这种分解唯一。 定理3.6：设A为对称正定矩阵，则存在三角分解A=LLT，其中L是非奇异下三角矩阵，且当限定L的对角线元素为正时，这种分解唯一。 证明：A对称正定，则A的各阶顺序主子式都为正 由定理2.5，