计算机数学课件

1、计算机数学,第十一章线性方程组,首页,计算机数学,基本要求、重点难点,11.1线性方程组的消元法,11.2线性方程组解的结构,11.3线性代数的应用实例,11.4演示与实验十,计算机数学,基本要求,理解线性方程组解的概念。 理解齐次线性方程组有非零解的充分必要条件及非齐次线性方程组有解的充分必要条件。 理解齐次线性方程组的基础解系的概念，掌握齐次线性方程组的基础解系和通解的求法。 理解非齐次线性方程组解的结构及通解的概念；掌握用行初等变换求线性方程组的通解的方法，会用特解及相应的导出组的基础解系表示非齐次线性方程组的通解,计算机数学,重点难点,重点： 线性方程组的解的理论与求解方法,计算机数学

2、,11.1 线性方程组的消元法,用消元法解线性方程组的具体做法是：对方程组反复施行初等变换，化为阶梯形方程组，然后从阶梯形方程组中看出原方程组是有惟一解，无穷多解，无解。这种方法称为高斯消元法,上一章中，我们研究了用克拉默法则和逆矩阵求线性方程组的解。它要求方程组必须是n个未知量和n个方程的线性方程组，而且系数行列式不等于零。对于未知量个数和线性方程组的个数不相等，或者相等但系数行列式为零时，上一章的两种方法将无能为力,计算机数学,设含有n个未知量、有m个方程式组成的方程组,其中系数a ij，常数b j 都是已知数，x i是未知量（也称为未知数）。当右端常数项b1,b2 , , b m不全为0

3、 时，称方程组（11.1.13）为非齐次线性方程组；当b1=b2= =b m= 0时，即,11.1.13,称为齐次线性方程组,2,计算机数学,由n个数k1, k2, , k n组成的一个有序数组（ k1, k2, , k n ），如果将它们依次代入方程组（11.1.13）中的x1,x2, ,x n后，（11.1.13）中的每个方程都变成恒等式，则称这个有序数组（k1, k2, , k n）为方程组（11.1.13）的一个解。显然由x1=0, x2=0, , x n=0 组成的有序数组（0, 0, , 0）是齐次线性方程组（2）的一个解，称之为齐次线性方程组（2）的无解，而当齐次线性方程组的未知

4、量取值不全为零时，称之为非零解,非齐次线性方程组（1）的矩阵表示形式为：可用矩阵形式表示为 AX = b,称A为方程组（1）的系数矩阵，X 为未知矩阵，B 为常数矩阵,计算机数学,将系数矩阵A 和常数矩阵B放在一起构成的矩阵,称为方程组（1）的增广矩阵。 齐次线性方程组（2）的矩阵表示形式为：AX = O (11.1.20),定理11.1线性方程组(11.1.13)有解的充要条件是系数矩阵的秩与增广矩阵的秩相等，即r(A)=r(B).对于方程组(11.1.14)中，若b=(0，0，0)T，则方程组为 AX=O， (11.1.20) 称为齐次线性方程组,计算机数学,定理11.2 齐次方程组(11

5、.1.20)一定有解：若r(A)=n则只有零解；它有非零解的充要条件是r(A)n,由上述定理可知，若m是系数矩阵的行数：(1) 当mn 时，r(A)mn，此时方程组(11.1.20)一定有非零解，即齐次方程中未知量的个数大于方程的个数就一定有非零解； (2) 当m=n 时，方程组(11.1.20)有非零解的充要条件是它的系数行列式det A=0 (3) 当m=n且r(A)=n 时，此时系数矩阵的行列式det A0，故方程组(11.1.20)只有零解； (4) 当mn 时，此时r(A)n，故存在方程组(11.1.20)的同解方程组，使“mn,计算机数学,在10.3.2中我们给出了两个特殊的矩阵列

6、向量和行向量，为了讨论方便用小写的希腊字母a，表示，而对于n元线性方程组(11.1.13)，每一个方程的系数都可以看成一个n维行向量，即 i=(ai1，ai2，a in) (i=1，2，m)，共有m个n维行向量a1，a2，am，叫做系数矩阵的行向量组。 每个未知数的系数构成一个列向量,共有n个列向量，称为系数矩阵m 维列向量组。相应地，方程组(11.1.13)的常数项也可以表为一个m 维列向量,可见，线性方程组(11.1.13)与n+1个m 维列向量组1，2， n，之间是一一对应的。可用向量组表示方程组。若方程组有解，即存在x1,x2,x n，使得=x11+x22+x n n成立，此时称是向量

7、组1,2, n的线性组合,11.2 线性方程组解的结构,计算机数学,11.2.1 向量的线性相关性,定义11.1 设a1，a2, m和都是n 维行(列)向量，若存在一组数1,2,am，使得=11+22+ m m， 则称向量是向量组1,2, m的线性组合或称可由向量组1,2, m 线性表出。 由此可得：方程组(11.1.13)有解的充要条件是可由向量组1，2， n 线性表出,例 11.2.1 求证任一n维向量=(a1，a2, n)是向量组1=(1，0，0)，2=(0，1，0)， n=(0，0，1)的线性组合。 证 事实上，令1=a1，2=a2， n=an，则有=(a1，a2，an)=a11+a2

8、2+a n n.即向量是向量组1，2， n的线性组合，或者说任意 n 维向量可由向量组1，2， n 线性表出。向量组1，2， n称为n维单位向量组,定义11.2 设n 维向量1,2, m，若存在一组不全为零的实数1,2, m，使11+22+mm =0 成立，则称向量组1,2, m线性相关。否则，称向量组1,2,am 线性无关,定义11.3 在向量组1,2, m中，若有r 个向量(rm)线性无关，而任意添加一个向量(r 个向量之外还有的话)都是线性相关，则称这 r 个向量构成的部分向量组称为原向量组的极大线性无关组，简称极大无关组,计算机数学,11.2.2 齐次线性方程组解的结构,设线性方程组(

9、11.1.13)，写成矩阵形式如(11.1.14)的形式，当b1=b2=b m=0，即b为零向量时，为齐次方程组(11.1.20)。也称为方程组(11.1.14)的导出方程组。方程组的解是一个列向量X=(x1，x2，x n) T，称为方程组的解向量。 设齐次方程组(11.1.20)有非零解，则它的解有下述性质,计算机数学,定义11.4 设1，2， s是方程组(11.1.20)的一组解向量，并且： (1) 1，2， s 线性无关； (2) 方程组(11.1.20)的任一解向量都可由向量组1，2， s 线性表出。 则称1，2， s 是线性方程组(11.1.20)的一个基础解系,定理11.3 若齐次

10、线性方程组(11.1.20)的系数矩阵A的秩r n (r 0)，那么方程组(11.1.20)有基础解系，且基础解系所含解向量的个数等于n-r,计算机数学,对于齐次线性方程组，其向量方程形式为：Ax=O，它的解向量可用通式表示为：=k11+k22+.+k n-rn-r，(其右端的1,2,n-r都是解向量：若取k1=1，其余的k为0，即可看出1为解向量，.。) 故我们可以说，Ax=0的解向量为某n-r 个线性无关的解向量的线性组合。 注：这任意 n-r 个线性无关的解向量是齐次线性方程组解空间中的一个最大线性无关组。是解空间的一个基,计算机数学,11.2.3 非齐次线性方程组解的结构,齐次线性方程

11、组(11.1.13)的解与它的导出方程组(11.1.20)的解之间有密切的关系，具有以下两个性质,计算机数学,定理11.4 若*是方程组(11.1.13)的一个解，1，2， n-r 是它的导出组(11.1.20)的一个基础解系，则方程组(11.1.13)的全部解为=*+k11+k22+k n-r n-r，(11.2.5)。 其中k1，k2，k n-r 是任意实数,证 先证是方程组(11.1.13)的一个解。事实上，由于 A=A(*+k11+k22+k n-rn-r) = A*+k1A1+k2A2+k n-r An-r=b+O+O=b. 再证方程组(11.1.13)的任意解都可以用式(11.2.

12、5)表示。设是方程组(11.1.13)的任意一个解，则由性质11.3知，-*可由导出组的一个基础解系表出，即有 -*=k11+k22+k n-rn-r， 于是，=*+k11+k22+k n-rn-r,计算机数学,11.3 线性代数的应用实例,11.3.1 线性规划问题,例 11.3.1,计算机数学,线性规划问题的一般形式如下：设有n个变量x1,x2,x n满足,S 称为目标函数，式(11.3.3)称为约束条件。 引入新的非负变量(称为松弛变量)x3,x4,x5就可以使不等式组(11.3.2)变为一组等式。因为2x1+x2比 80 小，加上某个正数量x3，使得它们的和为 80。类似地，也可以使式

13、(11.3.2)的另外两式变为等式，于是，有,计算机数学,显然，满足式(11.3.4)和式(11.3.5)的解x i(i=1,2,3,4,5)中的x1,x2必定满足式(11.3.1)和式(11.3.2)，因此，求满足式(11.3.5)的解x i(i=1,2,3,4,5),使50 x1+30 x2+0 x3+0 x4+0 x5(即式(11.3.4)取最大值，其中的x1、x2就是原线性规划问题的解。因此，我们将公式(11.3.3)改写成等式形式，称为线性规划问题的标准形式，即,其中bi0 (i=1,2,m).满足公式(11.3.6)的x1,x2,x n称为线性规划问题的最优解，相应地 max S=

14、S0 称为该问题的最优值,计算机数学,11.3.2 线性规划问题的初等解法,如果把S 亦视为一个变量，公式(11.3.6)写为,计算机数学,11.4 演示与实验十,11.4.1 实验目的,1. 学习用Mathematica判定非齐次线性方程组解的存在性； 2. 学习用Mathematica求齐次线性方程组的基础解系和通解； 3. 学习用Mathematica求非齐次线性方程组的通解和特解,计算机数学,11.4.2 内容与步骤,1. 用Mathematica判定非齐次线性方程组解的存在性。 根据线性方程组解的存在性定理，只要求出系数矩阵和增广矩阵的秩，即可判定方程组的解是否存在。 求矩阵的秩，除可以用10.6.2中的方法外，还可以用下面的命令：n-LengthNullSpaceA其中，n是矩阵A的列数，LengthNullSpaceA是齐次线性方程组AX=O的基础解系所含解的个数,计算机数学,2. 用Math