function Two_stage(zflag,b,A,Aflag,c)clear;clcformat ratfprintf(n);zflag=input(请输入所求目标函数的状态(求最大值为 1，求最小值为-1):);b=input(请输入资源矩阵 b:);A=input(请输入约束矩阵 A:);Aflag=input(请输入约束方程的状态矩阵(小于号为-1，等于为 0，大于号为 1););c=input(请输入目标函数价值系数矩阵 c:);%zflag=1;%b=9 15 5;%A=5 3 1;-5 6 15;2 1 1;%c=10 15 12;%Aflag=-1 -1 1;%-标准化if zflag=-1c=-1*c;endnb=size(b);mA,nA=size(A);for i=1:nbif b(i)0f(i)=b(i)/A(i,j); elsef(i)=10000; %如果相除是负值，把值设为 10000times2=times2+1;end end if times2=mAfprintf(此问题的无最优解 n)Isend=1;break; %跳出 if m10for i=1:mA if A(i,j)*b(i)0f(i)=b(i)/A(i,j); else f(i)=10000; %如果相除是负值，把值设为 10000ftimes=ftimes+1;endendif ftimes=mAfprintf(此问题的无最优解n)Isend=1;elsem4,row=min(f);CB(row)=c(column);N(row)=column; B=;for i=1:length(N)B=B,A(:,N(i);endb=Bb; % end the for cycle,已改变了 b 数组的值A=BA;endelseIsend=1;X=zeros(1,nA);for i=1:length(N)X(N(i)=b(i); fprintf(此问题的最优解为：n) Xfprintf(目标函数的最优值为：n)Z=c*(X)disp(*);endend %end ifend%end while%*第二阶段迭代函数结束