用高斯消去法解线性方程组。C++程序及结果.docVIP

用高斯消去法解线性方程组。设有n元线性方程组如下： 算法： 对于k从0开始到n-2结束，进行以下步骤： 首先，第k行第k列的数不能为0，若为0，则与下面不为0的行交换。将k行所有数除以第k行第k列的数. 其次，进行消去： a[i][j]=a[i][j]-a[i][k]*a[k][j],j,i=k+1,k+2,…，n-1 b[i]=b[i]-a[i][k]*b[k],i=k+1,k+2,……,n-1 最后，回代过程 x[n-1]=b[n-1]/a[n-1][n-1] x[i]=b[i]-∑a[i][j]x[j],I=n-2,…1,0 定义一个矩阵类Matrix作为基类，然后由矩阵类派生出线性方程组类Linequ。 程序清单： #includeiostream #includecmath using namespace std; class Matrix //定义矩阵类 { public: Matrix(int dims=2) //构造函数 { index=dims; //保护数据赋值 MatrixA=new double[index*index]; //动态内存分配 } ~Matrix() {delete[] MatrixA;} //内存释放 void setMatrix(double *rmatr) //设置矩阵值 { for(int i=0;iindex*index;i++) {*(MatrixA+i)=rmatr[i]; } //矩阵成员赋初值 } void printM(); //显示矩阵 protected: int index; double* MatrixA; }; class Linequ:public Matrix //线性方程组类 { public: Linequ(int dims=2):Matrix(dims) { sums=new double[dims]; solu=new double[dims]; } ~ Linequ(); void setLinequ(double*a,double*b); //方程赋值 void printL(); //显示方程 int Solve(); //高斯消元法求解 void showX(); //输出方程的解 private: double *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() { delete[]sums; delete[]solu; } void Linequ::setLinequ(double *a,double *b) //方程赋值 { setMatrix(a); for(int i=0;iindex;i++) sums[i]=b[i]; } 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 sums[i]endl; } } void Linequ::showX() //输出方程的解 { coutThe Result is:endl; for(int