高斯消元

高斯消元

  所谓高斯消元就是求解线性方程组的一种方法， 学过线性代数的人应该对该方法比较清楚，其实与我们求解多元一次方程的方法类似，而高斯消元则是采用矩阵消元的方式求解线性方程组的解。
  实现的思想很简单，但是步骤却比较繁琐，我们可以用一个二维矩阵来存储线性方程组中的系数，而用一维矩阵存储常数项，最后按照线性代数中将增广矩阵化简为阶梯型矩阵的方式，来求解出线性方程组的解。

例题

球形空间产生器

#include<bits/stdc++.h>
using namespace std;
int main(){
    int n;
    scanf("%d",&n);
    double p[n+1][n];
    double c[20][20];
    double b[20];
    for(int i=0;i<n+1;i++){
        for(int j=0;j<n;j++){
            scanf("%lf",&p[i][j]);
        }
    }
    //生成线性方程组矩阵
    for(int i=0;i<n;i++){
        b[i]=0.0;
        for(int j=0;j<n;j++){
            //将二次方程组转化为线性方程组
            c[i][j]=2*(p[i][j]-p[i+1][j]);
            b[i]+=p[i][j]*p[i][j]-p[i+1][j]*p[i+1][j];
        }
    }
    //高斯消元
    for(int i=0;i<n;i++){
        //找到主元不位0的一行，交换到第i行
        for(int j=i;j<n;j++){
            if(fabs(c[j][i])>1e-8){//找到主元不为0的一行交换到第i行，用来消去其他行
                for(int k=0;k<n;k++){
                    swap(c[i][k],c[j][k]);
                }
                swap(b[i],b[j]);
                break;
            }
        }
        //高斯消元
        //用其他行将去第i行的倍数，使得第i列除了第i行以外都是0
        for(int j=0;j<n;j++){
            if(i==j)continue;
            double rate=c[j][i]/c[i][i];
            for(int k=i;k<n;k++){
                c[j][k]-=c[i][k]*rate;
            }
            b[j]-=b[i]*rate;
        }
    }
    //得到结果
    for(int i=0;i<n;i++){
        printf("%.3f ",b[i]/c[i][i]);
    }
    return 0;
}

  我们先按照圆心到球上每一点的距离相等的性质得到n+1个2次方程，此时我们并不能直接求解方程组， 因为不是线性方程组，但是通过相邻两个点之间做差我们可以将二次方程转化为线性方程组，此时我们可以用二维数组模拟矩阵，再按照化简增广矩阵的方式求解线性方程组即可得到结果。