矩阵消元

矩阵消元法(Elimination)是常用的方程组Ax=b求解方法,但是该如何利用矩阵变换的思想来理解消元法呢?

假设有方程组


对于上述方程组,系数矩阵A为,整个消元过程为


 在用矩阵变换的思想理解矩阵消元之前,关于矩阵乘法的理解很重要,矩阵乘法可以从行的角度理解,列的角度理解,也可从矩阵单个元素的角度理解(这是我们通常的理解)。如果从行的角度,则相乘是矩阵行的线性组合,为了方便书写,这里给出的例子是向量与矩阵的相乘,因为矩阵-矩阵的相乘最终也是分解为矩阵-向量的相乘来实现的,如


这个过程表示进行了行变换,也就是说,要对一个矩阵进行变换,需在该矩阵左侧乘以一个矩阵。

如果从列的角度,则相乘是矩阵列的线性组合,如



这个过程表示进行了列变换,也就是说,要对一个矩阵进行列变换,需在该矩阵右侧乘以一个矩阵。

因此如果有矩阵A和矩阵B相乘A*B,我们可以理解为A对B进行了行变换,也可以理解为B对A进行了列变换。

 

回到矩阵的消元,有什么样的矩阵可以实现将A的第一行乘以-3加到第二行去,保持第一行不变?由于消元法进行的是行变换,因此上面的过程应该是在A的左侧乘以一个未知阵,即

设该未知消元矩阵E21为

则有

由上式易得E21=,E21表示将矩阵第2行第1列的数消去。

同样有什么样的矩阵可以实现将A的第二行乘以-2加到第三行去,保持第一、二行不变?即


按照上面的方法,可得消元矩阵E32=

也就是说上述的整个矩阵消元过程可以表示为 E32*(E21*A)=U

由于消元法中有可能涉及行之间的交换,因此还有一类重要的矩阵完成交换矩阵中两行次序的功能,称为置换矩阵(Permutation),同样按照上面矩阵相乘的思路,如果要交换行,则

如果要交换列,则

即如果要进行列交换,置换矩阵应放在右侧,如果进行行交换,置换矩阵放在左侧。

  • 5
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
### 回答1: C语言写矩阵消元代码的步骤如下: 1、定义一个二维数组,用于存储矩阵中的元素; 2、用for循环遍历数组,将最大元素置于对角线上; 3、用for循环实现除法运算,将其他元素除以对角线元素,形成单位化矩阵; 4、用for循环实现减法运算,将其他元素减去对角线元素,将其余元素置零; 5、用for循环实现乘法运算,将结果矩阵的元素乘以原矩阵的元素,得到最终结果。 ### 回答2: 矩阵消元是一种线性代数中常用的操作,可以将一个矩阵通过一系列行变换,得到一个上三角矩阵或者是行最简形矩阵。以下是用C语言编写的一个简单矩阵消元代码的示例: ```c #include <stdio.h> #define SIZE 3 // 矩阵的大小 void printMatrix(float matrix[SIZE][SIZE+1]) { for(int i=0; i<SIZE; i++) { for(int j=0; j<SIZE+1; j++) { printf("%6.2f ", matrix[i][j]); } printf("\n"); } } void gaussianElimination(float matrix[SIZE][SIZE+1]) { for(int i=0; i<SIZE-1; i++) { // 寻找主元素 int max_row = i; float max_val = matrix[i][i]; for(int j=i+1; j<SIZE; j++) { if(matrix[j][i] > max_val) { max_row = j; max_val = matrix[j][i]; } } // 交换行 if(max_row != i) { for(int j=i; j<SIZE+1; j++) { float temp = matrix[i][j]; matrix[i][j] = matrix[max_row][j]; matrix[max_row][j] = temp; } } // 消元 for(int j=i+1; j<SIZE; j++) { float factor = matrix[j][i] / matrix[i][i]; for(int k=i; k<SIZE+1; k++) { matrix[j][k] -= factor * matrix[i][k]; } } } } void backSubstitution(float matrix[SIZE][SIZE+1]) { float solution[SIZE]; for(int i=SIZE-1; i>=0; i--) { solution[i] = matrix[i][SIZE]; for(int j=i+1; j<SIZE; j++) { solution[i] -= matrix[i][j] * solution[j]; } solution[i] /= matrix[i][i]; } // 输出解 printf("解为:\n"); for(int i=0; i<SIZE; i++) { printf("x%d = %6.2f\n", i+1, solution[i]); } } int main() { float matrix[SIZE][SIZE+1] = { {2, 1, -1, -8}, {-3, -1, 2, 1}, {-2, 1, 2, -3} }; printf("原始矩阵:\n"); printMatrix(matrix); gaussianElimination(matrix); printf("上三角矩阵:\n"); printMatrix(matrix); backSubstitution(matrix); return 0; } ``` 这段代码实现了对一个大小为3的矩阵进行消元操作。在`main`函数中,我们首先定义了一个大小为3x4的二维数组`matrix`,其中前3列是矩阵的系数矩阵,最后1列是等号右侧的常数矩阵。 然后,通过调用`printMatrix`函数,我们输出了原始矩阵的内容。 接下来,我们调用`gaussianElimination`函数来进行高斯消元,将矩阵转化为上三角形式。在该函数中,我们采用了主元选取策略,确保每次消元时使用绝对值最大的元素作为主元素。然后通过行交换和消元操作,将主元下方的元素都变为0。 最后,我们调用`backSubstitution`函数对上三角矩阵进行回代操作,求解方程组的解,并输出结果。 运行该代码,输出结果将包括原始矩阵、上三角矩阵和方程组的解。这就是使用C语言编写的一个简单的矩阵消元的示例代码。 ### 回答3: 矩阵消元是线性代数中的一种常用操作,可以将一个矩阵转化为其标准形式,简化后续的计算操作。下面是一个用C语言写的简单矩阵消元的代码: ```c #include <stdio.h> #define MAX_SIZE 100 void printMatrix(int n, int mat[][MAX_SIZE]) { for(int i=0; i<n; i++) { for(int j=0; j<n; j++) { printf("%d ", mat[i][j]); } printf("\n"); } } void rowOperation(int n, int mat[][MAX_SIZE], int row1, int row2, int scalar) { for(int i=0; i<n; i++) { mat[row2][i] += scalar * mat[row1][i]; } } void rowEchelonForm(int n, int mat[][MAX_SIZE]) { int lead = 0; for(int r=0; r<n; r++) { if(lead >= n) return; int i = r; while(mat[i][lead] == 0) { i++; if(i == n) { i = r; lead++; if(lead == n) return; } } // 交换第r行和第i行 for(int j=0; j<n; j++) { int temp = mat[r][j]; mat[r][j] = mat[i][j]; mat[i][j] = temp; } // 将第r行的首元素变为1 int scalar = mat[r][lead]; for(int j=0; j<n; j++) { mat[r][j] /= scalar; } // 消去其他行中的首元素 for(int j=0; j<n; j++) { if(j != r) { scalar = mat[j][lead]; rowOperation(n, mat, r, j, -scalar); } } lead++; } } int main() { int n; printf("请输入矩阵的维度:"); scanf("%d", &n); int matrix[MAX_SIZE][MAX_SIZE]; printf("请输入矩阵的元素:\n"); for(int i=0; i<n; i++) { for(int j=0; j<n; j++) { scanf("%d", &matrix[i][j]); } } printf("原始矩阵为:\n"); printMatrix(n, matrix); rowEchelonForm(n, matrix); printf("矩阵消元后的结果为:\n"); printMatrix(n, matrix); return 0; } ``` 以上代码实现了一个求矩阵消元的函数`rowEchelonForm`,可以将输入的矩阵转化为行阶梯形式。在`main`函数中,首先获取矩阵的维度,然后输入矩阵的元素,调用`rowEchelonForm`函数进行消元,并打印消元后的结果。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值