最近在学数值分析，正好学到求解线性方程组。就自己动手简单实现了一下。关于本算法的原理可以在《数值分析》第5版(李庆扬编)，对应于该书的P145页，详细讲解了公式。因本人时间有限，暂时不详细编辑公式，等空闲了再来重新补充。
简要的说明下该算法的应用吧，高斯消元法在线性代数那门课肯定学过了，对于一般简单的3、4阶线性方程组，还可以进行纸上的笔运算，但是当数目过多，算起来就比较吃力了，所以借助于计算机实现。如果很多的数学软件里都有矩阵的求解的函数，如mathematic和matlab，因本人对这些软件不熟悉，所以也没去去尝试使用。
对于非奇异矩阵并且主对角线上的元素不为0。就可以用高斯消元法来求解。
我直接用简单的C++编程方式，有点类C风格，并没有运用C++的奇淫技巧，由于本人技术不成熟，时间有限，就暂时用简单的实现。等空闲了，一定把奇淫技巧运用上去。
先贴代码吧！
#include
using namespace std;
const int n = 3;
void gaussin(double a[n][n], double b[n])
{
//判断能否用高斯消元法，如果矩阵主对角线上有0元素存在是不能用的
for (int i = 0; i < n; i++)
if (a[i][i] == 0)
{
cout << "can't use gaussin meathod" << endl;
return;
}
int i, j, k;
double c[n]; //存储初等行变换的系数，用于行的相减
//消元的整个过程如下，总共n-1次消元过程。
for (k = 0; k < n - 1; k++)
{
//求出第K次初等行变换的系数
for (i = k + 1; i < n; i++)
c[i] = a[i][k] / a[k][k];
//第K次的消元计算
for (i = k + 1; i < n; i++)
{
for (j = 0; j < n; j++)
{
a[i][j] = a[i][j] - c[i] * a[k][j];
}
b[i] = b[i] - c[i] * b[k];
}
}
//解的存储数组
double x[n];
//先计算出最后一个未知数；
x[n - 1] = b[n - 1] / a[n - 1][n - 1];
//求出每个未知数的值
for (i = n - 2; i >= 0; i--)
{
double sum = 0;
for (j = i + 1; j < n; j++)
{
sum += a[i][j] * x[j];
}
x[i] = (b[i] - sum) / a[i][i];
}
cout << " the solution of the equations is:" << endl;
cout << endl;
for (i = 0; i < n; i++)
cout <
}
1. 我们进行判断主对角线元素是否有0值，如果有0值就不能用这个方法了。因为这个简单的高斯消元法并没有进行行列交换，所以主队角线不能为0，否则那个行变换的大小就不能得到(除数为0.undefined).2.为初等行变换计算变换大小，存储在数组C[N] 3.进行消元过程。 4.外循环继续第二次重复(2)(3)步骤，直到K=N-1次为止。 5.进行解的倒推。分析一下该算法的不足之处吧，由于我本人能力有限，对于C++的特性还不是很熟悉
1 由于这个数组传参问题就是我的弱项，我的想法是由用户指定任意大小的数组，但这样就要多一个参数。
2 本人只是做了个double类型的矩阵运算。我的思路是可以将这些修改成函数模板，用模板支持各种数据类型的矩阵。如果矩阵是整数类型的，里面的运算可能会产生double型，进行转换赋值会丢失精度。而且一般来说double型的矩阵计算，由于机器运算的原因，会有计算机产生的误差，所以改成模板，似乎没有什么应用。3.对于空间和性能分析，该算法的局限性很大，首先由于该算法不能对主对角线上为0的元素进行计算，所以局限性很大，关于改进版的算法，接下来我会陆续编写求解方法。其次，对于该算法，我的实现过程用了2个存储n维数组来保存值，对其实可以存储在算矩阵本身，因为消元过程完成后，那个空间基不会再使用了，不加上时间的消耗 ，并不是一个理想的算法，接下里我会改进此算法。
“`
对后续结果产生影响，按理应该能直接存储在矩阵中。 时间分析：该算法时间复杂度是O(N^3)，所以运算很大，int main()
{
double a[3][3] = { 1,1,1,0,4,-1,2,-2,1 };
double b[3] = { 6,5，1}
gaussin(a, b);
return 0;
}
“`这是测试代码，结果运算结果是1,2,3。我时间有限，没有进行很多测试，希望大家能与我多探讨，相互学习和促进，做更多的分析和讨论！如有错误和不适之处，请各位指出！
