计算方法讲义课件 四 线性方程组.docVIP

第四章 线性方程组 许多科学技术问题可归结为有多个未知量x1, x2, …, xn的线性方程组：，这里aij为系数。方程组的矩阵形式为：AX = B，其中。 线性方程组的数值解法有直接法和迭代法两大类。直接法不考虑舍入误差，通过有限步骤四则运算即能求得线性方程组的解。用克莱姆法则来求解线性代数方程组并不实用。 4.1高斯消去法 高斯消去法思想是通过逐步消元，把方程组化为系数矩阵为三角形，用回代法解此三角形方程组。 三角形方程组 三角形方程组分：下三角和上三角两类。 若aii1 0，i = 1, 2,…, n，则下三角的解为：，k = 2, 3,…, n，称为前推。 若aii 1 0, i = 1, 2,…, n，则上三角解为，k = n-1, n-2,…, 1，称为回代。 高斯消去法的过程实例 例 求解。（I）乘（）后加到（II）上，（I）乘（）后加（III）上，可消去（II）、（III）中x1，。 （II）乘（）后加（III），得：，回代得x3 = 2， x2 = 8，x1 = -13。 将原方程组转化为：，回代公式为： 。高斯消去法分为消元与回代过程。算法复杂度为n3。 例 求解。精确解为：，。用消去法第一步以0.0001为主元，从第二个方程中消x1后可得：，回代可得x1 = 0.00，注意这里假定计算过程中有舍入误差（考虑只保留三位有效数字）。这个解误差很大。误差分析表明：顺序消元法在系数矩阵A为对称正定时，可保证对舍入误差的数值稳定性。 在列主元消去法中，把各方程中要消去的那个未知数的系数按绝对值最大值作为主元素。 例：求解。 第一步将4选为主元素，然后将主元行换到第一行得到 消元归结到： ，回代得。 例：求解。解释。 例：求解。解释。 计算步骤： 输入：阶数，增广矩阵， 输出： 1．对 （1）按列选主元：选取 使 ， （2）如果，交换的 第k行与底l行元素， （3）消元计算：… 2．回代。 4．2追赶法 在实际中常遇到以下形式方程组：称为三对角矩阵。这种方程组有效算法—追赶法。追赶法是高斯消去法的简化，分消元与回代两过程。将原系数矩阵转化为：，逐步回代，即可依次求出xn，xn-1，…，x1。 综合追赶的过程，得如下计算公式：，。 4．3 初等变换的应用 高斯消去法实际上就是初等变换。对行的初等变换有三种：（1）某行乘一非零系数，（2）两行位置对调，（3）某行乘系数后加到另一行。同样有对列的初等变换。 定义： n阶单位矩阵I经过一次初等变换得到的矩阵称为初等矩阵： 对换单位矩阵I的i, j两行，所得初等矩阵记为I(ri,rj)。 用非零数k乘单位矩阵I的第i行，所得初等矩阵记为I(kri)。 把单位矩阵I的第i行的k倍加到第j行上，所得初等矩阵记为I(kri+rj)。 举例说明三种初等变换阵。 定理：对一个矩阵A作一次初等行变换所得到的矩阵B，等于一个对应的m阶初等矩阵左乘矩阵A；对A作一次初等列变换所得到的矩阵C，等于一个对应的n阶初等矩阵右乘矩阵A。 验证：= 应用1：逆矩阵： 设A = (aij)n′n是非奇，|A| 1 0，由于AA-1 = I，求A-1的问题相当于解下列线性方程组：。 例 求逆。，逆为。 应用2：行列式值： 行列式任意一行(列)的元素乘以同一个数后，加到另一行（列）的对应元素上，其行列式的值不变；任意对换两行（列）的位置其值反号；三角矩阵的行列式之值等于其主对角元素的乘积。因此，可用高斯消去法将 |A|化成： 应用3：秩分解 令A为矩阵，，，用单位阵记录变换过程。 4.4 矩阵LU分解 用矩阵理论来分析高斯消去法，实现矩阵的三角分解。 定理：上三角乘上三角还是上三角；下三角乘下三角还是下三角。 定理：非奇异下（上）三角矩阵的逆还是下（上）三角。 证明的方法一初等变换法，二是矩阵分块、归纳法。 初等变换法：设是单位下三角，考虑，都是下三角，就是的逆。 归纳法：。 考虑线性方程组。假设解此方程组用高期消去法能够完成（不进行交换两行的初等变换）。由于对施行初等变换相当于用初等矩阵左乘，所以用高斯消去法可以将矩阵分解为：。令，则，其中为单位下三角矩阵，为上三角矩阵。 定理：（矩阵LU分解）(LU Decomposition)设A为n′n实矩阵，如果解AX=B用高斯消去法能够完成（限制不进行行的交换），则矩阵A可分解为单位下三角矩阵L与上三角矩阵U的乘积：A = LU，且这种分解是唯一的。 证明：证唯一性。设 A = L1 U1 = LU，其中L1，L1为单位下三角阵，U1，U为上三角阵。由于和存在，于是有，上式右端为上三角矩阵，左边为单位下三角阵，故为单位矩阵。即L1 = L，U1=U。 例如对于一个的矩阵，就有 定理：约化主元素（i = 1, 2, …, k）