juzhenqiuzhi.m

来自「计算不等式约束情况下二次规划问题」· M 代码 · 共 49 行

function t=JuZhenQiuZhi(A)
  [m,n]=size(A);
   t=0;
  if (m>n) 
      nn=n;
  else
      nn=m;
  end
    
  for k=1:nn
      q=0;
       for i=k:m
           for j=k:n
               d=abs(A(i,j));
               if (d>q)
                   q=d;
                   hi=i;
                   lj=j;
               end
           end
       end
       if (q+1==1)
           return ;
       end
       t=t+1;
       if (hi~=k)
           for j=k:n
               d=A(k,j);
               A(k,j)=A(hi,j);
               A(hi,j)=d;
           end
       end
       if (lj~=k)
           for i=k:m
               d=A(i,k);
               A(i,k)=A(i,lj);
               A(i,lj)=d;
           end
       end
       for i=k+1:m
           d=A(i,k)/A(k,k);
           for j=k+1:n
               A(i,j)=A(i,j)-d*A(k,j);
           end
       end
  end       
    return