main.m

以上是两个内点法程序

%本程序的功能是用对偶仿射尺度内点法求解线性规划问题
%C表示n费用系数列向量，A表示m*n维技术系数矩阵，b表示m维右端项列向量
%W表示m维对偶变量列向量，S表示n维松弛变量列向量
%pr表示求解精度
C=input('请输入费用系数向量C：');
A=input('请输入技术系数矩阵A：');
b=input('请输入右端项向量b：');
W=input('请输入初始可行内点W：');
S=input('请输入初始可行内点S：');
pr=input('请输入求解精度pr:');
%m，n分别表示矩阵A的行数和列数
[m,n]=size(A);
%D用于表示仿射尺度变换对角阵，用n维行向量表示n*n维对角矩阵以节约存储空间
%D_ni用于表示仿射尺度变换对角阵的逆，用行向量表示以节约存储空间
D=zeros(1,n);D_ni=zeros(1,n);f=0.99;
%调用最优性检测程序
test