QQ登录

只需要一步,快速开始

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

GAMS示例——运输问题 【代码及演算过程】

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

1253

主题

443

听众

-516

积分

复兴中华数学头子

  • TA的每日心情
    开心
    2011-9-26 17:31
  • 签到天数: 3 天

    [LV.2]偶尔看看I

    自我介绍
    数学中国网站(www.madio.cn)是目前中国最大的数学建模交流社区

    邮箱绑定达人 优秀斑竹奖 发帖功臣 元老勋章 新人进步奖 原创写作奖 最具活力勋章 风雨历程奖

    群组越狱吧

    群组湖南工业大学数学建模同盟会

    群组四川农业大学数学建模协会

    群组重庆交通大学数学建模协会

    群组中国矿业大学数学建模协会

    跳转到指定楼层
    #
    发表于 2009-11-14 04:27 |只看该作者 |正序浏览
    |招呼Ta 关注Ta |邮箱已经成功绑定

    GAMS示例

    9 L9 a8 c/ Q0 d9 f& T. G: P$ E7 {- n

    下面这个问题主要是用来举例说明GAMS是如何让您以一种自然的方式建立模型.GAMS能够处理大的多和高度复杂的问题.在这里只展示了GAMS的很少一部分的基本特征.

    9 _! R4 D4 I+ B( w ) R3 }( ^$ w% P, _ {$ {7 Q. ?+ v+ I# P8 h" b. X+ L$ l: ]7 b/ C4 e# s+ T b* R- r7 R* q0 j! ~3 O6 @0 D1 C; t
       代数描述
    % I: m5 W$ ` `( ]) H( E" M

    下面是问题的标准的代数描述,这是用来最小化把货物从2个工厂运输到3个市场的费用.受约束于供需约束.

    Z4 k0 E0 C& [: w) u

    指标:

    ( P& _/ `9 M# y" o( [9 ?7 l

     i=工厂(plants)
    j=市场(markets)

    * W6 o7 J2 g$ H: ^+ R

    给定的数据:

    ' S7 k6 ^ m# D

    =在工厂i日常供应量(批数cases)
    =在市场j的日常需求量(批数cases)
    =在工厂i和市场j之间的距离(千英里)
    =在工厂i和市场j之间每单位的运输费用($/批/千英里)

    $ c: D" @ z* t/ h( h# i7 ? ' F& d) r0 L) z; ?- ]5 p9 V5 h* L- f) E$ B& {3 O$ M+ r+ a$ F9 a5 c G# s9 I: s2 |" C" _; \8 H j( H: Y& \8 p3 ~' ~/ S3 B8 @1 R0 @+ A; u1 z9 t. Y$ S F5 K$ H+ L/ \ ^* L) D9 \! b& l# |% [ S( y8 v. X) F; O7 W4 g' c; M- w' p' }1 F" K) I6 J- O, T! D" @6 T$ M' e7 d2 g P0 V3 o4 d3 s; M. b* {3 `# G) b; H! W, O, G& ^3 h* J# p! R0 y. a' H4 n% K# }' q! M0 ]4 G( q: Q# t' ^/ w- h. O: E, R, t/ Y# s5 _9 t4 C0 I0 u8 x) D7 p9 e8 O& r( _$ u/ W2 t8 z: g; a; `. }% y/ } W( c: z( o8 L4 ~' h# P$ ?$ F* I" A) f5 d- v8 t# J! S S: I% F$ l9 @! o; c) Q) u0 ^$ b( h# P3 t# g0 p1 W- K% k+ F- R9 `1 ^) u( n% V5 R% v0 r/ P/ m8 f; _+ P; q @5 d! f; ] H0 ^, T# d$ {) z% N6 e1 n! D; {+ U5 P4 c7 [+ B* B- l0 X3 k: v, H6 n4 H8 p$ U5 Q( B% X( Y( K; c( L$ e
    距离  
      市场  
    工厂 New York Chicago Topeka 供应量
    Seattle 2.5 1.7 1.8 350
    San Diego 2.5 1.8 1.4 600
    需求量 325 300 275  
    - H- b/ u/ d" ^- t, s7 `" z. _

     F=每批每千英里的费用$

    $ `8 A8 }- x- s& ^3 N- U' ~

    决策变量:

    1 f" e9 I/ G0 g: G

    =日常从工厂i运输到市场j的总量(批数cases)
    这里 适用所有i,j

    4 \/ X. f6 `( G5 ^7 e" W; U

    约束:

    ' k% G+ b- q% O+ r6 m u1 H# U; N" R7 E

    在工厂i的供应量**(批数cases):适用所有的i
    在市场j的需求量:适用所有的j

    : a' \0 @+ f4 @% L6 z; J

    目标函数:

    / F) Q2 A1 v6 f. ]6 t E

    Minimize (千元)

    8 Z. p$ D0 t5 @; @! K& @0 Z( N1 u6 ` $ h$ _, A. K) z$ n! v7 P& A7 \3 p, g( h- T1 S' D3 d& a, L& p1 v* r! A# v/ O+ d" l% B- f
       GAMS模型
    / b: T. C( T1 |

    同样的模型在GAMS中建模.简练的代数描述使得模型高度紧凑,并带有逻辑结构.内部的文档,比如对参数的解释和测量的单位,使得模型很容易读懂.
    ; U8 O/ d; {, j3 h

    集合(Sets)

    }. F$ @6 P4 J/ ]0 M K

    9 N0 [1 @6 C# J( H' x

    GAMS让您以直接的方式指定指标:声明和命名集合(这里是I和J),并列举它们的元素.

    ) C2 X; }& o0 z# E, t7 U

    参数

    - D& W3 Y( D) V

    2 `' ~) u0 G6 _! q$ {

    这里的数据输入被作为指标参数A(I)和B(J),值简单的被列出.

    ; h! n T, f6 C

    GAMS让您可以在模型的任意位置放置解释性文本(以小写格式显示),当您在开发它时.您的注释自动被结合到输出报告中的合适位置.

    % L5 @8 ^/ o4 f9 ]) V

    表格

    ( M: J* u( x1 g3 P4 ~0 o


    数据同样能够以方便的表格形式输入.GAMS让您以数据的基本形式来输入数据-转换是特定的代数化的.

    ' F* ^: p! Z: ]- e3 j

    标量(Scalar)

    $ r2 q) h8 @% X3 P6 p


    常量能够被声明为标量,它的值是指定的. ' [1 ^! p5 h! ]/ z( R7 z" Y0 P* a: v

    数据处理

    ! O2 }3 W% ~1 B& \' k8 V9 s$ m {

    [5 B' {: R& p A

    当数据值要被计算前,您首先要声明参数(比如,给它一个符号,随意给它编个指标),然后给它一个代数公式.GAMS将自动进行计算.

    # Q! o5 g0 J0 J* f

    变量

    3 i0 O3 p+ Q1 p% @1 |7 Q1 W

    * z z* t9 x* h& U# O9 S& u5 j

    决策变量以代数的方式表达,带有特定的指标.从这种常见的形式,GAMS在域中生成变量的每个实例.

    o) o3 v8 [) T3 ~2 c

    变量可以被指定为下列类型:任意(FREE),正值(POSITIVE),负值(NEGATIVE),二元(BINARY),或者整数(INTEGER).默认是任意(FREE).

    / t# C* `" L/ o6 W! \

    目标变量(这里是z)仅被声明,没有指标(index).

    * R( ?; O# F4 l6 m

    方程式

    ' p% F% q- a5 w5 S/ L

    ' C, z% F! s6 j& E. t' A9 u

    目标函数和约束方程式首先被通过指定名字来声明.然后它们的概括的算术公式被声明.GAMS现在已经有了足够的信息(从上面的数据输入和从在方程式中指定的算术关系)来自动生成每个单独的约束声明-就像您能在下面的输出报告中看到的.

    0 B4 \ X1 J* ?. k* p8 ^1 n% q

    =E=表示'equal to' (等于)
    =L=表示 'less than or equal to' (小于或等于)
    =G=表示'greater than or equal to' (大于或等于)

    + o3 _- ]) |; ~7 U! ?- q

    模型声明

    2 v2 S& M3 E% ]

    ' m$ y) l: l, _" W/ R

    模型被指定了一个唯一的名字(这里是TRANSPORT),模型缔造者指定那个方程式应该被包含到这个特别的公式中.在这里我们指定了ALL,也就是说所有的方程式都是模型的一部分.就等于是MODEL TRANSPORT /COST, SUPPLY, DEMAND/ . 这个方程式选择使您能够在单个的GAMS输入文件中以公式表达不同的模型,基于相同或不同的给定数据.

    % ]/ N4 l3 E6 A& p

    求解声明

    * M; W t, m' I. _# r) L

    0 P# i. `0 L' I/ a, l

    求解声明(1)告诉GAMS那个模型要被求解,(2)选择要使用的求解器(在这里是LP求解器),(3)表明优化方向,或者是求最小值,或者是求最大值.(4)指定目标变量.

    4 w3 Y- m( H6 K5 _ # H8 H2 n9 ]& J% g( _# d7 @6 p4 C& i; U, }. @/ _; [" b. H. E3 u. ~& Q/ ^4 d7 |4 J" f* k K( W. s. W+ `3 \9 h, [2 }* R
       GAMS输出报告(部分摘录)
    3 x, a; Y9 _! R' M

    完整的GAMS输出报告比下面列出的部分摘录详细的多,包含了更多的帮助用于解释和诊断您的模型.甚至您能够修改输出格式来符合您的特定的需要.

    % N0 G* k* c# P' F+ r4 x0 ^) i- C

    方程式列表

    8 X: i/ O. W1 J1 N4 ^' ?

    $ C6 M( O0 v0 |, H( w

    方程式列表显示从在GAMS输入中指定的分区(block)生成的单独的约束.在GAMS中使用者可以以一种非常紧凑的形式写下被索引的方程式分区(block),这将产生大量的单个方程式.在我们的示例中,我们指定了3个方程式分区,生成了6个单独的方程式.

    # g; A( N$ h) R4 D! ]' w2 L

    列列表

    F8 c0 ~7 g) I0 T! o8 B

    ( R8 ?# r1 N8 L) A: T' G

    列列表提供信息到生成的单独的变量上.变量X(I,J)扩展出6个单独的变量.当许多变量从一个分区中被生成,默认的列表只显示最初的3个(用户可以修改).

    6 G3 T% ]6 ^' }2 V; V y: J# u7 C/ I0 e& }9 i( ~0 h- Q$ r2 k3 m1 t1 W- G1 D1 C5 m7 o# m* ?2 ^1 m- ^; t p3 F7 ]" d+ P' d, D* N0 t
       求解信息
    0 D5 e. _( ^' J# Y

    ' C+ r `8 U7 H; K6 i9 O; `% _

    求解声明将生成模型(单个方程式和对应到特定模型的变量的产物).首先一些关于生成的模型的统计表将会被显示:方程式数,变量和非零元素.

    * e% ?' K6 x( j

    在求解汇总信息部分,我们看到BDMLP被调用来求解这个模型.BDMLP经过4次反复,耗时0.18秒找到了这个问题的最优解.求解信息下列的消息来自求解器.

    , P7 A. l1 m5 O/ H

    解(Solution)

    $ f4 H' a% V2 z6 E% E9 I

    3 R) i4 d. o0 g& Q& n+ S7 p

    解被显示在这里.边际值(marginals)对应方程式的重复和变量减少的花费.

    ( L. p p# n$ w: h+ \

    写工具不需要学习一门其它的语言.在GAMS所有的数据处理,模型定义和报告编写都是在一个单独的环境中完成的.

    : k- Z/ e" p5 s1 r/ d$ y5 F ; t g$ J2 ?, M) F6 F1 o) v2 j- x3 T4 N4 c2 ?' x% b7 I: b' H( w" ]; e7 j* D" a1 p, h$ x3 L( I( x
       参考
    3 ^8 w' x8 k1 x3 h1 p1 ^

    Dantzig G. B., Linear Programming and Extensions, Princeton University Press,Princeton, New Jersey, 1963, Chapter 3-3.

    GAMS Example.pdf

    143.39 KB, 下载次数: 57, 下载积分: 体力 -2 点

    GAMS 示例 (PDF格式)下载

    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏1 支持支持2 反对反对0 微信微信
    板砖sy        

    0

    主题

    9

    听众

    5

    积分

    升级  0%

    该用户从未签到

    自我介绍
    系统工程专业
    入门看这个就足够了,但是一旦独立进行编程的话,会出现很多bug,还有错误代码,都不知道什么意思,大神有没有解决这些问题的方法,辛苦了
    6 B) h# f6 W! O% h2 G
    回复

    使用道具 举报

    3

    主题

    14

    听众

    27

    积分

    升级  23.16%

  • TA的每日心情
    开心
    2017-3-1 21:26
  • 签到天数: 1 天

    [LV.1]初来乍到

    回复

    使用道具 举报

    LYJA        

    20

    主题

    3

    听众

    375

    积分

    升级  25%

  • TA的每日心情
    开心
    2017-8-27 23:29
  • 签到天数: 61 天

    [LV.6]常住居民II

    社区QQ达人

    回复

    使用道具 举报

    鸣凤123        

    0

    主题

    12

    听众

    244

    积分

    升级  72%

  • TA的每日心情
    开心
    2016-1-29 08:03
  • 签到天数: 34 天

    [LV.5]常住居民I

    群组2014美赛MCMA题备战群

    回复

    使用道具 举报

    7

    主题

    9

    听众

    779

    积分

    升级  44.75%

  • TA的每日心情
    开心
    2016-5-19 15:19
  • 签到天数: 200 天

    [LV.7]常住居民III

    自我介绍
    我是一个勇于挑战自我的人。

    社区QQ达人 新人进步奖

    群组西安交大数学建模

    群组学术交流B

    群组第一期sas基础实训课堂

    群组认证人才交流群

    回复

    使用道具 举报

    0

    主题

    6

    听众

    30

    积分

    升级  26.32%

  • TA的每日心情
    开心
    2013-4-10 08:59
  • 签到天数: 13 天

    [LV.3]偶尔看看II

    群组东北三省联盟

    wsq200808 发表于 2012-12-9 22:55
    3 W0 W- n$ p! Q6 j. k! s9 I你的时间怎么是明天呢?

    0 ?$ b4 p7 Y1 l0 j我是外星人。。。嘎嘎嘎
    回复

    使用道具 举报

    wsq200808        

    0

    主题

    7

    听众

    10

    积分

    升级  5.26%

  • TA的每日心情
    奋斗
    2012-12-9 23:18
  • 签到天数: 1 天

    [LV.1]初来乍到

    自我介绍
    数学掌控世界,统计预测未来!
    邓浩考研 发表于 2012-10-10 09:47
    " j; i8 I1 g- w  g帖子很好啊,入门很好啊。
    8 r( M. V( I0 F4 z
    你的时间怎么是明天呢?

    点评

    邓浩考研  我是外星人。。。嘎嘎嘎  详情 回复 发表于 2012-12-10 10:08
    回复

    使用道具 举报

    0

    主题

    6

    听众

    30

    积分

    升级  26.32%

  • TA的每日心情
    开心
    2013-4-10 08:59
  • 签到天数: 13 天

    [LV.3]偶尔看看II

    群组东北三省联盟

    帖子很好啊,入门很好啊。

    点评

    wsq200808  你的时间怎么是明天呢?  详情 回复 发表于 2012-12-9 22:55
    回复

    使用道具 举报

    市_长 实名认证       

    0

    主题

    0

    听众

    5

    积分

    升级  0%

    该用户从未签到

    回复

    使用道具 举报

    0

    主题

    4

    听众

    30

    积分

    升级  26.32%

  • TA的每日心情
    难过
    2012-4-24 16:16
  • 签到天数: 3 天

    [LV.2]偶尔看看I

    自我介绍
    数学建模
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-7-22 01:11 , Processed in 0.498410 second(s), 111 queries .

    回顶部