由于矩阵经有限次初等变换秩不变,且行阶梯形矩阵的秩就等于非零行的行数.因此考虑将矩阵化为行简化阶梯矩阵,确定其非零行数,即为原矩阵的秩.
clear
A=[2,3,1,-3,-7;1,2,0,-2,-4;3,-2,8,3,0;2,-3,7,4,3]
B=rref(A) %求A的行简化阶梯矩阵B
[m,n]=size(B);
while B(m,:)==0 %求B的非零行数
m=m-1;
end
m
A =
2 3 1 -3 -7
1 2 0 -2 -4
3 -2 8 3 0
2 -3 7 4 3
B =
1 0 2 0 -2
0 1 -1 0 3
0 0 0 1 4
0 0 0 0 0
m = 3
即R(A)=3