1.矩阵运算
matrixA = 3.0 * matrixB //数乘
matrixA = vector * matrixB /matrixA = matrixB * vector //向量乘
matrixC = matrixA + / - / * matrixB //矩阵加、减、乘
resultM = (DenseMatrix)matrixA.PointwiseDivide(matrixB); //点乘
上述所有运算符都有对应的函数,也可是利用函数进行运算,如“+”可以写成 matrixC = (DenseMatrix)matrixA.Add(matrixB);或者matrixA.Add(matrixB, matrixC);
matrixB = matrixA.Inverse()/Transpose() //求逆和转置
2.求解线性方程组
对于一阶线性方程组
5*x + 2*y - 4*z = -7
3*x - 7*y + 6*z = 38
4*x + 1*y + 5*z = 43
可以如下求解:
var matrixA = new DenseMatrix(new[,] { { 5.00, 2.00, -4.00 }, { 3.00, -7.00, 6.00 }, { 4.00, 1.00, 5.00 } });
var vectorB = new DenseVector(new[] { -7.0, 38.0, 43.0 });
var resultX = matrixA.LU().Solve(vectorB);
或者
resultX = matrixA.QR().Solve(vectorB);
或者
matrixA.Svd(true).Solve(vectorB, resultX);
或者
matrixA.GramSchmidt().Solve(vectorB, resultX);
3.矩阵秩、行列式、trace和范数
matrix.Determinant()/Rank()/ConditionNumber()/Trace() //行列式/秩/条件数/trace
matrix.L1Norm()/L2Norm()/FrobeniusNorm()/InfinityNorm() //范数