yueshuji0.m

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

function w=YueShuJi0(A,Xo,B)   
% A 为约束方程系数矩阵，B为约束方程常量，
% Xo为可行点，W为有效约束集,输出W
 m=size(A,1);
 for i=1:m  
   a=A(i,:);    %a为某一约束函数系数
   b=B(i,:);    %b为某一约束函数常数阵
   f=a*Xo-b;    %为约束函数
    if (f==0)   w(1,i)=i; 
    else   w(1,i)=0;
    end
 end