test.m

来自「以上是两个内点法程序」· M 代码 · 共 28 行

%本程序的功能是用于最优性检验
%improve表示调用迭代改进程序
%i,j表示循环变量
for i=1:n
    D(1,i)=X(i);
    D_ni(1,i)=1./X(i);
end
for i=1:m
    for j=1:n
        Z_Z(i,j)=A(i,j)*D(1,j)^2;
    end
end
W=inv(Z_Z*A')*Z_Z*C;
R=C-A'*W;
G=0;
for i=1:n
    if R(i,1)<0
        G=G+1;
    end
end
if G==0&D*R<pr
    disp('原始近优解X为：');
    disp(X);
    disp('对偶近优解W为：');
    disp(W);
else
    improve
end