C++简单矩阵运算.docVIP

C++简单矩阵运算 ////////////////////////////////////////////////////////////////////////////// ////////////////////////////////Matrix.h//////////////////////////////////// ///////////////////////////////////////////////////////////////////////////// #pragma once #include stdafx.h using namespace std; class Matrix { private: int row, column; double* T; public: Matrix(){}; Matrix(int, int); Matrix(int, int, double*); ~Matrix(); int GetRow(); int GetColumn(); void MultiplyRow(int,double);//某行乘以一个倍数 void ChangeRow(int,int);//互换两行 void AddToRow(int,int,double);//将某行的倍数加到另一行 double GetElement(int,int); //根据行列号获取某个元素 void SetElement(int,int,double); //根据行列号设置某个元素的值 Matrix Inverse(); Matrix transpose(); Matrix operator * (Matrix ); void operator = (Matrix ); }; ////////////////////////////////////////////////////////////////////////////// ////////////////////////////////Matrix.cpp//////////////////////////////////// ///////////////////////////////////////////////////////////////////////////// #pragma once #include stdafx.h #include wmatrix.h #include cmath Matrix::Matrix(int r, int c) { if(r 0 c 0) { row = r; column = c; T = new double[r*c]; for(int i = 0; i r; i++) { for(int j = 0; j c; j++) { T[i*c + j] = 0; } } } else { row = -1; column = -1; } } Matrix::Matrix(int r, int c, double* d) { if(r 0 c 0) { row = r; column = c; T = new double[r*c]; for(int i = 0; i r; i++) { for(int j = 0; j c; j++) { T[i*c + j] = d[i*c + j]; } } } else { row = -1; column = -1; } } Matrix Matrix::Inverse() { Matrix nM(-1, -1); if(row 0 column 0) { return nM; } if(row != column) { return nM; } int i,j,k,n; double max,temp; n = row; //单位矩阵 double* e = new double[n * n]; double* old = new double[n * n]; for(i = 0