实验1 单纯形法求解线性规划.docx

文档介绍：: .
sgma=c-cb*A;h=-1;fori=1:nifsgma(i)>0h=1;enden: .
sgma=c-cb*A;h=-1;fori=1:nifsgma(i)>0h=1;endendendcb;b；fm=sum(cb*b);ifh==-1vv=0;fori=1:mifxx(i)<=m1flg='nofeasibel';xx=[];fm=[];b=[]；
实验1单纯形法求解线性规划
成绩
专业班级信息学号20102姓名报告日期.
实验类型：•验证性实验O综合性实验O设计性实验实验目的：进一步熟练掌握单纯形法求解线性规划。实验内容：单纯形法求解线性规划4个(题目自选)
实验原理线性规划单纯形法(线性规划解有四种情形，唯一最优解，无穷多个最解，无界解，无可行解)
1 实验步骤要求上机实验前先编写出程序代码编辑录入程序调试程序并记录调试过程中出现的问题及修改程序的过程经反复调试后，运行程序并验证程序运行是否正确。
2 记录运行时的输入和输出。
预习编写程序代码：
实验报告：根据实验情况和结果撰写并递交实验报告。
参考程序function[xx,b,fm,sgma,AA,flg]=myprgmh(ml,m,n,A,b,c)%单纯形法求解性规划函数。默认标准型人工变量在最前、剩余变量在后构成基本量；%ml人工变量的个数;m基变量的个数；n所有变量的个数；%A约束方程的系数矩阵;AA最优表中约束方程的系数矩阵；%b约束方程右端列向量;输出b级变量的值；%c目标函数的系数。cb基变量的系数%输出xx为基变量的下标；%fm输出目标函数的值；%flg表示解得四种情况；B0=A(:,1:m);%B0初始可行基矩阵
(单位矩阵)；cb=c(:,1:m);xx=1:m;%xx变量的下标；
sgma=c-(cb*BO)*A;%sgma检验数；
vv=1;
h=-l;
AA=[];
sta=ones(m,l);
end
fori=1:n
end
ifsgma(i)>0
ifvv〜=1
h=1;
ww=0;
end
fori=m+1:n
end
ifsgma(i)<0
whileh>0
ww=ww+1;
[msg,mk]=max(sgma);
end
fori=1:m
end
ifA(i,mk)>0
ifww==n-m
sta(i)=b(i)/A(i,mk);
flg='Thereisonlyonesolution';
else
AA=A;
sta(i)=10000;
else
end
flg='Therearemanysolutions';
end
AA=A;
[mst,mr]=min(sta);
end
ifmst==10000
end
flg='unboundedsolution';
end
fm=inf;
end
xx=[];
end
b=[];
End