1、目录一摘要 2二实验目的 2三实验内容 2四建立数学模型 3五实验原理 5六MALTAB 程序代码及注释 7七结果运行测试 13八心得与感悟 15 2一摘要：线性规划是运筹学中研究较早、发展较快、应用广泛、方法较成熟的一个重要分支,它是辅助人们进行科学管理的一种数学方法.研究线性约束条件下线性目标函数的极值问题的数学理论和方法,英文缩写 LP。自 1946 年 G.B.Dantizig 提出单纯形法以来，它一直是求解线性规划问题的最有效的数学方法之一。单纯形法的理论根据是：线性规划问题的可行域是 n 维向量空间 Rn 中的多面凸集,其最优值如果存在必在该凸集的某顶点处达到。顶点所对应的可行解称
2、为基本可行解。通过引入普通单纯形法，依次迭代并判断，逐步逼近，最后得到最优解。若不是，则按照一定法则转换到另一改进的基本可行解，再鉴别；若仍不是，则再转换，按此重复进行。因基本可行解的个数有限，故经有限次转换必能得出问题的最优解。如果问题无最优解也可用此法判别。关键字：线性规划，单纯形法，最优值，最优解二实验目的：1.加强学生分析问题能力，锻炼数学建模能力。2.了解并掌握 MATLAB 软件中的线性规划问题的编程、求解和分析。3.利用所学的 MALTAB 语言，完成对单纯形法问题的编程设计。三实验内容：某商场决定，营业员每周连续工作 5 天后连续休息 2 天，轮流休息，据统计，商场每天需要营业
3、员如下：星期一：300，二：300；三：350，四：400，五：480，六：600；日：500；（1）商场人力资源部应如何安排每天上班的人数才能使商场总的营业员最少（2）若商场可以雇佣临时工，上班时间同正式工，若正式工每天工资80，临时工每天 100，问商场是否应雇佣临时工及雇佣多少名？3四建立数学模型：从实际问题中建立数学模型一般有以下三个步骤： 1.根据影响所要达到目的的因素找到决策变量； 2.由决策变量和所在达到目的之间的函数关系确定目标函数； 3.由决策变量所受的限制条件确定决策变量所要满足的约束条件。 当我们得到的数学模型的目标函数为线性函数，约束条件为线性等式或不等式时称此数学模型
4、为线性规划模型。 线性规划问题的标准形式：由题可知，可设每天上班人数分别应为 x1,x2,x3,x4,x5,x6,x7;建立下列数学模型4将其转化为标准形式为：即价值向量约束矩阵右端向量5五实验原理：根据单纯形法的原理，在线性规划问题中，决策变量（控制变量）x1，x2，x n 的值称为一个解,满足所有的约束条件的解称为可行解。使目标函数达到最大值（或最小值）的可行解称为最优解。这样，一个或多个最优解能在整个由约束条件所确定的可行区域内使目标函数达到最大值（或最小值） 。求解线性规划问题的目的就是要找出最优解。最优解可能出现下列情况之一：存在着一个最优解；存在着无穷多个最优解；不存在最优解，这只
5、在三种情况下发生，即没有可行解或各项约束条件不阻止目标函数的值无限增大（或向负的方向无限增大） 。单纯形法的一般解题步骤可归纳如下： 把线性规划问题的约束方程组表达成典范型方程组，找出基本可行解作为初始基本可行解。 若基本可行解不存在，即约束条件有矛盾，则问题无解。 若基本可行解存在，从初始基本可行解作为起点，根据最优性条件和可行性条件，引入非基变量取代某一基变量，找出目标函数值更优的另一基本可行解。 按步骤 3 进行迭代,直到对应检验数满足最优性条件（这时目标函数值不能再改善） ，即得到问题的最优解。 若迭代过程中发现问题的目标函数值无界，则终止迭代。6流程图如下：対原问题增加 m 个人工变
6、量构造辅助问题判断辅助问题最优值 g=0无解，停止人工变量 xj 是否为非基变量把人工变量对应的列从单纯形表中去掉进行换基迭代得到新矩阵 BB 中是否有人工变量得到初始可行基 B求对应典式和检验数 判断 k0得到最优解进行换基迭代得到新基判断 Ak0问题无界否是否否是是否是7六MALTAB 程序代码及注释：function x,minf,flag,cpt=dcxsf(A,b,c)format rat %使数据可以以分数形式输出c=-c; %将目标函数系数向量加负号得到单纯形表第0 行m,n=size(A); %求约束矩阵的行数和列数m1=m; %保存下原来的行数s=eye(m); %生成秩为m
7、的单位矩阵A=A s; %将s矩阵添加到A矩阵右侧A=A b; %将b向量添加到A矩阵右侧g1=zeros(1,n); %生成一个1行n 列的零矩阵g1x=ones(1,m); %生成一个 1行m列元素全为1的矩阵g1=g1 -x; %将g1和-x合并，产生一个新的前n列为0 后m列为-1的单行矩阵g=0; %初始化一个单元素零矩阵g1=g1 g;%将单元素零矩阵添加到g1 右侧，生成人工向量的检验向量g1s1=n+m+1; %记录目前列数之和s2=zeros(1,m+1);%生成1行m+1列的零矩阵s2c=c s2;%将s2添加到c右侧A1=zeros(m,1);%生成一个m 行1列的零矩阵
8、A1for i1=1:m A1(i1,1)=i1+n;%基变量的数值存储区endfor i=1:mg1(1,:)= g1(1,:)+A(i,:); enddecide=find(g1(1,1:m+n)0); %寻找 g1 中大于零的数值列数存于decide 数组中while isempty(decide) %当decide不为空i=decide(1); %将列数赋给i text=find(A(1:m,i)0); %text存放该列中所有数值大于零的行数if isempty(text) %如果text为空则无解 flag=0;break;endmin=inf;%min初始化为无穷大for i1=
9、1:m % 当该列存在大于零的数值时 if A(i1,i)0 % 再进行查找endif g1(1,s1)0 % 无解情况flag=-1; endif g1(1,s1)=0 %置空矩阵g1(1,:)=;for i6=1:mA(:,n+1)=;endfor i6=1:mc(n+1)=;enddecide=find(A1(1:m,1)n); % 查找是否有人工变量if isempty(decide) while isempty(decide)x1=decide(1); % 存放的是人工变量的行数 text=find(A(x1,1:n)=0); % 该行的所有元素都不为零的列坐标if isempty(
10、text) %如果text为空decide(1)=;A1(x1,:)=;A(x1,:)=;m=m-1;elsei=text(1); % i 保存列数 A(x1,:)=A(x1,:)/A(x1,i);%进行单位化A1(x1,1)=i;c(1,:)=c(1,:)+(-1*c(1,i)*A(x1,:);for i1=1:mif i1=x1A(i1,:)=A(i1,:)+(-1*A(i1,i)*A(x1,:);%进行换基迭代9endenddecide(1)=;%赋值为空text(1)=;endendenddecide=find(c(1,1:n)0); %decide存放c中该行中所有数值大于零的列数
11、while isempty(decide)i=decide(1); % 将列数赋给itext=find(A(:,i)0); % 保存大于零的项的行数if isempty(text) % 为空的时候，则无解flag=0; %有可行解但无最优值，flag记为0c=c;A; %将c添加到A矩阵上方的到对应解的单纯形表cpt=c;x=zeros(n,1);for o=1:nfor o1=1:mif o=A1(o1,1)x(o,1)=A(o1,n+1);endendenddisp(该问题有可行解但没有最优解!)minf=-inf;%无最优值，将minf赋值为无穷小x=;%解向量为空return;endm
12、in=inf; %min初始化为无穷大for i1=1:mif A(i1,i)0 % 再进行查找endc=c;A;%得到对应解的单纯形表cpt=c;x=zeros(n,1);for o=1:nfor o1=1:mif o=A1(o1,1)x(o,1)=A(o1,n+1);endendendminf=c(1,n+1);%得到最优值flag=1 ; %有最优值 记录flag 为1end if flag=0 %有可行解但没有最优值情况disp(该问题有可行解但没有最优解!)x=; %赋值为空minf=-inf;%最优值输出为无穷小endif flag=1 %有可行解有最优值情况disp(该问题存在最
13、优解!)x=x; minf=minf;end;if flag=-1 %没有可行解情况disp(该问题没有可行解!)x=;%赋值为空cpt=;minf=;endend11七结果运行测试：输入测试数据： A=1 1 1 1 1 0 0 -1 0 0 0 0 0 0;0 1 1 1 1 1 0 0 -1 0 0 0 0 0;0 0 1 1 1 1 1 0 0 -1 0 0 0 0;1 0 0 1 1 1 1 0 0 0 -1 0 0 0;1 1 0 0 1 1 1 0 0 0 0 -1 0 0;1 1 1 0 0 1 1 0 0 0 0 0 -1 0;1 1 1 1 0 0 1 0 0 0 0 0
14、0 -1; b=300 300 350 400 480 600 550; c=1 1 1 1 1 1 1 0 0 0 0 0 0 0; x,minf,flag,cpt=dcxsf(A,b,c)运行结果为：该问题存在最优解!x =170 290/3 120 50/3 0 200/3 440/3 310/3 0 0 0 0 0 0 minf =1850/3 flag =1 cpt =Columns 1 through 6 120 0 0 0 -1/3 0 1 0 0 0 -1 0 0 0 1 0 -1 0 0 0 0 0 2/3 0 0 0 0 0 2/3 1 0 1 0 0 2/3 0 0 0
15、0 0 -5/3 0 0 0 0 1 2/3 0 Columns 7 through 12 0 0 -1/3 0 -1/3 0 0 0 0 1 -1 1 0 0 0 0 0 1 1 0 2/3 -1 2/3 -1 0 0 -1/3 0 -1/3 0 0 0 -1/3 0 2/3 -1 0 1 -2/3 1 -2/3 1 0 0 -1/3 0 -1/3 0 Columns 13 through 15 -1/3 -1/3 1850/3 -1 0 170 -1 0 120 2/3 -1/3 440/3 -1/3 2/3 200/3 2/3 -1/3 290/3 -2/3 -2/3 310/3 2/3
16、 -1/3 50/3 解得最优解为(170 97 120 17 0 67 146)最优值为 61713八心得与感悟：通过本学期学习 MALTAB 软件，初步了解并掌握了 MATLAB 软件中的线性规划问题的编程求解的方法。通过对编程问题的求解，我明白了掌握编程的重要性。我们可以把单纯形法的知识运用到企业中去，运用线性规划并配合计算机进行测算非常简便易行，几分钟就可以拿出最优方案，提高了企业决策的科学性和可靠性。最开始学习的时候，我真的是一片茫然，完全不知道从何处入手，如何把问题和代码联系到一起，后来经过不断地探索，慢慢的掌握 maltab 软件的应用。期间也遇到了许许多多的问题和错误，通过和同学讨论，上网百度，一点一点将错误改正，最终完成了代码。我也知道干什么都不是一件容易的事，坚持下来总会有收获的。
