数值分析4 高斯主元素消去法.docVIP

§2高斯主元素消去法 引例（小主元） 【】 Gauss消去法解方程组 解：clear all a=[0.0120 1.00 2.00;1.00 2.63 5.24;-2.00 1.04 4.57]; b=[1.00;2.00;3.00];x=a\b 方程组的三位有效数字的解： Gauss消去法求解(取三位有效数字)： 解出。 【注】1）设Ax=b,其中A为n阶非奇异矩阵，可以应用高斯消元法。 2）消元过程中，即使，用其作除数会导致计算中间结果数量级严重增长和舍入误差的累积、扩大，最后使得计算结果不可靠。 3)应避免采用绝对值很小的主元素；对一般的系数矩阵，最好保持乘数，因此，在高斯消去法中应引进选主元技巧，以便减少计算过程中舍入误差对求解的影响。 clear all a=[0.0120 1.00 2.00;1.00 2.63 5.24;-2.00 1.04 4.57]; b=[1.00;2.00;3.00];x_value=vpa(a\b,15) %10位有效数字的近似解 a=[a,b];eps=1e-6;[n,m]=size(a);Gauss,x=vpa(x,15) %对比高斯消去法的结果 一、列主元素消去法 基本思想：在每轮消元之前，选列主元素（绝对值最大的元素）,使乘数（即消元因子） 步骤：设已进行k-1轮消元，得矩阵 一落千丈 1 23 Ｓ１：选列主元素： ……………（1） Ｓ２：换行：如果→０，则方程组解不唯一，停止运算； 否则，如果i0=k， 则可进行下一轮消元； 如果，则ri0 rk，然后进行下一轮消元。 算法与例子 步1，定义数组A，B，X；输入A，B，EPS 步2，对k=1，2，…，n-1做（2.1）~(2.5)： （2.1）求i0使 ， （2.2）如果小于某小数EPS，则输出无唯一解信息，停机， （2.3）如果i0=k，则做（2.5）， （2.4）交换i0，k两行对应元素， （2.5）对i，j=k+1，…，n消元计算： 步3，回代求解： 步4，输出结果。 【】回代解出：。 程序:见column_gauss.m a=[0.0120 1.00 2.00;1.00 2.63 5.24;-2.00 1.04 4.57]; b=[1.00;2.00;3.00];eps=1e-6; 二、行主元素消去法 （即在第k行选主元素） 在第k行选主元素： ……………(2) 三、完全主元素消去法 ( 即选主元素的范围为k≤i，j≤n) 选主元素: ……………（3） 【注】完全选主元消去法解方程组，在选主元素时要化费较多的机器时间，行主元消去法与列主元消去法运算量大体相同，实际计算时，用列主元消去法即可。 四、列主元高斯—约当消去法（Gauss-Jordam） 高斯消去法解方程组自始至终仅仅是对A(k)的第k行下面的元素进行消元计算。现考虑一个修正方法，即消元计算对A(k)的第k行上面、下面的元素都进行消元计算，最后不需要回代即可求得方程组的解。这是高斯—约当消去法。再引进按列选主元，就是列主元高斯—约当消去法。 对于第k步 【】 clear all a=[0.0120 1.00 2.00;1.00 2.63 5.24;-2.00 1.04 4.57]; b=[1.00;2.00;3.00];a=[a,b] a=rref(a) 得方程组的三位有效数字的解： 【】 解： 所以， Matlab中解法 a=[1 2 3;2 4 5;3 5 6];a=[a,eye(3)];rref(a)