《用高斯消去法解线性方程组C++程序及结果》由会员分享，可在线阅读，更多相关《用高斯消去法解线性方程组C++程序及结果（5页珍藏版）》请在金锄头文库上搜索。
1、用高斯消去法解线性方程组。设有n 元线性方程组如下：a00x0+a01x1+a0,n-1xn-1=b0a10x0+a11x1+a1,n-1xn-1=b1An-1,0x0+an-1,1x1+an-1,n-1xn-1=bn-1算法：对于 k 从 0 开始到 n-2 结束，进行以下步骤：首先，第 k 行第 k 列的数不能为 0，若为 0，则与下面不为 0 的行交换。将 k 行所有数除以第 k 行第 k 列的数.其次，进行消去：aij=aij-aik*akj,j,i=k+1,k+2,，n-1bi=bi-aik*bk,i=k+1,k+2,n-1最后，回代过程xn-1=bn-1/an-1n-1xi=bi-
2、aijxj,I=n-2,1,0定义一个矩阵类 Matrix 作为基类，然后由矩阵类派生出线性方程组类Linequ。程序清单：#include#includeusing namespace std;class Matrix /定义矩阵类public:Matrix(int dims=2) /构造函数index=dims; /保护数据赋值MatrixA=new doubleindex*index; /动态内存分配Matrix()delete MatrixA; /内存释放void setMatrix(double *rmatr) /设置矩阵值for(int i=0;iindex*index;i+)*(
3、MatrixA+i)=rmatri; /矩阵成员赋初值void printM(); /显示矩阵protected:int index;double* MatrixA;class Linequ:public Matrix /线性方程组类public:Linequ(int dims=2):Matrix(dims)sums=new doubledims;solu=new doubledims; Linequ();void setLinequ(double*a,double*b); /void printL(); /int Solve(); /void showX(); /private:double
4、 *sums;double *solu;void Matrix:printM() /显示矩阵的元素coutThe Matrix is:endl;for(int i=0;iindex;i+)for(int j=0;jindex;j+)cout*(MatrixA+i*index+j) ;coutendl;Linequ:Linequ()deletesums;deletesolu;void Linequ:setLinequ(double *a,double *b) /setMatrix(a);for(int i=0;iindex;i+)sumsi=bi;方程赋值显示方程高斯消元法求解输出方程的解方程赋
5、值void Linequ:printL() /显示方程coutThe Line eqution is:endl;for(int i=0;iindex;i+)for(int j=0;jindex;j+)cout*(MatrixA+i*index+j) ;cout sumsiendl;void Linequ:showX() /coutThe Result is:endl;for(int i=0;iindex;i+)coutXi=soluiendl;int Linequ:Solve() /int l,k,i,j,is,p,q,m=0;double d,t;l=1;for(k=0;k=index-2;
6、k+) /d=;for(i=k;i=index-1;i+)for(j=k;jd)d=t;is=i;if(d+=l=0;if(l=0)coutfailendl;return(0);d=MatrixAk*index+k;if(d=0)m+;输出方程的解解线性方程组消去过程d=MatrixA(k+m)*index+k;is=k+m;if(is!=k)for(j=k;j=index-1;j+)p=k*index+j;q=is*index+j;t=MatrixAp;MatrixAp=MatrixAq;MatrixAq=t;t=sumsk;sumsk=sumsis;sumsis=t;d=MatrixAk*
7、index+k; for(j=k;j=index-1;j+) p=k*index+j;MatrixAp=MatrixAp/d; sumsk=sumsk/d; for(i=k+1;i=index-1;i+) for(j=k+1;j=index-1;j+) p=i*index+j; MatrixAp=MatrixAp-MatrixAi*index+k*MatrixAk*index+j; sumsi=sumsi-MatrixAi*index+k*sumsk; d=MatrixA(index-1)*index+index-1;if(fabs(d)+= coutfail=0;i-) t=; for(j=i+1;j=index-1;j+) t=t+MatrixAi*index+j*soluj; solui=sumsi-t;int main() /主函数double a=2,0, ,2,0, ,0,2, 0,2;double b4=0,3,3,0;Linequ equl(4);(a,b);();if()();elsecoutFailendl;return 0;运行结果：
