Matrix(C#实现).docx

更多资料using System;using System.Collections.Generic;using System.Linq;using System.Text;using System.IO;using System.Diagnostics;namespace WindowsFormsApplication1 public class Matrix public Matrix(int row) m_data = new doublerow, row; public Matrix(int row, int col) m_data = new doublerow, col; /复制构造函数 public Matrix(Matrix m) /有什么用 int row = m.Row; int col = m.Col; m_data = new doublerow, col; for (int i = 0; i row; i+) for (int j = 0; j col; j+) m_datai, j = mi, j; /*/ /* /分配方阵的大小 /对于已含有内存的矩阵，将清空数据 public void SetSize(int row) m_data = new doublerow,row; /分配矩阵的大小 /对于已含有内存的矩阵，将清空数据 public void SetSize(int row,int col) m_data = new doublerow,col; */ /unit matrix:设为单位阵 public void SetUnit() for (int i = 0; i m_data.GetLength(0); i+) for (int j = 0; j m_data.GetLength(1); j+) m_datai, j = (i = j) ? 1 : 0); /设置元素值 public void SetValue(double d) for (int i = 0; i m_data.GetLength(0); i+) for (int j = 0; j m_data.GetLength(1); j+) m_datai, j = d; / Value extraction：返中行数 public int Row get return m_data.GetLength(0); /返回列数 public int Col get return m_data.GetLength(1); /重载索引 /存取数据成员 public double thisint row, int col get return m_datarow, col; set m_datarow, col = value; /primary change /初等变换对调两行：rirj public Matrix Exchange(int i, int j) double temp; for (int k = 0; k Col; k+) temp = m_datai, k; m_datai, k = m_dataj, k; m_dataj, k = temp; return this; /初等变换第index 行乘以mul public Matrix Multiple(int index, double mul) for (int j = 0; j Col; j+) m_dataindex, j *= mul; return this; /初等变换 第src行乘以mul加到第index行 public Matrix MultipleAdd(int index, int src, double mul) for (int j = 0; j Col; j+) m_dataindex, j += m_datasrc, j * mul; return this; /transpose 转置 public Matrix Transpose() Matrix ret = new Matrix(Col, Row); for (int i = 0; i Row; i+) for (int j = 0; j Col; j+) retj, i = m_datai, j; return ret; /binary addition 矩阵加 public static Matrix operator +(Matrix lhs, Matrix rhs) if (lhs.Row != rhs.Row) /异常 System.Exception e = new Exception(相加的两个矩阵的行数不等); throw e; if (lhs.Col != rhs.Col) /异常 System.Exception e = new Exception(相加的两个矩阵的列数不等); throw e; int row = lhs.Row; int col = lhs.Col; Matrix ret = new Matrix(row, col); for (int i = 0; i row; i+) for (int j = 0; j col; j+) double d = lhsi, j + rhsi, j; reti, j = d; return ret; /binary subtraction 矩阵减 public static Matrix operator -(Matrix lhs, Matrix rhs) if (lhs.Row != rhs.Row) /异常 System.Exception e = new Exception(相减的两个矩阵的行数不等); throw e; if (lhs.Col != rhs.Col) /异常 System.Exception e = new Exception(相减的两个矩阵的列数不等); throw e; int row = lhs.Row; int col = lhs.Col; Matrix ret = new Matrix(row, col); for (int i = 0; i row; i+) for (int j = 0; j col; j+) double d = lhsi, j - rhsi, j; reti, j = d; return ret; /binary multiple 矩阵乘 public static Matrix operator *(Matrix lhs, Matrix rhs) if (lhs.Col != rhs.Row) /异常 System.Exception e = new Exception(相乘的两个矩阵的行列数不匹配); throw e; Matrix ret = new Matrix(lhs.Row, rhs.Col); double temp; for (int i = 0; i lhs.Row; i+) for (int j = 0; j rhs.Col; j+) temp = 0; for (int k = 0; k lhs.Col; k+) temp += lhsi, k * rhsk, j; reti, j = temp; return ret; /binary division 矩阵除 public static Matrix operator /(Matrix lhs, Matrix rhs) return lhs * rhs.Inverse(); /unary addition单目加 public static Matrix operator +(Matrix m) Matrix ret = new Matrix(m); return ret; /unary subtraction 单目减 public static Matrix operator -(Matrix m) Matrix ret = new Matrix(m); for (int i = 0; i ret.Row; i+) for (int j = 0; j ret.Col; j+) reti, j = -reti, j; return ret; /number multiple 数乘 public static Matrix operator *(double d, Matrix m) Matrix ret = new Matrix(m); for (int i = 0; i ret.Row; i+) for (int j = 0; j ret.Col; j+) reti, j *= d; return ret; /number division 数除 public static Matrix operator /(double d, Matrix m) return d * m.Inverse(); /功能：返回列主元素的行号 /参数：row为开始查找的行号 /说明：在行号row,Col)范围内查找第row列中绝对值最大的元素，返回所在行号 int Pivot(int row) int index = row; for (int i = row + 1; i m_dataindex, row) index = i; return index; /inversion 逆阵：使用矩阵的初等变换，列主元素消去法 public Matrix Inverse() if (Row != Col) /异常,非方阵 System.Exception e = new Exception(求逆的矩阵不是方阵); throw e; / StreamWriter sw = new StreamWriter(.annexclose_matrix.txt); Matrix tmp = new Matrix(this); Matrix ret = new Matrix(Row); /单位阵 ret.SetUnit(); int maxIndex; double dMul; for (int i = 0; i Row; i+) maxIndex = tmp.Pivot(i); if (tmp.m_datamaxIndex, i = 0) System.Exception e = new Exception(求逆的矩阵的行列式的值等于0,); throw e; if (maxIndex != i) /下三角阵中此列的最大值不在当前行，交换 tmp.Exchange(i, maxIndex); ret.Exchange(i, maxIndex); ret.Multiple(i, 1 / tmpi, i); tmp.Multiple(i, 1 / tmpi, i); for (int j = i + 1; j 0; i-) for (int j = i - 1; j = 0; j-) dMul = -tmpj, i / tmpi, i; tmp.MultipleAdd(j, i, dMul); ret.MultipleAdd(j, i, dMul); /end for / sw.WriteLine(tmp=rn + tmp); / sw.WriteLine(ret=rn + ret); / sw.WriteLine(*=rn + this * ret); / sw.Close(); return ret; /end Inverse #region /*/ /* /inversion 逆阵：使用矩阵的初等变换，列主元素消去法 public Matrix Inverse() if(Row != Col) /异常,非方阵 System.Exception e = new Exception(求逆的矩阵不是方阵); throw e; / StreamWriter sw = new StreamWriter(.annexmatrix_mul.txt); / / Matrix tmp = new Matrix(this); Matrix ret =new Matrix(Row); /单位阵 ret.SetUnit(); int maxIndex; double dMul; for(int i=0;iRow;i+) maxIndex = tmp.Pivot(i); if(tmp.m_datamaxIndex,i=0) System.Exception e = new Exception(求逆的矩阵的行列式的值等于0,); throw e; if(maxIndex != i) /下三角阵中此列的最大值不在当前行，交换 tmp.Exchange(i,maxIndex); ret.Exchange(i,maxIndex); ret.Multiple(i,1/tmpi,i); / /sw.WriteLine(nul t+tmpi,i+t+reti,i); / tmp.Multiple(i,1/tmpi,i); /sw.WriteLine(mmm t+tmpi,i+t+reti,i); sw.WriteLine(111111111 tmp=rn+tmp); for(int j=i+1;j0;i-) for(int j=i-1;j=0;j-) dMul = -tmpj,i; tmp.MultipleAdd(j,i,dMul); ret.MultipleAdd(j,