Excel-单纯形法（大M法）求解 直接求解与规划求解功能

单纯形法

线性规划的标准形式为：max CX

其中，B >= 0，矩阵C、X、A、B的阶分别为 1 * n、 n * 1 、m * n 、m * 1 。
单纯形的一般形式如表1所示：

其中, E≥0,(D E)是由(1)式中(A B)经初等行变换所得,(D E)中有m-r(A)行全为0(否则,(1)式没有可行解;通常,删除单纯形中全为0的行),D中基变量所在列组成r(A)阶单位子矩阵,这个单位子矩阵称作基,同一基变量所在行列的交叉处为1。通常,单纯形法由两个阶段组成:构造初始单纯形和换基。构造初始单纯形的一般步骤为:对(1)式中(AB)做初等行变换,当变换所得(D E)中E≥0且D有r( A)阶单位子矩阵时,按表1相应填写初始单纯形中的X"、C、D、E、F、G、H、T,其中,T=C-FD。按_上述步骤构造初始单纯形有一定 难度,为此，对单纯形法做相应改进的算法应运而生,下述大M法就是其中之一。

大M法

用单纯形法对下述(2)式求解称作大M单纯形
法,简称大M法:

其中,Y=(y1y2… ym)T是人工变量,N=(-M -M… -M),I是m阶单位矩阵，M充分大,其它部分与(1)式相同。大M法依据(2)式和(1)式的下述关系，把求解(1)式的问题转换为求解(2)式:①Xo是(1)式最优解的充分必要条件是:(X。0)是(2)式最优解。②如果(2)式没有最优解,则(])式也没有最优解。③如果(2)式存在形如(X。Yg)的最优解且Y≠0,则(1)式没有最优解。

换基由四部分组成:检验最优解选择人基变量、选择出基变量、换基。
最优解检验准则:当T1和T0中没有正数且基中人工变量均为0(所在行的约束值均为0)时,当前单纯形的可行基解是最优解(之一)。
入基准则:当T1 或T0中有正数时,选择T0中最大数(不要求为正数)所在列入基,或者说,所在列的变量入基。
出基准则:当变量x;入基时,如果x;的约束系数中没有正数,说明(1)式无上界,没有最优解。否则,x,入基到{h, =e/d,ld; >0,i=1 ~m}中最小数所在行,该行上当前基变量出基。
换基:按照新基对应填写新的单纯形。当变量xj入基到第i行时,G、F、F。的第i行相应改为x;和.xj的目标系数,对(D E)做初等行变换，使D的第j列中第i行为1、其它行为0,计算新的检验数。

例子

直接计算
1、标准化得到

2、在E1:J1、E2:J3、C4:C6、A4:B6、B4: B6、D4 ;D6、E4 :J6中依次填写约束变量、目标系数、基变量基变量的目标系数、约束值、约束矩阵。
A7填写目标值中M的系数,公式为= SUMPRODUCT( A4 :A6 , D4:D6)
B7填写目标值中常数项,公式为= SUMPRODUCT( B4 :B6 ,D4:D6)
E7填写检验数中M的系数,公式为= E$2 - SUMPRODUCT( $ A4: $ A6,E4:E6 )如图2所示,设置E7的条件格式,条件公式为=AND(E7 >0,E7 = $ K7)具体格式自定,醒目即可。把E7向右一直复制到J7。
E8填写检验数中常数项,公式为=IF(E7= $ K7,E $3 - SUMPRODUCT( $ B4:$ B6,E4:E6),"")设置E8的条件格式,条件公式为= AND(E8>0,E8= $ K8)把E8向右一直复制到J8。
K7填写检验数中M的最大系数,公式为= MAX( E7:J7)把K7复制到K8。
3、检验最优解。 检验数E7:J7中有正数,需要换基。
4、选择人基变量。 E8:J8 中只有G8这一个检验数,因此,x3人基。
5、选择出基变量。 K4 中填写公式= IF(G4 >0,D4/G4,"")并向下一直复制到K6。出基值K4: K6中最小正数在K4和K5,不妨选x|出基。
6、换基。 把A4:K8复制到A9:K13, C9、A9、B9中相应填写人基变量x及其目标系数。D9:D11中依次填人公式=D4、=D5- D9、=D6+ D9并向右一直复制到J列。然后,转回第3步。重复继续类似做。

最后得到的结果如下所示：
最终结果为x1=0，x2=0，x3=2，max=2

规划求解功能
先根据原式，填好系数、约束条件等。D5-D6的右边是等下要求的最优解具体值。

E2=@MMULT(A2:C2,E5:E7)，然后往下拉至E4。
E1=@MMULT(A5:C5,E5:E7) ,这是max。
填完后可能是这样的。

找到Excel的数据选项，最后边的规划求解功能。
如果你的excel版本是64位的话，可能会有一些小问题。建议使用32位的求解。

按下图选择，设置目标填的是E1；通过更改可变单元格是E5-E7。遵守约束条件点击右边的添加，逐个添加上。

选择求解方法选为：单纯线性规划。点击他右边的选项后，得到如下图所示。

然后确定，求解完成后得到如下结果。
x1=0，x2=0，x3=2，max=2

可根据需要选择，是否生成运算报告。
三个报告如下图所示：
运算结果报告：

敏感性报告：

极限性报告：

参考教程：兵团教育学院学报 2009年 第4期