buchangyueshu3.m

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

function [alpha,x1,w1]=BuChangYueShu3(A,X,B,P,W)
 m=size(A,1);   %求约束函数系数阵的行数
 temp=0;ak=1;mi=0;
 for i=1:m  
   if(W(1,i)==0)   %alpha
     a=A(i,:);
     b=B(i,:);
     t=a*P;
     if (t<0)
        temp=(b-a*X)/t;
        if (temp<ak)  %alpha<1
           mi=i;
           ak=temp;
        end
     end
   end
 end
 if (mi~=0)   %alpha<1,增加约束
     W(1,mi)=mi;
 end
 w1=W;
 alpha=ak;
 x1=X+alpha*P; %调整X