【文章內容簡介】 實驗內容：通過構造矩陣類，實現(xiàn)矩陣的定義，包括：矩陣的加減乘除，求矩陣的轉置、求逆等，求矩陣的行列式。主要代碼： class LJH_CMatrix{public: LJH_CMatrix()。 // 默認構造函數(shù) LJH_CMatrix(int row, int column)。 // 構造函數(shù)一 LJH_CMatrix(const LJH_CMatrixamp。 m)。 // 復制構造函數(shù) ~LJH_CMatrix()。 // 默認析構函數(shù) void input()。//矩陣輸入 void output()。 // 輸出該矩陣 LJH_CMatrix transpose()。 // 矩陣轉置 //////////////////////////////////////////////////////////// LJH_CMatrix yuzishi(int i,int j)。//求矩陣的第(i,j)的余子式double hanglieshi()。//求矩陣的行列式 LJH_CMatrix bansui()。//求矩陣的伴隨矩陣 LJH_CMatrix inverse()。 // 矩陣求逆(伴隨矩陣除以行列式) ////////////////////////////////////////////////////////////// LJH_CMatrix inv()。//矩陣求逆(用高斯約當法) LJH_CMatrix amp。 change(int k,int l)。//交換矩陣的第k行和第l行 int max_cloumn(int k)。//求矩陣第k列的最大行數(shù) ////////////////////////////////////////////////////////////// // 設置(i,j)的值 void setValue(int row, int column, double value) { _A[row][column] = value。 } double getValue(int row, int column) const { return _A[row][column]。 } // 設置行、列的值 void setRow(const int row) { _row = row。 } int getRow() const { return _row。 } void setColunm(const int column) { _column = column。 } int getColumn() const { return _column。 }public:// 成員變量 double** _A。 // 或用這個定義vectorvectordouble _A。 int _row, /*行*/ _column。 // 列}。 void LJH_CMatrix::input()//輸入函數(shù) { int i,j。 cout請輸入矩陣的行數(shù)：endl。 cin_row。 cout請輸入矩陣的列數(shù)：endl。 cin_column。 cout請輸入矩陣中的元素：endl。 for(i=0。i_row。i++) for(j=0。j_column。j++) { cin_A[i][j]。 } } void LJH_CMatrix::output()//輸出函數(shù) { int i,j。 for(i=0。i_row。i++) { for(j=0。j_column。j++) { cout_A[i][j] 。 } coutendl。 } } LJH_CMatrix LJH_CMatrix::transpose()//轉置函數(shù) { LJH_CMatrix tem。 =this_column。 =this_row。 int i,j。 for(i=0。i_row。i++) for(j=0。j_column。j++) { [j][i]=_A[i][j]。 } return tem。 } LJH_CMatrix LJH_CMatrix::yuzishi(int i,int j)//求矩陣的余子式 { LJH_CMatrix temp。 =this_row1。 =this_column1。 int m,n,k=0,l。 for(m=0。m_row。m++) { l=0。 for(n=0。n_column。n++) { if(m!=iamp。amp。n!=j) { [k][l]=_A[m][n]。 } if(n!=j) l++。 } if(m!=i) k++。 } return temp。 } double LJH_CMatrix::hanglieshi()//求矩陣的行列式 { if(_row!=_column) {cerr此矩陣無行列式endl。 } if(_row==1amp。amp。_column==1) return _A[0][0]。 else { int i。 double sum=0。 for(i=0。i_column。i++) { sum+=pow(1,i)*_A[0][i]*thisyuzishi(0,i).hanglieshi()。 } return sum。 } } LJH_CMatrix LJH_CMatrix::bansui()//求伴隨矩陣 { LJH_CMatrix temp。 =this_column。 =this_row。 int i,j。 for(i=0。i_row。i++) for(j=0。j_column。j++) { [i][j]=pow(1,i+j)*thisyuzishi(i,j).hanglieshi()。 } return temp。 } LJH_CMatrix LJH_CMatrix::inverse()//矩陣求逆 { LJH_CMatrix temp。 int n。 n=thishanglieshi()。 temp=thisbansui()。 int i,j。 for(i=0。i_row。i++) for(j=0。j_column。j++) { [i][j]/=n。 } return temp。 } LJH_CMatrix amp。 LJH_CMatrix::change(int k,int l)//交換矩陣的第k行和第l行 { int i。 double j。 for(i=0。i_column。i++) { j=_A[k][i]。 _A[k][i]=_A[l][i]。 _A[l][i]=j。 } return *this。 } int LJH_CMatrix::max_cloumn(int k)//求矩陣第k列中從第k個元素之后絕對值最大的行數(shù) { int m=k。 double max=fabs(_A[k][k])。 for(int i=k+1。i_column。i++) { if(fabs(_A[i][k])max) { max=fabs(_A[i][k])。 m=i。 } } return m。 } LJH_CMatrix LJH_CMatrix::inv()//矩陣求逆，通過行列變換 { int i,j,m。 LJH_CMatrix E1。 E1=*this。 if(this_row!=this_column) { cerr該矩陣不能求逆endl。 } else { if(()==0) {cerr該矩陣不可逆：endl。 } else {/////////把矩陣E賦值成單位陣 LJH_CMatrix E。//創(chuàng)建一個和當前方陣階數(shù)相同的單位矩陣 =。 =。 for(i=0。i。i++) for(j=0。j。j++) {[i][j]=0。} for(i=0。i。i++) [i][i]=1。 //化上三角陣 int i,j,hang。 for(i=0。i。i++)// { hang=(i)。 if(hang!=i) { (i,hang)。 (i,hang)。 } double xishu。 for(m=i+1。m。m++)//第m行 { xishu=[m][i]/[i][i]。 for(j=0。j。j++) { [m][j]=[m][j]xishu*[i][j]。 [m][j]=[m][j]xishu*[i][j]。}} ////上面為求上三角陣 } ///////////////////////////化對角陣 for(i=this_column1。i0。i)// { double xishu1。 for(m=0。mi。m++) { xishu1=[m][i]/[i][i]。 for(j=0。jthis_column。j++) { [m][j]=[m][j]xishu1*[i][j]。 [m][j]=[m][j]xishu1*[i][j]。 }}} ///////////////////////////////////矩陣單位化 double xishu3。 for(i=0。i。i++) { xishu3=1/[i][i]。 for(j=0。j。j++) { [i][j]*=xishu3。 [i][j]*=xishu3。}} return E。 }}}void main(){ LJH_CMatrix array1,array2,array3,array4,array5,array6,array7。 cout請輸入矩陣array1相關信息：。 ()。