Cramer法则

最新推荐文章于 2023-11-20 10:57:03 发布

贝桑不止学Python

最新推荐文章于 2023-11-20 10:57:03 发布

阅读量218

点赞数

分类专栏：数值分析学习文章标签：算法数据结构 c语言

本文链接：https://blog.csdn.net/weixin_66030644/article/details/129903030

版权

数值分析学习专栏收录该内容

17 篇文章 0 订阅

订阅专栏

该代码实现了CramersRule算法来解n维线性方程组，通过计算矩阵的行列式和替换矩阵的列来找出解。程序首先读取方程组的系数和常数项，然后计算行列式，如果行列式不为零，则使用Cramer公式求解并输出结果；否则，输出Fail.。

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

#include <stdio.h>
#include <math.h>
#include <stdlib.h>

double determinant(int n, double A[][25]);						
int Cramer(int n, double A[][25], double b[], double x[]);		


int main(){
	
	int n;					
	scanf("%d", &n);
	
	double A[25][25], b[25], x[25];

	for(int i = 0; i < n; i++){			
		for(int j = 0; j < n; j++){
			scanf("%lf", &A[i][j]);
		}
	}
	for(int i = 0;i < n; i++){
		scanf("%lf", &b[i]);
	}
	
  	int flag = Cramer(n, A, b, x);				
  
	if(flag){								
		for(int i = 0;	 i < n; i++){
			printf("%lf ", x[i]);
		}
	}
	else{									
		printf("Fail.");		
	}

} 

double determinant(int n, double A[][25]){
	double ret = 0;
	double B[25][25];
	if (n == 1) {
		ret = A[0][0];
	} else {
		for(int i=0; i<n; i++) {
			for (int j=0; j<i; j ++) {
				for (int k=1; k<n; k ++) {
					B[j][k-1] = A[j][k]; 
				}
			}
			for (int j=i; j<n; j ++) {
				for (int k=1; k<n; k ++) {
					B[j][k-1] = A[j+1][k]; 
				}
			}
			ret += A[i][0] * pow(-1, i+0) * determinant(n-1, B);
		}
	}
	return ret;
	
}

int Cramer(int n, double A[][25], double *b, double *x){
	
	double D = determinant(n, A);			
	double As[25][25];					
	
	if(D > 1e-15){							
		
		for (int i=0; i<n; i++) {
			for (int j=0; j<n; j++) {
				As[i][j] = A[i][j];
			}
		}					
		
		for(int k=0; k<n; k++){		
			for (int l=0; l<n; l++) {
				As[l][k] = b[l];
			}						
			x[k] = determinant(n, As) / D;
			for (int i=0; i<n; i++) {
				for (int j=0; j<n; j++) {
					As[i][j] = A[i][j];
				}
			}
		}
		
		return 1;						
	
    }
	
  	return 0;			

}