线性方程组的矩阵求解算法
如何利用数学中的矩阵进行线性方程组的求解
如何利用数学中的矩阵进行线性方程组的求解线性方程组在数学中具有重要的应用价值，求解线性方程组是数学中的基本问题之一。
矩阵是求解线性方程组的有力工具，能够简化计算过程并提高求解效率。
本文将介绍如何利用数学中的矩阵进行线性方程组的求解。
一、矩阵的定义和基本性质矩阵是由数个数按一定规则排列形成的矩形数组。
矩阵可以表示为一个大写字母加上两个下标，例如A，其中A是矩阵的名称，下标表示矩阵的行数和列数。
矩阵的加法和乘法是指对应元素的加法和乘法运算。
矩阵加法要求两个矩阵具有相同的行数和列数；矩阵乘法要求第一个矩阵的列数等于第二个矩阵的行数。
二、线性方程组和矩阵表示线性方程组是一组线性等式的集合。
一个线性方程组可以用矩阵表示，其中系数矩阵是一个m行n列的矩阵，m表示方程组的数量，n 表示未知数的数量；向量b是一个m行1列的矩阵，称为常数向量；向量x是一个n行1列的矩阵，称为未知向量。
线性方程组可以写成Ax=b的形式。
三、矩阵求解线性方程组的方法1. 列主元高斯消元法列主元高斯消元法是一种求解线性方程组的基本方法。
具体步骤如下：(1) 首先将线性方程组写成增广矩阵的形式[A|b]。
(2) 选择第一列中绝对值最大的元素作为主元所在行，将该行与第一行交换。
(3) 将第一行乘以一个系数，使得主元所在列的其他元素都变为0。
(4) 重复第二步和第三步，直到将整个矩阵化为上三角矩阵。
(5) 从最后一行开始，倒序回代求解线性方程组。
2. 矩阵逆的方法如果矩阵A可逆，则可以用逆矩阵来求解线性方程组。
逆矩阵是指与原矩阵相乘得到单位矩阵的矩阵。
具体步骤如下：(1) 首先求出矩阵A的逆矩阵A^(-1)。
(2) 将线性方程组写成矩阵形式Ax=b。
(3) 两边同时左乘A^(-1)，得到x=A^(-1)b。
3. 矩阵的LU分解LU分解是将矩阵A分解为两个矩阵L和U的乘积的过程。
L是一个下三角矩阵，U是一个上三角矩阵。
具体步骤如下：(1) 首先将矩阵A写成增广矩阵的形式[A|b]。
如何解决数学中的方程组与矩阵问题
如何解决数学中的方程组与矩阵问题在数学中，方程组与矩阵问题是常见且重要的内容，解决这些问题需要一定的方法和技巧。
本文将介绍几种解决数学中方程组与矩阵问题的方法，帮助读者更好地理解和应用。
一、高斯消元法高斯消元法是一种用于求解线性方程组的方法，它通过矩阵变换将方程组转化为一个更简单的形式，从而找到解。
下面以一个具体的例子来说明高斯消元法的步骤：假设有如下的方程组：(1) 2x + 3y - z = 7(2) x - y + z = 2(3) 3x - 4y + 2z = 4首先将方程组写成增广矩阵的形式：[ 2 3 -1 | 7 ][ 1 -1 1 | 2 ][ 3 -4 2 | 4 ]接下来，通过一系列的行变换，使矩阵变为上三角矩阵：[ 2 3 -1 | 7 ][ 0 -5 3 | -10 ][ 0 0 3 | -3 ]然后，从最后一行开始，依次求出未知数的值。
首先可以得到 z = -1，再依次代入前面的方程中，求解出 y = 2 和 x = 1。
因此，方程组的解为 x = 1，y = 2，z = -1。
高斯消元法可以帮助我们快速求解线性方程组，但在实际应用中，需要注意矩阵的可逆性和唯一解。
二、矩阵求逆在某些情况下，我们需要求解一个矩阵的逆矩阵，以便更便利地解决方程组或其他相关问题。
矩阵求逆的方法有多种，这里介绍其中一种常见的方法——伴随矩阵法。
对于一个 n 阶方阵 A，如果存在一个 n 阶方阵 B，使得 AB = BA = I，其中 I 是 n 阶单位矩阵，则称矩阵 B 为 A 的逆矩阵，记作 A^-1。
那么如何求解一个矩阵的逆矩阵呢？下面以一个 2 阶方阵为例来说明：首先，假设有一个 2 阶方阵 A：[ a b ][ c d ]如果 A 的行列式不等于 0，即 ad - bc ≠ 0，那么 A 的逆矩阵存在。
为了求解 A 的逆矩阵，我们可以按照以下步骤进行：1. 计算 A 的行列式的值 det(A) = ad - bc。
矩阵求方程的解
矩阵求方程的解
矩阵可以被用来求解线性方程组。
线性方程组可以表示为以下形式：
A * x = b
其中，A 是一个系数矩阵，x 是未知向量，b 是已知向量。
矩阵求解线性方程组主要有两种方法：逆矩阵法和高斯消元法。
1.逆矩阵法：如果矩阵A 是可逆的（即行列式不等于零），
则可以通过以下公式求解线性方程组的解：
x = A⁻¹ * b
其中，A⁻¹ 表示矩阵 A 的逆矩阵，* 表示矩阵的乘法运算。
2.高斯消元法：高斯消元法是通过变换线性方程组的形式，
将其转化为上三角形式或者简化行阶梯形式。
然后，可以
通过回代的方式求解线性方程组的解。
具体步骤如下：
•用初等行变换将矩阵A 转化为上三角形式（或简化行阶梯形式）。
•根据变换后的矩阵形式，可以直接得到解的结果或通过回代得到解。
需要注意的是，在实际应用中，矩阵方程的求解可能会遇到多解、无解或条件问题等情况。
因此，在使用矩阵求解线性方程组时，需要对方程组的性质进行仔细分析，并进行适当的处理。
线代矩阵求解题技巧
线代矩阵求解题技巧线性代数是数学中的一个重要分支，广泛应用于科学和工程学科中。
矩阵求解是线性代数中的一个基本概念，它是解线性方程组、求特征值和特征向量等问题的重要工具。
下面将介绍一些线性代数矩阵求解的基本技巧。
1. 高斯消元法高斯消元法是求解线性方程组的常用方法之一。
该方法的基本思想是通过矩阵变换将线性方程组化为上三角形方程组或者行最简形式，从而得到方程组的解。
高斯消元法具体步骤如下：（1）将线性方程组写成增广矩阵的形式；（2）选取一个主元（通常选取主对角线上的元素），并通过一个变换将该元素下面的所有元素置零；（3）对主元元素下面的行执行类似的操作，直到所有元素都变为零或者上三角矩阵形式；（4）回代求解未知数。
2. LU分解LU分解是将一个矩阵分解为下三角矩阵L和上三角矩阵U的乘积的方法。
这个方法通常用于解决多次使用相同矩阵求解线性方程组的场景。
LU分解的具体步骤如下：（1）设一个n阶方阵A，将其分解为A=LU；（2）通过高斯消元法将A化为上三角矩阵U；（3）构造下三角矩阵L，使得A=LU成立。
3. 矩阵的逆和伴随矩阵对于一个可逆矩阵A，可以通过求解逆矩阵来求解线性方程组。
设A为n阶可逆方阵，若存在一个n阶矩阵B，满足AB=BA=I，那么B称为A的逆矩阵，记作A^(-1)。
逆矩阵可以通过伴随矩阵来求解。
对于n阶矩阵A，它的伴随矩阵记作adj(A)，它的定义为adj(A)=det(A)·A^(-1)，其中det(A)是A的行列式。
逆矩阵的求解可以通过以下步骤：（1）求解矩阵A的行列式det(A)；（2）求解矩阵A的伴随矩阵adj(A)；（3）求解矩阵A的逆矩阵A^(-1)，即A^(-1)=adj(A)/det(A)。
4. 特征值和特征向量特征值和特征向量在矩阵求解中起着重要作用。
设A 是一个n阶方阵，若存在一个非零向量X，满足AX=kX，其中k为常数，则k为A的一个特征值，X为对应的特征向量。
矩阵的求解方法和技巧
矩阵的求解方法和技巧矩阵的求解是线性代数中的一个重要问题，涉及到矩阵的性质、运算和解析方法等多个方面。
下面将介绍一些矩阵求解的常用方法和技巧。
1. 高斯消元法：高斯消元法是一种常用的线性方程组求解方法，适用于任意大小的方阵。
该方法的基本思想是通过矩阵的初等行变换，将方程组化为行最简的形式，从而求解出未知数的值。
具体操作步骤如下：1) 将方程组转化为增广矩阵形式；2) 选择一个主元（通常选择第一列的第一个非零元素）；3) 将该主元所在的行除以主元得到1；4) 用主元所在行乘以矩阵的某一行，再与原行相减，使得该行的主元所在列的其他元素都为0；5) 选择下一个主元，重复步骤3和4，直至将方程组化为行最简的形式（即上三角形矩阵）；6) 回代求解每个未知数的值。
2. 克拉默法则：克拉默法则适用于求解n元线性方程组（n个方程、n 个未知数），它是一种基于行列式的方法。
具体操作步骤如下：1) 将方程组转化为增广矩阵形式；2) 求出系数矩阵的行列式D；3) 分别将方程组的等号右边替换为未知数列矩阵，并求出每个矩阵列的行列式Dj；4) 利用克拉默法则的公式，未知数xi的值等于Dj除以D的商。
克拉默法则的优点是理论简单，适用于少数方程未知数的求解，但对于大规模的方程组来说，计算量较大。
3. LU分解法：LU分解是将矩阵按照一定的规则分解为一个下三角矩阵L和一个上三角矩阵U的乘积。
LU分解法适用于求解一大类线性方程组，对于已经进行了LU分解的矩阵，可以节省计算量，提高计算效率。
具体操作步骤如下：1) 对矩阵进行LU分解，得到下三角矩阵L和上三角矩阵U；2) 利用前代法（也称为Ly=b法）求解方程Ly=b，求出向量y；3) 利用回代法（也称为Ux=y法）求解方程Ux=y，求出向量x。
4. 矩阵的逆：矩阵的逆是指如果一个方阵存在逆矩阵，那么它和它的逆矩阵相乘得到一个单位矩阵。
矩阵的逆可以用来求解线性方程组的解。
具体操作步骤如下：1) 对矩阵A进行LU分解；2) 利用前代法求解方程Ly=b，求出向量y；3) 利用回代法求解方程Ux=y，求出向量x；4) 得到矩阵的逆矩阵A^-1。
常见的线性代数求解方法
常见的线性代数求解方法
1.列主元消去法
列主元消去法是一种经典的求解线性方程组的方法。
它通过将
方程组转化为上三角矩阵的形式来求解。
这个方法的关键在于选取
主元的策略。
一种常见的选取主元的策略是选择当前列中绝对值最
大的元素作为主元，然后进行消去操作，直到将矩阵转化为上三角
矩阵。
2.高斯-约当消去法
高斯-约当消去法是另一种常见的线性方程组求解方法。
它通
过消去矩阵的下三角部分来将线性方程组转化为上三角矩阵的形式。
这个方法也需要选择主元，常见的选择策略是选取当前行中绝对值
最大的元素作为主元，然后进行消去操作。
3.LU分解法
LU分解法是将矩阵分解为一对矩阵的乘积的方法。
这个方法的思想是先将矩阵分解为一个下三角矩阵和一个上三角矩阵，然后通过求解上三角矩阵和下三角矩阵的两个方程组来求解原始的线性方程组。
4.Jacobi迭代法
Jacobi迭代法是一种迭代求解线性方程组的方法。
它通过将原始的线性方程组转化为一个对角矩阵和另一个矩阵的乘积的形式，然后通过迭代求解这个对角矩阵和另一个矩阵的方程组来逼近线性方程组的解。
5.Gauss-Seidel迭代法
Gauss-Seidel迭代法是另一种迭代求解线性方程组的方法。
它与Jacobi迭代法类似，但是在每一次迭代中，它使用前一次迭代得到的部分解来更新当前的解。
这个方法通常比Jacobi迭代法收敛得更快。
以上是一些常见的线性代数求解方法。
每种方法都有其特点和适用范围，我们可以根据具体情况选择合适的方法来求解线性方程组的问题。
矩阵的线性方程组解法
矩阵的线性方程组解法线性方程组是数学中的重要概念，它描述了一组线性方程之间的关系。
而求解线性方程组的方法之一就是利用矩阵的运算进行计算。
本文将介绍几种常见的矩阵解法，以帮助读者更好地理解线性方程组求解的过程。
一、高斯消元法高斯消元法是求解线性方程组的基本方法之一。
它通过矩阵的行变换来简化系数矩阵，并最终将线性方程组化简为上三角形式。
步骤如下：1. 构建增广矩阵：将系数矩阵和常数向量合并成一个增广矩阵。
2. 初等行变换：利用加减乘除的运算，将增广矩阵化为上三角矩阵。
3. 回代求解：从方程组的最后一行开始，依次求解每个变量。
二、矩阵的逆解法对于非奇异矩阵（可逆矩阵），可以利用矩阵的逆求解线性方程组。
设线性方程组为Ax=b，其中A为系数矩阵，x为未知向量，b为常数向量。
解法如下：1. 判断A是否可逆：计算矩阵A的行列式，若不为零，则A可逆。
2. 计算逆矩阵：利用伴随矩阵法或初等变换法，求解A的逆矩阵A^-1。
3. 求解线性方程组：利用逆矩阵的性质，有 x=A^-1b。
三、克拉默法则克拉默法则是一种求解线性方程组的特殊方法，它通过计算行列式的比值来求解每个未知数的值。
步骤如下：1. 列出增广矩阵：将线性方程组化为增广矩阵形式。
2. 计算行列式：利用增广矩阵的系数部分，计算系数矩阵A的行列式det(A)。
3. 计算未知数：利用克拉默法则，有 xi=det(Ai)/det(A)，其中Ai是用b替换第i列得到的矩阵。
四、LU分解法LU分解法是一种将矩阵A分解为下三角矩阵L和上三角矩阵U的方法。
通过LU分解后，可以利用前代法和回代法求解线性方程组。
步骤如下：1. 进行LU分解：将系数矩阵A分解为下三角矩阵L和上三角矩阵U，有 A=LU。
2. 利用前代法求解Ly=b：先解 Ly=b 得到y的值。
3. 利用回代法求解Ux=y：再解 Ux=y 得到x的值。
总结：本文介绍了矩阵的线性方程组解法，包括高斯消元法、矩阵的逆解法、克拉默法则和LU分解法。
矩阵的线性方程组解集求解
矩阵的线性方程组解集求解线性方程组是线性代数中的重要概念，而解线性方程组就是求解方程组中未知数的解集。
在矩阵的线性方程组中，我们利用矩阵的运算和变换来求解线性方程组的解集。
本文将介绍矩阵的线性方程组求解的基本方法和步骤。
首先，我们来回顾一下线性方程组的定义：线性方程组是由多个线性方程组成的集合，其中每个方程都是线性的。
线性方程组的一般形式可以表示为：a1x1 + a2x2 + ... + anxn = b其中，a1, a2, ..., an 是系数，x1, x2, ..., xn 是未知数，b 是常数。
对于一个含有 m 个方程和 n 个未知数的线性方程组，可以使用矩阵的形式来表示：AX = B其中，A 是一个 m×n 矩阵，X 是一个 n×1 矩阵（列向量），B 是一个 m×1 矩阵（列向量）。
在这个形式下，我们的目标是求解 X 的取值。
下面，我们将介绍两种常见的矩阵的线性方程组求解方法：高斯消元法和矩阵的逆。
1. 高斯消元法高斯消元法是一种基本的矩阵求解方法，其基本思想是通过矩阵的初等行变换将线性方程组转化为上三角形式，从而求解未知数的值。
具体步骤如下：（1）将线性方程组的系数矩阵 A 与常数矩阵 B 合并为增广矩阵[A|B]。
（2）利用矩阵的初等行变换，将增广矩阵化为上三角形式。
（3）反向替换，从最后一行开始，求解每一个未知数的值。
（4）得到线性方程组的解集。
2. 矩阵的逆矩阵的逆是线性方程组求解的另一种方法。
对于方阵 A，如果存在一个方阵 B，使得 A×B = B×A = I，其中 I 是单位矩阵，则称矩阵 A 是可逆的，B 是 A 的逆矩阵。
利用矩阵的逆矩阵，我们可以通过以下方式求解线性方程组。
具体步骤如下：（1）对于矩阵 A，若 A 可逆，则将方程组 AX = B 两边同时左乘A 的逆矩阵 A^(-1)，得到 X = A^(-1)B。
（2）计算矩阵 A 的逆矩阵 A^(-1)。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
线性方程组的矩阵求解算法摘要 线性方程组的矩阵求解算法,只需在约当消元法的基础上,再对方程组的增广矩阵的行最简形进行行(列)删除和增加行,交换行等运算即可得到方程组的解,并且这种方法既可求解有唯一解的方程组.因而算法简单,易于实现.关键词 线性方程组;解向量;解法;约当消元法1 矩阵求解算法设有线性方程组m n A X b ⨯=,其增广矩阵())(1,m n A A b ⨯+=,算法的步骤如下: 第一步:利用约当消元法,把增广矩阵A 化为行最简形,设行最简形为()1m n B ⨯+.若()t i (),r A r =则方程组无解;否则设(),r A R =并执行以下步骤;第二步:删除B 中的所有零行和每一行第一个非零元素(这个非零元素一定是1)所在的列,得到矩阵()1,r n r D ⨯-+并记录每行的第一个非零元所在的列标,放在一维数组()1,,t r L 中,如第i 行的第一个非零元在第j 列,则()t i j =;第三步:构造矩阵()1m n r D H F ⨯-+⎛⎫= ⎪⎝⎭,其中 ()()1100001000010n r n r F -⨯-+-⎛⎫⎪-⎪= ⎪⎪-⎝⎭L L L L L L L L第四步:对矩阵H 中的行作交换运算:把H 中的第i 行(,1,1,i r r =-L 即从第r 行开始直到第一行)依次与其下一行交换,使之成为第()t i 行,交换运算结果后的矩阵记为G ,则G 中的前n r -个n 维列向量即为方程组的一个基础解系,最后一列向量即为方程组的一个特解; 第五步:写出方程组的通解.2 算法证明先证一个特殊情形,增广矩阵A 的行最简形矩阵B 的左上角为一r 阶的单位矩阵,即第i 行的第一个非零元的列标为i ,即()()1t i i i r =≤≤,所以设B 为1,1112,122,1100010001000r n r n r r rn r c c d c c d B c c d +++⎛⎫ ⎪ ⎪ ⎪=⎪ ⎪ ⎪ ⎪ ⎪⎝⎭L L L L LL L L L L L L L LL L L L L L L L LL 则1,1112,122,1r n r n r r rnr c c d c c d D c c d +++⎛⎫ ⎪ ⎪= ⎪ ⎪ ⎪⎝⎭L L L L L L L由上述算法可得H 为1,11,1112,12,222,1,2100001000010r r n r r n r r r r rn r c c c d c c c d c c c d H ++++++⎛⎫⎪ ⎪ ⎪⎪ ⎪=⎪-⎪- ⎪ ⎪ ⎪ ⎪-⎝⎭LL L L L L L L L L L L L L L L 由于()()1t i i i r =≤≤,故从H 得到G 时,H 中的行不需交换位置,即.G H = 那么矩阵B 的增广矩阵的线性方程组为111,111,222,112,1,.r r n n r r n nr r r r rn n x d c x c x x d c x c x x d c c x +++++=--⎧⎪=--⎪⎨⎪⎪=---⎩L L L L L L L L L L L令1,12,11,1100r r c c cr r α++⎛⎫ ⎪ ⎪ ⎪ ⎪+ ⎪= ⎪- ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭M M , 1,22,22,2011r r c c cr r α++⎛⎫ ⎪ ⎪ ⎪ ⎪+ ⎪= ⎪ ⎪- ⎪ ⎪ ⎪ ⎪-⎝⎭M M , ,L 12,001n n rn n r c c c α-⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪-⎝⎭M M 12000r d d d η⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭M M可以验证1,2,n r ααα-L 是方程组(1)所对应的齐次线性方程组的解,η是方程组(1)的特解,又12,,n r ααα-L 的后n r -个分量构成的向量组100-⎛⎫ ⎪ ⎪ ⎪ ⎪⎝⎭M ,001,,.01⎛⎫⎛⎫ ⎪ ⎪- ⎪ ⎪ ⎪ ⎪ ⎪ ⎪-⎝⎭⎝⎭LM M 线性无关,把它扩充成维向量组后也线性无关,所以12,,,n r ααα-L 线性无关,又因为()r A r =,所以方程组(1)的基础解系中有()n r A n r -=-个向量,因此12,,,n r ααα-L 即为方程组(1)的基础解系,特殊情形得证.对于行最简形矩阵B 为一般情形时,可以通过若干次列交换把它变形为上述特殊情形,但是,列交换将会导致最后结果中对应未知数的次序混乱,即在进行第i 列与第j 列的交换后,最后结果中i x 与j x 次序也就被交换了,因此,在这过程中,必须记住所进行的一切列交换,以便在最后结果中恢复,但若使用本矩阵求解算法,则可避免上述麻烦,为了叙述方便,还是只证一种特殊情形.设 121,2112,222,210000100010*******0r n r n r r rn r c c c d c c d B c c d +++⎛⎫ ⎪ ⎪ ⎪⎪= ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭L L L L L L L O L L L L L L L L L L L L L L L L L L L L LLL即()()()11,12,t t i i i r ==+≤≤则121,2112,222,200r n r n r r rnr c c c d c c d D c c d +++⎛⎫ ⎪⎪= ⎪ ⎪ ⎪⎝⎭L L L LL L L L,121,2112,222,200100001000010r n r n r r rn r c c c d c c d c c d H +++⎛⎫⎪ ⎪ ⎪⎪⎪= ⎪-⎪- ⎪ ⎪ ⎪ ⎪-⎝⎭LL L L L L L L L L L L L L L L , 121,2112,222,210000001000010r n r n r r rn r c c c d c c d G c c d +++⎛⎫ ⎪- ⎪ ⎪ ⎪⎪= ⎪ ⎪- ⎪ ⎪ ⎪ ⎪-⎝⎭LL L L LL L L L L L LL L L L 现在证明G 的前n r -个列向量是B 所对应的方程的基础解系,G 的最后一列是该方程组的特解,把矩阵B 的第2列依次与第3列,第4列,第列交换,得到矩阵'B121,2112,222',21000100001r n r n r r rnr c c c d c c d B c c d +++⎛⎫ ⎪⎪= ⎪ ⎪ ⎪⎝⎭L L L L L L L L L LL L L LL设矩阵所对应的方程组的解向量为12(,,,)T n x x x K ,'B 所对应的方程组的解向量为12(,,,)T n y y y K ,则有112132122,,,,,,,r r r r r n n x y x y x y x y x y x y ++++======L L即若1,2(,,)T n y y y L 是'B 所对应的方程组的解向量,则112,,2(,,,,,)T r r r n y y y y y y ++L L 是矩阵C 所对应的方程组的解向量,而由上述所证的特殊情形,'B 所对应的方程组的基础解系和一个特解分别为12'100,100c α⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪= ⎪- ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭M M 1,22,2,2'2010r r r r c c c α+++⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪= ⎪ ⎪- ⎪ ⎪ ⎪ ⎪⎝⎭M M , ,L 12'001n n rn n r c c c α-⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪-⎝⎭M M , 2000r d d d η⎛⎫⎪ ⎪ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭M M 由此可得矩阵所对应的方程组的基础解系和特解为110000c α⎛⎫ ⎪- ⎪ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭M M , 1,22,22,2010r r r r c c c α+++⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪= ⎪ ⎪- ⎪ ⎪ ⎪ ⎪⎝⎭M M , L , 12001n n n r rn c c c α-⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪-⎝⎭M M , 12000r d d d η⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭M M 而12,,,n r ααα-L ,η即为G 的列向量组,这一情形得证若B 为起它任意情形,只要重复上上述证明过程,即可得到证明.3 举例例 设有线性方程组12456712345671234567123456712345672322232612422436292551062411242x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x -+-+-=⎧⎪--+-++=⎪⎪-+-++-=⎨⎪--+-+-=⎪--+-+-=⎪⎩求其通解.解方程组的增广矩阵A 为121312212321611241121243629255106121411242A ---⎛⎫ ⎪--- ⎪ ⎪=--- ⎪---- ⎪ ⎪----⎝⎭A 的行最简形矩阵B 为12013122001381201391000104480000000000000B ---⎛⎫⎪-- ⎪ ⎪=--⎪ ⎪ ⎪ ⎪⎝⎭划掉B 中的最后两个零行和每行的第一个非零元所在的第一列,第三列,第四列,得矩阵D ,并且()()()11,23,34t t t ===2312208120139100448D ⎛⎫---⎪ ⎪=-⎪ ⎪--⎝⎭构造矩阵H231220812013910044810000010000010000010H ---⎛⎫⎪- ⎪ ⎪--⎪= ⎪-⎪- ⎪ ⎪- ⎪ ⎪-⎝⎭由于()34t =,所以应把中第3行依次与其后的行交换,使之成为第4行,然后因为()23t =,所以把H 中第2行依次与其后的行交换,使之成为第3行最后因()11t =,故第1行不需与任何行交换,这样变得到矩阵G ,23122100000812013910044810000010000010G ---⎛⎫⎪- ⎪ ⎪-⎪=-- ⎪⎪- ⎪ ⎪- ⎪ ⎪-⎝⎭所以方程组的通解为1234231220001028100913100484000100010000001k k k k ξ--⎛⎫⎛⎫⎛⎫-⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪- ⎪ ⎪ ⎪ ⎪⎪ ⎪ ⎪ ⎪- ⎪⎪ ⎪ ⎪ ⎪ ⎪ ⎪=+++-+- ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪- ⎪ ⎪ ⎪ ⎪⎪ ⎪ ⎪ ⎪ ⎪⎪- ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭-⎝⎭⎝⎭⎝⎭ 4.算法分析事实上,本算法是约当消元法的推广,因为若()()r A r A n ==时,最简形矩阵B的前n列为n阶单位矩阵,所以由B得D时,D为1n⨯矩阵,且为B的最后一列所构造成的矩阵,由D构造H时,不断增加行,由H得到G时,不需交换行,即G H D==,因而方程组的解向量为G,这也是约当消元法的结果也就是说约当()()==消元法是本r A r A n算法当时的特殊情形,由于本算法的所有加法和乘法都在把增广矩阵化为行最简形矩阵的着一过程中,所以有以下结论:1)算法的计算量与约当消元法的计算量相等;2)算法所需的存贮空间略多于约当消元法所需的存贮空间;3)在求方程组的通解时,其稳定性与精度和约当消元法的完全一致.另外,由于本算法从输入方程组到输出通解(或唯一解),中间的所有运算都是对矩阵进行的,所以算法简单,容易在计算机上实现,当然,由于本算法包含约当消元法,因而它除了有与约当消元法相同的缺点以外,它还有一个缺点:有时需要移动大量的元素,特别是当未知数的个数与方程的个数都很大时,元素的移动量可能更大.总之,本算法在约当消元法的基础上,不需增加乘法和加法运算,即可得到方程组的通解,因而本算法有一定的适用价值.参考文献[1]徐士良计算机常用算法[M] 北京: 清华大学出版社,1995.12[2]同济大学线性代数[M] 北京: 高等教育出版社, 2002.1[3]邓建中等计算方法[M] 西安: 西安交通大学出版社,2001.8[4]刘仲奎高等代数[M] 北京: 高等教育出版社,2003.8[5]浙江大学线性代数[M] 北京: 科学出版社,2001作者简介:王雪娇(1984.10.03) ,女,现就读于陇东学院数学系04级专科(1)班。