QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3041|回复: 1
打印 上一主题 下一主题

线性规划(一):运输问题 (产销平衡) & 指派问题

[复制链接]
字体大小: 正常 放大

907

主题

66

听众

17万

积分

  • TA的每日心情
    开心
    2023-3-15 17:49
  • 签到天数: 224 天

    [LV.7]常住居民III

    社区QQ达人 邮箱绑定达人 元老勋章 发帖功臣 新人进步奖 优秀斑竹奖 金点子奖 原创写作奖 最具活力勋章 助人为乐奖 风雨历程奖

    跳转到指定楼层
    1#
    发表于 2021-10-29 11:50 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    1.可以转化为线性规划的问题

    很多看起来不是线性规划的问题也可以通过变换变成线性规划的问题来解决。下面举几个例子:

    2. 运输问题(产销平衡)  :康—希表上作业法

    8 Q4 [- Q+ v8 @6 u4 B0 k
    3. 指派问题的数学模型

    - _7 {# G- Y$ m9 ]  f+ G
    上述指派问题的可行解可以用一个矩阵表示,其每行每列均有且只有一个元素为 1,其余元素均为 0;可以用 1,...,n 中的一个置换表示。 问题中的变量只能取 0 或 1,从而是一个 0-1 规划问题。一般的 0-1 规划问题求解 极为困难。但指派问题并不难解,其约束方程组的系数矩阵十分特殊(被称为全单位模矩阵,其各阶非零子式均为正负1 ),其非负可行解的分量只能取0或1,故约束 =0或1 可改写为 而不改变其解。此时指派问题被转化为一个特殊的运输问题,其中m=n, & n1 N! J$ D6 a1 S% D" w' V- F$ r, L
    求解指派问题的匈牙利算法

    有时问题会稍复杂一些。

    例 9 求解系数矩阵C 的指派问题

    1 n& I& v6 Q. ]- Z3 A% d& \

    " Z: g0 }4 ?- x" K# ~! Z& t: Q

    : C2 g" r- T0 z- s" i
    4. 指派问题的计算机求解整数规划问题的求解可以使用 Lingo 等专用软件。对于一般的整数规划问题,无法 直接利用 Matlab 的函数,必须利用 Matlab 编程实现分枝定界解法和割平面解法。但对 于指派问题等 0−1整数规划问题,可以直接利用 Matlab 的函数 bintprog 进行求解。
    ' T0 Z  \9 g! z

    * s3 y4 }$ E) l2 X$ p9 R) T
    解:编写 Matlab 程序如下: 求得最优值为 21,最优指派方案为 .

    0 W" p5 @+ N" {6 f# Jc=[3 8 2 10 3;8 7 2 9 7;6 4 2 7 5    * ^& @- z- V+ ^! o, {
       8 4 2 3 5;9 10 6 9 10]; ; b3 t$ _9 ?, [' ?: O, C4 J. U
    c=c(;
      ~; Z) o( ]6 M  J7 p1 x6 Va=zeros(10,25);
    5 s+ A; c  O) z: [" Efor i=1:5   
    ' _. X* B1 [3 H2 R- j" }4 m    a(i,(i-1)*5+1:5*i)=1;    * s8 z% J9 E. P) w& Q/ F' x# ]
        a(5+i,i:5:25)=1;
    ( l1 ^+ ^+ M- Z) [- G: Vend
      r# A4 u& H% e, g  ~b=ones(10,1); 2 A7 s6 ^2 e" u
    [x,y]=bintprog(c,[],[],a,b);
    3 V2 a8 T7 T' K' Kx=reshape(x,[5,5]),y ; y& {& H7 A3 q7 ^; S

      @0 C- C* r5 _& z$ |8 A9 }求解的 LINGO 程序如下:
    2 C8 K( {( S7 A9 tmodel: 1 K& v5 ?2 |. N6 O" I
    sets:
    $ i  [4 {5 e, cvar/1..5/; 1 J4 @, x) P) I. v9 m; G: C# i7 C! X* {
    link(var,var):c,x; ; E) V4 T) t6 f8 N$ {
    endsets 3 ]: b# R2 [& K& Q/ M5 e( k1 G
    data:
    0 O3 M% `# s% ~: l8 L- H! jc=3 8 2 10 3   
    - n& K# J* a$ I* T  8 7 2 9 7   
    + j6 K" C( C6 o5 c  6 4 2 7 5   
    ) ]5 n8 K8 ^; E, Y% P( ~  8 4 2 3 5   
    4 A; }: D) l# p( w5 t6 F  9 10 6 9 10;
    " }$ j, l; Q  T; ~0 Benddata 4 e+ c, w: N. \1 r
    min=@sum(link:c*x); . W1 ~- ^. j  t& i% g" Y: R4 W+ B  E" Q
    @for(var(i)sum(var(j):x(i,j))=1);
    0 R. Q" Z0 b/ q3 l@for(var(j)sum(var(i):x(i,j))=1); * X0 d2 L' y& h+ Y. \
    @for(linkbin(x)); 6 a* Y0 R) m. }
    end : g! x. s# n' P4 k
    - ]8 a. Y& [6 e" q% s( W& x
    8 V. S4 d9 X+ ^* W' q: e
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    sjlxdn        

    1

    主题

    2

    听众

    155

    积分

    升级  27.5%

  • TA的每日心情
    无聊
    2022-2-19 17:40
  • 签到天数: 30 天

    [LV.5]常住居民I

    国际赛参赛者

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-5-26 05:14 , Processed in 0.381645 second(s), 57 queries .

    回顶部