main.m

以上是两个内点法程序

%本程序的功能是利用原始仿射尺度算法求解线性规划问题
%线性规划问题条件及初始可行内点输入
%费用系数C为n维列向量，技术系数A为m*n维矩阵
%决策变量X为n维列向量，右端项b为m为列向量，求解精度pr为一个大于零的小正实数
C=input('请输入费用系数向量C：');
A=input('请输入技术系数矩阵A：');
b=input('请输入右端项向量b：');
X=input('请输入初始可行内点：');
pr=input('请输入求解精度pr:');
%为了节省存储空间，对角阵D用行向量表示，D_ni表示对角阵D的逆，并表示成行向量
%f为步长因子
[m,n]=size(A);
D=zeros(1,n);D_ni=zeros(1,n);f=0.99;
%test表示调用最优性检测程序
test