C#使用MathNet库进行矩阵运算

    Math.net计算库是专为.NET开发的开源科学计算库。包含了大部分矩阵运算、线性运算等功能。

1.安装

   项目 - 右键 - 管理NuGet程序包

2.矩阵定义

//矩阵的定义和初始化
var matrix1 = new DenseMatrix(3);     //3维方阵
var matrix2 = new DenseMatrix(2, 2);  //2维方阵
var matrix3 = new DenseMatrix(2, 3);  //2×3矩阵
var matrix4 = new DenseMatrix(3, 2);  //3×2矩阵
var matrix5 = DenseMatrix.CreateIdentity(3);  //单位矩阵
double[,] x = { { 1.0, 2.0,3.0 },
                { 3.0, 4.0,5.0 },
                { 3.0, 6.0,5.0 } };
var matrix6 =  DenseMatrix.OfArray(x);

3.矩阵操作

double[,] x = { { 1.0, 2.0,3.0 },
                { 3.0, 4.0,5.0 },
                { 3.0, 6.0,5.0 } };
var matrix = DenseMatrix.OfArray(x);
var submatrix = matrix.SubMatrix(1, 2, 1,2); //取从第1行开始的2行，第1列开始的2列 子矩阵
var row = matrix.Row(1, 1, 2);                //取从第1行第1列开始的2个行元素 
var column = matrix.Column(1, 1, 2);          //取从第1列第1行开始的2个列元素      
var data1 = matrix.ToRowArrays();             //矩阵变为行向量
var data2 = matrix.ToColumnArrays();          //矩阵变为列向量
var data3 = matrix.Diagonal();                //取矩阵的对角线元素向量

4.矩阵运算

double[,] x = { { 1.0, 2.0,3.0 },
                { 3.0, 4.0,5.0 },
                { 3.0, 6.0,5.0 } };
var m1 = DenseMatrix.OfArray(x);
var m2 = DenseMatrix.CreateIdentity(3);  //单位矩阵

var result = m1 * m2;
bool flag = result.Equals(m1);   //true
result = m1 + m2;
flag = result.Equals(m1);        //false