两阶段算法求解线性规划问题(运筹学实验).docVIP

两阶段算法求解线性规划问题(运筹学实验) %此函数是使用两阶段算法求解线性规划问题 function [x,minf,optmatrx,flag]=linp(A,c,b) %x可行解 minf最优解 oprmatrx最优解对应的单纯形表 flag判断解的情况的标志 fprintf(n 请输入约束条件的系数矩阵！n); A=input(A=); fprintf(n 请输入最右端的列向量!n); b=input(b=);%b为列向量 fprintf(n 请输入目标函数的系数矩阵!n); c=input(c=); format rat %可以让结果用分数输出 [m,n]=size(A); [p,q]=size(b); for i=1p if b(i)0 A(i,)=A(i,)-1; b(i)=b(i)-1; end end %第一阶段辅助问题的求解 auxchart=zeros(m+2,m+n+1);%确定第一张辅助单纯形表，并将其置0 auxchart(1,)=[-c,zeros(1,m+1)];%给表的第一行赋值 auxchart(2,)=[zeros(1,n),-1ones(1,m),0];%给表的第二行赋值 auxchart([3m+2],)=[A,eye(m,m),b];%将变量系数和右端向量填到表中 basvar=[n+1n+m];%基变量是辅助变量 %使基变量的检验数为零 for i=3m+2 auxchart(2,)=auxchart(2,)+auxchart(i,); end %进行辅助问题的换基迭代 checknum1=find(auxchart(2,[1m+n])0);%checknum1向量中记录检验数大于零的位置 while ~isempty(checknum1) column1=checknum1(1);%取第一个大于零的数 postivecheck1=auxchart([3m+2],column1);%postivecheck1为记录检验数大于零的列向量 rightvector1=auxchart([3m+2],m+n+1);%rightvector1为记录右端向量的矩阵 temp1=find(postivecheck10);%在postivecheck1向量中找大于零的位置 if isempty(temp1) fprintf(n 此问题不存在最优解!n); x=[] minf=[] optmatrx=[] flag=-1 return end min1=Inf;%初始化最小值为无穷 for i=3m+2 if postivecheck1(i-2)0rightvector1(i-2)postivecheck1(i-2)min1%找最小比 min1=rightvector1(i-2)postivecheck1(i-2); row1=i;%记录位置 end end auxchart(row1,)=auxchart(row1,)auxchart(row1,column1);%化为一 for i=1m+2 %进行初等行变换 if i~=row1 auxchart(i,)=auxchart(i,)+(-1auxchart(i,column1)auxchart(row1,)); end end basvar(row1-2)=column1;%出基和进基 checknum1=find(auxchart(2,[1m+n])0);%重新找检验向量中大于零的数的位置 end %如果辅助问题的最优解大于零，则原问题无可行解 if auxchart(2,m+n+1)0 fprintf(n 此问题不存在最优解!n); x=[] minf=[] optmatrx=[] flag=-1 return end judge=find(basvarn); %如果辅助问题的最优解等于零，且基变量中无辅助变量，则直接删除表中第二行和人工变量 if auxchart(2,m+n+1)==0isempty(judge) if isempty(judge) auxchart(2,)=[]; auxchart(,[n+1n+m])=[]; end %