2022年二阶段法matlab矩阵描述 .pdf

《2022年二阶段法matlab矩阵描述 .pdf》由会员分享，可在线阅读，更多相关《2022年二阶段法matlab矩阵描述 .pdf（11页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、functionTwo_stage(zflag,b,A,Aflag,c) clear;clc format rat fprintf(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

2、=-1 -1 1; %-标准化ifzflag=-1 c=-1*c; end nb=size(b); mA,nA=size(A); fori=1:nb if b(i)0 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 11 页 - - - - - - - - - b(i)=-1*b(i); A(i,:)=-1*A(i,:); end end num1=find(Aflag=-1); num2=find(Aflag=0); num3=find(Aflag=1); N=; c1

3、=zeros(1,length(c); %构造第一阶段目标函数系数矩阵fori=1:length(num1) A(num1(i),nA+1)=1; c(nA+1)=0; c1(nA+1)=0; nA=nA+1; N(i)=nA; end fori=1:length(num3) A(num3(i),nA+1)=-1; A(num3(i),nA+2)=1; c(nA+1)=0; c1(nA+1)=0; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 11 页 - - - -

4、- - - - - c1(nA+2)=-1; nA=nA+2; N(end+1)=nA; end fori=1:length(num2) A(num2(i),nA+1)=1; c1(nA+1)=-1; nA=nA+1; N(end+1)=nA; end %以下部分初始化基变量的下标值矩阵N,基变量的系数矩阵CB %- CB=zeros(1,mA); mA,nA=size(A); fori=1:mA CB(i)=c1(N(i) ; %确定基变量的系数end if (length(num2)+length(num3)=0 A,b,c1,N,CB=first_stage(A,b,c1,N,CB);

5、A,b,c,N,CB=second_stage(A,b,c,c1,N,CB); else fprintf(标准化后初始单纯形表为：); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 11 页 - - - - - - - - - c p=CB,N,b,A disp(*);A,b,c,N,CB,t1=second_stage_diedai(A,b,c,N,CB);%只调用第二阶段迭代过程end % 单纯形法第一阶段函数%= function A,b,c1,N,CB=firs

6、t_stage(A,b,c1,N,CB) k=1; % k 是times=0; Isend=0; temp1= 第;temp2=次后单纯形表为： ; b=b; whileIsend=0 if (times=0) fprintf(标准化后初始单纯形表为：); else 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 11 页 - - - - - - - - - fprintf(%s%d%s,temp1,times,temp2); end times=times+1; mA,

7、nA=size(A); t=zeros(1,mA); for j=1:nA %start for t(j)=c1(j)-CB*A(:,j); end %End for c1 p=CB,N,b,A t m1,j=max(t); %确定导入基，导入基的列下标为col k=max(t); if m10 f(i)=b(i)/A(i,j); else f(i)=10000; % 如果 相除 是负 值， 把值 设 为10000 times2=times2+1; end end if times2=mA fprintf(此问题的无最优解 n) Isend=1; break; %跳出 if m10 fori=

8、1:mA if A(i,j)*b(i)0 f(i)=b(i)/A(i,j); else f(i)=10000; %如果相除是负值，把值设为10000 ftimes=ftimes+1; end end ifftimes=mA fprintf(此问题的无最优解 n) Isend=1; else m4,row=min(f); CB(row)=c(column); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 11 页 - - - - - - - - - N(row)=colu

9、mn; B=; fori=1:length(N) B=B,A(:,N(i); end b=Bb; % end the for cycle,已改变了 b 数组的值A=BA; end else Isend=1; X=zeros(1,nA); fori=1:length(N) X(N(i)=b(i); fprintf(此问题的最优解为： n) X fprintf(目标函数的最优值为： n) Z=c*(X) disp(*);end end %end if end%end while 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 11 页 - - - - - - - - - %*第二阶 段迭 代函数结束名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 11 页 - - - - - - - - -