wuyueshu0.m

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

function [flag,x1]=Wuyueshu0(A,B,G,D)
 m=size(A,1); %无约束时求最优解
 x=inv(G)*(-D);
 flag=1;
 for i=1:m     
   a=A(i,:);    %a为某一约束函数系数
   b=B(i,:);    %b为某一约束函数常数阵
   f=a*x-b;     %为约束函数
   if (f<0) 
      flag=0; break;
   end
 end
 if (flag==1)
     x1=x;
 else
     x1=[];
 end