面向?qū)ο蟪绦蛟O(shè)計(jì)課程設(shè)計(jì)報(bào)告(編輯修改稿)
2025-02-14 23:44 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 for (int i=0。 iindexm。 i++) { for (int j=0。 jindexn。 j++) { [i*() +j] = [i*+j]。 } for (int k=indexn+1。 k。k++) { [i *()+k1] = [i*+k]。 } } for (int m=indexm+1。 m。 m++) { for (int j=0。j。j++) { [(m1)*()+j]=[m*+j]。 } for (int k=indexn+1。 k。k++) { [(m1)*()+k1]=[m*+k]。 } } return adj。}double Det(Matrix amp。a) //遞歸求行列式{ double det = 0。 if ( != ) { cout不是方陣，沒(méi)有行列式!endl。 cout求行列式退出endl。 } if ( == 1) { det = [0]。 return det。 } else { for (int i = 0。 i 。 i++) { if (i % 2 == 0) det += [i * ] * Det(Adjunct(a, i, 0))。 else det = [i * ] * Det(Adjunct(a, i, 0))。 } } return det。}Matrix Inv(Matrix amp。a) //求矩陣的逆 { Matrix temp 。 =。 =。 = new double[ * ]。 double det = Det(a)。 //矩陣的逆 = 伴隨矩陣 / 行列式 if (det == 0) //如果行列式的值為0，則沒(méi)有逆 { cout此矩陣沒(méi)有逆!endl。 cout求矩陣逆退出!。 exit(0)。 } for (int i=0。 i。i++) { for (int j=0。 j。j++) { if ((i+j) % 2 == 0) [i*+j]=Det(Adjunct(a,i,j))/det。 else [i*+j]=Det(Adjunct(a,i,j))/det。 } } return temp。} 求矩陣對(duì)應(yīng)的行列式值和求矩陣(方陣)逆設(shè)計(jì)模塊測(cè)試情況注意：本程序輸出的逆矩陣是先輸出列的，如上圖所示。按照正常的顯示是先輸出行的。上圖中的a矩陣的逆矩陣應(yīng)看成: 2 1 附錄：程序代碼//，Matrix類(lèi)的定義includeiostreamusing namespace std。typedef double Datatype。 class Matrix{public: Matrix()。 //無(wú)參構(gòu)造函數(shù) friend void In(Matrix amp。)。 //設(shè)置為友元的輸入函數(shù) friend void Out(Matrix amp。)。 //設(shè)置為友元的輸出函數(shù) Matrix operator +(Matrix amp。)。 //加法重載函數(shù)Matrix operator (Matrix amp。)。 //減法重載函數(shù) Matrix operator *(Matrix amp。)。 //乘法重載函數(shù) friend Matrix change(Matrix amp。)。 //轉(zhuǎn)置函數(shù) friend Matrix Adjunct
研究報(bào)告相關(guān)推薦
鄂ICP備17016276號(hào)-1