给定矩阵A，如何用matlab编程求解出所有的最高阶非零子式？

A=[  1 -2  2 3;
    -2  4 -1 3;
    -1  2  0 3;
     0  6  2 3;
     2 -6  3 4];
 [R,C] = size(A);
 Rank = rank(A);
 
 count = 0;
 
 for i = 1 : R-Rank+1
     for j = 1 : C-Rank+1
         B = A(i:i+Rank-1,j:j+Rank-1);
        if rank(B)>0
            B
            count = count + 1;
        end
     end
 end
 
 count

非常感谢

其实有问题，秩如果是k，2楼的程序算的只是连续k行连续k列非零子式，实际上需要算的是任意k行k列的非零子式。
if rank(B)>0 改为 if det(B)~=0，B后加一行det(B)。

[ 本帖最后由 gaoy1 于 2013-7-10 00:30 编辑 ]

确实有问题 感谢你的建议

修改了一下 这样应该可以：

A=[  1 -2  2 3;
     -2  4 -1 3;
     -1  2  0 3;
      0  6  2 3;
      2 -6  3 4];
 [R,C] = size(A);
 Rank = rank(A);

 count = 0;

 for i = 1 : R-Rank+1
      for j = 1 : C-Rank+1
         B = A(i:i+Rank-1,j:j+Rank-1);
         if rank(B) == Rank
             B
             count = count + 1;
         end
      end
 end

 count