jichujuzhen.m

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

function  [aa,bb]=JiChuJuZhen(A,B)
    mm=JuZhenQiuZhi(A);
    [m,n]=size(A);
    [m1,n1]=size(B); 
    for k=1:mm 
        i=1;
        while(i<=m)
            for j=1:n
                gg(k,j)=A(i,j);
            end
                s=JuZhenQiuZhi(gg);
            if (s==k)
                hh(1,k)=i;
                break;
            else
                i=i+1;
            end
        end
    end
   
    for i=1:mm
        k=hh(1,i);
        for j=1:n
            a1(i,j)=A(k,j);
        end
    end
        
    for i=1:mm
        k=hh(1,i);
        for j=1:n1
            b1(i,j)=B(k,j);
        end
    end
    aa=a1;
    bb=b1;