f.m

来自「有兴趣学系matlab的朋友」· M 代码 · 共 19 行

function y=f(D11,D12,A1,b1)
t=cputime;
y=[];
n=0;
bequ=zeros(15,1);
D1=0*D11+1*D12;
for i=1:10
    for j=i+1:10
        for  k=j+1:10
            aequ=[zeros(5,5*(i-1)),eye(5),zeros(5,50-5*i);zeros(5,5*(j-1)),eye(5),zeros(5,50-5*j);zeros(5,5*(k-1)),eye(5),zeros(5,50-5*k)];
            [X,FVAL,FLAG]=linprog(-D1,A1,b1,aequ,bequ,zeros(50,1),[]);
            n=n+1;
            d2=D12'*X;
            d1=D11'*X;
            y=[y,[n;i;j;k;d1;d2;d1+d2;FLAG;X]]; %  y=[y;X]
        end
    end
end
t=cputime-t