【文章內(nèi)容簡(jiǎn)介】 _m。 public int M { get { return _m。 } set { _m = value。 } } private int _n。 public int N { get { return _n。 } set { _n = value。 } } private double[,] _value。 /// summary /// 構(gòu)造函數(shù) /// /summary /// param name=m/param /// param name=n/param public Matrix(int m, int n) { _m = m。 _n = n。 _value = new double[m, n]。 for (int i = 0。 i _m。 i++) { for (int j = 0。 j _n。 j++) { _value[i, j] = 。 } } } /// summary /// 定義的二維索引器，用于獲取矩陣的元素值 /// /summary /// param name=m/param /// param name=n/param /// returns/returns public double this[int m, int n] { get { if (m = _m || n = _n) { throw new ApplicationException(下標(biāo)越界)。 } return _value[m, n]。 } set { if (m = _m || n = _n) { throw new ApplicationException(下標(biāo)越界)。 } _value[m, n] = value。 } } public double MaxValue { get { double max = 1。 for (int i = 0。 i _m。 i++) { for (int j = 0。 j _n。 j++) { if ((_value[i, j]) max) { max = (_value[i, j])。 } } } return max。 } } /// summary /// 運(yùn)算符重載，用于矩陣相加 /// /summary /// param name=a/param /// param name=b/param /// returns返回一個(gè)Matrix類型的對(duì)象/returns static public Matrix operator +(Matrix a, Matrix b) { if ( == amp。amp。 == ) { Matrix c = new Matrix(, )。 for (int i = 0。 i 。 i++) { for (int j = 0。 j 。 j++) { c[i, j] = a[i, j] + b[i, j]。 } } return c。 } else { throw new ApplicationException(矩陣結(jié)構(gòu)不同，無(wú)法相加！)。 } } /// summary /// 運(yùn)算符重載，用于矩陣相減 /// /summary /// param name=a/param /// param name=b/param /// returns返回一個(gè)Matrix類型的對(duì)象/returns static public Matrix operator (Matrix a, Matrix b) { if ( == amp。amp。 == ) { Matrix c = new Matrix(, )。 for (int i = 0。 i 。 i++) { for (int j = 0。 j 。 j++) { c[i, j] = a[i, j] b[i, j]。 } } return c。 } else { throw new ApplicationException(矩陣結(jié)構(gòu)不同，無(wú)法相減)。 } } /// summary /// 運(yùn)算符重載，用于矩陣相乘 /// /summary /// param name=a/param /// param name=b/param /// returns返回一個(gè)Matrix類型的對(duì)象/returns static public Matrix operator *(Matrix a, Matrix b) { if ( == ) { Matrix c = new Matrix(, )。 for (int i = 0。 i 。 i++) { for (int j = 0。 j 。 j++) { for (int k = 0。 k 。 k++) { c[i, j] += a[i, k] * b[k, j]。 } } } return c。 } else { throw new ApplicationException(無(wú)法相乘！)。 } } /// summary /// 靜態(tài)函數(shù)，用于矩陣求逆 /// /summary /// param name=b/param /// returns返回一個(gè)Matrix對(duì)象/returns static public Matrix Rinv(Matrix b) { //復(fù)制矩陣 Matrix a = new Matrix(, )。 for (int ii = 0。 ii 。 ii++) { for (int jj = 0。 jj 。 jj++) { a[ii, jj] = b[ii, jj]。 } } if ( != ) { throw new ApplicationException(只有方陣才能求逆！)。 } int[] _is = new int[]。 int[] _js = new int[]。 int i, j, k。 int n = 。 double d, p。 //判斷是否為奇異陣 for (k = 0。 k = n 1。 k++) { d = 。 for (i = k。 i = n 1。 i++) { for (j = k。 j = n 1。 j++) { //l=i*n+j。 p = (a[i, j])。 if (p d) { d = p。 _is[k] = i。 _js[k] = j。 } } } if (d + == ) { throw new ApplicationException(當(dāng)前矩陣為奇異陣，無(wú)法求逆！)。 } if (_is[k] != k) { for (j = 0。 j = n 1。 j++) { p = a[k, j]。 a[k, j] = a[_is[k], j]。 a[_is[k], j] = p。 }