.
高斯消元法解线性方程组
在工程技术和工程管理中有许多问题经常可以归结为线性方程组类型的数学模型,这些模型中方程和未知量个数常常有多个,而且方程个数与未知量个数也不一定相同。那么这样的线性方程组是否有解呢?如果有解,解是否唯一?若解不唯一,解的结构如何呢?这就是下面要讨论的问题。
一、线性方程组
设含有n个未知量、有m个方程式组成的方程组
?a11x1?a12x2???a1nxn?b1?ax?ax???ax?b?2112222nn2 ? (3.1)
?????????am1x1?am2x2???amnxn?bm其中系数aij,常数bj都是已知数,xi是未知量(也称为未知数)。当右端常数项b1, 称方程组(3.1)为非齐次线性方程组;当b1=b2= … =bm= b2, …, bm不全为0时,0时,即
?a11x1?a12x2???a1nxn?0?ax?ax???ax?0?2112222nn ? (3.2)
?????????am1x1?am2x2???amnxn?0称为齐次线性方程组。
由n个数k1, k2, …, kn组成的一个有序数组(k1, k2, …, kn),如果将它们依次代入方程组(3.1)中的x1, x2, …, xn后,(3.1)中的每个方程都变成恒等式,则称这个有序数组(k1, k2, …, kn)为方程组(3.1)的一个解。显然由x1=0, x2=0, …, xn=0组成的有序数组(0, 0, …, 0)是齐次线性方程组(3.2)的一个解,称之为齐次线性方程组(3.2)的零解,而当齐次线性方程组的未知量取值不全为零时,称之为非零解。
(利用矩阵来讨论线性方程组的解的情况或求线性方程组的解是很方便的。因此,我们先给出线性方程组的矩阵表示形式。) 非齐次线性方程组(3.1)的矩阵表示形式为: AX = B
其中
?x1??b1??a11a12?a1n??x??b??aa22?a2n?221?,X = ??,B = ?2? A = ???????????????????xaa?a?n??bn?m2mn??m1称A为方程组(3.1)的系数矩阵,X为未知矩阵,B为常数矩阵。将系数矩阵A和常数矩阵B放在一起构成的矩阵
.
.
?a11a12?a1nb1??a?a?ab21222n2? [AB]=??????????aa?abm2mnm??m1称为方程组(3.1)的增广矩阵。
齐次线性方程组(3.2)的矩阵表示形式为:AX = O
二、高斯消元法
(下面介绍利用矩阵求解方程组的方法,那么矩阵初等行变换会不会改变方程组的解呢?我们先看一个定理。)
定理3.1 若用初等行变换将增广矩阵[AB]化为[CD],则AX = B与CX = D是同解方程组。
证 由定理3.1可知,存在初等矩阵P1, P2, …, Pk,使 Pk…P2P1(AB) = (CD) 记Pk…P2P1 = P,则P可逆,即P?1存在。 设X1为方程组A X = B的解,即 AX1 = B 在上式两边左乘P,得 P AX1 = PB 即 CX1= D 说明X1也是方程组C X = D的解。反之,设X2为方程组C X = D的解,即 CX2= D 在上式两边左乘P?1,得 P?1CX2= P?1D 即 AX2 = B 说明X2也是方程组AX = B的解。 因此,方程组A X = B与C X = D的解相同,即它们是同解方程组。(证毕)
(由定理3.1可知,求方程组(3.1)的解,可以利用初等行变换将其增广矩阵[AB]化简。又有第二章定理2.10可知,通过初等行变换可以将[AB]化成阶梯形矩阵。因此,我们得到了求解线性方程组(3.1)的一般方法:)
用初等行变换将方程组(3.1)的增广矩阵[AB]化成阶梯形矩阵,再写出该阶梯形矩阵所对应的方程组,逐步回代,求出方程组的解。因为它们为同解方程组,所以也就得到了原方程组(3.1)的解。这种方法被称为高斯消元法,
(下面举例说明用消元法求一般线性方程组解的方法和步骤。)
.
.
?x1?x2?2x3?x4??1?x?5x?3x?2x?0?1234 例1 解线性方程组 ? (3.3)
3x?x?x?4x?21234????2x1?2x2?x3?x4?1 解 先写出增广矩阵[AB],再用初等行变换将其逐步化成阶梯形矩阵,即
1?2?1?1?②?①(?1)?11?2?1?1??1?1?③?①(?3)?04?1?11?5?3?20?①2??④? ?????[AB]=??3?11?0?4742?75??????221?1104?3?3?1?????11?2?1?1??11?2?1?1??04?1?11?④?③(1)?04?1?11?③?②④?②(?1)3????? ??????????006?00666?66?????00?2?2?200000????
上述四个增广矩阵所表示的四个线性方程组是同解方程组,最后一个增广矩阵表示的线性方程组为
?x1?x2?2x3?x4??1?4x2?x3?x4?1 ??6x3?6x4?6?1将最后一个方程乘,再将x4项移至等号的右端,得
6x3??x4?1
将其代入第二个方程,解得
x2?12
再将x2,x3代入第一个方程组,解得
x1??x4?12
因此,方程组(3.3)的解为
?x1??x4?12? ?x2?12 (3.4)
?x??x?14?3其中x4可以任意取值。
由于未知量x4的取值是任意实数,故方程组(3.3)的解有无穷多个。由此可知,表示式(3.4)表示了方程组(3.3)的所有解。表示式(3.4)中等号右端的未知量x4称为自由未知量,用自由未知量表示其它未知量的表示式(3.4)称为方程组(3.3)的一般解,当表示式(3.4)中的未知量x4取定一个值(如x4=1),得
11到方程组(3.3)的一个解(如x1??,x2?,x3?0,x4?1),称之为方程
22组(3.3)的特解。
注意,自由未知量的选取不是唯一的,如例1也可以将x3取作自由未知量。
.