实验2 单纯形法求解线性规划.docVIP

实验2 单纯形法求解线性规划 成绩 专业班级 信息123班 学号 201212030329 姓名 于晓敏 报告日期 . 实验类型：●验证性实验 ○综合性实验 ○设计性实验 实验目的：进一步熟练掌握单纯形法求解线性规划。 实验内容：单纯形法求解线性规划4个（题目自选） 实验原理 首先要找到一个初始基本可行解，求出对应的检验数，判断其是否是最优解，如果是就停止计算；否则，就进行迭代找到另一个能使得目标函数值更优的基本可行解，然后再判断其是否是最优解，如此反复进行下去，直到找到最优解或者判断线性规划问题无解为止（线性规划解有四种情形，唯一最优解，无穷多个最解，无界解，无可行解）。 实验步骤 1 要求上机实验前先编写出程序代码 2 编辑录入程序 3 调试程序并记录调试过程中出现的问题及修改程序的过程 4 经反复调试后，运行程序并验证程序运行是否正确。 5 记录运行时的输入和输出。 预习编写程序代码： 实验报告：根据实验情况和结果撰写并递交实验报告。 实验总结： 参考程序 线性规划1 方程 线性规划2 方程： 线性规划3 方程： 线性规划4 方程： 单纯形法matlab程序-ssimplex 求解标准型线性规划: A是单纯初始表，包括最后一行是初始的检验数，最后一列是资源向量b N是初始的基变量的下标 输出变量sol是最优解，其中松弛变量（或剩余变量）可能不为0 输出变量val是最优目标值，kk是迭代次数 一、建立M文件ssimplex function [sol,val,kk]=ssimplex(A,N) [mA,nA]=size(A); kk=0; flag=1; while flag kk=kk+1; if A(mA,:)=0 flag=0; sol=zeros(1,nA-1); for i=1:mA-1 sol(N(i))=A(i,nA); end val=-A(mA,nA); else for i=1:nA-1 if A(mA,i)0A(1:mA-1,i)=0 disp(have infinite solution!); flag=0; break; end end if flag temp=0; for i=1:nA-1 if A(mA,i)temp temp=A(mA,i); inb=i; end end sita=zeros(1,mA-1); for i=1:mA-1 if A(i,inb)0 sita(i)=A(i,nA)/A(i,inb); end end temp=inf; for i=1:mA-1 if sita(i)0sita(i)temp temp=sita(i); outb=i; end end for i=1:mA-1 if i==outb N(i)=inb; end end A(outb,:)=A(outb,:)/A(outb,inb); for i=1:mA if i~=outb A(i,:)=A(i,:)-A(outb,:)*A(i,inb); end end end end end 一、规划输入：A=[1 2 1 0 0 8;4 0 0 1 0 16;0 4 0 0 1 12;2 3 0 0 0 0] A = 1 2 1 0 0 8 4 0 0 1 0 16 0 4 0 0 1 12 2 3 0 0 0 0 初始的基变量的下标N N=[3 4 5] N = 3 4 5 运行程序： [sol,val,kk]=ssimplex(A,N) sol = 4 2 0 0 4 val = 14 kk = 4 二、线性规划2 加入松驰变量，化为标准型，得到单纯初始表A, A=[5 10 1 0 0 -50;2 1 0 1 0 1;1 4 0 0 1 -4;1 3 0 0 0 0] A = 5 10 1 0 0 -50 2 1 0 1