那时，线性代数刚学到行列式。突然有种冲动想做个行列式计算器。于是就有了下面这个n阶行列式计算器。(●—●)
运行效果
思路
其实用的是最简单的方式: 按行列式的第一行展开。
算法用的是递归。说的简单的点，类似数列的递推。如果: n阶行列式的值是a(n)，那么按照第一行展开，会得到余子式，而余子式是n-1阶的行列式，那么可以记为a(n-1)，其计算方法是和计算a(n)一样的。所以可以构造一个函数f，计算f(n)要用到f(n-1)，同理一直推下去，直到f(2)是一个固定的值。这个值会返回代入f(3)并计算出f(3)的值，如此直到计算出f(n)的值。
使用工具
- 开发工具: Visual Studio 2015
- 开发语言: C#
下面贴出核心代码
public double[,] cofactor(double[,] x, int n, int a)
{
double[,] result = new double[n - 1, n - 1];
for (int j = 0; j < a; j++)
{
for (int i = 1; i < n; i++)
{
result[i - 1, j] = x[i, j];
}
}
for (int j = a + 1; j < n; j++)
{
for (int i = 1; i < n; i++)
{
result[i - 1, j - 1] = x[i, j];
}
}
return result;
}
以上是求一个余子式的函数。最终返回n-1阶数组的结果。
其中第一行的参数中：x为原矩阵，n为阶数，a为按第一行展开到了第几个数。
下面是求行列式的函数，调用到了Pow(求幂)和上面的求余子式的函数。
public double det(double[,] x, int n)
{
double result = 0;
if (n == 2)
{
return x[0, 0] * x[1, 1] - x[0, 1] * x[1, 0];
}
else
{
for (int i = 0; i < n; i++)
{
result = result + x[0, i] * Math.Pow(-1, i) * det(cofactor(x, n, i), n - 1); //按第一行展开，降阶算
}
return result;
}
}
再加上输入和输出，就是一个完整的行列式计算器。
代码并不长，也没用到什么高级的东西。
仅供参考，欢迎指出问题。
©Fing
2015-10-21
