线性代数:矩阵消元

本节是网易公开课上的麻省理工大学线性代数课程第二节:矩阵消元 的学习笔记。

消元法

本节将讨论 消元法 来求解线性方程组。

线性方程组的系数矩阵如果是一个好的矩阵(good matrix)可以使用消元法来得到结果,否则无法得到结果。消元完成后,使用 回代 就可以得到方程组的解。

需要求解的方程组如下:

 x + 2y + z = 2
3x + 8y + z = 12
     4y + z = 2

方程组使用矩阵可以表示为 Ax = b

130284111xyz=2122

思考一个问题,方程一乘以多少,然后用方程二减去它能够将方程二的系数 x 消去。用矩阵语言描述即系数矩阵 A 的第一行乘以多少,然后用第二行减去它能够使得第二行第一列的元素为0。这里将系数矩阵 A 的第一个元素(1)称为主元,与方程中的 x 系数对应。 乘以的这个数称为消元乘数。

很简单,我们将第一行乘以3,然后用第二行减去它,就能够将 第二行第一列的元素变为0。

130284111row2=row23row1100224121

这里我们将矩阵的第二行第一列的元素((2,1))变为了0,所以使用 (2,1) 作为这一步的代号。下一步我们需要将第三行第一列变为0,代号为 (3,1)。很幸运的是,我们发现它已经是0了。以上的操作都是针对主元一来进行的,即保留第一个方程的系数 x ,将第二个和第三个方程的系数 x 都消去。

接下来我们要操作的是主元二,即保留第一个和第二个方程的系数 y ,将第三个方程的系数 y 都消去。对应到矩阵中的操作为将第三行第二列的元素变为0,代号为 (3,2)。操作过程为将第二行乘以多少,然后用第三行减去它,将第三行第二列元素变为0。答案很明显,消元乘数为2。

100224121row3=row32row2100
  • 4
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
### 回答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`函数进行消元,并打印消元后的结果。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值