improverref.m

来自「有兴趣学系matlab的朋友」· M 代码 · 共 20 行

function A=improverref(A)
%理论不严密，仅可用于低阶解，且要验证解。
[m n]=size(A);
A=rref(A);
i=1;
j=1;
while j<=m
    i=1;
     while i<=10000
         if any(mod(A(j,1:n)*i,1))==0
            A(j,1:n)=mod(A(j,1:n)*i,2);
             break;
             %         else if i==10000 disp('error!'); end
         end
         i=i+1;
     end
     j=j+1;
 end