认识LP
线性规划(Linear Programming) 特指目标函数和约束条件皆为线性的最优化问题.
- 目标函数: 多个变量形成的函数
- 约束条件: 由多个等式/不等式形成的约束条件
- 线性规划: 在线性约束条件下,目标函数求极值的问题
- 可行解: 满足线性约束条件下的解
- 可行域:所有可行解构成的集合
- 最优解: 使目标函数取得极值的可行解
线性
个人觉得最好理解是用向量了. 就是元素满足 加法和数乘 的形式
- \(f(a+b) = f(a)+f(b)\)
- \(f(ca) = c f(a), c为常数\)
当然要理解上面两个等式可能需要去理解向量空间, 线性变换这些内容,嗯, 反正我自己已经懂了, 有时间可以分享.
定义模型的步骤
前提一定是线性的哈
- 确定决策变量
- 确定线性目标函数, 求max 或 min
- 确定线性约束条件
- 写出数学模型
case1
球队运作
需求:
需要补充7名球员, 每名球员有攻击值和防守值, 希望7名球员的进攻值大于500, 防守值大于400, 且要尽可能省钱.
战力值 | 防守值 | 价格(万) | |
进攻型 | 90 | 60 | 1000 |
平衡性 | 80 | 80 | 800 |
防守型 | 40 | 95 | 500 |
求解:
定义决策变量: 假设补充进攻型,防守型,平衡性各a,b,c名, 总价格为y万元,即
(s.t 即 subject to "受制于")
\(min \ y = 1000a + 800b + 500c\)
s.t.
\(90a + 80b +40c >= 500 \\ 60a+ 80b + 95c >= 400 \\a + b+ c = 7 \\ a,b,c >=0\)
case2
采购方案
需求:
作为采购经理,有2000元经费, 需采购单价为50元的桌子若干和单价20元的椅子若干.
- 桌椅总数尽可能多
- 椅子数量不少于桌子, 且不多于桌子的1.5倍
求解:
定义决策变量: 购买x1张桌子, x2把椅子, 总数为y.
\(max \ y = x1 + x2\)
s.t.
\(50x1 + 20x2 <= 2000 \\ 1.5x1 >= x2 \\ x1<= x2 \\ x1,x2 >=0\)
几何求解
- 约束条件的交集构成可行域
- 最优解即平行移动目标函数, 使其在可行域上达到截距最大
- 而最优解, 就是再交集域的顶点, 而无需在内部考虑.(不用严格用什么向量,几何, 什么定比分店,证明,肯定是边界的顶点上嘛)
证明最优解在边界顶点
假设平面三角形域顶点分别为x1, x2, x3, 最优解x0在三角形内, 过顶点x1,和x0的直线与底边 x2-x3交于点x4.
通过中学学的定点分比, 对x0作分解
$x_0 = $
LP的标准形式
\(min \ c^Tx \\ s.t. \ Ax <=b \\ x >= 0\)
- x, c, 表列向量, \(c^T\)是行向量, \(c^Tx\)表示线性组合
- \(Ax=b\)
- if 目标函数是求max, 则 - max 即转为了求min
- if 约束条件有 >= , 则 - (>=) 即转为了 <=
将case2 转为标准型
\(min \ y = -x1 + -x2\)
s.t.
\(50x1 + 20x2 <= 2000 \\ -1.5x1+x2 <= 0 \\ x1-x2<= 0 \\ x1,x2 >=0\)
将case2转为松弛型
- 松弛型: 用等式约束来等价表述不等式约束
- 松弛变量度量了等式约束与原不等式约束直接的松弛或差别
- 其实就是为了求解方便呗
$min y = -x1-x2 $
s.t.
\(50x1 + 20x2 + a1 = 2000 \\ -1.5x1+x2 + a2 =0 \\ x1-x2+a3=0 \\ x1,x2,a1,a2,a3 >=0 \\ 其中a1,a2,a3为松弛变量\)
(附) 证明最优解在边界顶点
假设平面三角形域顶点分别为x1, x2, x3, 最优解x0在三角形内, 过顶点x1,和x0的直线与底边 x2-x3交于点x4.
通过中学学的定比分点, 对x0作分解
\(x_0 = \lambda_1 x_1 + (1-\lambda_1) x_4, 其中\lambda_1 = ||x0-x1||/||x0-x4||\)
\(x4 = \lambda_2 x2 + (1- \lambda_2) x3, 其中\lambda_2 = c(定点分比值))\)
即:
\(x_0 = \lambda_1 x1 + (1-\lambda_1) \lambda_2 x2 + (1-\lambda_1)(1-\lambda_2) x3\)
其中 \(\lambda_1 + (1-\lambda_1) \lambda_2 + (1-\lambda_1)(1-\lambda_2) = 1\)
假设 \(c^tx1 > =c^tx2 >= c^tx3, 根据大前提即:\\ c^tx0 >= c^tx1 >= c^tx2 > =c^tx3\)
即:
$c^tx_0 =x_0 = \lambda_1 c^t x1 + (1-\lambda_1) \lambda_2 c^t x2 + (1-\lambda_1)(1-\lambda_2) c^t x3 $
\(>=c^tx_0 =x_0 = \lambda_1 c^t x1 + (1-\lambda_1) \lambda_2 c^t x1 + (1-\lambda_1)(1-\lambda_2) c^t x1\)
$ = c^tx1$
即说明最优解并不是 x0, 而是顶点x1, 不在内部哦.
耐心和恒心, 总会获得回报的.