Question

C++ 高斯消元法解方程组，3阶即可。

求一份可以通过运行的代码，尽可能用简单的语句，因为才开始学，谢谢。
我也蛋疼啊，matlab可以算的，老板要求用C++编。。。唉。。蛋都碎了。。

如果谁有兴趣的话，这个问题也可以回答哦，也是我提的。不过目前还没人回答的。。
http://zhidao.baidu.com/question/266464294.html

Answer 1

计算方法留的上机作业，难度比前几次提高数倍。感觉从来还没写过这么复杂的数学题程序- -!

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

#define N 10 //矩阵大小范围

/*
* 使用已经求出的x，向前计算x（供getx()调用）
* float a[][] 系数矩阵
* float x[] 方程组解
* int i 解的序号
* int n 矩阵大小
* return 公式中需要的和
*/
float getm(float a[N][N], float x[N], int i, int n)
{
float m = 0;
int r;
for(r=i+1; r<n; r++)
{
m += a[i][r] * x[r];
}
return m;
}

/*
* 解方程组，计算x
* float a[][] 系数矩阵
* float b[] 右端项
* float x[] 方程组解
* int i 解的序号
* int n 矩阵大小
* return 方程组的第i个解
*/
float getx(float a[N][N], float b[N], float x[N], int i, int n)
{
float result;
if(i==n-1) //计算最后一个x的值
result = float(b[n-1]/a[n-1][n-1]);
else //计算其他x值（对于公式中的求和部分，需要调用getm()函数）
result = float((b[i]-getm(a,x,i,n))/a[i][i]);

return result;
}

void main()
{
//float a[N][N] = ,,};
//float b[N] = ;
float a[N][N]; //系数矩阵
float b[N]; //右端项
float x[N]; //方程组解
int i,j,k;
int n=N; //矩阵大小

/*用户手工输入矩阵*/
printf("请输入系数矩阵的大小：");
scanf("%d", &n);
printf("请连续输入矩阵值：");
for(i=0; i<n; i++)
{
for(j=0; j<n; j++)
scanf("%f", &a[i][j]);
}
printf("请输入右端项：");
for(i=0; i<n; i++)
{
scanf("%f", &b[i]);
}

/*显示原始矩阵*/
printf("\n原始矩阵\n");
for(i=0; i<n; i++)
{
for(j=0; j<n; j++)
printf("%f ",a[i][j]);
printf("\t|\t%f\n",b[i]);
}
printf("\n\n");

/*进行高斯消去*/
for(j=0; j<n-1; j++)
{
for(i=j+1; i<n; i++)
{
float m = (float)(a[i][j]/a[j][j]);
for(k=j; k<n; k++)
{
a[i][k] = a[i][k]-m*a[j][k];
}
b[i] = b[i]-m*b[j];
}
}

/*显示处理后矩阵*/
printf("高斯消去后矩阵\n");
for(i=0; i<n; i++)
{
for(j=0; j<n; j++)
printf("%f ",a[i][j]);
printf("\t|\t%f\n",b[i]);
}

/*回代方式解方程组*/
for(i=n-1; i>=0; i--)
{
x[i] = getx(a,b,x,i,n);
}

/*显示方程组解*/
printf("\n\n方程组解\n");
for(i=0; i<n; i++)
{
printf("x%d = %f\n", i+1,x[i]);
}
}

VC6.0下测试通过

===================================================================

Answer 2

原来写过一个，你看看能不能用：
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
float *ColPivot(float *c,int n);
void main()
{
int i,j;
float *x;
int n=4;
float a[4][5]={4.,-1.,0.,2.,-1.0,-1,4,-1,0,-7.0,0,-1,4,-1,9.0,2,0,-1,4.0,0.0};
x=ColPivot(a[0],n);
for (i=0;i<n;i++) printf("x[%d]=%f\n",i,x[i]);
getch();
}

float *ColPivot(float *c,int n)
{
int i,j,t,k;
float *x,p;
x=(float *)malloc(n*sizeof(float));
for(i=0;i<=n-2;i++)
{
k=i;
for(j=i+1;j<=n-1;j++)
if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i)))) k=j;
if(k!=i)
for(j=i;j<=n;j++)
{
p=*(c+i*(n+1)+j);
*(c+i*(n+1)+j)=*(c+k*(n+1)+j);
*(c+k*(n+1)+j)=p;
}
for(j=i+1;j<=n-1;j++)
{
p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
for(t=i;t<=n;t++)
*(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
}
}
for(i=n-1;i>=0;i--)
{
for(j=n-1;j>=i+1;j--)
(*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
}
return x;
}

追问

不行的哦。。。
很多个找不到标识符

追答

我编译了一下，计算正确。怎会找不到标识符？

追问

怪了也。。我直接把你的所有原版拷贝进去的，我用的visual studio2008，加上了它提示的要加的那个include，然后编译，不成功，提示少标识符，我打算贴上来的，但是超过字数了。

追答

程序你是否copy 错？再见一个：
http://zhidao.baidu.com/question/74759870.html

追问

我拷贝你给的网址里面的那个程序，按vs2008的要求添加上了#include "stdafx.h"，但是仍然出错，说是getch和malloc找不到标识符。

追答

stdafx.h -- MFC 的头文件
malloc 试 malloc.h
getch -- 删去。system("pause"); 也可以

Answer 3

这个是数值分析的知识，，
以前编过，机子重装后没了，，编太费事，希望你遇到高手

追问

呵呵，也谢谢了哟·~~

Answer 4

这个用矩阵可能会简单得多，虽然我学了线性代数、、但把它编出来、、
呵呵、我有待加强