QQ登录

只需要一步,快速开始

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

我是新手,MATLAB接触不久。有一题不会,求教求教

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

1

主题

8

听众

66

积分

升级  64.21%

  • TA的每日心情
    奋斗
    2013-9-13 10:09
  • 签到天数: 17 天

    [LV.4]偶尔看看III

    自我介绍
    我一直在路上。努力~

    群组Matlab讨论组

    跳转到指定楼层
    1#
    发表于 2013-7-15 16:44 |只看该作者 |正序浏览
    |招呼Ta 关注Ta
    对于如下线性规划问题(有3n个决策变量(x,r,s)和2n个约束)- b) ^. A5 l" Z
    main(-xn)1 }, g: Q0 |& f! H$ E$ o
    st
    4 X7 S, F3 F) Z% j  u4 G4x1-4r1=1;
    ( v, }5 e3 e. Y9 L! w" w' w$ V3 nx1+s1=1;% Y2 H6 ]* N5 A
    4xj-x(j-1)-4rj=0;(j=2,3,4......)
    9 n3 U3 j2 j! w0 z4 H4xj+x(j-1)+4sj=0;(j=2,3,4......)% `" ?. R7 `  F1 ~- b/ h+ P
    xj,rj,sj>=0;(j=2,3,4......)- {5 d" g# Q: u
    ' C; h# F; E* _
    请分别对n的不同取值(如n=2,10,50等),利用MATLAB软件求解上述规划。
    , N; s$ M5 ]. S$ k/ D其中j,n是下标 。
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏1 支持支持0 反对反对0 微信微信

    0

    主题

    11

    听众

    23

    积分

    升级  18.95%

  • TA的每日心情
    擦汗
    2014-9-14 19:04
  • 签到天数: 8 天

    [LV.3]偶尔看看II

    自我介绍
    喜爱数学建模
    回复

    使用道具 举报

    0

    主题

    7

    听众

    37

    积分

    升级  33.68%

  • TA的每日心情
    开心
    2014-1-29 21:12
  • 签到天数: 11 天

    [LV.3]偶尔看看II

    自我介绍
    一名医学生,对数学建模很感兴趣,目前正在参加培训,还望各位大神多多指教
    回复

    使用道具 举报

    玖玖        

    1

    主题

    8

    听众

    66

    积分

    升级  64.21%

  • TA的每日心情
    奋斗
    2013-9-13 10:09
  • 签到天数: 17 天

    [LV.4]偶尔看看III

    自我介绍
    我一直在路上。努力~

    群组Matlab讨论组

        syms n
    6 {$ U4 u+ K3 B# \9 H( mn=input('input n please.(press enter)n=')% G* H; M$ I' K& E
    if(n==1)( F' x7 N. `( j0 u0 |' ?! @
        a=[4,-4,0;1,0,1];% X" i) I/ j$ M+ l
        b=[1;1];
    * V0 k0 W3 D: }0 D3 g& T5 B( r    v=zeros(1,3);
    - k% A, g7 K( P    c=[-1;0;0];& |- |  Z) M# v* _
        [x,f]=linprog(c,[],[],a,b,v)
    + K: J! j% \* E7 b: ^+ k$ Q5 welse2 J% T' E3 k* h* X2 s  t) x# o0 S- f
        a=zeros(2*n,3*n);( S/ f  _* H0 u: J! o% t* \. M$ j
        for j=2:n" i$ E: v7 i- H: N8 R/ G0 I
                       x=zeros(1,n);
    3 L3 J' S8 e4 ~: `                   r=zeros(1,n);
    ; p- J8 f* s$ r2 }) K# y                   s=zeros(1,n);
    5 l0 j6 P7 H! L# p2 d  u) m                   x(j)=4;2 n7 `/ a3 \+ ^. r' b
                       x(j-1)=-1;
    ' m2 \  F' l& V' F" b8 ~                   r(j)=-4;
    8 p1 y6 z! Z* |: l+ z8 O" g" N                   a((j+1),=[x,r,s];2 x* l- D" E6 N3 I4 v
                       x(j)=4;6 ]" i4 G6 T. Z
                       x(j-1)=1;5 p/ u6 @5 A+ b7 R8 w1 w
                       r(j)=0;
    8 e; x' y. ^! \# O                   s(j)=4;# d, @2 j2 G3 M7 q7 ~2 ?* M7 |
                       a((j+n),=[x,r,s];
    : w$ U; P$ ^$ v" A4 h0 U    end
    # K8 y3 y0 O2 Y3 ]" @0 [( C9 c7 \/ Q0 _               v=zeros(1,3*n);
    ) n+ e. ?5 O" V. X              c1=-ones(n,1);
    * {5 ]8 t$ H) }9 z" s( j              c2=zeros(2*n,1);
    2 b4 @  ~& Z2 o$ Y              c=[c1;c2];
    ' @: x4 Y. f9 p               x=zeros(1,n);/ @1 e& z% W6 O/ b" i% y/ N
                   r=zeros(1,n);! C" \" t6 t5 [' }  X
                   s=zeros(1,n);
      k3 s6 J! O: E9 K# `               x(1)=4;
    ' o4 Q; Q% X& O" u5 g) Y               r(1)=-4;2 G+ G7 V/ k! A' W) |9 U
                   a(1,=[x,r,s];
    % B) q" l$ z/ Q) |5 G5 H& @               x(1)=1;; t2 U$ C% M7 D! N0 q0 g1 f
                   r(1)=0;
    0 I* {  ^+ G! _& ]               s(1)=1;
    : z1 G: b$ w' i4 e2 w% s               a(2,=[x,r,s];
    : o4 s  f8 J* T0 |. D               b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
    1 x5 F$ p2 F4 N7 w) z" n% D! I( d               [x,f]=linprog(c,[],[],a,b,v)
    - _! @: q$ I7 r6 i- ]2 D  k  Y9 send
    回复

    使用道具 举报

    玖玖        

    1

    主题

    8

    听众

    66

    积分

    升级  64.21%

  • TA的每日心情
    奋斗
    2013-9-13 10:09
  • 签到天数: 17 天

    [LV.4]偶尔看看III

    自我介绍
    我一直在路上。努力~

    群组Matlab讨论组

    回复

    使用道具 举报

    玖玖        

    1

    主题

    8

    听众

    66

    积分

    升级  64.21%

  • TA的每日心情
    奋斗
    2013-9-13 10:09
  • 签到天数: 17 天

    [LV.4]偶尔看看III

    自我介绍
    我一直在路上。努力~

    群组Matlab讨论组

    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-8-4 15:06 , Processed in 0.432120 second(s), 77 queries .

    回顶部