example3_12.m

来自「Matlab 6.X辅助优化计算与设计----源代码」· M 代码 · 共 21 行

n=7;
m=10;
p=0.8;		
    % 生成满足KKT条件的随机数据
    A=rand(m,n);
    c=rand(n,1);
d=rand(n,1);
G=diag(d);
ydes=rand(m,1)-p;
act=( ydes>p );		
ydes(act)=0*ydes(act);
eq=( ydes<0 );		
xdes=(A'*ydes-c)./d;
res=rand(m,1);
res(~act)=0*res(~act);
b=A*xdes-res;
prt=0;
     
disp('test of minqdef.m');
[x,y,ier]=minqdef(c,G,A,b,eq,prt);