高斯消元法解线性方程组
在工程技术和工程管理中有许多问题经常可以归结为线性方程组类型的数学模型,这些模型中方程和未知量个数常常有多个,而且方程个数与未知量个数也不一定相同。那么这样的线性方程组是否有解呢?如果有解,解是否唯一?若解不唯一,解的结构如何呢?这就是下面要讨论的问题。
一、线性方程组
设含有n个未知量、有m个方程式组成的方程组
?a11x1+a12x2+"+a1nxn=b1?ax+ax+"+ax=b?2112222nn2 ?(3.1) """"""???am1x1+am2x2+"+amnxn=bm
(也称为未知数)。当右端常数项b1, 其中系数aij,常数bj都是已知数,xi是未知量
称方程组(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,
(0, 0, …, 0)是齐次线性方程组(3.2)的一个解,x2=0, …, xn=0组成的有序数组
称之为齐次线性方程组(3.2)的零解,而当齐次线性方程组的未知量取值不全为零时,称之为非零解。
(利用矩阵来讨论线性方程组的解的情况或求线性方程组的解是很方便的。因此,我们先给出线性方程组的矩阵表示形式。)
非齐次线性方程组(3.1)的矩阵表示形式为:
AX = B
其中
?x1??b1??a11a12"a1n??x??b??aa22"a2n?212?,X = ??,B = ?2? A = ??#??#??""""???????aaax"m2mn??m1?n??bn?
称A为方程组(3.1)的系数矩阵,X为未知矩阵,B为常数矩阵。将系数矩阵A和常数矩阵B放在一起构成的矩阵
1
?a11a12"a1nb1???abaa"212222n? [AB]=??"""""???baa"am2mnm??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)的解。这种方法被称为高斯消元法,
(下面举例说明用消元法求一般线性方程组解的方法和步骤。)
2
?x1+x2?2x3?x4=?1?x+5x?3x?2x=0?1234 例1 解线性方程组 ? (3.3) 23xxx4x?++=1234????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项移至等号的右端,得 6
x3=?x4+1
将其代入第二个方程,解得
x2=2
再将x2,x3代入第一个方程组,解得
x1=?x4+2
因此,方程组(3.3)的解为
?x1=?x4+2? (3.4) ?x2=2
?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取作自由未知量。 3
如果将表示式(3.4)中的自由未知量x4取一任意常数k,即令x4= k,那么方程组(3.3)的一般解为
?x1=?k+2?x=2?2 ,其中k为任意常数。 ?x=?k+1?3
??x4=k
用矩阵形式表示为
?x1???k+2???1??2??x????0??2?22?=k??+?? (3.5) ??=??x3???k+1???1??1?????????kx??1??0??4??
其中k为任意常数。称表示式(3.5)为方程组(3.3)的全部解。
(用消元法解线性方程组的过程中,当增广矩阵经过初等行变换化成阶梯形矩阵后,要写出相应的方程组,然后再用回代的方法求出解。如果用矩阵将回代的过程表示出来,我们可以发现,这个过程实际上就是对阶梯形矩阵进一步简化,使其最终化成一个特殊的矩阵,从这个特殊矩阵中,就可以直接解出或“读出”方程组的解。例如,)对例1中的阶梯形矩阵进一步化简,
?11?2?1?1?③1?11011?6?04?1?11?①+③2?04002?+③??②? ???→??006?00111?66?????00??000?00000?
?10012?1②?01002?4①+②(?1)? ????→??00111????00000?
上述矩阵对应的方程组为
?x1+x4=2??x2=2
?x+x=14?3
将此方程组中含x4的项移到等号的右端,就得到原方程组(3.3)的一般解,
?x1=?x4+2? (3.4) ?x2=2
?x=?x+14?3
其中x4可以任意取值。
?x1+2x2?3x3=4?2x+3x?5x=7?123 例2 解线性方程组 ? x+x?x=439923?1
??2x1+5x2?8x3=8
4
解 利用初等行变换,将方程组的增广矩阵[AB]化成阶梯阵,再求解。即
?12?34??12?34??0?11?1??23?57?? ?→?[AB]=??0?53?7??43?99???????01202588????
?12?34??12?34??0?11?1??0?11?1??→?? →??00?2?2??0011???????00110000????
?1003??1207??0102??0102?? ?→? →??0011??0011?????00000000????
一般解为
?x1=3??x2=2 ?x=1?3
?x1+x2+x3=1? 例3 解线性方程组 ??x1+2x2?4x3=2 ?2x+5x?x=3123?
解 利用初等行变换,将方程组的增广矩阵[AB]化成阶梯阵,再求解。即
?1111??1111??→?03?33? [AB]=??12?42?????????25?13??03?31?
1??111? 03?33 →??????000?2?
阶梯形矩阵的第三行“0, 0, 0, -2”所表示的方程为:0x1+0x2+0x3=?2,由该方程可知,无论x1,x2,x3取何值,都不能满足这个方程。因此,原方程组无解。
三、线性方程组的解的判定
前面介绍了用高斯消元法解线性方程组的方法,通过例题可知,线性方程组的解的情况有三种:无穷多解、唯一解和无解。从求解过程可以看出,方程组(3.1)是否有解,关键在于增广矩阵[A B]化成阶梯非零行的行数与系数矩阵A化成阶梯形矩阵后非零行的行数是否相等。因此,线性方程组是否有解,就可以用其系数矩阵和增广矩阵的秩来描述了。
5
(2) 用初等行变换将增广矩阵化成阶梯阵,即
2?3?11??1?12?3?11????1?12?01?1?4?7?→ … →?? [A B ]=??2?31?0006?0??????31250000????
因为 r(AB)=r(A)=2<n(= 3),所以方程组有无穷多解。
(3) 用初等行变换将增广矩阵化成阶梯形矩阵,即
2?3?11??1?12?3?11????1?11?01?2?4?7?→… →?? [A B ]=??2?31?00?70?6??????31250000????
因为 r(AB)=r(A)= 3 = n,所以方程组有唯一解。
例5 判别下列齐次方程组是否有非零解? (机动)
?x1+3x2?7x3?8x4=0?2x+5x+4x+4x=0?1234 ?37230?x?x?x?x=1234???x1+4x2?12x3?16x4=0
解 用初等行变换将系数矩阵化成阶梯形矩阵,即
?7?8?3?13?7?8??1???0?118?254420? ?→?A =??02?23?27???3?7?2?3???????14?12?160158????
?13?7?8??13?7?8??0?11820??0?11820??→?? →??001313??001313?????000?1001312????
因为 r(A)= 4 = n,所以齐次方程组只有零解。
7
www.99jianzhu.com/包含内容:建筑图纸、PDF/word/ppt 流程,表格,案例,最新,免费下载,施工方案、工程书籍、建筑论文、合同表格、标准规范、CAD图纸等内容。