Exemple #1
0
        /// <summary>
        /// 求逆
        /// </summary>
        /// <param name="Mat"></param>
        /// <returns></returns>
        public Matrix Inverses(Matrix Mat)
        {
            int n = Mat.RowCount;

            Matrix LU = new Matrix(n, n), Inverse = new Matrix(n, n);
            Vector b = new Vector(n), Indx = new Vector(n);

            if (Mat.ColCount != Mat.RowCount)
            {
                throw new ArgumentException("维数不符合要求!");
            }

            // LU decomposition
            LU = Mat;
            MatrixUtilNew.LU_Decomp(LU, Indx);

            // Solve Ax=b for  unit vectors b_1..b_n

            for (int j = 0; j < n; j++)
            {
                b[j] = 1.0;                            // Set b to j-th unit vector
                MatrixUtilNew.LU_BackSub(LU, Indx, b); // Solve Ax=b
                Inverse.SetCol(j, b);                  // Copy result
            }
            ;

            return(Inverse);
        }
Exemple #2
0
 /// <summary>
 /// Transposition
 /// </summary>
 /// <returns></returns>
 public Matrix Transpose()
 {
     return(MatrixUtilNew.Transpose(this));
 }