基于拉普拉斯公式的行列式递归函数去循环优化咨询

我正尝试基于拉普拉斯公式构建一个计算行列式的递归函数，代码如下：

function y=detm(M)
n=length(M);
if n==2
    y=M(1,1).*M(2,2)-(M(1,2).*M(2,1));
else
    for i=1:n
        y=sum(((-1)^(1+n).*M(1,n).*detm(M([1:n-1,n+1:end],[1:n-1,n+1:end]))));
    end
end

目前的问题是，由于需要遍历1到n列，我无法找到去除for循环的方法。请问能否在当前函数内实现该需求，还是需要创建另一个递归函数？

嘿，我看到你在尝试用拉普拉斯公式写递归的行列式计算函数，而且想去掉for循环——其实完全可以在当前函数内实现，不用额外写另一个递归函数，不过先得修正你代码里的一个小问题：你在else块的for循环里误用了n作为列索引（应该是i），而且每次循环都重新计算sum，这会导致结果错误。

用向量化操作替代for循环的实现方案 #

我们可以利用Matlab的arrayfun函数来批量处理每一列的递归计算，同时结合数组操作来生成符号因子和子矩阵索引，完全去掉for循环。修改后的代码如下：

function y = detm(M)
    n = length(M);
    if n == 2
        % 2阶矩阵的行列式直接计算
        y = M(1,1)*M(2,2) - M(1,2)*M(2,1);
    else
        % 生成1到n的列索引数组
        col_indices = 1:n;
        % 计算每个列对应的拉普拉斯符号因子 (-1)^(1+i)
        sign_factors = (-1).^(1 + col_indices);
        % 提取矩阵第一行的所有元素
        first_row_elements = M(1, :);
        
        % 使用arrayfun遍历每个列i，计算去掉第i列和第一行后的子矩阵行列式
        sub_det_values = arrayfun(@(i) detm(M(2:end, col_indices ~= i)), col_indices);
        
        % 按照拉普拉斯公式求和得到最终行列式
        y = sum(sign_factors .* first_row_elements .* sub_det_values);
    end
end

代码解释 #

1. 符号因子生成：(-1).^(1 + col_indices)会生成一个和列数长度相同的数组，每个元素对应拉普拉斯展开中第i列的符号。
2. 子矩阵索引：col_indices ~= i会生成一个逻辑数组，标记出除了第i列之外的所有列，结合M(2:end, ...)就得到了去掉第一行和第i列的子矩阵。
3. 批量递归计算：arrayfun会自动遍历每个列索引i，调用detm计算对应子矩阵的行列式，返回一个包含所有子行列式的数组。
4. 求和运算：最后把符号因子、第一行元素、子行列式三者相乘后求和，就得到了原矩阵的行列式。

这样修改后，既去掉了for循环，又保证了拉普拉斯公式的正确实现，而且不需要额外的递归函数。

需要注意的是，虽然这个实现是学习拉普拉斯展开的好例子，但Matlab内置的det函数经过高度优化，实际使用中效率会比这个递归实现高很多哦。

内容的提问来源于stack exchange，提问作者gbox