C# 实现矩阵运算：乘法、数乘、求逆与随机数生成

class Matrix
    {
        public struct Complex
        {
            public double reg;public double img;
        }
        
        public static double[,] MatrixMult(double[,] matrix1, double[,] matrix2)
        
            //matrix1是m*n矩阵，matrix2是n*p矩阵，则result是m*p矩阵
        { 
            int m = matrix1.GetLength(0), n = matrix2.GetLength(0), p = matrix2.GetLength(1);
            double[,] result = new double[m,p];
            //矩阵乘法：c[i,j]=Sigma(k=1→n,a[i,k]*b[k,j])
            for (int i = 0; i < m; i++)
            {
                for (int j = 0; j < p; j++)
                {
                    //对乘加法则
                    for (int k = 0; k < n; k++)
                    {
                        result[i,j] += (matrix1[i,k] * matrix2[k,j]);
                    }
                }
            }
            return result;
        }
        ////end矩阵相乘
        public static double[,] MatrixNumberMult(double[,] dMatrix, double Num)
        {

            int m = dMatrix.GetLength(0);
            int n = dMatrix.GetLength(1);
            for (int i = 0; i < m; i++)
                for (int j = 0; j < n; j++)
                    dMatrix[i, j] = Num * dMatrix[i, j];
            return dMatrix;
        }

        ///矩阵求逆

        private static double MatrixValue(double[,] MatrixList, int Level)  //求得|A| 如果为0 说明不可逆
        {
            //计算行列式的方法
            //   a1 a2 a3
            //  b1 b2 b3
            //  c1 c2 c3
            // 结果为 a1・b2・c3+b1・c2・a3+c1・a2・b3-a3・b2・c1-b3・c2・a1-c3・a2・b1(注意对角线就容易记住了）
            double[,] dMatrix = new double[Level, Level];   //定义二维数组，行列数相同
            for (int i = 0; i < Level; i++)
                for (int j = 0; j < Level; j++)
                    dMatrix[i, j] = MatrixList[i, j];     //将参数的值，付给定义的数组

            double c, x;
            int k = 1;
            for (int i = 0, j = 0; i < Level && j < Level; i++, j++)
            {
                if (dMatrix[i, j] == 0)   //判断对角线上的数据是否为0
                {
                    int m = i;
                    for (; dMatrix[m, j] == 0; m++) ;  //如果对角线上数据为0，从该数据开始依次往后判断是否为0