24小时热门版块排行榜    

查看: 602  |  回复: 0

candybear

铁虫 (初入文坛)


[交流] 【求助】求助 运输问题gams模型纠错

最近用gams编写了一段网络运输问题的程序,虽然程序没有编译错误,可程序运输的结果错误,请求高手指点。模型的代码如下:
c   set of all cells
           /c1 * c14/
       cp  set of all cells
           /c1 * c14/
       cr(c)  set of source cells
           /c1,c5,c9/
       cs(c)  set of sink cells
           /c14/
       ch(c)  set of common cells
           /c1 * c13/
       t   sets of discrete time intervals
           /0 * 10 /
       tfirst(t) intial time(t=0)
           /0/
       tend(t) t=8
           /10/
AI(cp,c) relationship from  cp to c
               /c1.c2
                c2.c3
                c2.c6
                c3.c4
                c4.c14
                c5.c6
                c6.c7
                c6.c10
                c7.c8
                c7.c12
                c7.c13
                c8.c14
                c9.c10
                c10.c11
                c11.c14
                c12.c4
                c13.c11/  ;

       time.. z =e= sum((ch,t),x(ch,t));
       xx(c,t+1)..  x(c,t+1) =e= x(c,t)+sum(cp$AI(cp,c),yi(cp,c,t))-sum(cp$AO(c,cp),yo(c,cp,t))+d(c,t);
       xfirst(c,tfirst).. x(c,tfirst)=e=0;
       yofirst(c,cp,tfirst).. yo(c,cp,tfirst)=e=0;
       yifirst(cp,c,tfirst).. yi(cp,c,tfirst)=e=0 ;
       qi(c,t).. sum(cp$AI(cp,c),yi(cp,c,t))=l=q(c,t);
       qo(c,t).. sum(cp$AO(c,cp),yo(c,cp,t))=l=q(c,t);
       nn(c,t)..  sum(cp$AI(cp,c),yi(cp,c,t))+x(c,t)=l= n(c,t);
       yx(c,t)..  sum(cp$AO(c,cp),yo(c,cp,t))=l= x(c,t);
       xend(ch,tend).. x(ch,tend)=e=0;

display      x.l
AI为网络连接矩阵表达,AO与AI相同,只为方便表达约束条件中的输出流量。

程序的运行结果为:

              1           2           3
c1       27.000      15.000       3.000
c5       15.000       3.000
c9       32.000      20.000       8.000

但恳请 高手帮忙指点 此模型的错误之处  多谢
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 candybear 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 8/400 2026-09-01 12:58 by ssyjh
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 15/750 2026-09-01 12:23 by icm639
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
信息提示
请填处理意见