- 在线时间
- 791 小时
- 最后登录
- 2022-11-28
- 注册时间
- 2017-6-12
- 听众数
- 15
- 收听数
- 0
- 能力
- 120 分
- 体力
- 36467 点
- 威望
- 11 点
- 阅读权限
- 255
- 积分
- 13901
- 相册
- 0
- 日志
- 0
- 记录
- 1
- 帖子
- 616
- 主题
- 542
- 精华
- 12
- 分享
- 0
- 好友
- 225
TA的每日心情 | 开心 2020-11-14 17:15 |
|---|
签到天数: 74 天 [LV.6]常住居民II
 群组: 2019美赛冲刺课程 群组: 站长地区赛培训 群组: 2019考研数学 桃子老师 群组: 2018教师培训(呼伦贝 群组: 2019考研数学 站长系列 |
问题描述: 某市消防中心同时接到了三处火警电话。根据当前的火势,三处火警地点分 别需要 2 辆、2 辆和 3 辆消防车前往灭火。三处火警地点的损失将依赖于消防车到达的及时程度:记 为第 j 辆消防车到达火警地点i的时间,则三处火警地点的损失分别为![]()
q1 ~ n* J/ F$ r- T; Q
# g, F" Y; N( d2 ?9 q
! o4 o) z& m: u2 U( _ ; 。目前可供消防中心调度的消防车正好有 7辆,分别属于三个消防站(可用消防车数量分别为 3 辆、2 辆、2 辆)。消防车从三个消 防站到三个火警地点所需要的时间如表 6 所示。应如何调度消防车,才能使总损失最 小? 2 b8 A' k# m/ @8 |# _
7 }' ~! C) c+ x$ x, A
![]()
. e/ `" w* t8 t% Y4 u6 j) r
1 ~7 P" j, D/ q$ [问题二:如果三处火警地点的损失分别为 ;,调度方案是否需要改变?# Z; O, ~1 D. k
" u* C. ^1 j! S- H
- q% A# s I0 \! g; E, f3 u
# f$ e8 A. ?7 |' l1 m2 q, S3 x6 q
(1)问题分析
) G# G, v3 O$ }3 H1 W
( ?5 K4 w8 ?$ W% M( \, \本题考虑的是为每个火警地点分配消防车的问题,初步看来与线性规划中经典的运输问题有些类似。本题的问题可以看成是指派问题和运输问题的一种变形,我们下面首先把它变成一个运输问题建模求解。6 z5 [2 T7 g" h) H/ A9 Z
9 S% b \( T# J: d+ c( d! R
(2)决策变量& {. u& a0 ~9 U$ w4 o, Z/ U1 z
* X* Z7 H8 s# g8 ~/ J
为了用运输问题建模求解,我们很自然地把 3 个消防站看成供应点。如果直接把 3 个火警地点看成需求点,我们却不能很方便地描述消防车到达的先后次序,因此难以确 定损失的大小。下面我们把 7 辆车的需求分别看成 7 个需求点(分别对应于到达时间 . 用 表示消防站i是否向第 j 个需求点派车(1 表示派车,0表示不派车),则共有 21 个 0−1 变量., u: g5 ~8 N4 w8 |
, }3 N3 Y5 Z! P8 D7 ?# ]: v) }# e(3)模型建立# j9 x' U( k5 p* Z+ g
) F# D" F* ^ y' M0 g- a8 M: d
题目中给出的损失函数都是消防车到达时间的线性函数,所以由所给数据进行简 单的计算可知,如果消防站 1 向第 6 个需求点派车(即消防站 1 向火警地点 3 派车但该 消防车是到达火警地点 3 的第二辆车),则由此引起的损失为 72 98 = × 。同理计算,可以得到损失矩阵如表 7 所示(元素分别记为 )。 2 E$ e; t: V( W q
8 I! i& @- I, P- ?) z$ B4 `! j
3 R ?+ B. f3 T; l- F+ T- v
![]()
! s9 j, N; J6 F3 Q9 a( S
9 U, O9 l# J. ~- u' _' X![]()
& q/ k2 L7 ^+ R( j: Z. J* E) v- v
于是,使总损失最小的决策目标为/ |5 G4 X- N& Q ]# f0 y
! h8 d* H2 M) u: ]+ d Q. M ( 1 )
! X) J o$ v0 n+ s2 E& G. s: S% w7 _3 P5 H
约束条件:
1 P0 d" c4 g" @& v& X V& {5 V
& m2 h' ~) `2 P约束条件有两类,一类是消防站拥有的消防车的数量限制,另一类是 各需求点对消防车的需求量限制。
6 O- C" H" ^/ }0 Y4 ~记 ( i=1,2,3 )为第i个消防站拥有消防车的数量,则消防站拥有的消防车的数量限制可以表示为 ! v+ N4 ~% c6 T# b
( 2 )
$ J- a3 q) R* c3 v; ^ R4 L. E: D; r' o
各需求点对消防车的需求量限制可以表示为 ' b' ^; k, z; \9 c3 ?/ A3 j
4 s! `0 ~; n. q6 }+ U ( 3 )) C! b; ~$ u1 L1 l! ~( S
# T) ?, \" W2 L" z8 o
(4)模型求解 的lingo代码, q) Y$ w7 P2 I$ |7 E
% O+ a2 s# Y- S/ T. @& D$ S/ x
MODEL:
+ X |6 N9 c! L3 H2 r3 R& gTITLE 消防车问题;
% |' d% X8 _, K( eSETS: 4 B* V! x# v: k3 ?: p1 E) _
supply/1..3/:b;
. I( I8 D; M0 Y) Sneed/1..7/; * [6 M5 }: i- T# @& k8 M
links(supply,need):c,x; % m f+ f) R! M. X( f
ENDSETS 9 r' r% y! K) L& x
[OBJ]Min=@sum(links:c*x);
/ M. F( T# H" M+ g@FOR(supply(i):@sum(need(j):x(i,j))=b(i)); 0 H4 Y8 ~2 ^+ s9 o- G
@FOR(need(j):@sum(supply(i):x(i,j))=1);
0 J/ ~; X% j8 c4 E1 hDATA: 7 J2 W; ^7 ?7 Q+ k- J- |0 `5 h
b=3,2,2; , @ `8 V2 R' f. x# q9 o
c=36,24,49,21,81,72,45 3 F/ _ f( l4 n* x& W' g: j
30,20,56,24,99,88,55
" s _7 \5 J# M: t& N; Y 36,24,63,27,90,80,50; 4 {- a- c9 q' w; f+ m
ENDDATA ) `) r* U" _# V6 q: P2 R. o
END
8 Y2 X9 s" ]) Y0 e1 ~9 o, P2 ^求得结果为,消防站 1 应向火警地点 2 派 1 辆车,向火警地点 3 派 2 辆车;消防站 2 应向火警地点 1 派 2 辆车;消防站 3 应向火警地点 2、3 各派 1 辆车。最小总损失 为 329。
) l# {* V, H( L6 f, |7 j
& {+ }8 ?: Y2 ~(5)讨论
0 E& u9 a8 b7 k' u5 `( d3 D W# t9 A1 x
: G# B* D9 h( S' z! Z1)这个问题本质上仍然和经典的运输问题类似,可以把每辆车到达火场看做需求点,消防站看做供应点。在上面模型中,我们虽然假设 为 0− 1变量,但求解时是采用线性规划求解的,也就是说没有加上 为 0− 1 变量或整数变量的限制条件,但求解得到的结果中 正好是 0−1 变量。这一结果不是偶然的,而是运输问题特有的一种性质.
% n7 I4 Y$ j0 ^6 n! y4 O( U: u7 Q, c1 ^' r5 |6 @
2 )在上面模型中,没有考虑消防车到达各火警地点的先后次序约束,但得到的结果正好满足所有的先后次序约束,这一结果不是必然的,而只是巧合。如对例题后半部 分的情形,结果就不是这样了。显然,此时只需要修改损失矩阵如表 8 所示(所示(元素仍然分别记为 )
$ k0 w6 Y% c, z1 k6 Z8 G: q- C4 t) w# @3 i7 k$ `6 A
![]()
^# n, `2 w: w! k4 T p, D4 d4 U- N' F" b
此时重新将式( 1 )-( 3 ) 构成的线性规划模型输入 LINGO 求解,可以得到新的最优解 : 其它变量为 0(最小总损失仍为 329)
: h* h T0 b1 e2 k4 T% p/ }& @: i5 h( M n7 A' y+ z0 W
实际上,损失矩阵中只是 1、2 列交换了位置,3、4 列交换了位置,5、7 列 交换了位置,因此不用重新求解就可以直接看出以上新的最优解。
" k/ z, Y; L0 q8 f7 R* k" L6 u0 F5 [* R3 V9 o0 @+ P( h% V
但是,以上新的最优解却是不符合实际情况的。例如, 表明火警地点2 的第一辆消防车来自消防站 3,第二辆消防车来自消防站 1,但这是不合理的,因为 火警地点 2 与消费站 3 有 9min 的距离,大于与消防站 1 的 7min 的距离。分配给火警 地点 3 的消防车也有类似的不合理问题。为了解决这一问题,我们必须考虑消防车到达 各火警地点的先后次序约束,也就是说必须在简单的运输问题模型中增加一些新的约 束,以保证以上的不合理问题不再出现。& K: Y5 R7 o, z% h2 v6 Y& C
; W5 g9 h. y/ M. J. l: q
首先考虑火警地点 2。由于消防站 1 的消防车到达所需时间(7min)小于消防站 2 的消防车到达所需时间(8 分钟),并都小于消防站 3 的消防车到达所需时间(9 分钟), 因此火警地点 2 的第二辆消防车如果来自消防站 1,则火警地点 2 的第 1 辆消防车也一 定来自消防站 1;火警地点 2 的第 2 辆消防车如果来自消防站 2,则火警地点 2 的第 1 辆消防车一定来自消防站 1 或 2。因此,必须增加以下约束: ( 4 )
# O; i( {: ~+ D+ Y+ k& D h
# q5 _1 a8 s b7 b1 N同理,对火警地点 1,必须增加以下约束: ( 5 )
5 F& h0 M V! K+ W; D1 w" }+ v9 }$ z8 V+ m) X$ G# y% a. u7 y# i$ j2 f
对火警地点 3,必须增加以下约束: ( 6 )$ j# g& j0 t- I, D# L. c) l9 [) T- p
; {$ H8 ^( [: ~: t3 t4 ^, z/ h! v重新将式(1)~(6)构成的整数规划模型( 是 1 0− 变量)输入 LINGO 软件如下:3 s/ r T- u( R3 I; J# v" B
$ Z2 t) M$ t. r% U3 u; X7 {( B/ gMODEL:
5 k- ^- `, x4 \4 _' VTITLE 消防车问题; ; v& t5 C+ A! S: p- l
SETS: / u7 d7 L+ l% o- r X8 s
supply/1..3/:b; ( V2 z" `. s7 |5 ^- D J
need/1..7/; 5 L/ X! F& j3 r
links(supply,need):c,x;
}- l9 A6 C" MENDSETS , w, C7 q: `7 A e7 g
[OBJ]Min=@sum(links:c*x);
/ o9 ]/ [' L0 v9 ~: F9 V@FOR(supply(i):@sum(need(j):x(i,j))=b(i)); ( X) u* b2 }* Z( d2 Z& {1 {2 w+ H
@FOR(need(j):@sum(supply(i):x(i,j))=1);
% ?: m! ~6 y( c/ k/ [! n6 S( Rx(1,4)<x(1,3); 4 F( y$ Y. L9 v8 L
x(2,4)<x(1,3)+x(2,3); . s |* z% _+ Y& ]) h- @$ v$ ?
x(2,2)<x(2,1);
% [. }+ F# S' i, G D: ~5 ~ l" vx(1,6)<x(1,5); + C2 D% `2 {; b
x(1,7)<x(1,6); 8 Z L" u7 O* B" q& M. O- e f
x(3,6)<x(1,5)+x(3,5);
4 u+ {7 Y7 W' O2*x(3,7)<x(1,5)+x(1,6)+x(3,5)+x(3,6);
! U2 r& C0 G+ Z: G9 _2 C@for(links:@bin(x)); 6 l% {4 E( \: X( U2 Z8 z
DATA:& @: Y+ Y. {* x y7 H8 z
b=3,2,2;
0 L& a" W8 Y6 b5 f; Zc= 24 36 21 49 45 72 81 : d/ a& ?8 t1 o/ H
20 30 24 56 55 88 99
3 Y- z5 A4 V% S/ I& O; o1 p 24 36 27 63 50 80 90; 8 N: i! F7 U: w& }3 R
ENDDATA ' O) A& f- ]9 A7 I# E. W. D
END & G2 m5 Q9 N0 \2 J# b( l5 x+ I
求解可以得到: ,其它变量为 0(最小总损失仍为 335)。也就是说,消防站 1 应向火警地点 2 派 2 辆车,向火警地点 3 派 1 辆车;消防站 2 应向火警地点 1 派 2 辆车;消防站 3 应向火警地点 3 派 2 辆车。经过检 验可以发现,此时的派车方案是合理的。 z. l- I$ b- x* s, l- }
————————————————2 Z/ M3 c0 Y2 @
版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
- l4 m) X1 h2 \+ n* ^/ K% n原文链接:https://blog.csdn.net/qq_29831163/java/article/details/893881726 G' i! ?( x- f' X8 X
; t' L/ I" s8 a) ?: V1 ]- s4 n
! Z4 x7 u' o; b6 C( m$ C# i/ c2 B! C# |4 I
|
zan
|