高斯消元法解方程Ax=b的C源代码实现

最新推荐文章于 2023-03-23 18:40:52 发布

silent_rye

最新推荐文章于 2023-03-23 18:40:52 发布

阅读量1k

点赞数 1

文章标签： c语言线性方程组

本文链接：https://blog.csdn.net/silent_rye/article/details/50483910

版权

这篇博客介绍了如何用C语言实现高斯消元法来解决线性方程组的问题。代码中采用了列选主元策略，并通过两个辅助数组跟踪行变换状态和顺序。在回代过程中，利用这些信息逆序计算最终解。用户需输入方阵的维度n，随后依次输入A矩阵的n行，最后输入向量b。

摘要由CSDN通过智能技术生成

采用列选主元方式实现，过程中维护了两个数组，分别记录改行是否已经选取主元，与选取主元行数的顺序。回代时利用选取顺序逆序得到最终结果。输入时先输入矩阵A的维数（方阵）n，再分别输入n个行向量；最后输入b。代码如下（由于使用了_Bool类型数组，需添加编译选项-std=c99，或者你可以直接改成int类型的：

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
double Abs(double x) { return (x > 0) ? x : -x; }
int* GaussE(double **A, double *b, int size);
double* backSub(double **A, double *b, int *order, int size);
int findMaxLin(double **A, _Bool *flag, int col, int size);
int main(void) {
    int size;
    double **A;
    double *b;
    double *slt;
    int *order;
    int cnt;
    scanf("%d", &size);
    getchar();
    A = (double **)malloc(sizeof</