附录1：线规划的阅读材料.docVIP

附录1 线性规划的阅读材料 一、问题的提出 首先看几个例子。 例 1 某工厂生产A,B,C三种产品，每件产品所消耗的材料、工时以及能够获取的利润如下表所示。 产品 利润（元/件） 材料（kg/件） 工时（h/件） A 7 4 4 B 3 4 2 C 6 5 3 已知该厂每天的材料消耗不得超过600kg，工时不得超过1400h。问每天生产A,B,C三种产品各多少件可使利润最大？ 解 设 分别是 A,B,C 三种产品每天的生产量，由于材料的限制，三种产品的材料总消耗量不能超过600kg，这可用不等式表示为 同样，由工时限制条件，可得不等式 该厂的目标是使利润最大，用表示该厂每天的利润数，则 由于产品的生产量不可能是负的，所以还有不等式 综上所述，就得到这个问题的数学模型 max s.t 例 2 已知A,B两地生产某种商品，其产量分别不多于300t和600t，运往销地甲、乙、丙三处，销量分别不得少于100t、300t和400t。从产地到销地的每吨商品运输费用如下表所示。 产 地 销 地 甲 乙 丙 A 2 4 3 B 5 3 4 现在希望制订一种商品生产和调运方案，使总运费最小。 解 设从产地A运到甲、乙、丙三地的商品分别是和吨，从产地B运到甲、乙、丙三地的商品分别是和吨。 由产地A运到甲、乙、丙三地的商品总量即为A地的商品产量，应该不多于300t，于是得到不等式 同理，由产地B运到甲、乙、丙三地的商品总量即为B地的商品产量，应该不多于600t，即 甲地的商品销量就是A,B两地运到甲地的商品总量，应该不少于100t，则有 同理，对于乙地和丙地，分别有 目标是使运费最小，设总运费是，则 此外，由于所有的变量显然不能为负数，所以还要加上非负性的约束条件，于是可得此问题的数学模型 min s.t 以上两例中数学模型的共同特点是，求一组变量的值，这组变量要满足一定的约束条件，此约束条件可以用一组线性不等式或线性方程组来表示。同时还要使某个指标达到最小或最大，而这一指标也是这组变量的线性函数。具有这种特点的问题，称为线性规划问题。 二、线性规划问题的数学模型 线性规划问题的数学模型如下 min cx (1) s.t Ax ? b 或 min cx s.t Ax ? b (2) vl ? x ? vu 其中x为维未知向量，c为目标系数，约束系数A为一个的矩阵，b为维列向量。（2）式中的vl与vu都是维常数向量。 三、线性规划问题的解法 使用MATLAB中的函数LP可求解线性规划问题，其一般形式如下： 调用格式1 [x,lambda,how]=LP(c,A,b) 功能 求解线性规划问题（1）。返回值中，x为最优解，lambda为相应的最优Lagrange乘子，how为一个文本串（如’infeasible’），用来描述求解状况。 调用格式2 x=LP(c,A,b,vl,vu) 功能 求解线性规划问题（2）。 调用格式3 x=LP(c,A,b,vl,vu,x0) 功能 求解线性规划问题（2）,其中x0为初始点。 调用格式4 x=LP(c,A,b,vl,vu,x0,) 功能 求解线性规划问题（2）,参数指出，由A和b所确定的约 束条件中前个是等式。 调用格式5 x=LP(c,A,b,vl,vu,x0,,-1) 功能 求解线性规划问题（2）,其中最后一个参数-1表明，当目标函数无界或者无可行解时，执行过程本应给出的警告信息被关掉。 例 3 求解线性规划问题 min s.t 解 在MATLAB命令窗口中操作如下： c=[2；2]； A=[-1 1；-1 2]； b=[-1；0]； x=lp（c，A，b，[0 0]，[]） x= 1 0 例 4 求解线性规划问题 min s.t 解 在MATLAB命令窗口中操作如下： c=[2；2]； A=[1 -1；1 2]； b=[-1；-2]； x=lp（c，A，b，[0 0]，[]） Warning：The constraints are overly stringent; There is no feas