QQ登录

只需要一步,快速开始

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

[问题求助] 请求一个整数规划的matlab程序

[复制链接]
字体大小: 正常 放大
huiwen     中国数模人才认证   

4

主题

3

听众

258

积分

升级  79%

  • TA的每日心情
    开心
    2014-12-3 16:59
  • 签到天数: 84 天

    [LV.6]常住居民II

    群组Matlab讨论组

    群组西安交大数学建模

    群组2011年第一期数学建模

    跳转到指定楼层
    1#
    发表于 2011-8-25 16:51 |只看该作者 |正序浏览
    |招呼Ta 关注Ta
    最近一直在研究整数规划,可是连着看了两个相关编程,用MATLAB实践时总是出错!所以,求助于各位大侠们能否提供一个整数规划的程序。或者,能否分析一下面的程序哪里有问题?谢谢了啊!
    1. function [x,y]=IntLp(f,G,h,Geq,heq,lb,ub,x,id,options)! I\" a( N$ z( x- D- E
    2. global upper opt c x0 A b Aeq beq ID options;7 J/ C/ B* D; x8 O- o6 K4 a# k+ x
    3. if nargin<10,* T% {/ S- l) Q$ o* L
    4.     options =optimset({});
      3 l+ G: R' e, J' }& C
    5.     options.Display='off';8 N$ @  v( O. o
    6.     options.LargeScale='off';
      $ I) O  [/ \7 }; w
    7. end) [  n6 s9 N6 t* _  e! x+ d+ F
    8. if nargin<9,
      % m2 k; Y7 z! q
    9.     id=ones(size(f));, [3 L: s  }5 f+ C* b# w) {: [
    10. end
      / g\" |# u1 l, d: ]1 c6 Q, i8 V
    11. if nargin<8,$ E' _/ ]! T3 Q
    12.     x=[];\" @) `1 U8 ~& X9 i
    13. end$ f3 D* F3 U7 ?$ ?, T
    14. if nargin<7 | isempty(ub),
      6 `1 B2 }, M( _
    15.     ub =inf*ones(size(f));6 D4 N0 I5 B7 g  l$ l\" {4 b4 j5 Y4 B
    16. end
      ( h2 a$ O& Z+ i
    17. if nargin<6 | isempty(lb),) f% p3 ~  a. @/ ^/ I- s. T4 z
    18.     lb=zeros(size(f));4 x& W( ]& C5 d' P4 P1 {9 F. g- F
    19. end/ B4 j4 j+ G8 a
    20. if nargin<5,( B* D3 {0 `; B
    21.     heq=[];
      . ?5 \4 V1 S8 {9 ^$ h
    22. end6 \0 r; F8 m  p\" \- w) P3 O/ S8 g  H
    23. if nargin<4,
      , D8 T5 ?0 P8 e5 E- G
    24.     Geq=[];
      - |$ h5 a( ^; u% `. X9 _* j
    25. end  W* o( r+ g9 {& o
    26. upper=inf;c=f;x0=x;: ]) Q\" D' B0 m\" Y
    27. A=G;5 z0 `- }6 T$ e( h. t% Q! e. K
    28. b=h;Aeq=Geq;
      2 }( }: X9 B$ T0 p
    29. beq=heq;ID=id;
      ' ^/ ?* k8 B0 z
    30. ftemp=IntL_P(lb(:),ub(:));
      7 z& {; r1 n+ w; }) f
    31. %以下为子函数
      2 {- \# G  }8 B  \
    32. function ftemp=IntL_P(vlb,vub)
      $ }/ x+ \9 _5 \: ]0 T
    33. global upper opt c x0 A b Aeq beq ID options;/ z1 g# G# ]8 t! |
    34. [x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);. ]6 f' {6 w* A3 |) r1 |\" [
    35. if how<0
      4 O5 g; L' e) r8 N% U
    36.     return;
      4 h: f# o# M* V' Y7 R. [+ d: Q
    37. end
      : Z; i8 r% g+ U# @2 ^0 {. s
    38. if ftemp-upper>0.00005
      , }) h' J8 \/ ~7 [
    39.     return;
      4 [6 g1 a5 ?+ D$ D2 H
    40. end
      5 \$ Y3 M& R/ ^- B  r\" j! q% b. @
    41. if max(abs(x.*ID-round(x.*ID)))<0.00005( t3 Q+ u: R7 X1 \! o
    42.     if upper-ftemp>0.00005
        h4 _0 a- t; @\" W' _! O8 k4 \
    43.         opt=x';. A3 M+ s$ q( O8 ~% E
    44.         upper=ftemp;
      8 v8 g! u& i; B2 L* f, K
    45.         return;1 G9 n, Y4 X6 `% k! g$ V
    46.     else
      * I% p, `- J% f  b
    47.         opt=[opt;x'];
      ' M' D( ^$ T, V: d3 D5 v; L
    48.         return;
      5 b- W9 ?- Q1 [0 ~8 h% l. R% Y; v
    49.     end
      7 f! J. J; m. O
    50. end
      / m0 a0 [# `\" H
    51. notintx=find(abs(x-round(x))<=0.00005);8 v; K. D+ E4 ~) a7 ~8 N) q* z6 [- |
    52. intx=fix(x);9 g: j% O- i& B4 H
    53. tempvlb=vlb;tempvub=vub;* N9 \; v+ |' ]3 q) o& ]4 \7 l
    54. if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
      4 _7 o7 x/ R9 ]% b) H\" w4 v
    55.     tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;8 u9 j; V5 E# R2 i2 F
    56.     ftemp=IntLP(tempvlb,vub);3 l+ M& K\" S) \' I7 d3 [
    57. end4 e6 p- z) i: Q2 ~: w
    58. if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
      7 G8 |! D; Y# ]- L9 ?; I/ Y! M
    59.     tempvub(notintx(1,1),1)=intx(notintx(1,1),1);. D* B$ E. W/ t$ Q3 p
    60.     ftemp=IntL_P(vlb,tempvub);
      + w3 n' b6 g0 |8 I: ^\" ]# ~
    61. end
    复制代码
    zan
    转播转播0 分享淘帖0 分享分享1 收藏收藏0 支持支持0 反对反对0 微信微信
    13#
    无效楼层,该帖已经被删除

    0

    主题

    13

    听众

    473

    积分

    升级  57.67%

  • TA的每日心情
    奋斗
    2016-7-31 18:28
  • 签到天数: 72 天

    [LV.6]常住居民II

    自我介绍
    本科生

    社区QQ达人

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

    回复

    使用道具 举报

    NICE--LLL        

    0

    主题

    6

    听众

    98

    积分

    升级  97.89%

  • TA的每日心情
    奋斗
    2014-7-1 21:02
  • 签到天数: 31 天

    [LV.5]常住居民I

    自我介绍
    WU
    回复

    使用道具 举报

    17

    主题

    7

    听众

    216

    积分

    升级  58%

  • TA的每日心情
    无聊
    2014-6-11 16:12
  • 签到天数: 43 天

    [LV.5]常住居民I

    自我介绍
    河南大学2010级学生,国际经济与贸易

    新人进步奖

    群组数学建摸协会

    群组小草的客厅

    回复

    使用道具 举报

    4

    主题

    4

    听众

    1133

    积分

    升级  13.3%

  • TA的每日心情
    开心
    2013-5-14 11:03
  • 签到天数: 167 天

    [LV.7]常住居民III

    群组数学建摸协会

    群组数学建模

    群组西安交大数学建模

    群组2011建模讨论组

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

    回复

    使用道具 举报

    4

    主题

    4

    听众

    1133

    积分

    升级  13.3%

  • TA的每日心情
    开心
    2013-5-14 11:03
  • 签到天数: 167 天

    [LV.7]常住居民III

    群组数学建摸协会

    群组数学建模

    群组西安交大数学建模

    群组2011建模讨论组

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

    回复

    使用道具 举报

    4

    主题

    4

    听众

    1133

    积分

    升级  13.3%

  • TA的每日心情
    开心
    2013-5-14 11:03
  • 签到天数: 167 天

    [LV.7]常住居民III

    群组数学建摸协会

    群组数学建模

    群组西安交大数学建模

    群组2011建模讨论组

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

    回复

    使用道具 举报

    xueyues3        

    0

    主题

    0

    听众

    4

    积分

    升级  80%

    该用户从未签到

    路过!!!!!!& {5 E2 C. a% n1 q4 x# Y$ D
    2 e6 m0 F- }& U% x

    , ?' p; o. m( l  e2 H; m" Z3 J
    - C# ~2 U! G, @1 H0 e, j- y; J1 x- ^3 }. X1 E5 j
    ; S6 X1 O- V9 ]

    + S* o7 f6 d# A" G9 L4 m/ v. \: k1 z, h

    5 c- \2 J- K' P) k: q; y& e2 D4 n( [" `! a7 }+ \5 d( Z! e2 D+ m

    0 o- D, k4 ?# b( A9 K
    2 P0 G( g# i' Z/ i! R+ U. S4 O, [4 X- n- R0 u  {
    高密招聘网 gmzhaopin.com
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-4-15 07:59 , Processed in 0.531937 second(s), 100 queries .

    回顶部