24小时热门版块排行榜    

查看: 1299  |  回复: 1

qinjuan1983

新虫 (初入文坛)

[求助] lingo 求解二元优化问题,怎么总是没有可行解呢,应该是有的~多谢 已有1人参与

程序如下:
N,T为变量;
max=p*(@exp(-k*N))*dPN-A/T-V*dPN*(@exp(ct*T)-1)/T/ct-h*dPN*((@exp(ct*T)-1)/ct-T)/T/ct+p*IE*dPN*(M-N-T/2);
dPN=D0-a*p+b*c*@exp(c*N);
T+N<=M;
ct=0.1;a=2;b=2;c=0.1;p=20;k=0.2;A=20;v=10;D0=1000;Ie=0.01;M=1;h=2;
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NicoleLas

木虫 (正式写手)

【答案】应助回帖

修改后的代码:
CODE:
max=p*(@exp(-k*N))*dPN
   -CapA/T-V*dPN*(@exp(ct*T)-1)/T/ct
     -h*dPN*((@exp(ct*T)-1)/ct-T)/T/ct
      +p*IE*dPN*(M-N-T/2);

dPN=D0-CapA*p+b*c*@exp(c*N);

T+N<=M;

ct=0.1;a=2;b=2;c=0.1;p=20;k=0.2;CapA=20;v=10;D0=1000;Ie=0.01;M=1;h=2;

求解结果:
  Global optimal solution found.
  Objective value:                         5844.24124892
  Objective bound:                       5844.24131379
  Infeasibilities:                            0.00000000000
  Extended solver steps:                          2363
  Total solver iterations:                        232119
  Elapsed runtime seconds:                    12.00
  Model Class:                                       NLP
  Total variables:                         3
  Nonlinear variables:                  3
  Integer variables:                      0
  Total constraints:                      3
  Nonlinear constraints:               2
  Total nonzeros:                        7
  Nonlinear nonzeros:                 4

                                     Variable                Value             Reduced Cost
                                            P        20.0000000000            0.00000000000
                                            K       0.200000000000            0.00000000000
                                            N        0.00000000000            2520.64061804
                                          DPN        600.200000000            0.00000000000
                                         CAPA        20.0000000000            0.00000000000
                                            T       0.143666483948            0.00000000000
                                            V        10.0000000000            0.00000000000
                                           CT       0.100000000000            0.00000000000
                                            H        2.00000000000            0.00000000000
                                           IE       0.100000000000E-01        0.00000000000
                                            M        1.00000000000            0.00000000000
                                           D0        1000.00000000            0.00000000000
                                            B        2.00000000000            0.00000000000
                                            C       0.100000000000            0.00000000000
                                            A        2.00000000000            0.00000000000

                                          Row         Slack or Surplus           Dual Price
                                            1        5844.24124892            1.00000000000
                                            2        0.00000000000            9.96909790547
                                            3       0.856333516052            0.00000000000
                                            4        0.00000000000           -439.459146366
                                            5        0.00000000000            0.00000000000
                                            6        0.00000000000           0.996909786833
                                            7        0.00000000000            19.9381958109
                                            8        0.00000000000            406.388901255
                                            9        0.00000000000           0.248104333878E-05
                                           10        0.00000000000           -206.342521329
                                           11        0.00000000000           -604.532150019
                                           12        0.00000000000            9.96909790547
                                           13        0.00000000000            11141.7137658
                                           14        0.00000000000            120.040002480
                                           15        0.00000000000           -43.3215224445
Be water, my friend.
2楼2014-10-17 18:57:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qinjuan1983 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 07:09 by 4wMiSEwB6436
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +24 医学老男孩 2026-08-13 50/2500 2026-08-15 06:38 by 学员8dgXkO
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +11 Tide man 2026-08-10 12/600 2026-08-15 02:12 by home3163
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 7/350 2026-08-15 00:40 by 4wMiSEwB6436
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见