线性目标规划(线性目标规划、图解法、单纯形法)

线性目标规划


请添加图片描述

数学模型
  • 普通线性规划的不足

    • 不能处理多目标的优化问题
    • 不允许约束资源有丝毫超差

请添加图片描述

  • 实际上工厂在作决策时,还要考虑市场等一系列其他条件:

    • 根据市场信息,产品Ⅰ的销售量有下降的趋势,故考虑产品Ⅰ的产量不大于产品Ⅱ
    • 超过计划供应的原材料时,需用高价采购,会使成本大幅度增加
    • 应尽可能充分利用设备台时,但不希望加班
    • 应尽可能达到并超过计划利润指标56元

    因此考虑产品决策时,为多目标决策问题。在原材料供应受严格限制的基础上考虑:首先是产品Ⅱ的产量不低于产品I的产量;其次是充分利用设备有效台时,不加班;再次是利润额不小于56元。求决策方案。

  • 不考虑市场等其他条件时,模型及最优解
    max ⁡ z = 8 x 1 + 10 x 2 { 2 x 1 + x 2 ≤ 56 x 1 + 2 x 2 ≤ 70 x 1 , x 2 ≥ 0 ⟹ 最优解 x ∗ = ( 4 , 3 ) T ⟹ 最优值 z ∗ = 62 \begin{array}{l}\max z=8x_1+10x_2 \\ \left\{\begin{array}{c} 2x_1+x_2\le 56\\ x_1+2x_2\le 70\\ x_1,x_2\ge 0 \end{array}\right.\end{array} \begin{array}{l} \overset{\text{最优解}}{\Longrightarrow} x^*=(4,3)^T\\ \overset{\text{最优值}}{\Longrightarrow} z^*=62 \end{array} maxz=8x1+10x22x1+x256x1+2x270x1,x20最优解x=(4,3)T最优值z=62

  • 目标规划数学模型概念

    • 正负偏差变量 d + ,   d − d^+,~d^- d+, d
      • d + d^+ d+:决策值超过目标值的部分
      • d − d^- d:决策值未达到目标值的部分
      • 决策值不可能同时超过目标值又没有达到目标值,即 d + ⋅ d − = 0 d^+\cdot d^-=0 d+d=0
    • 绝对(硬)约束和目标(软)约束
      • 绝对约束:红线,底线,严格遵守,不可逾越
      • 目标约束:可在小范围内调节(具体问题具体分析,但总归不大,类似于“邻域”的感觉)
    • 优先因子:决策者在要求达到各个目标时,有主次轻重的不同。第一位达到的目标赋予优先因子 P 1 P_1 P1,次位的目标赋予优先因子 P 2 P_2 P2 ⋯ \cdots ,且满足 P k ≫ P k + 1 , k = 1 , 2 , ⋯   , K P_k\gg P_{k+1},k=1,2,\cdots,K PkPk+1,k=1,2,,K
    • 权系数:若要区别相同优先因子的两个目标的差别,可以分别赋予不同的权系数 ω j \omega_{j} ωj
  • 考虑多目标规划

    • 约束条件

      1. 看绝对约束: 2 x 1 + x 2 ≤ 11 2x_1+x_2\le 11 2x1+x211,这部分不进行处理

      2. 对目标约束按优先级从高到低一次判断,按正常情况列出方程,一般视角下建模
        P 1 : x 1 ≤ x 2 P 2 : x 1 + 2 x 2 ≤ 10 P 3 : 8 x 1 + 10 x 2 ≥ 56 \begin{array}{cl} &P_1:x_1\le x_2\\ &P_2:x_1+2x_2\le 10\\ &P_3:8x_1+10x_2\ge 56 \end{array} P1:x1x2P2:x1+2x210P3:8x1+10x256

      3. 转化为目标约束:变量部分+偏差=常数部分
        KaTeX parse error: Expected '}', got '&' at position 17: … \begin{array}{&̲&} &(x_1-x_2)+…

      4. 最终约束条件
        { 2 x 1 + x 2 ≤ 11 x 1 − x 2 + d 1 − − d 1 + = 0 x 1 + 2 x 2 + d 2 − − d 2 + = 10 8 x 1 + 10 x 2 + d 3 − − d 3 + = 56 x 1 , x 2 ≥ 0 , d i − ≥ 0 , d i + ≥ 0 \left\{\begin{array}{l} 2x_1+x_2\le 11\\ x_1-x_2+d_1^--d_1^+=0\\ x_1+2x_2+d_2^--d_2^+=10\\ 8x_1+10x_2+d_3^--d_3^+=56\\ x_1,x_2\ge 0,d_i^-\ge 0,d_i^+\ge 0 \end{array}\right. 2x1+x211x1x2+d1d1+=0x1+2x2+d2d2+=108x1+10x2+d3d3+=56x1,x20,di0,di+0

    • 目标函数,三种基本形式 min ⁡ z = f ( d + , d − ) \min z=f(d^+,d^-) minz=f(d+,d)

      • 要求恰好达到目标值,正负偏差变量都要尽可能小
        min ⁡ z = f ( d + + d − ) \min z=f(d^++d^-) minz=f(d++d)

      • 要求不超过目标值,即允许达不到目标值,就是正偏差尽可能小
        min ⁡ z = f ( d + ) \min z=f(d^+) minz=f(d+)

      • 要求超过目标值,即超过量不限,但是负偏差变量要尽可能小
        min ⁡ z = f ( d − ) \min z=f(d^-) minz=f(d)

      • 原问题
        min ⁡ z = P 1 d 1 + + P 2 ( d 2 − + d 2 + ) + P 3 d 3 − \min z=P_1d_1^++P_2(d_2^-+d_2^+)+P_3d_3^- minz=P1d1++P2(d2+d2+)+P3d3
        注: P 1 P_1 P1:产品Ⅱ的产量不低于产品I的产量

        注: P 2 P_2 P2:其次是充分利用设备有效台时,不加班

        注: P 3 P_3 P3:再次是利润额不小于56元

  • 一般数学模型
    min ⁡ z = ∑ l = 1 L P l ∑ k = 1 K ( ω l k − d k − + ω l k + d k + ) { ∑ j = 1 n c k j x j + d k − − d k + = g k , k = 1 , ⋯   , K ∑ j = 1 n a i j x j ≤ ( = , ≥ ) b i , i = 1 , ⋯   , m x j ≥ 0 , j = 1 , ⋯   , n d k − , d k + ≥ 0 , k = 1 , ⋯   , K \begin{array}{l} \min z=\sum\limits_{l=1}^LP_l\sum\limits_{k=1}^K(\omega_{lk}^-d_k^-+\omega_{lk}^+d_k^+)\\ \left\{\begin{array}{l} \sum\limits_{j=1}^nc_{kj}x_j+d_k^--d_k^+=g_k,k=1,\cdots,K\\ \sum\limits_{j=1}^na_{ij}x_j\le(=,\ge)b_i,i=1,\cdots,m\\ x_j\ge 0,j=1,\cdots,n\\ d_k^-,d_k^+\ge 0,k=1,\cdots,K \end{array}\right. \end{array} minz=l=1LPlk=1K(ωlkdk+ωlk+dk+)j=1nckjxj+dkdk+=gk,k=1,,Kj=1naijxj(=,)bi,i=1,,mxj0,j=1,,ndk,dk+0,k=1,,K

图解法
  • 主要针对二维情况

  • 画图原则

    • 按优先级从高到低画(含权系数)
    • 画一条线判断一条线

    请添加图片描述

    注意思考对问题进行建模的过程,参考视频

    请添加图片描述

    一定要把所有偏差变量讨论完

解目标规划的单纯形法
  • 求出的解为满意解。计算方法与普通单纯形法和大M法完全一致,只需要满足 P 1 ≫ P 2 ≫ ⋯ P_1\gg P_2\gg \cdots P1P2即可

  • e.g. 对于 min ⁡ z = P 1 d 1 + + P 2 ( d 2 − + d 2 + ) + P 3 d 3 − , s . t . { 2 x 1 + x 2 ≤ 11 x 1 − x 2 + d 1 − − d 1 + = 0 x 1 + 2 x 2 + d 2 − − d 2 + = 10 8 x 1 + 10 x 2 + d 3 − − d 3 + = 56 x 1 , x 2 ≥ 0 , d i − ≥ 0 , d i + ≥ 0 \min z=P_1d_1^++P_2(d_2^-+d_2^+)+P_3d_3^-,s.t.\left\{\begin{array}{l} 2x_1+x_2\le 11\\ x_1-x_2+d_1^--d_1^+=0\\ x_1+2x_2+d_2^--d_2^+=10\\ 8x_1+10x_2+d_3^--d_3^+=56\\ x_1,x_2\ge 0,d_i^-\ge 0,d_i^+\ge 0 \end{array}\right. minz=P1d1++P2(d2+d2+)+P3d3,s.t.2x1+x211x1x2+d1d1+=0x1+2x2+d2d2+=108x1+10x2+d3d3+=56x1,x20,di0,di+0

    1. 化标准型
      { 2 x 1 + x 2 + x 3 = 11 x 1 − x 2 + d 1 − − d 1 + = 0 x 1 + 2 x 2 + d 2 − − d 2 + = 10 8 x 1 + 10 x 2 + d 3 − − d 3 + = 56 x 1 , x 2 ≥ 0 , d i − ≥ 0 , d i + ≥ 0 \left\{\begin{array}{l} 2x_1+x_2+x_3=11\\ x_1-x_2+d_1^--d_1^+=0\\ x_1+2x_2+d_2^--d_2^+=10\\ 8x_1+10x_2+d_3^--d_3^+=56\\ x_1,x_2\ge 0,d_i^-\ge 0,d_i^+\ge 0 \end{array}\right. 2x1+x2+x3=11x1x2+d1d1+=0x1+2x2+d2d2+=108x1+10x2+d3d3+=56x1,x20,di0,di+0

    2. 画单纯形表

请添加图片描述
请添加图片描述
请添加图片描述
请添加图片描述

注:在表3中非基变量 d 3 + d_3^+ d3+的 检验数为0,表示存在多重满意解,将其作为换入变量可以计算 θ \theta θ得到 d 1 − d_1^- d1为换出变量,从而求得两个满意解 x ∗ = ( 2 , 4 ) T x^*=(2,4)^T x=(2,4)T x ∗ = ( 10 / 3 , 10 / 3 ) T x^*=(10/3,10/3)^T x=(10/3,10/3)T。此外,这两点的凸线性组合都是满意解

  • 3
    点赞
  • 22
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值