提出背景:
对于一般线性规划,AAA中未必刚好有一个mmm阶单位矩阵,因此没有现成的初始基可行解
两阶段方法
引入人工变量,构造辅助线性规划
在上述辅助线性规划执行单纯形法,知道找到辅助线性规划的最优解g∗g^*g∗
①若g∗>0g^*>0g∗>0,则原线性规划无解
②若g∗=0g^*=0g∗=0
- 基变量全在xix_ixi中,则基可行解就是原规划的初始基可行解
- 基变量不全在xix_ixi中,如yky_kyk仍是基变量,则检查单纯形表的第kkk行,若对应原线性规划的半段系数全是0,则可以删掉第kkk行,否则进行换基操作(无论正负,b值应该为0才对)
上述过程完成了第一阶段任务,第二阶段按照上述的初始基可行解,重新计算单纯形表格的判别数和函数值,然后执行单纯形法。
大M法
在约束中增加人工变量Xa=(y1,…,ym)TX_a=(y_1,\dots,y_m)^TXa=(y1,…,ym)T,同时在目标函数中加上惩罚项MeTXaMe^TX_aMeTXa,如下:
用单纯形法求解上述线性规划
- 上述线性规划有最优解,且Xa∗=0X_a^* = 0Xa∗=0时,原线性规划有最优解为X∗X^*X∗
- 上述线性规划有最优解,且Xa∗≠0X_a^* \neq 0Xa∗=0时,原线性规划无最优解