【原创】开源Math.NET基础数学类库使用(15)C#计算矩阵行列式

               本博客所有文章分类的总目录：【总目录】本博客博文总目录-实时更新 

开源Math.NET基础数学类库使用总目录：【目录】开源Math.NET基础数学类库使用总目录

　　上个月对Math.NET的基本使用进行了介绍，主要内容有矩阵，向量的相关操作，解析数据格式，数值积分，数据统计，相关函数，求解线性方程组以及随机数发生器的相关内容。这个月接着深入发掘Math.NET的各种功能，并对源代码进行分析，使得大家可以尽可能的使用Math.NET在.NET平台下轻易的开发数学计算相关的，或者可以将其中的源码快速移植到自己的系统中去(有时候并不需要所有的功能，只需要其中的部分功能代码)，今天要介绍的是Math.NET中利用C#计算矩阵行列式的功能。

本文原文地址：http://www.cnblogs.com/asxinyu/p/4304289.html

1.行列式概念与性质

　　行列式是关于方阵的元素所定义的一种运算，其运算的结果是一个数，称为方阵的行列式值，简称为方阵的行列式。

     

　　行列式的概念最初是伴随着方程组的求解而发展起来的。行列式的提出可以追溯到十七世纪，最初的雏形由日本数学家关孝和与德国数学家戈特弗里德·莱布尼茨各自独立得出，时间大致相同。日本数学家关孝和提出来的，他在1683年写了一部名为解伏题之法的著作，意思是“解行列式问题的方法”，书中对行列式的概念和它的展开已经有了清楚的叙述。欧洲第一个提出行列式概念的是德国数学家，微积分学奠基人之一莱布尼茨。

　　行列式是线性代数里面的一个基本概念，我们可以从其定义和性质中了解一下其作用：

1.行列式与它的转置行列式相等；

2.互换行列式的两行(列)，行列式变号；

3.行列式的某一行(列)的所有的元素都乘以同一数k，等于用数k乘此行列式；

4.行列式如果有两行(列)元素成比例，则此行列式等于零；

5.若行列式的某一列(行)的元素都是两数之和，则这个行列式是对应两个行列式的和；

6.把行列式的某一列(行)的各元素乘以同一数然后加到另一列(行)对应的元素上去，行列式不变

2.Math.NET计算行列式的实现

　　Math.NET在对行列式的计算过程中，只是把其作为矩阵计算的一个小部分功能，作为属性添加在各个矩阵分解算法的抽象和实现类中。因为Math.NET中矩阵的泛型类型的相关实现主要是支持Double和Float类型，所以基本上与泛型相关的类都实现了2个版本，在实际使用时再进行调用。而矩阵分解算法如：Cholesky，LU，QR，Svd等都有一个抽象泛型基类。在这些抽象类中都定义好了矩阵分解相关的计算，如行列式，方程求解等功能，然后对类进行继承，如Cholesky分解算法，的抽象 基类：

 1 internal abstract class Cholesky : Cholesky<double>
 2 {
 3     protected Cholesky(Matrix<double> factor)
 4         : base(factor)
 5     {
 6     }
 7 
 8     /// <summary>
 9     /// Gets the determinant of the matrix for which the Cholesky matrix was computed.
10     /// </summary>
11     public override double Determinant
12     {
13         get
14         {
15             var det = 1.0;
16             for (var j = 0; j < Factor.RowCount; j++)
17             {
18                 var d = Factor.At(j, j);
19                 det *= d*d;
20