[自然科学]线性方程组的直接解法1.ppt

yfnie@nwpu.edu.cn 线性方程组的直接解法 概述 高斯消去法 矩阵分解及其在解线性方程组中的应用 矩阵的条件数和方程组的性态 三角分解（2） 紧凑格式 * * 1.1 数值解法的必要性 求： 的解 的值，根据克莱姆(Gramer)法则可 表示为两个行列式之比： 1 概述 如： ,若用每秒万亿次（1012）浮点乘法运算的计算机（当前国内运算速度最快）全天工作，完成这些计算约需30年。若使用一般的个人电脑，每秒完成十亿次（109）浮点乘法运算，则完成这些计算约需3万年。 计算一个 阶行列式需要做 个乘法，求解上述方程共做 次乘除法。 研究数值方法的必要性 1、直接法： 只包含有限次四则运算。假定在计算过程中不产生舍入误差，计算结果是原方程组的精确解。 2、迭代法： 基于一定的递推公式建立逼近方程组近似解的序列。收敛性是其前提，然后还有收敛速度以及误差估计等问题。 Remark：由于运算过程中舍入误差的确存在，实际上直接方法得到的解也是方程组的近始解。 1.2 线性代数方程组的解法分类 设有线性方程组 为非奇异矩阵. 或写为矩阵形式 ，其中 2 高斯消去法 2.1 高斯顺序消去法 若 令 　　 ,用 乘 第一个方程加到第 个方程 ，并保留第一式，得 Step1： Gauss顺序消去法（续） 乘除法运算次数 若 令 ， 乘第二个方程加到第i个方程 ，得 用 Step2: Gauss顺序消去法 记为 （续2） 乘除法运算次数 设为 第k-1步消元完成后，有 （续3） 用- 来乘以第k个方程，加到第i个方程，并保留第k个方程， 得： （i=k+1,k+2,…n） 若 ,令 , （i=k+1,k+2,…n) Step k： Step K （续4） (i,j=k+1,…,n) （i=k+1,…,n） 乘除法运算次数 做完n-1步，原方程可化为同解的上三角方程组。 Gauss顺序消去法 （续5） 回代过程： （i=n-1,n-2,…2，1） 回代乘除运算次数： 在Gauss顺序消去法的消去过程中，若将右端列向量视为方程组A的第n＋1列，直接对该增广矩阵A进行行初等变换，将其变换为上三角矩阵，从而回代求解得原方程组的解。 计算量 Gauss顺序消去法消去过程所需的乘除运算次数为 总的乘除运算次数： 取n＝20，则N＝3060 在消去过程中采用Gauss-Jordan消去法，即将方程组化为对角形方程组，进而化为单位阵，则右端列向量就化为方程组的解向量。 该方法不需回代过程，但总的计算量为 n3/2+n2-n/2， 比Gauss顺序消去法有所增加。 Gauss－Jordan消去法 消去过程能够实施的条件是 回代过程条件增加条件 命题1 定义 Gauss顺序消去法的消元条件 高斯顺序消去法仅对原增广矩阵作行初等变换，故 命题证明 续 定理2： 若系数矩阵A的各阶顺序主子式均不为零则可以用高斯顺序消去法求解方程组A X=b。 定理1：高斯顺序消去法求解方程组A X=b的消元过程能够实施的充要条件是系数矩阵A的前n-1个顺序主子式均不为零。 ①用Gauss顺序消去法求解线性方程组的条件强于解的存在唯一性条件 斯顺序消去法能进行下去，但 或相对于 比较小时，计算产生的舍入误差将导致计算结果误差急剧增大，计算解与真解相差甚远，即该方法不稳定。 Gauss顺序消去法的不足 ② ，所有主元不为零时高 小主元是不稳定的根源，需要采用“选主元”技术，即选取绝对值较大的元素作为主元。 就是原来的 ,然后 2.2 选主元的高斯消去法 列主元消去法 增广矩阵 在第一列中选取绝对值最大的元素,设 = 调换第一行与第p行，这时的 再执行消去法的第一步 再考虑 右下角矩阵，在第二列选取绝对值最大的元素作为主元素,将其所在行与第二行交换，消元。依此类推直至将方程组化成上三角方程组，再进行回代就可求得解。 列主元消去法（续） 优点：只要系数矩阵非奇异，消元过程就能进行，并且主元相对较大，方法稳定好。 在系数矩阵中选取绝对值最大的元素作为主元素，如 然后交换增广矩阵的第一