数学建模社区-数学中国
标题:
线性规划（二）：运输问题 (产销平衡) & 指派问题
[打印本页]
作者:
浅夏110
时间:
2020-6-6 10:08
标题:
线性规划（二）：运输问题 (产销平衡) & 指派问题
1.可以转化为线性规划的问题
5 R3 ~1 w [; I& x- R+ Z" ^: T6 A
很多看起来不是线性规划的问题也可以通过变换变成线性规划的问题来解决。下面举几个例子：
& G! y3 u$ }- `- `4 ~
+ \; J7 T/ G0 C$ j8 Q* r
; H1 q' T$ v& H. C8 t1 \' o1 P" I, g
7 c$ \5 z! _3 f! f1 R" e4 i
8 K& l, s) D4 Z+ ~4 d) m3 D% l6 m. H
- h* [, Q2 i2 P/ l3 K
' Q8 j5 B( s A; B* T
2. 运输问题(产销平衡) ：康—希表上作业法
" Q" d$ ~) u0 U% A' X) d$ o
! M6 g( d6 U1 i9 r) b. X
) ]3 L3 t/ E: l( H. w
8 p5 I( [3 B' V" Z/ ~4 x _
3. 指派问题的数学模型
9 _( E. E% g) G S1 K8 P( ?
& v! p4 T- @9 t, n
0 z3 x) C# h3 a. x: P# E1 L3 Q
8 R- o' ?/ [, C X/ {
) h' d) [9 S/ p; e" T* ~" L& ~/ d
上述指派问题的可行解可以用一个矩阵表示，其每行每列均有且只有一个元素为 1，其余元素均为 0；可以用 1,...,n 中的一个置换表示。 问题中的变量只能取 0 或 1，从而是一个 0-1 规划问题。一般的 0-1 规划问题求解 极为困难。但指派问题并不难解，其约束方程组的系数矩阵十分特殊（被称为全单位模矩阵，其各阶非零子式均为 ），其非负可行解的分量只能取0或1，故约束 =0或1 可改写为 而不改变其解。此时指派问题被转化为一个特殊的运输问题，其中m=n, .
8 @. }) `: e! K7 H l* r1 E$ Z/ A! D+ R2 j5 T
3 p" {# K3 E9 E9 p
求解指派问题的匈牙利算法
: I$ |3 }& Q# X8 Q/ {/ ~, v4 ^
9 D% ~; W9 ~, f; A+ l; O1 V) L
0 j: G! E# b1 z8 |! E
0 ?8 _; t- h E8 B1 _, s2 k& F, q
; K ?& a2 ~9 h9 w3 ~# k$ L, s+ U
有时问题会稍复杂一些。
4 ]) O" v# m# J0 n) R2 a; H+ e
% J$ |; [5 k/ @* N, s
例 9 求解系数矩阵C 的指派问题
- V, L1 }8 t# l4 q
' Y) V+ {3 s+ P& G! e) M0 w- F5 K& H# m
5 S V7 h' ~9 V8 H* _5 y
' t# q% y6 Y0 R$ i
7 J2 }/ O" P, Z: P
: X1 _ J! X8 @& {- o
/ r7 K% r2 H+ x+ j; S1 e
/ x% }- O! H: k, j3 E
& M' ~6 f) d6 l8 B' V. |* n" @
* D/ u5 U' j3 z; O3 R
4. 指派问题的计算机求解
5 P( W4 z) |# s, h1 K% t' M
整数规划问题的求解可以使用 Lingo 等专用软件。对于一般的整数规划问题，无法 直接利用 Matlab 的函数，必须利用 Matlab 编程实现分枝定界解法和割平面解法。但对 于指派问题等 0−1整数规划问题，可以直接利用 Matlab 的函数 bintprog 进行求解。
1 W' j! J* O- Z: X
% _2 M8 Y! [+ G0 N
& w( x1 y: C) W7 H7 x
& w9 ?1 f9 O" T5 c3 x( m: t
解：编写 Matlab 程序如下： 求得最优值为 21，最优指派方案为
. u9 N" \% S! M2 \! g1 O+ N
% G% ?4 r; U. S5 L6 n2 V( B q& H
c=[3 8 2 10 3;8 7 2 9 7;6 4 2 7 5
: f ~$ t5 v0 r7 O- J8 Y
8 4 2 3 5;9 10 6 9 10];
% Q; v' l& \8 A7 w& K5 r& I8 e
c=c(
;
; S1 H8 P- I. R
a=zeros(10,25);
: v3 v4 H" {9 V1 h7 x4 m8 d% n
for i=1:5
1 L" X' ?" ^& ~, S& i
a(i,(i-1)*5+1:5*i)=1;
?& Q' e; H1 O
a(5+i,i:5:25)=1;
! `3 E( g. I/ e/ v2 x
end
. Q9 @; G, U$ P
b=ones(10,1);
4 }$ i8 N; _ }5 x6 t
[x,y]=bintprog(c,[],[],a,b);
& |9 o' H2 i; }2 b/ s, {0 c( J
x=reshape(x,[5,5]),y
+ Q1 k' a* i4 q, T8 j( u
# \7 o2 X- @8 u$ |+ d2 |8 r
求解的 LINGO 程序如下：
6 d; ~ ~! `. d$ Z6 E. s0 l* A$ K
+ v, F8 h$ l) e s
) F' b" v; J, J% ]% E2 I
model:
# y4 |( m+ J' W4 S
sets:
7 Q, y6 _4 {9 b" v+ s! R9 h
var/1..5/;
7 p7 o9 b0 y5 v( N: p$ K6 r
link(var,var):c,x;
9 B! C/ i: j' H; B! D8 Z
endsets
5 a; s% g' u2 L
data:
& n1 l5 C$ U' Q' X1 M6 W
c=3 8 2 10 3
8 {& ?2 j: T, q, i O2 e
8 7 2 9 7
1 e' u+ `0 E2 o4 t$ N
6 4 2 7 5
! C/ D* c% X4 F" K, A2 S: p7 R
8 4 2 3 5
4 q: Q/ v1 b2 V7 c& H' C& W
9 10 6 9 10;
e* N: w' q' }$ p
enddata
& H; b* @/ {) [5 h, }! [
min=@sum(link:c*x);
, D" b( M! ~5 e+ O6 v$ O0 w, B
@for(var(i)
sum(var(j):x(i,j))=1);
) }# b& \0 M' J& N+ n+ P1 E
@for(var(j)
sum(var(i):x(i,j))=1);
; L: Y* m5 h" d- y2 C: k
@for(link
bin(x));
) B7 T2 |& ^) s# ]
end
- d2 G9 o0 z- ?5 O2 l) }
6 [- b, O/ J" ^1 b+ Z0 p$ F
————————————————
4 l5 `* w# H( z/ D
版权声明：本文为CSDN博主「wamg潇潇」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
2 Y& s; m3 [: p1 x# i* x) S9 C
原文链接：https://blog.csdn.net/qq_29831163/article/details/88894966
; [) \' e0 I1 U; e% c
f# e: P0 e+ U8 Q
3 Z7 w7 }4 m9 U- q
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
8 K& l, s) D4 Z+ ~4 d) m3 D% l6 m. H
' t# q% y6 Y0 R$ i
; ; S1 H8 P- I. R
sum(var(j):x(i,j))=1); ) }# b& \0 M' J& N+ n+ P1 Esum(var(i):x(i,j))=1); ; L: Y* m5 h" d- y2 C: kbin(x));