球形空间产生器[高斯消元]

最新推荐文章于 2022-10-14 23:29:30 发布

原创最新推荐文章于 2022-10-14 23:29:30 发布 · 243 阅读

0 ·

CC 4.0 BY-SA版权

高斯消元专栏收录该内容

4 篇文章

订阅专栏

本文介绍了一种使用高斯消元法求解线性方程组的C++实现方法。通过交换行和消元操作，将矩阵转化为上三角形，从而求得方程组的解。代码中详细展示了如何读取系数和常数项，进行高斯消元，并输出结果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

传送门

每一行与第一行联立 , 然后高斯消元

#include<bits/stdc++.h>
#define N 15
using namespace std;
int n;
double c[N][N],f[N];
double gauss(){
	for(int i=1;i<=n;i++){
		for(int j=i;j<=n;j++){
			if(fabs(c[j][i]) > 1e-8){
				for(int k=1;k<=n+1;k++) swap(c[i][k],c[j][k]);
			}
		}
		for(int j=1;j<=n;j++){
			if(i==j) continue;
			double now = c[j][i] / c[i][i];
			for(int k=i;k<=n+1;k++) c[j][k] -= c[i][k] * now;
		}
	}
}
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%lf",&f[i]);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			double x;scanf("%lf",&x);
			c[i][j] = 2 * (x-f[j]);
			c[i][n+1] += x * x - (f[j] * f[j]);
		}
	}
	gauss();
	for(int i=1;i<=n;i++) printf("%0.3lf ",c[i][n+1]/c[i][i]);
	return 0;
}