c语言计算线性方程

已于 2024-08-14 17:36:25 修改
阅读量293 收藏 2
点赞数 13
文章标签： c语言 开发语言
于 2024-08-14 17:36:04 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/cuitlse/article/details/141196274
版权

手动解线性方程组有两种方法。

1.通过消元求解

2.通过矩阵行列式等来求解

我们现在要用c语言来实现用行列式求解线性方程组。在求解之前我们首先要学习行列式相关知识，如果不了解行列式请先去百度行列式进行学习，记得回来哦~

可以发现行列式的数字排法很像二维数组的排法，于是我们用二维数组保存行列式内容。所以我们只需要把输入的数据存入二维数组内。

因为数组的大小在定义的时候必须是确定的，而我们可能要算二元一次方程组，我们有四个数据，也可能是三元的我们有九个数据，于是二元的需要定义一个a[2][2]，三元的我们要定义一个a[3][3]。这不是一个确定的数，所以我们要才用一种特殊的方式(动态内存分配)。这样就能产生一个动态的数组。

代码如下：

//从输入的数据中存到数组内
double* getArrFromDet(int n)
{
	double num;
	double* p = (double*)calloc(1, sizeof(double) * n * n);
	if (!p)return -1;
	int i = 0;
	while (n * n != i)
	{
		scanf("%lf", &num);
		*(p + i) = num;
		i++;
	}
	return p;
}

现在我们剩下的问题就是计算数组里面的行列式的值。

我们学过计算行列式有一种方法，叫做拉普拉斯定理。拉普拉斯定理，计算降阶行列式的一种方法。该定理断言：在n阶行列式D=|aij| 中，任意取定k行（列），1≤k≤n-1，由这k行（列）的元素所构成的一切k阶子式与其代数余子式的乘积的和等于行列式D的值。此展式称为拉普拉斯展式。

于是我们就运用这个定理，我们运用第一行的各元素乘以对于的各元素的代数余子式。代数余子式又是一个行列式继续用这个方法，一直到最后只剩下一个就返还这个的值。这像不像是个递归，那我们就用递归来解决这个问题吧

哦~对了好像没告诉你们代数余子式怎么求。现在我们来讲一讲他的求法。

只需要把选定的这个元素的行，和列全部去除，组成的行列式就是余子式，然后代数余子式就是有±号，比如说在第1行第1列的元素的代数余子式就是1×余子式，而第1行第2列的元素就是-1×余子式，是不是能发现什么规律，就是(-1)∧(行＋列)×行列式就可以了。好了这个懂了

就来看代码吧。

递归结束条件，只剩下一个元素在行列式中。

//将数组进行行列计算
double calDet(double *arr,int n)
{
	int i = 0;
    double sum = 0;
	if (n == 1)
	{
		return *arr;
	}
	else
	{
		while (i < n)
		{
			double* arrTemp = (double*)calloc(1, sizeof(double) * (n - 1) * (n - 1));
			double* temp = arrTemp;
			//去除正在递归的该行列的
			for (int k =1; k < n; k++)
			{
				for (int j = 0; j < n; j++)
				{
					if (j != i)
					{
						*temp = *((arr + k * n) + j);
						temp++;
					}
				}
			}
			sum += pow(-1,i)*(*(arr + i) * calDet(arrTemp, n - 1));//递归调用
			i++;
 
		}
	}
	return sum;
}

现在我们会求行列式了，那问题只剩下一个了：计算多个行列式的值并保存下来然后再运用克拉默法则运算出方程组的解。具体代码如下：

//求解线性方程组
void solveEqu(int n)
{
	double value[127] = { 0 };
	//获取d的值
	printf("d=");
	value[0] = calDet(getArrFromDet(n), n);
	int i = 1;
	while (i != n + 1)
	{
		printf("d%d=", i);
		value[i] = calDet(getArrFromDet(n), n);
		i++;
	}
	//输出方程的解
	i = 1;
	while (i != n + 1)
	{
		if (value[0] == 0 && value[i] != 0)
		{
			printf("方程无解");
			break;
		}
		else
		{
 
			printf("x%d = %lf\n", i, value[i] / value[0]);
			i++;
		}
	}
}

这样我们就能解出一个线性方程组。

接下来让我们看看这个简单的主函数吧

int main()
{
	int n;
	printf("输入方程组的元数");
	scanf("%d", &n);
    solveEqu(n);
	while (1);
	return 0;
}

 对了,再加上主函数

#include<stdio.h>

在用用高斯消元解法

完整代码如下：

#include <stdio.h>
#include <math.h>
#define MAX 20

void swap(double* a,double*b){
	double tmp = *a;
	*a = *b;
	*b = tmp;
}

void swapArray(double a[MAX], double b[MAX], int n){
	int i;
	for(i=0;i<n;i++)
		swap(&a[i],&b[i]);	
}

double computeDet(double array[MAX][MAX], int n){
    // 通过初等行变换，再计算对角元素乘积，来计算行列式的值

    int i,j,k;
    int sign=0;             //行列式交换一次需要改变符号，此变量记录交换次数 
    double sum=1.0;         //结果
    double tmp;             //暂存乘积因子
    double zero=1e-6;       //两浮点数差距小于1e-6视为相等

    for(i=0;i<n-1;i++){
        k=1;
        while(fabs(array[i][i]-0.0)<zero && i+k<n){
            swapArray(array[i],array[i+(k++)],n);
			sign++;
        }

        for(j=i+1;j<n;j++){
			if(array[j][i]==0.0)			//如为0则那一行不用化简 
				continue;
			else{
				tmp = -(double)array[j][i]/array[i][i];		//保存乘积因子 
				for(k=i;k<n;k++)
					array[j][k] += (tmp*array[i][k]);
			}
		}
    }

    for(i=0;i<n;i++)
		sum *= array[i][i];
	if(sign%2!=0)				//交换偶数次符仍为正 
		sum *= -1;
    return sum;
}

void rowTrans(double array[MAX][MAX], int n){
    // 做初等行变换
    int i,j,k;
    double tmp;             //暂存乘积因子
    double zero=1e-6;       //两浮点数差距小于1e-6视为相等

    for(i=0;i<n-1;i++){
        k=1;
        while(fabs(array[i][i]-0.0)<zero && i+k<n){
            swapArray(array[i],array[i+(k++)],n);
        }

        for(j=i+1;j<n;j++){
			if(array[j][i]==0.0)			//如为0则那一行不用化简 
				continue;
			else{
				tmp = -(double)array[j][i]/array[i][i];		//保存乘积因子 
				for(k=i;k<n;k++)
					array[j][k] += (tmp*array[i][k]);
			}
		}
    }

}

void rowTrans_Ab(double array[MAX][MAX], int n){
    // 做A、b同时做初等行变换
    int i,j,k;
    double tmp;             //暂存乘积因子
    double zero=1e-6;       //两浮点数差距小于1e-6视为相等

    for(i=0;i<n;i++){
        k=1;
        while(fabs(array[i][i]-0.0)<zero && i+k<n){
            swapArray(array[i],array[i+(k++)],n+1);
        }

        for(j=i+1;j<n;j++){
			if(array[j][i]==0.0)			//如为0则那一行不用化简 
				continue;
			else{
				tmp = -(double)array[j][i]/array[i][i];		//保存乘积因子 
				for(k=i;k<=n;k++)
					array[j][k] += (tmp*array[i][k]);
			}
		}
    }

}

void printMatrix(double array[MAX][MAX], int n){
    int i,j;
    for(i=0;i<n;i++){
        for(j=0;j<n;j++){
            printf("%lf",array[i][j]);
            if(j==n-1){
                printf("\n");
            }else{
                printf("\t");
            }
        }
    }
}

void printArray(double x[MAX], int n){
    int i;
    for(i=0;i<n;i++){
        printf("%lf",x[i]);
        if(i==n-1){
            printf("\n");
        }else{
            printf("\t");
        }
    }
}

void copyMatrix(double a[MAX][MAX], double b[MAX][MAX], int n){
    int i,j;
    for(i=0;i<n;i++){
        for(j=0;j<n;i++){
            b[i][j]=a[i][j];
        }
    }
}

void copyArray(double a[MAX], double b[MAX], int n){
    int i;
    for(i=0;i<n;i++){
        b[i]=a[i];
    }
}

void solveLowerEquations(double A[MAX][MAX], double x[MAX], double b[MAX], int n){
    //求解系数矩阵为下三角矩阵的线性方程组
    int i,j;
    copyArray(b,x,n);
    for(i=0;i<n;i++){
        for(j=0;j<i;j++){
            x[i]-=A[i][j]*x[j];
        }
        x[i]/=A[i][i];
    }
}

void solveUpperEquations(double A[MAX][MAX], double x[MAX], double b[MAX], int n){
    //求解系数矩阵为上三角矩阵的线性方程组
    int i,j;
    copyArray(b,x,n);
    for(i=n-1;i>=0;i--){
        for(j=i+1;j<n;j++){
            x[i]-=A[i][j]*x[j];
        }
        x[i]/=A[i][i];
    }
}

//解一般方程
int main(){
    int i,j;
    int n;				//阶数 
    double tmp;         //临时变量
    double x[MAX];      //未知数矩阵
    double b[MAX];      //常系数矩阵
    double A[MAX][MAX]; //系数矩阵


    //数据初始化
    printf("Please enter the dimensions of the coefficient matrix:(n)");
	scanf("%d",&n);

    printf("Please enter the value of coefficient matrix:(A)\n");
	for(i=0;i<n;i++){
		for(j=0;j<n;j++)
			scanf("%lf",&A[i][j]);
	}

    printf("Please enter the value of the constant term matrix:(b)\n");
    for(i=0;i<n;i++){
        scanf("%lf",&b[i]);
    }

	//初始化x和Ab
    for(i=0;i<n;i++){
        x[i]=0.0;
        A[i][n]=b[i];
    }

    //对矩阵A做初等行变换
    rowTrans_Ab(A,n);
    
    printf("rowTrans:\n");
	printMatrix(A,n);
    
    //取变化后的b
    for(i=0;i<n;i++){
        b[i]=A[i][n];
    }

    //解上三角矩阵
    solveUpperEquations(A,x,b,n);

    //输出结果
    printf("The result are:\n");
    printArray(x,n);
    return 0;
}

还有Cholesky法求正定对称矩阵的线性方程组

完整代码如下：

#include <stdio.h>
#include <math.h>
#define MAX 20

void swap(double* a,double*b){
	double tmp = *a;
	*a = *b;
	*b = tmp;
}

void swapArray(double a[MAX], double b[MAX], int n){
	int i;
	for(i=0;i<n;i++)
		swap(&a[i],&b[i]);	
}

double computeDet(double array[MAX][MAX], int n){
    // 通过初等行变换，再计算对角元素乘积，来计算行列式的值

    int i,j,k;
    int sign=0;             //行列式交换一次需要改变符号，此变量记录交换次数 
    double sum=1.0;         //结果
    double tmp;             //暂存乘积因子
    double zero=1e-6;       //两浮点数差距小于1e-6视为相等

    for(i=0;i<n-1;i++){
        k=1;
        while(fabs(array[i][i]-0.0)<zero && i+k<n){
            swapArray(array[i],array[i+(k++)],n);
			sign++;
        }

        for(j=i+1;j<n;j++){
			if(array[j][i]==0.0)			//如为0则那一行不用化简 
				continue;
			else{
				tmp = -(double)array[j][i]/array[i][i];		//保存乘积因子 
				for(k=i;k<n;k++)
					array[j][k] += (tmp*array[i][k]);
			}
		}
    }

    for(i=0;i<n;i++)
		sum *= array[i][i];
	if(sign%2!=0)				//交换偶数次符仍为正 
		sum *= -1;
    return sum;

}

void rowTrans(double array[MAX][MAX], int n){
    // 做初等行变换
    int i,j,k;
    double tmp;             //暂存乘积因子
    double zero=1e-6;       //两浮点数差距小于1e-6视为相等

    for(i=0;i<n-1;i++){
        k=1;
        while(fabs(array[i][i]-0.0)<zero && i+k<n){
            swapArray(array[i],array[i+(k++)],n);
        }

        for(j=i+1;j<n;j++){
			if(array[j][i]==0.0)			//如为0则那一行不用化简 
				continue;
			else{
				tmp = -(double)array[j][i]/array[i][i];		//保存乘积因子 
				for(k=i;k<n;k++)
					array[j][k] += (tmp*array[i][k]);
			}
		}
    }

}

void rowTrans_Ab(double array[MAX][MAX], int n){
    // 做A、b同时做初等行变换
    int i,j,k;
    double tmp;             //暂存乘积因子
    double zero=1e-6;       //两浮点数差距小于1e-6视为相等

    for(i=0;i<n;i++){
        k=1;
        while(fabs(array[i][i]-0.0)<zero && i+k<n){
            swapArray(array[i],array[i+(k++)],n+1);
        }

        for(j=i+1;j<n;j++){
			if(array[j][i]==0.0)			//如为0则那一行不用化简 
				continue;
			else{
				tmp = -(double)array[j][i]/array[i][i];		//保存乘积因子 
				for(k=i;k<=n;k++)
					array[j][k] += (tmp*array[i][k]);
			}
		}
    }

}

void printMatrix(double array[MAX][MAX], int n){
    int i,j;
    for(i=0;i<n;i++){
        for(j=0;j<n;j++){
            printf("%lf",array[i][j]);
            if(j==n-1){
                printf("\n");
            }else{
                printf("\t");
            }
        }
    }
}

void printArray(double x[MAX], int n){
    int i;
    for(i=0;i<n;i++){
        printf("%lf",x[i]);
        if(i==n-1){
            printf("\n");
        }else{
            printf("\t");
        }
    }
}

void copyMatrix(double a[MAX][MAX], double b[MAX][MAX], int n){
    int i,j;
    for(i=0;i<n;i++){
        for(j=0;j<n;i++){
            b[i][j]=a[i][j];
        }
    }
}

void copyArray(double a[MAX], double b[MAX], int n){
    int i;
    for(i=0;i<n;i++){
        b[i]=a[i];
    }
}

void transpose(double a[MAX][MAX], int n){
    int i,j;
    for(i=0;i<n-1;i++){
        for(j=i+1;j<n;j++){
            swap(&a[i][j], &a[j][i]);
        }
    }
}

void solveLowerEquations(double A[MAX][MAX], double x[MAX], double b[MAX], int n){
    //求解系数矩阵为下三角矩阵的线性方程组
    int i,j;
    copyArray(b,x,n);
    for(i=0;i<n;i++){
        for(j=0;j<i;j++){
            x[i]-=A[i][j]*x[j];
        }
        x[i]/=A[i][i];
    }
}

void solveUpperEquations(double A[MAX][MAX], double x[MAX], double b[MAX], int n){
    //求解系数矩阵为上三角矩阵的线性方程组
    int i,j;
    copyArray(b,x,n);
    for(i=n-1;i>=0;i--){
        for(j=i+1;j<n;j++){
            x[i]-=A[i][j]*x[j];
        }
        x[i]/=A[i][i];
    }
}


int main(){
    int i,j;
    int n;				//阶数 
    double tmp;         //临时变量
    double x[MAX];      //未知数矩阵
    double A[MAX][MAX];  //系数矩阵
    double b[MAX];       //常系数矩阵
	
    //数据初始化
    printf("Please enter the dimensions of the coefficient matrix:(n)");
	scanf("%d",&n);

    printf("Please enter the value of coefficient matrix:(A)\n");
	for(i=0;i<n;i++){
		for(j=0;j<n;j++)
			scanf("%lf",&A[i][j]);
	}

    printf("Please enter the value of the constant term matrix:(b)\n");
    for(i=0;i<n;i++){
        scanf("%lf",&b[i]);
    }

    for(i=0;i<n;i++){
        x[i]=0.0;
    }

    //对矩阵A做初等行变换
    rowTrans(A,n);
    printf("rowTrans:\n");
	printMatrix(A,n);

    //处理主元，得到L.T(上三角矩阵)
    for(i=0;i<n;i++){
        if(A[i][i]<0){
            for(j=i;j<n;j++){
                A[i][j]*=-1;
            }
        }
        if(A[i][i]!=1){
            tmp=sqrt(A[i][i]);
            for(j=i;j<n;j++){
                A[i][j]/=tmp;
            }
        }
        
    }
    printf("L.T:\n");
    printMatrix(A,n);

    //求得L（下三角矩阵），求解第一步方程Ly=b，求得y
    transpose(A,n);
    printf("L:\n");
    printMatrix(A,n);

    solveLowerEquations(A,x,b,n);
    printf("y:\n");
    printArray(x,n);

    //求得L.T（上三角矩阵），求解第二步方程L.T*x=y，求得x
    transpose(A,n);
    copyArray(x,b,n);//x中暂存的y赋给b
    solveUpperEquations(A,x,b,n);

    printf("The result are:\n");
    printArray(x,n);
    return 0;
    //期望结果:(1,-1,0,2,1,-1,0,2)
}

cuitlse
关注
  • 13
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
C语言计算线性方程组
m0_57614828的博客
03-10 7229
今天让我们来用c解线性方程组。 手动解线性方程组有两种方法。 1.通过消元求解 2.通过矩阵行列式等来求解 我们现在要用c语言来实现用行列式求解线性方程组。在求解之前我们首先要学习行列式相关知识，如果不了解行列式请先去学习行列式 ...
C语言解线性方程组
qq_44199379的博客
11-22 5832
解线性方程组一般解法的C语言实现
解线性方程组 c语言,C语言解线性方程组
weixin_39828457的博客
05-19 795
《C语言解线性方程组》由会员分享，可在线阅读，更多相关《C语言解线性方程组(6页珍藏版)》请在人人文库网上搜索。1、实验一、用 C 语言解线性方程组一、实验名称 ：用 C 语言解线性方程组。二、实验原理 ：用初等行变化的方法将线性方程组的增广矩阵化简成最简型； 最简型中可以看出线性方程组解的情况。三、实验工具 ：计算机和 Microsoft Visual C+ 6.0 软件。四、实验原程序代码如下...
C语言一元线性回归方程
qq_39901989的博客
11-22 801
C语言一元线性回归方程
用C语言实现线性方程组的数值解法
qq_43479228的博客
09-09 3900
简介： 数值求解方法主要有以下两条途径： 直接法：利用Gauss消元或矩阵分解，通过有限次运算可求出精确解。 迭代法：构造迭代格式，产生迭代序列，通过无限次迭代过程求解。有限次截断得近似解。 其中，直接解法可采用Gauss消元法、Doolittle分解法（三角分解法、LU分解法）、Crout分解法（与Doolittle类似）、Cholesky分解法（平方根法）。 标注： 此笔记中的四种消元方法已用...
c语言编程非线性方程求解,c语言计算机编程三种方法求解非线性方程
weixin_36076907的博客
05-20 1301
c语言计算机编程三种方法求解非线性方程 本 科 专 业 学 年 论 文题 目：非线性方程求解比较姓 名： 何 娟 专 业： 计算机科学技术系 班 级： 08 级本科(2)班 指 导 老 师： 刘 晓 娜 完成日期： 2010 年 11 月 21 日题 目：非线性方程求解比较摘 要本文给出了三种求解非线性方程的方法，分别是二分法，牛顿迭代法，割弦法。二分法巧妙地利用插值得到的点以及有根区间中点这两点...
c语言非线性方程二分法,二分法求非线性方程根（C语言）
weixin_39756696的博客
05-19 726
this.p={ m:2,b:2,loftPermalink:'',id:'fks_087070085095085075085087083095085086087068093080080065',blogTitle:'二分法求非线性方程根(C语言)',blogAbstract:'二分法，C语言实现',blogTag:'数值分析',blogUrl:'blog/static/106796622007...
c语言非线性方程根的求解,实验三_非线性方程求根
weixin_39640265的博客
05-22 819
实验三 非线性方程求根专业班级： 姓名： 学号：一、实验目的熟悉C语言编程；学习非线性方程求根的方法及程序设计算法二、实验题目迭代函数对收敛性的影响用迭代法求方程的根f(x) 2x3 0x 1方案1：化方程为等价方程 x ( x ) 2取初值x0 0 ，迭代10次。3x 2x 1 (x) 方案2：化 f ( x ) 0 ...
C语言实现解线性方程组AX=b
4-5月更新计算机程序设计的一些实例，之后再更新计算机基础，机器智能相关的一些内容。
04-10 568
/解线性方程组原理是高斯消元，实际上只要用每列第一个不为0的数将所在列其他不为0数都化为0，int temp, temp1;//记录非零元素所在行数。//最后第k列不为0的数对应bk即为解xk。#define N 3 //5*5为例。//暂不考虑无解和无穷多解求基础解系情形。//解决n*n，同时有唯一解情况。
C语言二分法计算非线性方程的解
程序必须为人类阅读而编写，并且仅仅碰巧可以让机器执行
01-11 1140
文章目录输入输出代码 输入 输入为给定的 x 的范围以及要求精度，例如求解非线性方程 x3−x−1=0x^3 - x - 1 = 0x3−x−1=0 在 [1,2][1, 2][1,2] 区间内的解，其中精度取 10−510^{-5}10−5。 注：这里初始给的区间范围不能太宽，并且应当保证方程的根就在区间内。 输出 方程在给定精度下的近似解 1.32472。 代码 #include <iostream> using namespace std; double function(doub
C语言解线性方程的四种方法
09-04
在C语言中，解线性方程通常涉及矩阵运算，这对于理解和编程有一定的挑战。...总之，解线性方程在C语言中的实现涉及矩阵运算和数值计算，需要理解各种方法的原理，并关注数值稳定性、错误处理和计算效率。
C语言解线性方程组C语言解线性方程组.doc
07-02
C语言解线性方程组C语言解线性方程组
C语言解线性方程的四种方法共7页.pdf.zip
10-29
在C语言中，解线性方程是数值分析的一个重要组成部分。线性方程通常以矩阵形式表示，即Ax=b，其中A是一个系数矩阵，x是未知数向量，b是常数向量。本资料“C语言解线性方程的四种方法共7页”可能涵盖了以下几种常见的...
C语言程序设计实现解非齐次线性方程
04-13
采用阶梯矩阵找出非齐次线性方程组的增广矩阵的秩，用大小为未知量个数的双向栈存栈储自由未知量与非自由未知量，并给出在微机上运行的模拟人工解题的C语言计算程序。 1.2 要求 明确课程设计的目的，能根据课程设计...
计算方法实验二非线性方程求根实验报告.docx
05-30
非线性方程求根是计算方法中的一个重要课题，它涉及到数学、计算机科学以及工程等多个领域。本实验报告主要探讨了两种求解非线性方程根的方法：二分法和Newton迭代法，并通过C语言编程实现这两种方法，旨在使学生...
初学嵌入式-C语言常犯错误详解
2301_78622258的博客
08-09 1624
C语言常见错误详解
图数据库框架及其支持的开发语言和应用场景
极客栈
08-12 379
这些图数据库框架通过提供多种编程语言的支持，使得开发者能够灵活选择适合自己项目的语言进行开发。同时，社区和第三方开发者也会提供各种语言的库和工具，进一步丰富了图数据库的开发生态。
财务收支管理系统(自动统计，智能查询，).xlsm
08-13
工资表，财务报表，对账表，付款申请，财务报告，费用支出表，财务收支 适用人群：学习不同技术领域的小白或进阶学习者；可作为毕设项目、课程设计、大作业、工程实训或初期项目立项。
核电站建设和运营 头豹词条报告系列.pdf
最新发布
08-13
电子元件 电子行业 行业分析 数据分析 数据报告 行业报告
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值